Genome-wide association studies (GWASs) have identified more than 150 loci associated with blood lipid and cholesterol levels; however, the functional and molecular mechanisms for many associations are unknown. We examined the functional regulatory effects of candidate variants at the GALNT2 locus associated with high-density lipoprotein cholesterol (HDL-C). Fine-mapping and conditional analyses in the METSIM study identified a single locus harboring 25 noncoding variants (r(2) > 0.7 with the lead GWAS variants) strongly associated with total cholesterol in medium-sized HDL (e.g., rs17315646, p = 3.5 × 10(-12)). We used luciferase reporter assays in HepG2 cells to test all 25 variants for allelic differences in regulatory enhancer activity. rs2281721 showed allelic differences in transcriptional activity (75-fold [T] versus 27-fold [C] more than the empty-vector control), as did a separate 780-bp segment containing rs4846913, rs2144300, and rs6143660 (49-fold [AT(-) haplotype] versus 16-fold [CC(+) haplotype] more). Using electrophoretic mobility shift assays, we observed differential CEBPB binding to rs4846913, and we confirmed this binding in a native chromatin context by performing chromatin-immunoprecipitation (ChIP) assays in HepG2 and Huh-7 cell lines of differing genotypes. Additionally, sequence reads in HepG2 DNase-I-hypersensitivity and CEBPB ChIP-seq signals spanning rs4846913 showed significant allelic imbalance. Allelic-expression-imbalance assays performed with RNA from primary human hepatocyte samples and expression-quantitative-trait-locus (eQTL) data in human subcutaneous adipose tissue samples confirmed that alleles associated with increased HDL-C are associated with a modest increase in GALNT2 expression. Together, these data suggest that at least rs4846913 and rs2281721 play key roles in influencing GALNT2 expression at this HDL-C locus.

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.

Interindividual variation in the composition of the human gut microbiome was examined in relation to demographic and anthropometric traits, and to changes in dietary saturated fat intake and protein source. One hundred nine healthy men and women aged 21 to 65, with BMIs of 18 to 36, were randomized, after a two-week baseline diet, to high (15% total energy [E])- or low (7%E)-saturated-fat groups and randomly received three diets (four weeks each) in which the protein source (25%E) was mainly red meat (beef, pork) (12%E), white meat (chicken, turkey) (12%E), and nonmeat sources (nuts, beans, soy) (16%E). Taxonomic characterization using 16S ribosomal DNA was performed on fecal samples collected at each diet completion. Interindividual differences in age, body fat (%), height, ethnicity, sex, and alpha diversity (Shannon) were all significant factors, and most samples clustered by participant in the PCoA ordination. The dietary interventions did not significantly alter the overall microbiome community in ordination space, but there was an effect on taxon abundance levels. Saturated fat had a greater effect than protein source on taxon differential abundance, but protein source had a significant effect once the fat influence was removed. Higher alpha diversity predicted lower beta diversity between the experimental and baseline diets, indicating greater resistance to change in people with higher microbiome diversity. Our results suggest that interindividual differences outweighed the influence of these specific dietary changes on the microbiome and that moderate changes in saturated fat level and protein source correspond to modest changes in the microbiome.IMPORTANCE The microbiome has proven to influence health and disease, but how combinations of external factors affect the microbiome is relatively unknown. Diet can cause changes, but this is usually achieved by altering macronutrient ratios and has not focused on dietary protein source or saturated fat intake levels. In addition, each individual's unique microbiome profile can be an important factor during studies, and it has even been shown to affect therapeutic outcomes. We show here that the effects of individual differences outweighed the effect of experimental diets and that protein source is less influential than saturated fat level. This suggests that fat and protein composition, separate from macronutrient ratio and carbohydrate composition, is an important consideration in dietary studies.

