Content and Performance of the MiniMUGA Genotyping Array: A New Tool To Improve Rigor and Reproducibility in Mouse Research

The mutant mouse resource and research center (MMRRC) consortium: the US-based public mouse repository system

Now in its 25th year, the Mutant Mouse Resource and Research Center (MMRRC) consortium continues to serve the United States and international biomedical scientific community as a public repository and distribution archive of laboratory mouse models of human disease for research. Supported by the National Institutes of Health (NIH), the MMRRC consists of 4 regionally distributed and dedicated vivaria, offices, and specialized laboratory facilities and an Informatics Coordination and Service Center (ICSC). The overarching purpose of the MMRRC is to facilitate groundbreaking biomedical research by offering an extensive repertoire of mutant mice that are essential for advancing the understanding of human physiology and disease. The function of the MMRRC is to identify, acquire, evaluate, characterize, cryopreserve, and distribute mutant mouse strains to qualified biomedical investigators around the nation and the globe. Mouse strains accepted from the research community are held to the highest scientific standards to optimize reproducibility and enhance scientific rigor and transparency. All submitted strains are thoroughly reviewed, documented, and validated using extensive scientific quality control measures. In addition, the MMRRC conducts resource-related research on cryopreservation, mouse genetics, environmental conditions, and other topics that enhance operations of the MMRRC. Today, the MMRRC maintains an archive of mice, cryopreserved embryos and sperm, embryonic stem (ES) cell lines, and murine hybridomas for nearly 65,000 alleles. Since its inception, the MMRRC has fulfilled more than 20,000 orders from 13,651 scientists at 8441 institutions worldwide. The MMRRC also provides numerous services to assist researchers, including scientific consultation, technical assistance, genetic assays, microbiome analysis, analytical phenotyping, pathology, cryorecovery, husbandry, breeding and colony management, infectious disease surveillance, and disease modeling. The ICSC coordinates MMRRC operations, interacts with researchers, and manages the website (mmrrc.org) and online catalogue. Researchers benefit from an expansive list of well-defined mouse models of disease that meet the highest scientific standards while submitting investigators benefit by having their mouse strains cryopreserved, protected, and distributed in compliance with NIH policies.

Genetic background affects the strength of crossover interference in house mice

Meiotic recombination is required for faithful chromosome segregation in most sexually reproducing organisms and shapes the distribution of genetic variation in populations. Both the overall rate and the spatial distribution of crossovers vary within and between species. Adjacent crossovers on the same chromosome tend to be spaced more evenly than expected at random, a phenomenon known as crossover interference. Although interference has been observed in many taxa, the factors that influence the strength of interference are not well understood. We used house mice (Mus musculus), a well-established model system for understanding recombination, to study the effects of genetics and age on recombination rate and interference in the male germline. We analyzed crossover positions in 503 progeny from reciprocal F1 hybrids between inbred strains representing the three major subspecies of house mice. Consistent with previous studies, autosomal alleles from M. m. musculus tend to increase recombination rate, while inheriting a M. m. musculus X chromosome decreases recombination rate. Old males transmit an average of 0.6 more crossovers per meiosis (5.0%) than young males, though the effect varies across genetic backgrounds. We show that the strength of crossover interference depends on genotype, providing a rare demonstration that interference evolves over short timescales. Differences between reciprocal F1s suggest that X-linked factors modulate the strength of interference. Our findings motivate additional comparisons of interference among recently diverged species and further examination of the role of paternal age in determining the number and positioning of crossovers.

Single nucleotide polymorphism (SNP) characterisation of mouse inbred strains bred at MRC-National Institute for Medical Research

Inbred mouse strains have long proved useful as tools for biomedical research. They remove the effects of genetic background as an experimental variable. Within all mouse colonies, genetic drift is a recognised phenomenon and monitoring and documenting changes is important for experimental design and consistency. This communication documents the initial characterisation through SNP analysis of the inbred mouse strains bred and used at the time at the Medical Research Council National Institute for Medical Research (MRC-NIMR), Mill Hill, now The Crick Institute. These inbred strains were part of the foundation colonies for the many genetically modified mouse strains made at Mill Hill. We found small genetic changes in four of the nine inbred strains. Although phenotypic differences have not yet been found between the NIMR and the correspondent parental strains, I cannot discard that these may arise or have already arisen. This work has also authenticated the 129/SvJEvNimr-Gpi1c strain that was widely used at MRC-NIMR for gene targeting. All these inbred strains have been renamed according to The International Committee on Standardized Genetic Nomenclature for Mice.

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.

A Minimal Hybrid Sterility Genome Assembled by Chromosome Swapping Between Mouse Subspecies (Mus musculus)

Hybrid sterility is a reproductive isolation barrier between diverging taxa securing the early steps of speciation. Hybrid sterility is ubiquitous in the animal and plant kingdoms, but its genetic control is poorly understood. In our previous studies, we have uncovered the sterility of hybrids between musculus and domesticus subspecies of the house mouse, which is controlled by the Prdm9 gene, the X-linked Hstx2 locus, and subspecific heterozygosity for genetic background. To further investigate this form of genic-driven chromosomal sterility, we constructed a simplified hybrid sterility model within the genome of the domesticus subspecies by swapping domesticus autosomes with their homologous partners from the musculus subspecies. We show that the "sterility" allelic combination of Prdm9 and Hstx2 can be activated by a musculus/domesticus heterozygosity of as few as two autosomes, Chromosome 17 (Chr 17) and Chr 18 and is further enhanced when another heterosubspecific autosomal pair is present, whereas it has no effect on meiotic progression in the pure domesticus genome. In addition, we identify a new X-linked hybrid sterility locus, Hstx3, at the centromeric end of Chr X, which modulates the incompatibility between Prdm9 and Hstx2. These results further support our concept of chromosomal hybrid sterility based on evolutionarily accumulated divergence between homologous sequences. Based on these and previous results, we believe that future studies should include more information on the mutual recognition of homologous chromosomes at or before the first meiotic prophase in interspecific hybrids, as this may serve as a general reproductive isolation checkpoint in mice and other species.

Adam19 Deficiency Impacts Pulmonary Function: Human GWAS Follow-up in a Mouse Knockout Model

PURPOSE

Over 550 loci have been associated with human pulmonary function in genome-wide association studies (GWAS); however, the causal role of most remains uncertain. Single nucleotide polymorphisms in a disintegrin and metalloprotease domain 19 (ADAM19) are consistently related to pulmonary function in GWAS. Thus, we used a mouse model to investigate the causal link between Adam19 and pulmonary function.

METHODS

We created an Adam19 knockout (KO) mouse model and validated the gene targeting using RNA-Seq and RT-qPCR. Mouse body composition was assessed using dual-energy X-ray absorptiometry. Mouse lung function was measured using flexiVent.

RESULTS

Contrary to prior publications, the KO was not neonatal lethal. KO mice had lower body weight and shorter tibial length than wild-type (WT) mice. Their body composition revealed lower soft weight, fat weight, and bone mineral content. Adam19 KO had decreased baseline respiratory system elastance, minute work of breathing, tissue damping, tissue elastance, and forced expiratory flow at 50% forced vital capacity but higher FEV0.1 and FVC. Adam19 KO had attenuated tissue damping and tissue elastance in response to methacholine following LPS exposure. Adam19 KO also exhibited attenuated neutrophil extravasation into the airway after LPS administration compared to WT. RNA-Seq analysis of KO and WT lungs identified several differentially expressed genes (Cd300lg, Kpna2, and Pttg1) implicated in lung biology and pathogenesis. Gene set enrichment analysis identified negative enrichment for TNF pathways.

CONCLUSION

Our murine findings support a causal role of ADAM19, implicated in human GWAS, in regulating pulmonary function.

Genetic background affects the strength of crossover interference in house mice

