Search

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. [ABSTRACT FROM AUTHOR]

The development of peanut allergy is due to a combination of genetic and environmental factors, although specific genes have proven difficult to identify. Previously, we reported that peanut-sensitized Collaborative Cross strain CC027/GeniUnc (CC027) mice develop anaphylaxis upon oral challenge to peanut, in contrast to C3H/HeJ (C3H) mice. This study aimed to determine the genetic basis of orally induced anaphylaxis to peanut in CC027 mice. A genetic mapping population between CC027 and C3H mice was designed to identify the genetic factors that drive oral anaphylaxis. A total of 356 CC027xC3H backcrossed mice were generated, sensitized to peanut, then challenged to peanut by oral gavage. Anaphylaxis and peanut-specific IgE were quantified for all mice. T-cell phenotyping was conducted on CC027 mice and 5 additional Collaborative Cross strains. Anaphylaxis to peanut was absent in 77% of backcrossed mice, with 19% showing moderate anaphylaxis and 4% having severe anaphylaxis. There were 8 genetic loci associated with variation in response to peanut challenge—6 associated with anaphylaxis (temperature decrease) and 2 associated with peanut-specific IgE levels. There were 2 major loci that impacted multiple aspects of the severity of acute anaphylaxis, at which the CC027 allele was associated with worse outcome. At one of these loci, CC027 has a private genetic variant in the Themis gene. Consistent with described functions of Themis , we found that CC027 mice have more immature T cells with fewer CD8+, CD4+, and CD4+CD25+CD127− regulatory T cells. Our results demonstrate a key role for Themis in the orally reactive CC027 mouse model of peanut allergy. [ABSTRACT FROM AUTHOR]

Objectives: This study combines two innovative mouse models in a major gene discovery project to assess the influence of host genetics on asbestos related disease (ARD). Conventional genetics studies provided evidence that some susceptibility to mesothelioma is genetic. However, the identification of host modifier genes, the roles they may play, and whether they contribute to disease susceptibility remain unknown. Here we report a study designed to rapidly identify genes associated with mesothelioma susceptibility by combining the Collaborative Cross (CC) resource with the well-characterised MexTAg mesothelioma mouse model. Methods: The CC is a powerful mouse resource that harnesses over 90% of common genetic variation in the mouse species, allowing rapid identification of genes mediating complex traits. MexTAg mice rapidly, uniformly, and predictably develop mesothelioma, but only after asbestos exposure. To assess the influence of host genetics on ARD, we crossed 72 genetically distinct CC mouse strains with MexTAg mice and exposed the resulting CC-MexTAg (CCMT) progeny to asbestos and monitored them for traits including overall survival, the time to ARD onset (latency), the time between ARD onset and euthanasia (disease progression) and ascites volume. We identified phenotype-specific modifier genes associated with these traits and we validated the role of human orthologues in asbestos-induced carcinogenesis using human mesothelioma datasets. Results: We generated 72 genetically distinct CCMT strains and exposed their progeny (2,562 in total) to asbestos. Reflecting the genetic diversity of the CC, there was considerable variation in overall survival and disease latency. Surprisingly, however, there was no variation in disease progression, demonstrating that host genetic factors do have a significant influence during disease latency but have a limited role once disease is established. Quantitative trait loci (QTL) affecting ARD survival/latency were identified on chromosomes 6, 12 and X. Of the 97-protein coding candidate modifier genes that spanned these QTL, eight genes (CPED1, ORS1, NDUFA1, HS1BP3, IL13RA1, LSM8, TES and TSPAN12) were found to significantly affect outcome in both CCMT and human mesothelioma datasets. Conclusion: Host genetic factors affect susceptibility to development of asbestos associated disease. However, following mesothelioma establishment, genetic variation in molecular or immunological mechanisms did not affect disease progression. Identification of multiple candidate modifier genes and their human homologues with known associations in other advanced stage or metastatic cancers highlights the complexity of ARD and may provide a pathway to identify novel therapeutic targets. [ABSTRACT FROM AUTHOR]

Coronaviruses have caused three severe epidemics since the start of the 21st century: SARS, MERS and COVID-19. The severity of the ongoing COVID-19 pandemic and increasing likelihood of future coronavirus outbreaks motivates greater understanding of factors leading to severe coronavirus disease. We screened ten strains from the Collaborative Cross mouse genetic reference panel and identified strains CC006/TauUnc (CC006) and CC044/Unc (CC044) as coronavirus-susceptible and resistant, respectively, as indicated by variable weight loss and lung congestion scores four days post-infection. We generated a genetic mapping population of 755 CC006xCC044 F2 mice and exposed the mice to one of three genetically distinct mouse-adapted coronaviruses: clade 1a SARS-CoV MA15 (n=391), clade 1b SARS-CoV-2 MA10 (n=274), and clade 2 HKU3-CoV MA (n=90). Quantitative trait loci (QTL) mapping in SARS-CoV MA15- and SARS-CoV-2 MA10-infected F2 mice identified genetic loci associated with disease severity. Specifically, we identified seven loci associated with variation in outcome following infection with either virus, including one, HrS43, that is present in both groups. Three of these QTL, including HrS43, were also associated with HKU3-CoV MA outcome. HrS43 overlaps with a QTL previously reported by our lab that is associated with SARS-CoV MA15 outcome in CC011xCC074 F2 mice and is also syntenic with a human chromosomal region associated with severe COVID-19 outcomes in humans GWAS. The results reported here provide: (a) additional support for the involvement of this locus in SARS-CoV MA15 infection, (b) the first conclusive evidence that this locus is associated with susceptibility across the Sarbecovirus subgenus, and (c) demonstration of the relevance of mouse models in the study of coronavirus disease susceptibility in humans.

BackgroundObesity is a serious disease with a complex etiology characterized by overaccumulation of adiposity resulting in detrimental health outcomes. Given the liver's critical role in the biological processes that attenuate adiposity accumulation, elucidating the influence of genetics and dietary patterns on hepatic gene expression is fundamental for improving methods of obesity prevention and treatment. To determine how genetics and diet impact obesity development, mice from 22 strains of the genetically diverse recombinant inbred Collaborative Cross (CC) mouse panel were challenged to either a high-protein or high-fat high-sucrose diet, followed by extensive phenotyping and analysis of hepatic gene expression.ResultsOver 1000 genes differentially expressed by perturbed dietary macronutrient composition were enriched for biological processes related to metabolic pathways. Additionally, over 9000 genes were differentially expressed by strain and enriched for biological process involved in cell adhesion and signaling. Weighted gene co-expression network analysis identified multiple gene clusters (modules) associated with body fat % whose average expression levels were influenced by both dietary macronutrient composition and genetics. Each module was enriched for distinct types of biological functions.ConclusionsGenetic background affected hepatic gene expression in the CC overall, but diet macronutrient differences also altered expression of a specific subset of genes. Changes in macronutrient composition altered gene expression related to metabolic processes, while genetic background heavily influenced a broad range of cellular functions and processes irrespective of adiposity. Understanding the individual role of macronutrient composition, genetics, and their interaction is critical to developing therapeutic strategies and policy recommendations for precision nutrition.

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

Collateral blood flow varies greatly among humans for reasons that remain unclear, resulting in significant differences in ischemic tissue damage. A similarly large variation has also been found in mice that is caused by genetic background-dependent differences in the extent of collateral formation, termed collaterogenesis—a unique angiogenic process that occurs during development and determines collateral number and diameter in the adult. Previous studies have identified several quantitative trait loci (QTL) linked to this variation. However, understanding has been hampered by the use of closely related inbred strains that do not model the wide genetic variation present in the "outbred" human population. The Collaborative Cross (CC) multiparent mouse genetic reference panel was developed to address this limitation. Herein we measured the number and average diameter of cerebral collaterals in 60 CC strains, their 8 founder strains, 8 F1 crosses of CC strains selected for abundant versus sparse collaterals, and 2 intercross populations created from the latter. Collateral number evidenced 47-fold variation among the 60 CC strains, with 14% having poor, 25% poor-to-intermediate, 47% intermediate-to-good, and 13% good collateral abundance, that was associated with large differences in post-stroke infarct volume. Collateral number in skeletal muscle and intestine of selected high- and low-collateral strains evidenced the same relative abundance as in brain. Genome-wide mapping demonstrated that collateral abundance is a highly polymorphic trait. Subsequent analysis identified: 6 novel QTL circumscribing 28 high-priority candidate genes harboring putative loss-of-function polymorphisms (SNPs) associated with low collateral number; 335 predicted-deleterious SNPs present in their human orthologs; and 32 genes associated with vascular development but lacking protein coding variants. Six additional suggestive QTL (LOD > 4.5) were also identified in CC-wide QTL mapping. This study provides a comprehensive set of candidate genes for future investigations aimed at identifying signaling proteins within the collaterogenesis pathway whose variants potentially underlie genetic-dependent collateral insufficiency in brain and other tissues. [ABSTRACT FROM AUTHOR]

