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

Differential neuronal functions of LNX1 and LNX2 revealed by behavioural analysis in single and double knockout mice

BACKGROUND

Ligand of NUMB protein-X 1 (LNX1) and LNX2 proteins are closely related PDZ domain-containing E3 ubiquitin ligases that interact with and potentially modulate numerous synaptic and neurodevelopmentally important proteins. While both LNX1 and LNX2 are expressed in neurons, it is noteworthy that neuronal LNX1 isoforms lack the catalytic domain responsible for ubiquitination of substrates. Thus, the shared interaction partners of LNX1 and LNX2 might be differentially regulated by these proteins, with LNX1 acting as a stabilizing scaffold while LNX2 may promote their ubiquitination and degradation. Despite the identification of many LNX interacting proteins and substrates, our understanding of the distinct in vivo functions of LNX1 and LNX2 remains very incomplete.

RESULTS

We previously reported that mice lacking both LNX1 in the central nervous system and LNX2 globally exhibit decreased anxiety-related behaviour. Here we significantly extend this work by examining anxiety-related and risk-taking behaviours in Lnx1-/- and Lnx2-/- single knockout animals for the first time and by analysing previously unexplored aspects of behaviour in both single and double knockout animals. While the absence of both LNX1 and LNX2 contributes to the decreased anxiety-related behaviour of double knockout animals in the open field and elevated plus maze tests, the elimination of LNX2 plays a more prominent role in altered behaviour in the dark-light emergence test and wire beam bridge risk-taking paradigms. By contrast, Lnx knockout mice of all genotypes were indistinguishable from wildtype animals in the marble burying, stress-induced hyperthermia and novel object recognition tests. Analysis of the ultrasonic vocalizations of pups following maternal separation revealed significant differences in call properties and vocal repertoire for Lnx1-/- and Lnx1-/-;Lnx2-/- double knockout animals. Finally, decreased body weight previously noted in double knockout animals could be attributed largely to Lnx1 gene knockout.

CONCLUSIONS

These results identify specific roles of LNX1 and LNX2 proteins in modulating distinct aspects of anxiety and risk-taking behaviour and social communication in mice. They also reveal an unexpected role for neuronally expressed LNX1 isoforms in determining body weight. These novel insights into the differential neuronal functions of LNX1 and LNX2 proteins provide a foundation for mechanistic studies of these phenomena.

A Similar Mutation in the AAUU-Rich Elements of the Mouse TNF Gene Results in a Distinct Ileocolitic Phenotype: A New Strain of TNF-Overexpressing Mice

BACKGROUND

Tumor necrosis factor (TNF) is a pleiotropic cytokine that plays a critical role in the pathogenesis of immune-mediated diseases including inflammatory bowel disease (IBD). The stability of its mRNA transcript, determined in part by destabilizing sequences in its AAUU repeats (ARE) gene region, is an important regulator of its tissue and systemic levels. A deletion in the ARE region of the gene resulted in IBD and arthritis in mice and pigs, supporting a critical role for the cytokine in human IBD and several human arthritides. A mutation in the same area of the mouse genome by Genentech scientists (T.Y., M.K.) resulted in a similar but not identical phenotype.

METHODS

Here, we compare histopathological, cellular, and molecular features of the strains and propose reasons for their distinct phenotypes. First, while homozygous TNFΔARE mice develop severe arthritis and die after weaning, homozygous Genentech TNFΔARE (ΔG/ΔG) mice have normal lifespans, and males are often fertile.

RESULTS

We found that while the ileitic phenotype had peaked at 12 weeks of age in all mice, colitis progressed mostly after 20 weeks of age in heterozygous mice. Their variably penetrant arthritic phenotype progressed mostly after 20 weeks, also in heterozygous mice from both strains. There was expansion of central memory T and B cells in lymphoid organs of TNF-overproducing strains and their transcriptional profile shared well-known pathogenetic pathways with human IBD. Finally, we found differences in the mutated sequences within the ARE regions of the TNF gene and in their microbiota composition and genetic background. These differences likely explain their phenotypic differences.

CONCLUSIONS

In summary, we describe a different strain of TNF-overproducing mice with an overlapping, yet not identical phenotype, which may have differential applications than the original strain.

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.