Meiotic recombination is required for faithful chromosome segregation in most sexually reproducing organisms and shapes the distribution of genetic variation in populations. Both the overall rate and the spatial distribution of crossovers vary within and between species. Adjacent crossovers on the same chromosome tend to be spaced more evenly than expected at random, a phenomenon known as crossover interference. Although interference has been observed in many taxa, the factors that influence the strength of interference are not well understood. We used house mice (Mus musculus), a well-established model system for understanding recombination, to study the effects of genetics and age on recombination rate and interference in the male germline. We analyzed crossover positions in 503 progeny from reciprocal F1 hybrids between inbred strains representing the three major subspecies of house mice. Consistent with previous studies, autosomal alleles from M. m. musculus tend to increase recombination rate, while inheriting a M. m. musculus X chromosome decreases recombination rate. Old males transmit an average of 0.6 more crossovers per meiosis (5.0%) than young males, though the effect varies across genetic backgrounds. We show that the strength of crossover interference depends on genotype, providing a rare demonstration that interference evolves over short timescales. Differences between reciprocal F1s suggest that X-linked factors modulate the strength of interference. Our findings motivate additional comparisons of interference among recently diverged species and further examination of the role of paternal age in determining the number and positioning of crossovers.

Sarbecovirus disease susceptibility is conserved across viral and host models

Coronaviruses have caused three severe epidemics since the start of the 21st century: SARS, MERS and COVID-19. The severity of the ongoing COVID-19 pandemic and increasing likelihood of future coronavirus outbreaks motivates greater understanding of factors leading to severe coronavirus disease. We screened ten strains from the Collaborative Cross mouse genetic reference panel and identified strains CC006/TauUnc (CC006) and CC044/Unc (CC044) as coronavirus-susceptible and resistant, respectively, as indicated by variable weight loss and lung congestion scores four days post-infection. We generated a genetic mapping population of 755 CC006xCC044 F2 mice and exposed the mice to one of three genetically distinct mouse-adapted coronaviruses: clade 1a SARS-CoV MA15 (n=391), clade 1b SARS-CoV-2 MA10 (n=274), and clade 2 HKU3-CoV MA (n=90). Quantitative trait loci (QTL) mapping in SARS-CoV MA15- and SARS-CoV-2 MA10-infected F2 mice identified genetic loci associated with disease severity. Specifically, we identified seven loci associated with variation in outcome following infection with either virus, including one, HrS43, that is present in both groups. Three of these QTL, including HrS43, were also associated with HKU3-CoV MA outcome. HrS43 overlaps with a QTL previously reported by our lab that is associated with SARS-CoV MA15 outcome in CC011xCC074 F2 mice and is also syntenic with a human chromosomal region associated with severe COVID-19 outcomes in humans GWAS. The results reported here provide: (a) additional support for the involvement of this locus in SARS-CoV MA15 infection, (b) the first conclusive evidence that this locus is associated with susceptibility across the Sarbecovirus subgenus, and (c) demonstration of the relevance of mouse models in the study of coronavirus disease susceptibility in humans.

Identification of a genetic region linked to tolerance to MRSA infection using Collaborative Cross mice

Staphylococcus aureus (S. aureus) colonizes humans asymptomatically but can also cause opportunistic infections, ranging from mild skin infections to severe life-threatening conditions. Resistance and tolerance are two ways a host can survive an infection. Resistance is limiting the pathogen burden, while tolerance is limiting the health impact of a given pathogen burden. In previous work, we established that collaborative cross (CC) mouse line CC061 is highly susceptible to Methicillin-resistant S. aureus infection (MRSA, USA300), while CC024 is tolerant. To identify host genes involved in tolerance after S. aureus infection, we crossed CC061 mice and CC024 mice to generate F1 and F2 populations. Survival after MRSA infection in the F1 and F2 generations was 65% and 55% and followed a complex dominant inheritance pattern for the CC024 increased survival phenotype. Colonization in F2 animals was more extreme than in their parents, suggesting successful segregation of genetic factors. We identified a Quantitative Trait Locus (QTL) peak on chromosome 7 for survival and weight change after infection. In this QTL, the WSB/EiJ (WSB) allele was present in CC024 mice and contributed to their MRSA tolerant phenotype. Two genes, C5ar1 and C5ar2, have high-impact variants in this region. C5ar1 and C5ar2 are receptors for the complement factor C5a, an anaphylatoxin that can trigger a massive immune response by binding to these receptors. We hypothesize that C5a may have altered binding to variant receptors in CC024 mice, reducing damage caused by the cytokine storm and resulting in the ability to tolerate a higher pathogen burden and longer survival.

A mutation in Themis contributes to anaphylaxis severity following oral peanut challenge in CC027 mice

BACKGROUND

The development of peanut allergy is due to a combination of genetic and environmental factors, although specific genes have proven difficult to identify. Previously, we reported that peanut-sensitized Collaborative Cross strain CC027/GeniUnc (CC027) mice develop anaphylaxis upon oral challenge to peanut, in contrast to C3H/HeJ (C3H) mice.

OBJECTIVE

This study aimed to determine the genetic basis of orally induced anaphylaxis to peanut in CC027 mice.

METHODS

A genetic mapping population between CC027 and C3H mice was designed to identify the genetic factors that drive oral anaphylaxis. A total of 356 CC027xC3H backcrossed mice were generated, sensitized to peanut, then challenged to peanut by oral gavage. Anaphylaxis and peanut-specific IgE were quantified for all mice. T-cell phenotyping was conducted on CC027 mice and 5 additional Collaborative Cross strains.

RESULTS

Anaphylaxis to peanut was absent in 77% of backcrossed mice, with 19% showing moderate anaphylaxis and 4% having severe anaphylaxis. There were 8 genetic loci associated with variation in response to peanut challenge-6 associated with anaphylaxis (temperature decrease) and 2 associated with peanut-specific IgE levels. There were 2 major loci that impacted multiple aspects of the severity of acute anaphylaxis, at which the CC027 allele was associated with worse outcome. At one of these loci, CC027 has a private genetic variant in the Themis gene. Consistent with described functions of Themis, we found that CC027 mice have more immature T cells with fewer CD8+, CD4+, and CD4+CD25+CD127- regulatory T cells.

CONCLUSIONS

Our results demonstrate a key role for Themis in the orally reactive CC027 mouse model of peanut allergy.

Improving laboratory animal genetic reporting: LAG-R guidelines

The biomedical research community addresses reproducibility challenges in animal studies through standardized nomenclature, improved experimental design, transparent reporting, data sharing, and centralized repositories. The ARRIVE guidelines outline documentation standards for laboratory animals in experiments, but genetic information is often incomplete. To remedy this, we propose the Laboratory Animal Genetic Reporting (LAG-R) framework. LAG-R aims to document animals' genetic makeup in scientific publications, providing essential details for replication and appropriate model use. While verifying complete genetic compositions may be impractical, better reporting and validation efforts enhance reliability of research. LAG-R standardization will bolster reproducibility, peer review, and overall scientific rigor.

ATP-dependent citrate lyase Drives Left Ventricular Dysfunction by Metabolic Remodeling of the Heart

Background

Metabolic remodeling is a hallmark of the failing heart. Oncometabolic stress during cancer increases the activity and abundance of the ATP-dependent citrate lyase (ACL, Acly ), which promotes histone acetylation and cardiac adaptation. ACL is critical for the de novo synthesis of lipids, but how these metabolic alterations contribute to cardiac structural and functional changes remains unclear.

Methods

We utilized human heart tissue samples from healthy donor hearts and patients with hypertrophic cardiomyopathy. Further, we used CRISPR/Cas9 gene editing to inactivate Acly in cardiomyocytes of MyH6-Cas9 mice. In vivo, positron emission tomography and ex vivo stable isotope tracer labeling were used to quantify metabolic flux changes in response to the loss of ACL. We conducted a multi-omics analysis using RNA-sequencing and mass spectrometry-based metabolomics and proteomics. Experimental data were integrated into computational modeling using the metabolic network CardioNet to identify significantly dysregulated metabolic processes at a systems level.