Objectives: This study combines two innovative mouse models in a major gene discovery project to assess the influence of host genetics on asbestos related disease (ARD). Conventional genetics studies provided evidence that some susceptibility to mesothelioma is genetic. However, the identification of host modifier genes, the roles they may play, and whether they contribute to disease susceptibility remain unknown. Here we report a study designed to rapidly identify genes associated with mesothelioma susceptibility by combining the Collaborative Cross (CC) resource with the well-characterised MexTAg mesothelioma mouse model.Methods: The CC is a powerful mouse resource that harnesses over 90% of common genetic variation in the mouse species, allowing rapid identification of genes mediating complex traits. MexTAg mice rapidly, uniformly, and predictably develop mesothelioma, but only after asbestos exposure. To assess the influence of host genetics on ARD, we crossed 72 genetically distinct CC mouse strains with MexTAg mice and exposed the resulting CC-MexTAg (CCMT) progeny to asbestos and monitored them for traits including overall survival, the time to ARD onset (latency), the time between ARD onset and euthanasia (disease progression) and ascites volume. We identified phenotype-specific modifier genes associated with these traits and we validated the role of human orthologues in asbestos-induced carcinogenesis using human mesothelioma datasets.Results: We generated 72 genetically distinct CCMT strains and exposed their progeny (2,562 in total) to asbestos. Reflecting the genetic diversity of the CC, there was considerable variation in overall survival and disease latency. Surprisingly, however, there was no variation in disease progression, demonstrating that host genetic factors do have a significant influence during disease latency but have a limited role once disease is established. Quantitative trait loci (QTL) affecting ARD survival/latency were identified on chromosomes 6, 12 and X. Of the 97-protein coding candidate modifier genes that spanned these QTL, eight genes (CPED1, ORS1, NDUFA1, HS1BP3, IL13RA1, LSM8, TES and TSPAN12) were found to significantly affect outcome in both CCMT and human mesothelioma datasets.Conclusion: Host genetic factors affect susceptibility to development of asbestos associated disease. However, following mesothelioma establishment, genetic variation in molecular or immunological mechanisms did not affect disease progression. Identification of multiple candidate modifier genes and their human homologues with known associations in other advanced stage or metastatic cancers highlights the complexity of ARD and may provide a pathway to identify novel therapeutic targets.

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

Blood lead (Pb) level (BLL) is a commonly used biomarker to evaluate associations with health effects. However, interventions to reduce the adverse effects of Pb require relating BLL to external exposure. Moreover, risk mitigation actions need to ensure protection of more susceptible individuals with a greater tendency to accumulate Pb. Because little data is available to quantify inter-individual variability in biokinetics of Pb, we investigated the influence of genetics and diet on BLL in the genetically diverse Collaborative Cross (CC) mouse population. Adult female mice from 49 CC strains received either a standard mouse chow or a chow mimicking the American diet while being provided water ad libitum with 1000 ppm Pb for 4 weeks. In both arms of the study, inter-strain variability was observed; however, in American diet-fed animals, the BLL was greater and more variable. Importantly, the degree of variation in BLL among strains on the American diet was greater (2.3) than the default variability estimate (1.6) used in setting the regulatory standards. Genetic analysis identified suggestive diet-associated haplotypes that were associated with variation in BLL, largely contributed by the PWK/PhJ strain. This study quantified the variation in BLL that is due to genetic background, diet, and their interactions, and observed that it may be greater than that assumed for current regulatory standards for Pb in drinking water. Moreover, this work highlights the need of characterizing inter-individual variation in BLL to ensure adequate public health interventions aimed at reducing human health risks from Pb. [ABSTRACT FROM AUTHOR]

Abstract Background Type 2 diabetes (T2D) is an adult‐onset and obese form of diabetes caused by an interplay between genetic, epigenetic, and environmental components. Here, we have assessed a cohort of 11 genetically different collaborative cross (CC) mouse lines comprised of both sexes for T2D and obesity developments in response to oral infection and high‐fat diet (HFD) challenges. Methods Mice were fed with either the HFD or the standard chow diet (control group) for 12 weeks starting at the age of 8 weeks. At week 5 of the experiment, half of the mice of each diet group were infected with Porphyromonas gingivalis and Fusobacterium nucleatum bacteria strains. Throughout the 12‐week experimental period, body weight (BW) was recorded biweekly, and intraperitoneal glucose tolerance tests were performed at weeks 6 and 12 of the experiment to evaluate the glucose tolerance status of mice. Results Statistical analysis has shown the significance of phenotypic variations between the CC lines, which have different genetic backgrounds and sex effects in different experimental groups. The heritability of the studied phenotypes was estimated and ranged between 0.45 and 0.85. We applied machine learning methods to make an early call for T2D and its prognosis. The results showed that classification with random forest could reach the highest accuracy classification (ACC = 0.91) when all the attributes were used. Conclusion Using sex, diet, infection status, initial BW, and area under the curve (AUC) at week 6, we could classify the final phenotypes/outcomes at the end stage of the experiment (at 12 weeks).

Abstract Background Phosphorylation of proteins is a key step in the regulation of many cellular processes including activation of enzymes and signaling cascades. The abundance of a phosphorylated peptide (phosphopeptide) is determined by the abundance of its parent protein and the proportion of target sites that are phosphorylated. Results We quantified phosphopeptides, proteins, and transcripts in heart, liver, and kidney tissue samples of mice from 58 strains of the Collaborative Cross strain panel. We mapped ~700 phosphorylation quantitative trait loci (phQTL) across the three tissues and applied genetic mediation analysis to identify causal drivers of phosphorylation. We identified kinases, phosphatases, cytokines, and other factors, including both known and potentially novel interactions between target proteins and genes that regulate site-specific phosphorylation. Our analysis highlights multiple targets of pyruvate dehydrogenase kinase 1 (PDK1), a regulator of mitochondrial function that shows reduced activity in the NZO/HILtJ mouse, a polygenic model of obesity and type 2 diabetes. Conclusions Together, this integrative multi-omics analysis in genetically diverse CC strains provides a powerful tool to identify regulators of protein phosphorylation. The data generated in this study provides a resource for further exploration.

Defined as chronic excessive accumulation of adiposity, obesity results from long-term imbalance between energy intake and expenditure. The mechanisms behind how caloric imbalance occurs are complex and influenced by numerous biological and environmental factors, especially genetics, and diet. Population-based diet recommendations have had limited success partly due to the wide variation in physiological responses across individuals when they consume the same diet. Thus, it is necessary to broaden our understanding of how individual genetics and diet interact relative to the development of obesity for improving weight loss treatment. To determine how consumption of diets with different macronutrient composition alter adiposity and other obesity-related traits in a genetically diverse population, we analyzed body composition, metabolic rate, clinical blood chemistries, and circulating metabolites in 22 strains of mice from the Collaborative Cross (CC), a highly diverse recombinant inbred mouse population, before and after 8 weeks of feeding either a high protein or high fat high sucrose diet. At both baseline and post-diet, adiposity and other obesity-related traits exhibited a broad range of phenotypic variation based on CC strain; diet-induced changes in adiposity and other traits also depended largely on CC strain. In addition to estimating heritability at baseline, we also quantified the effect size of diet for each trait, which varied by trait and experimental diet. Our findings identified CC strains prone to developing obesity, demonstrate the genotypic and phenotypic diversity of the CC for studying complex traits, and highlight the importance of accounting for genetic differences when making dietary recommendations.

