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

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.

A susceptibility gene signature for ERBB2-driven mammary tumour development and metastasis in collaborative cross mice

BACKGROUND

Deeper insights into ERBB2-driven cancers are essential to develop new treatment approaches for ERBB2+ breast cancers (BCs). We employed the Collaborative Cross (CC) mouse model to unearth genetic factors underpinning Erbb2-driven mammary tumour development and metastasis.

METHODS

732 F1 hybrid female mice between FVB/N MMTV-Erbb2 and 30 CC strains were monitored for mammary tumour phenotypes. GWAS pinpointed SNPs that influence various tumour phenotypes. Multivariate analyses and models were used to construct the polygenic score and to develop a mouse tumour susceptibility gene signature (mTSGS), where the corresponding human ortholog was identified and designated as hTSGS. The importance and clinical value of hTSGS in human BC was evaluated using public datasets, encompassing TCGA, METABRIC, GSE96058, and I-SPY2 cohorts. The predictive power of mTSGS for response to chemotherapy was validated in vivo using genetically diverse MMTV-Erbb2 mice.

FINDINGS

Distinct variances in tumour onset, multiplicity, and metastatic patterns were observed in F1-hybrid female mice between FVB/N MMTV-Erbb2 and 30 CC strains. Besides lung metastasis, liver and kidney metastases emerged in specific CC strains. GWAS identified specific SNPs significantly associated with tumour onset, multiplicity, lung metastasis, and liver metastasis. Multivariate analyses flagged SNPs in 20 genes (Stx6, Ramp1, Traf3ip1, Nckap5, Pfkfb2, Trmt1l, Rprd1b, Rer1, Sepsecs, Rhobtb1, Tsen15, Abcc3, Arid5b, Tnr, Dock2, Tti1, Fam81a, Oxr1, Plxna2, and Tbc1d31) independently tied to various tumour characteristics, designated as a mTSGS. hTSGS scores (hTSGSS) based on their transcriptional level showed prognostic values, superseding clinical factors and PAM50 subtype across multiple human BC cohorts, and predicted pathological complete response independent of and superior to MammaPrint score in I-SPY2 study. The power of mTSGS score for predicting chemotherapy response was further validated in an in vivo mouse MMTV-Erbb2 model, showing that, like findings in human patients, mouse tumours with low mTSGS scores were most likely to respond to treatment.

INTERPRETATION

Our investigation has unveiled many new genes predisposing individuals to ERBB2-driven cancer. Translational findings indicate that hTSGS holds promise as a biomarker for refining treatment strategies for patients with BC.

FUNDING

The U.S. Department of Defense (DoD) Breast Cancer Research Program (BCRP) (BC190820), United States; MCIN/AEI/10.13039/501100011039 (PID2020-118527RB-I00, PDC2021-121735-I00), the "European Union Next Generation EU/PRTR," the Regional Government of Castile and León (CSI144P20), European Union.

Characterization of Collaborative Cross mouse founder strain CAST/EiJ as a novel model for lethal COVID-19

Mutations in SARS-CoV-2 variants of concern (VOCs) have expanded the viral host range beyond primates, and a limited range of other mammals, to mice, affording the opportunity to exploit genetically diverse mouse panels to model the broad range of responses to infection in patient populations. Here we surveyed responses to VOC infection in genetically diverse Collaborative Cross (CC) founder strains. Infection of wild-derived CC founder strains produced a broad range of viral burden, disease susceptibility and survival, whereas most other strains were resistant to disease despite measurable lung viral titers. In particular, CAST/EiJ, a wild-derived strain, developed high lung viral burdens, more severe lung pathology than seen in other CC strains, and a dysregulated cytokine profile resulting in morbidity and mortality. These inbred mouse strains may serve as a valuable platform to evaluate therapeutic countermeasures against severe COVID-19 and other coronavirus pandemics in the future.

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.

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.

Immunogenicity of an adjuvanted broadly active influenza vaccine in immunocompromised and diverse populations

