Genetic Regulation of Adipose Gene Expression and Cardio-Metabolic Traits.

Subcutaneous adipose tissue stores excess lipids and maintains energy balance. We performed expression quantitative trait locus (eQTL) analyses by using abdominal subcutaneous adipose tissue of 770 extensively phenotyped participants of the METSIM study. We identified cis-eQTLs for 12,400 genes at a 1% false-discovery rate. Among an approximately 680 known genome-wide association study (GWAS) loci for cardio-metabolic traits, we identified 140 coincident cis-eQTLs at 109 GWAS loci, including 93 eQTLs not previously described. At 49 of these 140 eQTLs, gene expression was nominally associated (p < 0.05) with levels of the GWAS trait. The size of our dataset enabled identification of five loci associated (p < 5 × 10-8) with at least five genes located >5 Mb away. These trans-eQTL signals confirmed and extended the previously reported KLF14-mediated network to 55 target genes, validated the CIITA regulation of class II MHC genes, and identified ZNF800 as a candidate master regulator. Finally, we observed similar expression-clinical trait correlations of genes associated with GWAS loci in both humans and a panel of genetically diverse mice. These results provide candidate genes for further investigation of their potential roles in adipose biology and in regulating cardio-metabolic traits.

A novel approach to measuring macrophage-specific reverse cholesterol transport in vivo in humans.

Reverse cholesterol transport (RCT) is thought to be an atheroprotective function of HDL, and macrophage-specific RCT in mice is inversely associated with atherosclerosis. We developed a novel method using 3H-cholesterol nanoparticles to selectively trace macrophage-specific RCT in vivo in humans. Use of 3H-cholesterol nanoparticles was initially tested in mice to assess the distribution of tracer and response to interventions known to increase RCT. Thirty healthy subjects received 3H-cholesterol nanoparticles intravenously, followed by blood and stool sample collection. Tracer counts were assessed in plasma, nonHDL, HDL, and fecal fractions. Data were analyzed by using multicompartmental modeling. Administration of 3H-cholesterol nanoparticles preferentially labeled macrophages of the reticuloendothelial system in mice, and counts were increased in mice treated with a liver X receptor agonist or reconstituted HDL, as compared with controls. In humans, tracer disappeared from plasma rapidly after injection of nanoparticles, followed by reappearance in HDL and nonHDL fractions. Counts present as free cholesterol increased rapidly and linearly in the first 240 min after nadir; counts in cholesteryl ester increased steadily over time. Estimates of fractional transfer rates of key RCT steps were obtained. These results support the use of 3H-cholesterol nanoparticles as a feasible approach for the measurement of macrophage RCT in vivo in humans.

Discovery of Highly Potent Liver X Receptor ß Agonists.

Introducing a uniquely substituted phenyl sulfone into a series of biphenyl imidazole liver X receptor (LXR) agonists afforded a dramatic potency improvement for induction of ATP binding cassette transporters, ABCA1 and ABCG1, in human whole blood. The agonist series demonstrated robust LXRß activity (>70%) with low partial LXRα agonist activity (<25%) in cell assays, providing a window between desired blood cell ABCG1 gene induction in cynomolgus monkeys and modest elevation of plasma triglycerides for agonist 15. The addition of polarity to the phenyl sulfone also reduced binding to the plasma protein, human α-1-acid glycoprotein. Agonist 15 was selected for clinical development based on the favorable combination of in vitro properties, excellent pharmacokinetic parameters, and a favorable lipid profile.

Beneficial and Adverse Effects of an LXR Agonist on Human Lipid and Lipoprotein Metabolism and Circulating Neutrophils.

The development of LXR agonists for the treatment of coronary artery disease has been challenged by undesirable properties in animal models. Here we show the effects of an LXR agonist on lipid and lipoprotein metabolism and neutrophils in human subjects. BMS-852927, a novel LXRß-selective compound, had favorable profiles in animal models with a wide therapeutic index in cynomolgus monkeys and mice. In healthy subjects and hypercholesterolemic patients, reverse cholesterol transport pathways were induced similarly to that in animal models. However, increased plasma and hepatic TG, plasma LDL-C, apoB, apoE, and CETP and decreased circulating neutrophils were also evident. Furthermore, similar increases in LDL-C were observed in normocholesterolemic subjects and statin-treated patients. The primate model markedly underestimated human lipogenic responses and did not predict human neutrophil effects. These studies demonstrate both beneficial and adverse LXR agonist clinical responses and emphasize the importance of further translational research in this area.

Loss of Function of GALNT2 Lowers High-Density Lipoproteins in Humans, Nonhuman Primates, and Rodents.