The gut‐brain axis is increasingly recognized as an important pathway involved in cocaine use disorder. Microbial products of the murine gut have been shown to affect striatal gene expression, and depletion of the microbiome by antibiotic treatment alters cocaine‐induced behavioral sensitization in C57BL/6J male mice. Some reports suggest that cocaine‐induced behavioral sensitization is correlated with drug self‐administration behavior in mice. Here, we profile the composition of the naïve microbiome and its response to cocaine sensitization in two collaborative cross (CC) strains. These strains display extremely divergent behavioral responses to cocaine sensitization. A high‐responding strain, CC004/TauUncJ (CC04), has a gut microbiome that contains a greater amount of Lactobacillus than the cocaine‐nonresponsive strain CC041/TauUncJ (CC41). The gut microbiome of CC41 is characterized by an abundance of Eisenbergella, Robinsonella and Ruminococcus. In response to cocaine, CC04 has an increased Barnsiella population, while the gut microbiome of CC41 displays no significant changes. PICRUSt functional analysis of the functional potential of the gut microbiome in CC04 shows a significant number of potential gut‐brain modules altered after exposure to cocaine, specifically those encoding for tryptophan synthesis, glutamine metabolism, and menaquinone synthesis (vitamin K2). Depletion of the microbiome by antibiotic treatment revealed an altered cocaine‐sensitization response following antibiotics in female CC04 mice. Depleting the microbiome by antibiotic treatment in males revealed increased infusions for CC04 during a cocaine intravenous self‐administration dose–response curve. Together these data suggest that genetic differences in cocaine‐related behaviors may involve the microbiome. [ABSTRACT FROM AUTHOR]

ObjectiveThe Collaborative Cross (CC) is a mouse population model with diverse and reproducible genetic backgrounds used to identify novel disease models and genes that contribute to human disease. Since spontaneous tumour susceptibility in CC mice remains unexplored, we assessed tumour incidence and spectrum.DesignWe monitored 293 mice from 18 CC strains for tumour development. Genetic association analysis and RNA sequencing were used to identify susceptibility loci and candidate genes. We analysed genomes of patients with gastric cancer to evaluate the relevance of genes identified in the CC mouse model and measured the expression levels of ISG15 by immunohistochemical staining using a gastric adenocarcinoma tissue microarray. Association of gene expression with overall survival (OS) was assessed by Kaplan-Meier analysis.ResultsCC mice displayed a wide range in the incidence and types of spontaneous tumours. More than 40% of CC036 mice developed gastric tumours within 1 year. Genetic association analysis identified Nfκb1 as a candidate susceptibility gene, while RNA sequencing analysis of non-tumour gastric tissues from CC036 mice showed significantly higher expression of inflammatory response genes. In human gastric cancers, the majority of human orthologues of the 166 mouse genes were preferentially altered by amplification or deletion and were significantly associated with OS. Higher expression of the CC036 inflammatory response gene signature is associated with poor OS. Finally, ISG15 protein is elevated in gastric adenocarcinomas and correlated with shortened patient OS.ConclusionsCC strains exhibit tremendous variation in tumour susceptibility, and we present CC036 as a spontaneous laboratory mouse model for studying human gastric tumourigenesis.

Abstract Background Zika virus (ZIKV) infection at postnatal or adult age can lead to neurological disorders associated with cognitive defects. Yet, how mature neurons respond to ZIKV remains substantially unexplored. Methods The impact of ZIKV infection on mature neurons and microglia was analyzed at the molecular and cellular levels, in vitro using immunocompetent primary cultured neurons and microglia, and in vivo in the brain of adult immunocompetent mice following intracranial ZIKV inoculation. We have used C57BL/6 and the genetically diverse Collaborative Cross mouse strains, displaying a broad range of susceptibility to ZIKV infection, to question the correlation between the effects induced by ZIKV infection on neurons and microglia and the in vivo susceptibility to ZIKV. Results As a result of a delayed induction of interferon beta (IFNB) expression and response, infected neurons displayed an inability to stop ZIKV replication, a trait that was further increased in neurons from susceptible mice. Alongside with an enhanced expression of ZIKV RNA, we observed in vivo, in the brain of susceptible mice, an increased level of active Iba1-expressing microglial cells occasionally engulfing neurons and displaying a gene expression profile close to the molecular signature of disease-associated microglia (DAM). In vivo as well as in vitro, only neurons and not microglial cells were identified as infected, raising the question of the mechanisms underlying microglia activation following brain ZIKV infection. Treatment of primary cultured microglia with conditioned media from ZIKV-infected neurons demonstrated that type-I interferons (IFNs-I) secreted by neurons late after infection activate non-infected microglial cells. In addition, ZIKV infection induced pathological phosphorylation of Tau (pTau) protein, a hallmark of neurodegenerative tauopathies, in vitro and in vivo with clusters of neurons displaying pTau surrounded by active microglial cells. Conclusions We show that ZIKV-infected mature neurons display an inability to stop viral replication in link with a delayed IFNB expression and response, while signaling microglia for activation through IFNs-I secreted at late times post-infection. In the brain of ZIKV-infected susceptible mice, uninfected microglial cells adopt an active morphology and a DAM expression profile, surrounding and sometimes engulfing neurons while ZIKV-infected neurons accumulate pTau, overall reflecting a tauopathy-like phenotype.

Neurological dysfunction following viral infection varies among individuals, largely due to differences in their genetic backgrounds. Gait patterns, which can be evaluated using measures of coordination, balance, posture, muscle function, step-to-step variability, and other factors, are also influenced by genetic background. Accordingly, to some extent gait can be characteristic of an individual, even prior to changes in neurological function. Because neuromuscular aspects of gait are under a certain degree of genetic control, the hypothesis tested was that gait parameters could be predictive of neuromuscular dysfunction following viral infection. The Collaborative Cross (CC) mouse resource was utilized to model genetically diverse populations and the DigiGait treadmill system used to provide quantitative and objective measurements of 131 gait parameters in 142 mice from 23 CC and SJL/J strains. DigiGait measurements were taken prior to infection with the neurotropic virus Theiler's Murine Encephalomyelitis Virus (TMEV). Neurological phenotypes were recorded over 90 days post-infection (d.p.i.), and the cumulative frequency of the observation of these phenotypes was statistically associated with discrete baseline DigiGait measurements. These associations represented spatial and postural aspects of gait influenced by the 90 d.p.i. phenotype score. Furthermore, associations were found between these gait parameters with sex and outcomes considered to show resistance, resilience, or susceptibility to severe neurological symptoms after long-term infection. For example, higher pre-infection measurement values for the Paw Drag parameter corresponded with greater disease severity at 90 d.p.i. Quantitative trait loci significantly associated with these DigiGait parameters revealed potential relationships between 28 differentially expressed genes (DEGs) and different aspects of gait influenced by viral infection. Thus, these potential candidate genes and genetic variations may be predictive of long-term neurological dysfunction. Overall, these findings demonstrate the predictive/prognostic value of quantitative and objective pre-infection DigiGait measurements for viral-induced neuromuscular dysfunction. [ABSTRACT FROM AUTHOR]

BackgroundEarly dietary introduction of peanut has shown efficacy in clinical trials and driven pediatric recommendations for early introduction of peanut to children with heightened allergy risk worldwide. Unfortunately, tolerance is not induced in every case, and a subset of patients are allergic prior to introduction. Here we assess peanut allergic sensitization and oral tolerance in genetically diverse mouse strains.ObjectiveWe aimed to determine whether environmental adjuvant-driven airway sensitization and oral tolerance to peanut could be induced in various genetically diverse mouse strains.MethodsC57BL/6J and 12 Collaborative Cross (CC) mouse strains were fed regular chow or ad libitum peanut butter to induce tolerance. Tolerance was tested by attempting to sensitize mice via intratracheal exposure to peanut and lipopolysaccharide (LPS), followed by intraperitoneal peanut challenge. Peanut-specific immunoglobulins and peanut-induced anaphylaxis were assessed.ResultsWithout oral peanut feeding, most CC strains (11/12) and C57BL/6J induced peanut-specific IgE and IgG1 following airway exposure to peanut and LPS. With oral peanut feeding none of the CC strains nor C57BL/6J mice became sensitized to peanut or experienced anaphylaxis following peanut challenge.ConclusionAllergic sensitization and oral tolerance to peanut can be achieved across a range of genetically diverse mice. Notably, the same strains that became allergic via airway sensitization were tolerized by feeding high doses of peanut butter before sensitization, suggesting that the order and route of peanut exposure are critical for determining the allergic fate.

