Collaborative Cross mice have diverse phenotypic responses to infection with Methicillin-resistant Staphylococcus aureus USA300.

Staphylococcus aureus (S. aureus) is an opportunistic pathogen causing diseases ranging from mild skin infections to life threatening conditions, including endocarditis, pneumonia, and sepsis. To identify host genes modulating this host-pathogen interaction, we infected 25 Collaborative Cross (CC) mouse strains with methicillin-resistant S. aureus (MRSA) and monitored disease progression for seven days using a surgically implanted telemetry system. CC strains varied widely in their response to intravenous MRSA infection. We identified eight 'susceptible' CC strains with high bacterial load, tissue damage, and reduced survival. Among the surviving strains, six with minimal colonization were classified as 'resistant', while the remaining six tolerated higher organ colonization ('tolerant'). The kidney was the most heavily colonized organ, but liver, spleen and lung colonization were better correlated with reduced survival. Resistant strains had higher pre-infection circulating neutrophils and lower post-infection tissue damage compared to susceptible and tolerant strains. We identified four CC strains with sexual dimorphism: all females survived the study period while all males met our euthanasia criteria earlier. In these CC strains, males had more baseline circulating monocytes and red blood cells. We identified several CC strains that may be useful as new models for endocarditis, myocarditis, pneumonia, and resistance to MRSA infection. Quantitative Trait Locus (QTL) analysis identified two significant loci, on Chromosomes 18 and 3, involved in early susceptibility and late survival after infection. We prioritized Npc1 and Ifi44l genes as the strongest candidates influencing survival using variant analysis and mRNA expression data from kidneys within these intervals.

Genetic and microbial determinants of azoxymethane-induced colorectal tumor susceptibility in Collaborative Cross mice and their implication in human cancer.

The insights into interactions between host genetics and gut microbiome (GM) in colorectal tumor susceptibility (CTS) remains lacking. We used Collaborative Cross mouse population model to identify genetic and microbial determinants of Azoxymethane-induced CTS. We identified 4417 CTS-associated single nucleotide polymorphisms (SNPs) containing 334 genes that were transcriptionally altered in human colorectal cancers (CRCs) and consistently clustered independent human CRC cohorts into two subgroups with different prognosis. We discovered a set of genera in early-life associated with CTS and defined a 16-genus signature that accurately predicted CTS, the majority of which were correlated with human CRCs. We identified 547 SNPs associated with abundances of these genera. Mediation analysis revealed GM as mediators partially exerting the effect of SNP UNC3869242 within Duox2 on CTS. Intestine cell-specific depletion of Duox2 altered GM composition and contribution of Duox2 depletion to CTS was significantly influenced by GM. Our findings provide potential novel targets for personalized CRC prevention and treatment.

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

BACKGROUND: Aspergillus fumigatus (Af) is one of the most ubiquitous fungi and its infection potency is suggested to be strongly controlled by the host genetic background. The aim of this study was to search for candidate genes associated with host susceptibility to Aspergillus fumigatus (Af) using an RNAseq approach in CC lines and hepatic gene expression. METHODS: We studied 31 male mice from 25 CC lines at 8 weeks old; the mice were infected with Af. Liver tissues were extracted from these mice 5 days post-infection, and next-generation RNA-sequencing (RNAseq) was performed. The GENE-E analysis platform was used to generate a clustered heat map matrix. RESULTS: Significant variation in body weight changes between CC lines was observed. Hepatic gene expression revealed 12 top prioritized candidate genes differentially expressed in resistant versus susceptible mice based on body weight changes. Interestingly, three candidate genes are located within genomic intervals of the previously mapped quantitative trait loci (QTL), including Gm16270 and Stox1 on chromosome 10 and Gm11033 on chromosome 8. CONCLUSIONS: Our findings emphasize the CC mouse model's power in fine mapping the genetic components underlying susceptibility towards Af. As a next step, eQTL analysis will be performed for our RNA-Seq data. Suggested candidate genes from our study will be further assessed with a human cohort with aspergillosis.

Genetic variation influences the skeletal response to hindlimb unloading in the eight founder strains of the diversity outbred mouse population.