Influenza virus outbreaks are a major burden worldwide each year. Current vaccination strategies are inadequate due to antigenic drift/shift of the virus and the elicitation of low immune responses. The use of computationally optimized broadly reactive antigen (COBRA) hemagglutinin (HA) immunogens subvert the constantly mutating viruses; however, they are poorly immunogenic on their own. To increase the immunogenicity of subunit vaccines such as this, adjuvants can be delivered with the vaccine. For example, agonists of the stimulator of interferon genes (STING) have proven efficacy as vaccine adjuvants. However, their use in high-risk populations most vulnerable to influenza virus infection has not been closely examined. Here, we utilize a vaccine platform consisting of acetalated dextran microparticles loaded with COBRA HA and the STING agonist cyclic GMP-AMP. We examine the immunogenicity of this platform in mouse models of obesity, aging, and chemotherapy-induced immunosuppression. Further, we examine vaccine efficacy in collaborative cross mice, a genetically diverse population that mimics human genetic heterogeneity. Overall, this vaccine platform had variable efficacy in these populations supporting work to better tailor adjuvants to specific populations.

Host Genetic Background Influences BCG-Induced Antibodies Cross-Reactive to SARS-CoV-2 Spike Protein

Mycobacterium bovis Bacillus Calmette-Guérin (BCG) protects against childhood tuberculosis; and unlike most vaccines, BCG broadly impacts immunity to other pathogens and even some cancers. Early in the COVID-19 pandemic, epidemiological studies identified a protective association between BCG vaccination and outcomes of SARS-CoV-2, but the associations in later studies were inconsistent. We sought possible reasons and noticed the study populations often lived in the same country. Since individuals from the same regions can share common ancestors, we hypothesized that genetic background could influence associations between BCG and SARS-CoV-2. To explore this hypothesis in a controlled environment, we performed a pilot study using Diversity Outbred mice. First, we identified amino acid sequences shared by BCG and SARS-CoV-2 spike protein. Next, we tested for IgG reactive to spike protein from BCG-vaccinated mice. Sera from some, but not all, BCG-vaccinated Diversity Outbred mice contained higher levels of IgG cross-reactive to SARS-CoV-2 spike protein than sera from BCG-vaccinated C57BL/6J inbred mice and unvaccinated mice. Although larger experimental studies are needed to obtain mechanistic insight, these findings suggest that genetic background may be an important variable contributing to different associations observed in human randomized clinical trials evaluating BCG vaccination on SARS-CoV-2 and COVID-19.

Mechanistic insights into bone remodelling dysregulation by human viral pathogens

Bone-related diseases (osteopathologies) associated with human virus infections have increased around the globe. Recent findings have highlighted the intricate interplay between viral infection, the host immune system and the bone remodelling process. Viral infections can disrupt bone homeostasis, contributing to conditions such as arthritis and soft tissue calcifications. Osteopathologies can occur after arbovirus infections such as chikungunya virus, dengue virus and Zika virus, as well as respiratory viruses, such as severe acute respiratory syndrome coronavirus 2 and enteroviruses such as Coxsackievirus B. Here we explore how human viruses dysregulate bone homeostasis, detailing viral factors, molecular mechanisms, host immune response changes and bone remodelling that ultimately result in osteopathologies. We highlight model systems and technologies to advance mechanistic understanding of viral-mediated bone alterations. Finally, we propose potential prophylactic and therapeutic strategies, introduce 'osteovirology' as a research field highlighting the underestimated roles of viruses in bone-related diseases, and discuss research avenues for further investigation.

Genetic architecture of the acute and persistent immune cell response after radiation exposure

Hematologic toxicity is a common side effect of multimodal cancer therapy. Nearly all animal studies investigating the causes of radiotherapy-induced hematologic toxicity use inbred strains with limited genetic diversity and do not reflect the diverse responses observed in humans. We used the population-based Collaborative Cross (CC) mouse resource to investigate the genetic architecture of the acute and persistent immune response after radiation exposure by measuring 22 immune parameters in 1,720 CC mice representing 35 strains. We determined relative acute and persistent radiation resistance scores at the individual strain level considering contributions from all immune parameters. Genome-wide association analysis identified quantitative trait loci associated with baseline and radiation responses. A cross-species radiation resistance score predicted recurrence-free survival in medulloblastoma patients. We present a community resource of immune parameters and genome-wide association analyses before and after radiation exposure for future investigations of the contributions of host genetics on radiosensitivity.

