Search

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

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.

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.

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

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

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

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.

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.

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 T cell receptor engagement. Whether bystander T cell activation is also an important defense mechanism in the mucosa is poorly understood. Here, we investigated whether 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-nonspecific 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 contributes a fast and effective innate-like response to infection in mucosal tissue.

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.

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.

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.

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

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.

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.

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.

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

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.

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.

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.

Summary: 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. : Noll et al. use the Collaborative Cross, a mouse genetic reference population, to map genetic loci associated with variation in the humoral response to influenza virus infection. Cross-dataset comparison shows that mapped loci are important for antibody response to multiple pathogens, and candidate genes with likely translational relevance are identified. Keywords: Collaborative Cross, humoral immunity, influenza, influenza virus, antibody, host genetics, genetic architecture, genetic mapping, genetic reference population, complex trait

ABSTRACT Host genetics plays an important role in determining the outcome of Mycobacterium tuberculosis infection. We previously found that Collaborative Cross (CC) mouse strains differ in their susceptibility to M. tuberculosis and that the CC042/GeniUnc (CC042) strain suffered from a rapidly progressive disease and failed to produce the protective cytokine gamma interferon (IFN-γ) in the lung. Here, we used parallel genetic and immunological approaches to investigate the basis of CC042 mouse susceptibility. Using a population derived from a CC001/Unc (CC001) × CC042 intercross, we mapped four quantitative trait loci (QTL) underlying tuberculosis immunophenotypes (Tip1 to Tip4). These included QTL that were associated with bacterial burden, IFN-γ production following infection, and an IFN-γ-independent mechanism of bacterial control. Further immunological characterization revealed that CC042 animals recruited relatively few antigen-specific T cells to the lung and that these T cells failed to express the integrin alpha L (αL; i.e., CD11a), which contributes to T cell activation and migration. These defects could be explained by a CC042 private variant in the Itgal gene, which encodes CD11a and is found within the Tip2 interval. This 15-bp deletion leads to aberrant mRNA splicing and is predicted to result in a truncated protein product. The ItgalCC042 genotype was associated with all measured disease traits, indicating that this variant is a major determinant of susceptibility in CC042 mice. The combined effect of functionally distinct Tip variants likely explains the profound susceptibility of CC042 mice and highlights the multigenic nature of tuberculosis control in the Collaborative Cross. IMPORTANCE The variable outcome of Mycobacterium tuberculosis infection observed in natural populations is difficult to model in genetically homogeneous small-animal models. The newly developed Collaborative Cross (CC) represents a reproducible panel of genetically diverse mice that display a broad range of phenotypic responses to infection. We explored the genetic basis of this variation, focusing on a CC line that is highly susceptible to M. tuberculosis infection. This study identified multiple quantitative trait loci associated with bacterial control and cytokine production, including one that is caused by a novel loss-of-function mutation in the Itgal gene, which is necessary for T cell recruitment to the infected lung. These studies verify the multigenic control of mycobacterial disease in the CC panel, identify genetic loci controlling diverse aspects of pathogenesis, and highlight the utility of the CC resource.

Influenza A virus (IAV) is a respiratory pathogen that causes substantial morbidity and mortality during both seasonal and pandemic outbreaks. Infection outcomes in unexposed populations are affected by host genetics, but the host genetic architecture is not well understood. Here, we obtain a broad view of how heritable factors affect a mouse model of response to IAV infection using an 8 × 8 diallel of the eight inbred founder strains of the Collaborative Cross (CC). Expanding on a prior statistical framework for modeling treatment response in diallels, we explore how a range of heritable effects modify acute host response to IAV through 4 d postinfection. Heritable effects in aggregate explained ∼57% of the variance in IAV-induced weight loss. Much of this was attributable to a pattern of additive effects that became more prominent through day 4 postinfection and was consistent with previous reports of antiinfluenza myxovirus resistance 1 (Mx1) polymorphisms segregating between these strains; these additive effects largely recapitulated haplotype effects observed at the Mx1 locus in a previous study of the incipient CC, and are also replicated here in a CC recombinant intercross population. Genetic dominance of protective Mx1 haplotypes was observed to differ by subspecies of origin: relative to the domesticus null Mx1 allele, musculus acts dominantly whereas castaneus acts additively. After controlling for Mx1, heritable effects, though less distinct, accounted for ∼34% of the phenotypic variance. Implications for future mapping studies are discussed.

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.

Host genetic variation is known to contribute to differential pathogenesis following infection. Mouse models allow direct assessment of host genetic factors responsible for susceptibility to Severe Acute Respiratory Syndrome coronavirus (SARS-CoV). Based on an assessment of early stage lines from the Collaborative Cross mouse multi-parent population, we identified two lines showing highly divergent susceptibilities to SARS-CoV: the resistant CC003/Unc and the susceptible CC053/Unc. We generated 264 F2 mice between these strains, and infected them with SARS-CoV. Weight loss, pulmonary hemorrhage, and viral load were all highly correlated disease phenotypes. We identified a quantitative trait locus of major effect on chromosome 18 (27.1–58.6 Mb) which affected weight loss, viral titer and hemorrhage. Additionally, each of these three phenotypes had distinct quantitative trait loci [Chr 9 (weight loss), Chrs 7 and 12 (virus titer), and Chr 15 (hemorrhage)]. We identified Ticam2, an adaptor protein in the TLR signaling pathways, as a candidate driving differential disease at the Chr 18 locus. Ticam2−/− mice were highly susceptible to SARS-CoV infection, exhibiting increased weight loss and more pulmonary hemorrhage than control mice. These results indicate a critical role for Ticam2 in SARS-CoV disease, and highlight the importance of host genetic variation in disease responses.

