Using the emerging Collaborative Cross to probe the immune system

Variations in the germinal center response revealed by genetically diverse mouse strains

The humoral response is complex and involves multiple cellular populations and signaling pathways. Bacterial and viral infections, as well as immunization regimens, can trigger this type of response, promoting the formation of microanatomical cellular structures called germinal centers (GCs). GCs formed in secondary lymphoid organs support the differentiation of high-affinity plasma cells and memory B cells. There is growing evidence that the quality of the humoral response is influenced by genetic variants. Using 12 genetically divergent mouse strains, we assessed the impact of genetics on GC cellular traits. At steady state, in the spleen, lymph nodes and Peyer's patches, we quantified GC B cells, plasma cells and follicular helper T cells. These traits were also quantified in the spleen of mice following immunization with a foreign antigen, namely, sheep red blood cells, in addition to the number and size of GCs. We observed both strain- and organ-specific variations in cell type abundance, as well as for GC number and size. Moreover, we find that some of these traits are highly heritable. Importantly, the results of this study inform on the impact of genetic diversity in shaping the GC response and identify the traits that are the most impacted by genetic background.

Analysis of hepatic lentiviral vector transduction; implications for preclinical studies and clinical gene therapy protocols

Lentiviral vector-transduced T-cells were approved by the FDA as gene therapy anti-cancer medications. Little is known about the host genetic variation effects on the safety and efficacy of the lentiviral vector gene delivery system. To narrow this knowledge-gap, we characterized hepatic gene delivery by lentiviral vectors across the Collaborative Cross (CC) mouse genetic reference population. For 24 weeks, we periodically measured hepatic luciferase expression from lentiviral vectors in 41 CC mouse strains. Hepatic and splenic vector copy numbers were determined. We report that CC mouse strains showed highly diverse outcomes following lentiviral gene delivery. For the first time, moderate correlation between mouse strain-specific sleeping patterns and transduction efficiency was observed. We associated two quantitative trait loci (QTLs) with intra-strain variations in transduction phenotypes, which mechanistically relates to the phenomenon of metastable epialleles. An additional QTL was associated with the kinetics of hepatic transgene expression. Genes comprised in the above QTLs are potential targets to personalize gene therapy protocols. Importantly, we identified two mouse strains that open new directions in characterizing continuous viral vector silencing and HIV latency. Our findings suggest that wide-range patient-specific outcomes of viral vector-based gene therapy should be expected. Thus, novel escalating dose-based clinical protocols should be considered.

Quantitative trait loci mapping provides insights into the genetic regulation of dendritic cell numbers in mouse tissues

To explore the influence of genetics on homeostatic regulation of dendritic cell (DC) numbers, we present a screen of DCs and their progenitors in lymphoid and non-lymphoid tissues in Collaborative Cross (CC) and Diversity Outbred (DO) mice. We report 30 and 71 loci with logarithm of the odds (LOD) scores >8.18 and ranging from 6.67 to 8.19, respectively. The analysis reveals the highly polygenic and pleiotropic architecture of this complex trait, including many of the previously identified genetic regulators of DC development and maturation. Two SNPs in genes potentially underlying variation in DC homeostasis, a splice variant in Gramd4 (rs235532740) and a missense variant in Orai3 (rs216659754), are confirmed by gene editing using CRISPR-Cas9. Gramd4 is a central regulator of DC homeostasis that impacts the entire DC lineage, and Orai3 regulates cDC2 numbers in tissues. Overall, the data reveal a large number of candidate genes regulating DC homeostasis in vivo.

Identifying genetic susceptibility to Aspergillus fumigatus infection using collaborative cross mice and RNA-Seq approach

BACKGROUND

Aspergillus fumigatus (Af) is one of the most ubiquitous fungi and its infection potency is suggested to be strongly controlled by the host genetic background. The aim of this study was to search for candidate genes associated with host susceptibility to Aspergillus fumigatus (Af) using an RNAseq approach in CC lines and hepatic gene expression.

METHODS

We studied 31 male mice from 25 CC lines at 8 weeks old; the mice were infected with Af. Liver tissues were extracted from these mice 5 days post-infection, and next-generation RNA-sequencing (RNAseq) was performed. The GENE-E analysis platform was used to generate a clustered heat map matrix.

RESULTS

Significant variation in body weight changes between CC lines was observed. Hepatic gene expression revealed 12 top prioritized candidate genes differentially expressed in resistant versus susceptible mice based on body weight changes. Interestingly, three candidate genes are located within genomic intervals of the previously mapped quantitative trait loci (QTL), including Gm16270 and Stox1 on chromosome 10 and Gm11033 on chromosome 8.

CONCLUSIONS

Our findings emphasize the CC mouse model's power in fine mapping the genetic components underlying susceptibility towards Af. As a next step, eQTL analysis will be performed for our RNA-Seq data. Suggested candidate genes from our study will be further assessed with a human cohort with aspergillosis.

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

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

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

Inflammation in response to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection drives severity of coronavirus disease 2019 (COVID-19) and is influenced by host genetics. To understand mechanisms of inflammation, animal models that reflect genetic diversity and clinical outcomes observed in humans are needed. We report a mouse panel comprising the genetically diverse Collaborative Cross (CC) founder strains crossed to human ACE2 transgenic mice (K18-hACE2) that confers susceptibility to SARS-CoV-2. Infection of CC x K18-hACE2 resulted in a spectrum of survival, viral replication kinetics, and immune profiles. Importantly, in contrast to the K18-hACE2 model, early type I interferon (IFN-I) and regulated proinflammatory responses were required for control of SARS-CoV-2 replication in PWK x K18-hACE2 mice that were highly resistant to disease. Thus, virus dynamics and inflammation observed in COVID-19 can be modeled in diverse mouse strains that provide a genetically tractable platform for understanding anti-coronavirus immunity.

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

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

Using the Collaborative Cross and Diversity Outbred Mice in Immunology

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

Elucidating Mechanisms of Tolerance to Salmonella Typhimurium across Long-Term Infections Using the Collaborative Cross

