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.

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 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.

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).

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.

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.

Which mouse multiparental population is right for your study? The Collaborative Cross inbred strains, their F1 hybrids, or the Diversity Outbred population

Multiparental populations (MPPs) encompass greater genetic diversity than traditional experimental crosses of two inbred strains, enabling broader surveys of genetic variation underlying complex traits. Two such mouse MPPs are the Collaborative Cross (CC) inbred panel and the Diversity Outbred (DO) population, which are descended from the same eight inbred strains. Additionally, the F1 intercrosses of CC strains (CC-RIX) have been used and enable study designs with replicate outbred mice. Genetic analyses commonly used by researchers to investigate complex traits in these populations include characterizing how heritable a trait is, i.e. its heritability, and mapping its underlying genetic loci, i.e. its quantitative trait loci (QTLs). Here we evaluate the relative merits of these populations for these tasks through simulation, as well as provide recommendations for performing the quantitative genetic analyses. We find that sample populations that include replicate animals, as possible with the CC and CC-RIX, provide more efficient and precise estimates of heritability. We report QTL mapping power curves for the CC, CC-RIX, and DO across a range of QTL effect sizes and polygenic backgrounds for samples of 174 and 500 mice. The utility of replicate animals in the CC and CC-RIX for mapping QTLs rapidly decreased as traits became more polygenic. Only large sample populations of 500 DO mice were well-powered to detect smaller effect loci (7.5-10%) for highly complex traits (80% polygenic background). All results were generated with our R package musppr, which we developed to simulate data from these MPPs and evaluate genetic analyses from user-provided genotypes.

Microbial glutamate metabolism predicts intravenous cocaine self-administration in diversity outbred mice

The gut microbiome is thought to play a critical role in the onset and development of psychiatric disorders, including depression and substance use disorder (SUD). To test the hypothesis that the microbiome affects addiction predisposing behaviors and cocaine intravenous self-administration (IVSA) and to identify specific microbes involved in the relationship, we performed 16S rRNA gene sequencing on feces from 228 diversity outbred mice. Twelve open field measures, two light-dark assay measures, one hole board and novelty place preference measure significantly differed between mice that acquired cocaine IVSA (ACQ) and those that failed to acquire IVSA (FACQ). We found that ACQ mice are more active and exploratory and display decreased fear than FACQ mice. The microbial abundances that differentiated ACQ from FACQ mice were an increased abundance of Barnesiella, Ruminococcus, and Robinsoniella and decreased Clostridium IV in ACQ mice. There was a sex-specific correlation between ACQ and microbial abundance, a reduced Lactobacillus abundance in ACQ male mice, and a decreased Blautia abundance in female ACQ mice. The abundance of Robinsoniella was correlated, and Clostridium IV inversely correlated with the number of doses of cocaine self-administered during acquisition. Functional analysis of the microbiome composition of a subset of mice suggested that gut-brain modules encoding glutamate metabolism genes are associated with the propensity to self-administer cocaine. These findings establish associations between the microbiome composition and glutamate metabolic potential and the ability to acquire cocaine IVSA thus indicating the potential translational impact of targeting the gut microbiome or microbial metabolites for treatment of SUD. This article is part of the Special Issue on "Microbiome & the Brain: Mechanisms & Maladies".

Effect of an obesogenic high-fat and high-sucrose diet on hepatic gene expression signatures in male Collaborative Cross mice

