A natural mutator allele shapes mutation spectrum variation in mice.

Although germline mutation rates and spectra can vary within and between species, common genetic modifiers of the mutation rate have not been identified in nature1. Here we searched for loci that influence germline mutagenesis using a uniquely powerful resource: a panel of recombinant inbred mouse lines known as the BXD, descended from the laboratory strains C57BL/6J (B haplotype) and DBA/2J (D haplotype). Each BXD lineage has been maintained by brother-sister mating in the near absence of natural selection, accumulating de novo mutations for up to 50 years on a known genetic background that is a unique linear mosaic of B and D haplotypes2. We show that mice inheriting D haplotypes at a quantitative trait locus on chromosome 4 accumulate C>A germline mutations at a 50% higher rate than those inheriting B haplotypes, primarily owing to the activity of a C>A-dominated mutational signature known as SBS18. The B and D quantitative trait locus haplotypes encode different alleles of Mutyh, a DNA repair gene that underlies the heritable cancer predisposition syndrome that causes colorectal tumors with a high SBS18 mutation load3,4. Both B and D Mutyh alleles are present in wild populations of Mus musculus domesticus, providing evidence that common genetic variation modulates germline mutagenesis in a model mammalian species.

GenomeMUSter mouse genetic variation service enables multitrait, multipopulation data integration and analysis.

Hundreds of inbred mouse strains and intercross populations have been used to characterize the function of genetic variants that contribute to disease. Thousands of disease-relevant traits have been characterized in mice and made publicly available. New strains and populations including consomics, the collaborative cross, expanded BXD, and inbred wild-derived strains add to existing complex disease mouse models, mapping populations, and sensitized backgrounds for engineered mutations. The genome sequences of inbred strains, along with dense genotypes from others, enable integrated analysis of trait-variant associations across populations, but these analyses are hampered by the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense variant resource by harmonizing multiple data sets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extendable to other model organisms. The result is a web- and programmatically accessible data service called GenomeMUSter, comprising single-nucleotide variants covering 657 strains at 106.8 million segregating sites. Interoperation with phenotype databases, analytic tools, and other resources enable a wealth of applications, including multitrait, multipopulation meta-analysis. We show this in cross-species comparisons of type 2 diabetes and substance use disorder meta-analyses, leveraging mouse data to characterize the likely role of human variant effects in disease. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

Systems genetics of metabolism: the use of the BXD murine reference panel for multiscalar integration of traits.

Metabolic homeostasis is achieved by complex molecular and cellular networks that differ significantly among individuals and are difficult to model with genetically engineered lines of mice optimized to study single gene function. Here, we systematically acquired metabolic phenotypes by using the EUMODIC EMPReSS protocols across a large panel of isogenic but diverse strains of mice (BXD type) to study the genetic control of metabolism. We generated and analyzed 140 classical phenotypes and deposited these in an open-access web service for systems genetics (www.genenetwork.org). Heritability, influence of sex, and genetic modifiers of traits were examined singly and jointly by using quantitative-trait locus (QTL) and expression QTL-mapping methods. Traits and networks were linked to loci encompassing both known variants and novel candidate genes, including alkaline phosphatase (ALPL), here linked to hypophosphatasia. The assembled and curated phenotypes provide key resources and exemplars that can be used to dissect complex metabolic traits and disorders.

A novel quantitative trait locus implicates Msh3 in the propensity for genome-wide short tandem repeat expansions in mice.

Short tandem repeats (STRs) are a class of rapidly mutating genetic elements typically characterized by repeated units of 1-6 bp. We leveraged whole-genome sequencing data for 152 recombinant inbred (RI) strains from the BXD family of mice to map loci that modulate genome-wide patterns of new mutations arising during parent-to-offspring transmission at STRs. We defined quantitative phenotypes describing the numbers and types of germline STR mutations in each strain and performed quantitative trait locus (QTL) analyses for each of these phenotypes. We identified a locus on Chromosome 13 at which strains inheriting the C57BL/6J (B) haplotype have a higher rate of STR expansions than those inheriting the DBA/2J (D) haplotype. The strongest candidate gene in this locus is Msh3, a known modifier of STR stability in cancer and at pathogenic repeat expansions in mice and humans, as well as a current drug target against Huntington's disease. The D haplotype at this locus harbors a cluster of variants near the 5' end of Msh3, including multiple missense variants near the DNA mismatch recognition domain. In contrast, the B haplotype contains a unique retrotransposon insertion. The rate of expansion covaries positively with Msh3 expression-with higher expression from the B haplotype. Finally, detailed analysis of mutation patterns showed that strains carrying the B allele have higher expansion rates, but slightly lower overall total mutation rates, compared with those with the D allele, particularly at tetranucleotide repeats. Our results suggest an important role for inherited variants in Msh3 in modulating genome-wide patterns of germline mutations at STRs.

