Genetic Regulation of Adipose Gene Expression and Cardio-Metabolic Traits

Identification and functional analysis of glycemic trait loci in the China Health and Nutrition Survey

To identify genetic contributions to type 2 diabetes (T2D) and related glycemic traits (fasting glucose, fasting insulin, and HbA1c), we conducted genome-wide association analyses (GWAS) in up to 7,178 Chinese subjects from nine provinces in the China Health and Nutrition Survey (CHNS). We examined patterns of population structure within CHNS and found that allele frequencies differed across provinces, consistent with genetic drift and population substructure. We further validated 32 previously described T2D- and glycemic trait-loci, including G6PC2 and SIX3-SIX2 associated with fasting glucose. At G6PC2, we replicated a known fasting glucose-associated variant (rs34177044) and identified a second signal (rs2232326), a low-frequency (4%), probably damaging missense variant (S324P). A variant within the lead fasting glucose-associated signal at SIX3-SIX2 co-localized with pancreatic islet expression quantitative trait loci (eQTL) for SIX3, SIX2, and three noncoding transcripts. To identify variants functionally responsible for the fasting glucose association at SIX3-SIX2, we tested five candidate variants for allelic differences in regulatory function. The rs12712928-C allele, associated with higher fasting glucose and lower transcript expression level, showed lower transcriptional activity in reporter assays and increased binding to GABP compared to the rs12712928-G, suggesting that rs12712928-C contributes to elevated fasting glucose levels by disrupting an islet enhancer, resulting in reduced gene expression. Taken together, these analyses identified multiple loci associated with glycemic traits across China, and suggest a regulatory mechanism at the SIX3-SIX2 fasting glucose GWAS locus.

Regulatory variants at KLF14 influence type 2 diabetes risk via a female-specific effect on adipocyte size and body composition

Individual risk of type 2 diabetes (T2D) is modified by perturbations to the mass, distribution and function of adipose tissue. To investigate the mechanisms underlying these associations, we explored the molecular, cellular and whole-body effects of T2D-associated alleles near KLF14. We show that KLF14 diabetes-risk alleles act in adipose tissue to reduce KLF14 expression and modulate, in trans, the expression of 385 genes. We demonstrate, in human cellular studies, that reduced KLF14 expression increases pre-adipocyte proliferation but disrupts lipogenesis, and in mice, that adipose tissue-specific deletion of Klf14 partially recapitulates the human phenotype of insulin resistance, dyslipidemia and T2D. We show that carriers of the KLF14 T2D risk allele shift body fat from gynoid stores to abdominal stores and display a marked increase in adipocyte cell size, and that these effects on fat distribution, and the T2D association, are female specific. The metabolic risk associated with variation at this imprinted locus depends on the sex both of the subject and of the parent from whom the risk allele derives.

The TMAO-Producing Enzyme Flavin-Containing Monooxygenase 3 Regulates Obesity and the Beiging of White Adipose Tissue

Emerging evidence suggests that microbes resident in the human intestine represent a key environmental factor contributing to obesity-associated disorders. Here, we demonstrate that the gut microbiota-initiated trimethylamine N-oxide (TMAO)-generating pathway is linked to obesity and energy metabolism. In multiple clinical cohorts, systemic levels of TMAO were observed to strongly associate with type 2 diabetes. In addition, circulating TMAO levels were associated with obesity traits in the different inbred strains represented in the Hybrid Mouse Diversity Panel. Further, antisense oligonucleotide-mediated knockdown or genetic deletion of the TMAO-producing enzyme flavin-containing monooxygenase 3 (FMO3) conferred protection against obesity in mice. Complimentary mouse and human studies indicate a negative regulatory role for FMO3 in the beiging of white adipose tissue. Collectively, our studies reveal a link between the TMAO-producing enzyme FMO3 and obesity and the beiging of white adipose tissue.

Transcriptional Regulatory Mechanisms in Adipose and Muscle Tissue Associated with Composite Glucometabolic Phenotypes

OBJECTIVE

Tissue-specific gene expression is associated with individual metabolic measures. However, these measures may not reflect the true but latent underlying biological phenotype. This study reports gene expression associations with multidimensional glucometabolic characterizations of obesity, glucose homeostasis, and lipid traits.

METHODS

Factor analysis was computed by using orthogonal rotation to construct composite phenotypes (CPs) from 23 traits in 256 African Americans without diabetes. Genome-wide transcript expression data from adipose and muscle were tested for association with CPs, and expression quantitative trait loci (eQTLs) were identified by associations between cis-acting single-nucleotide polymorphisms (SNPs) and gene expression.