During disuse, mechanical unloading causes extensive bone loss, decreasing bone volume and strength. Variations in bone mass and risk of osteoporosis are influenced by genetics; however, it remains unclear how genetic variation affects the skeletal response to unloading. We previously found that genetic variation affects the musculoskeletal response to 3 weeks of immobilization in the 8 Jackson Laboratory J:DO founder strains: C57Bl/6J, A/J, 129S1/SvImJ, NOD/ShiLtJ, NZO/HlLtJ, CAST/EiJ, PWK/PhJ, and WSB/EiJ. Hindlimb unloading (HLU) is the best model for simulating local and systemic contributors of disuse and therefore may have a greater impact on bones than immobilization. We hypothesized that genetic variation would affect the response to HLU across the eight founder strains. Mice of each founder strain were placed in HLU for 3 weeks, and the femurs and tibias were analyzed. There were significant HLU and mouse strain interactions on body weight, femur trabecular BV/TV, and femur ultimate force. This indicates that unloading only caused significant catabolic effects in some mouse strains. C57BL/6 J mice were most affected by unloading while other strains were more protected. There were significant HLU and mouse strain interactions on gene expression of genes encoding bone metabolism genes in the tibia. This indicates that unloading only caused significant effects on bone metabolism genes in some mouse strains. Different mouse strains respond to HLU differently, and this can be explained by genetic differences. These results suggest the outbred J:DO mice will be a powerful model for examining the effects of genetics on the skeletal response to HLU.

Mapping novel QTL and fine mapping of previously identified QTL associated with glucose tolerance using the collaborative cross mice.

A chronic metabolic illness, type 2 diabetes (T2D) is a polygenic and multifactorial complicated disease. With an estimated 463 million persons aged 20 to 79 having diabetes, the number is expected to rise to 700 million by 2045, creating a significant worldwide health burden. Polygenic variants of diabetes are influenced by environmental variables. T2D is regarded as a silent illness that can advance for years before being diagnosed. Finding genetic markers for T2D and metabolic syndrome in groups with similar environmental exposure is therefore essential to understanding the mechanism of such complex characteristic illnesses. So herein, we demonstrated the exclusive use of the collaborative cross (CC) mouse reference population to identify novel quantitative trait loci (QTL) and, subsequently, suggested genes associated with host glucose tolerance in response to a high-fat diet. In this study, we used 539 mice from 60 different CC lines. The diabetogenic effect in response to high-fat dietary challenge was measured by the three-hour intraperitoneal glucose tolerance test (IPGTT) test after 12 weeks of dietary challenge. Data analysis was performed using a statistical software package IBM SPSS Statistic 23. Afterward, blood glucose concentration at the specific and between different time points during the IPGTT assay and the total area under the curve (AUC0-180) of the glucose clearance was computed and utilized as a marker for the presence and severity of diabetes. The observed AUC0-180 averages for males and females were 51,267.5 and 36,537.5 mg/dL, respectively, representing a 1.4-fold difference in favor of females with lower AUC0-180 indicating adequate glucose clearance. The AUC0-180 mean differences between the sexes within each specific CC line varied widely within the CC population. A total of 46 QTL associated with the different studied phenotypes, designated as T2DSL and its number, for Type 2 Diabetes Specific Locus and its number, were identified during our study, among which 19 QTL were not previously mapped. The genomic interval of the remaining 27 QTL previously reported, were fine mapped in our study. The genomic positions of 40 of the mapped QTL overlapped (clustered) on 11 different peaks or close genomic positions, while the remaining 6 QTL were unique. Further, our study showed a complex pattern of haplotype effects of the founders, with the wild-derived strains (mainly PWK) playing a significant role in the increase of AUC values.

Identification of genetic drivers of plasma lipoprotein size in the Diversity Outbred mouse population.

Despite great progress in understanding lipoprotein physiology, there is still much to be learned about the genetic drivers of lipoprotein abundance, composition, and function. We used ion mobility spectrometry to survey 16 plasma lipoprotein subfractions in 500 Diversity Outbred mice maintained on a Western-style diet. We identified 21 quantitative trait loci (QTL) affecting lipoprotein abundance. To refine the QTL and link them to disease risk in humans, we asked if the human homologs of genes located at each QTL were associated with lipid traits in human genome-wide association studies. Integration of mouse QTL with human genome-wide association studies yielded candidate gene drivers for 18 of the 21 QTL. This approach enabled us to nominate the gene encoding the neutral ceramidase, Asah2, as a novel candidate driver at a QTL on chromosome 19 for large HDL particles (HDL-2b). To experimentally validate Asah2, we surveyed lipoproteins in Asah2-/- mice. Compared to wild-type mice, female Asah2-/- mice showed an increase in several lipoproteins, including HDL. Our results provide insights into the genetic regulation of circulating lipoproteins, as well as mechanisms by which lipoprotein subfractions may affect cardiovascular disease risk in humans.