Results

Here, we show that in mice, ACL drives metabolic adaptation in the heart to sustain contractile function, histone acetylation, and lipid modulation. Notably, we show that loss of ACL increases glucose oxidation while maintaining fatty acid oxidation. Ex vivo isotope tracing experiments revealed a reduced efflux of glucose-derived citrate from the mitochondria into the cytosol, confirming that citrate is required for reductive metabolism in the heart. We demonstrate that YAP inactivation facilitates ACL deficiency. Computational flux analysis and integrative multi-omics analysis indicate that loss of ACL induces alternative isocitrate dehydrogenase 1 flux to compensate.

Conclusions

This study mechanistically delineates how cardiac metabolism compensates for suppressed citrate metabolism in response to ACL loss and uncovers metabolic vulnerabilities in the heart.

Mapping of susceptibility loci for Ebola virus pathogenesis in mice

Ebola virus (EBOV), a major global health concern, causes severe, often fatal EBOV disease (EVD) in humans. Host genetic variation plays a critical role, yet the identity of host susceptibility loci in mammals remains unknown. Using genetic reference populations, we generate an F2 mapping cohort to identify host susceptibility loci that regulate EVD. While disease-resistant mice display minimal pathogenesis, susceptible mice display severe liver pathology consistent with EVD-like disease and transcriptional signatures associated with inflammatory and liver metabolic processes. A significant quantitative trait locus (QTL) for virus RNA load in blood is identified in chromosome (chr)8, and a severe clinical disease and mortality QTL is mapped to chr7, which includes the Trim5 locus. Using knockout mice, we validate the Trim5 locus as one potential driver of liver failure and mortality after infection. The identification of susceptibility loci provides insight into molecular genetic mechanisms regulating EVD progression and severity, potentially informing therapeutics and vaccination strategies.

Complementary and Inducible creERT2 Mouse Models for Functional Evaluation of Endothelial Cell Subtypes in the Bone Marrow

In the adult bone marrow (BM), endothelial cells (ECs) are an integral component of the hematopoietic stem cell (HSC)-supportive niche, which modulates HSC activity by producing secreted and membrane-bound paracrine signals. Within the BM, distinct vascular arteriole, transitional, and sinusoidal EC subtypes display unique paracrine expression profiles and create anatomically-discrete microenvironments. However, the relative contributions of vascular endothelial subtypes in supporting hematopoiesis is unclear. Moreover, constitutive expression and off-target activity of currently available endothelial-specific and endothelial-subtype-specific murine cre lines potentially confound data analysis and interpretation. To address this, we describe two tamoxifen-inducible cre-expressing lines, Vegfr3-creERT2 and Cx40-creERT2, that efficiently label sinusoidal/transitional and arteriole endothelium respectively in adult marrow, without off-target activity in hematopoietic or perivascular cells. Utilizing an established mouse model in which cre-dependent recombination constitutively-activates MAPK signaling within adult endothelium, we identify arteriole ECs as the driver of MAPK-mediated hematopoietic dysfunction. These results define complementary tamoxifen-inducible creERT2-expressing mouse lines that label functionally-discrete and non-overlapping sinusoidal/transitional and arteriole EC populations in the adult BM, providing a robust toolset to investigate the differential contributions of vascular subtypes in maintaining hematopoietic homeostasis.

Adam19 Deficiency Impacts Pulmonary Function: Human GWAS Follow-up in Mouse

Purpose

Over 550 loci have been associated with human pulmonary function in genome-wide association studies (GWAS); however, the causal role of most remains uncertain. Single nucleotide polymorphisms in a disintegrin and metalloprotease domain 19 (ADAM19) are consistently related to pulmonary function in GWAS. Thus, we used a mouse model to investigate the causal link between Adam19 and pulmonary function.

Methods

We created an Adam19 knockout (KO) mouse model and validated the gene targeting using RNA-Seq and RT-qPCR. Contrary to prior publications, the KO was not neonatal lethal. Thus, we phenotyped the Adam19 KO.

Results

KO mice had lower body weight and shorter tibial length than wild type (WT). Dual-energy X-ray Absorptiometry indicated lower soft weight, fat weight, and bone mineral content in KO mice. In lung function analyses using flexiVent, compared to WT, Adam19 KO had decreased baseline respiratory system elastance, minute work of breathing, tissue damping, tissue elastance, and forced expiratory flow at 50% forced vital capacity but higher FEV0.1 and FVC. Adam19 KO had attenuated tissue damping and tissue elastance in response to methacholine following LPS exposure. Adam19 KO also exhibited attenuated neutrophil extravasation into the airway after LPS administration compared to WT. RNA-Seq analysis of KO and WT lungs identified several differentially expressed genes (Cd300lg, Kpna2, and Pttg1) implicated in lung biology and pathogenesis. Gene set enrichment analysis identified negative enrichment for TNF pathways.

Conclusion

Our murine findings support a causal role of ADAM19, implicated in human GWAS, in regulating pulmonary function.

Contributions of mouse genetic strain background to age-related phenotypes in physically active HET3 mice

We assessed aging hallmarks in skin, muscle, and adipose in the genetically diverse HET3 mouse, and generated a broad dataset comparing these to individual animal diagnostic SNPs from the 4 founding inbred strains of the HET3 line. For middle- and old-aged HET3 mice, we provided running wheel exercise to ensure our observations were not purely representative of sedentary animals, but age-related phenotypes were not improved with running wheel activity. Adipose tissue fibrosis, peripheral neuropathy, and loss of neuromuscular junction integrity were consistent phenotypes in older-aged HET3 mice regardless of physical activity, but aspects of these phenotypes were moderated by the SNP% contributions of the founding strains for the HET3 line. Taken together, the genetic contribution of founder strain SNPs moderated age-related phenotypes in skin and muscle innervation and were dependent on biological sex and chronological age. However, there was not a single founder strain (BALB/cJ, C57BL/6J, C3H/HeJ, DBA/2J) that appeared to drive more protection or disease-risk across aging in this mouse line, but genetic diversity in general was more protective.

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.

Host Genetic Variation Impacts SARS-CoV-2 Vaccination Response in the Diversity Outbred Mouse Population

The COVID-19 pandemic led to the rapid and worldwide development of highly effective vaccines against SARS-CoV-2. However, there is significant individual-to-individual variation in vaccine efficacy due to factors including viral variants, host age, immune status, environmental and host genetic factors. Understanding those determinants driving this variation may inform the development of more broadly protective vaccine strategies. While host genetic factors are known to impact vaccine efficacy for respiratory pathogens such as influenza and tuberculosis, the impact of host genetic variation on vaccine efficacy against COVID-19 is not well understood. To model the impact of host genetic variation on SARS-CoV-2 vaccine efficacy, while controlling for the impact of non-genetic factors, we used the Diversity Outbred (DO) mouse model. We found that DO mice immunized against SARS-CoV-2 exhibited high levels of variation in vaccine-induced neutralizing antibody responses. While the majority of the vaccinated mice were protected from virus-induced disease, similar to human populations, we observed vaccine breakthrough in a subset of mice. Importantly, we found that this variation in neutralizing antibody, virus-induced disease, and viral titer is heritable, indicating that the DO serves as a useful model system for studying the contribution of genetic variation of both vaccines and disease outcomes.

Large differences in collateral blood vessel abundance among individuals arise from multiple genetic variants

