Common forms of atherosclerosis involve multiple genetic and environmental factors. While human genome-wide association studies have identified numerous loci contributing to coronary artery disease and its risk factors, these studies are unable to control environmental factors or examine detailed molecular traits in relevant tissues. We now report a study of natural variations contributing to atherosclerosis and related traits in over 100 inbred strains of mice from the Hybrid Mouse Diversity Panel (HMDP). The mice were made hyperlipidemic by transgenic expression of human apolipoprotein E-Leiden (APOE-Leiden) and human cholesteryl ester transfer protein (CETP). The mice were examined for lesion size and morphology as well as plasma lipid, insulin and glucose levels, and blood cell profiles. A subset of mice was studied for plasma levels of metabolites and cytokines. We also measured global transcript levels in aorta and liver. Finally, the uptake of acetylated LDL by macrophages from HMDP mice was quantitatively examined. Loci contributing to the traits were mapped using association analysis, and relationships among traits were examined using correlation and statistical modeling. A number of conclusions emerged. First, relationships among atherosclerosis and the risk factors in mice resemble those found in humans. Second, a number of trait-loci were identified, including some overlapping with previous human and mouse studies. Third, gene expression data enabled enrichment analysis of pathways contributing to atherosclerosis and prioritization of candidate genes at associated loci in both mice and humans. Fourth, the data provided a number of mechanistic inferences; for example, we detected no association between macrophage uptake of acetylated LDL and atherosclerosis. Fifth, broad sense heritability for atherosclerosis was much larger than narrow sense heritability, indicating an important role for gene-by-gene interactions. Sixth, stepwise linear regression showed that the combined variations in plasma metabolites, including LDL/VLDL-cholesterol, trimethylamine N-oxide (TMAO), arginine, glucose and insulin, account for approximately 30 to 40% of the variation in atherosclerotic lesion area. Overall, our data provide a rich resource for studies of complex interactions underlying atherosclerosis.

Malignant phyllodes tumor is a rare breast malignancy with sarcomatous overgrowth and with limited effective treatment options for recurrent and metastatic cases. Recent clinical trials indicated a potential for anti-angiogenic, anti-EGFR and immunotherapeutic approaches for patients with sarcomas, which led us to investigate these and other targetable pathways in malignant phyllodes tumor of the breast. Thirty-six malignant phyllodes tumors (including 8 metastatic tumors with two cases having matched primary and metastatic tumors) were profiled using gene sequencing, gene copy number analysis, whole genome expression, and protein expression. Whole genome expression analysis demonstrated consistent over-expression of genes involved in angiogenesis including VEGFA, Angiopoietin-2, VCAM1, PDGFRA, and PTTG1. EGFR protein overexpression was observed in 26/27 (96%) of cases with amplification of the EGFR gene in 8/24 (33%) cases. Two EGFR mutations were identified including EGFRvIII and a presumed pathogenic V774M mutation, respectively. The most common pathogenic mutations included TP53 (50%) and PIK3CA (15%). Cases with matched primary and metastatic tumors harbored identical mutations in both sites (PIK3CA/KRAS and RB1 gene mutations, respectively). Tumor expression of PD-L1 immunoregulatory protein was observed in 3/22 (14%) of cases. Overexpression of molecular biomarkers of increased angiogenesis, EGFR and immune checkpoints provides novel targeted therapy options in malignant phyllodes tumors of the breast.

