Genetic loci regulate Sarbecovirus pathogenesis: A comparison across mice and humans

Coronavirus (CoV) cause considerable morbidity and mortality in humans and other mammals, as evidenced by the emergence of Severe Acute Respiratory CoV (SARS-CoV) in 2003, Middle East Respiratory CoV (MERS-CoV) in 2012, and SARS-CoV-2 in 2019. Although poorly characterized, natural genetic variation in human and other mammals modulate virus pathogenesis, as reflected by the spectrum of clinical outcomes ranging from asymptomatic infections to lethal disease. Using multiple human epidemic and zoonotic Sarbecoviruses, coupled with murine Collaborative Cross genetic reference populations, we identify several dozen quantitative trait loci that regulate SARS-like group-2B CoV pathogenesis and replication. Under a Chr4 QTL, we deleted a candidate interferon stimulated gene, Trim14 which resulted in enhanced SARS-CoV titers and clinical disease, suggesting an antiviral role during infection. Importantly, about 60 % of the murine QTL encode susceptibility genes identified as priority candidates from human genome-wide association studies (GWAS) studies after SARS-CoV-2 infection, suggesting that similar selective forces have targeted analogous genes and pathways to regulate Sarbecovirus disease across diverse mammalian hosts. These studies provide an experimental platform in rodents to investigate the molecular-genetic mechanisms by which potential cross mammalian susceptibility loci and genes regulate type-specific and cross-SARS-like group 2B CoV replication, immunity, and pathogenesis in rodent models. Our study also provides a paradigm for identifying susceptibility loci for other highly heterogeneous and virulent viruses that sporadically emerge from zoonotic reservoirs to plague human and animal populations.

Adjuvant-dependent impact of inactivated SARS-CoV-2 vaccines during heterologous infection by a SARS-related coronavirus

Whole virus-based inactivated SARS-CoV-2 vaccines adjuvanted with aluminum hydroxide have been critical to the COVID-19 pandemic response. Although these vaccines are protective against homologous coronavirus infection, the emergence of novel variants and the presence of large zoonotic reservoirs harboring novel heterologous coronaviruses provide significant opportunities for vaccine breakthrough, which raises the risk of adverse outcomes like vaccine-associated enhanced respiratory disease. Here, we use a female mouse model of coronavirus disease to evaluate inactivated vaccine performance against either homologous challenge with SARS-CoV-2 or heterologous challenge with a bat-derived coronavirus that represents a potential emerging disease threat. We show that inactivated SARS-CoV-2 vaccines adjuvanted with aluminum hydroxide can cause enhanced respiratory disease during heterologous infection, while use of an alternative adjuvant does not drive disease and promotes heterologous viral clearance. In this work, we highlight the impact of adjuvant selection on inactivated vaccine safety and efficacy against heterologous coronavirus infection.

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 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.

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.

U-CAN-seq: A Universal Competition Assay by Nanopore Sequencing

RNA viruses quickly evolve subtle genotypic changes that can have major impacts on viral fitness and host range, with potential consequences for human health. It is therefore important to understand the evolutionary fitness of novel viral variants relative to well-studied genotypes of epidemic viruses. Competition assays are an effective and rigorous system with which to assess the relative fitness of viral genotypes. However, it is challenging to quickly and cheaply distinguish and quantify fitness differences between very similar viral genotypes. Here, we describe a protocol for using reverse transcription PCR in combination with commercial nanopore sequencing services to perform competition assays on untagged RNA viruses. Our assay, called the Universal Competition Assay by Nanopore Sequencing (U-CAN-seq), is relatively cheap and highly sensitive. We used a well-studied N24A mutation in the chikungunya virus (CHIKV) nsp3 gene to confirm that we could detect a competitive disadvantage using U-CAN-seq. We also used this approach to show that mutations to the CHIKV 5' conserved sequence element that disrupt sequence but not structure did not affect the fitness of CHIKV. However, similar mutations to an adjacent CHIKV stem loop (SL3) did cause a fitness disadvantage compared to wild-type CHIKV, suggesting that structure-independent, primary sequence determinants in this loop play an important role in CHIKV biology. Our novel findings illustrate the utility of the U-CAN-seq competition assay.

Unique immune profiles in collaborative cross mice linked to survival and viral clearance upon infection

The response to infection is generally heterogeneous and diverse, with some individuals remaining asymptomatic while others present with severe disease or a diverse range of symptoms. Here, we address the role of host genetics on immune phenotypes and clinical outcomes following viral infection by studying genetically diverse mice from the Collaborative Cross (CC), allowing for use of a small animal model with controlled genetic diversity while maintaining genetic replicates. We demonstrate variation by deeply profiling a broad range of innate and adaptive immune cell phenotypes at steady-state in 63 genetically distinct CC mouse strains and link baseline immune signatures with virologic and clinical disease outcomes following infection of mice with herpes simplex virus 2 (HSV-2) or severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This work serves as a resource for CC strain selection based on steady-state immune phenotypes or disease presentation upon viral infection, and further, points to possible pre-infection immune correlates of survival and early viral clearance upon infection.

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.

Host genetic variation guides hepacivirus clearance, chronicity, and liver fibrosis in mice

BACKGROUND AIMS

Human genetic variation is thought to guide the outcome of HCV infection, but model systems within which to dissect these host genetic mechanisms are limited. Norway rat hepacivirus, closely related to HCV, causes chronic liver infection in rats but causes acute self-limiting hepatitis in typical strains of laboratory mice, which resolves in 2 weeks. The Collaborative Cross (CC) is a robust mouse genetics resource comprised of a panel of recombinant inbred strains, which model the complexity of the human genome and provide a system within which to understand diseases driven by complex allelic variation.