Collateral blood flow varies greatly among humans for reasons that remain unclear, resulting in significant differences in ischemic tissue damage. A similarly large variation has also been found in mice that is caused by genetic background-dependent differences in the extent of collateral formation, termed collaterogenesis-a unique angiogenic process that occurs during development and determines collateral number and diameter in the adult. Previous studies have identified several quantitative trait loci (QTL) linked to this variation. However, understanding has been hampered by the use of closely related inbred strains that do not model the wide genetic variation present in the "outbred" human population. The Collaborative Cross (CC) multiparent mouse genetic reference panel was developed to address this limitation. Herein we measured the number and average diameter of cerebral collaterals in 60 CC strains, their 8 founder strains, 8 F1 crosses of CC strains selected for abundant versus sparse collaterals, and 2 intercross populations created from the latter. Collateral number evidenced 47-fold variation among the 60 CC strains, with 14% having poor, 25% poor-to-intermediate, 47% intermediate-to-good, and 13% good collateral abundance, that was associated with large differences in post-stroke infarct volume. Collateral number in skeletal muscle and intestine of selected high- and low-collateral strains evidenced the same relative abundance as in brain. Genome-wide mapping demonstrated that collateral abundance is a highly polymorphic trait. Subsequent analysis identified: 6 novel QTL circumscribing 28 high-priority candidate genes harboring putative loss-of-function polymorphisms (SNPs) associated with low collateral number; 335 predicted-deleterious SNPs present in their human orthologs; and 32 genes associated with vascular development but lacking protein coding variants. Six additional suggestive QTL (LOD > 4.5) were also identified in CC-wide QTL mapping. This study provides a comprehensive set of candidate genes for future investigations aimed at identifying signaling proteins within the collaterogenesis pathway whose variants potentially underlie genetic-dependent collateral insufficiency in brain and other tissues.

A mutation in Themis contributes to peanut-induced oral anaphylaxis in CC027 mice

Background

The development of peanut allergy is due to a combination of genetic and environmental factors, although specific genes have proven difficult to identify. Previously, we reported that peanut-sensitized CC027/GeniUnc (CC027) mice develop anaphylaxis upon oral challenge to peanut, unlike C3H/HeJ (C3H) mice.

Objective

To determine the genetic basis of orally-induced anaphylaxis to peanut in CC027 mice.

Methods

A genetic mapping population between CC027 and C3H mice was designed to identify the genetic factors that drive oral anaphylaxis. A total of 356 CC027xC3H backcrossed mice were generated, sensitized to peanut, then challenged to peanut by oral gavage. Anaphylaxis and peanut-specific IgE were quantified for all mice. T-cell phenotyping was conducted on CC027 and five additional CC strains.

Results

Anaphylaxis to peanut was absent in 77% of backcrossed mice, with 19% showing moderate anaphylaxis, and 4% having severe anaphylaxis. A total of eight genetic loci were associated with variation in response to peanut challenge, six associated with anaphylaxis (temperature decrease) and two associated with peanut-specific IgE levels. There were two major loci that impacted multiple aspects of the severity of acute anaphylaxis, at which the CC027 allele was associated with worse outcome. At one of these loci, CC027 has a private genetic variant in the Themis (thymocyte-expressed molecule involved in selection) gene. Consistent with Themis' described functions, we found that CC027 have more immature T cells with fewer CD8+, CD4+, and CD4+CD25+CD127- regulatory T cells.

Conclusion

Our results demonstrate a key role for Themis in the orally-reactive CC027 mouse model of peanut allergy.

Pinna Edge Biopsy of 7 and 21 Day Old C57BL/6 Mice as a Method for Identification and Genotyping

Identifying and genotyping mice prior to weaning can be useful for mouse colony management. Mice of an undesired genotype can be identified prior to weaning and removed from further study, resulting in a reduction of housing costs, and labor time. We hypothesized that a pinna edge biopsy (PEB) performed by removing a portion of its edge with scissors is a reliable method for identifying and genotyping mice on postnatal day (PND) 7 consistent with PND 21, weaned mice. The pinnae of 54 C57BL/NCrl6 mice were biopsied on PND 7, and another 54 were biopsied on PND 21. Nine pinna patterns were tested. The accuracy of pattern identification was assessed on PND 7, 14, 21, 30, and 63. The mean times were compared for performing the biopsy on PND 7 and PND 21 mice, and the average time taken to identify the patterns were determined. Weight, milk spot presence, pup rejection, morbidity, and mortality were examined at various time points. During the biopsy, bleeding of the pinna, urination, vocalization, and flinching were assessed. No significant differences were detected in DNA quality, relative DNA quantity, genotyping reliability, or body weight (P ≥ 0.05) between mice biopsied on PND 7 and PND 21. Flinching at the time of PEB was significantly higher in PND 21 mice as compared with PND 7 mice (P < 0.00001). Pinna pattern identification accuracy for mice biopsied on PND 7 and PND 21 were 96% and 98%, respectively. This study validates the use of PEB for simultaneous identification and genotyping of PND 7 mice.

A single dominant locus restricts retrovirus replication in YBR/Ei mice

Differential responses to viral infections are influenced by the genetic makeup of the host. Studies of resistance to retroviruses in human populations are complicated due to the inability to conduct proof-of-principle studies. Inbred mouse lines, which have a range of susceptible phenotypes to retroviruses, are an ideal tool to identify and characterize mechanisms of resistance and define their genetic underpinnings. YBR/Ei mice become infected with Mouse Mammary Tumor Virus, a mucosally transmitted murine retrovirus, but eliminate the virus from their pedigrees. Virus elimination correlates with a lack of virus-specific neonatal oral tolerance, which is a major mechanism for blocking the anti-virus response in susceptible mice. Virus control is unrelated to virus-neutralizing antibodies, cytotoxic CD8+ T cells, NK cells, and NK T cells, which are the best characterized mechanisms of resistance to retroviruses. We identified a single, dominant locus that controls the resistance mechanism, which we provisionally named attenuation of virus titers (Avt) and mapped to the distal region of chromosome 18. IMPORTANCE Elucidation of the mechanism that mediates resistance to retroviruses is of fundamental importance to human health, as it will ultimately lead to knowledge of the genetic differences among individuals in susceptibility to microbial infections.

Altered phenotypes due to genetic interaction between the mouse phosphoinositide biosynthesis genes Fig4 and Pip4k2c

Loss-of-function mutations of FIG4 are responsible for neurological disorders in human and mouse that result from reduced abundance of the signaling lipid PI(3,5)P2. In contrast, loss-of-function mutations of the phosphoinositide kinase PIP4K2C result in elevated abundance of PI(3,5)P2. These opposing effects on PI(3,5)P2 suggested that we might be able to compensate for deficiency of FIG4 by reducing expression of PIP4K2C. To test this hypothesis in a whole animal model, we generated triallelic mice with genotype Fig 4-/-, Pip4k2c+/-; these mice are null for Fig 4 and haploinsufficient for Pip4k2c. The neonatal lethality of Fig 4 null mice in the C57BL/6J strain background was rescued by reduced expression of Pip4k2c. The lysosome enlargement characteristic of Fig 4 null cells was also reduced by heterozygous loss of Pip4k2c. The data demonstrate interaction between these two genes, and suggest that inhibition of the kinase PIPK4C2 could be a target for treatment of FIG4 deficiency disorders such as Charcot-Marie-Tooth Type 4J and Yunis-Varón Syndrome.

Large differences in collateral blood vessel abundance among individuals arise from multiple genetic variants

Collateral blood flow varies greatly among humans for reasons that remain unclear, resulting in significant differences in ischemic tissue damage. A similarly large variation has also been found in mice that is caused by genetic background-dependent differences in the extent of collateral formation, termed collaterogenesis-a unique angiogenic process that occurs during development and determines collateral number and diameter in the adult. Previous studies have identified several quantitative trait loci (QTL) linked to this variation. However, understanding has been hampered by the use of closely related inbred strains that do not model the wide genetic variation present in the "outbred" human population. The Collaborative Cross (CC) multiparent mouse genetic reference panel was developed to address this limitation. Herein we measured the number and average diameter of cerebral collaterals in 60 CC strains, their 8 founder strains, 8 F1 crosses of CC strains selected for abundant versus sparse collaterals, and 2 intercross populations created from the latter. Collateral number evidenced 47-fold variation among the 60 CC strains, with 14% having poor, 25% poor-to-intermediate, 47% intermediate-to-good, and 13% good collateral abundance, that was associated with large differences in post-stroke infarct volume. Genome-wide mapping demonstrated that collateral abundance is a highly polymorphic trait. Subsequent analysis identified: 6 novel QTL circumscribing 28 high-priority candidate genes harboring putative loss-of-function polymorphisms (SNPs) associated with low collateral number; 335 predicted-deleterious SNPs present in their human orthologs; and 32 genes associated with vascular development but lacking protein coding variants. This study provides a comprehensive set of candidate genes for future investigations aimed at identifying signaling proteins within the collaterogenesis pathway whose variants potentially underlie genetic-dependent collateral insufficiency in brain and other tissues.

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.

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

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

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.