Understanding the molecular mechanisms underlying resistance and tolerance to pathogen infection may present the opportunity to develop novel interventions. Resistance is the absence of clinical disease with a low pathogen burden, while tolerance is minimal clinical disease with a high pathogen burden. Salmonella is a worldwide health concern. We studied 18 strains of collaborative cross mice that survive acute Salmonella Typhimurium (STm) infections. We infected these strains orally and monitored them for 3 weeks. Five strains cleared STm (resistant), six strains maintained a bacterial load and survived (tolerant), while seven strains survived >7 days but succumbed to infection within the study period and were called “delayed susceptible.” Tolerant strains were colonized in the Peyer’s patches, mesenteric lymph node, spleen, and liver, while resistant strains had significantly reduced bacterial colonization. Tolerant strains had lower preinfection core body temperatures and had disrupted circadian patterns of body temperature postinfection sooner than other strains. Tolerant strains had higher circulating total white blood cells than resistant strains, driven by increased numbers of neutrophils. Tolerant strains had more severe tissue damage and higher circulating levels of monocyte chemoattractant protein 1 (MCP-1) and interferon gamma (IFN-γ), but lower levels of epithelial neutrophil-activating protein 78 (ENA-78) than resistant strains. Quantitative trait locus (QTL) analysis revealed one significant association and six suggestive associations. Gene expression analysis identified 22 genes that are differentially regulated in tolerant versus resistant animals that overlapped these QTLs. Fibrinogen genes (Fga, Fgb, and Fgg) were found across the QTL, RNA, and top canonical pathways, making them the best candidate genes for differentiating tolerance and resistance.

Improving mouse models for the study of Alzheimer's disease

Alzheimer's disease (AD) is a complex neurodegenerative disease whose risk is influenced by genetic and environmental factors. Although a number of pathological hallmarks have been extensively studied over the last several decades, a complete picture of disease initiation and progression remains unclear. We now understand that numerous cell types and systems are involved in AD pathogenesis, and that this cellular profile may present differently for each individual, making the creation of relevant mouse models challenging. However, with increasingly diverse data made available by genome-wide association studies, we can identify and examine new genes and pathways involved in genetic risk for AD, many of which involve vascular health and inflammation. When developing mouse models, it is critical to assess (1) an aging timeline that represents onset and progression in humans, (2) genetic variants and context, (3) environmental factors present in human populations that result in both neuropathological and functional changes-themes that we address in this chapter.

A diversity outbred F1 mouse model identifies host-intrinsic genetic regulators of response to immune checkpoint inhibitors

Immune checkpoint inhibitors (ICI) have improved outcomes for a variety of malignancies; however, many patients fail to benefit. While tumor-intrinsic mechanisms are likely involved in therapy resistance, it is unclear to what extent host genetic background influences response. To investigate this, we utilized the Diversity Outbred (DO) and Collaborative Cross (CC) mouse models. DO mice are an outbred stock generated by crossbreeding eight inbred founder strains, and CC mice are recombinant inbred mice generated from the same eight founders. We generated 207 DOB6F1 mice representing 48 DO dams and demonstrated that these mice reliably accept the C57BL/6-syngeneic B16F0 tumor and that host genetic background influences response to ICI. Genetic linkage analysis from 142 mice identified multiple regions including one within chromosome 13 that associated with therapeutic response. We utilized 6 CC strains bearing the positive (NZO) or negative (C57BL/6) driver genotype in this locus. We found that 2/3 of predicted responder CCB6F1 crosses show reproducible ICI response. The chromosome 13 locus contains the murine prolactin family, which is a known immunomodulating cytokine associated with various autoimmune disorders. To directly test whether prolactin influences ICI response rates, we implanted inbred C57BL/6 mice with subcutaneous slow-release prolactin pellets to induce mild hyperprolactinemia. Prolactin augmented ICI response against B16F0, with increased CD8 infiltration and 5/8 mice exhibiting slowed tumor growth relative to controls. This study highlights the role of host genetics in ICI response and supports the use of F1 crosses in the DO and CC mouse populations as powerful cancer immunotherapy models.

Genetic background influences survival of infections with Salmonella enterica serovar Typhimurium in the Collaborative Cross

Salmonella infections typically cause self-limiting gastroenteritis, but in some individuals these bacteria can spread systemically and cause disseminated disease. Salmonella Typhimurium (STm), which causes severe systemic disease in most inbred mice, has been used as a model for disseminated disease. To screen for new infection phenotypes across a range of host genetics, we orally infected 32 Collaborative Cross (CC) mouse strains with STm and monitored their disease progression for seven days by telemetry. Our data revealed a broad range of phenotypes across CC strains in many parameters including survival, bacterial colonization, tissue damage, complete blood counts (CBC), and serum cytokines. Eighteen CC strains survived to day 7, while fourteen susceptible strains succumbed to infection before day 7. Several CC strains had sex differences in survival and colonization. Surviving strains had lower pre-infection baseline temperatures and were less active during their daily active period. Core body temperature disruptions were detected earlier after STm infection than activity disruptions, making temperature a better detector of illness. All CC strains had STm in spleen and liver, but susceptible strains were more highly colonized. Tissue damage was weakly negatively correlated to survival. We identified loci associated with survival on Chromosomes (Chr) 1, 2, 4, 7. Polymorphisms in Ncf2 and Slc11a1, known to reduce survival in mice after STm infections, are located in the Chr 1 interval, and the Chr 7 association overlaps with a previously identified QTL peak called Ses2. We identified two new genetic regions on Chr 2 and 4 associated with susceptibility to STm infection. Our data reveal the diversity of responses to STm infection across a range of host genetics and identified new candidate regions for survival of STm infection.

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

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

ONE SENTENCE SUMMARY

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

Variation in Host Resistance to Blastomyces dermatitidis: Potential Use of Genetic Reference Panels and Advances in Immunophenotyping of Diverse Mouse Strains

Host genetic determinants that underpin variation in susceptibility to systemic fungal infection are poorly understood. Genes responsible for complex traits can be identified by correlating variation in phenotype with allele in founder strains of wild mice with known genetic variation, assembled in genetic reference panels. In this work, we describe wide natural variation in both primary and acquired resistance to experimental pulmonary blastomycosis in eight founder strains, including 129, A/J, BL/6, CAST, NOD, NZO, PWK, and WSB of the Collaborative Cross collection, and the inbred DBA strain. These differences in susceptibility across strains were accompanied by sharp differences in the accumulation and function of immune cells in the lungs. Immune perturbations were mapped by identifying reagents that phenotypically mark immune cell populations in the distinct strains of mice. In particular, we uncovered marked differences between BL/6 and DBA/2 mouse strains in the development of acquired resistance. Our findings highlight the potential value in using genetic reference panels of mice, and particularly the BXD (recombinant inbred strains of mice from a cross of C57BL/6J and DBA/2J mice) collection harboring a cross between resistant BL/6 and susceptible DBA/2 mice, for unveiling genes linked with host resistance to fungal infection. IMPORTANCE Host genetic variation significantly impacts vulnerability to infectious diseases. While host variation in susceptibility to fungal infection with dimorphic fungi has long been recognized, genes that underpin this variation are poorly understood. We used a collection of seven mouse strains that represent nearly 90% of the genetic variation in mice to identify genetic variability among the strains in resistance to pulmonary infection with the dimorphic fungus Blastomyces dermatitidis. We analyzed differences between the strains in innate resistance by infecting naive mice and in acquired resistance by infecting vaccinated mice. We identified extreme variations in both innate and acquired resistance among the strains. In particular, we found sharp differences between C57BL/6 and DBA/2 strains in the ability to acquire vaccine-induced resistance. We also identified commercial reagents that allowed the phenotyping of immune cells from this strain collection of mice. Because there are additional mice harboring a genetic cross of the C57BL/6 and DBA/2 strains (BXD collection), such mice will permit future investigations to identify the genes that underlie differences in the ability to acquire resistance to infection.

