Protein coding variation in the J:ARC and J:DO outbred laboratory mouse stocks provides a molecular basis for distinct research applications.

Outbred laboratory mice (Mus musculus) are readily available and have high fecundity, making them a popular choice in biomedical research, especially toxicological and pharmacological applications. Direct high throughput genome sequencing (HTS) of these widely used research animals is an important genetic quality control measure that enhances research reproducibility. HTS data have been used to confirm the common origin of outbred stocks and to molecularly define distinct outbred populations. But these data have also revealed unexpected population structure and homozygosity in some populations; genetic features that emerge when outbred stocks are not properly maintained. We used exome sequencing to discover and interrogate protein-coding variation in a newly established population of Swiss-derived outbred stock (J:ARC) that is closely related to other, commonly used CD-1 outbred populations. We used these data to describe the genetic architecture of the J:ARC population including heterozygosity, minor allele frequency, LD decay, and we defined novel, protein-coding sequence variation. These data reveal the expected genetic architecture for a properly maintained outbred stock and provide a basis for the on-going genetic quality control. We also compared these data to protein-coding variation found in a multiparent outbred stock, the Diversity Outbred (J:DO). We found that the more recently derived, multiparent outbred stock has significantly higher interindividual variability, greater overall genetic variation, higher heterozygosity, and fewer novel variants than the Swiss-derived J:ARC stock. However, among the novel variants found in the J:DO stock, significantly more are predicted to be protein-damaging. The fact that individuals from this population can tolerate a higher load of potentially damaging variants highlights the buffering effects of allelic diversity and the differing selective pressures in these stocks. While both outbred stocks offer significant individual heterozygosity, our data provide a molecular basis for their intended applications, where the J:DO are best suited for studies requiring maximum, population-level genetic diversity and power for mapping, while the J:ARC are best suited as a general-purpose outbred stock with robust fecundity, relatively low allelic diversity, and less potential for extreme phenotypic variability.

A bipartite element with allele-specific functions safeguards DNA methylation imprints at the Dlk1-Dio3 locus.

Loss of imprinting (LOI) results in severe developmental defects, but the mechanisms preventing LOI remain incompletely understood. Here, we dissect the functional components of the imprinting control region of the essential Dlk1-Dio3 locus (called IG-DMR) in pluripotent stem cells. We demonstrate that the IG-DMR consists of two antagonistic elements: a paternally methylated CpG island that prevents recruitment of TET dioxygenases and a maternally unmethylated non-canonical enhancer that ensures expression of the Gtl2 lncRNA by counteracting de novo DNA methyltransferases. Genetic or epigenetic editing of these elements leads to distinct LOI phenotypes with characteristic alternations of allele-specific gene expression, DNA methylation, and 3D chromatin topology. Although repression of the Gtl2 promoter results in dysregulated imprinting, the stability of LOI phenotypes depends on the IG-DMR, suggesting a functional hierarchy. These findings establish the IG-DMR as a bipartite control element that maintains imprinting by allele-specific restriction of the DNA (de)methylation machinery.

Mouse Genetic Reference Populations: Cellular Platforms for Integrative Systems Genetics.

Interrogation of disease-relevant cellular and molecular traits exhibited by genetically diverse cell populations enables in vitro systems genetics approaches for uncovering the basic properties of cellular function and identity. Primary cells, stem cells, and organoids derived from genetically diverse mouse strains, such as Collaborative Cross and Diversity Outbred populations, offer the opportunity for parallel in vitro/in vivo screening. These panels provide genetic resolution for variant discovery and functional characterization, as well as disease modeling and in vivo validation capabilities. Here we review mouse cellular systems genetics approaches for characterizing the influence of genetic variation on signaling networks and phenotypic diversity, and we discuss approaches for data integration and cross-species validation.

Context-Dependent Requirement of Euchromatic Histone Methyltransferase Activity during Reprogramming to Pluripotency.

Methylation of histone 3 at lysine 9 (H3K9) constitutes a roadblock for cellular reprogramming. Interference with methyltransferases or activation of demethylases by the cofactor ascorbic acid (AA) facilitates the derivation of induced pluripotent stem cells (iPSCs), but possible interactions between specific methyltransferases and AA treatment remain insufficiently explored. We show that chemical inhibition of the methyltransferases EHMT1 and EHMT2 counteracts iPSC formation in an enhanced reprogramming system in the presence of AA, an effect that is dependent on EHMT1. EHMT inhibition during enhanced reprogramming is associated with rapid loss of H3K9 dimethylation, inefficient downregulation of somatic genes, and failed mesenchymal-to-epithelial transition. Furthermore, transient EHMT inhibition during reprogramming yields iPSCs that fail to efficiently give rise to viable mice upon blastocyst injection. Our observations establish novel functions of H3K9 methyltransferases and suggest that a functional balance between AA-stimulated enzymes and EHMTs supports efficient and less error-prone iPSC reprogramming to pluripotency.