Identification of collaborative cross mouse strains permissive to Salmonella enterica serovar Typhi infection

Salmonella enterica serovar Typhi is the causative agent of typhoid fever restricted to humans and does not replicate in commonly used inbred mice. Genetic variation in humans is far greater and more complex than that in a single inbred strain of mice. The Collaborative Cross (CC) is a large panel of recombinant inbred strains which has a wider range of genetic diversity than laboratory inbred mouse strains. We found that the CC003/Unc and CC053/Unc strains are permissive to intraperitoneal but not oral route of S. Typhi infection and show histopathological changes characteristic of human typhoid. These CC strains are immunocompetent, and immunization induces antigen-specific responses that can kill S. Typhi in vitro and control S. Typhi in vivo. Our results indicate that CC003/Unc and CC053/Unc strains can help identify the genetic basis for typhoid susceptibility, S. Typhi virulence mechanism(s) in vivo, and serve as a preclinical mammalian model system to identify effective vaccines and therapeutics strategies.

Genetic and Molecular Quality Control of Genetically Engineered Mice

Genetically engineered mice are used as avatars to understand mammalian gene function and develop therapies for human disease. During genetic modification, unintended changes can occur, and these changes may result in misassigned gene-phenotype relationships leading to incorrect or incomplete experimental interpretations. The types of unintended changes that may occur depend on the allele type being made and the genetic engineering approach used. Here we broadly categorize allele types as deletions, insertions, base changes, and transgenes derived from engineered embryonic stem (ES) cells or edited mouse embryos. However, the methods we describe can be adapted to other allele types and engineering strategies. We describe the sources and consequ ences of common unintended changes and best practices for detecting both intended and unintended changes by screening and genetic and molecular quality control (QC) of chimeras, founders, and their progeny. Employing these practices, along with careful allele design and good colony management, will increase the chance that investigations using genetically engineered mice will produce high-quality reproducible results, to enable a robust understanding of gene function, human disease etiology, and therapeutic development.

Genetic quality: a complex issue for experimental study reproducibility

Laboratory animal research involving mice, requires consideration of many factors to be controlled. Genetic quality is one factor that is often overlooked but is essential for the generation of reproducible experimental results. Whether experimental research involves inbred mice, spontaneous mutant, or genetically modified strains, exercising genetic quality through careful breeding, good recordkeeping, and prudent quality control steps such as validation of the presence of mutations and verification of the genetic background, will help ensure that experimental results are accurate and that reference controls are representative for the particular experiment. In this review paper, we will discuss various techniques used for the generation of genetically altered mice, and the different aspects to be considered regarding genetic quality, including inbred strains and substrains used, quality check controls during and after genetic manipulation and breeding. We also provide examples for when to use the different techniques and considerations on genetic quality checks. Further, we emphasize on the importance of establishing an in-house genetic quality program.

A Multitrait Locus Regulates Sarbecovirus Pathogenesis

Infectious diseases have shaped the human population genetic structure, and genetic variation influences the susceptibility to many viral diseases. However, a variety of challenges have made the implementation of traditional human Genome-wide Association Studies (GWAS) approaches to study these infectious outcomes challenging. In contrast, mouse models of infectious diseases provide an experimental control and precision, which facilitates analyses and mechanistic studies of the role of genetic variation on infection. Here we use a genetic mapping cross between two distinct Collaborative Cross mouse strains with respect to severe acute respiratory syndrome coronavirus (SARS-CoV) disease outcomes. We find several loci control differential disease outcome for a variety of traits in the context of SARS-CoV infection. Importantly, we identify a locus on mouse chromosome 9 that shows conserved synteny with a human GWAS locus for SARS-CoV-2 severe disease. We follow-up and confirm a role for this locus, and identify two candidate genes, CCR9 and CXCR6, that both play a key role in regulating the severity of SARS-CoV, SARS-CoV-2, and a distantly related bat sarbecovirus disease outcomes. As such we provide a template for using experimental mouse crosses to identify and characterize multitrait loci that regulate pathogenic infectious outcomes across species. IMPORTANCE Host genetic variation is an important determinant that predicts disease outcomes following infection. In the setting of highly pathogenic coronavirus infections genetic determinants underlying host susceptibility and mortality remain unclear. To elucidate the role of host genetic variation on sarbecovirus pathogenesis and disease outcomes, we utilized the Collaborative Cross (CC) mouse genetic reference population as a model to identify susceptibility alleles to SARS-CoV and SARS-CoV-2 infections. Our findings reveal that a multitrait loci found in chromosome 9 is an important regulator of sarbecovirus pathogenesis in mice. Within this locus, we identified and validated CCR9 and CXCR6 as important regulators of host disease outcomes. Specifically, both CCR9 and CXCR6 are protective against severe SARS-CoV, SARS-CoV-2, and SARS-related HKU3 virus disease in mice. This chromosome 9 multitrait locus may be important to help identify genes that regulate coronavirus disease outcomes in humans.

stuart: an R package for the curation of SNP genotypes from experimental crosses

Genetic mapping in 2-generation crosses requires genotyping, usually performed with single nucleotide polymorphism markers arrays which provide high-density genetic information. However, genetic analysis on raw genotypes can lead to spurious or unreliable results due to defective single nucleotide polymorphism assays or wrong genotype interpretation. Here, we introduce stuart, an open-source R package, which analyzes raw genotyping data to filter single nucleotide polymorphism markers based on informativeness, Mendelian inheritance pattern, and consistency with parental genotypes. The functions of this package provide a curation pipeline and formatting adequate for genetic analysis with the R/qtl package. stuart is available with detailed documentation from https://gitlab.pasteur.fr/mouselab/stuart/.

Genic and Chromosomal Components of Prdm9-Driven Hybrid Male Sterility in Mice (Mus musculus)

Hybrid sterility contributes to speciation by preventing gene flow between related taxa. Prdm9, the first and only hybrid male sterility (HMS) gene known in vertebrates, predetermines the sites of recombination between homologous chromosomes and their synapsis in early meiotic prophase. The asymmetric binding of PRDM9 to heterosubspecific homologs of Mus m. musculus x Mus m. domesticus F1 hybrids and increase of PRDM9-independent DNA double-strand break (DSB) hotspots results in difficult to repair DSBs, incomplete synapsis of homologous chromosomes and meiotic arrest at the first meiotic prophase. Here we show that Prdm9 behaves as a major HMS gene in mice outside the Mus m. musculus x Mus m. domesticus F1 hybrids, in the genomes composed of Mus m. castaneus and Mus m. musculus chromosomes segregating on the Mus m. domesticus background. The Prdm9cst/dom2 (castaneus/domesticus) allelic combination secures meiotic synapsis, testes weight and sperm count within physiological limits, while the Prdm9msc1/dom2 (musculus/domesticus) males show a range of fertility impairment. Out of five quantitative trait loci contributing to the Prdm9msc1/dom2-related infertility, four control either meiotic synapsis or fertility phenotypes and one controls both, synapsis and fertility. Whole-genome genotyping of individual chromosomes showed preferential involvement of nonrecombinant musculus chromosomes in asynapsis in accordance with the chromosomal character of HMS. Moreover, we show that the overall asynapsis rate can be estimated solely from the genotype of individual males by scoring the effect of nonrecombinant musculus chromosomes. Prdm9-controlled HMS represents an example of genetic architecture of HMS consisting of genic and chromosomal components.