Divergent Genetic Regulation of Nitric Oxide Production between C57BL/6J and Wild-Derived PWD/PhJ Mice Controls Postactivation Mitochondrial Metabolism, Cell Survival, and Bacterial Resistance in Dendritic Cells

Dendritic cell (DC) activation is characterized by sustained commitment to glycolysis that is a requirement for survival in DC subsets that express inducible NO synthase (Nos2) due to NO-mediated inhibition of mitochondrial respiration. This phenomenon primarily has been studied in DCs from the classic laboratory inbred mouse strain C57BL/6J (B6) mice, where DCs experience a loss of mitochondrial function due to NO accumulation. To assess the conservation of NO-driven metabolic regulation in DCs, we compared B6 mice to the wild-derived genetically divergent PWD/PhJ (PWD) strain. We show preserved mitochondrial respiration and enhanced postactivation survival due to attenuated NO production in LPS-stimulated PWD DCs phenocopying human monocyte-derived DCs. To genetically map this phenotype, we used a congenic mouse strain (B6.PWD-Chr11.2) that carries a PWD-derived portion of chromosome 11, including Nos2, on a B6 background. B6.PWD-Chr11.2 DCs show preserved mitochondrial function and produce lower NO levels than B6 DCs. We demonstrate that activated B6.PWD-Chr11.2 DCs maintain mitochondrial respiration and TCA cycle carbon flux, compared with B6 DCs. However, reduced NO production by the PWD Nos2 allele results in impaired cellular control of Listeria monocytogenes replication. These studies establish a natural genetic model for restrained endogenous NO production to investigate the contribution of NO in regulating the interplay between DC metabolism and immune function. These findings suggest that reported differences between human and murine DCs may be an artifact of the limited genetic diversity of the mouse models used, underscoring the need for mouse genetic diversity in immunology research.

Genetic and immunological contributors to virus-induced paralysis

Infection by a single virus can evoke diverse immune responses, resulting in different neurological outcomes, depending on the host's genetic background. To study heterogenous viral response, we use Theiler's Murine Encephalomyelitis Virus (TMEV) to model virally induced neurological phenotypes and immune responses in Collaborative Cross (CC) mice. The CC resource consists of genetically distinct and reproducible mouse lines, thus providing a population model with genetic heterogeneity similar to humans. We examined different CC strains for the effect of chronic stage TMEV-induced immune responses on neurological outcomes throughout 90 days post infection (dpi), with a particular focus on limb paralysis, by measuring serum levels of 23 different cytokines and chemokines. Each CC strain demonstrated a unique set of immune responses, regardless of presence or absence of TMEV RNA. Using stepwise regression, significant associations were identified between IL-1α, RANTES, and paralysis frequency scores. To better understand these interactions, we evaluated multiple aspects of the different CC genetic backgrounds, including haplotypes of genomic regions previously linked with TMEV pathogenesis and viral clearance or persistence, individual cytokine levels, and TMEV-relevant gene expression. These results demonstrate how loci previously associated with TMEV outcomes provide incomplete information regarding TMEV-induced paralysis in the CC strains. Overall, these findings provide insight into the complex roles of immune response in the pathogenesis of virus-associated neurological diseases influenced by host genetic background.

Utilizing preclinical models of genetic diversity to improve translation of phytochemical activities from rodents to humans and inform personalized nutrition

Mouse models are an essential tool in different areas of research, including nutrition and phytochemical research. Traditional inbred mouse models have allowed the discovery of therapeutical targets and mechanisms of action and expanded our knowledge of health and disease. However, these models lack the genetic variability typically found in human populations, which hinders the translatability of the results found in mice to humans. The development of genetically diverse mouse models, such as the collaborative cross (CC) or the diversity outbred (DO) models, has been a useful tool to overcome this obstacle in many fields, such as cancer, immunology and toxicology. However, these tools have not yet been widely adopted in the field of phytochemical research. As demonstrated in other disciplines, use of CC and DO models has the potential to provide invaluable insights for translation of phytochemicals from rodents to humans, which are desperately needed given the challenges and numerous failed clinical trials in this field. These models may prove informative for personalized use of phytochemicals in humans, including: predicting interindividual variability in phytochemical bioavailability and efficacy, identifying genetic loci or genes governing response to phytochemicals, identifying phytochemical mechanisms of action and therapeutic targets, and understanding the impact of genetic variability on individual response to phytochemicals. Such insights would prove invaluable for personalized implementation of phytochemicals in humans. This review will focus on the current work performed with genetically diverse mouse populations, and the research opportunities and advantages that these models can offer to phytochemical research.

Conservation and divergence in gene regulation between mouse and human immune cells deserves equal emphasis

Model organisms such as mice are important for basic research and serve as valuable tools in preclinical translational studies. A challenge with translating findings from mice to humans is identifying and separating evolutionarily conserved mechanisms in the immune system from those diverging between species. A significant emphasis has been placed on defining conserved gene regulation principles, with divergent mechanisms often overlooked. We put forward the perspective that both conserved and divergent mechanisms that regulate gene expression programs are of equal importance. With recent advances and availability of datasets, immunologists should take a closer look at the role for genetic diversity in altering gene expression programs between mouse and human immune cells.

Human-relevant mechanisms and risk factors for TAK-875-Induced liver injury identified via a gene pathway-based approach in Collaborative Cross mice

Development of TAK-875 was discontinued when a small number of serious drug-induced liver injury (DILI) cases were observed in Phase 3 clinical trials. Subsequent studies have identified hepatocellular oxidative stress, mitochondrial dysfunction, altered bile acid homeostasis, and immune response as mechanisms of TAK-875 DILI and the contribution of genetic risk factors in oxidative response and mitochondrial pathways to the toxicity susceptibility observed in patients. We tested the hypothesis that a novel preclinical approach based on gene pathway analysis in the livers of Collaborative Cross mice could be used to identify human-relevant mechanisms of toxicity and genetic risk factors at the level of the hepatocyte as reported in a human genome-wide association study. Eight (8) male mice (4 matched pairs) from each of 45 Collaborative Cross lines were treated with a single oral (gavage) dose of either vehicle or 600 mg/kg TAK-875. As expected, liver injury was not detected histologically and few changes in plasma biomarkers of hepatotoxicity were observed. However, gene expression profiling in the liver identified hundreds of transcripts responsive to TAK-875 treatment across all strains reflecting alterations in immune response and bile acid homeostasis and the interaction of treatment and strain reflecting oxidative stress and mitochondrial dysfunction. Fold-change expression values were then used to develop pathway-based phenotypes for genetic mapping which identified candidate risk factor genes for TAK-875 toxicity susceptibility at the level of the hepatocyte. Taken together, these findings support our hypothesis that a gene pathway-based approach using Collaborative Cross mice could inform sensitive strains, human-relevant mechanisms of toxicity, and genetic risk factors for TAK-875 DILI. This novel preclinical approach may be helpful in understanding, predicting, and ultimately preventing clinical DILI for other drugs.

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

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