APPROACH RESULTS

We infected a panel of CC strains with Norway rat hepacivirus and identified several that failed to clear the virus after 4 weeks. Strains displayed an array of virologic phenotypes ranging from delayed clearance (CC046) to chronicity (CC071, CC080) with viremia for at least 10 months. Body weight loss, hepatocyte infection frequency, viral evolution, T-cell recruitment to the liver, liver inflammation, and the capacity to develop liver fibrosis varied among infected CC strains.

CONCLUSIONS

These models recapitulate many aspects of HCV infection in humans and demonstrate that host genetic variation affects a multitude of viruses and host phenotypes. These models can be used to better understand the molecular mechanisms that drive hepacivirus clearance and chronicity, the virus and host interactions that promote chronic disease manifestations like liver fibrosis, therapeutic and vaccine performance, and how these factors are affected by host genetic variation.

Mice with FVB-derived sequence on chromosome 17 succumb to disseminated virus infection due to aberrant NK cell and T cell responses

Zoonotic arenavirus infections can result in viral hemorrhagic disease, characterized by platelet loss, petechia, and multi-organ injury. The mechanisms governing these outcomes are likely impacted by virus strain and infection dose, as well as an individual's genetic background and immune constitution. To better understand the processes leading to severe pathogenesis, we compared two strains of inbred mice, C57BL/6J (B6) and FVB/NJ (FVB), that have diametrically opposed outcomes during disseminated lymphocytic choriomeningitis virus (LCMV) infection. Infection caused minimal pathogenesis in B6 mice, whereas FVB mice developed acute hepatitis and perished due, in part, to aberrant NK cell and T cell responses. Susceptible mice showed an outgrowth of cytolytic CD4+ T cells and loss of Treg cells. B6 congenic mice with the FVB allele at a 25Mb locus on chromosome 17 recapitulated FVB pathogenesis upon infection. A locus containing a limited number of variants in immune-related genes greatly impacts survival during infection.

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.

Host genetic background is a barrier to broadly effective vaccine-mediated protection against tuberculosis

Heterogeneity in human immune responses is difficult to model in standard laboratory mice. To understand how host variation affects Bacillus Calmette Guerin-induced (BCG-induced) immunity against Mycobacterium tuberculosis, we studied 24 unique collaborative cross (CC) mouse strains, which differ primarily in the genes and alleles they inherit from founder strains. The CC strains were vaccinated with or without BCG and challenged with aerosolized M. tuberculosis. Since BCG protects only half of the CC strains tested, we concluded that host genetics has a major influence on BCG-induced immunity against M. tuberculosis infection, making it an important barrier to vaccine-mediated protection. Importantly, BCG efficacy is dissociable from inherent susceptibility to tuberculosis (TB). T cell immunity was extensively characterized to identify components associated with protection that were stimulated by BCG and recalled after M. tuberculosis infection. Although considerable diversity is observed, BCG has little impact on the composition of T cells in the lung after infection. Instead, variability is largely shaped by host genetics. BCG-elicited protection against TB correlated with changes in immune function. Thus, CC mice can be used to define correlates of protection and to identify vaccine strategies that protect a larger fraction of genetically diverse individuals instead of optimizing protection for a single genotype.

Immune resilience despite inflammatory stress promotes longevity and favorable health outcomes including resistance to infection

Some people remain healthier throughout life than others but the underlying reasons are poorly understood. Here we hypothesize this advantage is attributable in part to optimal immune resilience (IR), defined as the capacity to preserve and/or rapidly restore immune functions that promote disease resistance (immunocompetence) and control inflammation in infectious diseases as well as other causes of inflammatory stress. We gauge IR levels with two distinct peripheral blood metrics that quantify the balance between (i) CD8+ and CD4+ T-cell levels and (ii) gene expression signatures tracking longevity-associated immunocompetence and mortality-associated inflammation. Profiles of IR metrics in ~48,500 individuals collectively indicate that some persons resist degradation of IR both during aging and when challenged with varied inflammatory stressors. With this resistance, preservation of optimal IR tracked (i) a lower risk of HIV acquisition, AIDS development, symptomatic influenza infection, and recurrent skin cancer; (ii) survival during COVID-19 and sepsis; and (iii) longevity. IR degradation is potentially reversible by decreasing inflammatory stress. Overall, we show that optimal IR is a trait observed across the age spectrum, more common in females, and aligned with a specific immunocompetence-inflammation balance linked to favorable immunity-dependent health outcomes. IR metrics and mechanisms have utility both as biomarkers for measuring immune health and for improving health outcomes.

Neuroinvasive Flavivirus Pathogenesis Is Restricted by Host Genetic Factors in Collaborative Cross Mice, Independently of Oas1b