Merged magnetic resonance and light sheet microscopy of the whole mouse brain.

We have developed workflows to align 3D magnetic resonance histology (MRH) of the mouse brain with light sheet microscopy (LSM) and 3D delineations of the same specimen. We start with MRH of the brain in the skull with gradient echo and diffusion tensor imaging (DTI) at 15 µm isotropic resolution which is ~ 1,000 times higher than that of most preclinical MRI. Connectomes are generated with superresolution tract density images of ~5 µm. Brains are cleared, stained for selected proteins, and imaged by LSM at 1.8 µm/pixel. LSM data are registered into the reference MRH space with labels derived from the ABA common coordinate framework. The result is a high-dimensional integrated volume with registration (HiDiver) with alignment precision better than 50 µm. Throughput is sufficiently high that HiDiver is being used in quantitative studies of the impact of gene variants and aging on mouse brain cytoarchitecture and connectomics.

Metabolic networks of longevity.

Molecular and cellular networks implicated in aging depend on a multitude of proteins that collectively mount adaptive and contingent metabolic responses to environmental challenges. Here, we discuss the intimate links between metabolic regulation and longevity and outline new approaches for analyzing and manipulating such links to promote human health span.

Systems genetics in the rat HXB/BXH family identifies Tti2 as a pleiotropic quantitative trait gene for adult hippocampal neurogenesis and serum glucose.

Neurogenesis in the adult hippocampus contributes to learning and memory in the healthy brain but is dysregulated in metabolic and neurodegenerative diseases. The molecular relationships between neural stem cell activity, adult neurogenesis, and global metabolism are largely unknown. Here we applied unbiased systems genetics methods to quantify genetic covariation among adult neurogenesis and metabolic phenotypes in peripheral tissues of a genetically diverse family of rat strains, derived from a cross between the spontaneously hypertensive (SHR/OlaIpcv) strain and Brown Norway (BN-Lx/Cub). The HXB/BXH family is a very well established model to dissect genetic variants that modulate metabolic and cardiovascular diseases and we have accumulated deep phenome and transcriptome data in a FAIR-compliant resource for systematic and integrative analyses. Here we measured rates of precursor cell proliferation, survival of new neurons, and gene expression in the hippocampus of the entire HXB/BXH family, including both parents. These data were combined with published metabolic phenotypes to detect a neurometabolic quantitative trait locus (QTL) for serum glucose and neuronal survival on Chromosome 16: 62.1-66.3 Mb. We subsequently fine-mapped the key phenotype to a locus that includes the Telo2-interacting protein 2 gene (Tti2)-a chaperone that modulates the activity and stability of PIKK kinases. To verify the hypothesis that differences in neurogenesis and glucose levels are caused by a polymorphism in Tti2, we generated a targeted frameshift mutation on the SHR/OlaIpcv background. Heterozygous SHR-Tti2+/- mutants had lower rates of hippocampal neurogenesis and hallmarks of dysglycemia compared to wild-type littermates. Our findings highlight Tti2 as a causal genetic link between glucose metabolism and structural brain plasticity. In humans, more than 800 genomic variants are linked to TTI2 expression, seven of which have associations to protein and blood stem cell factor concentrations, blood pressure and frontotemporal dementia.

Facilitating Complex Trait Analysis via Reduced Complexity Crosses.

