Altered macronutrient composition and genetics influence the complex transcriptional network associated with adiposity in the Collaborative Cross.

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

Obesogenic and diabetic effects of CD44 in mice are sexually dimorphic and dependent on genetic background.

INTRODUCTION: CD44 is a candidate gene for obesity and diabetes development and may be a critical mediator of a systemic inflammation associated with obesity and diabetes. METHODS: We investigated the relationship of CD44 with obesity in CD44-deficient mice challenged with a high-fat diet. RESULTS: In mice fed a diet high in fat, cholesterol, and sucrose for 12 weeks fat mass accumulation was reduced in CD44-deficient mice bred onto both a C57BL/6J and the naturally TLR deficient C3H/HeJ background. Reduced fat mass could not be attributed to lower food intake or an increase in energy expenditure as measured by indirect calorimetry. However, we observed a 40-60% lower mRNA expression of the inflammation markers, F4/80, CD11b, TNF-α, and CD14, in adipose tissue of CD44-deficient mice on the C57BL/6J background but not the C3H/HeJ background, perhaps indicating that alternative factors may be affecting adiposity in this model. Measures of hepatic steatosis and insulin sensitivity were improved in CD44-deficient mice on a C57BL/6J but not in the C3H/HeJ mice. These results were highly sexually dimorphic as there were no detectable effects of CD44 inactivation in female mice on a C57BL/6 J or C3H/HeJ background. CONCLUSION: CD44 was associated with adiposity, liver fat, and glucose in male mice. However, the effects of CD44 on obesity may be independent of TLR4 signaling.

Genetic architecture modulates diet-induced hepatic mRNA and miRNA expression profiles in Diversity Outbred mice.

Genetic approaches in model organisms have consistently demonstrated that molecular traits such as gene expression are under genetic regulation, similar to clinical traits. The resulting expression quantitative trait loci (eQTL) have revolutionized our understanding of genetic regulation and identified numerous candidate genes for clinically relevant traits. More recently, these analyses have been extended to other molecular traits such as protein abundance, metabolite levels, and miRNA expression. Here, we performed global hepatic eQTL and microRNA expression quantitative trait loci (mirQTL) analysis in a population of Diversity Outbred mice fed two different diets. We identified several key features of eQTL and mirQTL, namely differences in the mode of genetic regulation (cis or trans) between mRNA and miRNA. Approximately 50% of mirQTL are regulated by a trans-acting factor, compared to ∼25% of eQTL. We note differences in the heritability of mRNA and miRNA expression and variance explained by each eQTL or mirQTL. In general, cis-acting variants affecting mRNA or miRNA expression explain more phenotypic variance than trans-acting variants. Finally, we investigated the effect of diet on the genetic architecture of eQTL and mirQTL, highlighting the critical effects of environment on both eQTL and mirQTL. Overall, these data underscore the complex genetic regulation of two well-characterized RNA classes (mRNA and miRNA) that have critical roles in the regulation of clinical traits and disease susceptibility.

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.

Genetic Architecture Modulates Diet-Induced Hepatic mRNA and miRNA Expression Profiles in Diversity Outbred Mice.

Genetic approaches in model organisms have consistently demonstrated that molecular traits such as gene expression are under genetic regulation, similar to clinical traits. The resulting expression quantitative trait loci (eQTL) have revolutionized our understanding of genetic regulation and identified numerous candidate genes for clinically relevant traits. More recently, these analyses have been extended to other molecular traits such as protein abundance, metabolite levels, and miRNA expression. Here, we performed global hepatic eQTL and microRNA expression quantitative trait loci (mirQTL) analysis in a population of Diversity Outbred mice fed two different diets. We identified several key features of eQTL and mirQTL, namely differences in the mode of genetic regulation (cis or trans) between mRNA and miRNA. Approximately 50% of mirQTL are regulated by a trans-acting factor, compared to ∼25% of eQTL. We note differences in the heritability of mRNA and miRNA expression and variance explained by each eQTL or mirQTL. In general, cis-acting variants affecting mRNA or miRNA expression explain more phenotypic variance than trans-acting variants. Lastly, we investigated the effect of diet on the genetic architecture of eQTL and mirQTL, highlighting the critical effects of environment on both eQTL and mirQTL. Overall, these data underscore the complex genetic regulation of two well-characterized RNA classes (mRNA and miRNA) that have critical roles in the regulation of clinical traits and disease susceptibility.

