Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explain one-fifth of heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ~2,000, ~3,700 and ~9,500 SNPs explained ~21%, ~24% and ~29% of phenotypic variance. Furthermore, all common variants together captured the majority (60%) of heritability. The 697 variants clustered in 423 loci enriched for genes, pathways, and tissue-types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/beta-catenin, and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, we conducted genome-wide association meta-analyses of waist and hip circumference-related traits in up to 224,459 individuals. We identified 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (WHRadjBMI) and an additional 19 loci newly associated with related waist and hip circumference measures (P<5×10−8). Twenty of the 49 WHRadjBMI loci showed significant sexual dimorphism, 19 of which displayed a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation, and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.

Sequence-based variation in gene expression is a key driver of disease risk. Common variants regulating expression in cis have been mapped in many expression quantitative trait locus (eQTL) studies, typically in single tissues from unrelated individuals. Here, we present a comprehensive analysis of gene expression across multiple tissues conducted in a large set of mono- and dizygotic twins that allows systematic dissection of genetic (cis and trans) and non-genetic effects on gene expression. Using identity-by-descent estimates, we show that at least 40% of the total heritable cis effect on expression cannot be accounted for by common cis variants, a finding that reveals the contribution of low-frequency and rare regulatory variants with respect to both transcriptional regulation and complex trait susceptibility. We show that a substantial proportion of gene expression heritability is trans to the structural gene, and we identify several replicating trans variants that act predominantly in a tissue-restricted manner and may regulate the transcription of many genes

While there have been studies exploring regulatory variation in one or more tissues, the complexity of tissue-specificity in multiple primary tissues is not yet well understood. We explore in depth the role of cis-regulatory variation in three human tissues: lymphoblastoid cell lines (LCL), skin, and fat. The samples (156 LCL, 160 skin, 166 fat) were derived simultaneously from a subset of well-phenotyped healthy female twins of the MuTHER resource. We discover an abundance of cis-eQTLs in each tissue similar to previous estimates (858 or 4.7% of genes). In addition, we apply factor analysis (FA) to remove effects of latent variables, thus more than doubling the number of our discoveries (1,822 eQTL genes). The unique study design (Matched Co-Twin Analysis—MCTA) permits immediate replication of eQTLs using co-twins (93%–98%) and validation of the considerable gain in eQTL discovery after FA correction. We highlight the challenges of comparing eQTLs between tissues. After verifying previous significance threshold-based estimates of tissue-specificity, we show their limitations given their dependency on statistical power. We propose that continuous estimates of the proportion of tissue-shared signals and direct comparison of the magnitude of effect on the fold change in expression are essential properties that jointly provide a biologically realistic view of tissue-specificity. Under this framework we demonstrate that 30% of eQTLs are shared among the three tissues studied, while another 29% appear exclusively tissue-specific. However, even among the shared eQTLs, a substantial proportion (10%–20%) have significant differences in the magnitude of fold change between genotypic classes across tissues. Our results underline the need to account for the complexity of eQTL tissue-specificity in an effort to assess consequences of such variants for complex traits.

Epigenetic modifications such as DNA methylation play a key role in gene regulation and disease susceptibility. However, little is known about the genome-wide frequency, localization, and function of methylation variation and how it is regulated by genetic and environmental factors. We utilized the Multiple Tissue Human Expression Resource (MuTHER) and generated Illumina 450K adipose methylome data from 648 twins. We found that individual CpGs had low variance and that variability was suppressed in promoters. We noted that DNA methylation variation was highly heritable (h(2)median = 0.34) and that shared environmental effects correlated with metabolic phenotype-associated CpGs. Analysis of methylation quantitative-trait loci (metQTL) revealed that 28% of CpGs were associated with nearby SNPs, and when overlapping them with adipose expression quantitative-trait loci (eQTL) from the same individuals, we found that 6% of the loci played a role in regulating both gene expression and DNA methylation. These associations were bidirectional, but there were pronounced negative associations for promoter CpGs. Integration of metQTL with adipose reference epigenomes and disease associations revealed significant enrichment of metQTL overlapping metabolic-trait or disease loci in enhancers (the strongest effects were for high-density lipoprotein cholesterol and body mass index [BMI]). We followed up with the BMI SNP rs713586, a cg01884057 metQTL that overlaps an enhancer upstream of ADCY3, and used bisulphite sequencing to refine this region. Our results showed widespread population invariability yet sequence dependence on adipose DNA methylation but that incorporating maps of regulatory elements aid in linking CpG variation to gene regulation and disease risk in a tissue-dependent manner.