Genetically diverse inbred strains are frequently used in quantitative trait mapping to identify sequence variants underlying trait variation. Poor locus resolution and high genetic complexity impede variant discovery. As a solution, we explore reduced complexity crosses (RCCs) between phenotypically divergent, yet genetically similar, rodent substrains. RCCs accelerate functional variant discovery via decreasing the number of segregating variants by orders of magnitude. The simplified genetic architecture of RCCs often permit immediate identification of causal variants or rapid fine-mapping of broad loci to smaller intervals. Whole-genome sequences of substrains make RCCs possible by supporting the development of array- and targeted sequencing-based genotyping platforms, coupled with rapid genome editing for variant validation. In summary, RCCs enhance discovery-based genetics of complex traits.

Increased development of T-bet+CD11c+ B cells predisposes to lupus in females: Analysis in BXD2 mouse and genetic crosses.

Cardinal features of lupus include elevated B cell activation and autoantibody production with a female sex preponderance. We quantified interactions of sex and genetic variation on the development of autoimmune B-cell phenotypes and autoantibodies in the BXD2 murine model of lupus using a cohort of backcrossed progeny (BXD2 x C57BL/6J) x BXD2. Sex was the key factor leading to increased total IgG, IgG2b, and autoantibodies. The percentage of T-bet+CD11c+ IgD+ activated naive B cells (aNAV) was higher in females and was associated with increased T-bet+CD11c+ IgD- age-related B cells, Fas+GL7+ germinal center B cells, Cxcr5-Icos+ peripheral T-helper cells, and Cxcr5+Icos+ follicular T-helper cells. IFN-ß was elevated in females. Variation in aNAV cells was mapped to Chr 7 in a locus that shows significant interactions between the female sex and heterozygous B/D variant. Our results suggest that activation of naive B cells forms the basis for the female-predominant development of autoantibodies in lupus-susceptible BXD2 mice.

Current Advancements in Mouse Models of Retinal Disease.

The field of retinal degenerative (RDs) disease study has been in a state of exponential growth from discovering the underlying genetic components of such diseases as age-related macular degeneration (AMD) and retinitis pigmentosa (RP) to the first gene therapy developed and approved for human Leber congenital amaurosis. However, a source for high-fidelity animal models of these complex, multifactorial, and/or polygenic diseases is a need that has yet to be fulfilled. While models for AMD and RP do exist, they often require aging the animals for a year or more, feeding special diets, or introduction of external modulators such as exposure to cigarette smoke. Currently, work is being done to uncover high-fidelity naturally occurring models of these retinal diseases with the hope and intent of providing the vision community the tools it needs to better understand, treat, and, one day, cure the patients suffering from these devastating afflictions.

Highlights from the Era of Open Source Web-Based Tools.

High digital connectivity and a focus on reproducibility are contributing to an open science revolution in neuroscience. Repositories and platforms have emerged across the whole spectrum of subdisciplines, paving the way for a paradigm shift in the way we share, analyze, and reuse vast amounts of data collected across many laboratories. Here, we describe how open access web-based tools are changing the landscape and culture of neuroscience, highlighting six free resources that span subdisciplines from behavior to whole-brain mapping, circuits, neurons, and gene variants.

Identification of cyclin D1 as a major modulator of 3-nitropropionic acid-induced striatal neurodegeneration.

Mitochondria dysfunction occurs in the aging brain as well as in several neurodegenerative disorders and predisposes neuronal cells to enhanced sensitivity to neurotoxins. 3-nitropropionic acid (3-NP) is a naturally occurring plant and fungal neurotoxin that causes neurodegeneration predominantly in the striatum by irreversibly inhibiting the tricarboxylic acid respiratory chain enzyme, succinate dehydrogenase (SDH), the main constituent of the mitochondria respiratory chain complex II. Significantly, although 3-NP-induced inhibition of SDH occurs in all brain regions, neurodegeneration occurs primarily and almost exclusively in the striatum for reasons still not understood. In rodents, 3-NP-induced striatal neurodegeneration depends on the strain background suggesting that genetic differences among genotypes modulate toxicant variability and mechanisms that underlie 3-NP-induced neuronal cell death. Using the large BXD family of recombinant inbred (RI) strains we demonstrate that variants in Ccnd1 - the gene encoding cyclin D1 - of the DBA/2 J parent underlie the resistance to 3-NP-induced striatal neurodegeneration. In contrast, the Ccnd1 variant inherited from the widely used C57BL/6 J parental strain confers sensitivity. Given that cellular stress triggers induction of cyclin D1 expression followed by cell-cycle re-entry and consequent neuronal cell death, we sought to determine if the C57BL/6 J and DBA/2 J Ccnd1 variants are differentially modulated in response to 3-NP. We confirm that 3-NP induces cyclin D1 expression in striatal neuronal cells of C57BL/6 J, but this response is blunted in the DBA/2 J. We further show that striatal-specific alternative processing of a highly conserved 3'UTR negative regulatory region of Ccnd1 co-segregates with the C57BL/6 J parental Ccnd1 allele in BXD strains and that its differential processing accounts for sensitivity or resistance to 3-NP. Our results indicate that naturally occurring Ccnd1 variants may play a role in the variability observed in neurodegenerative disorders involving mitochondria complex II dysfunction and point to cyclin D1 as a possible therapeutic target.