Background: Zika virus (ZIKV) infection at postnatal or adult age can lead to neurological disorders associated with cognitive defects. Yet, how mature neurons respond to ZIKV remains substantially unexplored. Methods: The impact of ZIKV infection on mature neurons and microglia was analyzed at the molecular and cellular levels, in vitro using immunocompetent primary cultured neurons and microglia, and in vivo in the brain of adult immunocompetent mice following intracranial ZIKV inoculation. We have used C57BL/6 and the genetically diverse Collaborative Cross mouse strains, displaying a broad range of susceptibility to ZIKV infection, to question the correlation between the effects induced by ZIKV infection on neurons and microglia and the in vivo susceptibility to ZIKV. Results: As a result of a delayed induction of interferon beta (IFNB) expression and response, infected neurons displayed an inability to stop ZIKV replication, a trait that was further increased in neurons from susceptible mice. Alongside with an enhanced expression of ZIKV RNA, we observed in vivo, in the brain of susceptible mice, an increased level of active Iba1-expressing microglial cells occasionally engulfing neurons and displaying a gene expression profile close to the molecular signature of disease-associated microglia (DAM). In vivo as well as in vitro, only neurons and not microglial cells were identified as infected, raising the question of the mechanisms underlying microglia activation following brain ZIKV infection. Treatment of primary cultured microglia with conditioned media from ZIKV-infected neurons demonstrated that type-I interferons (IFNs-I) secreted by neurons late after infection activate non-infected microglial cells. In addition, ZIKV infection induced pathological phosphorylation of Tau (pTau) protein, a hallmark of neurodegenerative tauopathies, in vitro and in vivo with clusters of neurons displaying pTau surrounded by active microglial cells. Conclusions: We show that ZIKV-infected mature neurons display an inability to stop viral replication in link with a delayed IFNB expression and response, while signaling microglia for activation through IFNs-I secreted at late times post-infection. In the brain of ZIKV-infected susceptible mice, uninfected microglial cells adopt an active morphology and a DAM expression profile, surrounding and sometimes engulfing neurons while ZIKV-infected neurons accumulate pTau, overall reflecting a tauopathy-like phenotype. [ABSTRACT FROM AUTHOR]

Non-alcoholic fatty liver disease (NAFLD) is a progressive liver disease, affecting 38% of adults globally. If left untreated, NAFLD may progress to more advanced forms of the disease, including non-alcoholic steatohepatitis (NASH), liver cirrhosis, and fibrosis. Early NAFLD detection is critical to prevent disease progression. Using an obesogenic high-fat and high-sucrose (HF/HS) diet, we characterized the progression of NAFLD in male and female Collaborative Cross CC042 mice after 20-, 40-, and 60-week intervals of chronic HF/HS diet feeding. The incidence and severity of liver steatosis, inflammation, and fibrosis increased in both sexes over time, with male mice progressing to a NASH-like disease state faster than female mice, as indicated by earlier and more pronounced changes in liver steatosis. Histopathological indication of macrovesicular steatosis and gene expression changes of key lipid metabolism genes were found to be elevated in both sexes after 20 weeks of HF/HS diet. Measurement of circulating markers of inflammation (CXCL10 and TNF-α), histopathological analysis of immune cell infiltrates, and gene expression changes in inflammation-related genes indicated significant liver inflammation after 40 and 60 weeks of HF/HS diet exposure in both sexes. Liver fibrosis, as assessed by Picosirius red and Masson's trichrome staining and changes in expression of key fibrosis related genes indicated significant changes after 40 and 60 weeks of HF/HS diet exposure. In conclusion, we present a preclinical animal model of dietary NAFLD progression, which recapitulates human pathophysiological and pathomorphological changes, that could be used to better understand the progression of NAFLD and support development of new therapeutics. • Chronically feeding mice a high-fat and high-sucrose diet induced changes consistent with clinical NAFLD. • The model exhibits histological indication of steatosis, inflammation, and fibrosis in both sexes. • Gene expression changes support NASH-like histopathological progression. • This preclinical NASH model could be used for future drug therapy studies. [ABSTRACT FROM AUTHOR]

Coronaviruses have caused three severe epidemics since the start of the 21st century: SARS, MERS and COVID-19. The severity of the ongoing COVID-19 pandemic and increasing likelihood of future coronavirus outbreaks motivates greater understanding of factors leading to severe coronavirus disease. We screened ten strains from the Collaborative Cross mouse genetic reference panel and identified strains CC006/TauUnc (CC006) and CC044/Unc (CC044) as coronavirus-susceptible and resistant, respectively, as indicated by variable weight loss and lung congestion scores four days post-infection. We generated a genetic mapping population of 755 CC006xCC044 F2 mice and exposed the mice to one of three genetically distinct mouse-adapted coronaviruses: clade 1a SARS-CoV MA15 (n=391), clade 1b SARS-CoV-2 MA10 (n=274), and clade 2 HKU3-CoV MA (n=90). Quantitative trait loci (QTL) mapping in SARS-CoV MA15- and SARS-CoV-2 MA10-infected F2 mice identified genetic loci associated with disease severity. Specifically, we identified seven loci associated with variation in outcome following infection with either virus, including one, HrS43 , that is present in both groups. Three of these QTL, including HrS43 , were also associated with HKU3-CoV MA outcome. HrS43 overlaps with a QTL previously reported by our lab that is associated with SARS-CoV MA15 outcome in CC011xCC074 F2 mice and is also syntenic with a human chromosomal region associated with severe COVID-19 outcomes in humans GWAS. The results reported here provide: (a) additional support for the involvement of this locus in SARS-CoV MA15 infection, (b) the first conclusive evidence that this locus is associated with susceptibility across the Sarbecovirus subgenus, and (c) demonstration of the relevance of mouse models in the study of coronavirus disease susceptibility in humans. [Display omitted] [ABSTRACT FROM AUTHOR]

Cocaine use and overdose deaths attributed to cocaine have increased significantly in the United States in the last 10 years. Despite the prevalence of cocaine use disorder (CUD) and the personal and societal problems it presents, there are currently no approved pharmaceutical treatments. The absence of treatment options is due, in part, to our lack of knowledge about the etiology of CUDs. There is ample evidence that genetics plays a role in increasing CUD risk but thus far, very few risk genes have been identified in human studies. Genetic studies in mice have been extremely useful for identifying genetic loci and genes, but have been limited to very few genetic backgrounds, leaving substantial phenotypic, and genetic diversity unexplored. Herein we report the measurement of cocaine-induced behavioral sensitization using a 19-day protocol that captures baseline locomotor activity, initial locomotor response to an acute exposure to cocaine and locomotor sensitization across 5 exposures to the drug. These behaviors were measured in 51 genetically diverse Collaborative Cross (CC) strains along with their inbred founder strains. The CC was generated by crossing eight genetically diverse inbred strains such that each inbred CC strain has genetic contributions from each of the founder strains. Inbred CC mice are infinitely reproducible and provide a stable, yet diverse genetic platform on which to study the genetic architecture and genetic correlations among phenotypes. We have identified significant differences in cocaine locomotor sensitivity and behavioral sensitization across the panel of CC strains and their founders. We have established relationships among cocaine sensitization behaviors and identified extreme responding strains that can be used in future studies aimed at understanding the genetic, biological, and pharmacological mechanisms that drive addiction-related behaviors. Finally, we have determined that these behaviors exhibit relatively robust heritability making them amenable to future genetic mapping studies to identify addiction risk genes and genetic pathways that can be studied as potential targets for the development of novel therapeutics. [ABSTRACT FROM AUTHOR]

Brain and skull tissues interact through molecular signalling and mechanical forces during head development, leading to a strong correlation between the neurocranium and the external brain surface. Therefore, when brain tissue is unavailable, neurocranial endocasts are often used to approximate brain size and shape. Evolutionary changes in brain morphology may have resulted in secondary changes to neurocranial morphology, but the developmental and genetic processes underlying this relationship are not well understood. Using automated phenotyping methods, we quantified the genetic basis of endocast variation across large genetically varied populations of laboratory mice in two ways: (1) to determine the contributions of various genetic factors to neurocranial form and (2) to help clarify whether a neurocranial variation is based on genetic variation that primarily impacts bone development or on genetic variation that primarily impacts brain development, leading to secondary changes in bone morphology. Our results indicate that endocast size is highly heritable and is primarily determined by additive genetic factors. In addition, a non‐additive inbreeding effect led to founder strains with lower neurocranial size, but relatively large brains compared to skull size; suggesting stronger canalization of brain size and/or a general allometric effect. Within an outbred sample of mice, we identified a locus on mouse chromosome 1 that is significantly associated with variation in several positively correlated endocast size measures. Because the protein‐coding genes at this locus have been previously associated with brain development and not with bone development, we propose that genetic variation at this locus leads primarily to variation in brain volume that secondarily leads to changes in neurocranial globularity. We identify a strain‐specific missense mutation within Akt3 that is a strong causal candidate for this genetic effect. Whilst it is not appropriate to generalize our hypothesis for this single locus to all other loci that also contribute to the complex trait of neurocranial skull morphology, our results further reveal the genetic basis of neurocranial variation and highlight the importance of the mechanical influence of brain growth in determining skull morphology. [ABSTRACT FROM AUTHOR]