Host Genetic Background Effect on Body Weight Changes Influenced by Heterozygous Smad4 Knockout Using Collaborative Cross Mouse Population.

Obesity and its attendant conditions have become major health problems worldwide, and obesity is currently ranked as the fifth most common cause of death globally. Complex environmental and genetic factors are causes of the current obesity epidemic. Diet, lifestyle, chemical exposure, and other confounding factors are difficult to manage in humans. The mice model is helpful in researching genetic BW gain because genetic and environmental risk factors can be controlled in mice. Studies in mouse strains with various genetic backgrounds and established genetic structures provide unparalleled opportunities to find and analyze trait-related genomic loci. In this study, we used the Collaborative Cross (CC), a large panel of recombinant inbred mouse strains, to present a predictive study using heterozygous Smad4 knockout profiles of CC mice to understand and effectively identify predispositions to body weight gain. Male C57Bl/6J Smad4+/- mice were mated with female mice from 10 different CC lines to create F1 mice (Smad4+/-x CC). Body weight (BW) was measured weekly until week 16 and then monthly until the end of the study (week 48). The heritability (H2) of the assessed traits was estimated and presented. Comparative analysis of various machine learning algorithms for predicting the BW changes and genotype of mice was conducted. Our data showed that the body weight records of F1 mice with different CC lines differed between wild-type and mutant Smad4 mice during the experiment. Genetic background affects weight gain and some lines gained more weight in the presence of heterozygous Smad4 knockout, while others gained less, but, in general, the mutation caused overweight mice, except for a few lines. In both control and mutant groups, female %BW had a higher heritability (H2) value than males. Additionally, both sexes with wild-type genotypes showed higher heritability values than the mutant group. Logistic regression provides the most accurate mouse genotype predictions using machine learning. We plan to validate the proposed method on more CC lines and mice per line to expand the literature on machine learning for BW prediction.

Susceptibility to Zika virus in a Collaborative Cross mouse strain is induced by Irf3 deficiency in vitro but requires other variants in vivo.

Zika virus (ZIKV) is a Flavivirus responsible for recent epidemics in Pacific Islands and in the Americas. In humans, the consequences of ZIKV infection range from asymptomatic infection to severe neurological disease such as Guillain-Barré syndrome or fetal neurodevelopmental defects, suggesting, among other factors, the influence of host genetic variants. We previously reported similar diverse outcomes of ZIKV infection in mice of the Collaborative Cross (CC), a collection of inbred strains with large genetic diversity. CC071/TauUnc (CC071) was the most susceptible CC strain with severe symptoms and lethality. Notably, CC071 has been recently reported to be also susceptible to other flaviviruses including dengue virus, Powassan virus, West Nile virus, and to Rift Valley fever virus. To identify the genetic origin of this broad susceptibility, we investigated ZIKV replication in mouse embryonic fibroblasts (MEFs) from CC071 and two resistant strains. CC071 showed uncontrolled ZIKV replication associated with delayed induction of type-I interferons (IFN-I). Genetic analysis identified a mutation in the Irf3 gene specific to the CC071 strain which prevents the protein phosphorylation required to activate interferon beta transcription. We demonstrated that this mutation induces the same defective IFN-I response and uncontrolled viral replication in MEFs as an Irf3 knock-out allele. By contrast, we also showed that Irf3 deficiency did not induce the high plasma viral load and clinical severity observed in CC071 mice and that susceptibility alleles at other genes, not associated with the IFN-I response, are required. Our results provide new insight into the in vitro and in vivo roles of Irf3, and into the genetic complexity of host responses to flaviviruses.

Gene expression genetics of the striatum of Diversity Outbred mice.