Zhx2 Is a Candidate Gene Underlying Oxymorphone Metabolite Brain Concentration Associated with State-Dependent Oxycodone Reward

Understanding the pharmacogenomics of opioid metabolism and behavior is vital to therapeutic success, as mutations can dramatically alter therapeutic efficacy and addiction liability. We found robust, sex-dependent BALB/c substrain differences in oxycodone behaviors and whole brain concentration of oxycodone metabolites. BALB/cJ females showed robust state-dependent oxycodone reward learning as measured via conditioned place preference when compared with the closely related BALB/cByJ substrain. Accordingly, BALB/cJ females also showed a robust increase in brain concentration of the inactive metabolite noroxycodone and the active metabolite oxymorphone compared with BALB/cByJ mice. Oxymorphone is a highly potent, full agonist at the mu opioid receptor that could enhance drug-induced interoception and state-dependent oxycodone reward learning. Quantitative trait locus (QTL) mapping in a BALB/c F2 reduced complexity cross revealed one major QTL on chromosome 15 underlying brain oxymorphone concentration that explained 32% of the female variance. BALB/cJ and BALB/cByJ differ by fewer than 10,000 variants, which can greatly facilitate candidate gene/variant identification. Hippocampal and striatal cis-expression QTL (eQTL) and exon-level eQTL analysis identified Zhx2, a candidate gene coding for a transcriptional repressor with a private BALB/cJ retroviral insertion that reduces Zhx2 expression and sex-dependent dysregulation of cytochrome P450 enzymes. Whole brain proteomics corroborated the Zhx2 eQTL and identified upregulated CYP2D11 that could increase brain oxymorphone in BALB/cJ females. To summarize, Zhx2 is a highly promising candidate gene underlying brain oxycodone metabolite levels. Future studies will validate Zhx2 and its site of action using reciprocal gene editing and tissue-specific viral manipulations in BALB/c substrains. SIGNIFICANCE STATEMENT: Our findings show that genetic variation can result in sex-specific alterations in whole brain concentration of a bioactive opioid metabolite after oxycodone administration, reinforcing the need for sex as a biological factor in pharmacogenomic studies. The cooccurrence of female-specific increased oxymorphone and state-dependent reward learning suggests that this minor yet potent and efficacious metabolite of oxycodone could increase opioid interoception and drug-cue associative learning of opioid reward, which has implications for cue-induced relapse of drug-seeking behavior and for precision pharmacogenetics.

The lupus susceptibility allele DRB1*03:01 encodes a disease-driving epitope

The HLA-DRB1*03:01 allele is a major genetic risk factor in systemic lupus erythematosus (SLE), but the mechanistic basis of the association is unclear. Here we show that in the presence of interferon gamma (IFN-γ), a short DRB1*03:01-encoded allelic epitope activates a characteristic lupus transcriptome in mouse and human macrophages. It also triggers a cascade of SLE-associated cellular aberrations, including endoplasmic reticulum stress, unfolded protein response, mitochondrial dysfunction, necroptotic cell death, and production of pro-inflammatory cytokines. Parenteral administration of IFN-γ to naïve DRB1*03:01 transgenic mice causes increased serum levels of anti-double stranded DNA antibodies, glomerular immune complex deposition and histopathological renal changes that resemble human lupus nephritis. This study provides evidence for a noncanonical, antigen presentation-independent mechanism of HLA-disease association in SLE and could lay new foundations for our understanding of key molecular mechanisms that trigger and propagate this devastating autoimmune disease.

Extracellular superoxide dismutase (EC-SOD) R213G variant reduces mitochondrial ROS and preserves mitochondrial function in bleomycin-induced lung injury: EC-SOD R213G variant and intracellular redox regulation

Extracellular superoxide dismutase (EC-SOD) is highly expressed in the lung and vasculature. A common human single nucleotide polymorphism (SNP) in the matrix binding region of EC-SOD leads to a single amino acid substitution, R213G, and alters EC-SOD tissue binding affinity. The change in tissue binding affinity redistributes EC-SOD from tissue to extracellular fluids. Mice (R213G mice) expressing a knock-in of this EC-SOD SNP exhibit elevated plasma and reduced lung EC-SOD content and activity and are protected against bleomycin-induced lung injury and inflammation. It is unknown how the redistribution of EC-SOD alters site-specific redox-regulated molecules relevant for protection. In this study, we tested the hypothesis that the change in the local EC-SOD content would influence not only the extracellular redox microenvironment where EC-SOD is localized but also protect the intracellular redox status of the lung. Mice were treated with bleomycin and harvested 7 days post-treatment. Superoxide levels, measured by electron paramagnetic resonance (EPR), were lower in plasma and Bronchoalveolar lavage fluid (BALF) cells in R213G mice compared to wild-type (WT) mice, while lung cellular superoxide levels in R213G mice were not elevated post-bleomycin compared to WT mice despite low lung EC-SOD levels. Lung glutathione redox potential (EhGSSG), determined by HPLC and fluorescence, was more oxidized in WT compared to R213G mice. In R213G mice, lung mitochondrial oxidative stress was reduced shown by mitochondrial superoxide level measured by EPR in lung and the resistance to bleomycin-induced cardiolipin oxidation. Bleomycin treatment suppressed mitochondrial respiration in WT mice. Mitochondrial function was impaired at baseline in R213G mice but did not exhibit further suppression in respiration post-bleomycin. Collectively, the results indicate that R213G variant preserves intracellular redox state and protects mitochondrial function in the setting of bleomycin-induced inflammation.

QTL-mapping in the obese Berlin Fat Mouse identifies additional candidate genes for obesity and fatty liver disease

The Berlin Fat Mouse Inbred line (BFMI) is a model for obesity and the metabolic syndrome. This study aimed to identify genetic variants associated with liver weight, liver triglycerides, and body weight using the obese BFMI sub-line BFMI861-S1. BFMI861-S1 mice are insulin resistant and store ectopic fat in the liver. In generation 10, 58 males and 65 females of the advanced intercross line (AIL) BFMI861-S1xB6N were phenotyped under a standard diet over 20 weeks. QTL analysis was performed after genotyping with the MiniMUGA Genotyping Array. Whole-genome sequencing and gene expression data of the parental lines was used for the prioritization of positional candidate genes. Three QTLs associated with liver weight, body weight, and subcutaneous adipose tissue (scAT) weight were identified. A highly significant QTL on chromosome (Chr) 1 (157–168 Mb) showed an association with liver weight. A QTL for body weight at 20 weeks was found on Chr 3 (34.1–40 Mb) overlapping with a QTL for scAT weight. In a multiple QTL mapping approach, an additional QTL affecting body weight at 16 weeks was identified on Chr 6 (9.5–26.1 Mb). Considering sequence variants and expression differences, Sec16b and Astn1 were prioritized as top positional candidate genes for the liver weight QTL on Chr 1; Met and Ica1 for the body weight QTL on Chr 6. Interestingly, all top candidate genes have previously been linked with metabolic traits. This study shows once more the power of an advanced intercross line for fine mapping. QTL mapping combined with a detailed prioritization approach allowed us to identify additional and plausible candidate genes linked to metabolic traits in the BFMI861-S1xB6N AIL. By reidentifying known candidate genes in a different crossing population the causal link with specific traits is underlined and additional evidence is given for further investigations.

Selective retinal ganglion cell loss and optic neuropathy in a humanized mouse model of familial dysautonomia

