Subspecific origin and haplotype diversity in the laboratory mouse

Phylogenomic analyses revealed a new lineage of house mouse (Mus musculus) in Gyirong Basin of Xizang Autonomous Region, China

In the present study, we collected 20 individuals and 12 individuals of wild mice from the Gyirong Basin of Xizang Autonomous Region in China and Sudurpashchim in Nepal. Phylogeny and genetic structure inferred from different types of genomic markers suggest that these samples all belong to Mus musculus, among which individuals from Gyirong Basin represent a new genomic lineage (named as M. m. gyirongus), and samples from Sudurpashchim represent an intermediate population between the central population and M. m. castaneus. M. m. gyirongus, along with M. m. domesticus and M. m. musculus, differentiated from the central population compactly during ∼ 272,000-251,000 years ago in the interglacial period. Three lineages all experienced continuous population decline before ∼ 70,000 years ago. Then, they underwent population fluctuations at different periods that might have been impacted by climate changes, migration history, and human activities. Genes related to the structure and function of neural synapses, reproduction and development, regulation of cell cycle and carcinogenesis, and immune response have undergone positive selection in the genome of M. m. gyirongus. The discovery of M. m. gyirongus not only helps us to better understand the evolutionary history of M. musculus, but also provides new regional resources for breeding novel laboratory mouse strains.

Intact, recombinant, and spliced forms of endogenous mouse mammary tumor viruses in inbred and wild mice

Endogenous retroviruses (ERVs) are chromosomally integrated viral copies that represent relics of past infections. Analysis of the sequenced genomes of 17 mouse strains, Mus musculus subspecies, and Mus spretus identified 29 ERVs of mouse mammary tumor viruses (MMTVs), termed Mtvs. The 15 laboratory mouse Mtvs are each present in multiple strains reflecting their common breeding history; most predate the development of inbred strains and were likely acquired by Mus musculus domesticus progenitors but have no orthologs in wild mice, whereas four, including the intact Mtv1, were likely endogenized more recently. One of the 14 Mtvs found in wild mice was distributed over a broad geographic range in southeast Asia. Most Mtvs are full-length, with multiple open reading frames, but Mtvs from many wild mice have an unusual envelope deletion corresponding to an intron of the viral rem accessory gene, suggesting its derivation from spliced MMTV cDNAs. These deleted envs have open reading frames, are found in globally distributed mice, and show subspecies-specific sequence variation consistent with their recurrent generation. The highly variable MMTV sag gene, responsible for resistance to exogenous infection, exhibits evidence of recombination as well as positive selection, consistent with its role in antiviral defense. In contrast, the spread of Mtvs in Mus musculus populations is not marked by an active arms race pitting the MMTV envelope against its cellular receptor. Thus, the acquisition of potentially disease-inducing Mtvs is a recent and ongoing process in Mus accompanied by recombination, positive selection, and a recurrent envelope deletion.

IMPORTANCE

Endogenous retroviruses (ERVs) are copies of viral genomes inserted into host chromosomes, producing a fossil record of past infections and virus-host co-adaptations. ERVs of mouse mammary tumor viruses (Mtvs) were found in all common laboratory strains, all Mus musculus subspecies, and a sister species, Mus spretus. Most laboratory mouse Mtvs predate inbred strain origins and were acquired by M. musculus domesticus, but although widely shared among strains, none of these were found in wild mice. Among wild mouse Mtvs, only one showed a broad geographic distribution. All M. musculus subspecies carry Mtvs with a large envelope deletion corresponding to the processed mRNA for the viral rem gene; such Mtvs likely derive from spliced viral mRNA. The Mtv sag gene responsible for resistance to exogenous infection is under purifying selection and has been subject to recombination, whereas the Mtv envelope and its cellular receptor show no evidence of genetic conflicts.

Phenogenomic resources immortalized in a panel of wild-derived strains of five species of house mice

The house mouse, Mus musculus, is a widely used animal model in biomedical research, with classical laboratory strains (CLS) being the most frequently employed. However, the limited genetic variability in CLS hinders their applicability in evolutionary studies. Wild-derived strains (WDS), on the other hand, provide a suitable resource for such investigations. This study quantifies genetic and phenotypic data of 101 WDS representing 5 species, 3 subspecies, and 8 natural Y consomic strains and compares them with CLS. Genetic variability was estimated using whole mtDNA sequences, the Prdm9 gene, and copy number variation at two sex chromosome-linked genes. WDS exhibit a large natural variation with up to 2173 polymorphic sites in mitogenomes, whereas CLS display 92 sites. Moreover, while CLS have two Prdm9 alleles, WDS harbour 46 different alleles. Although CLS resemble M. m. domesticus and M. m. musculus WDS, they differ from them in 10 and 14 out of 16 phenotypic traits, respectively. The results suggest that WDS can be a useful tool in evolutionary and biomedical studies with great potential for medical applications.

Gene-by-environment Interactions and Adaptive Body Size Variation in Mice From the Americas

The relationship between genotype and phenotype is often mediated by the environment. Moreover, gene-by-environment (GxE) interactions can contribute to variation in phenotypes and fitness. In the last 500 yr, house mice have invaded the Americas. Despite their short residence time, there is evidence of rapid climate adaptation, including shifts in body size and aspects of metabolism with latitude. Previous selection scans have identified candidate genes for metabolic adaptation. However, environmental variation in diet as well as GxE interactions likely impact body mass variation in wild populations. Here, we investigated the role of the environment and GxE interactions in shaping adaptive phenotypic variation. Using new locally adapted inbred strains from North and South America, we evaluated response to a high-fat diet, finding that sex, strain, diet, and the interaction between strain and diet contributed significantly to variation in body size. We also found that the transcriptional response to diet is largely strain-specific, indicating that GxE interactions affecting gene expression are pervasive. Next, we used crosses between strains from contrasting climates to characterize gene expression regulatory divergence on a standard diet and on a high-fat diet. We found that gene regulatory divergence is often condition-specific, particularly for trans-acting changes. Finally, we found evidence for lineage-specific selection on cis-regulatory variation involved in diverse processes, including lipid metabolism. Overlap with scans for selection identified candidate genes for environmental adaptation with diet-specific effects. Together, our results underscore the importance of environmental variation and GxE interactions in shaping adaptive variation in complex traits.

Machine vision-based frailty assessment for genetically diverse mice

Frailty indexes (FIs) capture health status in humans and model organisms. To accelerate our understanding of biological aging and carry out scalable interventional studies, high-throughput approaches are necessary. We previously introduced a machine vision-based visual frailty index (vFI) that uses mouse behavior in the open field to assess frailty using C57BL/6J (B6J) data. Aging trajectories are highly genetic and are frequently modeled in genetically diverse animals. In order to extend the vFI to genetically diverse mouse populations, we collect frailty and behavior data on a large cohort of aged Diversity Outbred (DO) mice. Combined with previous data, this represents one of the largest video-based aging behavior datasets to date. Using these data, we build accurate predictive models of frailty, chronological age, and even the proportion of life lived. The extension of automated and objective frailty assessment tools to genetically diverse mice will enable better modeling of aging mechanisms and enable high-throughput interventional aging studies.

Recent genetic drift in the co-diversified gut bacterial symbionts of laboratory mice