Brain transcriptional variation is a heritable trait that mediates complex behaviors, including addiction. Expression quantitative trait locus (eQTL) mapping reveals genomic regions harboring genetic variants that influence transcript abundance. In this study, we profiled transcript abundance in the striatum of 386 Diversity Outbred (J:DO) mice of both sexes using RNA-Seq. All mice were characterized using a behavioral battery of widely-used exploratory and risk-taking assays prior to transcriptional profiling. We performed eQTL mapping, incorporated the results into a browser-based eQTL viewer, and deposited co-expression network members in GeneWeaver. The eQTL viewer allows researchers to query specific genes to obtain allelic effect plots, analyze SNP associations, assess gene expression correlations, and apply mediation analysis to evaluate whether the regulatory variant is acting through the expression of another gene. GeneWeaver allows multi-species comparison of gene sets using statistical and combinatorial tools. This data resource allows users to find genetic variants that regulate differentially expressed transcripts and place them in the context of other studies of striatal gene expression and function in addiction-related behavior.

Transcriptional signatures measured in whole blood correlate with protection against tuberculosis in inbred and outbred mice.

Although BCG has been used for almost 100 years to immunize against Mycobacterium tuberculosis, TB remains a global public health threat. Numerous clinical trials are underway studying novel vaccine candidates and strategies to improve or replace BCG, but vaccine development still lacks a well-defined set of immune correlates to predict vaccine-induced protection against tuberculosis. This study aimed to address this gap by examining transcriptional responses to BCG vaccination in C57BL/6 inbred mice, coupled with protection studies using Diversity Outbred mice. We evaluated relative gene expression in blood obtained from vaccinated mice, because blood is easily accessible, and data can be translated to human studies. We first determined that the average peak time after vaccination is 14 days for gene expression of a small subset of immune-related genes in inbred mice. We then performed global transcriptomic analyses using whole blood samples obtained two weeks after mice were vaccinated with BCG. Using comparative bioinformatic analyses and qRT-PCR validation, we developed a working correlate panel of 18 genes that were highly correlated with administration of BCG but not heat-killed BCG. We then tested this gene panel using BCG-vaccinated Diversity Outbred mice and revealed associations between the expression of a subset of genes and disease outcomes after aerosol challenge with M. tuberculosis. These data therefore demonstrate that blood-based transcriptional immune correlates measured within a few weeks after vaccination can be derived to predict protection against M. tuberculosis, even in outbred populations.

Lipidomic QTL in Diversity Outbred mice identifies a novel function for α/ß hydrolase domain 2 (Abhd2) as an enzyme that metabolizes phosphatidylcholine and cardiolipin.

We and others have previously shown that genetic association can be used to make causal connections between gene loci and small molecules measured by mass spectrometry in the bloodstream and in tissues. We identified a locus on mouse chromosome 7 where several phospholipids in liver showed strong genetic association to distinct gene loci. In this study, we integrated gene expression data with genetic association data to identify a single gene at the chromosome 7 locus as the driver of the phospholipid phenotypes. The gene encodes α/ß-hydrolase domain 2 (Abhd2), one of 23 members of the ABHD gene family. We validated this observation by measuring lipids in a mouse with a whole-body deletion of Abhd2. The Abhd2KO mice had a significant increase in liver levels of phosphatidylcholine and phosphatidylethanolamine. Unexpectedly, we also found a decrease in two key mitochondrial lipids, cardiolipin and phosphatidylglycerol, in male Abhd2KO mice. These data suggest that Abhd2 plays a role in the synthesis, turnover, or remodeling of liver phospholipids.

Single-Cell Transcriptomics of Bone Marrow Stromal Cells in Diversity Outbred Mice: A Model for Population-Level scRNA-Seq Studies.