Cocaine use and overdose deaths attributed to cocaine have increased significantly in the United States in the last 10 years. Despite the prevalence of cocaine use disorder (CUD) and the personal and societal problems it presents, there are currently no approved pharmaceutical treatments. The absence of treatment options is due, in part, to our lack of knowledge about the etiology of CUDs. There is ample evidence that genetics plays a role in increasing CUD risk but thus far, very few risk genes have been identified in human studies. Genetic studies in mice have been extremely useful for identifying genetic loci and genes, but have been limited to very few genetic backgrounds, leaving substantial phenotypic, and genetic diversity unexplored. Herein we report the measurement of cocaine-induced behavioral sensitization using a 19-day protocol that captures baseline locomotor activity, initial locomotor response to an acute exposure to cocaine and locomotor sensitization across 5 exposures to the drug. These behaviors were measured in 51 genetically diverse Collaborative Cross (CC) strains along with their inbred founder strains. The CC was generated by crossing eight genetically diverse inbred strains such that each inbred CC strain has genetic contributions from each of the founder strains. Inbred CC mice are infinitely reproducible and provide a stable, yet diverse genetic platform on which to study the genetic architecture and genetic correlations among phenotypes. We have identified significant differences in cocaine locomotor sensitivity and behavioral sensitization across the panel of CC strains and their founders. We have established relationships among cocaine sensitization behaviors and identified extreme responding strains that can be used in future studies aimed at understanding the genetic, biological, and pharmacological mechanisms that drive addiction-related behaviors. Finally, we have determined that these behaviors exhibit relatively robust heritability making them amenable to future genetic mapping studies to identify addiction risk genes and genetic pathways that can be studied as potential targets for the development of novel therapeutics.

ABSTRACT Infectious diseases have shaped the human population genetic structure, and genetic variation influences the susceptibility to many viral diseases. However, a variety of challenges have made the implementation of traditional human Genome-wide Association Studies (GWAS) approaches to study these infectious outcomes challenging. In contrast, mouse models of infectious diseases provide an experimental control and precision, which facilitates analyses and mechanistic studies of the role of genetic variation on infection. Here we use a genetic mapping cross between two distinct Collaborative Cross mouse strains with respect to severe acute respiratory syndrome coronavirus (SARS-CoV) disease outcomes. We find several loci control differential disease outcome for a variety of traits in the context of SARS-CoV infection. Importantly, we identify a locus on mouse chromosome 9 that shows conserved synteny with a human GWAS locus for SARS-CoV-2 severe disease. We follow-up and confirm a role for this locus, and identify two candidate genes, CCR9 and CXCR6, that both play a key role in regulating the severity of SARS-CoV, SARS-CoV-2, and a distantly related bat sarbecovirus disease outcomes. As such we provide a template for using experimental mouse crosses to identify and characterize multitrait loci that regulate pathogenic infectious outcomes across species. IMPORTANCE Host genetic variation is an important determinant that predicts disease outcomes following infection. In the setting of highly pathogenic coronavirus infections genetic determinants underlying host susceptibility and mortality remain unclear. To elucidate the role of host genetic variation on sarbecovirus pathogenesis and disease outcomes, we utilized the Collaborative Cross (CC) mouse genetic reference population as a model to identify susceptibility alleles to SARS-CoV and SARS-CoV-2 infections. Our findings reveal that a multitrait loci found in chromosome 9 is an important regulator of sarbecovirus pathogenesis in mice. Within this locus, we identified and validated CCR9 and CXCR6 as important regulators of host disease outcomes. Specifically, both CCR9 and CXCR6 are protective against severe SARS-CoV, SARS-CoV-2, and SARS-related HKU3 virus disease in mice. This chromosome 9 multitrait locus may be important to help identify genes that regulate coronavirus disease outcomes in humans.

ABSTRACT 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. IMPORTANCE To survive a bacterial infection, an infected host can display resistance or tolerance. Resistance is indicated by a decrease in pathogen load, while for tolerance a high pathogen load is accompanied by minimal disease. We infected genetically diverse mice with Salmonella Typhimurium for 21 days and discovered new phenotypes for disease outcome (delayed susceptible, tolerant, and resistant). Tolerant strains showed the lowest preinfection core body temperatures and the most rapid disruption in circadian patterns of body temperature postinfection. Tolerant strains had higher circulating neutrophils and higher circulating levels of MCP-1 and IFN-γ, but lower levels of ENA-78 than did resistant strains, in addition to more severe tissue damage. QTL analysis revealed multiple associated regions, and gene expression analysis identified 22 genes that are differentially regulated in tolerant versus resistant animals in these regions. Fibrinogen genes (Fga, Fgb, and Fgg) were found across the QTL, RNA, and the top canonical pathways, suggesting a role in tolerance.

Abstract Background Familial adenomatous polyposis is an inherited genetic disease, characterized by colorectal polyps. It is caused by inactivating mutations in the Adenomatous polyposis coli (Apc) gene. Mice carrying a nonsense mutation in the Apc gene at R850, which is designated Apc Min/+ (Multiple intestinal neoplasia), develop intestinal adenomas. Several genetic modifier loci of Min (Mom) were previously mapped, but so far, most of the underlying genes have not been identified. To identify novel modifier loci associated with Apc Min/+ , we performed quantitative trait loci (QTL) analysis for polyp development using 49 F1 crosses between different Collaborative Cross (CC) lines and C57BL/6 J-Apc Min/+ mice. The CC population is a genetic reference panel of recombinant inbred lines, each line independently descended from eight genetically diverse founder strains. C57BL/6 J-Apc Min/+ males were mated with females from 49 CC lines. F1 offspring were terminated at 23 weeks and polyp counts from three sub-regions (SB1–3) of small intestinal and colon were recorded. Results The number of polyps in all these sub-regions and colon varied significantly between the different CC lines. At 95% genome-wide significance, we mapped nine novel QTL for variation in polyp number, with distinct QTL associated with each intestinal sub-region. QTL confidence intervals varied in width between 2.63–17.79 Mb. We extracted all genes in the mapped QTL at 90 and 95% CI levels using the BioInfoMiner online platform to extract, significantly enriched pathways and key linker genes, that act as regulatory and orchestrators of the phenotypic landscape associated with the Apc Min/+ mutation. Conclusions Genomic structure of the CC lines has allowed us to identify novel modifiers and confirmed some of the previously mapped modifiers. Key genes involved mainly in metabolic and immunological processes were identified. Future steps in this analysis will be to identify regulatory elements – and possible epistatic effects – located in the mapped QTL.

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. [Display omitted] • Collaborative Cross and Diversity Outbred mice show dendritic cell (DC) frequency variation • QTL mapping reveals the highly polygenic and pleiotropic architecture of DC homeostasis • Gramd4 and Orai3 regulate DC frequency Oliveira et al. use quantitative trait locus (QTL) mapping in Diversity Outbred mice to find candidate genes linked to dendritic cell (DC) homeostasis. Site-directed mutagenesis using CRISPR-Cas9 confirms two candidate genes, Gramd4 and Orai3. Overall, the data represent a resource for interrogating the mechanisms governing DC homeostasis in tissues. [ABSTRACT FROM AUTHOR]

Neurological dysfunction following viral infection varies among individuals, largely due to differences in their genetic backgrounds. Gait patterns, which can be evaluated using measures of coordination, balance, posture, muscle function, step-to-step variability, and other factors, are also influenced by genetic background. Accordingly, to some extent gait can be characteristic of an individual, even prior to changes in neurological function. Because neuromuscular aspects of gait are under a certain degree of genetic control, the hypothesis tested was that gait parameters could be predictive of neuromuscular dysfunction following viral infection. The Collaborative Cross (CC) mouse resource was utilized to model genetically diverse populations and the DigiGait treadmill system used to provide quantitative and objective measurements of 131 gait parameters in 142 mice from 23 CC and SJL/J strains. DigiGait measurements were taken prior to infection with the neurotropic virus Theiler’s Murine Encephalomyelitis Virus (TMEV). Neurological phenotypes were recorded over 90 days post-infection (d.p.i.), and the cumulative frequency of the observation of these phenotypes was statistically associated with discrete baseline DigiGait measurements. These associations represented spatial and postural aspects of gait influenced by the 90 d.p.i. phenotype score. Furthermore, associations were found between these gait parameters with sex and outcomes considered to show resistance, resilience, or susceptibility to severe neurological symptoms after long-term infection. For example, higher pre-infection measurement values for the Paw Drag parameter corresponded with greater disease severity at 90 d.p.i. Quantitative trait loci significantly associated with these DigiGait parameters revealed potential relationships between 28 differentially expressed genes (DEGs) and different aspects of gait influenced by viral infection. Thus, these potential candidate genes and genetic variations may be predictive of long-term neurological dysfunction. Overall, these findings demonstrate the predictive/prognostic value of quantitative and objective pre-infection DigiGait measurements for viral-induced neuromuscular dysfunction.