Diverse CD8 T Cell Responses to Viral Infection Revealed by the Collaborative Cross

Enhanced host protection against re-infection requires generation of memory T cells of sufficient quantity and functional quality. Unlike well-studied inbred mice, T cell responses of diverse size and quality are generated following infection of humans and outbred mice. Thus, additional models are needed that accurately reflect variation in immune outcomes in genetically diverse populations and to uncover underlying genetic causes. The Collaborative Cross (CC), a large recombinant inbred panel of mice, is an ideal model in this pursuit for the high degree of genetic variation present, because it allows for assessment of genetic factors underlying unique phenotypes. Here, we advance the utility of the CC as a tool to analyze the immune response to viral infection. We describe variability in resting immune cell composition and adaptive immune responses generated among CC strains following systemic virus infection and reveal quantitative trait loci responsible for generation of CD62L+ memory CD8 T cells.

Martin et al. advance the use of the Collaborative Cross (CC) for studying adaptive immune responses. They demonstrate that the CC better models variation in T cell responses seen in outbred mice and humans and that it can uncover genes linked to generation of qualitatively distinct memory cells following infection.

Synergizing Mouse and Human Studies to Understand the Heterogeneity of Obesity

Obesity is routinely considered as a single disease state, which drives a "one-size-fits-all" approach to treatment. We recently convened the first annual University of North Carolina Interdisciplinary Nutrition Sciences Symposium to discuss the heterogeneity of obesity and the need for translational science to advance understanding of this heterogeneity. The symposium aimed to advance scientific rigor in translational studies from animal to human models with the goal of identifying underlying mechanisms and treatments. In this review, we discuss fundamental gaps in knowledge of the heterogeneity of obesity ranging from cellular to population perspectives. We also advocate approaches to overcoming limitations in the field. Examples include the use of contemporary mouse genetic reference population models such as the Collaborative Cross and Diversity Outbred mice that effectively model human genetic diversity and the use of translational models that integrate -omics and computational approaches from pre-clinical to clinical models of obesity. Finally, we suggest best scientific practices to ensure strong rigor that will allow investigators to delineate the sources of heterogeneity in the population with obesity. Collectively, we propose that it is critical to think of obesity as a heterogeneous disease with complex mechanisms and etiologies, requiring unique prevention and treatment strategies tailored to the individual.

Onchocerca volvulus bivalent subunit vaccine induces protective immunity in genetically diverse collaborative cross recombinant inbred intercross mice

This study tests the hypothesis that an Onchocerca volvulus vaccine, consisting of two recombinant antigens (Ov-103 and Ov-RAL-2) formulated with the combination-adjuvant Advax-2, can induce protective immunity in genetically diverse Collaborative Cross recombinant inbred intercross mice (CC-RIX). CC-RIX lines were immunized with the O. volvulus vaccine and challenged with third-stage larvae. Equal and significant reductions in parasite survival were observed in 7 of 8 CC-RIX lines. Innate protective immunity was seen in the single CC-RIX line that did not demonstrate protective adaptive immunity. Analysis of a wide array of immune factors showed that each line of mice have a unique set of immune responses to vaccination and challenge suggesting that the vaccine is polyfunctional, inducing different equally-protective sets of immune responses based on the genetic background of the immunized host. Vaccine efficacy in genetically diverse mice suggests that it will also be effective in genetically complex human populations.

FLT3 Ligand Is Dispensable for the Final Stage of Type 1 Conventional Dendritic Cell Differentiation

Conventional dendritic cells (cDCs) are comprised of two major subsets, type 1 cDC (cDC1) and type 2 cDC (cDC2). As each cDC subset differentially influences the nature of immune responses, we sought factors that would allow the manipulation of their relative abundance. Notably, cDC1 are less abundant than cDC2 in both lymphoid and nonlymphoid organs. We demonstrate that this bias is already apparent in bone marrow precommitted precursors. However, comparison of five common inbred strains revealed a disparity in precursor-product relationship, in which mice with fewer precursors to cDC1 had more cDC1. This disparity associated with contrasting variations in CD135 (FLT3) expression on cDC subsets. Hence, we characterized the response to FLT3 ligand during cDC1 and cDC2 lineage differentiation and find that although FLT3 ligand is required throughout cDC2 differentiation, it is surprisingly dispensable during late-stage cDC1 differentiation. Overall, we find that tight regulation of FLT3 ligand levels throughout cDC differentiation dictates the cDC1 to cDC2 ratio in lymphoid organs.

Genomic, microbial and environmental standardization in animal experimentation limiting immunological discovery

Background

The use of inbred mice housed under standardized environmental conditions has been critical in identifying immuno-pathological mechanisms in different infectious and inflammatory diseases as well as revealing new therapeutic targets for clinical trials. Unfortunately, only a small percentage of preclinical intervention studies using well-defined mouse models of disease have progressed to clinically-effective treatments in patients. The reasons for this lack of bench-to-bedside transition are not completely understood; however, emerging data suggest that genetic diversity and housing environment may greatly influence muring immunity and inflammation.

Results

Accumulating evidence suggests that certain immune responses and/or disease phenotypes observed in inbred mice may be quite different than those observed in their outbred counterparts. These differences have been thought to contribute to differing immune responses to foreign and/or auto-antigens in mice vs. humans. There is also a growing literature demonstrating that mice housed under specific pathogen free conditions possess an immature immune system that remarkably affects their ability to respond to pathogens and/or inflammation when compared with mice exposed to a more diverse spectrum of microorganisms. Furthermore, recent studies demonstrate that mice develop chronic cold stress when housed at standard animal care facility temperatures (i.e. 22–24 °C). These temperatures have been shown alter immune responses to foreign and auto-antigens when compared with mice housed at their thermo-neutral body temperature of 30–32 °C.

Conclusions

Exposure of genetically diverse mice to a spectrum of environmentally-relevant microorganisms at housing temperatures that approximate their thermo-neutral zone may improve the chances of identifying new and more potent therapeutics to treat infectious and inflammatory diseases.

Identification of Candidate Risk Factor Genes for Human Idelalisib Toxicity Using a Collaborative Cross Approach