Genome-wide association studies (GWASs) have advanced our understanding of the genetics of osteoporosis; however, the challenge has been converting associations to causal genes. Studies have utilized transcriptomics data to link disease-associated variants to genes, but few population transcriptomics data sets have been generated on bone at the single-cell level. To address this challenge, we profiled the transcriptomes of bone marrow-derived stromal cells (BMSCs) cultured under osteogenic conditions from five diversity outbred (DO) mice using single-cell RNA-seq (scRNA-seq). The goal of the study was to determine if BMSCs could serve as a model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells from large populations of mice to inform genetic studies. By enriching for mesenchymal lineage cells in vitro, coupled with pooling of multiple samples and downstream genotype deconvolution, we demonstrate the scalability of this model for population-level studies. We demonstrate that dissociation of BMSCs from a heavily mineralized matrix had little effect on viability or their transcriptomic signatures. Furthermore, we show that BMSCs cultured under osteogenic conditions are diverse and consist of cells with characteristics of mesenchymal progenitors, marrow adipogenic lineage precursors (MALPs), osteoblasts, osteocyte-like cells, and immune cells. Importantly, all cells were similar from a transcriptomic perspective to cells isolated in vivo. We employed scRNA-seq analytical tools to confirm the biological identity of profiled cell types. SCENIC was used to reconstruct gene regulatory networks (GRNs), and we observed that cell types show GRNs expected of osteogenic and pre-adipogenic lineage cells. Further, CELLECT analysis showed that osteoblasts, osteocyte-like cells, and MALPs captured a significant component of bone mineral density (BMD) heritability. Together, these data suggest that BMSCs cultured under osteogenic conditions coupled with scRNA-seq can be used as a scalable and biologically informative model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells in large populations. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

A study of the influence of genetic variance and sex on the density and thickness of the calvarial bone in collaborative cross mice.

BACKGROUND: Bone microarchitecture is affected by multiple genes, each having a small effect on the external appearance. It is thus challenging to characterize the genes and their specific effect on bone thickness and porosity. The purpose of this study was to assess the heritability and the genetic variation effect, as well as the sex effect on the calvarial bone thickness (Ca.Th) and calvarial porosity (%PoV) using the Collaborative Cross (CC) mouse population. METHODS: In the study we examined the parietal bones of 56 mice from 9 lines of CC mice. Morphometric parameters were evaluated using microcomputed tomography (µCT) and included Ca.Th and %PoV. We then evaluated heritability, genetic versus environmental variance and the sex effect for these parameters. RESULTS: Our morphometric analysis showed that Ca.Th and %PoV are both significantly different among the CC lines with a broad sense heritability of 0.78 and 0.90, respectively. The sex effect within the lines was significant in line IL111 and showed higher values of Ca.Th and %PoV in females compared to males. In line IL19 there was a borderline sex effect in Ca.Th in which males showed higher values than females. CONCLUSIONS: These results stress the complexity of sex and genotype interactions controlling Ca.Th and %PoV, as the skeletal sexual dimorphism was dependent on the genetic background. This study also shows that the CC population is a powerful tool for establishing the genetic effect on these traits.

Neuroinvasive Flavivirus Pathogenesis Is Restricted by Host Genetic Factors in Collaborative Cross Mice, Independently of Oas1b.

Powassan virus (POWV) is an emerging tick-borne flavivirus that causes neuroinvasive diseases, including encephalitis, meningitis, and paralysis. Similar to other neuroinvasive flaviviruses, such as West Nile virus (WNV) and Japanese encephalitis virus (JEV), POWV disease presentation is heterogeneous, and the factors influencing disease outcome are not fully understood. We used Collaborative Cross (CC) mice to assess the impact of host genetic factors on POWV pathogenesis. We infected a panel of Oas1b-null CC lines with POWV and observed a range of susceptibility, indicating that host factors other than the well-characterized flavivirus restriction factor Oas1b modulate POWV pathogenesis in CC mice. Among the Oas1b-null CC lines, we identified multiple highly susceptible lines (0% survival), including CC071 and CC015, and two resistant lines, CC045 and CC057 (>75% survival). The susceptibility phenotypes generally were concordant among neuroinvasive flaviviruses, although we did identify one line, CC006, that was specifically resistant to JEV, suggesting that both pan-flavivirus and virus-specific mechanisms contribute to susceptibility phenotypes in CC mice. We found that POWV replication was restricted in bone marrow-derived macrophages from CC045 and CC057 mice, suggesting that resistance could result from cell-intrinsic restriction of viral replication. Although serum viral loads at 2 days postinfection were equivalent between resistant and susceptible CC lines, clearance of POWV from the serum was significantly enhanced in CC045 mice. Furthermore, CC045 mice had significantly lower viral loads in the brain at 7 days postinfection than did CC071 mice, suggesting that reduced central nervous system (CNS) infection contributes to the resistant phenotype of CC045 mice. IMPORTANCE Neuroinvasive flaviviruses, such as WNV, JEV, and POWV, are transmitted to humans by mosquitoes or ticks and can cause neurologic diseases, such as encephalitis, meningitis, and paralysis, and they can result in death or long-term sequelae. Although potentially severe, neuroinvasive disease is a rare outcome of flavivirus infection. The factors that determine whether someone develops severe disease after a flavivirus infection are not fully understood, but host genetic differences in polymorphic antiviral response genes likely contribute to the outcome of infection. We evaluated a panel of genetically diverse mice and identified lines with distinct outcomes following infection with POWV. We found that resistance to POWV pathogenesis corresponded to reduced viral replication in macrophages, more rapid clearance of virus in peripheral tissues, and reduced viral infection in the brain. These susceptible and resistant mouse lines will provide a system for investigating the pathogenic mechanisms of POWV and identifying polymorphic host genes that contribute to resistance.