Functionally distinct regions of the locus Leishmania major response 15 control IgE or IFNγ level in addition to skin lesions

Leishmaniasis, a disease caused by parasites of Leishmania spp., endangers more than 1 billion people living in endemic countries and has three clinical forms: cutaneous, mucocutaneous, and visceral. Understanding of individual differences in susceptibility to infection and heterogeneity of its pathology is largely lacking. Different mouse strains show a broad and heterogeneous range of disease manifestations such as skin lesions, splenomegaly, hepatomegaly, and increased serum levels of immunoglobulin E and several cytokines. Genome-wide mapping of these strain differences detected more than 30 quantitative trait loci (QTLs) that control the response to Leishmania major. Some control different combinations of disease manifestations, but the nature of this heterogeneity is not yet clear. In this study, we analyzed the L. major response locus Lmr15 originally mapped in the strain CcS-9 which carries 12.5% of the genome of the resistant strain STS on the genetic background of the susceptible strain BALB/c. For this analysis, we used the advanced intercross line K3FV between the strains BALB/c and STS. We confirmed the previously detected loci Lmr15, Lmr18, Lmr24, and Lmr27 and performed genetic dissection of the effects of Lmr15 on chromosome 11. We prepared the interval-specific recombinant strains 6232HS1 and 6229FUD, carrying two STS-derived segments comprising the peak linkage of Lmr15 whose lengths were 6.32 and 17.4 Mbp, respectively, and analyzed their response to L. major infection. These experiments revealed at least two linked but functionally distinct chromosomal regions controlling IFNγ response and IgE response, respectively, in addition to the control of skin lesions. Bioinformatics and expression analysis identified the potential candidate gene Top3a. This finding further clarifies the genetic organization of factors relevant to understanding the differences in the individual risk of disease.

Using the collaborative cross to identify the role of host genetics in defining the murine gut microbiome

BACKGROUND

The human gut microbiota is a complex community comprised of trillions of bacteria and is critical for the digestion and absorption of nutrients. Bacterial communities of the intestinal microbiota influence the development of several conditions and diseases. We studied the effect of host genetics on gut microbial composition using Collaborative Cross (CC) mice. CC mice are a panel of mice that are genetically diverse across strains, but genetically identical within a given strain allowing repetition and deeper analysis than is possible with other collections of genetically diverse mice.

RESULTS

16S rRNA from the feces of 167 mice from 28 different CC strains was sequenced and analyzed using the Qiime2 pipeline. We observed a large variance in the bacterial composition across CC strains starting at the phylum level. Using bacterial composition data, we identified 17 significant Quantitative Trait Loci (QTL) linked to 14 genera on 9 different mouse chromosomes. Genes within these intervals were analyzed for significant association with pathways and the previously known human GWAS database using Enrichr analysis and Genecards database. Multiple host genes involved in obesity, glucose homeostasis, immunity, neurological diseases, and many other protein-coding genes located in these regions may play roles in determining the composition of the gut microbiota. A subset of these CC mice was infected with Salmonella Typhimurium. Using infection outcome data, an increase in abundance of genus Lachnospiraceae and decrease in genus Parasutterella correlated with positive health outcomes after infection. Machine learning classifiers accurately predicted the CC strain and the infection outcome using pre-infection bacterial composition data from the feces.

CONCLUSION

Our study supports the hypothesis that multiple host genes influence the gut microbiome composition and homeostasis, and that certain organisms may influence health outcomes after S. Typhimurium infection. Video Abstract.

Mouse models for immuno-oncology

Realizing the clinical promise of cancer immunotherapy is hindered by gaps in our knowledge of in vivo mechanisms underlying treatment response as well as treatment limiting toxicity. Preclinical in vivo model systems and technologies are required to address these knowledge gaps and to surmount the challenges faced in the clinical application of immunotherapy. Mice are commonly used for basic and translational research to support development and testing of immune interventions, including for cancer. Here, we discuss the advantages and the limitations of current models as well as future developments.