Laboratory mice (Mus musculus domesticus) harbor gut bacterial strains that are distinct from those of wild mice but whose evolutionary histories are unclear. Here, we show that laboratory mice have retained gut bacterial lineages that diversified in parallel (co-diversified) with rodent species for > 25 million years, but that laboratory-mouse gut microbiota (LGM) strains of these ancestral symbionts have experienced accelerated accumulation of genetic load during the past ~ 120 years of captivity. Compared to closely related wild-mouse gut microbiota (WGM) strains, co-diversified LGM strains displayed significantly faster genome-wide rates of nonsynonymous substitutions, indicating elevated genetic drift-a difference that was absent in non-co-diversified symbiont clades. Competition experiments in germ-free mice further indicated that LGM strains within co-diversified clades displayed significantly reduced fitness in vivo compared to WGM relatives to an extent not observed within non-co-diversified clades. Thus, stochastic processes (e.g., bottlenecks), not natural selection in the laboratory, have been the predominant evolutionary forces underlying divergence of co-diversified symbiont strains between laboratory and wild house mice. Our results show that gut bacterial lineages conserved in diverse rodent species have acquired novel mutational burdens in laboratory mice, providing an evolutionary rationale for restoring laboratory mice with wild gut bacterial strain diversity.

Genetic architecture of the red blood cell proteome in genetically diverse mice reveals central role of hemoglobin beta cysteine redox status in maintaining circulating glutathione pools

Red blood cells (RBCs) transport oxygen but accumulate oxidative damage over time, reducing function in vivo and during storage-critical for transfusions. To explore genetic influences on RBC resilience, we profiled proteins, metabolites, and lipids from fresh and stored RBCs obtained from 350 genetically diverse mice. Our analysis identified over 6,000 quantitative trait loci (QTL). Compared to other tissues, prevalence of trans genetic effects over cis reflects the absence of de novo protein synthesis in anucleated RBCs. QTL hotspots at Hbb, Hba, Mon1a, and storage-specific Steap3 linked ferroptosis to hemolysis. Proteasome components clustered at multiple loci, underscoring the importance of degrading oxidized proteins. Post-translational modifications (PTMs) mapped predominantly to hemoglobins, particularly cysteine residues. Loss of reactive C93 in humanized mice (HBB C93A) disrupted redox balance, affecting glutathione pools, protein glutathionylation, and redox PTMs. These findings highlight genetic regulation of RBC oxidation, with implications for transfusion biology and oxidative stress-dependent hemolytic disorders.

Spatiotemporal Diffusion, Colonization, and Antibody Responses in Susceptible C57BL/6J Mice Orally Infected with Toxoplasma gondii Cysts

Toxoplasma gondii is an obligate intracellular protozoan that infects humans and other mammals. The C57BL/6J mouse strain is regarded as an ideal model organism for studying T. gondii due to its susceptibility to T. gondii infection and its other advantages over other laboratory animals. However, systematic studies on the response dynamics of the susceptible C57BL/6J mice after oral infection with T. gondii cysts are lacking. To address this research gap, we investigated the spatiotemporal dynamics of infection, colonization, and antibody fluctuations in susceptible C57BL/6J mice orally infected with Type II T. gondii ME49 strain cysts. Mice were orally challenged with T. gondii cysts to examine the infection dynamics. Daily monitoring was conducted for 60 days post-infection (dpi) to assess animals' clinical signs and survival rates. The parasite burden in various organs was quantified using qPCR targeting the T. gondii B1 gene. The serum antibody responses were evaluated using ELISA. The cyst burden in the mouse brain was assessed via histology and immunofluorescence. T. gondii infection induced clinical symptoms in the mice, including fever and weight loss. T. gondii rapidly invaded the mice's small intestine, spleen, lungs, liver, and heart via the bloodstream within 1-5 dpi. T. gondii had breached the blood-brain barrier and colonized the brain by 7 dpi. The levels of Toxoplasma-specific IgG antibodies increased and stabilized for two months (until the experiment ended). Systemic parasite dissemination occurred rapidly, infiltrating most tissues and organs, leading to pronounced enteritis and multi-organ damage due to inflammation. The tachyzoites differentiated into bradyzoites when T. gondii infection progressed from the acute to the chronic phase in mice, forming tissue cysts in organs, including the muscles and brain. As a result, the predilection site of T. gondii in mice is the brain, which is where the cysts persisted for the host's lifetime and continuously induced meningitis. These findings provide valuable insights into the spatiotemporal diffusion, colonization, predilection sites, temporal antibody dynamics, pathogen detection methodologies, and histopathological changes in C57BL/6J mice following oral infection with T. gondii cysts. These insights are important for elucidating T. gondii's pathogenesis and host-T. gondii interaction.

Immunological roads diverged: mapping tuberculosis outcomes in mice

The journey from phenotypic observation to causal genetic mechanism is a long and challenging road. For pathogens like Mycobacterium tuberculosis (Mtb), which causes tuberculosis (TB), host-pathogen coevolution has spanned millennia, costing millions of human lives. Mammalian models can systematically recapitulate host genetic variation, producing a spectrum of disease outcomes. Leveraging genome sequences and deep phenotyping data from infected mouse genetic reference populations (GRPs), quantitative trait locus (QTL) mapping approaches have successfully identified host genomic regions associated with TB phenotypes. Here, we review the ongoing optimization of QTL mapping study design alongside advances in mouse GRPs. These next-generation resources and approaches have enabled identification of novel host-pathogen interactions governing one of the most prevalent infectious diseases in the world today.

Protocol for murine multi-tissue deep immunophenotyping using a 40-color full-spectrum flow cytometry panel

The innate and adaptive immune systems, though often studied separately, interact deeply and respond to stimuli simultaneously, with leukocytes displaying a range of pro- to anti-inflammatory phenotypes. This protocol details a procedure for characterizing murine innate and adaptive immune phenotypes using a 40-color full-spectral flow cytometry panel. We describe steps for organ collection, sample preparation, immunofluorescent staining, and acquisition to reproducibly and cost-effectively study tissue-resident leukocytes, their subpopulations, and inflammatory status in various organs.

Anapc5 and Anapc7 as genetic modifiers of KIF18A function in fertility and mitotic progression

The kinesin family member 18A (KIF18A) is an essential regulator of microtubule dynamics and chromosome alignment during mitosis. Functional dependency on KIF18A varies by cell type and genetic context but the heritable factors that influence this dependency remain unknown. To address this, we took advantage of the variable penetrance observed in different mouse strain backgrounds to screen for loci that modulate germ cell depletion in the absence of KIF18A. We found a significant association at a Chr5 locus where anaphase promoting complex subunits 5 (Anapc5) and 7 (Anapc7) were the top candidate genes. We found that both genes were differentially expressed in a sensitive strain background when compared to resistant strain background at key timepoints in gonadal development. We also identified a novel retroviral insertion in Anapc7 that may in part explain the observed expression differences. In cell line models, we found that depletion of KIF18A induced mitotic arrest, which was partially rescued by co-depletion of ANAPC7 (APC7) and exacerbated by co-depletion of ANAPC5 (APC5). These findings suggest that differential expression and activity of Anapc5 and Anapc7 may influence sensitivity to KIF18A depletion in germ cells and CIN cells, with potential implications for optimizing antineoplastic therapies.

Stem cell transcriptional profiles from mouse subspecies reveal cis-regulatory evolution at translation genes

A key goal of evolutionary genomics is to harness molecular data to draw inferences about selective forces that have acted on genomes. The field progresses in large part through the development of advanced molecular-evolution analysis methods. Here we explored the intersection between classical sequence-based tests for selection and an empirical expression-based approach, using stem cells from Mus musculus subspecies as a model. Using a test of directional, cis-regulatory evolution across genes in pathways, we discovered a unique program of induction of translation genes in stem cells of the Southeast Asian mouse M. m. castaneus relative to its sister taxa. We then mined population-genomic sequences to pursue underlying regulatory mechanisms for this expression divergence, finding robust evidence for alleles unique to M. m. castaneus at the upstream regions of the translation genes. We interpret our data under a model of changes in lineage-specific pressures across Mus musculus in stem cells with high translational capacity. Our findings underscore the rigor of integrating expression and sequence-based methods to generate hypotheses about evolutionary events from long ago.