Nonalcoholic fatty liver disease (NAFLD), the most prevalent chronic liver disease, is characterized by substantial variations in case-level severity. In this study, we used a genetically diverse Collaborative Cross (CC) mouse population model to analyze the global transcriptome and clarify the molecular mechanisms involved in hepatic fat accumulation that determine the level and severity of NAFLD. Twenty-four strains of male CC mice were maintained on a high-fat/high-sucrose (HF/HS) diet for 12 wk, and their hepatic gene expression profiles were determined by next-generation RNA sequencing. We found that the development of the nonalcoholic fatty liver (NAFL) phenotype in CC mice coincided with significant changes in the expression of hepatic genes at the population level, evidenced by the presence of 724 differentially expressed genes involved in lipid and carbohydrate metabolism, cell morphology, vitamin and mineral metabolism, energy production, and DNA replication, recombination, and repair. Importantly, expression of 68 of these genes strongly correlated with the extent of hepatic lipid accumulation in the overall population of HF/HS diet-fed male CC mice. Results of partial least squares (PLS) modeling showed that these derived hepatic gene expression signatures help to identify the individual mouse strains that are highly susceptible to the development of NAFLD induced by an HF/HS diet. These findings imply that gene expression profiling, combined with a PLS modeling approach, may be a useful tool to predict NAFLD severity in genetically diverse patient populations.NEW & NOTEWORTHY Feeding male Collaborative Cross mice an obesogenic diet allows modeling NAFLD at the population level. The development of NAFLD coincided with significant hepatic transcriptomic changes in this model. Genes (724) were differentially expressed and expression of 68 genes strongly correlated with the extent of hepatic lipid accumulation. Partial least squares modeling showed that derived hepatic gene expression signatures may help to identify individual mouse strains that are highly susceptible to the development of NAFLD.

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

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

Intramembranous bone regeneration in diversity outbred mice is heritable

There are over one million cases of failed bone repair in the U.S. annually, resulting in substantial patient morbidity and societal costs. Multiple candidate genes affecting bone traits such as bone mineral density have been identified in human subjects and animal models using genome-wide association studies (GWAS). This approach for understanding the genetic factors affecting bone repair is impractical in human subjects but could be performed in a model organism if there is sufficient variability and heritability in the bone regeneration response. Diversity Outbred (DO) mice, which have significant genetic diversity and have been used to examine multiple intact bone traits, would be an excellent possibility. Thus, we sought to evaluate the phenotypic distribution of bone regeneration, sex effects and heritability of intramembranous bone regeneration on day 7 following femoral marrow ablation in 47 12-week old DO mice (23 males, 24 females). Compared to a previous study using 4 inbred mouse strains, we found similar levels of variability in the amount of regenerated bone (coefficient of variation of 86 % v. 88 %) with approximately the same degree of heritability (0.42 v. 0.49). There was a trend toward more bone regeneration in males than females. The amount of regenerated bone was either weakly or not correlated with bone mass at intact sites, suggesting that the genetic factors responsible for bone regeneration and intact bone phenotypes are at least partially independent. In conclusion, we demonstrate that DO mice exhibit variation and heritability of intramembranous bone regeneration that will be suitable for future GWAS.

Non-alcoholic fatty liver disease-associated DNA methylation and gene expression alterations in the livers of Collaborative Cross mice fed an obesogenic high-fat and high-sucrose diet

Non-alcoholic fatty liver disease (NAFLD) is a highly prevalent chronic liver disease, and patient susceptibility to its onset and progression is influenced by several factors. In this study, we investigated whether altered hepatic DNA methylation in liver tissue correlates with the degree of severity of NAFLD-like liver injury induced by a high-fat and high-sucrose (HF/HS) diet in Collaborative Cross (CC) mice. Using genome-wide targeted bisulphite DNA methylation next-generation sequencing, we found that mice with different non-alcoholic fatty liver (NAFL) phenotypes could be distinguished by changes in hepatic DNA methylation profiles. Specifically, NAFL-prone male CC042 mice exhibited more prominent DNA methylation changes compared with male CC011 mice and female CC011 and CC042 mice that developed only a mild NAFL phenotype. Moreover, these mouse strains demonstrated different patterns of DNA methylation. While the HF/HS diet induced both DNA hypomethylation and DNA hypermethylation changes in all the mouse strains, the NAFL-prone male CC042 mice demonstrated a global predominance of DNA hypermethylation, whereas a more pronounced DNA hypomethylation pattern developed in the mild-NAFL phenotypic mice. In a targeted analysis of selected genes that contain differentially methylated regions (DMRs), we identified NAFL phenotype-associated differences in DNA methylation and gene expression of the Apoa4, Gls2, and Apom genes in severe NAFL-prone mice but not in mice with mild NAFL phenotypes. These changes in the expression of Apoa4 and Gls2 coincided with similar findings in a human in vitro cell model of diet-induced steatosis and in patients with NAFL. These results suggest that changes in the expression and DNA methylation status of these three genes may serve as a set of predictive markers for the development of NAFLD.