Distinct gene programs underpinning disease tolerance and resistance in influenza virus infection

When challenged with an invading pathogen, the host-defense response is engaged to eliminate the pathogen (resistance) and to maintain health in the presence of the pathogen (disease tolerance). However, the identification of distinct molecular programs underpinning disease tolerance and resistance remained obscure. We exploited transcriptional and physiological monitoring across 33 mouse strains, during in vivo influenza virus infection, to identify two host-defense gene programs-one is associated with hallmarks of disease tolerance and the other with hallmarks of resistance. Both programs constitute generic responses in multiple mouse and human cell types. Our study describes the organizational principles of these programs and validates Arhgdia as a regulator of disease-tolerance states in epithelial cells. We further reveal that the baseline disease-tolerance state in peritoneal macrophages is associated with the pathophysiological response to injury and infection. Our framework provides a paradigm for the understanding of disease tolerance and resistance at the molecular level.

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.

Using the Collaborative Cross and Diversity Outbred Mice in Immunology

The Collaborative Cross (CC) and the Diversity Outbred (DO) stock mouse panels are the most powerful murine genetics tools available to the genetics community. Together, they combine the strength of inbred animal models with the diversity of outbred populations. Using the 63 CC strains or a panel of DO mice, each derived from the same 8 parental mouse strains, researchers can map genetic contributions to exceptionally complex immunological and infectious disease traits that would require far greater powering if performed by genome-wide association studies (GWAS) in human populations. These tools allow genes to be studied in heterozygous and homozygous states and provide a platform to study epistasis between interacting loci. Most importantly, once a quantitative phenotype is investigated and quantitative trait loci are identified, confirmatory genetic studies can be performed, which is often problematic using the GWAS approach. In addition, novel stable mouse models for immune phenotypes are often derived from studies utilizing the DO and CC mice that can serve as stronger model systems than existing ones in the field. The CC/DO systems have contributed to the fields of cancer immunology, autoimmunity, vaccinology, infectious disease, allergy, tissue rejection, and tolerance but have thus far been greatly underutilized. In this article, we present a recent review of the field and point out key areas of immunology that are ripe for further investigation and awaiting new CC/DO research projects. We also highlight some of the strong computational tools that have been developed for analyzing CC/DO genetic and phenotypic data. Additionally, we have formed a centralized community on the CyVerse infrastructure where immunogeneticists can utilize those software tools, collaborate with groups across the world, and expand the use of the CC and DO systems for investigating immunogenetic phenomena. © 2022 Wiley Periodicals LLC.

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.

Standing on the shoulders of mice

While inbred mice have informed most of what we know about the immune system in the modern era, they have clear limitations with respect to their ability to be informative regarding genetic heterogeneity or microbial influences. They have also not been very predictive as models of human disease or vaccination results. Although there are concerted attempts to compensate for these flaws, the rapid rise of human studies, driven by both technical and conceptual advances, promises to fill in these gaps, as well as provide direct information about human diseases and vaccination responses. Work on human immunity has already provided important additional perspectives on basic immunology such as the importance of clonal deletion to self-tolerance, and while many challenges remain, it seems inevitable that "the human model" will continue to inform us about the immune system and even allow for the discovery of new mechanisms.

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.

Genetically diverse mouse models of SARS-CoV-2 infection reproduce clinical variation and cytokine responses in COVID-19

Inflammation in response to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection drives severity of coronavirus disease 2019 (COVID-19), with effective versus dysregulated responses influenced by host genetics. To understand mechanisms of inflammation, animal models that reflect genetic diversity and clinical outcomes observed in humans are needed. We report a mouse panel comprising the diverse genetic backgrounds of the Collaborative Cross founder strains crossed to K18-hACE2 transgenic mice that confers high susceptibility to SARS-CoV-2. Infection of CC x K18-hACE2 F1 progeny resulted in a spectrum of weight loss, survival, viral replication kinetics, histopathology, and cytokine profiles, some of which were sex-specific. Importantly, survival was associated with early type I interferon (IFN) expression and a phased proinflammatory response distinct from mice with severe disease. Thus, dynamics of inflammatory responses observed in COVID-19 can be modeled in diverse mouse strains that provide a genetically tractable platform for understanding antiviral immunity and evaluating countermeasures.