Powassan virus (POWV) is an emerging tick-borne flavivirus that causes neuroinvasive diseases, including encephalitis, meningitis, and paralysis. Similar to other neuroinvasive flaviviruses, such as West Nile virus (WNV) and Japanese encephalitis virus (JEV), POWV disease presentation is heterogeneous, and the factors influencing disease outcome are not fully understood. We used Collaborative Cross (CC) mice to assess the impact of host genetic factors on POWV pathogenesis. We infected a panel of Oas1b-null CC lines with POWV and observed a range of susceptibility, indicating that host factors other than the well-characterized flavivirus restriction factor Oas1b modulate POWV pathogenesis in CC mice. Among the Oas1b-null CC lines, we identified multiple highly susceptible lines (0% survival), including CC071 and CC015, and two resistant lines, CC045 and CC057 (>75% survival). The susceptibility phenotypes generally were concordant among neuroinvasive flaviviruses, although we did identify one line, CC006, that was specifically resistant to JEV, suggesting that both pan-flavivirus and virus-specific mechanisms contribute to susceptibility phenotypes in CC mice. We found that POWV replication was restricted in bone marrow-derived macrophages from CC045 and CC057 mice, suggesting that resistance could result from cell-intrinsic restriction of viral replication. Although serum viral loads at 2 days postinfection were equivalent between resistant and susceptible CC lines, clearance of POWV from the serum was significantly enhanced in CC045 mice. Furthermore, CC045 mice had significantly lower viral loads in the brain at 7 days postinfection than did CC071 mice, suggesting that reduced central nervous system (CNS) infection contributes to the resistant phenotype of CC045 mice. IMPORTANCE Neuroinvasive flaviviruses, such as WNV, JEV, and POWV, are transmitted to humans by mosquitoes or ticks and can cause neurologic diseases, such as encephalitis, meningitis, and paralysis, and they can result in death or long-term sequelae. Although potentially severe, neuroinvasive disease is a rare outcome of flavivirus infection. The factors that determine whether someone develops severe disease after a flavivirus infection are not fully understood, but host genetic differences in polymorphic antiviral response genes likely contribute to the outcome of infection. We evaluated a panel of genetically diverse mice and identified lines with distinct outcomes following infection with POWV. We found that resistance to POWV pathogenesis corresponded to reduced viral replication in macrophages, more rapid clearance of virus in peripheral tissues, and reduced viral infection in the brain. These susceptible and resistant mouse lines will provide a system for investigating the pathogenic mechanisms of POWV and identifying polymorphic host genes that contribute to resistance.

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.

Memory T cells possess an innate-like function in local protection from mucosal infection

Mucosal infections pose a significant global health burden. Antigen-specific tissue resident T cells are critical to maintaining barrier immunity. Previous studies in the context of systemic infection suggest that memory CD8 T cells may also provide innate-like protection against antigenically unrelated pathogens independent of TCR engagement. Whether "bystander T cell activation" is also an important defense mechanism in the mucosa is poorly understood. Here, we investigated if innate-like memory CD8 T cells could protect against a model mucosal virus infection, herpes simplex virus 2 (HSV-2). We found that immunization with an irrelevant antigen delayed disease progression from lethal HSV-2 challenge, suggesting that memory CD8 T cells may mediate protection despite the lack of antigen-specificity. Upon HSV-2 infection, we observed an early infiltration, rather than substantial local proliferation, of antigen-non-specific CD8 T cells, which became bystander-activated only within the infected mucosal tissue. Critically, we show that bystander-activated CD8 T cells are sufficient to reduce early viral burden after HSV-2 infection. Finally, local cytokine cues within the tissue microenvironment after infection were sufficient for bystander activation of mucosal tissue memory CD8 T cells from mice and humans. Altogether, our findings suggest that local bystander-activation of CD8 memory T cells contribute a fast and effective innate-like response to infection in mucosal tissue.

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.

Human ACE2 expression, a major tropism determinant for SARS-CoV-2, is regulated by upstream and intragenic elements

Angiotensin-converting enzyme 2 (ACE2), part of the renin-angiotensin system (RAS), serves as an entry point for SARS-CoV-2, leading to viral proliferation in permissive cell types. Using mouse lines in which the Ace2 locus has been humanized by syntenic replacement, we show that regulation of basal and interferon induced ACE2 expression, relative expression levels of different ACE2 transcripts, and sexual dimorphism in ACE2 expression are unique to each species, differ between tissues, and are determined by both intragenic and upstream promoter elements. Our results indicate that the higher levels of expression of ACE2 observed in the lungs of mice relative to humans may reflect the fact that the mouse promoter drives expression of ACE2 in populous airway club cells while the human promoter drives expression in alveolar type 2 (AT2) cells. In contrast to transgenic mice in which human ACE2 is expressed in ciliated cells under the control of the human FOXJ1 promoter, mice expressing ACE2 in club cells under the control of the endogenous Ace2 promoter show a robust immune response after infection with SARS-CoV-2, leading to rapid clearance of the virus. This supports a model in which differential expression of ACE2 determines which cell types in the lung are infected, and this in turn modulates the host response and outcome of COVID-19.

Forward genetic screen of homeostatic antibody levels in the Collaborative Cross identifies MBD1 as a novel regulator of B cell homeostasis

Variation in immune homeostasis, the state in which the immune system is maintained in the absence of stimulation, is highly variable across populations. This variation is attributed to both genetic and environmental factors. However, the identity and function of specific regulators have been difficult to identify in humans. We evaluated homeostatic antibody levels in the serum of the Collaborative Cross (CC) mouse genetic reference population. We found heritable variation in all antibody isotypes and subtypes measured. We identified 4 quantitative trait loci (QTL) associated with 3 IgG subtypes: IgG1, IgG2b, and IgG2c. While 3 of these QTL map to genome regions of known immunological significance (major histocompatibility and immunoglobulin heavy chain locus), Qih1 (associated with variation in IgG1) mapped to a novel locus on Chromosome 18. We further associated this locus with B cell proportions in the spleen and identify Methyl-CpG binding domain protein 1 under this locus as a novel regulator of homeostatic IgG1 levels in the serum and marginal zone B cells (MZB) in the spleen, consistent with a role in MZB differentiation to antibody secreting cells.

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.

