Unraveling the genetics of arsenic toxicity with cellular morphology QTL.

The health risks that arise from environmental exposures vary widely within and across human populations, and these differences are largely determined by genetic variation and gene-by-environment (gene-environment) interactions. However, risk assessment in laboratory mice typically involves isogenic strains and therefore, does not account for these known genetic effects. In this context, genetically heterogenous cell lines from laboratory mice are promising tools for population-based screening because they provide a way to introduce genetic variation in risk assessment without increasing animal use. Cell lines from genetic reference populations of laboratory mice offer genetic diversity, power for genetic mapping, and potentially, predictive value for in vivo experimentation in genetically matched individuals. To explore this further, we derived a panel of fibroblast lines from a genetic reference population of laboratory mice (the Diversity Outbred, DO). We then used high-content imaging to capture hundreds of cell morphology traits in cells exposed to the oxidative stress-inducing arsenic metabolite monomethylarsonous acid (MMAIII). We employed dose-response modeling to capture latent parameters of response and we then used these parameters to identify several hundred cell morphology quantitative trait loci (cmQTL). Response cmQTL encompass genes with established associations with cellular responses to arsenic exposure, including Abcc4 and Txnrd1, as well as novel gene candidates like Xrcc2. Moreover, baseline trait cmQTL highlight the influence of natural variation on fundamental aspects of nuclear morphology. We show that the natural variants influencing response include both coding and non-coding variation, and that cmQTL haplotypes can be used to predict response in orthogonal cell lines. Our study sheds light on the major molecular initiating events of oxidative stress that are under genetic regulation, including the NRF2-mediated antioxidant response, cellular detoxification pathways, DNA damage repair response, and cell death trajectories.

Genetic control of the pluripotency epigenome determines differentiation bias in mouse embryonic stem cells.

Genetically diverse pluripotent stem cells display varied, heritable responses to differentiation cues. Here, we harnessed these disparities through derivation of mouse embryonic stem cells from the BXD genetic reference panel, along with C57BL/6J (B6) and DBA/2J (D2) parental strains, to identify loci regulating cell state transitions. Upon transition to formative pluripotency, B6 stem cells quickly dissolved naïve networks adopting gene expression modules indicative of neuroectoderm lineages, whereas D2 retained aspects of naïve pluripotency. Spontaneous formation of embryoid bodies identified divergent differentiation where B6 showed a propensity toward neuroectoderm and D2 toward definitive endoderm. Genetic mapping identified major trans-acting loci co-regulating chromatin accessibility and gene expression in both naïve and formative pluripotency. These loci distally modulated occupancy of pluripotency factors at hundreds of regulatory elements. One trans-acting locus on Chr 12 primarily impacted chromatin accessibility in embryonic stem cells, while in epiblast-like cells, the same locus subsequently influenced expression of genes enriched for neurogenesis, suggesting early chromatin priming. These results demonstrate genetically determined biases in lineage commitment and identify major regulators of the pluripotency epigenome.

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.

The spontaneous mouse mutant low set ears (Lse) is caused by tandem duplication of Fgf3 and Fgf4.

The external ear develops from an organized convergence of ventrally migrating neural crest cells into the first and second branchial arches. Defects in external ear position are often symptomatic of complex syndromes such as Apert, Treacher-Collins, and Crouzon Syndrome. The low set ears (Lse) spontaneous mouse mutant is characterized by the dominant inheritance of a ventrally shifted external ear position and an abnormal external auditory meatus (EAM). We identified the causative mutation as a 148 Kb tandem duplication on Chromosome 7, which includes the entire coding sequences of Fgf3 and Fgf4. Duplications of FGF3 and FGF4 occur in 11q duplication syndrome in humans and are associated with craniofacial anomalies, among other features. Intercrosses of Lse-affected mice revealed perinatal lethality in homozygotes, and Lse/Lse embryos display additional phenotypes including polydactyly, abnormal eye morphology, and cleft secondary palate. The duplication results in increased Fgf3 and Fgf4 expression in the branchial arches and additional discrete domains in the developing embryo. This ectopic overexpression resulted in functional FGF signaling, demonstrated by increased Spry2 and Etv5 expression in overlapping domains of the developing arches. Finally, a genetic interaction between Fgf3/4 overexpression and Twist1, a regulator of skull suture development, resulted in perinatal lethality, cleft palate, and polydactyly in compound heterozygotes. These data indicate a role for Fgf3 and Fgf4 in external ear and palate development and provide a novel mouse model for further interrogation of the biological consequences of human FGF3/4 duplication.