A common step in the analysis of multiparent populations (MPPs) is genotype reconstruction: identifying the founder origin of haplotypes from dense marker data. This process often makes use of a probability model for the pattern of founder alleles along chromosomes, including the relative frequency of founder alleles and the probability of exchanges among them, which depend on a model for meiotic recombination and on the mating design for the population. While the precise experimental design used to generate the population may be used to derive a precise characterization of the model for exchanges among founder alleles, this can be tedious, particularly given the great variety of experimental designs that have been proposed. We describe an approximate model that can be applied for a variety of MPPs. We have implemented the approach in the R/qtl2 software, and we illustrate its use in applications to publicly available data on Diversity Outbred and Collaborative Cross mice. [ABSTRACT FROM AUTHOR]

Abstract Background The Collaborative Cross (CC) mouse population is a valuable resource to study the genetic basis of complex traits, such as obesity. Although the development of obesity is influenced by environmental factors, underlying genetic mechanisms play a crucial role in the response to these factors. The interplay between the genetic background and the gene expression pattern can provide further insight into this response, but we lack robust and easily reproducible workflows to integrate genomic and transcriptomic information in the CC mouse population. Results We established an automated and reproducible integrative workflow to analyse complex traits in the CC mouse genetic reference panel at the genomic and transcriptomic levels. We implemented the analytical workflow to assess the underlying genetic mechanisms of host susceptibility to diet induced obesity and integrated these results with diet induced changes in the hepatic gene expression of susceptible and resistant mice. Hepatic gene expression differs significantly between obese and non-obese mice, with a significant sex effect, where male and female mice exhibit different responses and coping mechanisms. Conclusion Integration of the data showed that different genes but similar pathways are involved in the genetic susceptibility and disturbed in diet induced obesity. Genetic mechanisms underlying susceptibility to high-fat diet induced obesity are different in female and male mice. The clear distinction we observed in the systemic response to the high-fat diet challenge and to obesity between male and female mice points to the need for further research into distinct sex-related mechanisms in metabolic disease.

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.

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. [ABSTRACT FROM AUTHOR]

In species with single-locus, chromosome-based mechanisms of sex determination, the laws of segregation predict an equal ratio of females to males at birth. Here, we show that departures from this Mendelian expectation are commonplace in the 8-way recombinant inbred Collaborative Cross (CC) mouse population. More than one-third of CC strains exhibit significant sex ratio distortion (SRD) at wean, with twice as many male-biased than female-biased strains. We show that these pervasive sex biases persist across multiple breeding environments, are stable over time, and are not mediated by random maternal effects. SRD exhibits a heritable component, but QTL mapping analyses fail to nominate any large effect loci. These findings, combined with the reported absence of sex ratio biases in the CC founder strains, suggest that SRD manifests from multilocus combinations of alleles only uncovered in recombined CC genomes. We explore several potential complex genetic mechanisms for SRD, including allelic interactions leading to sex-biased lethality, genetic sex reversal, chromosome drive mediated by sex-linked selfish elements, and incompatibilities between specific maternal and paternal genotypes. We show that no one mechanism offers a singular explanation for this population-wide SRD. Instead, our data present preliminary evidence for the action of distinct mechanisms of SRD at play in different strains. Taken together, our work exposes the pervasiveness of SRD in the CC population and nominates the CC as a powerful resource for investigating diverse genetic causes of biased sex chromosome transmission. [ABSTRACT FROM AUTHOR]

Drugs of abuse, including alcohol and stimulants like cocaine, produce effects that are subject to individual variability, and genetic variation accounts for at least a portion of those differences. Notably, research in both animal models and human subjects point toward reward sensitivity and impulsivity as being trait characteristics that predict relatively greater positive subjective responses to stimulant drugs. Here we describe use of the eight collaborative cross (CC) founder strains and 38 (reversal learning) or 10 (all other tests) CC strains to examine the heritability of reward sensitivity and impulsivity traits, as well as genetic correlations between these measures and existing addiction‐related phenotypes. Strains were all tested for activity in an open field and reward sensitivity (intake of chocolate BOOST®). Mice were then divided into two counterbalanced groups and underwent reversal learning (impulsive action and waiting impulsivity) or delay discounting (impulsive choice). CC and founder mice show significant heritability for impulsive action, impulsive choice, waiting impulsivity, locomotor activity, and reward sensitivity, with each impulsive phenotype determined to be non‐correlating, independent traits. This research was conducted within the broader, inter‐laboratory effort of the Center for Systems Neurogenetics of Addiction (CSNA) to characterize CC and DO mice for multiple, cocaine abuse related traits. These data will facilitate the discovery of genetic correlations between predictive traits, which will then guide discovery of genes and genetic variants that contribute to addictive behaviors. [ABSTRACT FROM AUTHOR]

Dopaminergic neurons (DA neurons) are controlled by multiple factors, many involved in neurological disease. Parkinson's disease motor symptoms are caused by the demise of nigral DA neurons, leading to loss of striatal dopamine (DA). Here, we measured DA concentration in the dorsal striatum of 32 members of Collaborative Cross (CC) family and their eight founder strains. Striatal DA varied greatly in founders, and differences were highly heritable in the inbred CC progeny. We identified a locus, containing 164 genes, linked to DA concentration in the dorsal striatum on chromosome X. We used RNAseq profiling of the ventral midbrain of two founders with substantial difference in striatal DA–C56BL/6 J and A/J—to highlight potential protein‐coding candidates modulating this trait. Among the five differentially expressed genes within the locus, we found that the gene coding for the collagen IV alpha 6 chain (Col4a6) was expressed nine times less in A/J than in C57BL/6J. Using single cell RNA‐seq data from developing human midbrain, we found that COL4A6 is highly expressed in radial glia‐like cells and neuronal progenitors, indicating a role in neuronal development. Collagen IV alpha‐6 chain (COL4A6) controls axogenesis in simple model organisms. Consistent with these findings, A/J mice had less striatal axonal branching than C57BL/6J mice. We tentatively conclude that DA concentration and axonal branching in dorsal striatum are modulated by COL4A6, possibly during development. Our study shows that genetic mapping based on an easily measured Central Nervous System (CNS) trait, using the CC population, combined with follow‐up observations, can parse heritability of such a trait, and nominate novel functions for commonly expressed proteins. [ABSTRACT FROM AUTHOR]

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.

Summary: Genetically diverse mouse populations are powerful tools for characterizing the regulation of the proteome and its relationship to whole-organism phenotypes. We used mass spectrometry to profile and quantify the abundance of 6,798 proteins in liver tissue from mice of both sexes across 58 Collaborative Cross (CC) inbred strains. We previously collected liver proteomics data from the related Diversity Outbred (DO) mice and their founder strains. We show concordance across the proteomics datasets despite being generated from separate experiments, allowing comparative analysis. We map protein abundance quantitative trait loci (pQTLs), identifying 1,087 local and 285 distal in the CC mice and 1,706 local and 414 distal in the DO mice. We find that regulatory effects on individual proteins are conserved across the mouse populations, in particular for local genetic variation and sex differences. In comparison, proteins that form complexes are often co-regulated, displaying varying genetic architectures, and overall show lower heritability and map fewer pQTLs. We have made this resource publicly available to enable quantitative analyses of the regulation of the proteome.