Using the Collaborative Cross and Diversity Outbred Mice in Immunology

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

Genetic diversity of collaborative cross mice enables identification of novel rift valley fever virus encephalitis model

Rift Valley fever (RVF) is an arboviral disease of humans and livestock responsible for severe economic and human health impacts. In humans, RVF spans a variety of clinical manifestations, ranging from an acute flu-like illness to severe forms of disease, including late-onset encephalitis. The large variations in human RVF disease are inadequately represented by current murine models, which overwhelmingly die of early-onset hepatitis. Existing mouse models of RVF encephalitis are either immunosuppressed, display an inconsistent phenotype, or develop encephalitis only when challenged via intranasal or aerosol exposure. In this study, the genetically defined recombinant inbred mouse resource known as the Collaborative Cross (CC) was used to identify mice with additional RVF disease phenotypes when challenged via a peripheral foot-pad route to mimic mosquito-bite exposure. Wild-type Rift Valley fever virus (RVFV) challenge of 20 CC strains revealed three distinct disease phenotypes: early-onset hepatitis, mixed phenotype, and late-onset encephalitis. Strain CC057/Unc, with the most divergent phenotype, which died of late-onset encephalitis at a median of 11 days post-infection, is the first mouse strain to develop consistent encephalitis following peripheral challenge. CC057/Unc mice were directly compared to C57BL/6 mice, which uniformly succumb to hepatitis within 2–4 days of infection. Encephalitic disease in CC057/Unc mice was characterized by high viral RNA loads in brain tissue, accompanied by clearance of viral RNA from the periphery, low ALT levels, lymphopenia, and neutrophilia. In contrast, C57BL/6 mice succumbed from hepatitis at 3 days post-infection with high viral RNA loads in the liver, viremia, high ALT levels, lymphopenia, and thrombocytopenia. The identification of a strain of CC mice as an RVFV encephalitis model will allow for future investigation into the pathogenesis and treatment of RVF encephalitic disease and indicates that genetic background makes a major contribution to RVF disease variation.

Genetic background influences the effect of thirdhand smoke exposure on anxiety and memory in Collaborative Cross mice

Growing evidence indicates that thirdhand smoke (THS) exposure induces many adverse health effects. However, it is unclear how THS exposure affects behavior and how host genetic background modulates phenotypic changes. Here we used the Collaborative Cross (CC) mouse population-based model to assess behavioral alterations immediately after THS exposure from 4 to 9 weeks of age. We first measured anxiety-like behavior in six strains using light/dark box combined with a custom multivariate mouse tracking system. We developed an anxiety risk scoring system based on anxiety-related traits and then evaluated the THS impact on them. THS exposure significantly decreased anxiety risk in CC019 (P = 0.002) and CC051 (P = 0.009), but increased anxiety risk in CC036 (P < 0.001), while the other three strains did not show significant changes in anxiety-related traits. Such differences were driven by female mice for the six measures of anxiety-like behavior. Memory potential was measured in the same cohort of mice using the passive avoidance assay. Both THS-exposed male and female CC019 mice displayed significant memory loss compared to controls while no significant changes were found in the other five strains. This study provides strong evidence that THS exposure leads to strain-dependent changes in anxiety-like behavior and memory, suggesting that host genetic variations play a critical role in individual susceptibility to THS-induced effects.

Hepatic transcriptional profile reveals the role of diet and genetic backgrounds on metabolic traits in female progenitor strains of the Collaborative Cross

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.

Diet-induced obesity in genetically diverse collaborative cross mouse founder strains reveals diverse phenotype response and amelioration by quercetin treatment in 129S1/SvImJ, PWK/EiJ, CAST/PhJ, and WSB/EiJ mice