High-Diversity Mouse Populations for Complex Traits.

Contemporary mouse genetic reference populations are a powerful platform to discover complex disease mechanisms. Advanced high-diversity mouse populations include the Collaborative Cross (CC) strains, Diversity Outbred (DO) stock, and their isogenic founder strains. When used in systems genetics and integrative genomics analyses, these populations efficiently harnesses known genetic variation for precise and contextualized identification of complex disease mechanisms. Extensive genetic, genomic, and phenotypic data are already available for these high-diversity mouse populations and a growing suite of data analysis tools have been developed to support research on diverse mice. This integrated resource can be used to discover and evaluate disease mechanisms relevant across species.

Large-scale discovery of mouse transgenic integration sites reveals frequent structural variation and insertional mutagenesis.

Transgenesis has been a mainstay of mouse genetics for over 30 yr, providing numerous models of human disease and critical genetic tools in widespread use today. Generated through the random integration of DNA fragments into the host genome, transgenesis can lead to insertional mutagenesis if a coding gene or an essential element is disrupted, and there is evidence that larger scale structural variation can accompany the integration. The insertion sites of only a tiny fraction of the thousands of transgenic lines in existence have been discovered and reported, due in part to limitations in the discovery tools. Targeted locus amplification (TLA) provides a robust and efficient means to identify both the insertion site and content of transgenes through deep sequencing of genomic loci linked to specific known transgene cassettes. Here, we report the first large-scale analysis of transgene insertion sites from 40 highly used transgenic mouse lines. We show that the transgenes disrupt the coding sequence of endogenous genes in half of the lines, frequently involving large deletions and/or structural variations at the insertion site. Furthermore, we identify a number of unexpected sequences in some of the transgenes, including undocumented cassettes and contaminating DNA fragments. We demonstrate that these transgene insertions can have phenotypic consequences, which could confound certain experiments, emphasizing the need for careful attention to control strategies. Together, these data show that transgenic alleles display a high rate of potentially confounding genetic events and highlight the need for careful characterization of each line to assure interpretable and reproducible experiments.

Biallelic mutations in the ferredoxin reductase gene cause novel mitochondriopathy with optic atrophy.

Iron-sulfur (Fe-S) clusters are ubiquitous cofactors essential to various cellular processes, including mitochondrial respiration, DNA repair, and iron homeostasis. A steadily increasing number of disorders are being associated with disrupted biogenesis of Fe-S clusters. Here, we conducted whole-exome sequencing of patients with optic atrophy and other neurological signs of mitochondriopathy and identified 17 individuals from 13 unrelated families with recessive mutations in FDXR, encoding the mitochondrial membrane-associated flavoprotein ferrodoxin reductase required for electron transport from NADPH to cytochrome P450. In vitro enzymatic assays in patient fibroblast cells showed deficient ferredoxin NADP reductase activity and mitochondrial dysfunction evidenced by low oxygen consumption rates (OCRs), complex activities, ATP production and increased reactive oxygen species (ROS). Such defects were rescued by overexpression of wild-type FDXR. Moreover, we found that mice carrying a spontaneous mutation allelic to the most common mutation found in patients displayed progressive gait abnormalities and vision loss, in addition to biochemical defects consistent with the major clinical features of the disease. Taken together, these data provide the first demonstration that germline, hypomorphic mutations in FDXR cause a novel mitochondriopathy and optic atrophy in humans.

Exome sequencing reveals pathogenic mutations in 91 strains of mice with Mendelian disorders.