A wide range of viruses cause neurological manifestations in their hosts. Infection by neurotropic viruses as well as the resulting immune response can irreversibly disrupt the complex structural and functional architecture of the brain, depending in part on host genetic background. The interaction between host genetic background, neurological response to viral infection, and subsequent clinical manifestations remains poorly understood. In the present study, we used the genetically diverse Collaborative Cross (CC) mouse resource to better understand how differences in genetic background drive clinical signs and neuropathological manifestations of acute Theiler’s murine encephalomyelitis virus (TMEV) infection. For the first time, we characterized variations of TMEV viral tropism and load based on host genetic background, and correlated viral load with microglial/macrophage activation. For five CC strains (CC002, CC023, CC027, CC057, and CC078) infected with TMEV, we compared clinical signs, lesion distribution, microglial/macrophage response, expression, and distribution of TMEV mRNA, and identified genetic loci relevant to the early acute (4 days post-infection [dpi]) and late acute (14 dpi) timepoints. We examined brain pathology to determine possible causes of strain-specific differences in clinical signs, and found that fields CA1 and CA2 of the hippocampal formation were especially targeted by TMEV across all strains. Using Iba-1 immunolabeling, we identified and characterized strain- and timepoint-specific variation in microglial/macrophage reactivity in the hippocampal formation. Because viral clearance can influence disease outcome, we used RNA in situ hybridization to quantify viral load and TMEV mRNA distribution at both timepoints. TMEV mRNA expression was broadly distributed in the hippocampal formation at 4 dpi in all strains but varied between radiating and clustered distribution depending on the CC strain. We found a positive correlation between microglial/macrophage reactivity and TMEV mRNA expression at 4 dpi. At 14 dpi, we observed a dramatic reduction in TMEV mRNA expression, and localization to the medial portion of field CA1 and field CA2. To better understand how host genetic background can influence pathological outcomes, we identified quantitative trait loci associated with frequency of lesions in a particular brain region and with microglial/macrophage reactivity. These QTL were located near several loci of interest: lysosomal trafficking regulator (Lyst) and nidogen 1 (Nid1), and transmembrane protein 106 B (Tmem106b). Together, these results provide a novel understanding about the influences of genetic variation on the acute neuropathological and immunopathological environment and viral load, which collectively lead to variable disease outcomes. Our findings reveal possible avenues for future investigation which may lead to more effective intervention strategies and treatment regimens.

Thoracic cancers pose a significant global health burden. Immune checkpoint blockade therapies have improved treatment outcomes, but durable responses remain limited. Understanding how the host immune system interacts with a developing tumor is essential for the rational development of improved treatments for thoracic malignancies. Recent technical advances have improved our understanding of the mutational burden of cancer cells and changes in cancer-specific gene expression, providing a detailed understanding of the complex biology underpinning tumor-host interactions. While there has been much focus on the genetic alterations associated with cancer cells and how they may impact treatment outcomes, how host genetics affects cancer development is also critical and will greatly determine treatment response. Genome-wide association studies (GWAS) have identified genetic variants associated with cancer predisposition. This approach has successfully identified host genetic risk factors associated with common thoracic cancers like lung cancer, but is less effective for rare cancers like malignant mesothelioma. To assess how host genetics impacts rare thoracic cancers, we used the Collaborative Cross (CC); a powerful murine genetic resource designed to maximize genetic diversity and rapidly identify genes associated with any biological trait. We are using the CC in conjunction with our asbestos-induced MexTAg mouse model, to identify host genes associated with mesothelioma development. Once genes that moderate tumor development and progression are known, human homologues can be identified and human datasets interrogated to validate their association with disease outcome. Furthermore, our CC−MexTAg animal model enables in-depth study of the tumor microenvironment, allowing the correlation of immune cell infiltration and gene expression signatures with disease development. This strategy provides a detailed picture of the underlying biological pathways associated with mesothelioma susceptibility and progression; knowledge that is crucial for the rational development of new diagnostic and therapeutic strategies. Here we discuss the influence of host genetics on developing an effective immune response to thoracic cancers. We highlight current knowledge gaps, and with a focus on mesothelioma, describe the development and application of the CC-MexTAg to overcome limitations and illustrate how the knowledge gained from this unique study will inform the rational design of future treatments of mesothelioma.

Background: Familial adenomatous polyposis is an inherited genetic disease, characterized by colorectal polyps. It is caused by inactivating mutations in the Adenomatous polyposis coli (Apc) gene. Mice carrying a nonsense mutation in the Apc gene at R850, which is designated ApcMin/+ (Multiple intestinal neoplasia), develop intestinal adenomas. Several genetic modifier loci of Min (Mom) were previously mapped, but so far, most of the underlying genes have not been identified. To identify novel modifier loci associated with ApcMin/+, we performed quantitative trait loci (QTL) analysis for polyp development using 49 F1 crosses between different Collaborative Cross (CC) lines and C57BL/6 J-ApcMin/+mice. The CC population is a genetic reference panel of recombinant inbred lines, each line independently descended from eight genetically diverse founder strains. C57BL/6 J-ApcMin/+ males were mated with females from 49 CC lines. F1 offspring were terminated at 23 weeks and polyp counts from three sub-regions (SB1–3) of small intestinal and colon were recorded. Results: The number of polyps in all these sub-regions and colon varied significantly between the different CC lines. At 95% genome-wide significance, we mapped nine novel QTL for variation in polyp number, with distinct QTL associated with each intestinal sub-region. QTL confidence intervals varied in width between 2.63–17.79 Mb. We extracted all genes in the mapped QTL at 90 and 95% CI levels using the BioInfoMiner online platform to extract, significantly enriched pathways and key linker genes, that act as regulatory and orchestrators of the phenotypic landscape associated with the ApcMin/+ mutation. Conclusions: Genomic structure of the CC lines has allowed us to identify novel modifiers and confirmed some of the previously mapped modifiers. Key genes involved mainly in metabolic and immunological processes were identified. Future steps in this analysis will be to identify regulatory elements – and possible epistatic effects – located in the mapped QTL. [ABSTRACT FROM AUTHOR]

Thoracic cancers pose a significant global health burden. Immune checkpoint blockade therapies have improved treatment outcomes, but durable responses remain limited. Understanding how the host immune system interacts with a developing tumor is essential for the rational development of improved treatments for thoracic malignancies. Recent technical advances have improved our understanding of the mutational burden of cancer cells and changes in cancer-specific gene expression, providing a detailed understanding of the complex biology underpinning tumor-host interactions. While there has been much focus on the genetic alterations associated with cancer cells and how they may impact treatment outcomes, how host genetics affects cancer development is also critical and will greatly determine treatment response. Genome-wide association studies (GWAS) have identified genetic variants associated with cancer predisposition. This approach has successfully identified host genetic risk factors associated with common thoracic cancers like lung cancer, but is less effective for rare cancers like malignant mesothelioma. To assess how host genetics impacts rare thoracic cancers, we used the Collaborative Cross (CC); a powerful murine genetic resource designed to maximize genetic diversity and rapidly identify genes associated with any biological trait. We are using the CC in conjunction with our asbestos-induced MexTAg mouse model, to identify host genes associated with mesothelioma development. Once genes that moderate tumor development and progression are known, human homologues can be identified and human datasets interrogated to validate their association with disease outcome. Furthermore, our CC−MexTAg animal model enables in-depth study of the tumor microenvironment, allowing the correlation of immune cell infiltration and gene expression signatures with disease development. This strategy provides a detailed picture of the underlying biological pathways associated with mesothelioma susceptibility and progression; knowledge that is crucial for the rational development of new diagnostic and therapeutic strategies. Here we discuss the influence of host genetics on developing an effective immune response to thoracic cancers. We highlight current knowledge gaps, and with a focus on mesothelioma, describe the development and application of the CC-MexTAg to overcome limitations and illustrate how the knowledge gained from this unique study will inform the rational design of future treatments of mesothelioma. [ABSTRACT FROM AUTHOR]

Specific fat distributions are risk factors for complex diseases, including coronary heart disease and obstructive sleep apnea. To demonstrate the utility of high-diversity mouse models for elucidating genetic associations, we describe the phenotyping and heritability of fat distributions within the five classical inbred and three wild-derived founder mouse strains of the Collaborative Cross and Diversity Outbred mice. Measurements of subcutaneous and internal fat volumes in the abdomen, thorax and neck, and fat volumes in the tongue and pericardium were obtained using magnetic resonance imaging in male mice from the A/J (n = 12), C57BL/6J (n = 17), 129S1/SvlmJ (n = 12), NOD/LtJ (n = 14), NZO/HILtJ (n = 12), CAST/EiJ (n= 14), PWK/PhJ (n = 12), and WSB/EiJ (n = 15) strains. Phenotypes were compared across strains using analysis of variance and heritability estimated as the proportion of phenotypic variability attributable to strain. Heritability ranged from 44 to 91% across traits, including >70% heritability of tongue fat. A majority of heritability estimates remained significant controlling for body weight, suggesting genetic influences independent of general obesity. Principal components analysis supports genetic influences on overall obesity and specific to increased pericardial and intra-neck fat. Thus, among the founder strains of the Collaborative Cross and Diversity Outbred mice, we observed significant heritability of subcutaneous and internal fat volumes in the neck, thorax and abdomen, pericardial fat volume and tongue fat volume, consistent with genetic architecture playing an important role in explaining trait variability. Findings pave the way for studies utilizing high-diversity mouse models to identify genes affecting fat distributions and, in turn, influencing risk for associated complex disorders. [ABSTRACT FROM AUTHOR]