Familial dysautonomia (FD) is an autosomal recessive neurodegenerative disease caused by a splicing mutation in the gene encoding Elongator complex protein 1 (ELP1, also known as IKBKAP). This mutation results in tissue-specific skipping of exon 20 with a corresponding reduction of ELP1 protein, predominantly in the central and peripheral nervous system. Although FD patients have a complex neurological phenotype caused by continuous depletion of sensory and autonomic neurons, progressive visual decline leading to blindness is one of the most problematic aspects of the disease, as it severely affects their quality of life. To better understand the disease mechanism as well as to test the in vivo efficacy of targeted therapies for FD, we have recently generated a novel phenotypic mouse model, TgFD9; IkbkapΔ20/flox. This mouse exhibits most of the clinical features of the disease and accurately recapitulates the tissue-specific splicing defect observed in FD patients. Driven by the dire need to develop therapies targeting retinal degeneration in FD, herein, we comprehensively characterized the progression of the retinal phenotype in this mouse, and we demonstrated that it is possible to correct ELP1 splicing defect in the retina using the splicing modulator compound (SMC) BPN-15477.

A Multitrait Locus Regulates Sarbecovirus Pathogenesis

Infectious diseases have shaped the human population genetic structure, and genetic variation influences the susceptibility to many viral diseases. However, a variety of challenges have made the implementation of traditional human Genome-wide Association Studies (GWAS) approaches to study these infectious outcomes challenging. In contrast, mouse models of infectious diseases provide an experimental control and precision, which facilitates analyses and mechanistic studies of the role of genetic variation on infection. Here we use a genetic mapping cross between two distinct Collaborative Cross mouse strains with respect to SARS-CoV disease outcomes. We find several loci control differential disease outcome for a variety of traits in the context of SARS-CoV infection. Importantly, we identify a locus on mouse Chromosome 9 that shows conserved synteny with a human GWAS locus for SARS-CoV-2 severe disease. We follow-up and confirm a role for this locus, and identify two candidate genes, CCR9 and CXCR6 that both play a key role in regulating the severity of SARS-CoV, SARS-CoV-2 and a distantly related bat sarbecovirus disease outcomes. As such we provide a template for using experimental mouse crosses to identify and characterize multitrait loci that regulate pathogenic infectious outcomes across species.

Cocaine-Induced Locomotor Activation Differs Across Inbred Mouse Substrains

Cocaine use disorders (CUD) are devastating for affected individuals and impose a significant societal burden, but there are currently no FDA-approved therapies. The development of novel and effective treatments has been hindered by substantial gaps in our knowledge about the etiology of these disorders. The risk for developing a CUD is influenced by genetics, the environment and complex interactions between the two. Identifying specific genes and environmental risk factors that increase CUD risk would provide an avenue for the development of novel treatments. Rodent models of addiction-relevant behaviors have been a valuable tool for studying the genetics of behavioral responses to drugs of abuse. Traditional genetic mapping using genetically and phenotypically divergent inbred mice has been successful in identifying numerous chromosomal regions that influence addiction-relevant behaviors, but these strategies rarely result in identification of the causal gene or genetic variant. To overcome this challenge, reduced complexity crosses (RCC) between closely related inbred mouse strains have been proposed as a method for rapidly identifying and validating functional variants. The RCC approach is dependent on identifying phenotypic differences between substrains. To date, however, the study of addiction-relevant behaviors has been limited to very few sets of substrains, mostly comprising the C57BL/6 lineage. The present study expands upon the current literature to assess cocaine-induced locomotor activation in 20 inbred mouse substrains representing six inbred strain lineages (A/J, BALB/c, FVB/N, C3H/He, DBA/2 and NOD) that were either bred in-house or supplied directly by a commercial vendor. To our knowledge, we are the first to identify significant differences in cocaine-induced locomotor response in several of these inbred substrains. The identification of substrain differences allows for the initiation of RCC populations to more rapidly identify specific genetic variants associated with acute cocaine response. The observation of behavioral profiles that differ between mice generated in-house and those that are vendor-supplied also presents an opportunity to investigate the influence of environmental factors on cocaine-induced locomotor activity.

Pierre CL, Williams A, Sohni A, Wilkinson MF, Gymrek M, Sebat J, Palmer AA. SNPs, short tandem repeats, and structural variants are responsible for differential gene expression across C57BL/6 and C57BL/10 substrains

Mouse substrains are an invaluable model for understanding disease. We compared C57BL/6J, which is the most commonly used inbred mouse strain, with eight C57BL/6 and five C57BL/10 closely related inbred substrains. Whole-genome sequencing and RNA-sequencing analysis yielded 352,631 SNPs, 109,096 indels, 150,344 short tandem repeats (STRs), 3,425 structural variants (SVs), and 2,826 differentially expressed genes (DE genes) among these 14 strains; 312,981 SNPs (89%) distinguished the B6 and B10 lineages. These SNPs were clustered into 28 short segments that are likely due to introgressed haplotypes rather than new mutations. Outside of these introgressed regions, we identified 53 SVs, protein-truncating SNPs, and frameshifting indels that were associated with DE genes. Our results can be used for both forward and reverse genetic approaches and illustrate how introgression and mutational processes give rise to differences among these widely used inbred substrains.

Natural rodent model of viral transmission reveals biological features of virus population dynamics

Emerging viruses threaten global health, but few experimental models can characterize the virus and host factors necessary for within- and cross-species transmission. Here, we leverage a model whereby pet store mice or rats-which harbor natural rodent pathogens-are cohoused with laboratory mice. This "dirty" mouse model offers a platform for studying acute transmission of viruses between and within hosts via natural mechanisms. We identified numerous viruses and other microbial species that transmit to cohoused mice, including prospective new members of the Coronaviridae, Astroviridae, Picornaviridae, and Narnaviridae families, and uncovered pathogen interactions that promote or prevent virus transmission. We also evaluated transmission dynamics of murine astroviruses during transmission and spread within a new host. Finally, by cohousing our laboratory mice with the bedding of pet store rats, we identified cross-species transmission of a rat astrovirus. Overall, this model system allows for the analysis of transmission of natural rodent viruses and is a platform to further characterize barriers to zoonosis.

Host-pathogen genetic interactions underlie tuberculosis susceptibility in genetically diverse mice

The outcome of an encounter with Mycobacterium tuberculosis (Mtb) depends on the pathogen’s ability to adapt to the variable immune pressures exerted by the host. Understanding this interplay has proven difficult, largely because experimentally tractable animal models do not recapitulate the heterogeneity of tuberculosis disease. We leveraged the genetically diverse Collaborative Cross (CC) mouse panel in conjunction with a library of Mtb mutants to create a resource for associating bacterial genetic requirements with host genetics and immunity. We report that CC strains vary dramatically in their susceptibility to infection and produce qualitatively distinct immune states. Global analysis of Mtb transposon mutant fitness (TnSeq) across the CC panel revealed that many virulence pathways are only required in specific host microenvironments, identifying a large fraction of the pathogen’s genome that has been maintained to ensure fitness in a diverse population. Both immunological and bacterial traits can be associated with genetic variants distributed across the mouse genome, making the CC a unique population for identifying specific host-pathogen genetic interactions that influence pathogenesis.

The Mutant Mouse Resource and Research Center (MMRRC): the NIH-supported National Public Repository and Distribution Archive of Mutant Mouse Models in the USA

The Mutant Mouse Resource and Research Center (MMRRC) Program is the pre-eminent public national mutant mouse repository and distribution archive in the USA, serving as a national resource of mutant mice available to the global scientific community for biomedical research. Established more than two decades ago with grants from the National Institutes of Health (NIH), the MMRRC Program supports a Consortium of regionally distributed and dedicated vivaria, laboratories, and offices (Centers) and an Informatics Coordination and Service Center (ICSC) at three academic teaching and research universities and one non-profit genetic research institution. The MMRRC Program accepts the submission of unique, scientifically rigorous, and experimentally valuable genetically altered and other mouse models donated by academic and commercial scientists and organizations for deposition, maintenance, preservation, and dissemination to scientists upon request. The four Centers maintain an archive of nearly 60,000 mutant alleles as live mice, frozen germplasm, and/or embryonic stem (ES) cells. Since its inception, the Centers have fulfilled 13,184 orders for mutant mouse models from 9591 scientists at 6626 institutions around the globe. Centers also provide numerous services that facilitate using mutant mouse models obtained from the MMRRC, including genetic assays, microbiome analysis, analytical phenotyping and pathology, cryorecovery, mouse husbandry, infectious disease surveillance and diagnosis, and disease modeling. The ICSC coordinates activities between the Centers, manages the website (mmrrc.org) and online catalog, and conducts communication, outreach, and education to the research community. Centers preserve, secure, and protect mutant mouse lines in perpetuity, promote rigor and reproducibility in scientific experiments using mice, provide experiential training and consultation in the responsible use of mice in research, and pursue cutting edge technologies to advance biomedical studies using mice to improve human health. Researchers benefit from an expansive list of well-defined mouse models of disease that meet the highest standards of rigor and reproducibility, while donating investigators benefit by having their mouse lines preserved, protected, and distributed in compliance with NIH policies.