The intersection of genome-wide association analyses with physiological and functional data indicates that variants regulating islet gene transcription influence type 2 diabetes (T2D) predisposition and glucose homeostasis. However, the specific genes through which these regulatory variants act remain poorly characterized. We generated expression quantitative trait locus (eQTL) data in 118 human islet samples using RNA-sequencing and high-density genotyping. We identified fourteen loci at which cis-exon-eQTL signals overlapped active islet chromatin signatures and were coincident with established T2D and/or glycemic trait associations. ‎At some, these data provide an experimental link between GWAS signals and biological candidates, such as DGKB and ADCY5. At others, the cis-signals implicate genes with no prior connection to islet biology, including WARS and ZMIZ1. At the ZMIZ1 locus, we show that perturbation of ZMIZ1 expression in human islets and beta-cells influences exocytosis and insulin secretion, highlighting a novel role for ZMIZ1 in the maintenance of glucose homeostasis. Together, these findings provide a significant advance in the mechanistic insights of T2D and glycemic trait association loci.

Aims/hypothesisMicroRNAs regulate a broad range of biological mechanisms. To investigate the relationship between microRNA expression and type 2 diabetes, we compared global microRNA expression in insulin target tissues from three inbred rat strains that differ in diabetes susceptibility.MethodsUsing microarrays, we measured the expression of 283 microRNAs in adipose, liver and muscle tissue from hyperglycaemic (Goto–Kakizaki), intermediate glycaemic (Wistar Kyoto) and normoglycaemic (Brown Norway) rats (n = 5 for each strain). Expression was compared across strains and validated using quantitative RT-PCR. Furthermore, microRNA expression variation in adipose tissue was investigated in 3T3-L1 adipocytes exposed to hyperglycaemic conditions.ResultsWe found 29 significantly differentiated microRNAs (padjusted < 0.05): nine in adipose tissue, 18 in liver and two in muscle. Of these, five microRNAs had expression patterns that correlated with the strain-specific glycaemic phenotype. MiR-222 (padjusted = 0.0005) and miR-27a (padjusted = 0.006) were upregulated in adipose tissue; miR-195 (padjusted = 0.006) and miR-103 (padjusted = 0.04) were upregulated in liver; and miR-10b (padjusted = 0.004) was downregulated in muscle. Exposure of 3T3-L1 adipocytes to increased glucose concentration upregulated the expression of miR-222 (p = 0.008), miR-27a (p = 0.02) and the previously reported miR-29a (p = 0.02). Predicted target genes of these differentially expressed microRNAs are involved in pathways relevant to type 2 diabetes.ConclusionThe expression patterns of miR-222, miR-27a, miR-195, miR-103 and miR-10b varied with hyperglycaemia, suggesting a role for these microRNAs in the pathophysiology of type 2 diabetes, as modelled by the Gyoto–Kakizaki rat. We observed similar patterns of expression of miR-222, miR-27a and miR-29a in adipocytes as a response to increased glucose levels, which supports our hypothesis that altered expression of microRNAs accompanies primary events related to the pathogenesis of type 2 diabetes.Electronic supplementary materialThe online version of this article (doi:10.1007/s00125-010-1667-2) contains supplementary material, which is available to authorised users.

The majority of genetic risk variants for type 2 diabetes (T2D) affect insulin secretion, but the mechanisms through which they influence pancreatic islet function remain largely unknown. We functionally characterized human islets to determine secretory, biophysical, and ultrastructural features in relation to genetic risk profiles in diabetic and nondiabetic donors. Islets from donors with T2D exhibited impaired insulin secretion, which was more pronounced in lean than obese diabetic donors. We assessed the impact of 14 disease susceptibility variants on measures of glucose sensing, exocytosis, and structure. Variants near TCF7L2 and ADRA2A were associated with reduced glucose-induced insulin secretion, whereas susceptibility variants near ADRA2A, KCNJ11, KCNQ1, and TCF7L2 were associated with reduced depolarization-evoked insulin exocytosis. KCNQ1, ADRA2A, KCNJ11, HHEX/IDE, and SLC2A2 variants affected granule docking. We combined our results to create a novel genetic risk score for β-cell dysfunction that includes aberrant granule docking, decreased Ca2+ sensitivity of exocytosis, and reduced insulin release. Individuals with a high risk score displayed an impaired response to intravenous glucose and deteriorating insulin secretion over time. Our results underscore the importance of defects in β-cell exocytosis in T2D and demonstrate the potential of cellular phenotypic characterization in the elucidation of complex genetic disorders.