Mice have provided critical mechanistic understandings of clinical traits underlying metabolic syndrome (MetSyn) and susceptibility to MetSyn in mice is known to vary among inbred strains. We investigated the diet- and strain-dependent effects on metabolic traits in the eight Collaborative Cross (CC) founder strains (A/J, C57BL/6J, 129S1/SvImJ, NOD/ShiLtJ, NZO/HILtJ, CAST/EiJ, PWK/PhJ, and WSB/EiJ). Liver transcriptomics analysis showed that both atherogenic diet and host genetics have profound effects on the liver transcriptome, which may be related to differences in metabolic traits observed between strains. We found strain differences in circulating trimethylamine N-oxide (TMAO) concentration and liver triglyceride content, both of which are traits associated with metabolic diseases. Using a network approach, we identified a module of transcripts associated with TMAO and liver triglyceride content, which was enriched in functional pathways. Interrogation of the module related to metabolic traits identified NADPH oxidase 4 (Nox4), a gene for a key enzyme in the production of reactive oxygen species, which showed a strong association with plasma TMAO and liver triglyceride. Interestingly, Nox4 was identified as the highest expressed in the C57BL/6J and NZO/HILtJ strains and the lowest expressed in the CAST/EiJ strain. Based on these results, we suggest that there may be genetic variation in the contribution of Nox4 to the regulation of plasma TMAO and liver triglyceride content. In summary, we show that liver transcriptomic analysis identified diet- or strain-specific pathways for metabolic traits in the Collaborative Cross (CC) founder strains. [ABSTRACT FROM AUTHOR]

A general Bayesian model, Diploffect, is described for estimating the effects of founder haplotypes at quantitative trait loci (QTL) detected in multiparental genetic populations; such populations include the Collaborative Cross (CC), Heterogeneous Socks (HS), and many others for which local genetic variation is well described by an underlying, usually probabilistically inferred, haplotype mosaic. Our aim is to provide a framework for coherent estimation of haplotype and diplotype (haplotype pair) effects that takes into account the following: uncertainty in haplotype composition for each individual; uncertainty arising from small sample sizes and infrequently observed haplotype combinations; possible effects of dominance (for noninbred subjects); genetic background; and that provides a means to incorporate data that may be incomplete or has a hierarchical structure. Using the results of a probabilistic haplotype reconstruction as prior information, we obtain posterior distributions at the QTL for both haplotype effects and haplotype composition. Two alternative computational approaches are supplied: a Markov chain Monte Carlo sampler and a procedure based on importance sampling of integrated nested Laplace approximations. Using simulations of QTL in the incipient CC (pre-CC) and Northport HS populations, we compare the accuracy of Diploffect, approximations to it, and more commonly used approaches based on Haley-Knott regression, describing trade-offs between these methods. We also estimate effects for three QTL previously identified in those populations, obtaining posterior intervals that describe how the phenotype might be affected by diplotype substitutions at the modeled locus.

Genetic variation between diverse mouse species is well-characterized, yet existing knowledge of the mouse transcriptome comes largely from one mouse strain (C57BL/6J). As such, it is unlikely to reflect the transcriptional complexity of the mouse species. Gene transcription is dynamic and condition-specific; therefore, to better understand the mouse transcriptional response to respiratory virus infection, we infected the eight founder strains of the Collaborative Cross with either influenza A virus or severe acute respiratory syndrome coronavirus and sequenced lung RNA samples at 2 and 4 days after infection. We found numerous instances of transcripts that were not present in the C57BL/6J reference annotation, indicating that a nontrivial proportion of the mouse genome is transcribed but poorly annotated. Of these novel transcripts, 2150 could be aligned to human or rat genomes, but not to existing mouse genomes, suggesting functionally conserved sequences not yet recorded in mouse genomes. We also found that respiratory virus infection induced differential expression of 4287 splicing junctions, resulting in strain-specific isoform expression. Of these, 59 were influenced by strain-specific mutations within 2 base pairs of key intron-exon boundaries, suggesting cis-regulated expression. Our results reveal the complexity of the transcriptional response to viral infection, previously undocumented genomic elements, and extensive diversity in the response across mouse strains. These findings identify hitherto unexplored transcriptional patterns and undocumented transcripts in genetically diverse mice. Host genetic variation drives the complexity and diversity of the host response by eliciting starkly different transcriptional profiles in response to a viral infection.

The Collaborative Cross (CC) is an emerging panel of recombinant inbred (RI) mouse strains. Each strain is genetically distinct but all descended from the same eight inbred founders. In 66 strains from incipient lines of the CC (pre-CC), as well as the 8 CC founders and some of their F1 offspring, we examined subsets of lymphocytes and antigen-presenting cells. We found significant variation among the founders, with even greater diversity in the pre-CC. Genome-wide association using inferred haplotypes detected highly significant loci controlling B-to-T cell ratio, CD8 T-cell numbers, CD11c and CD23 expression. Comparison of overall strain effects in the CC founders with strain effects at QTL in the pre-CC revealed sharp contrasts in the genetic architecture of two traits with significant loci: variation in CD23 can be explained largely by additive genetics at one locus, whereas variation in B-to-T ratio has a more complex etiology. For CD23, we found a strong QTL whose confidence interval contained the CD23 structural gene Fcer2a. Our data on the pre-CC demonstrate the utility of the CC for studying immunophenotypes and the value of integrating founder, CC and F1 data. The extreme immunophenotypes observed could have pleiotropic effects in other CC experiments.

Defined as chronic excessive accumulation of adiposity, obesity results from long-term imbalance between energy intake and expenditure. The mechanisms behind how caloric imbalance occurs are complex and influenced by numerous biological and environmental factors, especially genetics, and diet. Population-based diet recommendations have had limited success partly due to the wide variation in physiological responses across individuals when they consume the same diet. Thus, it is necessary to broaden our understanding of how individual genetics and diet interact relative to the development of obesity for improving weight loss treatment. To determine how consumption of diets with different macronutrient composition alter adiposity and other obesity-related traits in a genetically diverse population, we analyzed body composition, metabolic rate, clinical blood chemistries, and circulating metabolites in 22 strains of mice from the Collaborative Cross (CC), a highly diverse recombinant inbred mouse population, before and after 8 weeks of feeding either a high protein or high fat high sucrose diet. At both baseline and post-diet, adiposity and other obesity-related traits exhibited a broad range of phenotypic variation based on CC strain; diet-induced changes in adiposity and other traits also depended largely on CC strain. In addition to estimating heritability at baseline, we also quantified the effect size of diet for each trait, which varied by trait and experimental diet. Our findings identified CC strains prone to developing obesity, demonstrate the genotypic and phenotypic diversity of the CC for studying complex traits, and highlight the importance of accounting for genetic differences when making dietary recommendations.

Background: Interindividual variation in voluntary ethanol consumption and ethanol response is partially influenced by genetic variation. Discovery of the genes and allelic variants that affect these phenotypes may clarify the etiology and pathophysiology of problematic alcohol use, including alcohol use disorder. Genetically diverse mouse populations, which demonstrate heritable variation in ethanol consumption, can be utilized to discover the genes and gene networks that influence this trait. The Collaborative Cross (CC) recombinant inbred strains, Diversity Outbred (DO) population and their 8 founder strains are complementary mouse resources that capture substantial genetic diversity and can demonstrate expansive phenotypic variation in heritable traits. These populations may be utilized to discover candidate genes and gene networks that moderate ethanol consumption and other ethanol‐related traits. Methods: We characterized ethanol consumption, preference, and pharmacokinetics in the 8 founder strains and 10 CC strains in 12‐hour drinking sessions during the dark phase of the circadian cycle. Results: Ethanol consumption was substantially heritable, both early in ethanol access and over a chronic intermittent access schedule. Ethanol pharmacokinetics were also heritable; however, no association between strain‐level ethanol consumption and pharmacokinetics was detected. The PWK/PhJ strain was the highest drinking strain, with consumption substantially exceeding that of the C57BL/6J strain, which is commonly used as a model of "high" or "binge" drinking. Notably, we found strong evidence that sex moderated genetic effects on voluntary ethanol drinking. Conclusions: Collectively, this research serves as a foundation for expanded genetic study of ethanol consumption in the CC/DO and related populations. Moreover, we identified reference strains with extreme consumption phenotypes that effectively represent polygenic models of excessive ethanol use. [ABSTRACT FROM AUTHOR]