Idelalisib is a phosphatidylinositol 3-kinase inhibitor highly selective for the delta isoform that has shown good efficacy in treating chronic lymphocytic leukemia and follicular lymphoma. In clinical trials, however, idelalisib was associated with rare, but potentially serious liver and lung toxicities. In this study, we used the Collaborative Cross (CC) mouse population to identify genetic factors associated with the drug response that may inform risk management strategies for idelalisib in humans. Eight male mice (4 matched pairs) from 50 CC lines were treated once daily for 14 days by oral gavage with either vehicle or idelalisib at a dose selected to achieve clinically relevant peak plasma concentrations (150 mg/kg/day). The drug was well tolerated across all CC lines, and there were no observations of overt liver injury. Differences across CC lines were seen in drug concentration in plasma samples collected at the approximate Tmax on study Days 1, 7, and 14. There were also small but statistically significant treatment-induced alterations in plasma total bile acids and microRNA-122, and these may indicate early hepatocellular stress required for immune-mediated hepatotoxicity in humans. Idelalisib treatment further induced significant elevations in the total cell count of terminal bronchoalveolar lavage fluid, which may be analogous to pneumonitis observed in the clinic. Genetic mapping identified loci associated with interim plasma idelalisib concentration and the other 3 treatment-related endpoints. Thirteen priority candidate quantitative trait genes identified in CC mice may now guide interrogation of risk factors for adverse drug responses associated with idelalisib in humans.

The Diversity Outbred Mouse Population Is an Improved Animal Model of Vaccination against Tuberculosis That Reflects Heterogeneity of Protection

We vaccinated the Diversity Outbred (DO) population of mice with BCG, the only vaccine currently used to protect against tuberculosis, and then challenged them with M. tuberculosis by aerosol. We found that the BCG-vaccinated DO mouse population exhibited a wide range of outcomes, in which outcomes in individual mice ranged from minimal respiratory or systemic disease to fulminant disease and death. The breadth of these outcomes appears similar to the range seen in people, indicating that DO mice may serve as an improved small-animal model to study tuberculosis infection and immunity. Moreover, sophisticated tools are available for the use of these mice to map genes contributing to control of vaccination. Thus, the present studies provided an important new tool in the fight against tuberculosis.

Many studies of Mycobacterium tuberculosis infection and immunity have used mouse models. However, outcomes of vaccination and challenge with M. tuberculosis in inbred mouse strains do not reflect the full range of outcomes seen in people. Previous studies indicated that the novel Diversity Outbred (DO) mouse population exhibited a spectrum of outcomes after primary aerosol infection with M. tuberculosis. Here, we demonstrate the value of this novel mouse population for studies of vaccination against M. tuberculosis aerosol challenge. Using the only currently licensed tuberculosis vaccine, we found that the DO population readily controlled systemic Mycobacterium bovis BCG bacterial burdens and that BCG vaccination significantly improved survival across the DO population upon challenge with M. tuberculosis. Many individual DO mice that were vaccinated with BCG and then challenged with M. tuberculosis exhibited low bacterial burdens, low or even no systemic dissemination, little weight loss, and only minor lung pathology. In contrast, some BCG-vaccinated DO mice progressed quickly to fulminant disease upon M. tuberculosis challenge. Across the population, most of these disease parameters were at most modestly correlated with each other and were often discordant. This result suggests the need for a multiparameter metric to better characterize “disease” and “protection,” with closer similarity to the complex case definitions used in people. Taken together, these results demonstrate that DO mice provide a novel small-animal model of vaccination against tuberculosis that better reflects the wide spectrum of outcomes seen in people.

IMPORTANCE We vaccinated the Diversity Outbred (DO) population of mice with BCG, the only vaccine currently used to protect against tuberculosis, and then challenged them with M. tuberculosis by aerosol. We found that the BCG-vaccinated DO mouse population exhibited a wide range of outcomes, in which outcomes in individual mice ranged from minimal respiratory or systemic disease to fulminant disease and death. The breadth of these outcomes appears similar to the range seen in people, indicating that DO mice may serve as an improved small-animal model to study tuberculosis infection and immunity. Moreover, sophisticated tools are available for the use of these mice to map genes contributing to control of vaccination. Thus, the present studies provided an important new tool in the fight against tuberculosis.

A comprehensive and comparative phenotypic analysis of the collaborative founder strains identifies new and known phenotypes

The collaborative cross (CC) is a large panel of mouse-inbred lines derived from eight founder strains (NOD/ShiLtJ, NZO/HILtJ, A/J, C57BL/6J, 129S1/SvImJ, CAST/EiJ, PWK/PhJ, and WSB/EiJ). Here, we performed a comprehensive and comparative phenotyping screening to identify phenotypic differences and similarities between the eight founder strains. In total, more than 300 parameters including allergy, behavior, cardiovascular, clinical blood chemistry, dysmorphology, bone and cartilage, energy metabolism, eye and vision, immunology, lung function, neurology, nociception, and pathology were analyzed; in most traits from sixteen females and sixteen males. We identified over 270 parameters that were significantly different between strains. This study highlights the value of the founder and CC strains for phenotype-genotype associations of many genetic traits that are highly relevant to human diseases. All data described here are publicly available from the mouse phenome database for analyses and downloads.

Immune Relevant and Immune Deficient Mice: Options and Opportunities in Translational Research

In 1989 ILAR published a list and description of immunodeficient rodents used in research. Since then, advances in understanding of molecular mechanisms; recognition of genetic, epigenetic microbial, and other influences on immunity; and capabilities in manipulating genomes and microbiomes have increased options and opportunities for selecting mice and designing studies to answer important mechanistic and therapeutic questions. Despite numerous scientific breakthroughs that have benefitted from research in mice, there is debate about the relevance and predictive or translational value of research in mice. Reproducibility of results obtained from mice and other research models also is a well-publicized concern. This review summarizes resources to inform the selection and use of immune relevant mouse strains and stocks, aiming to improve the utility, validity, and reproducibility of research in mice. Immune sufficient genetic variations, immune relevant spontaneous mutations, immunodeficient and autoimmune phenotypes, and selected induced conditions are emphasized.

Determinants of QTL Mapping Power in the Realized Collaborative Cross

The Collaborative Cross (CC) is a mouse genetic reference population whose range of applications includes quantitative trait loci (QTL) mapping. The design of a CC QTL mapping study involves multiple decisions, including which and how many strains to use, and how many replicates per strain to phenotype, all viewed within the context of hypothesized QTL architecture. Until now, these decisions have been informed largely by early power analyses that were based on simulated, hypothetical CC genomes. Now that more than 50 CC strains are available and more than 70 CC genomes have been observed, it is possible to characterize power based on realized CC genomes. We report power analyses from extensive simulations and examine several key considerations: 1) the number of strains and biological replicates, 2) the QTL effect size, 3) the presence of population structure, and 4) the distribution of functionally distinct alleles among the founder strains at the QTL. We also provide general power estimates to aide in the design of future experiments. All analyses were conducted with our R package, SPARCC (Simulated Power Analysis in the Realized Collaborative Cross), developed for performing either large scale power analyses or those tailored to particular CC experiments.