Renal oncocytoma is a benign tumor with characteristic histologic findings. We describe an oncocytoma-like renal tumor with progression to high-grade oncocytic carcinoma and metastasis. A 74-year-old man with no family history of cancer presented with hematuria. Computed tomography showed an 11 cm heterogeneous multilobulated mass in the right kidney lower pole, enlarged aortocaval lymph nodes, and multiple lung nodules. In the nephrectomy specimen, approximately one third of the renal tumor histologically showed regions classic for benign oncocytoma transitioning to regions of high-grade carcinoma without sharp demarcation. With extensive genomic investigation using single nucleotide polymorphism-based array virtual karyotyping, multiregion sequencing, and expression array analysis, we were able to show a common lineage between the benign oncocytoma and high-grade oncocytic carcinoma regions in the tumor. We were also able to show karyotypic differences underlying this progression. The benign oncocytoma showed no chromosomal aberrations, whereas the high-grade oncocytic carcinoma showed loss of the 17p region housing FLCN (folliculin [Birt-Hogg-Dubé protein]), loss of 8p, and gain of 8q. Gene expression patterns supported dysregulation and activation of phosphoinositide 3-kinase (PI3K)/v-akt murine thymoma viral oncogene homolog (Akt), mitogen-activated protein kinase (MAPK)/extracellular-signal-regulated kinase (ERK), and mechanistic target of rapamycin (serine/threonine kinase) (mTOR) pathways in the high-grade oncocytic carcinoma regions. This was partly attributable to FLCN underexpression but further accentuated by overexpression of numerous genes on 8q. In the high-grade oncocytic carcinoma region, vascular endothelial growth factor A along with metalloproteinases matrix metallopeptidase 9 and matrix metallopeptidase 12 were overexpressed, facilitating angiogenesis and invasiveness. Genetic molecular testing provided evidence for the development of an aggressive oncocytic carcinoma from an oncocytoma, leading to aggressive targeted treatment but eventual death 39 months after the diagnosis.

We profiled and analyzed 283 metabolites representing eight major classes of molecules including Lipids, Carbohydrates, Amino Acids, Peptides, Xenobiotics, Vitamins and Cofactors, Energy Metabolism, and Nucleotides in mouse liver of 104 inbred and recombinant inbred strains. We find that metabolites exhibit a wide range of variation, as has been previously observed with metabolites in blood serum. Using genome-wide association analysis, we mapped 40% of the quantified metabolites to at least one locus in the genome and for 75% of the loci mapped we identified at least one candidate gene by local expression QTL analysis of the transcripts. Moreover, we validated 2 of 3 of the significant loci examined by adenoviral overexpression of the genes in mice. In our GWAS results, we find that at significant loci the peak markers explained on average between 20 and 40% of variation in the metabolites. Moreover, 39% of loci found to be regulating liver metabolites in mice were also found in human GWAS results for serum metabolites, providing support for similarity in genetic regulation of metabolites between mice and human. We also integrated the metabolomic data with transcriptomic and clinical phenotypic data to evaluate the extent of co-variation across various biological scales.

Metabolomics studies hold promise for the discovery of pathways linked to disease processes. Cardiovascular disease (CVD) represents the leading cause of death and morbidity worldwide. Here we used a metabolomics approach to generate unbiased small-molecule metabolic profiles in plasma that predict risk for CVD. Three metabolites of the dietary lipid phosphatidylcholine--choline, trimethylamine N-oxide (TMAO) and betaine--were identified and then shown to predict risk for CVD in an independent large clinical cohort. Dietary supplementation of mice with choline, TMAO or betaine promoted upregulation of multiple macrophage scavenger receptors linked to atherosclerosis, and supplementation with choline or TMAO promoted atherosclerosis. Studies using germ-free mice confirmed a critical role for dietary choline and gut flora in TMAO production, augmented macrophage cholesterol accumulation and foam cell formation. Suppression of intestinal microflora in atherosclerosis-prone mice inhibited dietary-choline-enhanced atherosclerosis. Genetic variations controlling expression of flavin monooxygenases, an enzymatic source of TMAO, segregated with atherosclerosis in hyperlipidaemic mice. Discovery of a relationship between gut-flora-dependent metabolism of dietary phosphatidylcholine and CVD pathogenesis provides opportunities for the development of new diagnostic tests and therapeutic approaches for atherosclerotic heart disease.