ONE SENTENCE SUMMARY

Genetically diverse mice display a broad spectrum of clinically relevant responses to SARS-CoV-2 infection, reflecting variability in COVID-19 disease.

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.

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.

Host genetic control of natural killer cell diversity revealed in the Collaborative Cross

Natural killer (NK) cells are innate effectors armed with cytotoxic and cytokine-secreting capacities whose spontaneous antitumor activity is key to numerous immunotherapeutic strategies. However, current mouse models fail to mirror the extensive immune system variation that exists in the human population which may impact on NK cell-based therapies. We performed a comprehensive profiling of NK cells in the Collaborative Cross (CC), a collection of novel recombinant inbred mouse strains whose genetic diversity matches that of humans, thereby providing a unique and highly diverse small animal model for the study of immune variation. We demonstrate that NK cells from CC strains displayed a breadth of phenotypic and functional variation reminiscent of that reported for humans with regards to cell numbers, key marker expression, and functional capacities. We took advantage of the vast genetic diversity of the CC and identified nine genomic loci through quantitative trait locus mapping driving these phenotypic variations. SNP haplotype patterns and variant effect analyses identified candidate genes associated with lung NK cell numbers, frequencies of CD94⁺ NK cells, and expression levels of NKp46. Thus, we demonstrate that the CC represents an outstanding resource to study NK cell diversity and its regulation by host genetics.

In-Depth Analysis of Genetic Variation Associated with Severe West Nile Viral Disease

West Nile virus (WNV) is a mosquito-borne virus which causes symptomatic disease in a minority of infected humans. To identify novel genetic variants associated with severe disease, we utilized data from an existing case-control study of WNV and included population controls for an expanded analysis. We conducted imputation and gene-gene interaction analysis in the largest and most comprehensive genetic study conducted to date for West Nile neuroinvasive disease (WNND). Within the imputed West Nile virus dataset (severe cases n = 381 and asymptomatic/mild controls = 441), we found novel loci within the MCF.2 Cell Line Derived Transforming Sequence Like (MCF2L) gene (rs9549655 and rs2297192) through the individual loci analyses, although none reached statistical significance. Incorporating population controls from the Wisconsin Longitudinal Study on Aging (n = 9012) did not identify additional novel variants, a possible reflection of the cohort's inclusion of individuals who could develop mild or severe WNV disease upon infection. Many of the top gene-gene interaction results were intergenic, with currently undefined biological roles, highlighting the need for further investigation into these regions and other identified gene targets in severe WNND. Further studies including larger sample sizes and more diverse populations reflective of those at risk are needed to fully understand the genetic architecture of severe WNDD and provide guidance on viable targets for therapeutic and vaccine development.

A Mouse-Adapted SARS-CoV-2 Induces Acute Lung Injury and Mortality in Standard Laboratory Mice

The SARS-CoV-2 pandemic has caused extreme human suffering and economic harm. We generated and characterized a new mouse-adapted SARS-CoV-2 virus that captures multiple aspects of severe COVID-19 disease in standard laboratory mice. This SARS-CoV-2 model exhibits the spectrum of morbidity and mortality of COVID-19 disease as well as aspects of host genetics, age, cellular tropisms, elevated Th1 cytokines, and loss of surfactant expression and pulmonary function linked to pathological features of acute lung injury (ALI) and acute respiratory distress syndrome (ARDS). This model can rapidly access existing mouse resources to elucidate the role of host genetics, underlying molecular mechanisms governing SARS-CoV-2 pathogenesis, and the protective or pathogenic immune responses related to disease severity. The model promises to provide a robust platform for studies of ALI and ARDS to evaluate vaccine and antiviral drug performance, including in the most vulnerable populations (i.e., the aged) using standard laboratory mice.