Search

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.

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.

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

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

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

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

1,3-Butadiene (BD) exposure is known to cause numerous adverse health effects, including cancer, in animals and humans. BD is metabolized to reactive epoxide intermediates, which are genotoxic, but it is not well know what other effects BD has on cellular metabolism. We examined the effects of exposure to BD on the mouse lung metabolome in the genetically heterogeneous collaborative cross outbred mouse model. Mice were exposed to 3 concentra-tions of BD for 10 days (2, 20, and 200 ppm), and lung tissues were analyzed using high-resolution mass spectrometry-based metabolomics. As compared to controls (0 ppm BD), BD had extensive effects on lung metabolism at all concentrations of exposure, including the lowest concentration of 2 ppm, as reflected by reprogramming of multiple metabolic pathways. Metabolites participating in glycolysis and the tricarboxylic acid cycle were elevated, with 8 out of 10 metabolites demonstrating a 2 to 8-fold increase, including the oncometabolite fumarate. Fatty acid levels, sphingosine, and sphinganine were decreased (2 to 8-fold), and fatty acyl-CoAs were significantly increased (16 to 31-fold), suggesting adjustments in lipid metabolism. Furthermore, metabolites involved in basic amino acid metabolism, steroid hormone metabolism, and nucleic acid metabolism were significantly altered. Overall, these changes mirror the metabolic alterations found in lung cancer cells, suggesting that very low doses of BD induce metabolic adaptations that may prevent or promote adverse health effects such as tumor formation., (Published by Elsevier B.V.)

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

Arsenic methyltransferase (AS3MT) is the key enzyme in the pathway for the methylation of inorganic arsenic (iAs), a potent human carcinogen and diabetogen. AS3MT converts iAs to mono- and dimethylated arsenic species (MAs, DMAs) that are excreted mainly in urine. Polymorphisms in AS3MT is a key genetic factor affecting iAs metabolism and toxicity. The present study examined the role of As3mt polymorphisms in the susceptibility to the diabetogenic effects of iAs exposure using two Collaborative Cross mouse strains, CC021/Unc and CC027/GeniUnc, carrying different As3mt haplotypes. Male mice from the two strains were exposed to iAs in drinking water (0, 0.1 or 50 ppm) for 11 weeks. Blood glucose and plasma insulin levels were measured after 6-h fasting and 15 min after i.p. injection of glucose. Body composition was determined using magnetic resonance imaging. To asses iAs metabolism, the concentrations of iAs, MAs and DMAs were measured in urine. The results show that CC021 mice, both iAs-exposed and controls, had higher body fat percentage, lower fasting blood glucose, higher fasting plasma insulin, and were more insulin resistant than their CC027 counterparts. iAs exposure had a minor effect on diabetes indicators and only in CC027 mice. Blood glucose levels 15 min after glucose injection were significantly higher in CC027 mice exposed to 0.1 ppm iAs than in control mice. No significant differences were found in the concentrations or proportions of arsenic species in urine of CC021 and CC027 mice at the same exposure level. These results suggest that the differences in As3mt haplotypes did not affect the profiles of iAs or its metabolites in mouse urine. The major differences in diabetes indicators were associated with the genetic backgrounds of CC021 and CC027 mice. The effects of iAs exposure, while minor, were genotype- and dose-dependent., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Endocrine disrupting chemicals (EDCs) are compounds that alter the structure and function of the endocrine system and may be contributing to disorders of the reproductive, metabolic, neuroendocrine and other complex systems. Typically, these outcomes cannot be modeled in cell-based or other simple systems necessitating the use of animal testing. Appropriate animal model selection is required to effectively recapitulate the human experience, including relevant dosing and windows of exposure, and ensure translational utility and reproducibility. While classical toxicology heavily relies on inbred rats and mice, and focuses on apical endpoints such as tumor formation or birth defects, EDC researchers have used a greater diversity of species to effectively model more subtle but significant outcomes such as changes in pubertal timing, mammary gland development, and social behaviors. Advances in genomics, neuroimaging and other tools are making a wider range of animal models more widely available to EDC researchers., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Cytochrome P450 2E1 (CYP2E1) metabolizes low molecular weight hydrophobic compounds, including 1,3-butadiene, which is converted by CYP2E1 to electrophilic epoxide metabolites that covalently modify cellular proteins and DNA. Previous CYP2E1 studies have mainly focused on the enzyme localized in the endoplasmic reticulum (erCYP2E1); however, active CYP2E1 has also been found in mitochondria (mtCYP2E1) and the distribution of CYP2E1 between organelles can influence an individual's response to exposure. Relatively few studies have focused on the contribution of mtCYP2E1 to activation of chemical toxicants. We hypothesized that CYP2E1 bioactivation of 1,3-butadiene within mitochondria adversely affects mitochondrial respiratory complexes I-IV. A population of Collaborative Cross mice was exposed to air (control) or 200ppm 1,3-butadiene. Subcellular fractions (mitochondria, DNA, and microsomes) were collected from frozen livers and CYP2E1 activity was measured in microsomes and mitochondria. Individual activities of mitochondrial respiratory complexes I-IV were measured using in vitro assays and purified mitochondrial fractions. In air- and 1,3-butadiene-exposed mouse samples, mtDNA copy numbers were assessed by RT-PCR, and mtDNA integrity was assessed through a PCR-based assay. No significant changes in mtDNA copy number or integrity were observed; however, there was a decrease in overall activity of mitochondrial respiratory complexes I, II, and IV after 1,3-butadiene exposure. Additionally, higher mtCYP2E1 (but not erCYP2E1) activity was correlated with decreased mitochondrial respiratory complex activity (in complexes I-IV) in the 1,3-butadiene-exposed (not control) animals. Together, these results represent the first in vivo link between mitochondrial CYP2E1 activity and mitochondrial toxicity., (Copyright © 2017 Elsevier B.V. All rights reserved.)