Identifying gene function and module connections by the integration of multispecies expression compendia.

The functions of many eukaryotic genes are still poorly understood. Here, we developed and validated a new method, termed GeneBridge, which is based on two linked approaches to impute gene function and bridge genes with biological processes. First, Gene-Module Association Determination (G-MAD) allows the annotation of gene function. Second, Module-Module Association Determination (M-MAD) allows predicting connectivity among modules. We applied the GeneBridge tools to large-scale multispecies expression compendia-1700 data sets with over 300,000 samples from human, mouse, rat, fly, worm, and yeast-collected in this study. G-MAD identifies novel functions of genes-for example, DDT in mitochondrial respiration and WDFY4 in T cell activation-and also suggests novel components for modules, such as for cholesterol biosynthesis. By applying G-MAD on data sets from respective tissues, tissue-specific functions of genes were identified-for instance, the roles of EHHADH in liver and kidney, as well as SLC6A1 in brain and liver. Using M-MAD, we identified a list of module-module associations, such as those between mitochondria and proteasome, mitochondria and histone demethylation, as well as ribosomes and lipid biosynthesis. The GeneBridge tools together with the expression compendia are available as an open resource, which will facilitate the identification of connections linking genes, modules, phenotypes, and diseases.

Gabra2 is a genetic modifier of Dravet syndrome in mice.

Pathogenic variants in epilepsy genes result in a spectrum of clinical severity. One source of phenotypic heterogeneity is modifier genes that affect expressivity of a primary pathogenic variant. Mouse epilepsy models also display varying degrees of clinical severity on different genetic backgrounds. Mice with heterozygous deletion of Scn1a (Scn1a+/-) model Dravet syndrome, a severe epilepsy most often caused by SCN1A haploinsufficiency. Scn1a+/- mice recapitulate features of Dravet syndrome, including spontaneous seizures, sudden death, and cognitive/behavioral deficits. Scn1a+/- mice maintained on the 129S6/SvEvTac (129) strain have normal lifespan and no spontaneous seizures. In contrast, admixture with C57BL/6J (B6) results in epilepsy and premature lethality. We previously mapped Dravet Survival Modifier loci (Dsm1-Dsm5) responsible for strain-dependent differences in survival. Gabra2, encoding the GABAA α2 subunit, was nominated as a candidate modifier at Dsm1. Direct measurement of GABAA receptors found lower abundance of α2-containing receptors in hippocampal synapses of B6 mice relative to 129. We also identified a B6-specific single nucleotide deletion within Gabra2 that lowers mRNA and protein by nearly 50%. Repair of this deletion reestablished normal levels of Gabra2 expression. In this study, we used B6 mice with a repaired Gabra2 allele to evaluate Gabra2 as a genetic modifier of severity in Scn1a+/- mice. Gabra2 repair restored transcript and protein expression, increased abundance of α2-containing GABAA receptors in hippocampal synapses, and rescued epilepsy phenotypes of Scn1a+/- mice. These findings validate Gabra2 as a genetic modifier of Dravet syndrome, and support enhancing function of α2-containing GABAA receptors as treatment strategy for Dravet syndrome.

Genetic differences in ethanol consumption: effects on iron, copper, and zinc regulation in mouse hippocampus.

