26
|
Schäfer A, Leist SR, Gralinski LE, Martinez DR, Winkler ES, Okuda K, Hawkins PE, Gully KL, Graham RL, Scobey DT, Bell TA, Hock P, Shaw GD, Loome JF, Madden EA, Anderson E, Baxter VK, Taft-Benz SA, Zweigart MR, May SR, Dong S, Clark M, Miller DR, Lynch RM, Heise MT, Tisch R, Boucher RC, Pardo Manuel de Villena F, Montgomery SA, Diamond MS, Ferris MT, Baric RS. A Multitrait Locus Regulates Sarbecovirus Pathogenesis. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2022. [PMID: 35677067 DOI: 10.1101/2022.06.01.494461] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
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 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.
Collapse
|
27
|
Beierle JA, Yao EJ, Goldstein SI, Scotellaro JL, Sena KD, Averin O, Moody DE, Reilly CA, Emili A, Peltz G, Ferris MT, Bryant CD. A reduced complexity cross between BALB/c substrains identifies Zhx2 as a candidate gene underlying oxycodone metabolite brain concentration and state‐dependent learning of opioid reward. FASEB J 2022. [DOI: 10.1096/fasebj.2022.36.s1.r4573] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
28
|
Gaines CH, Schoenrock SA, Farrington J, Lee DF, Aponte-Collazo LJ, Shaw GD, Miller DR, Ferris MT, Pardo-Manuel de Villena F, Tarantino LM. Cocaine-Induced Locomotor Activation Differs Across Inbred Mouse Substrains. Front Psychiatry 2022; 13:800245. [PMID: 35599758 PMCID: PMC9120424 DOI: 10.3389/fpsyt.2022.800245] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/22/2021] [Accepted: 03/22/2022] [Indexed: 11/13/2022] Open
Abstract
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.
Collapse
|
29
|
Cruz Cisneros MC, Parotti B, Noll K, Bell TA, Hock P, Heise MT, Ferris MT. Genetic analysis of differential responses to adjuvanted influenza vaccination. THE JOURNAL OF IMMUNOLOGY 2022. [DOI: 10.4049/jimmunol.208.supp.64.08] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Abstract
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.
Collapse
|
30
|
Smith CM, Baker RE, Proulx MK, Mishra BB, Long JE, Park SW, Lee HN, Kiritsy MC, Bellerose MM, Olive AJ, Murphy KC, Papavinasasundaram K, Boehm FJ, Reames CJ, Meade RK, Hampton BK, Linnertz CL, Shaw GD, Hock P, Bell TA, Ehrt S, Schnappinger D, Pardo-Manuel de Villena F, Ferris MT, Ioerger TR, Sassetti CM. Host-pathogen genetic interactions underlie tuberculosis susceptibility in genetically diverse mice. eLife 2022; 11:74419. [PMID: 35112666 PMCID: PMC8846590 DOI: 10.7554/elife.74419] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2021] [Accepted: 01/27/2022] [Indexed: 11/21/2022] Open
Abstract
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.
Collapse
|
31
|
Beierle JA, Yao EJ, Goldstein SI, Scotellaro JL, Sena KD, Linnertz CA, Willits AB, Kader L, Young EE, Peltz G, Emili A, Ferris MT, Bryant CD. Genetic basis of thermal nociceptive sensitivity and brain weight in a BALB/c reduced complexity cross. Mol Pain 2022; 18:17448069221079540. [PMID: 35088629 PMCID: PMC8891926 DOI: 10.1177/17448069221079540] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Accepted: 01/20/2022] [Indexed: 11/30/2022] Open
Abstract
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.
Collapse
|
32
|
Fay EJ, Balla KM, Roach SN, Shepherd FK, Putri DS, Wiggen TD, Goldstein SA, Pierson MJ, Ferris MT, Thefaine CE, Tucker A, Salnikov M, Cortez V, Compton SR, Kotenko SV, Hunter RC, Masopust D, Elde NC, Langlois RA. Natural rodent model of viral transmission reveals biological features of virus population dynamics. J Exp Med 2021; 219:212940. [PMID: 34958350 PMCID: PMC8713297 DOI: 10.1084/jem.20211220] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2021] [Revised: 11/05/2021] [Accepted: 12/08/2021] [Indexed: 12/21/2022] Open
Abstract
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.
Collapse
|
33
|
Keele GR, Zhang T, Pham DT, Vincent M, Bell TA, Hock P, Shaw GD, Paulo JA, Munger SC, Pardo-Manuel de Villena F, Ferris MT, Gygi SP, Churchill GA. Regulation of protein abundance in genetically diverse mouse populations. CELL GENOMICS 2021; 1:100003. [PMID: 36212994 PMCID: PMC9536773 DOI: 10.1016/j.xgen.2021.100003] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/05/2020] [Revised: 04/01/2021] [Accepted: 05/10/2021] [Indexed: 12/12/2022]
Abstract
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.
Collapse
|
34
|
Schäfer A, Gralinski LE, Leist SR, Winkler ES, Hampton BK, Mooney MA, Jensen KL, Graham RL, Agnihothram S, Jeng S, Chamberlin S, Bell TA, Scobey DT, VanBlargan LA, Thackray LB, Hock P, Miller DR, Shaw GD, de Villena FPM, McWeeney SK, Montgomery SA, Diamond MS, Heise MT, Menachery VD, Ferris MT, Baric RS. Common Mechanism of SARS-CoV and SARS-CoV-2 Pathogenesis across Species. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2021:2021.05.14.444205. [PMID: 34013261 PMCID: PMC8132217 DOI: 10.1101/2021.05.14.444205] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
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.
Collapse
|
35
|
Hampton BK, Jensen KL, Whitmore AC, Gralinski LE, Leist SR, Linnertz CL, Maurizio P, Menachery VD, Morrison CR, Noll KE, Plante KS, Shäfer A, Shaw GD, West A, de Villena FPM, Baric RS, Heise MT, Ferris MT. Genetic regulation of immune homeostatic lung leukocyte populations influences respiratory virus induced disease in collaborative cross mice. THE JOURNAL OF IMMUNOLOGY 2021. [DOI: 10.4049/jimmunol.206.supp.24.05] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
Abstract
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.
Collapse
|
36
|
Graham JB, Swarts JL, Edwards KR, Voss KM, Green R, Jeng S, Miller DR, Mooney MA, McWeeney SK, Ferris MT, Pardo-Manuel de Villena F, Gale M, Lund JM. Correlation of Regulatory T Cell Numbers with Disease Tolerance upon Virus Infection. Immunohorizons 2021; 5:157-169. [PMID: 33893179 PMCID: PMC8281504 DOI: 10.4049/immunohorizons.2100009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2021] [Accepted: 02/19/2021] [Indexed: 11/19/2022] Open
Abstract
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.
Collapse
|
37
|
Smeekens JM, Johnson-Weaver BT, Hinton AL, Azcarate-Peril MA, Moran TP, Immormino RM, Kesselring JR, Steinbach EC, Orgel KA, Staats HF, Burks AW, Mucha PJ, Ferris MT, Kulis MD. Fecal IgA, Antigen Absorption, and Gut Microbiome Composition Are Associated With Food Antigen Sensitization in Genetically Susceptible Mice. Front Immunol 2021; 11:599637. [PMID: 33542716 PMCID: PMC7850988 DOI: 10.3389/fimmu.2020.599637] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2020] [Accepted: 11/25/2020] [Indexed: 01/04/2023] Open
Abstract
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.
Collapse
|
38
|
Graham JB, Swarts JL, Leist SR, Schäfer A, Menachery VD, Gralinski LE, Jeng S, Miller DR, Mooney MA, McWeeney SK, Ferris MT, Pardo-Manuel de Villena F, Heise MT, Baric RS, Lund JM. Baseline T cell immune phenotypes predict virologic and disease control upon SARS-CoV infection in Collaborative Cross mice. PLoS Pathog 2021; 17:e1009287. [PMID: 33513210 PMCID: PMC7875398 DOI: 10.1371/journal.ppat.1009287] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2020] [Revised: 02/10/2021] [Accepted: 01/05/2021] [Indexed: 12/25/2022] Open
Abstract
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.
Collapse
|
39
|
Sigmon JS, Blanchard MW, Baric RS, Bell TA, Brennan J, Brockmann GA, Burks AW, Calabrese JM, Caron KM, Cheney RE, Ciavatta D, Conlon F, Darr DB, Faber J, Franklin C, Gershon TR, Gralinski L, Gu B, Gaines CH, Hagan RS, Heimsath EG, Heise MT, Hock P, Ideraabdullah F, Jennette JC, Kafri T, Kashfeen A, Kulis M, Kumar V, Linnertz C, Livraghi-Butrico A, Lloyd KCK, Lutz C, Lynch RM, Magnuson T, Matsushima GK, McMullan R, Miller DR, Mohlke KL, Moy SS, Murphy CEY, Najarian M, O'Brien L, Palmer AA, Philpot BD, Randell SH, Reinholdt L, Ren Y, Rockwood S, Rogala AR, Saraswatula A, Sassetti CM, Schisler JC, Schoenrock SA, Shaw GD, Shorter JR, Smith CM, St Pierre CL, Tarantino LM, Threadgill DW, Valdar W, Vilen BJ, Wardwell K, Whitmire JK, Williams L, Zylka MJ, Ferris MT, McMillan L, Manuel de Villena FP. Content and Performance of the MiniMUGA Genotyping Array: A New Tool To Improve Rigor and Reproducibility in Mouse Research. Genetics 2020; 216:905-930. [PMID: 33067325 PMCID: PMC7768238 DOI: 10.1534/genetics.120.303596] [Citation(s) in RCA: 42] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2020] [Accepted: 10/06/2020] [Indexed: 12/14/2022] Open
Abstract
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.
Collapse
|
40
|
Graham JB, Swarts JL, Leist SR, Schäfer A, Menachery VD, Gralinski LE, Jeng S, Miller DR, Mooney MA, McWeeney SK, Ferris MT, de Villena FPM, Heise MT, Baric RS, Lund JM. Baseline T cell immune phenotypes predict virologic and disease control upon SARS-CoV infection. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2020:2020.09.21.306837. [PMID: 32995791 PMCID: PMC7523117 DOI: 10.1101/2020.09.21.306837] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
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.
Collapse
|
41
|
Graham JB, Swarts JL, Menachery VD, Gralinski LE, Schäfer A, Plante KS, Morrison CR, Voss KM, Green R, Choonoo G, Jeng S, Miller DR, Mooney MA, McWeeney SK, Ferris MT, Pardo-Manuel de Villena F, Gale M, Heise MT, Baric RS, Lund JM. Immune Predictors of Mortality After Ribonucleic Acid Virus Infection. J Infect Dis 2020; 221:882-889. [PMID: 31621854 PMCID: PMC7107456 DOI: 10.1093/infdis/jiz531] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2019] [Accepted: 10/11/2019] [Indexed: 12/14/2022] Open
Abstract
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.
Collapse
|
42
|
Hampton BK, Jensen K, Whitmore AC, Morrison CR, Plante KS, Leist SR, Gralinski LE, Menachery VD, Schäfer A, de Villena FPM, Baric RS, Heise MT, Ferris MT. Genetic regulation of immune homeostatic lung leukocyte populations. THE JOURNAL OF IMMUNOLOGY 2020. [DOI: 10.4049/jimmunol.204.supp.235.4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Abstract
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.
Collapse
|
43
|
Noll KE, Whitmore AC, West A, McCarthy MK, Morrison CR, Plante KS, Hampton BK, Kollmus H, Pilzner C, Leist SR, Gralinski LE, Menachery VD, Schäfer A, Miller D, Shaw G, Mooney M, McWeeney S, Pardo-Manuel de Villena F, Schughart K, Morrison TE, Baric RS, Ferris MT, Heise MT. Complex Genetic Architecture Underlies Regulation of Influenza-A-Virus-Specific Antibody Responses in the Collaborative Cross. Cell Rep 2020; 31:107587. [PMID: 32348764 PMCID: PMC7195006 DOI: 10.1016/j.celrep.2020.107587] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2019] [Revised: 02/20/2020] [Accepted: 04/08/2020] [Indexed: 02/07/2023] Open
Abstract
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.
Collapse
|
44
|
Shorter JR, Najarian ML, Bell TA, Blanchard M, Ferris MT, Hock P, Kashfeen A, Kirchoff KE, Linnertz CL, Sigmon JS, Miller DR, McMillan L, Pardo-Manuel de Villena F. Whole Genome Sequencing and Progress Toward Full Inbreeding of the Mouse Collaborative Cross Population. G3 (BETHESDA, MD.) 2019; 9:1303-1311. [PMID: 30858237 PMCID: PMC6505143 DOI: 10.1534/g3.119.400039] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/28/2019] [Accepted: 03/08/2019] [Indexed: 12/20/2022]
Abstract
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.
Collapse
|
45
|
Menachery VD, McAnarney ET, Ferris MT. Age-dependent Susceptibility To Virus Infection Is Governed by Host Genetics. THE JOURNAL OF IMMUNOLOGY 2019. [DOI: 10.4049/jimmunol.202.supp.197.9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Abstract
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.
Collapse
|
46
|
Noll KE, Ferris MT, Heise MT. The Collaborative Cross: A Systems Genetics Resource for Studying Host-Pathogen Interactions. Cell Host Microbe 2019; 25:484-498. [PMID: 30974083 PMCID: PMC6494101 DOI: 10.1016/j.chom.2019.03.009] [Citation(s) in RCA: 45] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
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.
Collapse
|
47
|
Ferris MT, Hood DW. Host genetic regulation of immune-based and infectious diseases : Introduction to mammalian genome special issue: genetics of infectious disease. Mamm Genome 2019; 29:365-366. [PMID: 30171339 PMCID: PMC7088047 DOI: 10.1007/s00335-018-9779-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
48
|
Orgel K, Smeekens JM, Ye P, Fotsch L, Guo R, Miller DR, Pardo-Manuel de Villena F, Burks AW, Ferris MT, Kulis MD. Genetic diversity between mouse strains allows identification of the CC027/GeniUnc strain as an orally reactive model of peanut allergy. J Allergy Clin Immunol 2018; 143:1027-1037.e7. [PMID: 30342892 PMCID: PMC7252586 DOI: 10.1016/j.jaci.2018.10.009] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2017] [Revised: 09/17/2018] [Accepted: 10/01/2018] [Indexed: 02/08/2023]
Abstract
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 TH2 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.
Collapse
|
49
|
Graham JB, Swarts JL, Mooney M, Choonoo G, Jeng S, Miller DR, Ferris MT, McWeeney S, Lund JM. Extensive Homeostatic T Cell Phenotypic Variation within the Collaborative Cross. Cell Rep 2018; 21:2313-2325. [PMID: 29166619 PMCID: PMC5728448 DOI: 10.1016/j.celrep.2017.10.093] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2017] [Revised: 07/14/2017] [Accepted: 10/25/2017] [Indexed: 12/11/2022] Open
Abstract
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
Collapse
|
50
|
McMullan RC, Ferris MT, Bell TA, Menachery VD, Baric RS, Hua K, Pomp D, Smith‐Ryan AE, de Villena FP. CC002/Unc females are mouse models of exercise-induced paradoxical fat response. Physiol Rep 2018; 6:e13716. [PMID: 29924460 PMCID: PMC6009762 DOI: 10.14814/phy2.13716] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2018] [Accepted: 05/03/2018] [Indexed: 12/13/2022] Open
Abstract
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.
Collapse
|