Genetic diversity of collaborative cross mice enables identification of novel rift valley fever virus encephalitis model

Rift Valley fever (RVF) is an arboviral disease of humans and livestock responsible for severe economic and human health impacts. In humans, RVF spans a variety of clinical manifestations, ranging from an acute flu-like illness to severe forms of disease, including late-onset encephalitis. The large variations in human RVF disease are inadequately represented by current murine models, which overwhelmingly die of early-onset hepatitis. Existing mouse models of RVF encephalitis are either immunosuppressed, display an inconsistent phenotype, or develop encephalitis only when challenged via intranasal or aerosol exposure. In this study, the genetically defined recombinant inbred mouse resource known as the Collaborative Cross (CC) was used to identify mice with additional RVF disease phenotypes when challenged via a peripheral foot-pad route to mimic mosquito-bite exposure. Wild-type Rift Valley fever virus (RVFV) challenge of 20 CC strains revealed three distinct disease phenotypes: early-onset hepatitis, mixed phenotype, and late-onset encephalitis. Strain CC057/Unc, with the most divergent phenotype, which died of late-onset encephalitis at a median of 11 days post-infection, is the first mouse strain to develop consistent encephalitis following peripheral challenge. CC057/Unc mice were directly compared to C57BL/6 mice, which uniformly succumb to hepatitis within 2–4 days of infection. Encephalitic disease in CC057/Unc mice was characterized by high viral RNA loads in brain tissue, accompanied by clearance of viral RNA from the periphery, low ALT levels, lymphopenia, and neutrophilia. In contrast, C57BL/6 mice succumbed from hepatitis at 3 days post-infection with high viral RNA loads in the liver, viremia, high ALT levels, lymphopenia, and thrombocytopenia. The identification of a strain of CC mice as an RVFV encephalitis model will allow for future investigation into the pathogenesis and treatment of RVF encephalitic disease and indicates that genetic background makes a major contribution to RVF disease variation.

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.

Genetic analysis of differential responses to adjuvanted influenza vaccination

Each year the seasonal influenza vaccine is administered to protect individuals from infection, but vaccination often shows variable efficacy against influenza infections. While the inclusion of adjuvants can boost vaccine efficacy, individuals also vary in their response to adjuvants, and understanding the basis of this variation would provide insights into adjuvant mechanisms and advance vaccine development. Host genetics contribute to variation in vaccine response, but the specific genetic factors that impact immune responses to adjuvanted vaccination are not well understood. To better understand how the genetic background of individuals can modify vaccine and adjuvant-induced immune responses, we screened a set of Collaborative Cross (CC) mouse strains with immunization via inactivated influenza with or without alum. Across these genetically diverse CC strains, we observed wide variation in levels of influenza specific antibodies and isotypes (e.g., IgG1, IgG2). We further focused on two strains, CC002 and CC072, that had contrasting antibody responses to vaccine + alum. To understand the genetic factors driving these responses, we generated an F2 population to map quantitative trait loci (qtl) driving the differences between CC002 and CC072 in vaccine responses. Among the F2 population, variation across antibody levels were observed in mice vaccinated with the vaccine with and without alum that reflected the variation in the parent strains. QTL were mapped for variation in antibody levels for mice vaccinated without alum, further supporting a role for host genetics in regulating vaccine responses. Future analysis will focus on analyzing specific host factors modulating antibody responses to adjuvanted vaccination.

Supported by grants from NIH (U19 AI100625, 5T32 GM007092, 1T32 GM135128, 5R25 GM055336-19) and Burroughs Wellcome Fund GDEP to MCC. 

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.

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.

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.

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.

Common Mechanism of SARS-CoV and SARS-CoV-2 Pathogenesis across Species

Sarbecovirus (CoV) infections, including Severe Acute Respiratory CoV (SARS-CoV) and SARS-CoV-2, are considerable human threats. Human GWAS studies have recently identified loci associated with variation in SARS-CoV-2 susceptibility. However, genetically tractable models that reproduce human CoV disease outcomes are needed to mechanistically evaluate genetic determinants of CoV susceptibility. We used the Collaborative Cross (CC) and human GWAS datasets to elucidate host susceptibility loci that regulate CoV infections and to identify host quantitative trait loci that modulate severe CoV and pan-CoV disease outcomes including a major disease regulating loci including CCR9. CCR9 ablation resulted in enhanced titer, weight loss, respiratory dysfunction, mortality, and inflammation, providing mechanistic support in mitigating protection from severe SARS-CoV-2 pathogenesis across species. This study represents a comprehensive analysis of susceptibility loci for an entire genus of human pathogens conducted, identifies a large collection of susceptibility loci and candidate genes that regulate multiple aspects type-specific and cross-CoV pathogenesis, and also validates the paradigm of using the CC platform to identify common cross-species susceptibility loci and genes for newly emerging and pre-epidemic viruses.

Genetic regulation of immune homeostatic lung leukocyte populations influences respiratory virus induced disease in collaborative cross mice