Voluntary nicotine consumption and reward in a subset of diversity outbred founder strains

BACKGROUND

Nicotine is largely responsible for the initiation and maintenance of tobacco dependence and contributes to a global health problem.

AIMS

This study characterizes nicotine oral consumption and preference in male and female mice of several Diversity Outbred (DO) founder strains: C57BL/6J, A/J, 129S1/SvImJ, PWK/PhJ, NOD/ShiLtJ, and CAST/EiJ. It assesses the impact of nicotine concentration on intake and preference, the potential interaction of strain with sex, and estimates the degree of heritable variation in nicotine consumption.

METHODS

Two-bottle choice oral self-administration paradigm was used to assess nicotine intake, nicotine preference, and total fluid intake in male and female mice of each strain in a concentration-response manner. A conditioned place preference (CPP) test was performed to evaluate the rewarding and aversive effects of nicotine in certain strains after systemic administration of the drug.

RESULTS

The highest nicotine-consuming strain was found to be 129S1/SvlmJ, and the lowest nicotine-consuming strain was A/J. Strain differences in nicotine intake were not due to differences in bitter and sweet tastes as shown in the saccharine and quinine two-bottle choice tests. A/J strain showed no significant CPP for nicotine while the 129S1/SvImJ strain showed a significant CPP for nicotine and a higher preference when compared to the C57BL/6J strain. Heritability estimates of nicotine intake were sex dependent and concentration dependent.

CONCLUSIONS

Data support that nicotine consumption patterns are heritable with an influence of genotype in a voluntary oral self-administration paradigm. Results pave the way for future studies with the highly recombinant DO mice that might lead to the identification of novel genetic loci and genes influencing nicotine consumption.

Mapping DNA Methylation to Cardiac Pathologies Induced by Beta-Adrenergic Stimulation in a Large Panel of Mice

Background

Heart failure (HF) is a leading cause of morbidity and mortality worldwide, with over 18 million deaths annually. Despite extensive research, genetic and environmental factors contributing to HF remain complex and poorly understood. Recent studies suggest that epigenetic modifications, such as DNA methylation, may play a crucial role in regulating HF-associated phenotypes. In this study, we leverage the Hybrid Mouse Diversity Panel (HMDP), a cohort of over 100 inbred mouse strains, to investigate the role of DNA methylation in HF progression.

Objective

We aim to identify epigenetic modifications associated with HF by integrating DNA methylation data with gene expression and phenotypic traits. Using isoproterenol (ISO)-induced cardiac hypertrophy and failure in HMDP mice, we explore the relationship between methylation patterns and HF susceptibility.

Methods

We performed reduced representational bisulfite sequencing (RRBS) to capture DNA methylation at single-nucleotide resolution in the left ventricles of 90 HMDP mouse strains under both control and ISO-treated conditions. We identified differentially methylated regions (DMRs) and performed an epigenome-wide association study (EWAS) using the MACAU algorithm. We identified likely candidate genes within each locus through integration of our results with previously reported sequence variation, gene expression, and HF-related phenotypes. In vitro approaches were employed to validate key findings, including gene knockdown experiments in neonatal rat ventricular myocytes (NRVMs). We also examined the effects of preventing DNA methyltransferase activity on HF progression.

Results

Our EWAS identified 56 CpG loci significantly associated with HF phenotypes, including 18 loci where baseline DNA methylation predicted post-ISO HF progression. Key candidate genes, such as Prkag2, Anks1, and Mospd3, were identified based on their epigenetic regulation and association with HF traits. In vitro follow-up on a number of genes confirmed that knockdown of Anks1 and Mospd3 in NRVMs resulted in significant alterations in cell size and blunting of ISO-induced hypertrophy, demonstrating their functional relevance in HF pathology.Furthermore, treatment with the DNA methyltransferase inhibitor RG108 in ISO-treated BTBRT mice significantly reduced cardiac hypertrophy and preserved ejection fraction compared to mice only treated with ISO, highlighting the therapeutic potential of targeting DNA methylation in HF. Differential expression analysis revealed that RG108 treatment restored the expression of several methylation-sensitive genes, further supporting the role of epigenetic regulation in HF.

Conclusion

Our study demonstrates a clear interplay between DNA methylation, gene expression, and HF-associated phenotypes. We identified several novel epigenetic loci and candidate genes that contribute to HF progression, offering new insights into the molecular mechanisms of HF. These findings underscore the importance of epigenetic regulation in cardiac disease and suggest potential therapeutic strategies for modifying HF outcomes through targeting DNA methylation.

Machine vision based frailty assessment for genetically diverse mice

Frailty indexes (FIs) capture health status in humans and model organisms. To accelerate our understanding of biological aging and carry out scalable interventional studies, high-throughput approaches are necessary. We previously introduced a machine vision-based visual frailty index (vFI) that uses mouse behavior in the open field to assess frailty using C57BL/6J (B6J) data. Aging trajectories are highly genetic and are frequently modeled in genetically diverse animals. In order to extend the vFI to genetically diverse mouse populations, we collect frailty and behavior data on a large cohort of aged Diversity Outbred (DO) mice. Combined with previous data, this represents one of the largest video-based aging behavior datasets to date. Using these data, we build accurate predictive models of frailty, chronological age, and even the proportion of life lived. The extension of automated and objective frailty assessment tools to genetically diverse mice will enable better modeling of aging mechanisms and enable high-throughput interventional aging studies.

The updated mouse universal genotyping array bioinformatic pipeline improves genetic QC in laboratory mice

The MiniMUGA genotyping array is a popular tool for genetic quality control of laboratory mice and genotyping samples from most experimental crosses involving laboratory strains, particularly for reduced complexity crosses. The content of the production version of the MiniMUGA array is fixed; however, there is the opportunity to improve the array's performance and the associated report's usefulness by leveraging thousands of samples genotyped since the initial description of MiniMUGA. Here, we report our efforts to update and improve marker annotation, increase the number and the reliability of the consensus genotypes for classical inbred strains and substrains, and increase the number of constructs reliably detected with MiniMUGA. In addition, we have implemented key changes in the informatics pipeline to identify and quantify the contribution of specific genetic backgrounds to the makeup of a given sample, remove arbitrary thresholds, include the Y Chromosome and mitochondrial genome in the ideogram, and improve robust detection of the presence of commercially available substrains based on diagnostic alleles. Finally, we have updated the layout of the report to simplify the interpretation and completeness of the analysis and added a section summarizing the ideogram in table format. These changes will be of general interest to the mouse research community and will be instrumental in our goal of improving the rigor and reproducibility of mouse-based biomedical research.

Hyperglycemia and cognitive impairments anticipate the onset of an overt type 2 diabetes-like phenotype in TALLYHO/JngJ mice

Type 2 Diabetes mellitus (T2DM) is a metabolic disorder characterized by chronic hyperglycemia, resulting from deficits in insulin secretion, insulin action, or both. Whilst the role of insulin in the peripheral nervous system has been ascertained in countless studies, its role in the central nervous system (CNS) is emerging only recently. Brain insulin has been lately associated with brain disorders like Alzheimer's disease, obsessive compulsive disorder, and attention deficit hyperactivity disorder. Thus, understanding the role of insulin as a common risk factor for mental and somatic comorbidities may disclose novel preventative and therapeutic approaches. We evaluated general metabolism (glucose tolerance, insulin sensitivity, energy expenditure, lipid metabolism, and polydipsia) and cognitive capabilities (attention, cognitive flexibility, and memory), in adolescent, young adult, and adult male and female TALLYHO/JngJ mice (TH, previously reported to constitute a valid experimental model of T2DM due to impaired insulin signaling). Adult TH mice have also been studied for alterations in gut microbiota diversity and composition. While TH mice exhibited profound deficits in cognitive flexibility and altered glucose metabolism, we observed that these alterations emerged either much earlier (males) or independent of (females) a comprehensive constellation of symptoms, isomorphic to an overt T2DM-like phenotype (insulin resistance, polydipsia, higher energy expenditure, and altered lipid metabolism). We also observed significant sex-dependent alterations in gut microbiota alpha diversity and taxonomy in adult TH mice. Deficits in insulin signaling may represent a common risk factor for both T2DM and CNS-related deficits, which may stem from (partly) independent mechanisms.