Genetic basis of thermal nociceptive sensitivity and brain weight in a BALB/c reduced complexity cross

Thermal nociception involves the transmission of temperature-related noxious information from the periphery to the CNS and is a heritable trait that could predict transition to persistent pain. Rodent forward genetics complement human studies by controlling genetic complexity and environmental factors, analysis of end point tissue, and validation of variants on appropriate genetic backgrounds. Reduced complexity crosses between nearly identical inbred substrains with robust trait differences can greatly facilitate unbiased discovery of novel genes and variants. We found BALB/cByJ mice showed enhanced sensitivity on the 53.5°C hot plate and mechanical stimulation in the von Frey test compared to BALB/cJ mice and replicated decreased gross brain weight in BALB/cByJ versus BALB/cJ. We then identified a quantitative trait locus (QTL) on chromosome 13 for hot plate sensitivity (LOD = 10.7; p < 0.001; peak = 56 Mb) and a QTL for brain weight on chromosome 5 (LOD = 8.7; p < 0.001). Expression QTL mapping of brain tissues identified H2afy (56.07 Mb) as the top transcript with the strongest association at the hot plate locus (FDR = 0.0002) and spliceome analysis identified differential exon usage within H2afy associated with the same locus. Whole brain proteomics further supported decreased H2AFY expression could underlie enhanced hot plate sensitivity, and identified ACADS as a candidate for reduced brain weight. To summarize, a BALB/c reduced complexity cross combined with multiple-omics approaches facilitated identification of candidate genes underlying thermal nociception and brain weight. These substrains provide a powerful, reciprocal platform for future validation of candidate variants.

Laboratory Codes in Nomenclature and Scientific Communication (Advancing Organism Nomenclature in Scientific Communication to Improve Research Reporting and Reproducibility)

Laboratory registration codes, also known as laboratory codes or lab codes, are a key element in standardized laboratory animal and genetic nomenclature. As such they are critical to accurate scientific communication and to research reproducibility and integrity. The original committee on Mouse Genetic Nomenclature published nomenclature conventions for mice genetics in 1940, and then conventions for inbred strains in 1952. Unique designations were needed, and have been in use since the 1950s, for the sources of animals and substrains, for the laboratories that identified new alleles or mutations, and then for developers of transgenes and induced mutations. Current laboratory codes are typically a 2- to 4-letter acronym for an institution or an investigator. Unique codes are assigned from the International Laboratory Code Registry, which was developed and is maintained by ILAR in the National Academies (National Academies of Sciences Engineering and Medicine and previously National Academy of Sciences). As a resource for the global research community, the registry has been online since 1997. Since 2003 mouse and rat genetic and strain nomenclature rules have been reviewed and updated annually as a joint effort of the International Committee on Standardized Genetic Nomenclature for Mice and the Rat Genome and Nomenclature Committee. The current nomenclature conventions (particularly conventions for non-inbred animals) are applicable beyond rodents, although not widely adopted. Ongoing recognition, since at least the 1930s, of the research relevance of genetic backgrounds and origins of animals, and of spontaneous and induced genetic variants speaks to the need for broader application of standardized nomenclature for animals in research, particularly given the increasing numbers and complexities of genetically modified swine, nonhuman primates, fish, and other species.

Cardiac proteomics reveals sex chromosome-dependent differences between males and females that arise prior to gonad formation

Sex disparities in cardiac homeostasis and heart disease are well documented, with differences attributed to actions of sex hormones. However, studies have indicated sex chromosomes act outside of the gonads to function without mediation by gonadal hormones. Here, we performed transcriptional and proteomics profiling to define differences between male and female mouse hearts. We demonstrate, contrary to current dogma, cardiac sex disparities are controlled not only by sex hormones but also through a sex-chromosome mechanism. Using Turner syndrome (XO) and Klinefelter (XXY) models, we find the sex-chromosome pathway is established by X-linked gene dosage. We demonstrate cardiac sex disparities occur at the earliest stages of heart formation, a period before gonad formation. Using these datasets, we identify and define a role for alpha-1B-glycoprotein (A1BG), showing loss of A1BG leads to cardiac defects in females, but not males. These studies provide resources for studying sex-biased cardiac disease states.

Patterns and mechanisms of sex ratio distortion in the Collaborative Cross mouse mapping population

In species with single-locus, chromosome-based mechanisms of sex determination, the laws of segregation predict an equal ratio of females to males at birth. Here, we show that departures from this Mendelian expectation are commonplace in the 8-way recombinant inbred Collaborative Cross (CC) mouse population. More than one-third of CC strains exhibit significant sex ratio distortion (SRD) at wean, with twice as many male-biased than female-biased strains. We show that these pervasive sex biases persist across multiple breeding environments, are stable over time, and are not mediated by random maternal effects. SRD exhibits a heritable component, but QTL mapping analyses fail to nominate any large effect loci. These findings, combined with the reported absence of sex ratio biases in the CC founder strains, suggest that SRD manifests from multilocus combinations of alleles only uncovered in recombined CC genomes. We explore several potential complex genetic mechanisms for SRD, including allelic interactions leading to sex-biased lethality, genetic sex reversal, chromosome drive mediated by sex-linked selfish elements, and incompatibilities between specific maternal and paternal genotypes. We show that no one mechanism offers a singular explanation for this population-wide SRD. Instead, our data present preliminary evidence for the action of distinct mechanisms of SRD at play in different strains. Taken together, our work exposes the pervasiveness of SRD in the CC population and nominates the CC as a powerful resource for investigating diverse genetic causes of biased sex chromosome transmission.

Regulation of protein abundance in genetically diverse mouse populations

Genetically diverse mouse populations are powerful tools for characterizing the regulation of the proteome and its relationship to whole-organism phenotypes. We used mass spectrometry to profile and quantify the abundance of 6,798 proteins in liver tissue from mice of both sexes across 58 Collaborative Cross (CC) inbred strains. We previously collected liver proteomics data from the related Diversity Outbred (DO) mice and their founder strains. We show concordance across the proteomics datasets despite being generated from separate experiments, allowing comparative analysis. We map protein abundance quantitative trait loci (pQTLs), identifying 1,087 local and 285 distal in the CC mice and 1,706 local and 414 distal in the DO mice. We find that regulatory effects on individual proteins are conserved across the mouse populations, in particular for local genetic variation and sex differences. In comparison, proteins that form complexes are often co-regulated, displaying varying genetic architectures, and overall show lower heritability and map fewer pQTLs. We have made this resource publicly available to enable quantitative analyses of the regulation of the proteome.

Importing genetically altered animals: ensuring quality

The reproducibility of research using laboratory animals requires reliable management of their quality, in particular of their genetics, health and environment, all of which contribute to their phenotypes. The point at which these biological materials are transferred between researchers is particularly sensitive, as it may result in a loss of integrity of the animals and/or their documentation. Here, we describe the various aspects of laboratory animal quality that should be confirmed when sharing rodent research models. We also discuss how repositories of biological materials support the scientific community to ensure the continuity of the quality of laboratory animals. Both the concept of quality and the role of repositories themselves extend to all exchanges of biological materials and all networks that support the sharing of these reagents.