Integrative analysis of hepatic transcriptional profiles reveals genetic regulation of atherosclerosis in hyperlipidemic Diversity Outbred-F1 mice.

Atherogenesis is an insipidus but precipitating process leading to serious consequences of many cardiovascular diseases (CVD). Numerous genetic loci contributing to atherosclerosis have been identified in human genome-wide association studies, but these studies have limitations in the ability to control environmental factors and to decipher cause/effect relationships. To assess the power of hyperlipidemic Diversity Outbred (DO) mice in facilitating quantitative trait loci (QTL) analysis of complex traits, we generated a high-resolution genetic panel of atherosclerosis susceptible (DO-F1) mouse cohort by crossing 200 DO females with C57BL/6J males carrying two human genes: encoding apolipoprotein E3-Leiden and cholesterol ester transfer protein. We examined atherosclerotic traits including plasma lipids and glucose in the 235 female and 226 male progeny before and after 16 weeks of a high-fat/cholesterol diet, and aortic plaque size at 24 weeks. We also assessed the liver transcriptome using RNA-sequencing. Our QTL mapping for atherosclerotic traits identified one previously reported female-specific QTL on Chr10 with a narrower interval of 22.73 to 30.80 Mb, and one novel male-specific QTL at 31.89 to 40.25 Mb on Chr19. Liver transcription levels of several genes within each QTL were highly correlated with the atherogenic traits. A majority of these candidates have already known atherogenic potential in humans and/or mice, but integrative QTL, eQTL, and correlation analyses further pointed Ptprk as a major candidate of the Chr10 QTL, while Pten and Cyp2c67 of the Chr19 QTL in our DO-F1 cohort. Finally, through additional analyses of RNA-seq data we identified genetic regulation of hepatic transcription factors, including Nr1h3, contributes to atherogenesis in this cohort. Thus, an integrative approach using DO-F1 mice effectively validates the influence of genetic factors on atherosclerosis in DO mice and suggests an opportunity to discover therapeutics in the setting of hyperlipidemia.

Intravenous BCG Vaccination of Diversity Outbred Mice Results in Moderately Enhanced Protection against Challenge with Mycobacterium tuberculosis Compared to Intradermal Vaccination.

Tuberculosis is still the leading cause of death globally from any infectious disease, despite the widespread use of the live attenuated vaccine Bacille Calmette Guerin (BCG). While BCG has some efficacy against disseminated TB disease in children, protection wanes into adulthood resulting in over 1.8 million TB deaths per year. This has led to efforts to develop novel vaccine candidates that either replace or boost BCG, as well as to test novel delivery mechanisms to enhance BCG's efficacy. Traditional BCG vaccination is performed as an intradermal (ID) injection but delivering BCG by an alternate route may enhance the depth and breadth of protection. Previously, we demonstrated that phenotypically and genotypically disparate Diversity Outbred (DO) mice have heterogenous responses to M. tuberculosis challenge following intradermal BCG vaccination. Here, we utilize DO mice to examine BCG-induced protection when BCG is delivered systemically via intravenous (IV) administration. We find that DO mice vaccinated with IV BCG had a greater distribution of BCG throughout their organs compared to ID-vaccinated animals. However, compared to ID-vaccinated mice, M. tuberculosis burdens in lungs and spleens were not significantly reduced in animals vaccinated with BCG IV, nor was lung inflammation significantly altered. Nonetheless, DO mice that received BCG IV had increased survival over those vaccinated by the traditional ID route. Thus, our results suggest that delivering BCG by the alternate IV route enhances protection as detected in this diverse small animal model.