Spontaneously arising mouse mutations have served as the foundation for understanding gene function for more than 100 years. We have used exome sequencing in an effort to identify the causative mutations for 172 distinct, spontaneously arising mouse models of Mendelian disorders, including a broad range of clinically relevant phenotypes. To analyze the resulting data, we developed an analytics pipeline that is optimized for mouse exome data and a variation database that allows for reproducible, user-defined data mining as well as nomination of mutation candidates through knowledge-based integration of sample and variant data. Using these new tools, putative pathogenic mutations were identified for 91 (53%) of the strains in our study. Despite the increased power offered by potentially unlimited pedigrees and controlled breeding, about half of our exome cases remained unsolved. Using a combination of manual analyses of exome alignments and whole-genome sequencing, we provide evidence that a large fraction of unsolved exome cases have underlying structural mutations. This result directly informs efforts to investigate the similar proportion of apparently Mendelian human phenotypes that are recalcitrant to exome sequencing.

Biallelic mutations in FDXR cause neurodegeneration associated with inflammation.

Mitochondrial dysfunction lies behind many neurodegenerative disorders, owing largely to the intense energy requirements of most neurons. Such mitochondrial dysfunction may work through a variety of mechanisms, from direct disruption of the electron transport chain to abnormal mitochondrial biogenesis. Recently, we have identified biallelic mutations in the mitochondrial flavoprotein "ferredoxin reductase" (FDXR) gene as a novel cause of mitochondriopathy, peripheral neuropathy, and optic atrophy. In this report, we expand upon those results by describing two new cases of disease-causing FDXR variants in patients with variable severity of phenotypes, including evidence of an inflammatory response in brain autopsy. To investigate the underlying pathogenesis, we examined neurodegeneration in a mouse model. We found that Fdxr mutant mouse brain tissues share pathological changes similar to those seen in patient autopsy material, including increased astrocytes. Furthermore, we show that these abnormalities are associated with increased levels of markers for both neurodegeneration and gliosis, with the latter implying inflammation as a major factor in the pathology of Fdxr mutations. These data provide further insight into the pathogenic mechanism of FDXR-mediated central neuropathy, and suggest an avenue for mechanistic studies that will ultimately inform treatment.

Discovery and characterization of spontaneous mouse models of craniofacial dysmorphology.

Craniofacial abnormalities are among the most common features of human genetic syndromes and disorders. The etiology of these conditions is often complex, influenced by both genetic context and the environment. Frequently, craniofacial abnormalities present as part of a syndrome with clear comorbid phenotypes, providing additional insight into mechanisms of the causative gene or pathway. The mouse has been a key tool in our understanding of the genetic mechanisms of craniofacial development and disease, and can provide excellent models for human craniofacial abnormalities. While powerful genetic engineering tools in the mouse have contributed significantly our understanding of craniofacial development and dysmorphology, forward genetic approaches provide an unbiased means to identify new genes and pathways. Moreover, spontaneous mutations can occur on any number of genetic backgrounds, potentially revealing critical genes that require a specific genetic context. Here we report discovery and phenotyping of 43 craniofacial mouse models, derived primarily from a screen for spontaneous mutations in production colonies at the Jackson Laboratory. We identify the causative gene for 33 lines, including novel genes in pathways not previously connected to craniofacial development, and novel alleles of known genes that present with unique phenotypes. Together with our detailed characterization, this work provides a valuable gene discovery resource for the craniofacial community, and a rich source of mouse models for further investigation.

CRISPRtools: a flexible computational platform for performing CRISPR/Cas9 experiments in the mouse.