Survey of gene, lncRNA and transposon transcription patterns in four mouse organs highlights shared and organ-specific sex-biased regulation

BACKGROUND

Sex-biased gene regulation is the basis of sexual dimorphism in phenotypes and has been studied across different cell types and different developmental stages. However, sex-biased expression of transposable elements (TEs), which represent nearly half of the mammalian genome and have the potential of influencing genome integrity and regulation, remains underexplored.

RESULTS

We report a survey of gene, lncRNA, and TE expression in four organs from mice with different combinations of gonadal and genetic sex. The data show remarkable variability among organs with respect to the impact of gonadal sex on transcription with the strongest effects observed in the liver. In contrast, the X-chromosome dosage alone had a modest influence on sex-biased transcription across organs, albeit interaction between X-dosage and gonadal sex cannot be ruled out. The presence of the Y-chromosome influences TE, but not gene or lncRNA, expression in the liver. Notably, 90% of sex-biased TEs (sDETEs) reside in clusters. Moreover, 54% of these clusters overlap or reside less than 100 kb from sex-biased genes or lncRNAs, share the same sex bias, and also have higher expression levels than sDETE clusters that do not co-localize with other types of sex-biased transcripts. We test the heterochromatic sink hypothesis that predicts higher expression of TEs in XX individuals finding no evidence to support it.

CONCLUSIONS

Our data show that sex-biased expression of TEs varies among organs with the highest numbers of sDETEs found in the liver following trends observed for genes and lncRNAs. It is enhanced by proximity to other types of sex-biased transcripts.

An ALPK3 truncation variant causing autosomal dominant hypertrophic cardiomyopathy is partially rescued by mavacamten

The ALPK3 gene encodes alpha-protein kinase 3, a cardiac pseudo-kinase of unknown function. Heterozygous truncating variants (ALPK3tv) can cause dominant adult-onset hypertrophic cardiomyopathy (HCM). Here we confirm an excess of ALPK3tv in sarcomere-gene negative HCM patients. Moreover, we generated a novel knock-in mouse model carrying an ALPK3tv (K201X). Homozygous animals displayed hypertrophy and systolic dysfunction. Heterozygous animals demonstrated no obvious baseline; however, they had an aggravated hypertrophic response upon chronic adrenergic challenge. Isolated, unloaded cardiomyocytes from heterozygous and homozygous mice showed reduced basal sarcomere length with prolonged relaxation, whilst calcium transients showed increased diastolic calcium levels. Protein kinase A-mediated phosphorylation, including that of cardiac troponin I, was significantly decreased. In agreement with the cellular HCM phenotype, reduced ratios of myosin heads in the super-relaxed state were measured. Contractile and calcium handling defects were partly corrected by treatment with mavacamten, a novel myosin inhibitor. For the first time with a non-sarcomere HCM variant, we have demonstrated hallmark changes in cardiac contractility and calcium handling. Mavacamten is able to partially rescue the cellular phenotype, hence could be beneficial to HCM patients with ALPK3tv. Moreover, our data points at a potential role of ALPK3 as a modulator of protein kinase A signalling.

A conditional smoothened (smo) allele on an inbred C57BL/6J genetic background has a hypomorphic smo mutant phenotype

We have introduced the floxed allele of Smoothened (Smo) carried by the mouse line Smotm2Amc into the C57BL/6J strain by serial backcross. Recapitulation of the Smo null phenotype was confirmed by deleting the allele using E2a-cre and intercrossing heterozygous Smo ± mice. No homozygous mutant embryos were identified at E9.5, suggesting the null phenotype is at least as severe as that observed on a mixed genetic background. While healthy and fertile homozygous floxed mice were regularly obtained after intercrosses, their numbers at weaning were reduced relative to Mendelian expectation, suggesting the unrecombined allele is itself hypomorphic. This hypothesis is supported by characterization of transcription of the floxed allele, which revealed that its expression was variably reduced relative to wild-type Smo.

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.

Metabolism of inorganic arsenic in mice carrying the human AS3MT gene and fed folate deficient or folate supplemented diet