The oligoadenylate-synthetase (Oas) gene locus provides innate immune resistance to virus infection. In mouse models, variation in the Oas1b gene influences host susceptibility to flavivirus infection. However, the impact of Oas variation on overall innate immune programming and global gene expression among tissues and in different genetic backgrounds has not been defined. We examined how Oas1b acts in spleen and brain tissue to limit West Nile virus (WNV) susceptibility and disease across a range of genetic backgrounds. The laboratory founder strains of the mouse Collaborative Cross (CC) (A/J, C57BL/6J, 129S1/SvImJ, NOD/ShiLtJ, and NZO/HlLtJ) all encode a truncated, defective Oas1b, whereas the three wild-derived inbred founder strains (CAST/EiJ, PWK/PhJ, and WSB/EiJ) encode a full-length OAS1B protein. We assessed disease profiles and transcriptional signatures of F1 hybrids derived from these founder strains. F1 hybrids included wild-type Oas1b (F/F), homozygous null Oas1b (N/N), and heterozygous offspring of both parental combinations (F/N and N/F). These mice were challenged with WNV, and brain and spleen samples were harvested for global gene expression analysis. We found that the Oas1b haplotype played a role in WNV susceptibility and disease metrics, but the presence of a functional Oas1b allele in heterozygous offspring did not absolutely predict protection against disease. Our results indicate that Oas1b status as wild-type or truncated, and overall Oas1b gene dosage, link with novel innate immune gene signatures that impact specific biological pathways for the control of flavivirus infection and immunity through both Oas1b-dependent and independent processes.

The success of lentiviral vectors in curing fatal genetic and acquired diseases has opened a new era in human gene therapy. However, variability in the efficacy and safety of this therapeutic approach has been reported in human patients. Consequently, lentiviral-vector-based gene therapy is limited to incurable human diseases, with little understanding of the underlying causes of adverse effects and poor efficacy. To assess the role that host genetic variation has on efficacy of gene therapy, we characterized lentiviral-vector gene therapy within a set of 12 collaborative cross mouse strains. Lentiviral vectors carrying the firefly luciferase cDNA under the control of a liver-specific promoter were administered to female mice, with total-body and hepatic luciferase expression periodically monitored through 41 weeks post-vector administration. Vector copy number per diploid genome in mouse liver and spleen was determined at the end of this study. We identified major strain-specific contributions to overall success of transduction, vector biodistribution, maximum luciferase expression, and the kinetics of luciferase expression throughout the study. Our results highlight the importance of genetic variation on gene-therapeutic efficacy; provide new models with which to more rigorously assess gene therapy approaches; and suggest that redesigning preclinical studies of gene-therapy methodologies might be appropriate.

West Nile Virus (WNV) is a mosquito-transmitted virus from the Flaviviridae family that causes fever in 1 in 5 infected people. WNV can also become neuro-invasive and cross the blood-brain barrier leading to severe neurological symptoms in a subset of WNV infected individuals [1]. WNV neuro-invasion is believed to be influenced by a number of factors including host genetics. In order to explore these effects and recapitulate the complex immune genetic differences among individuals, we studied gene expression following WNV infection in the Collaborative Cross (CC) model. The CC is a mouse genetics resource composed of >70 independently bred, octo-parental recombinant inbred mouse lines [2]. To identify the individual host gene expression signatures influencing protection or susceptibility to WNV disease and WNV neuroinvasion, we used the nanostring nsolver platform to quantify gene expression in brain tissue isolated from WNV-infected CC mice at days 4, 7 and 12 post-infection [3]. This nanostring technology provided a high throughput, non-amplification based mRNA quantitation method to detect immune genes involved in neuro-invasion. Data was deposited into the Gene Expression Omnibus (GEO) under accession GSE85999.

Flaviviruses are hematophagous arthropod-viruses that pose global challenges to human health. Like Zika virus, West Nile Virus (WNV) is a flavivirus for which no approved vaccine exists [1]. The role host genetics play in early detection and response to WNV still remains largely unexplained. In order to capture the impact of genetic variation on innate immune responses, we studied gene expression following WNV infection using the collaborative cross (CC). The CC is a mouse genetics resource composed of hundreds of independently bred, octo-parental recombinant inbred mouse lines [2]. To accurately capture the host immune gene expression signatures of West Nile infection, we used the nanostring platform to evaluate expression in spleen tissue isolated from CC mice infected with WNV over a time course of 4, 7, and 12 days' post-infection [3]. Nanostring is a non-amplification based digital method to quantitate gene expression that uses color-coded molecular barcodes to detect hundreds of transcripts in a sample. Using this approach, we identified unique gene signatures in spleen tissue at days 4, 7, and 12 following WNV infection, which delineated distinct differences between asymptomatic and symptomatic CC lines. We also identified novel immune genes. Data was deposited into the Gene Expression Omnibus under accession GSE86000.