Human genetics studies have implicated GALNT2, encoding GalNAc-T2, as a regulator of high-density lipoprotein cholesterol (HDL-C) metabolism, but the mechanisms relating GALNT2 to HDL-C remain unclear. We investigated the impact of homozygous GALNT2 deficiency on HDL-C in humans and mammalian models. We identified two humans homozygous for loss-of-function mutations in GALNT2 who demonstrated low HDL-C. We also found that GALNT2 loss of function in mice, rats, and nonhuman primates decreased HDL-C. O-glycoproteomics studies of a human GALNT2-deficient subject validated ANGPTL3 and ApoC-III as GalNAc-T2 targets. Additional glycoproteomics in rodents identified targets influencing HDL-C, including phospholipid transfer protein (PLTP). GALNT2 deficiency reduced plasma PLTP activity in humans and rodents, and in mice this was rescued by reconstitution of hepatic Galnt2. We also found that GALNT2 GWAS SNPs associated with reduced HDL-C also correlate with lower hepatic GALNT2 expression. These results posit GALNT2 as a direct modulator of HDL metabolism across mammals.

Genetic Architecture of Atherosclerosis in Mice: A Systems Genetics Analysis of Common Inbred Strains.

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.

The genetic architecture of NAFLD among inbred strains of mice.

To identify genetic and environmental factors contributing to the pathogenesis of non-alcoholic fatty liver disease, we examined liver steatosis and related clinical and molecular traits in more than 100 unique inbred mouse strains, which were fed a diet rich in fat and carbohydrates. A >30-fold variation in hepatic TG accumulation was observed among the strains. Genome-wide association studies revealed three loci associated with hepatic TG accumulation. Utilizing transcriptomic data from the liver and adipose tissue, we identified several high-confidence candidate genes for hepatic steatosis, including Gde1, a glycerophosphodiester phosphodiesterase not previously implicated in triglyceride metabolism. We confirmed the role of Gde1 by in vivo hepatic over-expression and shRNA knockdown studies. We hypothesize that Gde1 expression increases TG production by contributing to the production of glycerol-3-phosphate. Our multi-level data, including transcript levels, metabolite levels, and gut microbiota composition, provide a framework for understanding genetic and environmental interactions underlying hepatic steatosis.

Genetic architecture of insulin resistance in the mouse.

Insulin resistance (IR) is a complex trait with multiple genetic and environmental components. Confounded by large differences between the sexes, environment, and disease pathology, the genetic basis of IR has been difficult to dissect. Here we examine IR and related traits in a diverse population of more than 100 unique male and female inbred mouse strains after feeding a diet rich in fat and refined carbohydrates. Our results show dramatic variation in IR among strains of mice and widespread differences between sexes that are dependent on genotype. We uncover more than 15 genome-wide significant loci and validate a gene, Agpat5, associated with IR. We also integrate plasma metabolite levels and global gene expression from liver and adipose tissue to identify metabolite quantitative trait loci (mQTL) and expression QTL (eQTL), respectively. Our results provide a resource for analysis of interactions between diet, sex, and genetic background in IR.

Pharmacological characterization of a novel liver X receptor agonist with partial LXRα activity and a favorable window in nonhuman primates.