Tuberculosis susceptibility in genetically diverse mice reveals functional diversity of neutrophils

Tuberculosis (TB) is a heterogenous disease in humans with individuals exhibiting a wide range of susceptibility. This heterogeneity is not captured by standard laboratory mouse lines. We used a new collection of 19 wild-derived inbred mouse lines collected from diverse geographic sites to identify novel phenotypes during Mycobacterium tuberculosis ( Mtb ) infection. Wild derived mice have heterogenous immune responses to infection that result in differential ability to control disease at early timepoints. Correlation analysis with multiple parameters including sex, weight, and cellular immune responses in the lungs revealed that enhanced control of infection is associated with increased numbers of CD4 T cells, CD8 T cells and B cells. Surprisingly, we did not observe strong correlations between IFN-γ production and control of infection. Although in most lines high neutrophils were associated with susceptibility, we identified a mouse line that harbors high neutrophils numbers yet controls infection. Using single-cell RNA sequencing, we identified a novel neutrophil signature associated with failure to control infection.

Stem cell transcriptional profiles from mouse subspecies reveal cis -regulatory evolution at translation genes

A key goal of evolutionary genomics is to harness molecular data to draw inferences about selective forces that have acted on genomes. The field progresses in large part through the development of advanced molecular-evolution analysis methods. Here we explored the intersection between classical sequence-based tests for selection and an empirical expression-based approach, using stem cells from Mus musculus subspecies as a model. Using a test of directional, cis -regulatory evolution across genes in pathways, we discovered a unique program of induction of translation genes in stem cells of the Southeast Asian mouse M. m. castaneus relative to its sister taxa. We then mined population-genomic sequences to pursue underlying regulatory mechanisms for this expression divergence, finding robust evidence for alleles unique to M. m. castaneus at the upstream regions of the translation genes. We interpret our data under a model of changes in lineage-specific pressures across Mus musculus in stem cells with high translational capacity. Our findings underscore the rigor of integrating expression and sequence-based methods to generate hypotheses about evolutionary events from long ago.

Insertion of short L1 sequences generates inter-strain histone acetylation differences in the mouse

BACKGROUND

Gene expression divergence between populations and between individuals can emerge from genetic variations within the genes and/or in the cis regulatory elements. Since epigenetic modifications regulate gene expression, it is conceivable that epigenetic variations in cis regulatory elements can also be a source of gene expression divergence.

RESULTS

In this study, we compared histone acetylation (namely, H3K9ac) profiles in two mouse strains of different subspecies origin, C57BL/6 J (B6) and MSM/Ms (MSM), as well as their F1 hybrids. This identified 319 regions of strain-specific acetylation, about half of which were observed between the alleles of F1 hybrids. While the allele-specific presence of the interferon regulatory factor 3 (IRF3) binding sequence was associated with allele-specific histone acetylation, we also revealed that B6-specific insertions of a short 3' fragment of LINE-1 (L1) retrotransposon occur within or proximal to MSM-specific acetylated regions. Furthermore, even in hyperacetylated domains, flanking regions of non-polymorphic 3' L1 fragments were hypoacetylated, suggesting a general activity of the 3' L1 fragment to induce hypoacetylation. Indeed, we confirmed the binding of the 3' region of L1 by three Krüppel-associated box domain-containing zinc finger proteins (KZFPs), which interact with histone deacetylases. These results suggest that even a short insertion of L1 would be excluded from gene- and acetylation-rich regions by natural selection. Finally, mRNA-seq analysis for F1 hybrids was carried out, which disclosed a link between allele-specific promoter/enhancer acetylation and gene expression.

CONCLUSIONS

This study disclosed a number of genetic changes that have changed the histone acetylation levels during the evolution of mouse subspecies, a part of which is associated with gene expression changes. Insertions of even a very short L1 fragment can decrease the acetylation level in their neighboring regions and thereby have been counter-selected in gene-rich regions, which may explain a long-standing mystery of discrete genomic distribution of LINEs and SINEs.

Inner ear morphology in wild versus laboratory house mice

The semicircular canals of the inner ear are involved in balance and velocity control. Being crucial to ensure efficient mobility, their morphology exhibits an evolutionary conservatism attributed to stabilizing selection. Release of selection in slow-moving animals has been argued to lead to morphological divergence and increased inter-individual variation. In its natural habitat, the house mouse Mus musculus moves in a tridimensional space where efficient balance is required. In contrast, laboratory mice in standard cages are severely restricted in their ability to move, which possibly reduces selection on the inner ear morphology. This effect was tested by comparing four groups of mice: several populations of wild mice trapped in commensal habitats in France; their second-generation laboratory offspring, to assess plastic effects related to breeding conditions; a standard laboratory strain (Swiss) that evolved for many generations in a regime of mobility reduction; and hybrids between wild offspring and Swiss mice. The morphology of the semicircular canals was quantified using a set of 3D landmarks and semi-landmarks analyzed using geometric morphometric protocols. Levels of inter-population, inter-individual (disparity) and intra-individual (asymmetry) variation were compared. All wild mice shared a similar inner ear morphology, in contrast to the important divergence of the Swiss strain. The release of selection in the laboratory strain obviously allowed for an important and rapid drift in the otherwise conserved structure. Shared traits between the inner ear of the lab strain and domestic pigs suggested a common response to mobility reduction in captivity. The lab-bred offspring of wild mice also differed from their wild relatives, suggesting plastic response related to maternal locomotory behavior, since inner ear morphology matures before birth in mammals. The signature observed in lab-bred wild mice and the lab strain was however not congruent, suggesting that plasticity did not participate to the divergence of the laboratory strain. However, contrary to the expectation, wild mice displayed slightly higher levels of inter-individual variation than laboratory mice, possibly due to the higher levels of genetic variance within and among wild populations compared to the lab strain. Differences in fluctuating asymmetry levels were detected, with the laboratory strain occasionally displaying higher asymmetry scores than its wild relatives. This suggests that there may indeed be a release of selection and/or a decrease in developmental stability in the laboratory strain.

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.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for over a century. However, laboratory mice capture only a subset of the genetic variation found in wild mouse populations, ultimately limiting the potential of classical inbred strains to uncover phenotype-associated variants and pathways. Wild mouse populations are reservoirs of genetic diversity that could facilitate the discovery of new functional and disease-associated alleles, but the scarcity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently developed, sequenced, and phenotyped a set of 11 inbred strains derived from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from one of five environmentally distinct locations across North and South America. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the GRCm39 mouse reference, with 42.5% of variants in the Nachman strain genomes absent from current classical inbred mouse strain panels. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels. These novel wild-derived inbred mouse strain resources are set to empower new discoveries in both basic and preclinical research.

The Updated Mouse Universal Genotyping Array Bioinformatic Pipeline Improves Genetic QC in Laboratory Mice