Genome editing using the CRISPR/Cas9 RNA-guided endonuclease system has rapidly become a driving force for discovery in modern biomedical research. This simple yet elegant system has been widely used to generate both loss-of-function alleles and precision knock-in mutations using single-stranded donor oligonucleotides. Our CRISPRtools platform supports both of these applications in order to facilitate the use of CRISPR/Cas9. While there are several tools that facilitate CRISPR/Cas9 design and screen for potential off-target sites, the process is typically performed sequentially on single genes, limiting scalability for large-scale programs. Here, the design principle underlying gene ablation is based upon using paired guides flanking a critical region/exon of interest to create deletions. Guide pairs are rank ordered based upon published efficiency scores and off-target analyses, and reported in a concise format for downstream implementation. The exon deletion strategy simplifies characterization of founder animals and is the strategy employed for the majority of knockouts in the mouse. In proof-of-principle experiments, the effectiveness of this approach is demonstrated using microinjection and electroporation to introduce CRISPR/Cas9 components into mouse zygotes to delete critical exons.

Mouse genomic variation and its effect on phenotypes and gene regulation.

We report genome sequences of 17 inbred strains of laboratory mice and identify almost ten times more variants than previously known. We use these genomes to explore the phylogenetic history of the laboratory mouse and to examine the functional consequences of allele-specific variation on transcript abundance, revealing that at least 12% of transcripts show a significant tissue-specific expression bias. By identifying candidate functional variants at 718 quantitative trait loci we show that the molecular nature of functional variants and their position relative to genes vary according to the effect size of the locus. These sequences provide a starting point for a new era in the functional analysis of a key model organism.

Kif18a is specifically required for mitotic progression during germ line development.

Genome integrity in the developing germ line is strictly required for fecundity. In proliferating somatic cells and in germ cells, there are mitotic checkpoint mechanisms that ensure accurate chromosome segregation and euploidy. There is growing evidence of mitotic cell cycle components that are uniquely required in the germ line to ensure genome integrity. We previously showed that the primary phenotype of germ cell deficient 2 (gcd2) mutant mice is infertility due to germ cell depletion during embryogenesis. Here we show that the underlying mutation is a mis-sense mutation, R308K, in the motor domain of the kinesin-8 family member, KIF18A, a protein that is expressed in a variety of proliferative tissues and is a key regulator of chromosome alignment during mitosis. Despite the conservative nature of the mutation, we show that its functional consequences are equivalent to KIF18A deficiency in HeLa cells. We also show that somatic cells progress through mitosis, despite having chromosome alignment defects, while germ cells with similar chromosome alignment defects undergo mitotic arrest and apoptosis. Our data provide evidence for differential requirements for chromosome alignment in germ and somatic cells and show that Kif18a is one of a growing number of genes that are specifically required for cell cycle progression in proliferating germ cells.

Nuclear localization of PRDM9 and its role in meiotic chromatin modifications and homologous synapsis.

Developmental progress of germ cells through meiotic phases is closely tied to ongoing meiotic recombination. In mammals, recombination preferentially occurs in genomic regions known as hotspots; the protein that activates these hotspots is PRDM9, containing a genetically variable zinc finger (ZNF) domain and a PR-SET domain with histone H3K4 trimethyltransferase activity. PRDM9 is required for fertility in mice, but little is known about its localization and developmental dynamics. Application of spermatogenic stage-specific markers demonstrates that PRDM9 accumulates in male germ cell nuclei at pre-leptonema to early leptonema but is no longer detectable in nuclei by late zygonema. By the pachytene stage, PRDM9-dependent histone H3K4 trimethyl marks on hotspots also disappear. PRDM9 localizes to nuclei concurrently with the deposition of meiotic cohesin complexes, but is not required for incorporation of cohesin complex proteins into chromosomal axial elements, or accumulation of normal numbers of RAD51 foci on meiotic chromatin by late zygonema. Germ cells lacking PRDM9 exhibit inefficient homology recognition and synapsis, with aberrant repair of meiotic DNA double-strand breaks and transcriptional abnormalities characteristic of meiotic silencing of unsynapsed chromatin. Together, these results on the developmental time course for nuclear localization of PRDM9 establish its direct window of function and demonstrate the independence of chromosome axial element formation from the concurrent PRDM9-mediated activation of recombination hotspots.

Dsp rul: a spontaneous mouse mutation in desmoplakin as a model of Carvajal-Huerta syndrome.