Genotyping microarrays are an important resource for genetic mapping, population genetics, and monitoring of the genetic integrity of laboratory stocks. We have developed the third generation of the Mouse Universal Genotyping Array (MUGA) series, GigaMUGA, a 143,259-probe Illumina Infinium II array for the house mouse (Mus musculus). The bulk of the content of GigaMUGA is optimized for genetic mapping in the Collaborative Cross and Diversity Outbred populations, and for substrain-level identification of laboratory mice. In addition to 141,090 single nucleotide polymorphism probes, GigaMUGA contains 2006 probes for copy number concentrated in structurally polymorphic regions of the mouse genome. The performance of the array is characterized in a set of 500 high-quality reference samples spanning laboratory inbred strains, recombinant inbred lines, outbred stocks, and wild-caught mice. GigaMUGA is highly informative across a wide range of genetically diverse samples, from laboratory substrains to other Mus species. In addition to describing the content and performance of the array, we provide detailed probe-level annotation and recommendations for quality control.

Infection with West Nile virus (WNV) leads to a range of disease outcomes, including chronic infection, though lack of a robust mouse model of chronic WNV infection has precluded identification of the immune events contributing to persistent infection. Using the Collaborative Cross, a population of recombinant inbred mouse strains with high levels of standing genetic variation, we have identified a mouse model of persistent WNV disease, with persistence of viral loads within the brain. Compared to lines exhibiting no disease or marked disease, the F1 cross CC(032x013)F1 displays a strong immunoregulatory signature upon infection that correlates with restraint of the WNV-directed cytolytic response. We hypothesize that this regulatory T cell response sufficiently restrains the immune response such that a chronic infection can be maintained in the CNS. Use of this new mouse model of chronic neuroinvasive virus will be critical in developing improved strategies to prevent prolonged disease in humans.

New systems genetics approaches are needed to rapidly identify host genes and genetic networks that regulate complex disease outcomes. Using genetically diverse animals from incipient lines of the Collaborative Cross mouse panel, we demonstrate a greatly expanded range of phenotypes relative to classical mouse models of SARS-CoV infection including lung pathology, weight loss and viral titer. Genetic mapping revealed several loci contributing to differential disease responses, including an 8.5Mb locus associated with vascular cuffing on chromosome 3 that contained 23 genes and 13 noncoding RNAs. Integrating phenotypic and genetic data narrowed this region to a single gene, Trim55, an E3 ubiquitin ligase with a role in muscle fiber maintenance. Lung pathology and transcriptomic data from mice genetically deficient in Trim55 were used to validate its role in SARS-CoV-induced vascular cuffing and inflammation. These data establish the Collaborative Cross platform as a powerful genetic resource for uncovering genetic contributions of complex traits in microbial disease severity, inflammation and virus replication in models of outbred populations.

ABSTRACT West Nile virus (WNV) is an emerging neuroinvasive flavivirus that now causes significant morbidity and mortality worldwide. The innate and adaptive immune responses to WNV infection have been well studied in C57BL/6J inbred mice, but this model lacks the variations in susceptibility, immunity, and outcome to WNV infection that are observed in humans, thus limiting its usefulness to understand the mechanisms of WNV infection and immunity dynamics. To build a model of WNV infection that captures human infection outcomes, we have used the Collaborative Cross (CC) mouse model. We show that this model, which recapitulates the genetic diversity of the human population, demonstrates diversity in susceptibility and outcomes of WNV infection observed in humans. Using multiple F1 crosses of CC mice, we identified a wide range of susceptibilities to infection, as demonstrated through differences in survival, clinical disease score, viral titer, and innate and adaptive immune responses in both peripheral tissues and the central nervous system. Additionally, we examined the Oas1b alleles in the CC mice and confirmed the previous finding that Oas1b plays a role in susceptibility to WNV; however, even within a given Oas1b allele status, we identified a wide range of strain-specific WNV-associated phenotypes. These results confirmed that the CC model is effective for identifying a repertoire of host genes involved in WNV resistance and susceptibility. The CC effectively models a wide range of WNV clinical, virologic, and immune phenotypes, thus overcoming the limitations of the traditional C57BL/6J model, allowing genetic and mechanistic studies of WNV infection and immunity in differently susceptible populations. IMPORTANCE Mouse models of West Nile virus infection have revealed important details regarding the innate and adaptive immune responses to this emerging viral infection. However, traditional mouse models lack the genetic diversity present in human populations and therefore limit our ability to study various disease outcomes and immunologic mechanisms subsequent to West Nile virus infection. In this study, we used the Collaborative Cross mouse model to more effectively model the wide range of clinical, virologic, and immune phenotypes present upon West Nile virus infection in humans.