A Susceptibility Locus on Chromosome 13 Profoundly Impacts the Stability of Genomic Imprinting in Mouse Pluripotent Stem Cells.

Cultured pluripotent cells accumulate detrimental chromatin alterations, including DNA methylation changes at imprinted genes known as loss of imprinting (LOI). Although the occurrence of LOI is considered a stochastic phenomenon, here we document a genetic determinant that segregates mouse pluripotent cells into stable and unstable cell lines. Unstable lines exhibit hypermethylation at Dlk1-Dio3 and other imprinted loci, in addition to impaired developmental potential. Stimulation of demethylases by ascorbic acid prevents LOI and loss of developmental potential. Susceptibility to LOI greatly differs between commonly used mouse strains, which we use to map a causal region on chromosome 13 with quantitative trait locus (QTL) analysis. Our observations identify a strong genetic determinant of locus-specific chromatin abnormalities in pluripotent cells and provide a non-invasive way to suppress them. This highlights the importance of considering genetics in conjunction with culture conditions for assuring the quality of pluripotent cells for biomedical applications.

Analysis of 3D genomic interactions identifies candidate host genes that transposable elements potentially regulate.

BACKGROUND: The organization of chromatin in the nucleus plays an essential role in gene regulation. About half of the mammalian genome comprises transposable elements. Given their repetitive nature, reads associated with these elements are generally discarded or randomly distributed among elements of the same type in genome-wide analyses. Thus, it is challenging to identify the activities and properties of individual transposons. As a result, we only have a partial understanding of how transposons contribute to chromatin folding and how they impact gene regulation. RESULTS: Using PCR and Capture-based chromosome conformation capture (3C) approaches, collectively called 4Tran, we take advantage of the repetitive nature of transposons to capture interactions from multiple copies of endogenous retrovirus (ERVs) in the human and mouse genomes. With 4Tran-PCR, reads are selectively mapped to unique regions in the genome. This enables the identification of transposable element interaction profiles for individual ERV families and integration events specific to particular genomes. With this approach, we demonstrate that transposons engage in long-range intra-chromosomal interactions guided by the separation of chromosomes into A and B compartments as well as topologically associated domains (TADs). In contrast to 4Tran-PCR, Capture-4Tran can uniquely identify both ends of an interaction that involve retroviral repeat sequences, providing a powerful tool for uncovering the individual transposable element insertions that interact with and potentially regulate target genes. CONCLUSIONS: 4Tran provides new insight into the manner in which transposons contribute to chromosome architecture and identifies target genes that transposable elements can potentially control.

Trim32 Deficiency Enhances Th2 Immunity and Predisposes to Features of Atopic Dermatitis.

Altered innate immunity is a feature of certain skin inflammatory diseases such as psoriasis and atopic dermatitis (AD). In this study, we provide evidence that deficiency in Trim32 (a tripartite motif [TRIM] protein with innate antiviral activity) contributes to a T helper type 2 biased response and predisposes to features of AD in mice. On treatment with the toll-like receptor 7 agonist imquimod (IMQ), Trim32 knockout mice displayed compromised psoriasiform phenotypes and defective T helper type 17 response. Instead, IMQ treatment of Trim32 knockout mice induced AD-like phenotypes with enhanced skin infiltration of eosinophils and mast cells, elevation of T helper type 2 cytokines/chemokines expression, and reduced expression of filaggrin protein expression. Furthermore, although the induction of phosphorylated Stat3 and RelA was compromised after IMQ treatment in the knockout mice, phosphorylated Stat6 was elevated. CC chemokine ligand 20 induction by tumor necrosis factor-α and IL-17A was reduced in Trim32-deficient keratinocytes, whereas CC chemokine ligand 5 induction by tumor necrosis factor-α and IL-4 was enhanced. In addition, Trim32 protein levels were elevated in mice treated with IMQ. Unlike Trim32 overexpression in psoriasis, TRIM32 levels were low in patients with AD. Based on Trim32 induction by IMQ, the lower levels of TRIM32 in AD skin compared with healthy control and psoriatic skin suggest a defective TRIM32 pathway in AD pathogenesis.

A reporter model to visualize imprinting stability at the Dlk1 locus during mouse development and in pluripotent cells.