One common characteristic of neurodegenerative diseases is dysregulation of iron, usually with observed increases in its concentration in various regions. Heavy alcohol consumption is believed to contribute to such iron dysregulation in the brain with accompanying dementia. To examine this effect and related genetic-based individual differences in an animal model, we subjected female mice from 12 BXD recombinant inbred strains to 16 weeks of alcohol consumption using the drinking in the dark (DID) method. Daily consumption was recorded and at the end of 16 weeks hippocampus tissues harvested. Concentrations of iron, copper and zinc were measured using X-ray fluorescence technology. The results showed that, DID increased iron overall across all strains, ranging from 3 to 68%. Copper and Zinc both decreased, ranging from 0.4-42 and 5-35% respectively. Analysis of variance revealed significant strain by treatment interactions for all three metals. Additionally, in the DID group, we observed strain differences in reduction of hippocampus mass. These findings are particularly interesting to us because high alcohol consumption in humans has been associated with neurodegeneration and dementia related to disruption of iron regulation. The findings of alcohol consumption associated decreases in copper and zinc are novel. The role of copper regulation and neurological function related to alcohol consumption is as yet largely unexplored. The role of zinc is better known as a neuromodulator in the hippocampus and appears to be protective against neurological damage. It would seem then, that the alcohol-related decrease in zinc in the hippocampus would be of concern and warrants further study.

Identifying the molecular systems that influence cognitive resilience to Alzheimer's disease in genetically diverse mice.

Individual differences in cognitive decline during normal aging and Alzheimer's disease (AD) are common, but the molecular mechanisms underlying these distinct outcomes are not fully understood. We utilized a combination of genetic, molecular, and behavioral data from a mouse population designed to model human variation in cognitive outcomes to search for the molecular mechanisms behind this population-wide variation. Specifically, we used a systems genetics approach to relate gene expression to cognitive outcomes during AD and normal aging. Statistical causal-inference Bayesian modeling was used to model systematic genetic perturbations matched with cognitive data that identified astrocyte and microglia molecular networks as drivers of cognitive resilience to AD. Using genetic mapping, we identified Fgf2 as a potential regulator of the astrocyte network associated with individual differences in short-term memory. We also identified several immune genes as regulators of a microglia network associated with individual differences in long-term memory, which was partly mediated by amyloid burden. Finally, significant overlap between mouse and two different human coexpression networks provided strong evidence of translational relevance for the genetically diverse AD-BXD panel as a model of late-onset AD. Together, this work identified two candidate molecular pathways enriched for microglia and astrocyte genes that serve as causal AD cognitive biomarkers, and provided a greater understanding of processes that modulate individual and population-wide differences in cognitive outcomes during AD.

Identifying modifier genes for hypertrophic cardiomyopathy.

BACKGROUND: Hypertrophic cardiomyopathy (HCM) severity greatly varies among patients even with the same HCM gene mutations. This variation is largely regulated by modifier gene(s), which, however, remain largely unknown. The current study is aimed to identify modifier genes using BXD strains, a large murine genetic reference population (GRP) derived from crosses between C57BL/6 J (B6) and D2 DBA/2 J (D2) mice. D2 mice natualy carrythe genetic basis and phenotypes of HCM. METHODS: Myocardial hypertrophy, the major phenotype of HCM, was determined by cardiomyocyte size on cardiac sections in 30 BXD strains, and their parental B6 and D2 strains and morphometric analysis was performed. Quantitative Trait Locus (QTL) mapping for cardiomyocyte sizes was conducted with WebQTL in GeneNetwork. Correlation of cardiomyocyte size and cardiac gene expression in BXDs accessed from GeneNetwork were evaluated. QTL candidate genes associated with cardiomyocyte sizes were prioritized based on the score system. RESULTS: Cardiomyocyte size varied significantly among BXD strains. Interval mapping on cardiomyocyte size data showed a significant QTL on chromosome (Chr) 2 at 66- 73.5 Mb and a suggestive QTL on Chr 5 at 20.9-39.7 Mb. Further score system revealed a high QTL score for Xirp2 in Chr 2. Xirp2 encodes xin actin-binding repeat containing 2, which is highly expressed in cardiac tissue and associate with cardiomyopathy and heart failure. In Chr5 QTL, Nos3, encoding nitric oxide synthase 3, received the highest score, which is significantly correlated with cardiomyocyte size. CONCLUSION: These results indicate that Xirp2 and Nos3 serve as novel candidate modifier genes for myocardial hypertrophy in HCM. These candidate genes will be validated in our future studies.