HDLs are nanoparticles with more than 80 associated proteins, phospholipids, cholesterol, and cholesteryl esters. The potential inverse relation of HDL to coronary artery disease (CAD) and the effects of HDL on myriad other inflammatory conditions warrant a better understanding of the genetic basis of the HDL proteome. We conducted a comprehensive genetic analysis of the regulation of the proteome of HDL isolated from a panel of 100 diverse inbred strains of mice (the hybrid mouse diversity panel) and examined protein composition and efflux capacity to identify novel factors that affect the HDL proteome. Genetic analysis revealed widely varied HDL protein levels across the strains. Some of this variation was explained by local cis-acting regulation, termed cis-protein quantitative trait loci (QTLs). Variations in apoA-II and apoC-3 affected the abundance of multiple HDL proteins, indicating a coordinated regulation. We identified modules of covarying proteins and defined a protein-protein interaction network that describes the protein composition of the naturally occurring subspecies of HDL in mice. Sterol efflux capacity varied up to 3-fold across the strains, and HDL proteins displayed distinct correlation patterns with macrophage and ABCA1-specific cholesterol efflux capacity and cholesterol exchange, suggesting that subspecies of HDL participate in discrete functions. The baseline and stimulated sterol efflux capacity phenotypes were associated with distinct QTLs with smaller effect size, suggesting a multigenetic regulation. Our results highlight the complexity of HDL particles by revealing the high degree of heterogeneity and intercorrelation, some of which is associated with functional variation, and support the concept that HDL-cholesterol alone is not an accurate measure of HDL's properties, such as protection against CAD.

We report the genetic analysis of a "humanized" hyperlipidemic mouse model for progressive nonalcoholic steatohepatitis (NASH) and fibrosis. Mice carrying transgenes for human apolipoprotein E*3-Leiden and cholesteryl ester transfer protein and fed a "Western" diet were studied on the genetic backgrounds of over 100 inbred mouse strains. The mice developed hepatic inflammation and fibrosis that was highly dependent on genetic background, with vast differences in the degree of fibrosis. Histological analysis showed features characteristic of human NASH, including macrovesicular steatosis, hepatocellular ballooning, inflammatory foci, and pericellular collagen deposition. Time course experiments indicated that while hepatic triglyceride levels increased steadily on the diet, hepatic fibrosis occurred at about 12 weeks. We found that the genetic variation predisposing to NASH and fibrosis differs markedly from that predisposing to simple steatosis, consistent with a multistep model in which distinct genetic factors are involved. Moreover, genome-wide association identified distinct genetic loci contributing to steatosis and NASH. Finally, we used hepatic expression data from the mouse panel and from 68 bariatric surgery patients with normal liver, steatosis, or NASH to identify enriched biological pathways. Conclusion: The pathways showed substantial overlap between our mouse model and the human disease.

Inter-tissue communication via secreted proteins has been established as a vital mechanism for proper physiologic homeostasis. Here, we report a bioinformatics framework using a mouse reference population, the Hybrid Mouse Diversity Panel (HMDP), which integrates global multi-tissue expression data and publicly available resources to identify and functionally annotate novel circuits of tissue-tissue communication. We validate this method by showing that we can identify known as well as novel endocrine factors responsible for communication between tissues. We further show the utility of this approach by identification and mechanistic characterization of two new endocrine factors. Adipose-derived Lipocalin-5 is shown to enhance skeletal muscle mitochondrial function, and liver-secreted Notum promotes browning of white adipose tissue, also known as "beiging." We demonstrate the general applicability of the method by providing in vivo evidence for three additional novel molecules mediating tissue-tissue interactions.