RESULTS

The factor analysis identified six CPs. CPs 1 through 6 individually explained 34%, 12%, 9%, 8%, 6%, and 5% of the variation in 23 glucometabolic traits studied. There were 3,994 and 929 CP-associated transcripts identified in adipose and muscle tissue, respectively; CP2 had the largest number of associated transcripts. Pathway analysis identified multiple canonical pathways from the CP-associated transcripts. In adipose and muscle, significant cis-eQTLs were identified for 558 and 164 CP-associated transcripts (q-value < 0.01), respectively.

CONCLUSIONS

Adipose and muscle transcripts comprehensively define pathways involved in regulating glucometabolic disorders. Cis-eQTLs for CP-associated genes may act as primary causal determinants of glucometabolic phenotypes by regulating transcription of key genes.

From Genotype to Phenotype: A Primer on the Functional Follow-up of Genome-Wide Association Studies in Cardiovascular Disease

Genome-wide association studies (GWASs) have implicated many human genomic loci in the development of complex traits. The loci identified by these studies are potentially involved in novel pathways that contribute to disease pathophysiology. However, eventual therapeutic targeting of these pathways relies on bridging the gap between genetic association and function, a task that first requires validation of causal genetic variants, casual genes, and directionality of effect. Executing this task requires basic knowledge of interpreting GWAS results and prioritizing candidates for further study, in addition to understanding the experimental methods available for evaluating candidate variants. Here we review the basic genetic principles of genome-wide association studies, the computational and experimental tools used for identifying causal variants and genes, and salient illustrative examples of how cardiovascular loci have undergone functional investigation.

Voluntary exposure to a toxin: the genetic influence on ethanol consumption

Ethyl alcohol is a toxin that, when consumed at high levels, produces organ damage and death. One way to prevent or ameliorate this damage in humans is to reduce the exposure of organs to alcohol by reducing alcohol ingestion. Both the propensity to consume large volumes of alcohol and the susceptibility of human organs to alcohol-induced damage exhibit a strong genetic influence. We have developed an integrative genetic/genomic approach to identify transcriptional networks that predispose complex traits, including propensity for alcohol consumption and propensity for alcohol-induced organ damage. In our approach, the phenotype is assessed in a panel of recombinant inbred (RI) rat strains, and quantitative trait locus (QTL) analysis is performed. Transcriptome data from tissues/organs of naïve RI rat strains are used to identify transcriptional networks using Weighted Gene Coexpression Network Analysis (WGCNA). Correlation of the first principal component of transcriptional coexpression modules with the phenotype across the rat strains, and overlap of QTLs for the phenotype and the QTLs for the coexpression modules (module eigengene QTL) provide the criteria for identification of the functionally related groups of genes that contribute to the phenotype (candidate modules). While we previously identified a brain transcriptional module whose QTL overlapped with a QTL for levels of alcohol consumption in HXB/BXH RI rat strains and 12 selected rat lines, this module did not account for all of the genetic variation in alcohol consumption. Our search for QTL overlap and correlation of coexpression modules with phenotype can, however, be applied to any organ in which the transcriptome has been measured, and this represents a holistic approach in the search for genetic contributors to complex traits. Previous work has implicated liver/brain interactions, particularly involving inflammatory/immune processes, as influencing alcohol consumption levels. We have now analyzed the liver transcriptome of the HXB/BXH RI rat panel in relation to the behavioral trait of alcohol consumption. We used RNA-Seq and microarray data to construct liver transcriptional networks, and identified a liver candidate transcriptional coexpression module that explained 24% of the genetic variance in voluntary alcohol consumption. The transcripts in this module focus attention on liver secretory products that influence inflammatory and immune signaling pathways. We propose that these liver secretory products can interact with brain mechanisms that affect alcohol consumption, and targeting these pathways provides a potential approach to reducing high levels of alcohol intake and also protecting the integrity of the liver and other organs.

HUGIn: Hi-C Unifying Genomic Interrogator

MOTIVATION

High throughput chromatin conformation capture (3C) technologies, such as Hi-C and ChIA-PET, have the potential to elucidate the functional roles of non-coding variants. However, most of published genome-wide unbiased chromatin organization studies have used cultured cell lines, limiting their generalizability.

RESULTS