Dissecting the Genetic Architecture of Cystatin C in Diversity Outbred Mice.

Plasma concentration of Cystatin C (CysC) level is a biomarker of glomerular filtration rate in the kidney. We use a Systems Genetics approach to investigate the genetic determinants of plasma CysC concentration. To do so we perform Quantitative Trait Loci (QTL) and expression QTL (eQTL) analysis of 120 Diversity Outbred (DO) female mice, 56 weeks of age. We performed network analysis of kidney gene expression to determine if the gene modules with common functions are associated with kidney biomarkers of chronic kidney diseases. Our data demonstrates that plasma concentrations and kidney mRNA levels of CysC are associated with genetic variation and are transcriptionally coregulated by immune genes. Specifically, Type-I interferon signaling genes are coexpressed with Cst3 mRNA levels and associated with CysC concentrations in plasma. Our findings demonstrate the complex control of CysC by genetic polymorphisms and inflammatory pathways.

Genetic Background Shapes Phenotypic Response to Diet for Adiposity in the Collaborative Cross.

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

microRNA-146a-5p association with the cardiometabolic disease risk factor TMAO.

Trimethylamine-N-oxide (TMAO), a microbial choline metabolism byproduct that is processed in the liver and excreted into circulation, is associated with increased atherosclerotic lesion formation and cardiovascular disease risk. Genetic regulators of TMAO levels are largely unknown. In the present study, we used 288 mice from a genetically heterogeneous mouse population [Diversity Outbred (DO)] to determine hepatic microRNA associations with TMAO in the context of an atherogenic diet. We also validated findings in two additional animal models of atherosclerosis: liver-specific insulin receptor knockout mice fed a chow diet (LIRKO) and African green monkeys fed high-fat/high-cholesterol diet. Small RNA-sequencing analysis in DO mice, LIRKO mice, and African green monkeys identified only one hepatic microRNA (miR-146a-5p) that is aberrantly expressed across all three models. Moreover, miR-146a-5p levels are associated with circulating TMAO after atherogenic diet in each of these models. We also performed high-resolution genetic mapping and identified a novel quantitative trait locus on Chromosome 12 for TMAO levels. This interval includes two genes, Numb and Dlst, which are inversely correlated with both miR-146a and TMAO and are predicted targets of miR-146a. Both of these genes have been validated as direct targets of miR-146a, though in other cellular contexts. This is the first report to our knowledge of a link between miR-146 and TMAO. Our findings suggest that miR-146-5p, as well as one or more genes at the Chromosome 12 QTL (possibly Numb or Dlst), is strongly linked to TMAO levels and likely involved in the control of atherosclerosis.

Leptin Signaling Mediates Obesity-Associated CSC Enrichment and EMT in Preclinical TNBC Models.

Obesity is associated with poor prognosis in triple-negative breast cancer (TNBC). Preclinical models of TNBC were used to test the hypothesis that increased leptin signaling drives obesity-associated TNBC development by promoting cancer stem cell (CSC) enrichment and/or epithelial-to-mesenchymal transition (EMT). MMTV-Wnt-1 transgenic mice, which develop spontaneous basal-like, triple-negative mammary tumors, received either a control diet (10% kcal from fat) or a diet-induced obesity regimen (DIO, 60% kcal from fat) for up to 42 weeks (n = 15/group). Mice were monitored for tumor development and euthanized when tumor diameter reached 1.5 cm. Tumoral gene expression was assessed via RNA sequencing (RNA-seq). DIO mice had greater body weight and percent body fat at termination than controls. DIO mice, versus controls, demonstrated reduced survival, increased systemic metabolic and inflammatory perturbations, upregulated tumoral CSC/EMT gene signature, elevated tumoral aldehyde dehydrogenase activity (a CSC marker), and greater leptin signaling. In cell culture experiments using TNBC cells (murine: E-Wnt and M-Wnt; human: MDA-MB-231), leptin enhanced mammosphere formation, and media supplemented with serum from DIO versus control mice increased cell viability, migration, invasion, and CSC- and EMT-related gene expression, including Foxc2, Twist2, Vim, Akt3, and Sox2 In E-Wnt cells, knockdown of leptin receptor ablated these procancer effects induced by DIO mouse serum. These findings indicate that increased leptin signaling is causally linked to obesity-associated TNBC development by promoting CSC enrichment and EMT.Implications: Leptin-associated signals impacting CSC and EMT may provide new targets and intervention strategies for decreasing TNBC burden in obese women. Mol Cancer Res; 16(5); 869-79. ©2018 AACR.