The Systems Genetics Resource (SGR) (http://systems.genetics.ucla.edu) is a new open-access web application and database that contains genotypes and clinical and intermediate phenotypes from both human and mouse studies. The mouse data include studies using crosses between specific inbred strains and studies using the Hybrid Mouse Diversity Panel. SGR is designed to assist researchers studying genes and pathways contributing to complex disease traits, including obesity, diabetes, atherosclerosis, heart failure, osteoporosis, and lipoprotein metabolism. Over the next few years, we hope to add data relevant to deafness, addiction, hepatic steatosis, toxin responses, and vascular injury. The intermediate phenotypes include expression array data for a variety of tissues and cultured cells, metabolite levels, and protein levels. Pre-computed tables of genetic loci controlling intermediate and clinical phenotypes, as well as phenotype correlations, are accessed via a user-friendly web interface. The web site includes detailed protocols for all of the studies. Data from published studies are freely available; unpublished studies have restricted access during their embargo period.

The osteoblast-lineage consists of cells at various stages of maturation that are essential for skeletal development, growth, and maintenance. Over the past decade, many of the signaling cascades that regulate this lineage have been elucidated; however, little is known of the networks that coordinate, modulate, and transmit these signals. Here, we identify a gene network specific to the osteoblast-lineage through the reconstruction of a bone co-expression network using microarray profiles collected on 96 Hybrid Mouse Diversity Panel (HMDP) inbred strains. Of the 21 modules that comprised the bone network, module 9 (M9) contained genes that were highly correlated with prototypical osteoblast maker genes and were more highly expressed in osteoblasts relative to other bone cells. In addition, the M9 contained many of the key genes that define the osteoblast-lineage, which together suggested that it was specific to this lineage. To use the M9 to identify novel osteoblast genes and highlight its biological relevance, we knocked-down the expression of its two most connected "hub" genes, Maged1 and Pard6g. Their perturbation altered both osteoblast proliferation and differentiation. Furthermore, we demonstrated the mice deficient in Maged1 had decreased bone mineral density (BMD). It was also discovered that a local expression quantitative trait locus (eQTL) regulating the Wnt signaling antagonist Sfrp1 was a key driver of the M9. We also show that the M9 is associated with BMD in the HMDP and is enriched for genes implicated in the regulation of human BMD through genome-wide association studies. In conclusion, we have identified a physiologically relevant gene network and used it to discover novel genes and regulatory mechanisms involved in the function of osteoblast-lineage cells. Our results highlight the power of harnessing natural genetic variation to generate co-expression networks that can be used to gain insight into the function of specific cell-types.

A segregating population of (C57BL/6J x DBA/2J)F2 intercross mice was studied for obesity-related traits and for global gene expression in liver. Quantitative trait locus analyses were applied to the subcutaneous fat-mass trait and all gene-expression data. These data were then used to identify gene sets that are differentially perturbed in lean and obese mice.

Systems-oriented genetic approaches that incorporate gene expression and genotype data are valuable in the quest for genetic regulatory loci underlying complex traits. Gene coexpression network analysis lends itself to identification of entire groups of differentially regulated genes-a highly relevant endeavor in finding the underpinnings of complex traits that are, by definition, polygenic in nature. Here we describe one such approach based on liver gene expression and genotype data from an F(2) mouse inter-cross utilizing weighted gene coexpression network analysis (WGCNA) of gene expression data to identify physiologically relevant modules. We describe two strategies: single-network analysis and differential network analysis. Single-network analysis reveals the presence of a physiologically interesting module that can be found in two distinct mouse crosses. Module quantitative trait loci (mQTLs) that perturb this module were discovered. In addition, we report a list of genetic drivers for this module. Differential network analysis reveals differences in connectivity and module structure between two networks based on the liver expression data of lean and obese mice. Functional annotation of these genes suggests a biological pathway involving epidermal growth factor (EGF). Our results demonstrate the utility of WGCNA in identifying genetic drivers and in finding genetic pathways represented by gene modules. These examples provide evidence that integration of network properties may well help chart the path across the gene-trait chasm.

Background The gut microbiota is modulated by a combination of diet, host genetics, and sex effects. The magnitude of these effects and interactions among them is important to understanding inter-individual variability in gut microbiota. In a previous study, mouse strain-specific responses to American and ketogenic diets were observed along with several QTL for metabolic traits. In the current study, we searched for genetic variants underlying differences in the gut microbiota in response to American and ketogenic diets, which are high in fat and vary in carbohydrate composition, between C57BL/6J (B6) and FVB/NJ (FVB) mouse strains. Results Genetic mapping of microbial features revealed 18 loci under the QTL model (i.e., marginal effects that are not specific to diet or sex), 12 loci under the QTL by diet model, and 1 locus under the QTL by sex model. Multiple metabolic and microbial features map to the distal part of Chr 1 and Chr 16 along with eigenvectors extracted from principal coordinate analysis of measures of β-diversity. Bilophila , Ruminiclostridium 9 , and Rikenella (Chr 1) were identified as sex and diet independent QTL candidate keystone organisms and Rikenelleceae RC9 Gut Group (Chr 16) was identified as a diet-specific, candidate keystone organism in confirmatory factor analyses of traits mapping to these regions. For many microbial features, irrespective of which QTL model was used, diet or the interaction between diet and a genotype were the strongest predictors of the abundance of each microbial trait. Sex, while important to the analyses, was not as strong of a predictor for microbial abundances. Conclusions These results demonstrate that sex, diet, and genetic background have different magnitudes of effects on inter-individual differences in gut microbiota. Therefore, Precision Nutrition through the integration of genetic variation, microbiota, and sex affecting microbiota variation will be important to predict response to diets varying in carbohydrate composition.

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.

Metabolic diseases are major public health issues worldwide and are responsible for disproportionately higher healthcare costs and increased complications of many diseases including SARS-CoV-2 infection. The Western Diet (WD) specifically is believed to be a major contributor to the global metabolic disease epidemic. In contrast, the Mediterranean diet (MeD), Ketogenic diet (KD), and Japanese diet (JD) are often considered beneficial for metabolic health. Yet, there is a growing appreciation that the effect of diet on metabolic health varies depending on several factors including host genetics. Additionally, poor metabolic health has also been attributed to altered gut microbial composition and/or function. To understand the complex relationship between host genetics, gut microbiota, and dietary patterns, we treated four widely used metabolically diverse inbred mouse strains (A/J, C57BL/6J, FVB/NJ, and NOD/ShiLtJ) with four human-relevant diets (MeD, JD, KD, WD), and a control mouse chow from 6 weeks to 30 weeks of age. We found that diet-induced alteration of gut microbiota (α-diversity, β-diversity, and abundance of several bacteria including Bifidobacterium, Ruminococcus, Turicibacter, Faecalibaculum, and Akkermansia) is significantly modified by host genetics. In addition, depending on the gut microbiota, the same diet could have different metabolic health effects. Our study also revealed that C57BL/6J mice are more susceptible to altered gut microbiota compared to other strains in this study indicating that host genetics is an important modulator of the diet-microbiota-metabolic health axis. Overall, our study demonstrated complex interactions between host genetics, gut microbiota, and diet on metabolic health; indicating the need to consider both host genetics and the gut microbiota in the development of new and more effective precision nutrition strategies to improve metabolic health.

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.

Our understanding of the genetic basis of learning and memory remains shrouded in mystery. To explore the genetic networks governing the biology of conditional fear, we used a systems genetics approach to analyze a hybrid mouse diversity panel (HMDP) with high mapping resolution.

The relationships between the levels of transcripts and the levels of the proteins they encode have not been examined comprehensively in mammals, although previous work in plants and yeast suggest a surprisingly modest correlation. We have examined this issue using a genetic approach in which natural variations were used to perturb both transcript levels and protein levels among inbred strains of mice. We quantified over 5,000 peptides and over 22,000 transcripts in livers of 97 inbred and recombinant inbred strains and focused on the 7,185 most heritable transcripts and 486 most reliable proteins. The transcript levels were quantified by microarray analysis in three replicates and the proteins were quantified by Liquid Chromatography-Mass Spectrometry using O(18)-reference-based isotope labeling approach. We show that the levels of transcripts and proteins correlate significantly for only about half of the genes tested, with an average correlation of 0.27, and the correlations of transcripts and proteins varied depending on the cellular location and biological function of the gene. We examined technical and biological factors that could contribute to the modest correlation. For example, differential splicing clearly affects the analyses for certain genes; but, based on deep sequencing, this does not substantially contribute to the overall estimate of the correlation. We also employed genome-wide association analyses to map loci controlling both transcript and protein levels. Surprisingly, little overlap was observed between the protein- and transcript-mapped loci. We have typed numerous clinically relevant traits among the strains, including adiposity, lipoprotein levels, and tissue parameters. Using correlation analysis, we found that a low number of clinical trait relationships are preserved between the protein and mRNA gene products and that the majority of such relationships are specific to either the protein levels or transcript levels. Surprisingly, transcript levels were more strongly correlated with clinical traits than protein levels. In light of the widespread use of high-throughput technologies in both clinical and basic research, the results presented have practical as well as basic implications.

The Hybrid Mouse Diversity Panel (HMDP) is a collection of approximately 100 well-characterized inbred strains of mice that can be used to analyze the genetic and environmental factors underlying complex traits. While not nearly as powerful for mapping genetic loci contributing to the traits as human genome-wide association studies, it has some important advantages. First, environmental factors can be controlled. Second, relevant tissues are accessible for global molecular phenotyping. Finally, because inbred strains are renewable, results from separate studies can be integrated. Thus far, the HMDP has been studied for traits relevant to obesity, diabetes, atherosclerosis, osteoporosis, heart failure, immune regulation, fatty liver disease, and host-gut microbiota interactions. High-throughput technologies have been used to examine the genomes, epigenomes, transcriptomes, proteomes, metabolomes, and microbiomes of the mice under various environmental conditions. All of the published data are available and can be readily used to formulate hypotheses about genes, pathways and interactions.

The genetic background of apolipoprotein E (apoE) deficient mice influences atherosclerotic plaque development. We previously reported three quantitative trait loci (QTL), Aath1-Aath3, that affect aortic arch atherosclerosis independently of those in the aortic root in a cross between C57BL6 apoEKO mice (B6-apoE) and 129S6 apoEKO mice (129-apoE). To gain further insight into genetic factors that influence atherosclerosis at different vascular locations, we analyzed 335 F2 mice from an intercross between 129-apoE and apoEKO mice on a DBA/2J genetic background (DBA-apoE). The extent of atherosclerosis in the aortic arch was very similar in the two parental strains. Nevertheless, a genome-wide scan identified two significant QTL for plaque size in the aortic arch: Aath4 on Chromosome (Chr) 2 at 137 Mb and Aath5 on Chr 10 at 51 Mb. The DBA alleles of Aath4 and Aath5 respectively confer susceptibility and resistance to aortic arch atherosclerosis over 129 alleles. Both QTL are also independent of those affecting plaque size at the aortic root. Genome analysis suggests that athero-susceptibility of Aath4 in DBA may be contributed by multiple genes, including Mertk and Cd93, that play roles in phagocytosis of apoptotic cells and modulate inflammation. A candidate gene for Aath5 is Stab2, the DBA allele of which is associated with 10 times higher plasma hyaluronan than the 129 allele. Overall, our identification of two new QTL that affect atherosclerosis in an aortic arch-specific manner further supports the involvement of distinct pathological processes at different vascular locations.

Atherosclerosis, the underlying cause of cardiovascular disease, results from both genetic and environmental factors.

Heritable epigenetic factors can contribute to complex disease etiology. Here we examine the contribution of DNA methylation to complex traits that are precursors to heart disease, diabetes, and osteoporosis. We profiled DNA methylation in the liver using bisulfite sequencing in 90 mouse inbred strains, genome-wide expression levels, proteomics, metabolomics, and 68 clinical traits and performed epigenome-wide association studies (EWAS). We found associations with numerous clinical traits including bone density, insulin resistance, expression, and protein and metabolite levels. A large proportion of associations were unique to EWAS and were not identified using GWAS. Methylation levels were regulated by genetics largely in cis, but we also found evidence of trans regulation, and we demonstrate that genetic variation in the methionine synthase reductase gene Mtrr affects methylation of hundreds of CpGs throughout the genome. Our results indicate that natural variation in methylation levels contributes to the etiology of complex clinical traits.

Obesity is a highly heritable disease driven by complex interactions between genetic and environmental factors. Human genome-wide association studies (GWAS) have identified a number of loci contributing to obesity; however, a major limitation of these studies is the inability to assess environmental interactions common to obesity. Using a systems genetics approach, we measured obesity traits, global gene expression, and gut microbiota composition in response to a high-fat/high-sucrose (HF/HS) diet of more than 100 inbred strains of mice. Here we show that HF/HS feeding promotes robust, strain-specific changes in obesity that are not accounted for by food intake and provide evidence for a genetically determined set point for obesity. GWAS analysis identified 11 genome-wide significant loci associated with obesity traits, several of which overlap with loci identified in human studies. We also show strong relationships between genotype and gut microbiota plasticity during HF/HS feeding and identify gut microbial phylotypes associated with obesity.