Genomic imprinting results in the monoallelic expression of genes that encode important regulators of growth and proliferation. Dysregulation of imprinted genes, such as those within the Dlk1-Dio3 locus, is associated with developmental syndromes and specific diseases. Our ability to interrogate causes of imprinting instability has been hindered by the absence of suitable model systems. Here, we describe a Dlk1 knock-in reporter mouse that enables single-cell visualization of allele-specific expression and prospective isolation of cells, simultaneously. We show that this 'imprinting reporter mouse' can be used to detect tissue-specific Dlk1 expression patterns in developing embryos. We also apply this system to pluripotent cell culture and demonstrate that it faithfully indicates DNA methylation changes induced upon cellular reprogramming. Finally, the reporter system reveals the role of elevated oxygen levels in eroding imprinted Dlk1 expression during prolonged culture and in vitro differentiation. The possibility to study allele-specific expression in different contexts makes our reporter system a useful tool to dissect the regulation of genomic imprinting in normal development and disease.

A Damage-Independent Role for 53BP1 that Impacts Break Order and Igh Architecture during Class Switch Recombination.

During class switch recombination (CSR), B cells replace the Igh Cµ or δ exons with another downstream constant region exon (CH), altering the antibody isotype. CSR occurs through the introduction of AID-mediated double-strand breaks (DSBs) in switch regions and subsequent ligation of broken ends. Here, we developed an assay to investigate the dynamics of DSB formation in individual cells. We demonstrate that the upstream switch region Sµ is first targeted during recombination and that the mechanism underlying this control relies on 53BP1. Surprisingly, regulation of break order occurs through residual binding of 53BP1 to chromatin before the introduction of damage and independent of its established role in DNA repair. Using chromosome conformation capture, we show that 53BP1 mediates changes in chromatin architecture that affect break order. Finally, our results explain how changes in Igh architecture in the absence of 53BP1 could promote inversional rearrangements that compromise CSR.

4C-ker: A Method to Reproducibly Identify Genome-Wide Interactions Captured by 4C-Seq Experiments.

4C-Seq has proven to be a powerful technique to identify genome-wide interactions with a single locus of interest (or "bait") that can be important for gene regulation. However, analysis of 4C-Seq data is complicated by the many biases inherent to the technique. An important consideration when dealing with 4C-Seq data is the differences in resolution of signal across the genome that result from differences in 3D distance separation from the bait. This leads to the highest signal in the region immediately surrounding the bait and increasingly lower signals in far-cis and trans. Another important aspect of 4C-Seq experiments is the resolution, which is greatly influenced by the choice of restriction enzyme and the frequency at which it can cut the genome. Thus, it is important that a 4C-Seq analysis method is flexible enough to analyze data generated using different enzymes and to identify interactions across the entire genome. Current methods for 4C-Seq analysis only identify interactions in regions near the bait or in regions located in far-cis and trans, but no method comprehensively analyzes 4C signals of different length scales. In addition, some methods also fail in experiments where chromatin fragments are generated using frequent cutter restriction enzymes. Here, we describe 4C-ker, a Hidden-Markov Model based pipeline that identifies regions throughout the genome that interact with the 4C bait locus. In addition, we incorporate methods for the identification of differential interactions in multiple 4C-seq datasets collected from different genotypes or experimental conditions. Adaptive window sizes are used to correct for differences in signal coverage in near-bait regions, far-cis and trans chromosomes. Using several datasets, we demonstrate that 4C-ker outperforms all existing 4C-Seq pipelines in its ability to reproducibly identify interaction domains at all genomic ranges with different resolution enzymes.

Dermal reflectivity determined by optical coherence tomography is an indicator of epidermal hyperplasia and dermal edema within inflamed skin.

Psoriasis is a common inflammatory skin disease resulting from genetic and environmental alterations of cutaneous immune responses. While numerous therapeutic targets involved in the immunopathogenesis of psoriasis have been identified, the in vivo dynamics of inflammation in psoriasis remain unclear. We undertook in vivo time course focus-tracked optical coherence tomography (OCT) imaging to noninvasively document cutaneous alterations in mouse skin treated topically with Imiquimod (IMQ), an established model of a psoriasis-like disease. Quantitative appraisal of dermal architectural changes was achieved through a two parameter fit of OCT axial scans in the dermis of the form A(x, y, z) = ρ(x, y)exp [-µ(x, y)z]. Ensemble averaging over 2000 axial scans per mouse in each treatment arm revealed no significant changes in the average dermal attenuation rate, <µ>, however the average local dermal reflectivity <ρ>, decreased significantly following 1, 3, and 6 days of IMQ treatment (p < 0.001) in comparison to vehicle-treated control mice. In contrast, epidermal and dermal thickness changes were only significant when comparing controls and 6-day IMQ treated mice. This suggests that dermal alterations, attributed to collagen fiber bundle enlargement, occur prior to epidermal thickness changes due to hyperplasia and dermal thickness changes due to edema. Dermal reflectivity positively correlated with epidermal hyperplasia (r(epi)(2) = 0.78) and dermal edema (r(derm)(2) = 0.86). Our results suggest that dermal reflectivity as measured by OCT can be utilized to quantify a psoriasis-like disease in mice, and thus has the potential to aid in the quantitative assessment of psoriasis in humans.