We developed a web browser, HUGIn, to visualize Hi-C data generated from 21 human primary tissues and cell lines. HUGIn enables assessment of chromatin contacts both constitutive across and specific to tissue(s) and/or cell line(s) at any genomic loci, including GWAS SNPs, eQTLs and cis-regulatory elements, facilitating the understanding of both GWAS and eQTL results and functional genomics data.

AVAILABILITY AND IMPLEMENTATION

HUGIn is available at http://yunliweb.its.unc.edu/HUGIn.

CONTACT

yunli@med.unc.edu or hum@ccf.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

HDL Cholesterol Metabolism and the Risk of CHD: New Insights from Human Genetics

PURPOSE OF REVIEW

Elevated high-density lipoprotein cholesterol levels in the blood (HDL-C) represent one of the strongest epidemiological surrogates for protection against coronary heart disease (CHD), but recent human genetic and pharmacological intervention studies have raised controversy about the causality of this relationship. Here, we review recent discoveries from human genome studies using new analytic tools as well as relevant animal studies that have both addressed, and in some cases, fueled this controversy.

RECENT FINDINGS

Methodologic developments in genotyping and sequencing, such as genome-wide association studies (GWAS), exome sequencing, and exome array genotyping, have been applied to the study of HDL-C and risk of CHD in large, multi-ethnic populations. Some of these efforts focused on population-wide variation in common variants have uncovered new polymorphisms at novel loci associated with HDL-C and, in some cases, CHD risk. Other efforts have discovered loss-of-function variants for the first time in genes previously implicated in HDL metabolism through common variant studies or animal models. These studies have allowed the genetic relationship between these pathways, HDL-C and CHD to be explored in humans for the first time through analysis tools such as Mendelian randomization. We explore these discoveries for selected key HDL-C genes CETP, LCAT, LIPG, SCARB1, and novel loci implicated from GWAS including GALNT2, KLF14, and TTC39B. Recent human genetics findings have identified new nodes regulating HDL metabolism while reshaping our current understanding of known candidate genes to HDL and CHD risk through the study of critical variants across model systems. Despite their effect on HDL-C, variants in many of the reviewed genes were found to lack any association with CHD. These data collectively indicate that HDL-C concentration, which represents a static picture of a very dynamic and heterogeneous metabolic milieu, is unlikely to be itself causally protective against CHD. In this context, human genetics represent an extremely valuable tool to further explore the biological mechanisms regulating HDL metabolism and investigate what role, if any, HDL plays in the pathogenesis of CHD.

Genetic regulation of adipose tissue transcript expression is involved in modulating serum triglyceride and HDL-cholesterol

Dyslipidemia is a major contributor to the increased cardiovascular disease and mortality associated with obesity and type 2 diabetes. We hypothesized that variation in expression of adipose tissue transcripts is associated with serum lipid concentrations in African Americans (AAs), and common genetic variants regulate expression levels of these transcripts. Fasting serum lipid levels, genome-wide transcript expression profiles of subcutaneous adipose tissue, and genome-wide SNP genotypes were analyzed in a cohort of non-diabetic AAs (N=250). Serum triglyceride (TRIG) and high density lipoprotein-cholesterol (HDL-C) levels were associated (FDR<0.01) with expression level of 1021 and 1875 adipose tissue transcripts, respectively, but none associated with total cholesterol or LDL-C levels. Serum HDL-C-associated transcripts were enriched for salient biological pathways, including branched-chain amino acid degradation, and oxidative phosphorylation. Genes in immuno-inflammatory pathways were activated among individuals with higher serum TRIG levels. We identified significant cis-regulatory SNPs (cis-eSNPs) for 449 serum lipid-associated transcripts in adipose tissue. The cis-eSNPs of 12 genes were nominally associated (p<0.001) with serum lipid level in genome wide association studies in Global Lipids Genetics Consortium (GLGC) cohorts. Allelic effect direction of cis-eSNPs on expression of MARCH2, BEST1 and TMEM258 matched with effect direction of these SNP alleles on serum TRIG or HDL-C levels in GLGC cohorts. These data suggest that expressions of serum lipid-associated transcripts in adipose tissue are dependent on common cis-eSNPs in African Americans. Thus, genetically-mediated transcriptional regulation in adipose tissue may play a role in reducing HDL-C and increasing TRIG in serum.

Association analyses of East Asian individuals and trans-ancestry analyses with European individuals reveal new loci associated with cholesterol and triglyceride levels