A Diallel of the Mouse Collaborative Cross Founders Reveals Strong Strain-Specific Maternal Effects on Litter Size

Reproductive success in the eight founder strains of the Collaborative Cross (CC) was measured using a diallel-mating scheme. Over a 48-month period we generated 4,448 litters, and provided 24,782 weaned pups for use in 16 different published experiments. We identified factors that affect the average litter size in a cross by estimating the overall contribution of parent-of-origin, heterosis, inbred, and epistatic effects using a Bayesian zero-truncated overdispersed Poisson mixed model. The phenotypic variance of litter size has a substantial contribution (82%) from unexplained and environmental sources, but no detectable effect of seasonality. Most of the explained variance was due to additive effects (9.2%) and parental sex (maternal vs. paternal strain; 5.8%), with epistasis accounting for 3.4%. Within the parental effects, the effect of the dam's strain explained more than the sire's strain (13.2% vs. 1.8%), and the dam's strain effects account for 74.2% of total variation explained. Dams from strains C57BL/6J and NOD/ShiLtJ increased the expected litter size by a mean of 1.66 and 1.79 pups, whereas dams from strains WSB/EiJ, PWK/PhJ, and CAST/EiJ reduced expected litter size by a mean of 1.51, 0.81, and 0.90 pups. Finally, there was no strong evidence for strain-specific effects on sex ratio distortion. Overall, these results demonstrate that strains vary substantially in their reproductive ability depending on their genetic background, and that litter size is largely determined by dam's strain rather than sire's strain effects, as expected. This analysis adds to our understanding of factors that influence litter size in mammals, and also helps to explain breeding successes and failures in the extinct lines and surviving CC strains.

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

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

Advanced Genetic Approaches in Discovery and Characterization of Genes Involved With Osteoporosis in Mouse and Human

Osteoporosis is a complex condition with contributions from, and interactions between, multiple genetic loci and environmental factors. This review summarizes key advances in the application of genetic approaches for the identification of osteoporosis susceptibility genes. Genome-wide linkage analysis (GWLA) is the classical approach for identification of genes that cause monogenic diseases; however, it has shown limited success for complex diseases like osteoporosis. In contrast, genome-wide association studies (GWAS) have successfully identified over 200 osteoporosis susceptibility loci with genome-wide significance, and have provided most of the candidate genes identified to date. Phenome-wide association studies (PheWAS) apply a phenotype-to-genotype approach which can be used to complement GWAS. PheWAS is capable of characterizing the association between osteoporosis and uncommon and rare genetic variants. Another alternative approach, whole genome sequencing (WGS), will enable the discovery of uncommon and rare genetic variants in osteoporosis. Meta-analysis with increasing statistical power can offer greater confidence in gene searching through the analysis of combined results across genetic studies. Recently, new approaches to gene discovery include animal phenotype based models such as the Collaborative Cross and ENU mutagenesis. Site-directed mutagenesis and genome editing tools such as CRISPR/Cas9, TALENs and ZNFs have been used in functional analysis of candidate genes in vitro and in vivo. These resources are revolutionizing the identification of osteoporosis susceptibility genes through the use of genetically defined inbred mouse libraries, which are screened for bone phenotypes that are then correlated with known genetic variation. Identification of osteoporosis-related susceptibility genes by genetic approaches enables further characterization of gene function in animal models, with the ultimate aim being the identification of novel therapeutic targets for osteoporosis.

Variable outcomes of human heart attack recapitulated in genetically diverse mice

Clinical variation in patient responses to myocardial infarction (MI) has been difficult to model in laboratory animals. To assess the genetic basis of variation in outcomes after heart attack, we characterized responses to acute MI in the Collaborative Cross (CC), a multi-parental panel of genetically diverse mouse strains. Striking differences in post-MI functional, morphological, and myocardial scar features were detected across 32 CC founder and recombinant inbred strains. Transcriptomic analyses revealed a plausible link between increased intrinsic cardiac oxidative phosphorylation levels and MI-induced heart failure. The emergence of significant quantitative trait loci for several post-MI traits indicates that utilizing CC strains is a valid approach for gene network discovery in cardiovascular disease, enabling more accurate clinical risk assessment and prediction.

Defining Memory CD8 T Cell

CD8 T cells comprising the memory pool display considerable heterogeneity, with individual cells differing in phenotype and function. This review will focus on our current understanding of heterogeneity within the antigen-specific memory CD8 T cell compartment and classifications of memory CD8 T cell subsets with defined and discrete functionalities. Recent data suggest that phenotype and/or function of numerically stable circulatory memory CD8 T cells are defined by the age of memory CD8 T cell (or time after initial antigen-encounter). In addition, history of antigen stimulations has a profound effect on memory CD8 T cell populations, suggesting that repeated infections (or vaccination) have the capacity to further shape the memory CD8 T cell pool. Finally, genetic background of hosts and history of exposure to diverse microorganisms likely contribute to the observed heterogeneity in the memory CD8 T cell compartment. Extending our tool box and exploring alternative mouse models (i.e., "dirty" and/or outbred mice) to encompass and better model diversity observed in humans will remain an important goal for the near future that will likely shed new light into the mechanisms that govern biology of memory CD8 T cells.

Giving the Genes a Shuffle: Using Natural Variation to Understand Host Genetic Contributions to Viral Infections

The laboratory mouse has proved an invaluable model to identify host factors that regulate the progression and outcome of virus-induced disease. The paradigm is to use single-gene knockouts in inbred mouse strains or genetic mapping studies using biparental mouse populations. However, genetic variation among these mouse strains is limited compared with the diversity seen in human populations. To address this disconnect, a multiparental mouse population has been developed to specifically dissect the multigenetic regulation of complex disease traits. The Collaborative Cross (CC) population of recombinant inbred mouse strains is a well-suited systems-genetics tool to identify susceptibility alleles that control viral and microbial infection outcomes and immune responses and to test the promise of personalized medicine.

Mapping novel genetic loci associated with female liver weight variations using Collaborative Cross mice

BACKGROUND

Liver weight is a complex trait, controlled by polygenic factors and differs within populations. Dissecting the genetic architecture underlying these variations will facilitate the search for key role candidate genes involved directly in the hepatomegaly process and indirectly involved in related diseases etiology.

METHODS