Studies of spontaneous mutations in mice have provided valuable disease models and important insights into the mechanisms of human disease. Ruffled (rul) is a new autosomal recessive mutation causing abnormal hair coat in mice. The rul allele arose spontaneously in the RB156Bnr/EiJ inbred mouse strain. In addition to an abnormal coat texture, we found diffuse epidermal blistering, abnormal electrocardiograms (ECGs), and ventricular fibrosis in mutant animals. Using high-throughput sequencing (HTS) we found a frameshift mutation at 38,288,978bp of chromosome 13 in the desmoplakin gene (Dsp). The predicted mutant protein is truncated at the c-terminus and missing the majority of the plakin repeat domain. The phenotypes found in Dsp(rul) mice closely model a rare human disorder, Carvajal-Huerta syndrome. Carvajal-Huerta syndrome (CHS) is a rare cardiocutaneous disorder that presents in humans with wooly hair, palmoplantar keratoderma and ventricular cardiomyopathy. CHS results from an autosomal recessive mutation on the 3' end of desmoplakin (DSP) truncating the full length protein. The Dsp(rul) mouse provides a new model to investigate the pathogenesis of CHS, as well as the underlying basic biology of the adhesion molecules coded by the desmosomal genes.

Discovery of transgene insertion sites by high throughput sequencing of mate pair libraries.

BACKGROUND: Transgenesis by random integration of a transgene into the genome of a zygote has become a reliable and powerful method for the creation of new mouse strains that express exogenous genes, including human disease genes, tissue specific reporter genes or genes that allow for tissue specific recombination. Nearly 6,500 transgenic alleles have been created by random integration in embryos over the last 30 years, but for the vast majority of these strains, the transgene insertion sites remain uncharacterized. RESULTS: To obtain a complete understanding of how insertion sites might contribute to phenotypic outcomes, to more cost effectively manage transgenic strains, and to fully understand mechanisms of instability in transgene expression, we've developed methodology and a scoring scheme for transgene insertion site discovery using high throughput sequencing data. CONCLUSIONS: Similar to other molecular approaches to transgene insertion site discovery, high-throughput sequencing of standard paired-end libraries is hindered by low signal to noise ratios. This problem is exacerbated when the transgene consists of sequences that are also present in the host genome. We've found that high throughput sequencing data from mate-pair libraries are more informative when compared to data from standard paired end libraries. We also show examples of the genomic regions that harbor transgenes, which have in common a preponderance of repetitive sequences.

Imputation of 3D genome structure by genetic-epigenetic interaction modeling in mice.

Gene expression is known to be affected by interactions between local genetic variation and DNA accessibility, with the latter organized into three-dimensional chromatin structures. Analyses of these interactions have previously been limited, obscuring their regulatory context, and the extent to which they occur throughout the genome. Here, we undertake a genome-scale analysis of these interactions in a genetically diverse population to systematically identify global genetic-epigenetic interaction, and reveal constraints imposed by chromatin structure. We establish the extent and structure of genotype-by-epigenotype interaction using embryonic stem cells derived from Diversity Outbred mice. This mouse population segregates millions of variants from eight inbred founders, enabling precision genetic mapping with extensive genotypic and phenotypic diversity. With 176 samples profiled for genotype, gene expression, and open chromatin, we used regression modeling to infer genetic-epigenetic interactions on a genome-wide scale. Our results demonstrate that statistical interactions between genetic variants and chromatin accessibility are common throughout the genome. We found that these interactions occur within the local area of the affected gene, and that this locality corresponds to topologically associated domains (TADs). The likelihood of interaction was most strongly defined by the three-dimensional (3D) domain structure rather than linear DNA sequence. We show that stable 3D genome structure is an effective tool to guide searches for regulatory elements and, conversely, that regulatory elements in genetically diverse populations provide a means to infer 3D genome structure. We confirmed this finding with CTCF ChIP-seq that revealed strain-specific binding in the inbred founder mice. In stem cells, open chromatin participating in the most significant regression models demonstrated an enrichment for developmental genes and the TAD-forming CTCF-binding complex, providing an opportunity for statistical inference of shifting TAD boundaries operating during early development. These findings provide evidence that genetic and epigenetic factors operate within the context of 3D chromatin structure.