Immune homeostasis is the state that the immune system maintains in the absence of insult. Perturbation of immune homeostasis can impact autoimmunity/allergy and adversely affect infectious responses. To date, much of the analysis of immune homeostasis has focused on systemic immunity, but it is also likely to be important in an organ specific manner. Since the lungs are a major site of environmental exposure and infection, we used the Collaborative Cross (CC) mouse genetic reference population to study the genetic regulation of the breadth of baseline immune cell populations in the lung and identify loci regulating these cells at homeostasis. We found that all 54 immune cell phenotypes measured showed strong genetic variation in cell type abundances. We identified 28 quantitative trait loci associated with variation in 24 immune cell populations or the relationship between cell populations. Further, we identified significant associations between 20 of these loci and responses to either influenza A virus (IAV) or Severe acute respiratory syndrome associated coronavirus (SARS-CoV) disease in the same strains of mice. Notably, a locus mapped for variation in Ly6C+ monocyte abundance was associated with SARS-CoV weight loss, as well as titer at days 2 and 4 post-infection. This locus is also associated with influenza-induced disease, as measured by weight loss post-infection. Our analysis highlights the strong genetic control of immune homeostasis, and the key role that immune homeostasis plays in contributing to downstream infectious responses. In particular, our analysis suggests that the abundance of a variety of lung leukocyte populations prior to infection could serve as predictors of immune responses to respiratory viruses.

Correlation of Regulatory T Cell Numbers with Disease Tolerance upon Virus Infection

The goal of a successful immune response is to clear the pathogen while sparing host tissues from damage associated with pathogen replication and active immunity. Regulatory T cells (Treg) have been implicated in maintaining this balance as they contribute both to the organization of immune responses as well as restriction of inflammation and immune activation to limit immunopathology. To determine if Treg abundance prior to pathogen encounter can be used to predict the success of an antiviral immune response, we used genetically diverse mice from the collaborative cross infected with West Nile virus (WNV). We identified collaborative cross lines with extreme Treg abundance at steady state, either high or low, and used mice with these extreme phenotypes to demonstrate that baseline Treg quantity predicted the magnitude of the CD8 T cell response to WNV infection, although higher numbers of baseline Tregs were associated with reduced CD8 T cell functionality in terms of TNF and granzyme B expression. Finally, we found that abundance of CD44⁺ Tregs in the spleen at steady state was correlated with an increased early viral load within the spleen without an association with clinical disease. Thus, we propose that Tregs participate in disease tolerance in the context of WNV infection by tuning an appropriately focused and balanced immune response to control the virus while at the same time minimizing immunopathology and clinical disease. We hypothesize that Tregs limit the antiviral CD8 T cell function to curb immunopathology at the expense of early viral control as an overall host survival strategy.

Fecal IgA, Antigen Absorption, and Gut Microbiome Composition Are Associated With Food Antigen Sensitization in Genetically Susceptible Mice

Food allergy is a potentially fatal disease affecting 8% of children and has become increasingly common in the past two decades. Despite the prevalence and severe nature of the disease, the mechanisms underlying sensitization remain to be further elucidated. The Collaborative Cross is a genetically diverse panel of inbred mice that were specifically developed to study the influence of genetics on complex diseases. Using this panel of mouse strains, we previously demonstrated CC027/GeniUnc mice, but not C3H/HeJ mice, develop peanut allergy after oral exposure to peanut in the absence of a Th2-skewing adjuvant. Here, we investigated factors associated with sensitization in CC027/GeniUnc mice following oral exposure to peanut, walnut, milk, or egg. CC027/GeniUnc mice mounted antigen-specific IgE responses to peanut, walnut and egg, but not milk, while C3H/HeJ mice were not sensitized to any antigen. Naïve CC027/GeniUnc mice had markedly lower total fecal IgA compared to C3H/HeJ, which was accompanied by stark differences in gut microbiome composition. Sensitized CC027/GeniUnc mice had significantly fewer CD3+ T cells but higher numbers of CXCR5+ B cells and T follicular helper cells in the mesenteric lymph nodes compared to C3H/HeJ mice, which is consistent with their relative immunoglobulin production. After oral challenge to the corresponding food, peanut- and walnut-sensitized CC027/GeniUnc mice experienced anaphylaxis, whereas mice exposed to milk and egg did not. Ara h 2 was detected in serum collected post-challenge from peanut-sensitized mice, indicating increased absorption of this allergen, while Bos d 5 and Gal d 2 were not detected in mice exposed to milk and egg, respectively. Machine learning on the change in gut microbiome composition as a result of food protein exposure identified a unique signature in CC027/GeniUnc mice that experienced anaphylaxis, including the depletion of Akkermansia. Overall, these results demonstrate several factors associated with enteral sensitization in CC027/GeniUnc mice, including diminished total fecal IgA, increased allergen absorption and altered gut microbiome composition. Furthermore, peanuts and tree nuts may have inherent properties distinct from milk and eggs that contribute to allergy.

Baseline T cell immune phenotypes predict virologic and disease control upon SARS-CoV infection in Collaborative Cross mice

The COVID-19 pandemic has revealed that infection with SARS-CoV-2 can result in a wide range of clinical outcomes in humans. An incomplete understanding of immune correlates of protection represents a major barrier to the design of vaccines and therapeutic approaches to prevent infection or limit disease. This deficit is largely due to the lack of prospectively collected, pre-infection samples from individuals that go on to become infected with SARS-CoV-2. Here, we utilized data from genetically diverse Collaborative Cross (CC) mice infected with SARS-CoV to determine whether baseline T cell signatures are associated with a lack of viral control and severe disease upon infection. SARS-CoV infection of CC mice results in a variety of viral load trajectories and disease outcomes. Overall, a dysregulated, pro-inflammatory signature of circulating T cells at baseline was associated with severe disease upon infection. Our study serves as proof of concept that circulating T cell signatures at baseline can predict clinical and virologic outcomes upon SARS-CoV infection. Identification of basal immune predictors in humans could allow for identification of individuals at highest risk of severe clinical and virologic outcomes upon infection, who may thus most benefit from available clinical interventions to restrict infection and disease.

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