Arsenic (+3 oxidation state) methyltransferase (AS3MT) catalyzes the S-adenosylmethionine (SAM)-dependent methylation of inorganic arsenic (iAs), yielding monomethyl‑arsenic (MAs) and dimethyl‑arsenic (DMAs) metabolites. The formation of DMAs in this pathway is considered a key mechanism for iAs detoxification. Availability of SAM for iAs methylation depends in part on dietary intake of folate. Results of population studies suggest that supplementation with folate stimulates iAs methylation, increasing DMAs and decreasing iAs and MAs proportions in urine and/or blood. The goal of the present study was to determine if folate intake affects methylation and clearance of iAs in a recently established mouse strain that expresses human AS3MT and exhibits a human-like pattern of iAs metabolism. The humanized male and female mice were fed folate-deficient (FD) or folate-supplemented (FS) diet for 6 weeks, followed by exposure to 0 ppb or 400 ppb iAs in drinking water for 5 weeks, while on the same types of diet. The concentrations and proportions of iAs, MAs and DMAs were determined in urine, liver, kidneys, and spleen. The diet-, sex- and dose-related differences were assessed by t-test or a non-parametric test; Bonferroni test was used to correct for multiple comparisons. In general, proportions of DMAs were greater and proportions of iAs were smaller in urine and tissues of FS mice as compared to FD mice. However, folate supplementation also increased MAs proportions. Notably, the folate intake had no effect on the concentrations of total arsenic either in the urine or the tissues. These results suggest that, similar to humans, folate supplementation stimulates iAs methylation in the humanized mice. However, the stimulation of iAs methylation is not associated with clearance of arsenic from tissues, possibly due to an inefficient conversion of MAs to DMAs.

Characterisation of sleep in a mouse model of CLN3 disease revealed sex-specific sleep disturbances

The neuronal ceroid lipofuscinoses (NCLs) are a group of recessively inherited neurodegenerative diseases characterizsed by lysosomal storage of fluorescent materials. CLN3 disease, or juvenile Batten disease, is the most common NCL that is caused by mutations in the Ceroid Lipofuscinosis, Neuronal 3 (CLN3) gene. Sleep disturbances are among the most common symptoms associated with CLN3 disease that deteriorate the patients' life quality, yet this is understudied and has not been delineated in animal models of the disease. The current study utilised PiezoSleep, a non-invasive, automated piezoelectric motion sensing system, to classify sleep and wakefulness in a Cln3Δex1-6/Δex1-6 (Cln3KO) mouse model and age- and sex-matched wild-type (WT) controls. The sleep-wake classification by PiezoSleep was found to be about 90% accurate when validated against simultaneous polysomnographic recordings including electroencephalography (EEG) and electromyography (EMG) in a small cohort of WT and Cln3KO mice. Our large cohort PiezoSleep study revealed sleep abnormalities during the light period in male Cln3KO mice compared with WT male mice, and more subtle differences in Cln3KO female mice in the dark period compared with WT female mice. Our characterisation of sleep in the Cln3KO mouse model aligns with sleep abnormalities seen in CLN3 disease patients and serves as a basis to continue examining sleep disturbances commonly reported for CLN3 disease and other NCLs. As the first animal model study capturing sleep disturbances in CLN3 disease, our work will facilitate future studies into the potential mechanism behind sleep disturbances associated with the disease and the potential treatment strategies.

Pathobiont-triggered induction of epithelial IDO1 drives regional susceptibility to Inflammatory Bowel Disease

The structure and function of the mammalian gut vary by region, yet why inflammatory diseases manifest in specific regions and not others remains unclear. We use a TNF-overexpressing Crohn's disease (CD) model (TnfΔARE/+), which typically presents in the terminal ileum (TI), to investigate how environmental factors interact with the host's immune susceptibility to drive region-specific disease. We identified Chlamydia muridarum, an intracellular bacterium and murine counterpart to the human sexually transmitted C. trachomatis, as necessary and sufficient to trigger disease manifestation in the ascending colon (AC), another common site of human CD. Disease manifestation in the AC depends on indoleamine 2,3-dioxygenase (IDO1) activity induced by hypersensitive surface secretory cells in genetically susceptible hosts. Single-cell and microbial analyses of human specimens also implicates this pathobiont-epithelial IDO1 pathway in patients with a history of CD in the AC. Our findings demonstrate that genetic and microbial factors can independently drive region-specific disease and provide a unique model to study CD specific to the AC.

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.

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.