An in vitro neurogenetics platform for precision disease modeling in the mouse.

The power and scope of disease modeling can be markedly enhanced through the incorporation of broad genetic diversity. The introduction of pathogenic mutations into a single inbred mouse strain sometimes fails to mimic human disease. We describe a cross-species precision disease modeling platform that exploits mouse genetic diversity to bridge cell-based modeling with whole organism analysis. We developed a universal protocol that permitted robust and reproducible neural differentiation of genetically diverse human and mouse pluripotent stem cell lines and then carried out a proof-of-concept study of the neurodevelopmental gene DYRK1A. Results in vitro reliably predicted the effects of genetic background on Dyrk1a loss-of-function phenotypes in vivo. Transcriptomic comparison of responsive and unresponsive strains identified molecular pathways conferring sensitivity or resilience to Dyrk1a1A loss and highlighted differential messenger RNA isoform usage as an important determinant of response. This cross-species strategy provides a powerful tool in the functional analysis of candidate disease variants identified through human genetic studies.

Mapping the genetic landscape establishing a tumor immune microenvironment favorable for anti-PD-1 response in mice and humans.

Identifying host genetic factors modulating immune checkpoint inhibitor (ICI) efficacy has been experimentally challenging because of variations in both host and tumor genomes, differences in the microbiome, and patient life exposures. Utilizing the Collaborative Cross (CC) multi-parent mouse genetic resource population, we developed an approach that fixes the tumor genomic configuration while varying host genetics. With this approach, we discovered that response to anti-PD-1 (aPD1) immunotherapy was significantly heritable in four distinct murine tumor models (H2 between 0.18-0.40). For the MC38 colorectal carcinoma system (H2 = 0.40), we mapped four significant ICI response quantitative trait loci (QTL) localized to mouse chromosomes (mChr) 5, 9, 15 and 17, and identified significant epistatic interactions between specific QTL pairs. Differentially expressed genes within these QTL were highly enriched for immune genes and pathways mediating allograft rejection and graft vs host disease. Using a cross species analytical approach, we found a core network of 48 genes within the four QTLs that showed significant prognostic value for overall survival in aPD1 treated human cohorts that outperformed all other existing validated immunotherapy biomarkers, especially in human tumors of the previously defined immune subtype 4. Functional blockade of two top candidate immune targets within the 48 gene network, GM-CSF and high affinity IL-2/IL-15 signaling, completely abrogated the MC38 tumor transcriptional response to aPD1 therapy in vivo. Thus, we have established a powerful cross species in vivo platform capable of uncovering host genetic factors that establish the tumor immune microenvironment configuration propitious for ICI response.

Cell morphology QTL reveal gene by environment interactions in a genetically diverse cell population.

Genetically heterogenous cell lines from laboratory mice are promising tools for population-based screening as they offer power for genetic mapping, and potentially, predictive value for in vivo experimentation in genetically matched individuals. To explore this further, we derived a panel of fibroblast lines from a genetic reference population of laboratory mice (the Diversity Outbred, DO). We then used high-content imaging to capture hundreds of cell morphology traits in cells exposed to the oxidative stress-inducing arsenic metabolite monomethylarsonous acid (MMAIII). We employed dose-response modeling to capture latent parameters of response and we then used these parameters to identify several hundred cell morphology quantitative trait loci (cmQTL). Response cmQTL encompass genes with established associations with cellular responses to arsenic exposure, including Abcc4 and Txnrd1, as well as novel gene candidates like Xrcc2. Moreover, baseline trait cmQTL highlight the influence of natural variation on fundamental aspects of nuclear morphology. We show that the natural variants influencing response include both coding and non-coding variation, and that cmQTL haplotypes can be used to predict response in orthogonal cell lines. Our study sheds light on the major molecular initiating events of oxidative stress that are under genetic regulation, including the NRF2-mediated antioxidant response, cellular detoxification pathways, DNA damage repair response, and cell death trajectories.