The laboratory mouse is the most widely used animal model for biomedical research, due in part to its well-annotated genome, wealth of genetic resources, and the ability to precisely manipulate its genome. Despite the importance of genetics for mouse research, genetic quality control (QC) is not standardized, in part due to the lack of cost-effective, informative, and robust platforms. Genotyping arrays are standard tools for mouse research and remain an attractive alternative even in the era of high-throughput whole-genome sequencing. Here, we describe the content and performance of a new iteration of the Mouse Universal Genotyping Array (MUGA), MiniMUGA, an array-based genetic QC platform with over 11,000 probes. In addition to robust discrimination between most classical and wild-derived laboratory strains, MiniMUGA was designed to contain features not available in other platforms: (1) chromosomal sex determination, (2) discrimination between substrains from multiple commercial vendors, (3) diagnostic SNPs for popular laboratory strains, (4) detection of constructs used in genetically engineered mice, and (5) an easy-to-interpret report summarizing these results. In-depth annotation of all probes should facilitate custom analyses by individual researchers. To determine the performance of MiniMUGA, we genotyped 6899 samples from a wide variety of genetic backgrounds. The performance of MiniMUGA compares favorably with three previous iterations of the MUGA family of arrays, both in discrimination capabilities and robustness. We have generated publicly available consensus genotypes for 241 inbred strains including classical, wild-derived, and recombinant inbred lines. Here, we also report the detection of a substantial number of XO and XXY individuals across a variety of sample types, new markers that expand the utility of reduced complexity crosses to genetic backgrounds other than C57BL/6, and the robust detection of 17 genetic constructs. We provide preliminary evidence that the array can be used to identify both partial sex chromosome duplication and mosaicism, and that diagnostic SNPs can be used to determine how long inbred mice have been bred independently from the relevant main stock. We conclude that MiniMUGA is a valuable platform for genetic QC, and an important new tool to increase the rigor and reproducibility of mouse research.

Baseline T cell immune phenotypes predict virologic and disease control upon SARS-CoV infection

The COVID-19 pandemic has revealed that infection with SARS-CoV-2 can result in a wide range of clinical outcomes in humans, from asymptomatic or mild disease to severe disease that can require mechanical ventilation. An incomplete understanding of immune correlates of protection represents a major barrier to the design of vaccines and therapeutic approaches to prevent infection or limit disease. This deficit is largely due to the lack of prospectively collected, pre-infection samples from indiviuals that go on to become infected with SARS-CoV-2. Here, we utilized data from a screen of genetically diverse mice from the Collaborative Cross (CC) infected with SARS-CoV to determine whether circulating baseline T cell signatures are associated with a lack of viral control and severe disease upon infection. SARS-CoV infection of CC mice results in a variety of viral load trajectories and disease outcomes. Further, early control of virus in the lung correlates with an increased abundance of activated CD4 and CD8 T cells and regulatory T cells prior to infections across strains. A basal propensity of T cells to express IFNg and IL17 over TNFa also correlated with early viral control. Overall, a dysregulated, pro-inflammatory signature of circulating T cells at baseline was associated with severe disease upon infection. While future studies of human samples prior to infection with SARS-CoV-2 are required, our studies in mice with SARS-CoV serve as proof of concept that circulating T cell signatures at baseline can predict clinical and virologic outcomes upon SARS-CoV infection. Identification of basal immune predictors in humans could allow for identification of individuals at highest risk of severe clinical and virologic outcomes upon infection, who may thus most benefit from available clinical interventions to restrict infection and disease.

SUMMARY

We used a screen of genetically diverse mice from the Collaborative Cross infected with mouse-adapted SARS-CoV in combination with comprehensive pre-infection immunophenotyping to identify baseline circulating immune correlates of severe virologic and clinical outcomes upon SARS-CoV infection.

Immune Predictors of Mortality After Ribonucleic Acid Virus Infection

BACKGROUND

Virus infections result in a range of clinical outcomes for the host, from asymptomatic to severe or even lethal disease. Despite global efforts to prevent and treat virus infections to limit morbidity and mortality, the continued emergence and re-emergence of new outbreaks as well as common infections such as influenza persist as a health threat. Challenges to the prevention of severe disease after virus infection include both a paucity of protective vaccines as well as the early identification of individuals with the highest risk that may require supportive treatment.

METHODS

We completed a screen of mice from the Collaborative Cross (CC) that we infected with influenza, severe acute respiratory syndrome-coronavirus, and West Nile virus.

RESULTS

The CC mice exhibited a range of disease manifestations upon infections, and we used this natural variation to identify strains with mortality after infection and strains exhibiting no mortality. We then used comprehensive preinfection immunophenotyping to identify global baseline immune correlates of protection from mortality to virus infection.

CONCLUSIONS

These data suggest that immune phenotypes might be leveraged to identify humans at highest risk of adverse clinical outcomes upon infection, who may most benefit from intensive clinical interventions, in addition to providing insight for rational vaccine design.

Genetic regulation of immune homeostatic lung leukocyte populations