Large-scale meta-analyses of genome-wide association studies (GWAS) have identified >175 loci associated with fasting cholesterol levels, including total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG). With differences in linkage disequilibrium (LD) structure and allele frequencies between ancestry groups, studies in additional large samples may detect new associations. We conducted staged GWAS meta-analyses in up to 69,414 East Asian individuals from 24 studies with participants from Japan, the Philippines, Korea, China, Singapore, and Taiwan. These meta-analyses identified (P < 5 × 10-8) three novel loci associated with HDL-C near CD163-APOBEC1 (P = 7.4 × 10-9), NCOA2 (P = 1.6 × 10-8), and NID2-PTGDR (P = 4.2 × 10-8), and one novel locus associated with TG near WDR11-FGFR2 (P = 2.7 × 10-10). Conditional analyses identified a second signal near CD163-APOBEC1. We then combined results from the East Asian meta-analysis with association results from up to 187,365 European individuals from the Global Lipids Genetics Consortium in a trans-ancestry meta-analysis. This analysis identified (log10Bayes Factor ≥6.1) eight additional novel lipid loci. Among the twelve total loci identified, the index variants at eight loci have demonstrated at least nominal significance with other metabolic traits in prior studies, and two loci exhibited coincident eQTLs (P < 1 × 10-5) in subcutaneous adipose tissue for BPTF and PDGFC. Taken together, these analyses identified multiple novel lipid loci, providing new potential therapeutic targets.

Trans-ancestry Fine Mapping and Molecular Assays Identify Regulatory Variants at the ANGPTL8 HDL-C GWAS Locus

Recent genome-wide association studies (GWAS) have identified variants associated with high-density lipoprotein cholesterol (HDL-C) located in or near the ANGPTL8 gene. Given the extensive sharing of GWAS loci across populations, we hypothesized that at least one shared variant at this locus affects HDL-C. The HDL-C–associated variants are coincident with expression quantitative trait loci for ANGPTL8 and DOCK6 in subcutaneous adipose tissue; however, only ANGPTL8 expression levels are associated with HDL-C levels. We identified a 400-bp promoter region of ANGPTL8 and enhancer regions within 5 kb that contribute to regulating expression in liver and adipose. To identify variants functionally responsible for the HDL-C association, we performed fine-mapping analyses and selected 13 candidate variants that overlap putative regulatory regions to test for allelic differences in regulatory function. Of these variants, rs12463177-G increased transcriptional activity (1.5-fold, P = 0.004) and showed differential protein binding. Six additional variants (rs17699089, rs200788077, rs56322906, rs3760782, rs737337, and rs3745683) showed evidence of allelic differences in transcriptional activity and/or protein binding. Taken together, these data suggest a regulatory mechanism at the ANGPTL8 HDL-C GWAS locus involving tissue-selective expression and at least one functional variant.

A genomic exploration identifies mechanisms that may explain adverse cardiovascular effects of COX-2 inhibitors

Cyclooxygenase-2 inhibitors (coxibs) are characterized by multiple molecular off-target effects and increased coronary artery disease (CAD) risk. Here, we systematically explored common variants of genes representing molecular targets of coxibs for association with CAD. Given a broad spectrum of pleiotropic effects of coxibs, our intention was to narrow potential mechanisms affecting CAD risk as we hypothesized that the affected genes may also display genomic signals of coronary disease risk. A Drug Gene Interaction Database search identified 47 gene products to be affected by coxibs. We traced association signals in 200-kb regions surrounding these genes in 84,813 CAD cases and 202,543 controls. Based on a threshold of 1 × 10⁻⁵ (Bonferroni correction for 3131 haplotype blocks), four gene loci yielded significant associations. The lead SNPs were rs7270354 (MMP9), rs4888383 (BCAR1), rs6905288 (VEGFA1), and rs556321 (CACNA1E). By additional genotyping, rs7270354 at MMP9 and rs4888383 at BCAR1 also reached the established GWAS threshold for genome-wide significance. The findings demonstrate overlap of genes affected by coxibs and those mediating CAD risk and points to further mechanisms, which are potentially responsible for coxib-associated CAD risk. The novel approach furthermore suggests that genetic studies may be useful to explore the clinical relevance of off-target drug effects.

Epistatic Networks Jointly Influence Phenotypes Related to Metabolic Disease and Gene Expression in Diversity Outbred Mice

In this study, Tyler et al. analyzed the complex genetic architecture of metabolic disease-related traits using the Diversity Outbred mouse population