Liver weight of 506 mice generated from 39 different Collaborative Cross (CC) lines with both sexes at age 20 weeks old was determined using an electronic balance. Genomic DNA of the CC lines was genotyped with high-density single nucleotide polymorphic markers.

RESULTS

Statistical analysis revealed a significant (P < 0.05) variation of liver weight between the CC lines, with broad sense heritability (H 2) of 0.32 and genetic coefficient of variation (CVG) of 0.28. Subsequently, quantitative trait locus (QTL) mapping was performed, and results showed a significant QTL only for females on chromosome 8 at genomic interval 88.61-93.38 Mb (4.77 Mb). Three suggestive QTL were mapped at chromosomes 4, 12 and 13. The four QTL were designated as LWL1-LWL4 referring to liver weight loci 1-4 on chromosomes 8, 4, 12 and 13, respectively.

CONCLUSION

To our knowledge, this report presents, for the first time, the utilization of the CC for mapping QTL associated with baseline liver weight in mice. Our findings demonstrate that liver weight is a complex trait controlled by multiple genetic factors that differ significantly between sexes.

The Collaborative Cross mouse model for dissecting genetic susceptibility to infectious diseases

Infectious diseases, also known as communicable diseases, refer to a full range of maladies caused by pathogen invasion to the host body. Host response towards an infectious pathogen varies between individuals, and can be defined by responses from asymptomatic to lethal. Host response to infectious pathogens is considered as a complex trait controlled by gene-gene (host-pathogen) and gene-environment interactions, leading to the extensive phenotypic variations between individuals. With the advancement of the human genome mapping approaches and tools, various genome-wide association studies (GWAS) were performed, aimed at mapping the genetic basis underlying host susceptibility towards infectious pathogens. In parallel, immense efforts were invested in enhancing the genetic mapping resolution and gene-cloning efficacy, using advanced mouse models including advanced intercross lines; outbred populations; consomic, congenic; and recombinant inbred lines. Notwithstanding the evident advances achieved using these mouse models, the genetic diversity was low and quantitative trait loci (QTL) mapping resolution was inadequate. Consequently, the Collaborative Cross (CC) mouse model was established by full-reciprocal mating of eight divergent founder strains of mice (A/J, C57BL/6J, 129S1/SvImJ, NOD/LtJ, NZO/HiLtJ, CAST/Ei, PWK/PhJ, and WSB/EiJ) generating a next-generation mouse genetic reference population (CC lines). Presently, the CC mouse model population comprises a set of about 200 recombinant inbred CC lines exhibiting a unique high genetic diversity and which are accessible for multidisciplinary studies. The CC mouse model efficacy was validated by various studies in our lab and others, accomplishing high-resolution (< 1 MB) QTL genomic mapping for a variety of complex traits, using about 50 CC lines (3-4 mice per line). Herein, we present a number of studies demonstrating the power of the CC mouse model, which has been utilized in our lab for mapping the genetic basis of host susceptibility to various infectious pathogens. These include Aspergillus fumigatus, Klebsiella pneumoniae, Porphyromonas gingivalis and Fusobacterium nucleatum (causing oral mixed infection), Pseudomonas aeruginosa, and the bacterial toxins Lipopolysaccharide and Lipoteichoic acid.

Extensive Homeostatic T Cell Phenotypic Variation within the Collaborative Cross

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

•

The Collaborative Cross models the phenotypic diversity observed in human immunity

•

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

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

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

Bayesian Diallel Analysis Reveals Mx1-Dependent and Mx1-Independent Effects on Response to Influenza A Virus in Mice

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.

Editor's Highlight: Candidate Risk Factors and Mechanisms for Tolvaptan-Induced Liver Injury Are Identified Using a Collaborative Cross Approach

Clinical trials of tolvaptan showed it to be a promising candidate for the treatment of Autosomal Dominant Polycystic Kidney Disease (ADPKD) but also revealed potential for idiosyncratic drug-induced liver injury (DILI) in this patient population. To identify risk factors and mechanisms underlying tolvaptan DILI, 8 mice in each of 45 strains of the genetically diverse Collaborative Cross (CC) mouse population were treated with a single oral dose of either tolvaptan or vehicle. Significant elevations in plasma alanine aminotransferase (ALT) were observed in tolvaptan-treated animals in 3 of the 45 strains. Genetic mapping coupled with transcriptomic analysis in the liver was used to identify several candidate susceptibility genes including epoxide hydrolase 2, interferon regulatory factor 3, and mitochondrial fission factor. Gene pathway analysis revealed that oxidative stress and immune response pathways were activated in response to tolvaptan treatment across all strains, but genes involved in regulation of bile acid homeostasis were most associated with tolvaptan-induced elevations in ALT. Secretory leukocyte peptidase inhibitor (Slpi) mRNA was also induced in the susceptible strains and was associated with increased plasma levels of Slpi protein, suggesting a potential serum marker for DILI susceptibility. In summary, tolvaptan induced signs of oxidative stress, mitochondrial dysfunction, and innate immune response in all strains, but variation in bile acid homeostasis was most associated with susceptibility to the liver response. This CC study has indicated potential mechanisms underlying tolvaptan DILI and biomarkers of susceptibility that may be useful in managing the risk of DILI in ADPKD patients.

Genetic and epigenetic determinants of inter-individual variability in responses to toxicants

It is well established that genetic variability has a major impact on susceptibility to common diseases, responses to drugs and toxicants, and influences disease-related outcomes. The appreciation that epigenetic marks also vary across the population is growing with more data becoming available from studies in humans and model organisms. In addition, the links between genetic variability, toxicity outcomes and epigenetics are being actively explored. Recent studies demonstrate that gene-by-environment interactions involve both chromatin states and transcriptional regulation, and that epigenetics provides important mechanistic clues to connect expression-related quantitative trait loci (QTL) and disease outcomes. However, studies of Gene×Environment×Epigenetics further extend the complexity of the experimental designs and create a challenge for selecting the most informative epigenetic readouts that can be feasibly performed to interrogate multiple individuals, exposures, tissue types and toxicity phenotypes. We propose that among the many possible epigenetic experimental methodologies, assessment of chromatin accessibility coupled with total RNA levels provides a cost-effective and comprehensive option to sufficiently characterize the complexity of epigenetic and regulatory activity in the context of understanding the inter-individual variability in responses to toxicants.

Inbred Strain Variant Database (ISVdb): A Repository for Probabilistically Informed Sequence Differences Among the Collaborative Cross Strains and Their Founders