Immune homeostasis is the state where the immune system maintains stability in the absence of insult. Much of the analysis of immune homeostasis has focused on systemic immunity, but it is also likely to be important in an organ specific manner. There is evidence that homeostatic immunity can affect subsequent responses to infection or vaccination. Since the lungs are a major site of infection, we used the Collaborative Cross (CC) mouse genetic reference population to study the genetic regulation of the breadth of baseline immune cell populations in the lung and identify loci regulating these cells at the steady state. We found that all immune cell populations measured showed strong genetic (i.e. strain-specific) variation in cell type abundances. We identified 12 quantitative trait loci (QTL) associated with variation in 12 immune cell populations or the relationships between cell populations. Given the role of various immune cells in the lungs during respiratory virus pathogenesis, we asked whether any of the mapped QTL correlated with influenza A virus (IAV) or Severe acute respiratory syndrome associated coronavirus (SARS-CoV) disease following infection in the same strains of mice. Notably, a locus we mapped for baseline abundance of CD8+ T cells in the lungs was associated with peak weight loss following IAV infection. Additionally, a locus mapped for variation in Ly6C+ monocyte/macrophage abundance was associated with SARS-CoV titer at days 2 and 4 post-infection. These data suggest that abundance of lung leukocyte populations prior to infection could serve as predictors of immune responses to respiratory viruses.

Complex Genetic Architecture Underlies Regulation of Influenza-A-Virus-Specific Antibody Responses in the Collaborative Cross

Host genetic factors play a fundamental role in regulating humoral immunity to viral infection, including influenza A virus (IAV). Here, we utilize the Collaborative Cross (CC), a mouse genetic reference population, to study genetic regulation of variation in antibody response following IAV infection. CC mice show significant heritable variation in the magnitude, kinetics, and composition of IAV-specific antibody response. We map 23 genetic loci associated with this variation. Analysis of a subset of these loci finds that they broadly affect the antibody response to IAV as well as other viruses. Candidate genes are identified based on predicted variant consequences and haplotype-specific expression patterns, and several show overlap with genes identified in human mapping studies. These findings demonstrate that the host antibody response to IAV infection is under complex genetic control and highlight the utility of the CC in modeling and identifying genetic factors with translational relevance to human health and disease.

Whole Genome Sequencing and Progress Toward Full Inbreeding of the Mouse Collaborative Cross Population

Two key features of recombinant inbred panels are well-characterized genomes and reproducibility. Here we report on the sequenced genomes of six additional Collaborative Cross (CC) strains and on inbreeding progress of 72 CC strains. We have previously reported on the sequences of 69 CC strains that were publicly available, bringing the total of CC strains with whole genome sequence up to 75. The sequencing of these six CC strains updates the efforts toward inbreeding undertaken by the UNC Systems Genetics Core. The timing reflects our competing mandates to release to the public as many CC strains as possible while achieving an acceptable level of inbreeding. The new six strains have a higher than average founder contribution from non-domesticus strains than the previously released CC strains. Five of the six strains also have high residual heterozygosity (>14%), which may be related to non-domesticus founder contributions. Finally, we report on updated estimates on residual heterozygosity across the entire CC population using a novel, simple and cost effective genotyping platform on three mice from each strain. We observe a reduction in residual heterozygosity across all previously released CC strains. We discuss the optimal use of different genetic resources available for the CC population.

Age-dependent Susceptibility To Virus Infection Is Governed by Host Genetics

While improved sanitation and medical care have reduced the burden of infectious disease, it remains the 4th leading cause of death in the elderly. Coupled with the expansion of the aged population over the next two decades, understanding the interplay between infection and the senescent host is critical. Both the microbiome and host genetics have been predicted to play critical roles in aging immunity. In this work, we explore the role of each in age-dependent susceptibility to virus infection. Severe acute respiratory syndrome coronavirus (SARS-CoV), the first outbreak virus of the 21st century, produces increased disease in aged patients as compared to the young. Importantly, age-dependent susceptibility is observed in mouse models of SARS-CoV infection. Our results indicate that despite microbiome variation between young and aged animals, limiting these differences had only a marginal effect on disease in the aged following infection. In contrast, age-dependent susceptibility was not fully conserved across genetically diverse Collaborative Cross (CC) mice. We found that while the majority of CC lines maintained increased disease as they aged, several lines showed no distinction between young and aged animals. In one CC line, we found that aged mice were more resistant to SARS-CoV infection than their younger equivalent. Notably, age-dependent susceptibility was not simply associated with higher viral loads driving augmented disease. Instead, the results suggest that host factors including tissue and inflammation responses govern disease in the senescent host. Together, the work highlights the importance of genetic elements in driving age-dependent disease following infection.

The Collaborative Cross: A Systems Genetics Resource for Studying Host-Pathogen Interactions

Host genetic variation has a major impact on infectious disease susceptibility. The study of pathogen resistance genes, largely aided by mouse models, has significantly advanced our understanding of infectious disease pathogenesis. The Collaborative Cross (CC), a newly developed multi-parental mouse genetic reference population, serves as a tractable model system to study how pathogens interact with genetically diverse populations. In this review, we summarize progress utilizing the CC as a platform to develop improved models of pathogen-induced disease and to map polymorphic host response loci associated with variation in susceptibility to pathogens.

Genetic diversity between mouse strains allows identification of the CC027/GeniUnc strain as an orally reactive model of peanut allergy

BACKGROUND