Refining risk estimates for lead in drinking water based on the impact of genetics and diet on blood lead levels using the Collaborative Cross mouse population.

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.

Unraveling the Host Genetic Background Effect on Internal Organ Weight Influenced by Obesity and Diabetes Using Collaborative Cross Mice.

Type 2 diabetes mellitus (T2DM) is a severe chronic epidemic that results from the body's improper usage of the hormone insulin. Globally, 700 million people are expected to have received a diabetes diagnosis by 2045, according to the International Diabetes Federation (IDF). Cancer and macro- and microvascular illnesses are only a few immediate and long-term issues it could lead to. T2DM accelerates the effect of organ weights by triggering a hyperinflammatory response in the body's organs, inhibiting tissue repair and resolving inflammation. Understanding how genetic variation translates into different clinical presentations may highlight the mechanisms through which dietary elements may initiate or accelerate inflammatory disease processes and suggest potential disease-prevention techniques. To address the host genetic background effect on the organ weight by utilizing the newly developed mouse model, the Collaborative Cross mice (CC). The study was conducted on 207 genetically different CC mice from 8 CC lines of both sexes. The experiment started with 8-week-old mice for 12 weeks. During this period, one group maintained a standard chow diet (CHD), while the other group maintained a high-fat diet (HFD). In addition, body weight was recorded bi-weekly, and at the end of the study, a glucose tolerance test, as well as tissue collection (liver, spleen, heart), were conducted. Our study observed a strong effect of HFD on blood glucose clearance among different CC lines. The HFD decreased the blood glucose clearance displayed by the significant Area Under Curve (AUC) values in both populations. In addition, variation in body weight changes among the different CC lines in response to HFD. The female liver weight significantly increased compared to males in the overall population when exposed to HFD. Moreover, males showed higher heritability values than females on the same diet. Regardless of the dietary challenge, the liver weight in the overall male population correlated positively with the final body weight. The liver weight results revealed that three different CC lines perform well under classification models. The regression results also varied among organs. Accordingly, the differences among these lines correspond to the genetic variance, and we suspect that some genetic factors invoke different body responses to HFD. Further investigations, such as quantitative trait loci (QTL) analysis and genomic studies, could find these genetic elements. These findings would prove critical factors for developing personalized medicine, as they could indicate future body responses to numerous situations early, thus preventing the development of complex diseases.

The microbial community dynamics of cocaine sensitization in two behaviorally divergent strains of collaborative cross mice.

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.

Thirdhand tobacco smoke exposure increases the genetic background-dependent risk of pan-tumor development in Collaborative Cross mice.

Increasing evidence has shown that thirdhand smoke (THS) exposure is likely to induce adverse health effects. An important knowledge gap remains in our understanding of THS exposure related to cancer risk in the human population. Population-based animal models are useful and powerful in investigating the interplay between host genetics and THS exposure on cancer risk. Here, we used the Collaborative Cross (CC) mouse population-based model system, which recapitulates the genetic and phenotypic diversity observed in the human population, to assess cancer risk after a short period of exposure, between 4 and 9 weeks of age. Eight CC strains (CC001, CC019, CC026, CC036, CC037, CC041, CC042 and CC051) were included in our study. We quantified pan-tumor incidence, tumor burden per mouse, organ tumor spectrum and tumor-free survival until 18 months of age. At the population level, we observed a significantly increased pan-tumor incidence and tumor burden per mouse in THS-treated mice as compared to the control (p = 3.04E-06). Lung and liver tissues exhibited the largest risk of undergoing tumorigenesis after THS exposure. Tumor-free survival was significantly reduced in THS-treated mice compared to control (p = 0.044). At the individual strain level, we observed a large variation in tumor incidence across the 8 CC strains. CC036 and CC041 exhibited a significant increase in pan-tumor incidence (p = 0.0084 and p = 0.000066, respectively) after THS exposure compared to control. We conclude that early-life THS exposure increases tumor development in CC mice and that host genetic background plays an important role in individual susceptibility to THS-induced tumorigenesis. Genetic background is an important factor that should be taken into account when determining human cancer risk of THS exposure.