Genetic studies of multidimensional phenotypes can potentially link genetic variation, gene expression, and physiological data to create multi-scale models of complex traits. The challenge of reducing these data to specific hypotheses has become increasingly acute with the advent of genome-scale data resources. Multi-parent populations derived from model organisms provide a resource for developing methods to understand this complexity. In this study, we simultaneously modeled body composition, serum biomarkers, and liver transcript abundances from 474 Diversity Outbred mice. This population contained both sexes and two dietary cohorts. Transcript data were reduced to functional gene modules with weighted gene coexpression network analysis (WGCNA), which were used as summary phenotypes representing enriched biological processes. These module phenotypes were jointly analyzed with body composition and serum biomarkers in a combined analysis of pleiotropy and epistasis (CAPE), which inferred networks of epistatic interactions between quantitative trait loci that affect one or more traits. This network frequently mapped interactions between alleles of different ancestries, providing evidence of both genetic synergy and redundancy between haplotypes. Furthermore, a number of loci interacted with sex and diet to yield sex-specific genetic effects and alleles that potentially protect individuals from the effects of a high-fat diet. Although the epistatic interactions explained small amounts of trait variance, the combination of directional interactions, allelic specificity, and high genomic resolution provided context to generate hypotheses for the roles of specific genes in complex traits. Our approach moves beyond the cataloging of single loci to infer genetic networks that map genetic etiology by simultaneously modeling all phenotypes.

A Low-Frequency Inactivating AKT2 Variant Enriched in the Finnish Population Is Associated With Fasting Insulin Levels and Type 2 Diabetes Risk

To identify novel coding association signals and facilitate characterization of mechanisms influencing glycemic traits and type 2 diabetes risk, we analyzed 109,215 variants derived from exome array genotyping together with an additional 390,225 variants from exome sequence in up to 39,339 normoglycemic individuals from five ancestry groups. We identified a novel association between the coding variant (p.Pro50Thr) in AKT2 and fasting plasma insulin (FI), a gene in which rare fully penetrant mutations are causal for monogenic glycemic disorders. The low-frequency allele is associated with a 12% increase in FI levels. This variant is present at 1.1% frequency in Finns but virtually absent in individuals from other ancestries. Carriers of the FI-increasing allele had increased 2-h insulin values, decreased insulin sensitivity, and increased risk of type 2 diabetes (odds ratio 1.05). In cellular studies, the AKT2-Thr50 protein exhibited a partial loss of function. We extend the allelic spectrum for coding variants in AKT2 associated with disorders of glucose homeostasis and demonstrate bidirectional effects of variants within the pleckstrin homology domain of AKT2.

Multi-omics approaches to disease

High-throughput technologies have revolutionized medical research. The advent of genotyping arrays enabled large-scale genome-wide association studies and methods for examining global transcript levels, which gave rise to the field of “integrative genetics”. Other omics technologies, such as proteomics and metabolomics, are now often incorporated into the everyday methodology of biological researchers. In this review, we provide an overview of such omics technologies and focus on methods for their integration across multiple omics layers. As compared to studies of a single omics type, multi-omics offers the opportunity to understand the flow of information that underlies disease.

The Metabolic Syndrome in Men study: a resource for studies of metabolic and cardiovascular diseases

The Metabolic Syndrome in Men (METSIM) study is a population-based study including 10,197 Finnish men examined in 2005-2010. The aim of the study is to investigate nongenetic and genetic factors associated with the risk of T2D and CVD, and with cardiovascular risk factors. The protocol includes a detailed phenotyping of the participants, an oral glucose tolerance test, fasting laboratory measurements including proton NMR measurements, mass spectometry metabolomics, adipose tissue biopsies from 1,400 participants, and a stool sample. In our ongoing follow-up study, we have, to date, reexamined 6,496 participants. Extensive genotyping and exome sequencing have been performed for essentially all METSIM participants, and >2,000 METSIM participants have been whole-genome sequenced. We have identified several nongenetic markers associated with the development of diabetes and cardiovascular events, and participated in several genetic association studies to identify gene variants associated with diabetes, hyperglycemia, and cardiovascular risk factors. The generation of a phenotype and genotype resource in the METSIM study allows us to proceed toward a "systems genetics" approach, which includes statistical methods to quantitate and integrate intermediate phenotypes, such as transcript, protein, or metabolite levels, to provide a global view of the molecular architecture of complex traits.