The MiniMUGA genotyping array is a popular tool for genetic QC of laboratory mice and genotyping of samples from most types of experimental crosses involving laboratory strains, particularly for reduced complexity crosses. The content of the production version of the MiniMUGA array is fixed; however, there is the opportunity to improve array's performance and the associated report's usefulness by leveraging thousands of samples genotyped since the initial description of MiniMUGA in 2020. Here we report our efforts to update and improve marker annotation, increase the number and the reliability of the consensus genotypes for inbred strains and increase the number of constructs that can reliably be detected with MiniMUGA. In addition, we have implemented key changes in the informatics pipeline to identify and quantify the contribution of specific genetic backgrounds to the makeup of a given sample, remove arbitrary thresholds, include the Y Chromosome and mitochondrial genome in the ideogram, and improve robust detection of the presence of commercially available substrains based on diagnostic alleles. Finally, we have made changes to the layout of the report, to simplify the interpretation and completeness of the analysis and added a table summarizing the ideogram. We believe that these changes will be of general interest to the mouse research community and will be instrumental in our goal of improving the rigor and reproducibility of mouse-based biomedical research.

Female behavior drives the formation of distinct social structures in C57BL/6J versus wild-derived outbred mice in field enclosures

BACKGROUND

Social behavior and social organization have major influences on individual health and fitness. Yet, biomedical research focuses on studying a few genotypes under impoverished social conditions. Understanding how lab conditions have modified social organizations of model organisms, such as lab mice, relative to natural populations is a missing link between socioecology and biomedical science.

RESULTS

Using a common garden design, we describe the formation of social structure in the well-studied laboratory mouse strain, C57BL/6J, in replicated mixed-sex populations over 10-day trials compared to control trials with wild-derived outbred house mice in outdoor field enclosures. We focus on three key features of mouse social systems: (i) territory establishment in males, (ii) female social relationships, and (iii) the social networks formed by the populations. Male territorial behaviors were similar but muted in C57 compared to wild-derived mice. Female C57 sharply differed from wild-derived females, showing little social bias toward cage mates and exploring substantially more of the enclosures compared to all other groups. Female behavior consistently generated denser social networks in C57 than in wild-derived mice.

CONCLUSIONS

C57 and wild-derived mice individually vary in their social and spatial behaviors which scale to shape overall social organization. The repeatable societies formed under field conditions highlights opportunities to experimentally study the interplay between society and individual biology using model organisms.

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.

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.

The Expansion of House Mouse Major Urinary Protein Genes Likely Did Not Facilitate Commensalism with Humans

Mouse wild-derived strains (WDSs) combine the advantages of classical laboratory stocks and wild animals, and thus appear to be promising tools for diverse biomedical and evolutionary studies. We employed 18 WDSs representing three non-synanthropic species (Mus spretus, Mus spicilegus, and M. macedonicus) and three house mouse subspecies (Mus musculus musculus, M. m. domesticus, M. m. castaneus), which are all important human commensals to explore whether the number of major urinary protein (MUP) genes and their final protein levels in urine are correlated with the level of commensalism. Contrary to expectations, the MUP copy number (CN) and protein excretion in the strains derived from M. m. castaneus, which is supposed to be the strongest commensal, were not significantly different from the non-commensal species. Regardless of an overall tendency for higher MUP amounts in taxa with a higher CN, there was no significant correlation at the strain level. Our study thus suggests that expansion of the Mup cluster, which appeared before the house mouse diversification, is unlikely to facilitate commensalism with humans in three house mouse subspecies. Finally, we found considerable variation among con(sub)specific WDSs, warning against generalisations of results based on a few strains.

Network analysis identifies strain-dependent response to tau and tau seeding-associated genes

Previous research demonstrated that genetic heterogeneity is a critical factor in modeling amyloid accumulation and other Alzheimer's disease phenotypes. However, it is unknown what mechanisms underlie these effects of genetic background on modeling tau aggregate-driven pathogenicity. In this study, we induced tau aggregation in wild-derived mice by expressing MAPT. To investigate the effect of genetic background on the action of tau aggregates, we performed RNA sequencing with brains of C57BL/6J, CAST/EiJ, PWK/PhJ, and WSB/EiJ mice (n = 64) and determined core transcriptional signature conserved in all genetic backgrounds and signature unique to wild-derived backgrounds. By measuring tau seeding activity using the cortex, we identified 19 key genes associated with tau seeding and amyloid response. Interestingly, microglial pathways were strongly associated with tau seeding activity in CAST/EiJ and PWK/PhJ backgrounds. Collectively, our study demonstrates that mouse genetic context affects tau-mediated alteration of transcriptome and tau seeding. The gene modules associated with tau seeding provide an important resource to better model tauopathy.

Novel regulators of islet function identified from genetic variation in mouse islet Ca2+ oscillations

Insufficient insulin secretion to meet metabolic demand results in diabetes. The intracellular flux of Ca2+ into β-cells triggers insulin release. Since genetics strongly influences variation in islet secretory responses, we surveyed islet Ca2+ dynamics in eight genetically diverse mouse strains. We found high strain variation in response to four conditions: (1) 8 mM glucose; (2) 8 mM glucose plus amino acids; (3) 8 mM glucose, amino acids, plus 10 nM glucose-dependent insulinotropic polypeptide (GIP); and (4) 2 mM glucose. These stimuli interrogate β-cell function, α- to β-cell signaling, and incretin responses. We then correlated components of the Ca2+ waveforms to islet protein abundances in the same strains used for the Ca2+ measurements. To focus on proteins relevant to human islet function, we identified human orthologues of correlated mouse proteins that are proximal to glycemic-associated single-nucleotide polymorphisms in human genome-wide association studies. Several orthologues have previously been shown to regulate insulin secretion (e.g. ABCC8, PCSK1, and GCK), supporting our mouse-to-human integration as a discovery platform. By integrating these data, we nominate novel regulators of islet Ca2+ oscillations and insulin secretion with potential relevance for human islet function. We also provide a resource for identifying appropriate mouse strains in which to study these regulators.

Identification of genomic loci regulating platelet plasminogen activator inhibitor-1 in mice

BACKGROUND

Plasminogen activator inhibitor-1 (PAI-1, Serpine1) is an important circulating fibrinolysis inhibitor. PAI-1 exists in 2 pools, packaged within platelet α-granules and freely circulating in plasma. Elevated plasma PAI-1 levels are associated with cardiovascular disease. However, little is known about the regulation of platelet PAI-1 (pPAI-1).

OBJECTIVES

We investigated the genetic control of pPAI-1 levels in mice and humans.

METHODS

We measured pPAI-1 antigen levels via enzyme-linked immunosorbent assay in platelets isolated from 10 inbred mouse strains, including LEWES/EiJ (LEWES) and C57BL/6J (B6). LEWES and B6 were crossed to produce the F1 generation, B6LEWESF1. B6LEWESF1 mice were intercrossed to produce B6LEWESF2 mice. These mice were subjected to genome-wide genetic marker genotyping followed by quantitative trait locus analysis to identify pPAI-1 regulatory loci.

RESULTS

We identified differences in pPAI-1 between several laboratory strains, with LEWES having pPAI-1 levels more than 10-fold higher than those in B6. Quantitative trait locus analysis of B6LEWESF2 offspring identified a major pPAI-1 regulatory locus on chromosome 5 from 136.1 to 137.6 Mb (logarithm of the odds score, 16.2). Significant pPAI-1 modifier loci on chromosomes 6 and 13 were also identified.

CONCLUSION

Identification of pPAI-1 genomic regulatory elements provides insights into platelet/megakaryocyte-specific and cell type-specific gene expression. This information can be used to design more precise therapeutic targets for diseases where PAI-1 plays a role.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for a century. However, laboratory mice capture a narrow subset of the genetic variation found in wild mouse populations. This consideration inherently restricts the scope of potential discovery in laboratory models and narrows the pool of potentially identified phenotype-associated variants and pathways. Wild mouse populations are reservoirs of predicted functional and disease-associated alleles, but the sparsity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently imported, sequenced, and phenotyped a set of 11 wild-derived inbred strains developed from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from five environmentally diverse locations across North and South America: Saratoga Springs, New York, USA; Gainesville, Florida, USA; Manaus, Brazil; Tucson, Arizona, USA; and Edmonton, Alberta, Canada. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the mouse reference assembly, with 42.5% of variants in the Nachman strain genomes absent from classical inbred mouse strains. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels alone. Taken together, our work introduces a novel wild-derived inbred mouse strain resource that will enable new discoveries in basic and preclinical research. These strains are currently available through The Jackson Laboratory Repository under laboratory code NachJ.