Significant evidence suggests protective effects of flavonoids against obesity in animal models, but these often do not translate to humans. One explanation for this disconnect is use of a few mouse strains (notably C57BL/6 J) in obesity studies. Obesity is a multifactorial disease. The underlying causes are not fully replicated by the high-fat C57BL/6 J model, despite phenotypic similarities. Furthermore, the impact of genetic factors on the activities of flavonoids is unknown. This study was designed to explore how diverse mouse strains respond to diet-induced obesity when fed a representative flavonoid. A subset of Collaborative Cross founder strains (males and females) were placed on dietary treatments (low-fat, high-fat, high-fat with quercetin, high-fat with quercetin and antibiotics) longitudinally. Diverse responses were observed across strains and sexes. Quercetin appeared to moderately blunt weight gain in male C57 and both sexes of 129S1/SvImJ mice, and slightly increased weight gain in female C57 mice. Surprisingly, quercetin dramatically blunted weight gain in male, but not female, PWK/PhJ mice. For female mice, quercetin blunted weight gain (relative to the high-fat phase) in CAST/PhJ, PWK/EiJ and WSB/EiJ mice compared to C57. Antibiotics did not generally result in loss of protective effects of quercetin. This highlights complex interactions between genetic factors, sex, obesity stimuli, and flavonoid intake, and the need to move away from single inbred mouse models to enhance translatability to diverse humans. These data justify use of genetically diverse Collaborative Cross and Diversity Outbred models which are emerging as invaluable tools in the field of personalized nutrition.

Facial shape and allometry quantitative trait locus intervals in the Diversity Outbred mouse are enriched for known skeletal and facial development genes

The biology of how faces are built and come to differ from one another is complex. Discovering normal variants that contribute to differences in facial morphology is one key to untangling this complexity, with important implications for medicine and evolutionary biology. This study maps quantitative trait loci (QTL) for skeletal facial shape using Diversity Outbred (DO) mice. The DO is a randomly outcrossed population with high heterozygosity that captures the allelic diversity of eight inbred mouse lines from three subspecies. The study uses a sample of 1147 DO animals (the largest sample yet employed for a shape QTL study in mouse), each characterized by 22 three-dimensional landmarks, 56,885 autosomal and X-chromosome markers, and sex and age classifiers. We identified 37 facial shape QTL across 20 shape principal components (PCs) using a mixed effects regression that accounts for kinship among observations. The QTL include some previously identified intervals as well as new regions that expand the list of potential targets for future experimental study. Three QTL characterized shape associations with size (allometry). Median support interval size was 3.5 Mb. Narrowing additional analysis to QTL for the five largest magnitude shape PCs, we found significant overrepresentation of genes with known roles in growth, skeletal and facial development, and sensory organ development. For most intervals, one or more of these genes lies within 0.25 Mb of the QTL's peak. QTL effect sizes were small, with none explaining more than 0.5% of facial shape variation. Thus, our results are consistent with a model of facial diversity that is influenced by key genes in skeletal and facial development and, simultaneously, is highly polygenic.

Whole Genome Sequencing and Progress Toward Full Inbreeding of the Mouse Collaborative Cross Population

Two key features of recombinant inbred panels are well-characterized genomes and reproducibility. Here we report on the sequenced genomes of six additional Collaborative Cross (CC) strains and on inbreeding progress of 72 CC strains. We have previously reported on the sequences of 69 CC strains that were publicly available, bringing the total of CC strains with whole genome sequence up to 75. The sequencing of these six CC strains updates the efforts toward inbreeding undertaken by the UNC Systems Genetics Core. The timing reflects our competing mandates to release to the public as many CC strains as possible while achieving an acceptable level of inbreeding. The new six strains have a higher than average founder contribution from non-domesticus strains than the previously released CC strains. Five of the six strains also have high residual heterozygosity (>14%), which may be related to non-domesticus founder contributions. Finally, we report on updated estimates on residual heterozygosity across the entire CC population using a novel, simple and cost effective genotyping platform on three mice from each strain. We observe a reduction in residual heterozygosity across all previously released CC strains. We discuss the optimal use of different genetic resources available for the CC population.