We describe a simple bioinformatics method for biomarker discovery that is based on the analysis of global transcript levels in a population of inbred mouse strains showing variation for disease-related traits. This method has advantages such as controlled environment and accessibility to heart and plasma tissue in the preclinical selection stage. We illustrate the approach by identifying candidate heart failure (HF) biomarkers by overlaying mouse transcriptome and clinical traits from 91 Hybrid Mouse Diversity Panel (HMDP) inbred strains and human HF transcriptome from the Myocardial Applied Genomics Network (MAGNet) consortium. We found that some of the top differentially expressed genes correlated with known human HF biomarkers, such as galectin-3 and tissue inhibitor of metalloproteinase 1. Using ELISA assays, we investigated one novel candidate, Glycoprotein NMB, in a mouse model of chronic β-adrenergic stimulation by isoproterenol (ISO) induced HF. We observed significantly lower GPNMB plasma levels in the ISO model compared to the control group (p-value = 0.007). In addition, we assessed GPNMB plasma levels among 389 HF cases and controls from the METabolic Syndrome In Men (METSIM) study. Lower levels of GPNMB were also observed in patients with HF from the METSIM study compared to non-HF controls (p-value < 0.0001). In summary, we have identified several candidate biomarkers for HF using the cardiac transcriptome data in a population of mice that may be directly relevant and applicable to human populations.

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.

The simplest definition of cis-eQTLs versus trans, refers to genetic variants that affect expression in an allele specific manner, with implications on underlying mechanism. Yet, due to technical limitations of expression microarrays, the vast majority of eQTL studies performed in the last decade used a genomic distance based definition as a surrogate for cis, therefore exploring local rather than cis-eQTLs.

We investigated the association of the levels of ketone bodies (KBs) with hyperglycemia and with 62 genetic risk variants regulating glucose levels or type 2 diabetes in the population-based Metabolic Syndrome in Men (METSIM) study, including 9,398 Finnish men without diabetes or newly diagnosed type 2 diabetes. Increasing fasting and 2-h plasma glucose levels were associated with elevated levels of acetoacetate (AcAc) and β-hydroxybutyrate (BHB). AcAc and BHB predicted an increase in the glucose area under the curve in an oral glucose tolerance test, and AcAc predicted the conversion to type 2 diabetes in a 5-year follow-up of the METSIM cohort. Impaired insulin secretion, but not insulin resistance, explained these findings. Of the 62 single nucleotide polymorphisms associated with the risk of type 2 diabetes or hyperglycemia, the glucose-increasing C allele of GCKR significantly associated with elevated levels of fasting BHB levels. Adipose tissue mRNA expression levels of genes involved in ketolysis were significantly associated with insulin sensitivity (Matsuda index). In conclusion, high levels of KBs predicted subsequent worsening of hyperglycemia, and a common variant of GCKR was significantly associated with BHB levels.

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.

We report a systems genetic analysis of high density lipoprotein (HDL) levels in an F2 intercross between inbred strains CAST/EiJ and C57BL/6J. We previously showed that there are dramatic differences in HDL metabolism in a cross between these strains, and we now report co-expression network analysis of HDL that integrates global expression data from liver and adipose with relevant metabolic traits. Using data from a total of 293 F2 intercross mice, we constructed weighted gene co-expression networks and identified modules (subnetworks) associated with HDL and clinical traits. These were examined for genes implicated in HDL levels based on large human genome-wide associations studies (GWAS) and examined with respect to conservation between tissue and sexes in a total of 9 data sets. We identify genes that are consistently ranked high by association with HDL across the 9 data sets. We focus in particular on two genes, Wfdc2 and Hdac3, that are located in close proximity to HDL QTL peaks where causal testing indicates that they may affect HDL. Our results provide a rich resource for studies of complex metabolic interactions involving HDL. This article is part of a Special Issue entitled Advances in High Density Lipoprotein Formation and Metabolism: A Tribute to John F. Oram (1945-2010).

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.

We have examined maternal mechanisms for adult-onset glucose intolerance, increased adiposity, and atherosclerosis using two mouse models for intrauterine growth restriction (IUGR): maternal protein restriction and hypercholesterolemia.

Ambient air pollution is associated with increased cardiovascular morbidity and mortality. We have found that exposure to ambient ultrafine particulate matter, highly enriched in redox cycling organic chemicals, promotes atherosclerosis in mice. We hypothesize that these pro-oxidative chemicals could synergize with oxidized lipid components generated in low-density lipoprotein particles to enhance vascular inflammation and atherosclerosis.

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.