Systems genetics identifies a co-regulated module of liver microRNAs associated with plasma LDL cholesterol in murine diet-induced dyslipidemia.

Chronically altered levels of circulating lipids, termed dyslipidemia, is a significant risk factor for a number of metabolic and cardiovascular morbidities. MicroRNAs (miRNAs) have emerged as important regulators of lipid balance, have been implicated in dyslipidemia, and have been proposed as candidate therapeutic targets in lipid-related disorders including atherosclerosis. A major limitation of most murine studies of miRNAs in lipid metabolic disorders is that they have been performed in just one (or very few) inbred strains, such as C57BL/6. Moreover, although individual miRNAs have been associated with lipid phenotypes, it is well understood that miRNAs likely work together in functional modules. To address these limitations, we implemented a systems genetics strategy using the Diversity Outbred (DO) mouse population. Specifically, we performed gene and miRNA expression profiling in the livers from ~300 genetically distinct DO mice after 18 wk on either a high-fat/high-cholesterol diet or a high-protein diet. Large-scale correlative analysis of these data with a wide range of cardio-metabolic end points revealed a co-regulated module of miRNAs significantly associated with circulating low-density lipoprotein cholesterol (LDL-C) levels. The hubs of this module were identified as miR-199a, miR-181b, miR-27a, miR-21_-_1, and miR-24. In sum, we demonstrate that a high-fat/high-cholesterol diet robustly rewires the miRNA regulatory network, and we identify a small group of co-regulated miRNAs that may exert coordinated effects to control circulating LDL-C.

Genetic network identifies novel pathways contributing to atherosclerosis susceptibility in the innominate artery.

BACKGROUND: Atherosclerosis, the underlying cause of cardiovascular disease, results from both genetic and environmental factors. METHODS: In the current study we take a systems-based approach using weighted gene co-expression analysis to identify a candidate pathway of genes related to atherosclerosis. Bioinformatic analyses are performed to identify candidate genes and interactions and several novel genes are characterized using in-vitro studies. RESULTS: We identify 1 coexpression module associated with innominate artery atherosclerosis that is also enriched for inflammatory and macrophage gene signatures. Using a series of bioinformatics analysis, we further prioritize the genes in this pathway and identify Cd44 as a critical mediator of the atherosclerosis. We validate our predictions generated by the network analysis using Cd44 knockout mice. CONCLUSION: These results indicate that alterations in Cd44 expression mediate inflammation through a complex transcriptional network involving a number of previously uncharacterized genes.

Responsiveness of cardiometabolic-related microbiota to diet is influenced by host genetics.

Intestinal microbial community structure is driven by host genetics in addition to environmental factors such as diet. In comparison with environmental influences, the effect of host genetics on intestinal microbiota, and how host-driven differences alter host metabolism is unclear. Additionally, the interaction between host genetics and diet, and the impact on the intestinal microbiome and possible down-stream effect on host metabolism is not fully understood, but represents another aspects of inter-individual variation in disease risk. The objectives of this study were to investigate how diet and genetic background shape microbial communities, and how these diet- and genetic-driven microbial differences relate to cardiometabolic phenotypes. To determine these effects, we used the 8 progenitor strains of the collaborative cross/diversity outbred mapping panels (C57BL/6J, A/J, NOD/ShiLtJ, NZO/HILtJ, WSB/EiJ, CAST/EiJ, PWK/PhJ, and 129S1/SvImJ). 16s rRNA profiling of enteric microbial communities in addition to the assessment of phenotypes central to cardiometabolic health was conducted under baseline nutritional conditions and in response to diets varying in atherogenic nutrient (fat, cholesterol, cholic acid) composition. These studies revealed strain-driven differences in enteric microbial communities which were retained with dietary intervention. Diet-strain interactions were seen for a core group of cardiometabolic-related microbial taxa. In conclusion, these studies highlight diet and genetically regulated cardiometabolic-related microbial taxa. Furthermore, we demonstrate the progenitor model is useful for nutrigenomic-based studies and screens seeking to investigate the interaction between genetic background and the phenotypic and microbial response to diet.