GenomeMUSter mouse genetic variation service enables multi-trait, multi-population data integration and analyses

Hundreds of inbred laboratory mouse strains and intercross populations have been used to functionalize 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 the Collaborative Cross, expanded BXD and inbred wild-derived strains add to set of complex disease mouse models, genetic mapping resources and sensitized backgrounds against which to evaluate engineered mutations. The genome sequences of many inbred strains, along with dense genotypes from others could allow integrated analysis of trait - variant associations across populations, but these analyses are not feasible due to the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense data resource by harmonizing multiple variant datasets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extensible to other model organism species. The result is a web- and programmatically-accessible data service called GenomeMUSter ( https://muster.jax.org ), comprising allelic data covering 657 strains at 106.8M segregating sites. Interoperation with phenotype databases, analytic tools and other resources enable a wealth of applications including multi-trait, multi-population meta-analysis. We demonstrate this in a cross-species comparison of the meta-analysis of Type 2 Diabetes and of substance use disorders, resulting in the more specific characterization of the role of human variant effects in light of mouse phenotype data. 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.

Seven naturally variant loci serve as genetic modifiers of Lamc2jeb induced non-Herlitz junctional Epidermolysis Bullosa in mice

Epidermolysis Bullosa (EB) is a group of rare genetic disorders that compromise the structural integrity of the skin such that blisters and subsequent erosions occur after minor trauma. While primary genetic risk of all subforms of EB adhere to Mendelian patterns of inheritance, their clinical presentations and severities can vary greatly, implying genetic modifiers. The Lamc2jeb mouse model of non-Herlitz junctional EB (JEB-nH) demonstrated that genetic modifiers can contribute substantially to the phenotypic variability of JEB and likely other forms of EB. The innocuous changes in an 'EB related gene', Col17a1, have shown it to be a dominant modifier of Lamc2jeb. This work identifies six additional Quantitative Trait Loci (QTL) that modify disease in Lamc2jeb/jeb mice. Three QTL include other known 'EB related genes', with the strongest modifier effect mapping to a region including the epidermal hemi-desmosomal structural gene dystonin (Dst-e/Bpag1-e). Three other QTL map to intervals devoid of known EB-associated genes. Of these, one contains the nuclear receptor coactivator Ppargc1a as its primary candidate and the others contain related genes Pparg and Igf1, suggesting modifier pathways. These results, demonstrating the potent disease modifying effects of normally innocuous genetic variants, greatly expand the landscape of genetic modifiers of EB and therapeutic approaches that may be applied.

Which mouse multiparental population is right for your study? The Collaborative Cross inbred strains, their F1 hybrids, or the Diversity Outbred population

Multiparental populations (MPPs) encompass greater genetic diversity than traditional experimental crosses of two inbred strains, enabling broader surveys of genetic variation underlying complex traits. Two such mouse MPPs are the Collaborative Cross (CC) inbred panel and the Diversity Outbred (DO) population, which are descended from the same eight inbred strains. Additionally, the F1 intercrosses of CC strains (CC-RIX) have been used and enable study designs with replicate outbred mice. Genetic analyses commonly used by researchers to investigate complex traits in these populations include characterizing how heritable a trait is, i.e. its heritability, and mapping its underlying genetic loci, i.e. its quantitative trait loci (QTLs). Here we evaluate the relative merits of these populations for these tasks through simulation, as well as provide recommendations for performing the quantitative genetic analyses. We find that sample populations that include replicate animals, as possible with the CC and CC-RIX, provide more efficient and precise estimates of heritability. We report QTL mapping power curves for the CC, CC-RIX, and DO across a range of QTL effect sizes and polygenic backgrounds for samples of 174 and 500 mice. The utility of replicate animals in the CC and CC-RIX for mapping QTLs rapidly decreased as traits became more polygenic. Only large sample populations of 500 DO mice were well-powered to detect smaller effect loci (7.5-10%) for highly complex traits (80% polygenic background). All results were generated with our R package musppr, which we developed to simulate data from these MPPs and evaluate genetic analyses from user-provided genotypes.

Multi-omics analysis identifies drivers of protein phosphorylation

BACKGROUND

Phosphorylation of proteins is a key step in the regulation of many cellular processes including activation of enzymes and signaling cascades. The abundance of a phosphorylated peptide (phosphopeptide) is determined by the abundance of its parent protein and the proportion of target sites that are phosphorylated.

RESULTS

We quantified phosphopeptides, proteins, and transcripts in heart, liver, and kidney tissue samples of mice from 58 strains of the Collaborative Cross strain panel. We mapped ~700 phosphorylation quantitative trait loci (phQTL) across the three tissues and applied genetic mediation analysis to identify causal drivers of phosphorylation. We identified kinases, phosphatases, cytokines, and other factors, including both known and potentially novel interactions between target proteins and genes that regulate site-specific phosphorylation. Our analysis highlights multiple targets of pyruvate dehydrogenase kinase 1 (PDK1), a regulator of mitochondrial function that shows reduced activity in the NZO/HILtJ mouse, a polygenic model of obesity and type 2 diabetes.

CONCLUSIONS

Together, this integrative multi-omics analysis in genetically diverse CC strains provides a powerful tool to identify regulators of protein phosphorylation. The data generated in this study provides a resource for further exploration.

Next-generation PFAS 6:2 fluorotelomer sulfonate reduces plaque formation in exposed white-footed mice

6:2 fluorotelomer sulfonate (6:2 FTS) has been used as a replacement for legacy per- and polyfluoroalkyl substances (PFAS). We assessed reproductive and developmental effects in a human-wildlife hybrid animal model based on the association of adverse effects linked to legacy PFAS with these sensitive life stages. In this study, white-footed mice were exposed orally to 0, 0.2, 1, 5, or 25 mg/kg-day 6:2 FTS for 112 days (4 weeks premating exposure plus at least 4 weeks mating exposure). Pregnancy and fertility indices were calculated, and litter production, total litter size, live litter size, stillbirths, litter loss, average pup weight, and pinna unfolding were assessed. Sex steroid and thyroid hormone serum levels were assessed. Body weight, histopathology, and immune function were also assessed in this study. Reproductive endpoints were not significantly altered in response to 6:2 FTS. Spleen weight increased in male mice dosed with 6:2 FTS. Immune function determined via a plaque-forming cell (PFC) assay was decreased in both male and female mice in the 2 highest doses. A low benchmark dose was calculated based on PFCs as the critical effect and was found to be 2.63 and 2.26 mg/kg-day 6:2 FTS in male and female mice, respectively. This study characterizes 6:2 FTS as being potentially immunotoxic with little evidence of effect on reproduction and development; furthermore, it models acceptable levels of exposure. These 2 pieces of information together will aid regulators in setting environmental exposure limits for this PFAS currently thought to be less toxic than other PFAS.

A genetic locus complements resistance to Bordetella pertussis-induced histamine sensitization

Histamine plays pivotal role in normal physiology and dysregulated production of histamine or signaling through histamine receptors (HRH) can promote pathology. Previously, we showed that Bordetella pertussis or pertussis toxin can induce histamine sensitization in laboratory inbred mice and is genetically controlled by Hrh1/HRH1. HRH1 allotypes differ at three amino acid residues with P263-V313-L331 and L263-M313-S331, imparting sensitization and resistance respectively. Unexpectedly, we found several wild-derived inbred strains that carry the resistant HRH1 allotype (L263-M313-S331) but exhibit histamine sensitization. This suggests the existence of a locus modifying pertussis-dependent histamine sensitization. Congenic mapping identified the location of this modifier locus on mouse chromosome 6 within a functional linkage disequilibrium domain encoding multiple loci controlling sensitization to histamine. We utilized interval-specific single-nucleotide polymorphism (SNP) based association testing across laboratory and wild-derived inbred mouse strains and functional prioritization analyses to identify candidate genes for this modifier locus. Atg7, Plxnd1, Tmcc1, Mkrn2, Il17re, Pparg, Lhfpl4, Vgll4, Rho and Syn2 are candidate genes within this modifier locus, which we named Bphse, enhancer of Bordetella pertussis induced histamine sensitization. Taken together, these results identify, using the evolutionarily significant diversity of wild-derived inbred mice, additional genetic mechanisms controlling histamine sensitization.

Reproductive genomics of the mouse: implications for human fertility and infertility

Genetic analyses of mammalian gametogenesis and fertility have the potential to inform about two important and interrelated clinical areas: infertility and contraception. Here, we address the genetics and genomics underlying gamete formation, productivity and function in the context of reproductive success in mammalian systems, primarily mouse and human. Although much is known about the specific genes and proteins required for meiotic processes and sperm function, we know relatively little about other gametic determinants of overall fertility, such as regulation of gamete numbers, duration of gamete production, and gamete selection and function in fertilization. As fertility is not a binary trait, attention is now appropriately focused on the oligogenic, quantitative aspects of reproduction. Multiparent mouse populations, created by complex crossing strategies, exhibit genetic diversity similar to human populations and will be valuable resources for genetic discovery, helping to overcome current limitations to our knowledge of mammalian reproductive genetics. Finally, we discuss how what we know about the genomics of reproduction can ultimately be brought to the clinic, informing our concepts of human fertility and infertility, and improving assisted reproductive technologies.

Network analysis reveals strain-dependent response to misfolded tau aggregates

Mouse genetic backgrounds have been shown to modulate amyloid accumulation and propagation of tau aggregates. Previous research into these effects has highlighted the importance of studying the impact of genetic heterogeneity on modeling Alzheimer's disease. However, it is unknown what mechanisms underly these effects of genetic background on modeling Alzheimer's disease, specifically tau aggregate-driven pathogenicity. In this study, we induced tau aggregation in wild-derived mice by expressing MAPT (P301L). To investigate the effect of genetic background on the action of tau aggregates, we performed RNA sequencing with brains of 6-month-old C57BL/6J, CAST/EiJ, PWK/PhJ, and WSB/EiJ mice (n=64). We also measured tau seeding activity in the cortex of these mice. We identified three gene signatures: core transcriptional signature, unique signature for each wild-derived genetic background, and tau seeding-associated signature. Our data suggest that microglial response to tau seeds is elevated in CAST/EiJ and PWK/PhJ mice. Together, our study provides the first evidence that mouse genetic context influences the seeding of tau.

SUMMARY

Seeding of tau predates the phosphorylation and spreading of tau aggregates. Acri and colleagues report transcriptomic responses to tau and elevated tau seeds in wild-derived mice. This paper creates a rich resource by combining genetics, tau biosensor assays, and transcriptomics.

Mouse Phenome Database: towards a more FAIR-compliant and TRUST-worthy data repository and tool suite for phenotypes and genotypes

The Mouse Phenome Database (MPD; https://phenome.jax.org; RRID:SCR_003212), supported by the US National Institutes of Health, is a Biomedical Data Repository listed in the Trans-NIH Biomedical Informatics Coordinating Committee registry. As an increasingly FAIR-compliant and TRUST-worthy data repository, MPD accepts phenotype and genotype data from mouse experiments and curates, organizes, integrates, archives, and distributes those data using community standards. Data are accompanied by rich metadata, including widely used ontologies and detailed protocols. Data are from all over the world and represent genetic, behavioral, morphological, and physiological disease-related characteristics in mice at baseline or those exposed to drugs or other treatments. MPD houses data from over 6000 strains and populations, representing many reproducible strain types and heterogenous populations such as the Diversity Outbred where each mouse is unique but can be genotyped throughout the genome. A suite of analysis tools is available to aggregate, visualize, and analyze these data within and across studies and populations in an increasingly traceable and reproducible manner. We have refined existing resources and developed new tools to continue to provide users with access to consistent, high-quality data that has translational relevance in a modernized infrastructure that enables interaction with a suite of bioinformatics analytic and data services.

Rapid Evolution of the Fine-scale Recombination Landscape in Wild House Mouse (Mus musculus) Populations

Meiotic recombination is an important evolutionary force and an essential meiotic process. In many species, recombination events concentrate into hotspots defined by the site-specific binding of PRMD9. Rapid evolution of Prdm9's zinc finger DNA-binding array leads to remarkably abrupt shifts in the genomic distribution of hotspots between species, but the question of how Prdm9 allelic variation shapes the landscape of recombination between populations remains less well understood. Wild house mice (Mus musculus) harbor exceptional Prdm9 diversity, with >150 alleles identified to date, and pose a particularly powerful system for addressing this open question. We employed a coalescent-based approach to construct broad- and fine-scale sex-averaged recombination maps from contemporary patterns of linkage disequilibrium in nine geographically isolated wild house mouse populations, including multiple populations from each of three subspecies. Comparing maps between wild mouse populations and subspecies reveals several themes. First, we report weak fine- and broad-scale recombination map conservation across subspecies and populations, with genetic divergence offering no clear prediction for recombination map divergence. Second, most hotspots are unique to one population, an outcome consistent with minimal sharing of Prdm9 alleles between surveyed populations. Finally, by contrasting aggregate hotspot activity on the X versus autosomes, we uncover evidence for population-specific differences in the degree and direction of sex dimorphism for recombination. Overall, our findings illuminate the variability of both the broad- and fine-scale recombination landscape in M. musculus and underscore the functional impact of Prdm9 allelic variation in wild mouse populations.

Taxonomic assessment of two wild house mouse subspecies using whole-genome sequencing

The house mouse species complex (Mus musculus) is comprised of three primary subspecies. A large number of secondary subspecies have also been suggested on the basis of divergent morphology and molecular variation at limited numbers of markers. While the phylogenetic relationships among the primary M. musculus subspecies are well-defined, relationships among secondary subspecies and between secondary and primary subspecies remain less clear. Here, we integrate de novo genome sequencing of museum-stored specimens of house mice from one secondary subspecies (M. m. bactrianus) and publicly available genome sequences of house mice previously characterized as M. m. helgolandicus, with whole genome sequences from diverse representatives of the three primary house mouse subspecies. We show that mice assigned to the secondary M. m. bactrianus and M. m. helgolandicus subspecies are not genetically differentiated from M. m. castaneus and M. m. domesticus, respectively. Overall, our work suggests that the M. m. bactrianus and M. m. helgolandicus subspecies are not well-justified taxonomic entities, emphasizing the importance of leveraging whole-genome sequence data to inform subspecies designations. Additionally, our investigation provides tailored experimental procedures for generating whole genome sequences from air-dried mouse skins, along with key genomic resources to inform future genomic studies of wild mouse diversity.

The contribution of sex chromosome conflict to disrupted spermatogenesis in hybrid house mice

Incompatibilities on the sex chromosomes are important in the evolution of hybrid male sterility, but the evolutionary forces underlying this phenomenon are unclear. House mice (Mus musculus) lineages have provided powerful models for understanding the genetic basis of hybrid male sterility. X chromosome-autosome interactions cause strong incompatibilities in M. musculus F1 hybrids, but variation in sterility phenotypes suggests a more complex genetic basis. In addition, XY chromosome conflict has resulted in rapid expansions of ampliconic genes with dosage-dependent expression that is essential to spermatogenesis. Here, we evaluated the contribution of XY lineage mismatch to male fertility and stage-specific gene expression in hybrid mice. We performed backcrosses between two house mouse subspecies to generate reciprocal Y-introgression strains and used these strains to test the effects of XY mismatch in hybrids. Our transcriptome analyses of sorted spermatid cells revealed widespread overexpression of the X chromosome in sterile F1 hybrids independent of Y chromosome subspecies origin. Thus, postmeiotic overexpression of the X chromosome in sterile F1 mouse hybrids is likely a downstream consequence of disrupted meiotic X-inactivation rather than XY gene copy number imbalance. Y chromosome introgression did result in subfertility phenotypes and disrupted expression of several autosomal genes in mice with an otherwise nonhybrid genomic background, suggesting that Y-linked incompatibilities contribute to reproductive barriers, but likely not as a direct consequence of XY conflict. Collectively, these findings suggest that rapid sex chromosome gene family evolution driven by genomic conflict has not resulted in strong male reproductive barriers between these subspecies of house mice.

Architecture of the outbred brown fat proteome defines regulators of metabolic physiology

Brown adipose tissue (BAT) regulates metabolic physiology. However, nearly all mechanistic studies of BAT protein function occur in a single inbred mouse strain, which has limited the understanding of generalizable mechanisms of BAT regulation over physiology. Here, we perform deep quantitative proteomics of BAT across a cohort of 163 genetically defined diversity outbred mice, a model that parallels the genetic and phenotypic variation found in humans. We leverage this diversity to define the functional architecture of the outbred BAT proteome, comprising 10,479 proteins. We assign co-operative functions to 2,578 proteins, enabling systematic discovery of regulators of BAT. We also identify 638 proteins that correlate with protection from, or sensitivity to, at least one parameter of metabolic disease. We use these findings to uncover SFXN5, LETMD1, and ATP1A2 as modulators of BAT thermogenesis or adiposity, and provide OPABAT as a resource for understanding the conserved mechanisms of BAT regulation over metabolic physiology.

A Locus on Chromosome 15 Contributes to Acute Ozone-induced Lung Injury in Collaborative Cross Mice

Ozone (O3)-induced respiratory toxicity varies considerably within the human population and across inbred mouse strains, indicative of gene-environment interactions (GxE). Though previous studies have identified several quantitative trait loci (QTL) and candidate genes underlying responses to O3 exposure, precise mechanisms of susceptibility remain incompletely described. We sought to update our understanding of the genetic architecture of O3 responsiveness using the Collaborative Cross (CC) recombinant inbred mouse panel. We evaluated hallmark O3-induced inflammation and injury phenotypes in 56 CC strains after exposure to filtered air or 2 ppm O3, and performed focused genetic analysis of variation in lung injury, as reflected by protein in lung lavage fluid. Strain-dependent responses to O3 were clear, and QTL mapping revealed two novel loci on Chr (Chromosomes) 10 (peak, 26.2 Mb; 80% confidence interval [CI], 24.6-43.6 Mb) and 15 (peak, 47.1 Mb; 80% CI, 40.2-54.9 Mb), the latter surpassing the 95% significance threshold. At the Chr 15 locus, C57BL/6J and CAST/EiJ founder haplotypes were associated with higher lung injury responses compared with all other CC founder haplotypes. With further statistical analysis and a weight of evidence approach, we delimited the Chr 15 QTL to an ∼2 Mb region containing 21 genes (10 protein coding) and nominated three candidate genes, namely Oxr1, Rspo2, and Angpt1. Gene and protein expression data further supported Oxr1 and Angpt1 as priority candidate genes. In summary, we have shown that O3-induced lung injury is modulated by genetic variation, identified two high priority candidate genes, and demonstrated the value of the CC for detecting GxE.

Whole-genome sequencing reveals the genetic mechanisms of domestication in classical inbred mice

Background

The laboratory mouse was domesticated from the wild house mouse. Understanding the genetics underlying domestication in laboratory mice, especially in the widely used classical inbred mice, is vital for studies using mouse models. However, the genetic mechanism of laboratory mouse domestication remains unknown due to lack of adequate genomic sequences of wild mice.

Results

We analyze the genetic relationships by whole-genome resequencing of 36 wild mice and 36 inbred strains. All classical inbred mice cluster together distinctly from wild and wild-derived inbred mice. Using nucleotide diversity analysis, Fst, and XP-CLR, we identify 339 positively selected genes that are closely associated with nervous system function. Approximately one third of these positively selected genes are highly expressed in brain tissues, and genetic mouse models of 125 genes in the positively selected genes exhibit abnormal behavioral or nervous system phenotypes. These positively selected genes show a higher ratio of differential expression between wild and classical inbred mice compared with all genes, especially in the hippocampus and frontal lobe. Using a mutant mouse model, we find that the SNP rs27900929 (T>C) in gene Astn2 significantly reduces the tameness of mice and modifies the ratio of the two Astn2 (a/b) isoforms.

Conclusion

Our study indicates that classical inbred mice experienced high selection pressure during domestication under laboratory conditions. The analysis shows the positively selected genes are closely associated with behavior and the nervous system in mice. Tameness may be related to the Astn2 mutation and regulated by the ratio of the two Astn2 (a/b) isoforms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02772-1.

Population structure and inbreeding in wild house mice (Mus musculus) at different geographic scales

House mice (Mus musculus) have spread globally as a result of their commensal relationship with humans. In the form of laboratory strains, both inbred and outbred, they are also among the most widely used model organisms in biomedical research. Although the general outlines of house mouse dispersal and population structure are well known, details have been obscured by either limited sample size or small numbers of markers. Here we examine ancestry, population structure, and inbreeding using SNP microarray genotypes in a cohort of 814 wild mice spanning five continents and all major subspecies of Mus, with a focus on M. m. domesticus. We find that the major axis of genetic variation in M. m. domesticus is a south-to-north gradient within Europe and the Mediterranean. The dominant ancestry component in North America, Australia, New Zealand, and various small offshore islands are of northern European origin. Next we show that inbreeding is surprisingly pervasive and highly variable, even between nearby populations. By inspecting the length distribution of homozygous segments in individual genomes, we find that inbreeding in commensal populations is mostly due to consanguinity. Our results offer new insight into the natural history of an important model organism for medicine and evolutionary biology.

Reproductive and developmental toxicity of perfluorooctane sulfonate (PFOS) in the white-footed mouse (Peromyscus leucopus)

Concerns about per- and polyfluoroalkyl substances (PFAS) stem from their ubiquitous presence in the environment, bioaccumulation, resistance to degradation, and toxicity. Previously, toxicity data relevant to ecological risk assessment has largely been aquatic, terrestrial invertebrates, or avian in origin. In this study, repeated oral exposures of perfluorooctane sulfonate (PFOS) were administered to white-footed mice (Peromyscus leucopus) to evaluate effects on reproduction and development. Prenatal exposure to high doses of PFOS caused neonatal mortality, though growth and development were unaffected by low doses. Additionally, parental (P) generation animals exhibited increased liver:body weight, increased hepatocyte cytoplasmic vacuolization, and decreased serum thyroxine (T4) levels. Total litter loss was selected as the protective critical effect in this study resulting in a benchmark dose low (BMDL) of 0.12 mg/kg-d PFOS. Importantly, PFOS exposure has been linked to reduced adult recruitment in myriad species and at similar thresholds to this study. Similarities in critical/toxicologic effects across taxa may add confidence in risk assessments at sites with multiple taxa or environments.