Variability and heritability of mouse brain structure: Microscopic MRI atlases and connectomes for diverse strains.

Genome-wide association studies have demonstrated significant links between human brain structure and common DNA variants. Similar studies with rodents have been challenging because of smaller brain volumes. Using high field MRI (9.4 T) and compressed sensing, we have achieved microscopic resolution and sufficiently high throughput for rodent population studies. We generated whole brain structural MRI and diffusion connectomes for four diverse isogenic lines of mice (C57BL/6J, DBA/2J, CAST/EiJ, and BTBR) at spatial resolution 20,000 times higher than human connectomes. We measured narrow sense heritability (h2) I.e. the fraction of variance explained by strains in a simple ANOVA model for volumes and scalar diffusion metrics, and estimates of residual technical error for 166 regions in each hemisphere and connectivity between the regions. Volumes of discrete brain regions had the highest mean heritability (0.71 ± 0.23 SD, n = 332), followed by fractional anisotropy (0.54 ± 0.26), radial diffusivity (0.34 ± 0.022), and axial diffusivity (0.28 ± 0.19). Connection profiles were statistically different in 280 of 322 nodes across all four strains. Nearly 150 of the connection profiles were statistically different between the C57BL/6J, DBA/2J, and CAST/EiJ lines. Microscopic whole brain MRI/DTI has allowed us to identify significant heritable phenotypes in brain volume, scalar DTI metrics, and quantitative connectomes.

Genetic variability of T cell responses in hypersensitivity pneumonitis identified using the BXD genetic reference panel.

Hypersensitivity pneumonitis (HP) is an interstitial lung disease that may progress to fibrosis and significant risk of death. HP develops following repeated exposures to inhaled environmental antigens; however, only a fraction of the exposed population develops the disease, suggesting that host genetics contribute to disease susceptibility. We used the BXD family of mice with the Saccharopolyspora rectivirgula (SR) model of HP to investigate the role of genetics in susceptibility to HP. The BXD family is derived from a B6 mother and a D2 father and has been used to map susceptibility loci to numerous diseases. B6, D2, and BXD progeny strains were exposed to SR for 3 wk, and the development of HP was monitored. The B6 and D2 strains developed alveolitis; however, the cellular composition was neutrophilic in the D2 strain and more lymphocytic in the B6 strain. Hematoxylin-eosin staining of lung sections revealed lymphoid aggregates in B6 lungs, whereas D2 lungs exhibited a neutrophilic infiltration. Twenty-eight BXD strains of mice were tested, and the results reveal significant heritable variation for numbers of CD4+ or CD8+ T cells in the air spaces. There was significant genetic variability for lymphoid aggregates and alveolar wall thickening. We mapped a significant quantitative trait locus (QTL) on chromosome 18 for CD8+CD69+ T cells that includes cadherin 2 (Cdh2), an excellent candidate gene associated with epithelial-mesenchymal transition, which is upregulated in lungs of strains with HP. These results demonstrate that the BXD family is a valuable and translationally relevant model to identify genes contributing to HP and to devise early and effective interventions.

Interleukin 17-producing T helper cells and interleukin 17 orchestrate autoreactive germinal center development in autoimmune BXD2 mice.

Interleukin 17 (IL-17) is a cytokine associated with inflammation, autoimmunity and defense against some bacteria. Here we show that IL-17 can promote autoimmune disease through a mechanism distinct from its proinflammatory effects. As compared with wild-type mice, autoimmune BXD2 mice express more IL-17 and show spontaneous development of germinal centers (GCs) before they increase production of pathogenic autoantibodies. We show that blocking IL-17 signaling disrupts CD4+ T cell and B cell interactions required for the formation of GCs and that mice lacking the IL-17 receptor have reduced GC B cell development and humoral responses. Production of IL-17 correlates with upregulated expression of the genes Rgs13 and Rgs16, which encode regulators of G-protein signaling, and results in suppression of the B cell chemotactic response to the chemokine CXCL12. These findings suggest a mechanism by which IL-17 drives autoimmune responses by promoting the formation of spontaneous GCs.