The Collaborative Cross (CC) is a panel of recently established multiparental recombinant inbred mouse strains. For the CC, as for any multiparental population (MPP), effective experimental design and analysis benefit from detailed knowledge of the genetic differences between strains. Such differences can be directly determined by sequencing, but until now whole-genome sequencing was not publicly available for individual CC strains. An alternative and complementary approach is to infer genetic differences by combining two pieces of information: probabilistic estimates of the CC haplotype mosaic from a custom genotyping array, and probabilistic variant calls from sequencing of the CC founders. The computation for this inference, especially when performed genome-wide, can be intricate and time-consuming, requiring the researcher to generate nontrivial and potentially error-prone scripts. To provide standardized, easy-to-access CC sequence information, we have developed the Inbred Strain Variant Database (ISVdb). The ISVdb provides, for all the exonic variants from the Sanger Institute mouse sequencing dataset, direct sequence information for CC founders and, critically, the imputed sequence information for CC strains. Notably, the ISVdb also: (1) provides predicted variant consequence metadata; (2) allows rapid simulation of F1 populations; and (3) preserves imputation uncertainty, which will allow imputed data to be refined in the future as additional sequencing and genotyping data are collected. The ISVdb information is housed in an SQL database and is easily accessible through a custom online interface (http://isvdb.unc.edu), reducing the analytic burden on any researcher using the CC.

Allelic Variation in the Toll-Like Receptor Adaptor Protein Ticam2 Contributes to SARS-Coronavirus Pathogenesis in Mice

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.

Characterization of Variability in Toxicokinetics and Toxicodynamics of Tetrachloroethylene Using the Collaborative Cross Mouse Population

BACKGROUND

Evaluation of interindividual variability is a challenging step in risk assessment. For most environmental pollutants, including perchloroethylene (PERC), experimental data are lacking, resulting in default assumptions being used to account for variability in toxicokinetics and toxicodynamics.

OBJECTIVE

We quantitatively examined the relationship between PERC toxicokinetics and toxicodynamics at the population level to test whether individuals with increased oxidative metabolism are be more sensitive to hepatotoxicity following PERC exposure.

METHODS

Male mice from 45 strains of the Collaborative Cross (CC) were orally administered a single dose of PERC (1,000 mg/kg) or vehicle (Alkamuls-EL620) and euthanized at various time points (n = 1/strain/time). Concentration–time profiles were generated for PERC and its primary oxidative metabolite trichloroacetate (TCA) in multiple tissues. Toxicodynamic phenotyping was also performed.

RESULTS

Significant variability among strains was observed in toxicokinetics of PERC and TCA in every tissue examined. Based on area under the curve (AUC), the range of liver TCA levels spanned nearly an order of magnitude (~8-fold). Expression of liver cytochrome P4502E1 did not correlate with TCA levels. Toxicodynamic phenotyping revealed an effect of PERC on bodyweight loss, induction of peroxisome proliferator activated receptor-alpha (PPARα)-regulated genes, and dysregulation of hepatic lipid homeostasis. Clustering was observed among a) liver levels of PERC, TCA, and triglycerides; b) TCA levels in liver and kidney; and c) TCA levels in serum, brain, fat, and lung.

CONCLUSIONS

Using the CC mouse population model, we have demonstrated a complex and highly variable relationship between PERC and TCA toxicokinetics and toxicodynamics at the population level. https://doi.org/10.1289/EHP788.

The Collaborative Cross Resource for Systems Genetics Research of Infectious Diseases

An increasing body of evidence highlights the role of host genetic variation in driving susceptibility to severe disease following pathogen infection. In order to fully appreciate the importance of host genetics on infection susceptibility and resulting disease, genetically variable experimental model systems should be employed. These systems allow for the identification, characterization, and mechanistic dissection of genetic variants that cause differential disease responses. Herein we discuss application of the Collaborative Cross (CC) panel of recombinant inbred strains to study viral pathogenesis, focusing on practical considerations for experimental design, assessment and analysis of disease responses within the CC, as well as some of the resources developed for the CC. Although the focus of this chapter is on viral pathogenesis, many of the methods presented within are applicable to studies of other pathogens, as well as to case-control designs in genetically diverse populations.

Humanized Immunoglobulin Mice: Models for HIV Vaccine Testing and Studying the Broadly Neutralizing Antibody Problem

A vaccine that can effectively prevent HIV-1 transmission remains paramount to ending the HIV pandemic, but to do so, will likely need to induce broadly neutralizing antibody (bnAb) responses. A major technical hurdle toward achieving this goal has been a shortage of animal models with the ability to systematically pinpoint roadblocks to bnAb induction and to rank vaccine strategies based on their ability to stimulate bnAb development. Over the past 6 years, immunoglobulin (Ig) knock-in (KI) technology has been leveraged to express bnAbs in mice, an approach that has enabled elucidation of various B-cell tolerance mechanisms limiting bnAb production and evaluation of strategies to circumvent such processes. From these studies, in conjunction with the wealth of information recently obtained regarding the evolutionary pathways and paratopes/epitopes of multiple bnAbs, it has become clear that the very features of bnAbs desired for their function will be problematic to elicit by traditional vaccine paradigms, necessitating more iterative testing of new vaccine concepts. To meet this need, novel bnAb KI models have now been engineered to express either inferred prerearranged V(D)J exons (or unrearranged germline V, D, or J segments that can be assembled into functional rearranged V(D)J exons) encoding predecessors of mature bnAbs. One encouraging approach that has materialized from studies using such newer models is sequential administration of immunogens designed to bind progressively more mature bnAb predecessors. In this review, insights into the regulation and induction of bnAbs based on the use of KI models will be discussed, as will new Ig KI approaches for higher-throughput production and/or altering expression of bnAbs in vivo, so as to further enable vaccine-guided bnAb induction studies.

Dissecting the Effect of Genetic Variation on the Hepatic Expression of Drug Disposition Genes across the Collaborative Cross Mouse Strains

A central challenge in pharmaceutical research is to investigate genetic variation in response to drugs. The Collaborative Cross (CC) mouse reference population is a promising model for pharmacogenomic studies because of its large amount of genetic variation, genetic reproducibility, and dense recombination sites. While the CC lines are phenotypically diverse, their genetic diversity in drug disposition processes, such as detoxification reactions, is still largely uncharacterized. Here we systematically measured RNA-sequencing expression profiles from livers of 29 CC lines under baseline conditions. We then leveraged a reference collection of metabolic biotransformation pathways to map potential relations between drugs and their underlying expression quantitative trait loci (eQTLs). By applying this approach on proximal eQTLs, including eQTLs acting on the overall expression of genes and on the expression of particular transcript isoforms, we were able to construct the organization of hepatic eQTL-drug connectivity across the CC population. The analysis revealed a substantial impact of genetic variation acting on drug biotransformation, allowed mapping of potential joint genetic effects in the context of individual drugs, and demonstrated crosstalk between drug metabolism and lipid metabolism. Our findings provide a resource for investigating drug disposition in the CC strains, and offer a new paradigm for integrating biotransformation reactions to corresponding variations in DNA sequences.