Improved animal models are needed to understand the genetic and environmental factors that contribute to food allergy.

OBJECTIVE

We sought to assess food allergy phenotypes in a genetically diverse collection of mice.

METHODS

We selected 16 Collaborative Cross (CC) mouse strains, as well as the classic inbred C57BL/6J, C3H/HeJ, and BALB/cJ strains, for screening. Female mice were sensitized to peanut intragastrically with or without cholera toxin and then challenged with peanut by means of oral gavage or intraperitoneal injection and assessed for anaphylaxis. Peanut-specific immunoglobulins, T-cell cytokines, regulatory T cells, mast cells, and basophils were quantified.

RESULTS

Eleven of the 16 CC strains had allergic reactions to intraperitoneal peanut challenge, whereas only CC027/GeniUnc mice reproducibly experienced severe symptoms after oral food challenge (OFC). CC027/GeniUnc, C3H/HeJ, and C57BL/6J mice all mounted a T_H2 response against peanut, leading to production of IL-4 and IgE, but only the CC027/GeniUnc mice reacted to OFC. Orally induced anaphylaxis in CC027/GeniUnc mice was correlated with serum levels of Ara h 2 in circulation but not with allergen-specific IgE or mucosal mast cell protease 1 levels, indicating systemic allergen absorption is important for anaphylaxis through the gastrointestinal tract. Furthermore, CC027/GeniUnc, but not C3H/HeJ or BALB/cJ, mice can be sensitized in the absence of cholera toxin and react on OFC to peanut.

CONCLUSIONS

We have identified and characterized CC027/GeniUnc mice as a strain that is genetically susceptible to peanut allergy and prone to severe reactions after OFC. More broadly, these findings demonstrate the untapped potential of the CC population in developing novel models for allergy research.

Extensive Homeostatic T Cell Phenotypic Variation within the Collaborative Cross

The Collaborative Cross (CC) is a panel of reproducible recombinant inbred mouse strains with high levels of standing genetic variation, affording an unprecedented opportunity to perform experiments in a small animal model containing controlled genetic diversity while allowing for genetic replicates. Here, we advance the utility of this unique mouse resource for immunology research because it allows for both examination and genetic dissection of mechanisms behind adaptive immune states in mice with distinct and defined genetic makeups. This approach is based on quantitative trait locus mapping: identifying genetically variant genome regions associated with phenotypic variance in traits of interest. Furthermore, the CC can be utilized for mouse model development; distinct strains have unique immunophenotypes and immune properties, making them suitable for research on particular diseases and infections. Here, we describe variations in cellular immune phenotypes across F1 crosses of CC strains and reveal quantitative trait loci responsible for several immune phenotypes.

•

The Collaborative Cross models the phenotypic diversity observed in human immunity

•

QTL mapping in the CC reveals candidate genes linked to T cell phenotypes

CC002/Unc females are mouse models of exercise-induced paradoxical fat response

Exercise results in beneficial health outcomes and protects against a variety of chronic diseases. However, U.S. exercise guidelines recommend identical exercise programs for everyone, despite individual variation in responses to these programs, including paradoxical fat gain. Experimental models of exercise-induced paradoxical outcomes may enable the dissection of underlying physiological mechanisms as well as the evaluation of potential interventions. Whereas several studies have identified individual mice exhibiting paradoxical fat gain following exercise, no systematic effort has been conducted to identify and characterize models of paradoxical response. Strains from the Collaborative Cross (CC) genetic reference population were used due to its high levels of genetic variation, its reproducible nature, and the observation that the CC is a rich source of novel disease models, to assess the impact genetic background has on exercise responses. We identified the strain CC002/Unc as an exercise-induced paradoxical fat response model in a controlled voluntary exercise study across multiple ages in female mice. We also found sex and genetic differences were consistent with this pattern in a study of forced exercise programs. These results provide a novel model for studies to determine the mechanisms behind paradoxical metabolic responses to exercise, and enable development of more rational personalized exercise recommendations based on factors such as age, sex, and genetic background.

Methyl-CpG binding domain protein-1 polymorphisms impact immune homeostatic levels of IgG

Host immunity plays a key role in protection against infection, however, there is significant individual-to-individual variation in the induction and quality of immune response that is driven in part by genetic factors. A growing body of evidence suggests that variation in immune status prior to specific insult may influence subsequent response to vaccination or pathogen challenge. Therefore, understanding the variation in immune status prior to specific insult is necessary for assessing variability in immune response to vaccination or pathogen challenge. The Collaborative Cross (CC) mouse population is a genetically diverse population that was designed to mimic the genetic and phenotypic diversity that is observed in humans and allows us to investigate mechanistically the contribution of genetics to complex phenotypes. Consequently, we asked whether there was variation in antibody levels in 57 CC strains prior to pathogen challenge and if that variation was associated with specific genetic factors. To assess this, we measured antibody levels for IgM, total IgG, IgG1, IgG2a, IgG2b, and IgG2c by ELISA and mapped quantitative trait loci (QTL) associated with variation in antibody concentrations. We observed substantial variation in all antibody isotypes measured. We also identified a suggestive (p=0.1) QTL on chromosome 18. Allele effects show that the C57Bl/6J allele at the locus of interest is linked to higher levels of total IgG. Subsequent analysis of genes under the QTL identified Methyl-CpG binding domain protein-1 (Mbd1) as the most promising candidate gene associated with total IgG variation. Ongoing studies are evaluating the role of Mbd1 in the regulation of total IgG production.