Liver X Receptors (LXRs) α and ß are nuclear hormone receptors that regulate multiple genes involved in reverse cholesterol transport (RCT) and are potential drug targets for atherosclerosis. However, full pan agonists also activate lipogenic genes, resulting in elevated plasma and hepatic lipids. We report the pharmacology of BMS-779788 [2-(2-(1-(2-chlorophenyl)-1-methylethyl)-1-(3'-(methylsulfonyl)-4-biphenylyl)-1H-imidazol-4-yl)-2-propanol], a potent partial LXR agonist with LXRß selectivity, which has an improved therapeutic window in the cynomolgus monkey compared with a full pan agonist. BMS-779788 induced LXR target genes in blood in vivo with an EC50 = 610 nM, a value similar to its in vitro blood gene induction potency. BMS-779788 was 29- and 12-fold less potent than the full agonist T0901317 in elevating plasma triglyceride and LDL cholesterol, respectively, with similar results for plasma cholesteryl ester transfer protein and apolipoprotein B. However, ABCA1 and ABCG1 mRNA inductions in blood, which are critical for RCT, were comparable. Increased liver triglyceride was observed after 7-day treatment with BMS-779788 at the highest dose tested and was nearly identical to the dose response for plasma triglyceride, consistent with the central role of liver LXR in these lipogenic effects. Dose-dependent increases in biliary cholesterol and decreases in phospholipid and bile acid occurred in BMS-779788-treated animals, similar to LXR agonist effects reported in mouse. In summary, BMS-779788, a partial LXRß selective agonist, has decreased lipogenic potential compared with a full pan agonist in cynomolgus monkeys, with similar potency in the induction of genes known to stimulate RCT. This provides support in nonhuman primates for improving LXR agonist therapeutic windows by limiting LXRα activity.

Liver X receptor (LXR) partial agonists: biaryl pyrazoles and imidazoles displaying a preference for LXRß.

A series of biaryl pyrazole and imidazole Liver X Receptor (LXR) partial agonists has been synthesized displaying LXRß selectivity. The LXRß selective partial agonist 18 was identified with potent induction of ATP binding transporters ABCA1 and ABCG1 in human whole blood (EC50=1.2µM, 55% efficacy). In mice 18 displayed peripheral induction of ABCA1 at 3 and 10mpk doses with no significant elevation of plasma or hepatic triglycerides at these doses, showing an improved profile compared to a full pan-agonist.

Genetic regulation of mouse liver metabolite levels.

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.

Genetic regulation of human adipose microRNA expression and its consequences for metabolic traits.

The genetics of messenger RNA (mRNA) expression has been extensively studied in humans and other organisms, but little is known about genetic factors contributing to microRNA (miRNA) expression. We examined natural variation of miRNA expression in adipose tissue in a population of 200 men who have been carefully characterized for metabolic syndrome (MetSyn) phenotypes as part of the Metabolic Syndrome in Men (METSIM) study. We genotyped the subjects using high-density single-nucleotide polymorphism microarrays and quantified the mRNA abundance using genome-wide expression arrays and miRNA abundance using next-generation sequencing. We reliably quantified 356 miRNA species that were expressed in human adipose tissue, a limited number of which made up most of the expressed miRNAs. We mapped the miRNA abundance as an expression quantitative trait and determined cis regulation of expression for nine of the miRNAs and of the processing of one miRNA (miR-28). The degree of genetic variation of miRNA expression was substantially less than that of mRNAs. For the majority of the miRNAs, genetic regulation of expression was independent of the expression of mRNA from which the miRNA is transcribed. We also showed that for 108 miRNAs, mapped reads displayed widespread variation from the canonical sequence. We found a total of 24 miRNAs to be significantly associated with MetSyn traits. We suggest a regulatory role for miR-204-5p which was predicted to inhibit acetyl coenzyme A carboxylase ß, a key fatty acid oxidation enzyme that has been shown to play a role in regulating body fat and insulin resistance in adipose tissue.

Genetic control of obesity and gut microbiota composition in response to high-fat, high-sucrose diet in mice.

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.

Unraveling inflammatory responses using systems genetics and gene-environment interactions in macrophages.

Many common diseases have an important inflammatory component mediated in part by macrophages. Here we used a systems genetics strategy to examine the role of common genetic variation in macrophage responses to inflammatory stimuli. We examined genome-wide transcript levels in macrophages from 92 strains of the Hybrid Mouse Diversity Panel. We exposed macrophages to control media, bacterial lipopolysaccharide (LPS), or oxidized phospholipids. We performed association mapping under each condition and identified several thousand expression quantitative trait loci (eQTL), gene-by-environment interactions, and eQTL "hot spots" that specifically control LPS responses. We used siRNA knockdown of candidate genes to validate an eQTL hot spot in chromosome 8 and identified the gene 2310061C15Rik as a regulator of inflammatory responses in macrophages. We have created a public database where the data presented here can be used as a resource for understanding many common inflammatory traits that are modeled in the mouse and for the dissection of regulatory relationships between genes.

Hybrid mouse diversity panel: a panel of inbred mouse strains suitable for analysis of complex genetic traits.

We have developed an association-based approach using classical inbred strains of mice in which we correct for population structure, which is very extensive in mice, using an efficient mixed-model algorithm. Our approach includes inbred parental strains as well as recombinant inbred strains in order to capture loci with effect sizes typical of complex traits in mice (in the range of 5% of total trait variance). Over the last few years, we have typed the hybrid mouse diversity panel (HMDP) strains for a variety of clinical traits as well as intermediate phenotypes and have shown that the HMDP has sufficient power to map genes for highly complex traits with resolution that is in most cases less than a megabase. In this essay, we review our experience with the HMDP, describe various ongoing projects, and discuss how the HMDP may fit into the larger picture of common diseases and different approaches.

Genome-wide study of gene variants associated with differential cardiovascular event reduction by pravastatin therapy.

Statin therapy reduces the risk of coronary heart disease (CHD), however, the person-to-person variability in response to statin therapy is not well understood. We have investigated the effect of genetic variation on the reduction of CHD events by pravastatin. First, we conducted a genome-wide association study of 682 CHD cases from the Cholesterol and Recurrent Events (CARE) trial and 383 CHD cases from the West of Scotland Coronary Prevention Study (WOSCOPS), two randomized, placebo-controlled studies of pravastatin. In a combined case-only analysis, 79 single nucleotide polymorphisms (SNPs) were associated with differential CHD event reduction by pravastatin according to genotype (P<0.0001), and these SNPs were analyzed in a second stage that included cases as well as non-cases from CARE and WOSCOPS and patients from the PROspective Study of Pravastatin in the Elderly at Risk/PHArmacogenomic study of Statins in the Elderly at risk for cardiovascular disease (PROSPER/PHASE), a randomized placebo controlled study of pravastatin in the elderly. We found that one of these SNPs (rs13279522) was associated with differential CHD event reduction by pravastatin therapy in all 3 studies: P = 0.002 in CARE, P = 0.01 in WOSCOPS, P = 0.002 in PROSPER/PHASE. In a combined analysis of CARE, WOSCOPS, and PROSPER/PHASE, the hazard ratio for CHD when comparing pravastatin with placebo decreased by a factor of 0.63 (95% CI: 0.52 to 0.75) for each extra copy of the minor allele (P = 4.8 × 10(-7)). This SNP is located in DnaJ homolog subfamily C member 5B (DNAJC5B) and merits investigation in additional randomized studies of pravastatin and other statins.

Comparative analysis of proteome and transcriptome variation in mouse.

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.

Genetic variants in the KIF6 region and coronary event reduction from statin therapy.

A single nucleotide polymorphism (SNP) in KIF6, a member of the KIF9 family of kinesins, is associated with differential coronary event reduction from statin therapy in four randomized controlled trials; this SNP (rs20455) is also associated with the risk for coronary heart disease (CHD) in multiple prospective studies. We investigated whether other common SNPs in the KIF6 region were associated with event reduction from statin therapy. Of the 170 SNPs in the KIF6 region investigated in the Cholesterol and Recurrent Events trial (CARE), 28 were associated with differential event reduction from statin therapy (P (interaction) < 01 in Caucasians, adjusted for age and sex) and were further investigated in the Pravastatin or Atorvastatin Evaluation and Infection Therapy-Thrombolysis In Myocardial Infarction 22 (PROVE IT-TIMI22) and West of Scotland Coronary Prevention Study (WOSCOPS). These analyses revealed that two SNPs (rs9462535 and rs9471077), in addition to rs20455, were associated with event reduction from statin therapy (P (interaction) < 0.1 in each of the three studies). The relative risk reduction ranged from 37 to 50% (P < 0.01) in carriers of the minor alleles of these SNPs and from -4 to 13% (P > 0.4) in non-carriers. These three SNPs are in high linkage disequilibrium with one another (r (2) > 0.84). Functional studies of these variants may help to understand the role of KIF6 in the pathogenesis of CHD and differential response to statin therapy.

Effect of pravastatin therapy on coronary events in carriers of the KIF6 719Arg allele from the cholesterol and recurrent events trial.

A previous genetic analysis of the Cholesterol and Recurrent Events (CARE) trial found that carriers of the 719Arg allele of the kinesin family member 6 gene (KIF6) (rs20455), but not noncarriers, received significant event reduction from pravastatin therapy. However, that previous analysis of CARE included only Caucasian patients and was limited to the myocardial infarction components of the primary end point. Therefore, the aim of this study was to investigate whether pravastatin therapy reduced primary end point events in KIF6 719Arg carriers and noncarriers, separately, in the combined ethnic groups of CARE. The effect of pravastatin therapy on primary end point events (fatal coronary event or nonfatal myocardial infarction) was investigated in Cox regression models that adjusted for population structure using either self-reported ethnicity or the principal components of genetic heterogeneity. After adjustment for age, gender, and self-reported ethnicity, pravastatin therapy reduced events in carriers of KIF6 719Arg (hazard ratio [HR] 0.63, 95% confidence interval [CI] 0.49 to 0.83) but not in noncarriers (HR 1.01, 95% CI 0.69 to 1.45) (p for interaction = 0.049). After adjustment for age, gender, traditional risk factors, and principal components, pravastatin therapy reduced events in carriers of 719Arg (HR 0.64, 95% CI 0.49 to 0.85) but not in noncarriers (HR 0.90, 95% CI 0.62 to 1.32) (p for interaction = 0.14). In conclusion, in an analysis that included CARE patients of all ethnic groups, pravastatin therapy significantly and substantially reduced primary end point events in carriers of the KIF6 719Arg allele but not in noncarriers.