Search

Objective: To identify genetic risk factors for incident cardiovascular disease (CVD) among people with type 2 diabetes (T2D)., Research Design and Methods: We conducted a multiancestry time-to-event genome-wide association study for incident CVD among people with T2D. We also tested 204 known coronary artery disease (CAD) variants for association with incident CVD., Results: Among 49,230 participants with T2D, 8,956 had incident CVD events (event rate 18.2%). We identified three novel genetic loci for incident CVD: rs147138607 (near CACNA1E/ZNF648, hazard ratio [HR] 1.23, P = 3.6 × 10-9), rs77142250 (near HS3ST1, HR 1.89, P = 9.9 × 10-9), and rs335407 (near TFB1M/NOX3, HR 1.25, P = 1.5 × 10-8). Among 204 known CAD loci, 5 were associated with incident CVD in T2D (multiple comparison-adjusted P < 0.00024, 0.05/204). A standardized polygenic score of these 204 variants was associated with incident CVD with HR 1.14 (P = 1.0 × 10-16)., Conclusions: The data point to novel and known genomic regions associated with incident CVD among individuals with T2D., (© 2024 by the American Diabetes Association.)

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes 1,2 and molecular mechanisms that are often specific to cell type 3,4 . Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P < 5 × 10 -8 ) that map to 611 loci, of which 145 loci are, to our knowledge, previously unreported. We define eight non-overlapping clusters of T2D signals that are characterized by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type-specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial cells and enteroendocrine cells. We build cluster-specific partitioned polygenic scores 5 in a further 279,552 individuals of diverse ancestry, including 30,288 cases of T2D, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned polygenic scores are associated with coronary artery disease, peripheral artery disease and end-stage diabetic nephropathy across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings show the value of integrating multi-ancestry genome-wide association study data with single-cell epigenomics to disentangle the aetiological heterogeneity that drives the development and progression of T2D. This might offer a route to optimize global access to genetically informed diabetes care., (© 2024. The Author(s).)

Objective: Clonal hematopoiesis of indeterminate potential (CHIP) is an aging-related accumulation of somatic mutations in hematopoietic stem cells, leading to clonal expansion. CHIP presence has been implicated in atherosclerotic coronary heart disease (CHD) and all-cause mortality, but its association with incident type 2 diabetes (T2D) is unknown. We hypothesized that CHIP is associated with elevated risk of T2D., Research Design and Methods: CHIP was derived from whole-genome sequencing of blood DNA in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine (TOPMed) prospective cohorts. We performed analysis for 17,637 participants from six cohorts, without prior T2D, cardiovascular disease, or cancer. We evaluated baseline CHIP versus no CHIP prevalence with incident T2D, including associations with DNMT3A, TET2, ASXL1, JAK2, and TP53 variants. We estimated multivariable-adjusted hazard ratios (HRs) and 95% CIs with adjustment for age, sex, BMI, smoking, alcohol, education, self-reported race/ethnicity, and combined cohorts' estimates via fixed-effects meta-analysis., Results: Mean (SD) age was 63.4 (11.5) years, 76% were female, and CHIP prevalence was 6.0% (n = 1,055) at baseline. T2D was diagnosed in n = 2,467 over mean follow-up of 9.8 years. Participants with CHIP had 23% (CI 1.04, 1.45) higher risk of T2D than those with no CHIP. Specifically, higher risk was for TET2 (HR 1.48; CI 1.05, 2.08) and ASXL1 (HR 1.76; CI 1.03, 2.99) mutations; DNMT3A was nonsignificant (HR 1.15; CI 0.93, 1.43). Statistical power was limited for JAK2 and TP53 analyses., Conclusions: CHIP was associated with higher incidence of T2D. CHIP mutations located on genes implicated in CHD and mortality were also related to T2D, suggesting shared aging-related pathology., (© 2023 by the American Diabetes Association.)

Background: Insulin resistance (IR) is a major risk factor for Alzheimer's disease (AD) dementia. The mechanisms by which IR predisposes to AD are not well-understood. Epigenetic studies may help identify molecular signatures of IR associated with AD, thus improving our understanding of the biological and regulatory mechanisms linking IR and AD., Methods: We conducted an epigenome-wide association study of IR, quantified using the homeostatic model assessment of IR (HOMA-IR) and adjusted for body mass index, in 3,167 participants from the Framingham Heart Study (FHS) without type 2 diabetes at the time of blood draw used for methylation measurement. We identified DNA methylation markers associated with IR at the genome-wide level accounting for multiple testing (P < 1.1 × 10 -7 ) and evaluated their association with neurological traits in participants from the FHS (N = 3040) and the Religious Orders Study/Memory and Aging Project (ROSMAP, N = 707). DNA methylation profiles were measured in blood (FHS) or dorsolateral prefrontal cortex (ROSMAP) using the Illumina HumanMethylation450 BeadChip. Linear regressions (ROSMAP) or mixed-effects models accounting for familial relatedness (FHS) adjusted for age, sex, cohort, self-reported race, batch, and cell type proportions were used to assess associations between DNA methylation and neurological traits accounting for multiple testing., Results: We confirmed the strong association of blood DNA methylation with IR at three loci (cg17901584-DHCR24, cg17058475-CPT1A, cg00574958-CPT1A, and cg06500161-ABCG1). In FHS, higher levels of blood DNA methylation at cg00574958 and cg17058475 were both associated with lower IR (P = 2.4 × 10 -11 and P = 9.0 × 10 -8 ), larger total brain volumes (P = 0.03 and P = 9.7 × 10 -4 ), and smaller log lateral ventricular volumes (P = 0.07 and P = 0.03). In ROSMAP, higher levels of brain DNA methylation at the same two CPT1A markers were associated with greater risk of cognitive impairment (P = 0.005 and P = 0.02) and higher AD-related indices (CERAD score: P = 5 × 10 -4 and 0.001; Braak stage: P = 0.004 and P = 0.01)., Conclusions: Our results suggest potentially distinct epigenetic regulatory mechanisms between peripheral blood and dorsolateral prefrontal cortex tissues underlying IR and AD at CPT1A locus., (© 2023. The Author(s).)

Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on 'around the clock' glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor (GLP1R), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification., (© 2023. The Author(s).)

Aims/hypothesis: The Latino population has been systematically underrepresented in large-scale genetic analyses, and previous studies have relied on the imputation of ungenotyped variants based on the 1000 Genomes (1000G) imputation panel, which results in suboptimal capture of low-frequency or Latino-enriched variants. The National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) released the largest multi-ancestry genotype reference panel representing a unique opportunity to analyse rare genetic variations in the Latino population. We hypothesise that a more comprehensive analysis of low/rare variation using the TOPMed panel would improve our knowledge of the genetics of type 2 diabetes in the Latino population., Methods: We evaluated the TOPMed imputation performance using genotyping array and whole-exome sequence data in six Latino cohorts. To evaluate the ability of TOPMed imputation to increase the number of identified loci, we performed a Latino type 2 diabetes genome-wide association study (GWAS) meta-analysis in 8150 individuals with type 2 diabetes and 10,735 control individuals and replicated the results in six additional cohorts including whole-genome sequence data from the All of Us cohort., Results: Compared with imputation with 1000G, the TOPMed panel improved the identification of rare and low-frequency variants. We identified 26 genome-wide significant signals including a novel variant (minor allele frequency 1.7%; OR 1.37, p=3.4 × 10 -9 ). A Latino-tailored polygenic score constructed from our data and GWAS data from East Asian and European populations improved the prediction accuracy in a Latino target dataset, explaining up to 7.6% of the type 2 diabetes risk variance., Conclusions/interpretation: Our results demonstrate the utility of TOPMed imputation for identifying low-frequency variants in understudied populations, leading to the discovery of novel disease associations and the improvement of polygenic scores., Data Availability: Full summary statistics are available through the Common Metabolic Diseases Knowledge Portal ( https://t2d.hugeamp.org/downloads.html ) and through the GWAS catalog ( https://www.ebi.ac.uk/gwas/ , accession ID: GCST90255648). Polygenic score (PS) weights for each ancestry are available via the PGS catalog ( https://www.pgscatalog.org , publication ID: PGP000445, scores IDs: PGS003443, PGS003444 and PGS003445)., (© 2023. The Author(s).)

Distinct tissue-specific mechanisms mediate insulin action in fasting and postprandial states. Previous genetic studies have largely focused on insulin resistance in the fasting state, where hepatic insulin action dominates. Here we studied genetic variants influencing insulin levels measured 2 h after a glucose challenge in >55,000 participants from three ancestry groups. We identified ten new loci (P < 5 × 10 -8 ) not previously associated with postchallenge insulin resistance, eight of which were shown to share their genetic architecture with type 2 diabetes in colocalization analyses. We investigated candidate genes at a subset of associated loci in cultured cells and identified nine candidate genes newly implicated in the expression or trafficking of GLUT4, the key glucose transporter in postprandial glucose uptake in muscle and fat. By focusing on postprandial insulin resistance, we highlighted the mechanisms of action at type 2 diabetes loci that are not adequately captured by studies of fasting glycemic traits., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Background: Whether genetics contribute to the rising prevalence of obesity or its cardiovascular consequences in today's obesogenic environment remains unclear. We sought to determine whether the effects of a higher aggregate genetic burden of obesity risk on body mass index (BMI) or cardiovascular disease (CVD) differed by birth year., Methods: We split the FHS (Framingham Heart Study) into 4 equally sized birth cohorts (birth year before 1932, 1932 to 1946, 1947 to 1959, and after 1960). We modeled a genetic predisposition to obesity using an additive genetic risk score (GRS) of 941 BMI-associated variants and tested for GRS-birth year interaction on log-BMI (outcome) when participants were around 50 years old (N=7693). We repeated the analysis using a GRS of 109 BMI-associated variants that increased CVD risk factors (type 2 diabetes, blood pressure, total cholesterol, and high-density lipoprotein) in addition to BMI. We then evaluated whether the effects of the BMI GRSs on CVD risk differed by birth cohort when participants were around 60 years old (N=5493)., Results: Compared with participants born before 1932 (mean age, 50.8 yrs [2.4]), those born after 1960 (mean age, 43.3 years [4.5]) had higher BMI (median, 25.4 [23.3-28.0] kg/m 2 versus 26.9 [interquartile range, 23.7-30.6] kg/m 2 ). The effect of the 941-variant BMI GRS on BMI and CVD risk was stronger in people who were born in later years (GRS-birth year interaction: P =0.0007 and P =0.04 respectively)., Conclusions: The significant GRS-birth year interactions indicate that common genetic variants have larger effects on middle-age BMI and CVD risk in people born more recently. These findings suggest that the increasingly obesogenic environment may amplify the impact of genetics on the risk of obesity and possibly its cardiovascular consequences.

Background: Considering relatives' health history in logistic regression for case-control genome-wide association studies (CC-GWAS) may provide new information that increases accuracy and power to detect disease associated genetic variants. We conducted simulations and analyzed type 2 diabetes (T2D) data from the Framingham Heart Study (FHS) to compare two methods, liability threshold model conditional on both case-control status and family history (LT-FH) and Fam-meta, which incorporate family history into CC-GWAS., Results: In our simulation scenario of trait with modest T2D heritability (h 2  = 0.28), variant minor allele frequency ranging from 1% to 50%, and 1% of phenotype variance explained by the genetic variants, Fam-meta had the highest overall power, while both methods incorporating family history were more powerful than CC-GWAS. All three methods had controlled type I error rates, while LT-FH was the most conservative with a lower-than-expected error rate. In addition, we observed a substantial increase in power of the two familial history methods compared to CC-GWAS when the prevalence of the phenotype increased with age. Furthermore, we showed that, when only the phenotypes of more distant relatives were available, Fam-meta still remained more powerful than CC-GWAS, confirming that leveraging disease history of both close and distant relatives can increase power of association analyses. Using FHS data, we confirmed the well-known association of TCF7L2 region with T2D at the genome-wide threshold of P-value < 5 × 10 -8 , and both familial history methods increased the significance of the region compared to CC-GWAS. We identified two loci at 5q35 (ADAMTS2) and 5q23 (PRR16), not previously reported for T2D using CC-GWAS and Fam-meta; both genes play a role in cardiovascular diseases. Additionally, CC-GWAS detected one more significant locus at 13q31 (GPC6) reported associated with T2D-related traits., Conclusions: Overall, LT-FH and Fam-meta had higher power than CC-GWAS in simulations, especially using phenotypes that were more prevalent in older age groups, and both methods detected known genetic variants with lower P-values in real data application, highlighting the benefits of including family history in genetic association studies., (© 2022. The Author(s).)

The genetic determinants of fasting glucose (FG) and fasting insulin (FI) have been studied mostly through genome arrays, resulting in over 100 associated variants. We extended this work with high-coverage whole genome sequencing analyses from fifteen cohorts in NHLBI's Trans-Omics for Precision Medicine (TOPMed) program. Over 23,000 non-diabetic individuals from five race-ethnicities/populations (African, Asian, European, Hispanic and Samoan) were included. Eight variants were significantly associated with FG or FI across previously identified regions MTNR1B, G6PC2, GCK, GCKR and FOXA2. We additionally characterize suggestive associations with FG or FI near previously identified SLC30A8, TCF7L2, and ADCY5 regions as well as APOB, PTPRT, and ROBO1. Functional annotation resources including the Diabetes Epigenome Atlas were compiled for each signal (chromatin states, annotation principal components, and others) to elucidate variant-to-function hypotheses. We provide a catalog of nucleotide-resolution genomic variation spanning intergenic and intronic regions creating a foundation for future sequencing-based investigations of glycemic traits., (© 2022. The Author(s).)

Objective: Tryptophan can be catabolised to various metabolites through host kynurenine and microbial indole pathways. We aimed to examine relationships of host and microbial tryptophan metabolites with incident type 2 diabetes (T2D), host genetics, diet and gut microbiota., Method: We analysed associations between circulating levels of 11 tryptophan metabolites and incident T2D in 9180 participants of diverse racial/ethnic backgrounds from five cohorts. We examined host genome-wide variants, dietary intake and gut microbiome associated with these metabolites., Results: Tryptophan, four kynurenine-pathway metabolites (kynurenine, kynurenate, xanthurenate and quinolinate) and indolelactate were positively associated with T2D risk, while indolepropionate was inversely associated with T2D risk. We identified multiple host genetic variants, dietary factors, gut bacteria and their potential interplay associated with these T2D-relaetd metabolites. Intakes of fibre-rich foods, but not protein/tryptophan-rich foods, were the dietary factors most strongly associated with tryptophan metabolites. The fibre-indolepropionate association was partially explained by indolepropionate-associated gut bacteria, mostly fibre-using Firmicutes . We identified a novel association between a host functional LCT variant (determining lactase persistence) and serum indolepropionate, which might be related to a host gene-diet interaction on gut Bifidobacterium , a probiotic bacterium significantly associated with indolepropionate independent of other fibre-related bacteria. Higher milk intake was associated with higher levels of gut Bifidobacterium and serum indolepropionate only among genetically lactase non-persistent individuals., Conclusion: Higher milk intake among lactase non-persistent individuals, and higher fibre intake were associated with a favourable profile of circulating tryptophan metabolites for T2D, potentially through the host-microbial cross-talk shifting tryptophan metabolism toward gut microbial indolepropionate production., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2022. No commercial re-use. See rights and permissions. Published by BMJ.)

We assembled an ancestrally diverse collection of genome-wide association studies (GWAS) of type 2 diabetes (T2D) in 180,834 affected individuals and 1,159,055 controls (48.9% non-European descent) through the Diabetes Meta-Analysis of Trans-Ethnic association studies (DIAMANTE) Consortium. Multi-ancestry GWAS meta-analysis identified 237 loci attaining stringent genome-wide significance (P < 5 × 10 -9 ), which were delineated to 338 distinct association signals. Fine-mapping of these signals was enhanced by the increased sample size and expanded population diversity of the multi-ancestry meta-analysis, which localized 54.4% of T2D associations to a single variant with >50% posterior probability. This improved fine-mapping enabled systematic assessment of candidate causal genes and molecular mechanisms through which T2D associations are mediated, laying the foundations for functional investigations. Multi-ancestry genetic risk scores enhanced transferability of T2D prediction across diverse populations. Our study provides a step toward more effective clinical translation of T2D GWAS to improve global health for all, irrespective of genetic background., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

A challenge in standard genetic studies is maintaining good power to detect associations, especially for low prevalent diseases and rare variants. The traditional methods are most powerful when evaluating the association between variants in balanced study designs. Without accounting for family correlation and unbalanced case-control ratio, these analyses could result in inflated type I error. One cost-effective solution to increase statistical power is exploitation of available family history (FH) that contains valuable information about disease heritability. Here, we develop methods to address the aforementioned type I error issues while providing optimal power to analyze aggregates of rare variants by incorporating additional information from FH. With enhanced power in these methods exploiting FH and accounting for relatedness and unbalanced designs, we successfully detect genes with suggestive associations with Alzheimer disease, dementia, and type 2 diabetes by using the exome chip data from the Framingham Heart Study., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Objective: Type 2 diabetes (T2D) has heterogeneous patient clinical characteristics and outcomes. In previous work, we investigated the genetic basis of this heterogeneity by clustering 94 T2D genetic loci using their associations with 47 diabetes-related traits and identified five clusters, termed β-cell, proinsulin, obesity, lipodystrophy, and liver/lipid. The relationship between these clusters and individual-level metabolic disease outcomes has not been assessed., Research Design and Methods: Here we constructed individual-level partitioned polygenic scores (pPS) for these five clusters in 12 studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (n = 454,193) and tested for cross-sectional association with T2D-related outcomes, including blood pressure, renal function, insulin use, age at T2D diagnosis, and coronary artery disease (CAD)., Results: Despite all clusters containing T2D risk-increasing alleles, they had differential associations with metabolic outcomes. Increased obesity and lipodystrophy cluster pPS, which had opposite directions of association with measures of adiposity, were both significantly associated with increased blood pressure and hypertension. The lipodystrophy and liver/lipid cluster pPS were each associated with CAD, with increasing and decreasing effects, respectively. An increased liver/lipid cluster pPS was also significantly associated with reduced renal function. The liver/lipid cluster includes known loci linked to liver lipid metabolism (e.g., GCKR, PNPLA3, and TM6SF2), and these findings suggest that cardiovascular disease risk and renal function may be impacted by these loci through their shared disease pathway., Conclusions: Our findings support that genetically driven pathways leading to T2D also predispose differentially to clinical outcomes., (© 2022 by the American Diabetes Association.)

Objective: LDL cholesterol (LDLc)-lowering drugs modestly increase body weight and type 2 diabetes risk, but the extent to which the diabetogenic effect of lowering LDLc is mediated through increased BMI is unknown., Research Design and Methods: We conducted summary-level univariable and multivariable Mendelian randomization (MR) analyses in 921,908 participants to investigate the effect of lowering LDLc on type 2 diabetes risk and the proportion of this effect mediated through BMI. We used data from 92,532 participants from 14 observational studies to replicate findings in individual-level MR analyses., Results: A 1-SD decrease in genetically predicted LDLc was associated with increased type 2 diabetes odds (odds ratio [OR] 1.12 [95% CI 1.01, 1.24]) and BMI (β = 0.07 SD units [95% CI 0.02, 0.12]) in univariable MR analyses. The multivariable MR analysis showed evidence of an indirect effect of lowering LDLc on type 2 diabetes through BMI (OR 1.04 [95% CI 1.01, 1.08]) with a proportion mediated of 38% of the total effect (P = 0.03). Total and indirect effect estimates were similar across a number of sensitivity analyses. Individual-level MR analyses confirmed the indirect effect of lowering LDLc on type 2 diabetes through BMI with an estimated proportion mediated of 8% (P = 0.04)., Conclusions: These findings suggest that the diabetogenic effect attributed to lowering LDLc is partially mediated through increased BMI. Our results could help advance understanding of adipose tissue and lipids in type 2 diabetes pathophysiology and inform strategies to reduce diabetes risk among individuals taking LDLc-lowering medications., (© 2021 by the American Diabetes Association.)

Hundreds of thousands of genetic variants have been reported to cause severe monogenic diseases, but the probability that a variant carrier develops the disease (termed penetrance) is unknown for virtually all of them. Additionally, the clinical utility of common polygenetic variation remains uncertain. Using exome sequencing from 77,184 adult individuals (38,618 multi-ancestral individuals from a type 2 diabetes case-control study and 38,566 participants from the UK Biobank, for whom genotype array data were also available), we apply clinical standard-of-care gene variant curation for eight monogenic metabolic conditions. Rare variants causing monogenic diabetes and dyslipidemias display effect sizes significantly larger than the top 1% of the corresponding polygenic scores. Nevertheless, penetrance estimates for monogenic variant carriers average 60% or lower for most conditions. We assess epidemiologic and genetic factors contributing to risk prediction in monogenic variant carriers, demonstrating that inclusion of polygenic variation significantly improves biomarker estimation for two monogenic dyslipidemias.

Smoking is a potentially causal behavioral risk factor for type 2 diabetes (T2D), but not all smokers develop T2D. It is unknown whether genetic factors partially explain this variation. We performed genome-environment-wide interaction studies to identify loci exhibiting potential interaction with baseline smoking status (ever vs. never) on incident T2D and fasting glucose (FG). Analyses were performed in participants of European (EA) and African ancestry (AA) separately. Discovery analyses were conducted using genotype data from the 50,000-single-nucleotide polymorphism (SNP) ITMAT-Broad-CARe (IBC) array in 5 cohorts from from the Candidate Gene Association Resource Consortium (n = 23,189). Replication was performed in up to 16 studies from the Cohorts for Heart Aging Research in Genomic Epidemiology Consortium (n = 74,584). In meta-analysis of discovery and replication estimates, 5 SNPs met at least one criterion for potential interaction with smoking on incident T2D at p<1x10-7 (adjusted for multiple hypothesis-testing with the IBC array). Two SNPs had significant joint effects in the overall model and significant main effects only in one smoking stratum: rs140637 (FBN1) in AA individuals had a significant main effect only among smokers, and rs1444261 (closest gene C2orf63) in EA individuals had a significant main effect only among nonsmokers. Three additional SNPs were identified as having potential interaction by exhibiting a significant main effects only in smokers: rs1801232 (CUBN) in AA individuals, rs12243326 (TCF7L2) in EA individuals, and rs4132670 (TCF7L2) in EA individuals. No SNP met significance for potential interaction with smoking on baseline FG. The identification of these loci provides evidence for genetic interactions with smoking exposure that may explain some of the heterogeneity in the association between smoking and T2D., Competing Interests: I have read the journal's policy and the authors of this manuscript have the following competing interests. Dr. Meigs currently has a research grant from GlaxoSmithKline and serves on a consultancy board for Interleukin Genetics. Dr. Florez has received consulting honoraria from Daiichi-Sankyo and AstraZeneca. Dr. Mike A. Nalls is supported by a consulting contract between Data Tecnica International LLC and the National Institute on Aging (NIA), National Institutes of Health (NIH), Bethesda, MD, USA. Dr. Nalls also consults for Illumina Inc., the Michael J. Fox Foundation, and the University of California Healthcare. DR. Jose C. Florez, Consulting honoraria from Janssen Pharmaceuticals and Goldfinch Bio, and speaking honorarium from Novo Nordisk The other authors declare no conflicts of interest. This does not alter our adherence to PLOS ONE policies on sharing data and materials. There are no patents, products in development or marketed products associated with this research to declare.

Epigenetic changes may contribute substantially to risks of diseases of aging. Previous studies reported seven methylation variable positions (MVPs) robustly associated with incident type 2 diabetes mellitus (T2DM). However, their causal roles in T2DM are unclear. In an incident T2DM case-cohort study nested within the population-based European Prospective Investigation into Cancer and Nutrition (EPIC)-Norfolk cohort, we used whole blood DNA collected at baseline, up to 11 years before T2DM onset, to investigate the role of methylation in the etiology of T2DM. We identified 15 novel MVPs with robust associations with incident T2DM and robustly confirmed three MVPs identified previously (near to TXNIP , ABCG1 , and SREBF1 ). All 18 MVPs showed directionally consistent associations with incident and prevalent T2DM in independent studies. Further conditional analyses suggested that the identified epigenetic signals appear related to T2DM via glucose and obesity-related pathways acting before the collection of baseline samples. We integrated genome-wide genetic data to identify methylation-associated quantitative trait loci robustly associated with 16 of the 18 MVPs and found one MVP, cg00574958 at CPT1A , with a possible direct causal role in T2DM. None of the implicated genes were previously highlighted by genetic association studies, suggesting that DNA methylation studies may reveal novel biological mechanisms involved in tissue responses to glycemia., (© 2019 by the American Diabetes Association.)

Scope: Insulin resistance (IR) and inflammation are hallmarks of type 2 diabetes (T2D). The nod-like receptor pyrin domain containing-3 (NLRP3) inflammasome is a metabolic sensor activated by saturated fatty acids (SFA) initiating IL-1β inflammation and IR. Interactions between SFA intake and NLRP3-related genetic variants may alter T2D risk factors., Methods: Meta-analyses of six Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (n = 19 005) tested interactions between SFA and NLRP3-related single-nucleotide polymorphisms (SNPs) and modulation of fasting insulin, fasting glucose, and homeostasis model assessment of insulin resistance., Results: SFA interacted with rs12143966, wherein each 1% increase in SFA intake increased insulin by 0.0063 IU mL -1 (SE ± 0.002, p = 0.001) per each major (G) allele copy. rs4925663, interacted with SFA (β ± SE = -0.0058 ± 0.002, p = 0.004) to increase insulin by 0.0058 IU mL -1 , per additional copy of the major (C) allele. Both associations are close to the significance threshold (p < 0.0001). rs4925663 causes a missense mutation affecting NLRP3 expression., Conclusion: Two NLRP3-related SNPs showed potential interaction with SFA to modulate fasting insulin. Greater dietary SFA intake accentuates T2D risk, which, subject to functional validation, may be further elaborated depending on NLRP3-related genetic variants., (© 2019 The Authors. Published by WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Despite existing reports on differential DNA methylation in type 2 diabetes (T2D) and obesity, our understanding of its functional relevance remains limited. Here we show the effect of differential methylation in the early phases of T2D pathology by a blood-based epigenome-wide association study of 4808 non-diabetic Europeans in the discovery phase and 11,750 individuals in the replication. We identify CpGs in LETM1, RBM20, IRS2, MAN2A2 and the 1q25.3 region associated with fasting insulin, and in FCRL6, SLAMF1, APOBEC3H and the 15q26.1 region with fasting glucose. In silico cross-omics analyses highlight the role of differential methylation in the crosstalk between the adaptive immune system and glucose homeostasis. The differential methylation explains at least 16.9% of the association between obesity and insulin. Our study sheds light on the biological interactions between genetic variants driving differential methylation and gene expression in the early pathogenesis of T2D.

Protein-coding genetic variants that strongly affect disease risk can yield relevant clues to disease pathogenesis. Here we report exome-sequencing analyses of 20,791 individuals with type 2 diabetes (T2D) and 24,440 non-diabetic control participants from 5 ancestries. We identify gene-level associations of rare variants (with minor allele frequencies of less than 0.5%) in 4 genes at exome-wide significance, including a series of more than 30 SLC30A8 alleles that conveys protection against T2D, and in 12 gene sets, including those corresponding to T2D drug targets (P = 6.1 × 10 -3 ) and candidate genes from knockout mice (P = 5.2 × 10 -3 ). Within our study, the strongest T2D gene-level signals for rare variants explain at most 25% of the heritability of the strongest common single-variant signals, and the gene-level effect sizes of the rare variants that we observed in established T2D drug targets will require 75,000-185,000 sequenced cases to achieve exome-wide significance. We propose a method to interpret these modest rare-variant associations and to incorporate these associations into future target or gene prioritization efforts.

We expanded GWAS discovery for type 2 diabetes (T2D) by combining data from 898,130 European-descent individuals (9% cases), after imputation to high-density reference panels. With these data, we (i) extend the inventory of T2D-risk variants (243 loci, 135 newly implicated in T2D predisposition, comprising 403 distinct association signals); (ii) enrich discovery of lower-frequency risk alleles (80 index variants with minor allele frequency <5%, 14 with estimated allelic odds ratio >2); (iii) substantially improve fine-mapping of causal variants (at 51 signals, one variant accounted for >80% posterior probability of association (PPA)); (iv) extend fine-mapping through integration of tissue-specific epigenomic information (islet regulatory annotations extend the number of variants with PPA >80% to 73); (v) highlight validated therapeutic targets (18 genes with associations attributable to coding variants); and (vi) demonstrate enhanced potential for clinical translation (genome-wide chip heritability explains 18% of T2D risk; individuals in the extremes of a T2D polygenic risk score differ more than ninefold in prevalence).

Aims/hypothesis: Identifying the metabolite profile of individuals with normal fasting glucose (NFG [<5.55 mmol/l]) who progressed to type 2 diabetes may give novel insights into early type 2 diabetes disease interception and detection., Methods: We conducted a population-based prospective study among 1150 Framingham Heart Study Offspring cohort participants, age 40-65 years, with NFG. Plasma metabolites were profiled by LC-MS/MS. Penalised regression models were used to select measured metabolites for type 2 diabetes incidence classification (training dataset) and to internally validate the discriminatory capability of selected metabolites beyond conventional type 2 diabetes risk factors (testing dataset)., Results: Over a follow-up period of 20 years, 95 individuals with NFG developed type 2 diabetes. Nineteen metabolites were selected repeatedly in the training dataset for type 2 diabetes incidence classification and were found to improve type 2 diabetes risk prediction beyond conventional type 2 diabetes risk factors (AUC was 0.81 for risk factors vs 0.90 for risk factors + metabolites, p = 1.1 × 10 -4 ). Using pathway enrichment analysis, the nitrogen metabolism pathway, which includes three prioritised metabolites (glycine, taurine and phenylalanine), was significantly enriched for association with type 2 diabetes risk at the false discovery rate of 5% (p = 0.047). In adjusted Cox proportional hazard models, the type 2 diabetes risk per 1 SD increase in glycine, taurine and phenylalanine was 0.65 (95% CI 0.54, 0.78), 0.73 (95% CI 0.59, 0.9) and 1.35 (95% CI 1.11, 1.65), respectively. Mendelian randomisation demonstrated a similar relationship for type 2 diabetes risk per 1 SD genetically increased glycine (OR 0.89 [95% CI 0.8, 0.99]) and phenylalanine (OR 1.6 [95% CI 1.08, 2.4])., Conclusions/interpretation: In individuals with NFG, information from a discrete set of 19 metabolites improved prediction of type 2 diabetes beyond conventional risk factors. In addition, the nitrogen metabolism pathway and its components emerged as a potential effector of earliest stages of type 2 diabetes pathophysiology.

We aggregated coding variant data for 81,412 type 2 diabetes cases and 370,832 controls of diverse ancestry, identifying 40 coding variant association signals (P < 2.2 × 10 -7 ); of these, 16 map outside known risk-associated loci. We make two important observations. First, only five of these signals are driven by low-frequency variants: even for these, effect sizes are modest (odds ratio ≤1.29). Second, when we used large-scale genome-wide association data to fine-map the associated variants in their regional context, accounting for the global enrichment of complex trait associations in coding sequence, compelling evidence for coding variant causality was obtained for only 16 signals. At 13 others, the associated coding variants clearly represent 'false leads' with potential to generate erroneous mechanistic inference. Coding variant associations offer a direct route to biological insight for complex diseases and identification of validated therapeutic targets; however, appropriate mechanistic inference requires careful specification of their causal contribution to disease predisposition.

A major challenge in evaluating the contribution of rare variants to complex disease is identifying enough copies of the rare alleles to permit informative statistical analysis. To investigate the contribution of rare variants to the risk of type 2 diabetes (T2D) and related traits, we performed deep whole-genome analysis of 1,034 members of 20 large Mexican-American families with high prevalence of T2D. If rare variants of large effect accounted for much of the diabetes risk in these families, our experiment was powered to detect association. Using gene expression data on 21,677 transcripts for 643 pedigree members, we identified evidence for large-effect rare-variant cis -expression quantitative trait loci that could not be detected in population studies, validating our approach. However, we did not identify any rare variants of large effect associated with T2D, or the related traits of fasting glucose and insulin, suggesting that large-effect rare variants account for only a modest fraction of the genetic risk of these traits in this sample of families. Reliable identification of large-effect rare variants will require larger samples of extended pedigrees or different study designs that further enrich for such variants., Competing Interests: Conflict of interest statement: S.E.L., J. Laramie, and R.G.T. were employees of Complete Genomics during this study. T.M.T. is an employee of Regeneron Pharmaceuticals. D.A. is an employee of Vertex Pharmaceuticals.

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

To characterize type 2 diabetes (T2D)-associated variation across the allele frequency spectrum, we conducted a meta-analysis of genome-wide association data from 26,676 T2D case and 132,532 control subjects of European ancestry after imputation using the 1000 Genomes multiethnic reference panel. Promising association signals were followed up in additional data sets (of 14,545 or 7,397 T2D case and 38,994 or 71,604 control subjects). We identified 13 novel T2D-associated loci ( P < 5 × 10 -8 ), including variants near the GLP2R , GIP , and HLA-DQA1 genes. Our analysis brought the total number of independent T2D associations to 128 distinct signals at 113 loci. Despite substantially increased sample size and more complete coverage of low-frequency variation, all novel associations were driven by common single nucleotide variants. Credible sets of potentially causal variants were generally larger than those based on imputation with earlier reference panels, consistent with resolution of causal signals to common risk haplotypes. Stratification of T2D-associated loci based on T2D-related quantitative trait associations revealed tissue-specific enrichment of regulatory annotations in pancreatic islet enhancers for loci influencing insulin secretion and in adipocytes, monocytes, and hepatocytes for insulin action-associated loci. These findings highlight the predominant role played by common variants of modest effect and the diversity of biological mechanisms influencing T2D pathophysiology., (© 2017 by the American Diabetes Association.)

Background: Glycated hemoglobin (HbA1c) is used to diagnose type 2 diabetes (T2D) and assess glycemic control in patients with diabetes. Previous genome-wide association studies (GWAS) have identified 18 HbA1c-associated genetic variants. These variants proved to be classifiable by their likely biological action as erythrocytic (also associated with erythrocyte traits) or glycemic (associated with other glucose-related traits). In this study, we tested the hypotheses that, in a very large scale GWAS, we would identify more genetic variants associated with HbA1c and that HbA1c variants implicated in erythrocytic biology would affect the diagnostic accuracy of HbA1c. We therefore expanded the number of HbA1c-associated loci and tested the effect of genetic risk-scores comprised of erythrocytic or glycemic variants on incident diabetes prediction and on prevalent diabetes screening performance. Throughout this multiancestry study, we kept a focus on interancestry differences in HbA1c genetics performance that might influence race-ancestry differences in health outcomes., Methods & Findings: Using genome-wide association meta-analyses in up to 159,940 individuals from 82 cohorts of European, African, East Asian, and South Asian ancestry, we identified 60 common genetic variants associated with HbA1c. We classified variants as implicated in glycemic, erythrocytic, or unclassified biology and tested whether additive genetic scores of erythrocytic variants (GS-E) or glycemic variants (GS-G) were associated with higher T2D incidence in multiethnic longitudinal cohorts (N = 33,241). Nineteen glycemic and 22 erythrocytic variants were associated with HbA1c at genome-wide significance. GS-G was associated with higher T2D risk (incidence OR = 1.05, 95% CI 1.04-1.06, per HbA1c-raising allele, p = 3 × 10-29); whereas GS-E was not (OR = 1.00, 95% CI 0.99-1.01, p = 0.60). In Europeans and Asians, erythrocytic variants in aggregate had only modest effects on the diagnostic accuracy of HbA1c. Yet, in African Americans, the X-linked G6PD G202A variant (T-allele frequency 11%) was associated with an absolute decrease in HbA1c of 0.81%-units (95% CI 0.66-0.96) per allele in hemizygous men, and 0.68%-units (95% CI 0.38-0.97) in homozygous women. The G6PD variant may cause approximately 2% (N = 0.65 million, 95% CI 0.55-0.74) of African American adults with T2D to remain undiagnosed when screened with HbA1c. Limitations include the smaller sample sizes for non-European ancestries and the inability to classify approximately one-third of the variants. Further studies in large multiethnic cohorts with HbA1c, glycemic, and erythrocytic traits are required to better determine the biological action of the unclassified variants., Conclusions: As G6PD deficiency can be clinically silent until illness strikes, we recommend investigation of the possible benefits of screening for the G6PD genotype along with using HbA1c to diagnose T2D in populations of African ancestry or groups where G6PD deficiency is common. Screening with direct glucose measurements, or genetically-informed HbA1c diagnostic thresholds in people with G6PD deficiency, may be required to avoid missed or delayed diagnoses.

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 ., (© 2017 by the American Diabetes Association.)

Objective: This study tested the hypothesis that genetically raised hyperglycemia increases coronary artery disease (CAD) risk separately from the risk conferred by type 2 diabetes as a whole., Research Design and Methods: We conducted a Mendelian randomization (MR) analysis using summary-level statistics from the largest published meta-analyses of genome-wide association studies (GWAS) for fasting glucose (FG) ( n = 133,010 participants free of diabetes) and CAD ( n = 63,746 case subjects and 130,681 control subjects) of predominantly European ancestry. FG-increasing variants associated with type 2 diabetes from the largest GWAS for type 2 diabetes were excluded. Variants with pleiotropic effects on other CAD risk factors (blood lipids, blood pressure, and obesity) were excluded using summary-level data from the largest published GWAS. Data from the Framingham Heart Study were used to validate the MR instrument and to build an FG genetic risk score (GRS)., Results: In an instrumental variable analysis comprising 12 FG-raising variants, a 1 mmol/L increase in FG revealed an effect-size estimate of 1.43 CAD odds (95% CI 1.14-1.79). The association was preserved after excluding variants for heterogeneity and pleiotropic effects on other CAD risk factors (odds ratio [OR] 1.33 [95% CI 1.02-1.73]). The 12 FG-increasing variants did not significantly increase type 2 diabetes risk (OR 1.05 [95% CI 0.91-1.23]), and its prevalence was constant across FG GRS quintiles ( P = 0.72)., Conclusions: Our data support that genetic predisposition to hyperglycemia raises the odds of CAD separately from type 2 diabetes and other CAD risk factors. These findings suggest that modulating glycemia may provide cardiovascular benefit., (© 2017 by the American Diabetes Association.)

Background: Obesity and diabetes family history are the two strongest risk factors for type 2 diabetes (T2D). Prior work shows that an individual's obesity risk is associated with obesity in social contacts, but whether T2D risk follows similar patterns is unknown., Objective: We aimed to estimate the relationship between obesity or diabetes in an individual's social contacts and his/her T2D risk. We hypothesized that obesity and diabetes in social contacts would increase an individual's T2D risk., Design: This was a retrospective analysis of the community-based Framingham Offspring Study (FOS)., Participants: FOS participants with T2D status, height and weight, and at least one social contact were eligible for this study (n = 4797 at Exam 1). Participants' interpersonal ties, cardiometabolic and demographic variables were available at eight exams from 1971 to 2008, and a T2D additive polygenic risk score was measured at the fifth exam., Main Measures: Primary exposures were T2D (fasting glucose ≥ 7 mmol/L or taking diabetes medications) and obesity status (BMI ≥ 30 kg/m(2)) of social contacts at a prior exam. Primary outcome was incident T2D in participants., Key Results: Incident T2D was associated with having a social contact with diabetes (OR 1.32, p = 0.004) or with obesity (OR 1.21, p = 0.004). In stratified analyses, incident T2D was associated with diabetes in siblings (OR 1.64, p = 0.001) and obesity in spouses (OR 1.54, p = 0.0004). The associations between diabetes and obesity in social contacts and an individual's incident diabetes risk were stronger in individuals with a high diabetes genetic risk score., Conclusions: T2D and obesity in social contacts, particularly siblings and spouses, were associated with an individual's risk of incident diabetes even after accounting for parental T2D history. Assessing risk factors in an individual's siblings and spouses can inform T2D risk; furthermore, social network based lifestyle interventions involving spouses and siblings might be a novel T2D prevention approach., Competing Interests: Compliance with Ethical Standards Funding SR is supported by NIH National Research Service Award T32HP10251 and the Ryoichi Sasakawa Fellowship Fund. JBM and JD are supported by NIH RO1DK78616, and JBM is also supported by K24 DK080140. We gratefully acknowledge the FHS study participants and staff for their contributions. This work was partially supported by the National Heart, Lung and Blood Institute's Framingham Heart Study (Contract No. N01-HC-25195) and its contract with Affymetrix, Inc for genotyping services (Contract No. N02-HL-6-4278). Prior Presentations Elements of this manuscript were presented (1) orally at the Society of General Internal Medicine Annual Meeting on 24 April 2015 in Toronto, Canada and (2) as a poster at the American Diabetes Association Scientific Sessions on 7 June 2015 in Boston, MA. Conflict of Interest The authors have no conflicts of interest to disclose pertaining to this work.

The genetic architecture of common traits, including the number, frequency, and effect sizes of inherited variants that contribute to individual risk, has been long debated. Genome-wide association studies have identified scores of common variants associated with type 2 diabetes, but in aggregate, these explain only a fraction of the heritability of this disease. Here, to test the hypothesis that lower-frequency variants explain much of the remainder, the GoT2D and T2D-GENES consortia performed whole-genome sequencing in 2,657 European individuals with and without diabetes, and exome sequencing in 12,940 individuals from five ancestry groups. To increase statistical power, we expanded the sample size via genotyping and imputation in a further 111,548 subjects. Variants associated with type 2 diabetes after sequencing were overwhelmingly common and most fell within regions previously identified by genome-wide association studies. Comprehensive enumeration of sequence variation is necessary to identify functional alleles that provide important clues to disease pathophysiology, but large-scale sequencing does not support the idea that lower-frequency variants have a major role in predisposition to type 2 diabetes.

Knowledge of the genetic basis of the type 2 diabetes (T2D)-related quantitative traits fasting glucose (FG) and insulin (FI) in African ancestry (AA) individuals has been limited. In non-diabetic subjects of AA (n = 20,209) and European ancestry (EA; n = 57,292), we performed trans-ethnic (AA+EA) fine-mapping of 54 established EA FG or FI loci with detailed functional annotation, assessed their relevance in AA individuals, and sought previously undescribed loci through trans-ethnic (AA+EA) meta-analysis. We narrowed credible sets of variants driving association signals for 22/54 EA-associated loci; 18/22 credible sets overlapped with active islet-specific enhancers or transcription factor (TF) binding sites, and 21/22 contained at least one TF motif. Of the 54 EA-associated loci, 23 were shared between EA and AA. Replication with an additional 10,096 AA individuals identified two previously undescribed FI loci, chrX FAM133A (rs213676) and chr5 PELO (rs6450057). Trans-ethnic analyses with regulatory annotation illuminate the genetic architecture of glycemic traits and suggest gene regulation as a target to advance precision medicine for T2D. Our approach to utilize state-of-the-art functional annotation and implement trans-ethnic association analysis for discovery and fine-mapping offers a framework for further follow-up and characterization of GWAS signals of complex trait loci., (Copyright © 2016 American Society of Human Genetics. All rights reserved.)

To gain insight into potential regulatory mechanisms through which the effects of variants at four established type 2 diabetes (T2D) susceptibility loci (CDKAL1, CDKN2A-B, IGF2BP2 and KCNQ1) are mediated, we undertook transancestral fine-mapping in 22 086 cases and 42 539 controls of East Asian, European, South Asian, African American and Mexican American descent. Through high-density imputation and conditional analyses, we identified seven distinct association signals at these four loci, each with allelic effects on T2D susceptibility that were homogenous across ancestry groups. By leveraging differences in the structure of linkage disequilibrium between diverse populations, and increased sample size, we localised the variants most likely to drive each distinct association signal. We demonstrated that integration of these genetic fine-mapping data with genomic annotation can highlight potential causal regulatory elements in T2D-relevant tissues. These analyses provide insight into the mechanisms through which T2D association signals are mediated, and suggest future routes to understanding the biology of specific disease susceptibility loci., (© The Author 2016. Published by Oxford University Press.)

Objective: Type 2 diabetes (T2D) is associated with increased mortality in ethnically diverse populations, although the extent to which this association is genetically determined is unknown. We sought to determine whether T2D-related genetic variants predicted all-cause mortality, even after accounting for BMI, in the Third National Health and Nutrition Examination Survey., Research Design and Methods: We modeled mortality risk using a genetic risk score (GRS) from a weighted sum of risk alleles at 38 T2D-related single nucleotide polymorphisms. In age-, sex-, and BMI-adjusted logistic regression models, accounting for the complex survey design, we tested the association with mortality in 6,501 participants. We repeated the analysis within ethnicities (2,528 non-Hispanic white [NHW], 1,979 non-Hispanic black [NHB], and 1,994 Mexican American [MA]) and within BMI categories (<25, 25-30, and ≥30 kg/m(2)). Significance was set at P < 0.05., Results: Over 17 years, 1,556 participants died. GRS was associated with mortality risk (OR 1.04 [95% CI 1.00-1.07] per T2D-associated risk allele, P = 0.05). Within ethnicities, GRS was positively associated with mortality risk in NHW and NHB, but not in MA (0.95 [0.90-1.01], P = 0.07). The negative trend in MA was largely driven by those with BMI <25 kg/m(2) (0.91 [0.82-1.00]). In NHW, the positive association was strongest among those with BMI ≥30 kg/m(2) (1.07 [1.02-1.12])., Conclusions: In the U.S., a higher T2D genetic risk was associated with increased mortality risk, especially among obese NHW. The underlying genetic basis for mortality likely involves complex interactions with factors related to ethnicity, T2D, and body weight., (© 2016 by the American Diabetes Association. Readers may use this article as long as the work is properly cited, the use is educational and not for profit, and the work is not altered.)

We performed fine mapping of 39 established type 2 diabetes (T2D) loci in 27,206 cases and 57,574 controls of European ancestry. We identified 49 distinct association signals at these loci, including five mapping in or near KCNQ1. 'Credible sets' of the variants most likely to drive each distinct signal mapped predominantly to noncoding sequence, implying that association with T2D is mediated through gene regulation. Credible set variants were enriched for overlap with FOXA2 chromatin immunoprecipitation binding sites in human islet and liver cells, including at MTNR1B, where fine mapping implicated rs10830963 as driving T2D association. We confirmed that the T2D risk allele for this SNP increases FOXA2-bound enhancer activity in islet- and liver-derived cells. We observed allele-specific differences in NEUROD1 binding in islet-derived cells, consistent with evidence that the T2D risk allele increases islet MTNR1B expression. Our study demonstrates how integration of genetic and genomic information can define molecular mechanisms through which variants underlying association signals exert their effects on disease.

Background: Interleukin-18 (IL-18) is a pleiotropic cytokine centrally involved in the cytokine cascade with complex immunomodulatory functions in innate and acquired immunity. Circulating IL-18 concentrations are associated with type 2 diabetes mellitus, cardiovascular events, and diverse inflammatory and autoimmune disorders., Methods and Results: To identify causal variants affecting circulating IL-18 concentrations, we applied various omics and molecular biology approaches. By genome-wide association study, we confirmed association of IL-18 levels with a single nucleotide polymorphism in the untranslated exon 2 of the inflammasome component NLRC4 (NLR family, caspase recruitment domain-containing 4) gene on chromosome 2 (rs385076; P=2.4 × 10(-45)). Subsequent molecular analyses by gene expression analysis and reporter gene assays indicated an effect of rs385076 on NLRC4 expression and differential isoform usage by modulating binding of the transcription factor PU.1., Conclusions: Our study provides evidence for the functional causality of single nucleotide polymorphism rs385076 within the NLRC4 gene in relation to IL-18 activation., (© 2015 American Heart Association, Inc.)

Background: Type 2 diabetes mellitus (T2D) and cardiovascular disease share risk factors and subclinical atherosclerosis (SCA) predicts events in those with and without diabetes mellitus. T2D genetic risk may predict both T2D and SCA. We hypothesized that greater T2D genetic risk is associated with higher extent of SCA., Methods and Results: In a cross-sectional analysis, including ≤9210 European Americans, 3773 African Americans, 1446 Hispanic Americans, and 773 Chinese Americans without known cardiovascular disease and enrolled in the Framingham Heart Study, Coronary Artery Risk Development in Young Adults, Multi-Ethnic Study of Atherosclerosis, and Genetic Epidemiology Network of Arteriopathy studies, we tested a 62 T2D-loci genetic risk score for association with measures of SCA, including coronary artery or abdominal aortic calcium score, common and internal carotid artery intima-media thickness, and ankle-brachial index. We used ancestry-stratified linear regression models, with random effects accounting for family relatedness when appropriate, applying a genetic-only (adjusted for sex) and a full SCA risk factors-adjusted model (significance, P<0.01=0.05/5, number of traits analyzed). An inverse association with coronary artery calcium score in Multi-Ethnic Study of Atherosclerosis Europeans (fully-adjusted P=0.004) and with common carotid artery intima-media thickness in the Framingham Heart Study (P=0.009) was not confirmed in other study cohorts, either separately or in meta-analysis. Secondary analyses showed no consistent associations with β-cell and insulin resistance genetic risk sub-scores in the Framingham Heart Study and in the Coronary Artery Risk Development in Young Adults., Conclusions: SCA does not have a major genetic component linked to a burden of 62 T2D loci identified by large genome-wide association studies. A shared T2D-SCA genetic basis, if any, might become apparent from better functional information about both T2D and cardiovascular disease risk loci., (© 2015 American Heart Association, Inc.)

Aims/hypothesis: Type 2 diabetes mellitus in parents is a strong determinant of diabetes risk in their offspring. We hypothesise that offspring diabetes risk associated with parental diabetes is mediated by metabolic risk factors., Methods: We studied initially non-diabetic participants of the Framingham Offspring Study. Metabolic risk was estimated using beta cell corrected insulin response (CIR), HOMA-IR or a count of metabolic syndrome components (metabolic syndrome score [MSS]). Dietary risk and physical activity were estimated using questionnaire responses. Genetic risk score (GRS) was estimated as the count of 62 type 2 diabetes risk alleles. The outcome of incident diabetes in offspring was examined across levels of parental diabetes exposure, accounting for sibling correlation and adjusting for age, sex and putative mediators. The proportion mediated was estimated by comparing regression coefficients for parental diabetes with (β adj) and without (β unadj) adjustments for CIR, HOMA-IR, MSS and GRS (percentage mediated = 1 - β adj / β unadj)., Results: Metabolic factors mediated 11% of offspring diabetes risk associated with parental diabetes, corresponding to a reduction in OR per diabetic parent from 2.13 to 1.96. GRS mediated 9% of risk, corresponding to a reduction in OR per diabetic parent from 2.13 to 1.99., Conclusions/interpretation: Metabolic risk factors partially mediated offspring type 2 diabetes risk conferred by parental diabetes to a similar magnitude as genetic risk. However, a substantial proportion of offspring diabetes risk associated with parental diabetes remains unexplained by metabolic factors, genetic risk, diet and physical activity, suggesting that important familial influences on diabetes risk remain undiscovered.

Fasting glucose and insulin are intermediate traits for type 2 diabetes. Here we explore the role of coding variation on these traits by analysis of variants on the HumanExome BeadChip in 60,564 non-diabetic individuals and in 16,491 T2D cases and 81,877 controls. We identify a novel association of a low-frequency nonsynonymous SNV in GLP1R (A316T; rs10305492; MAF=1.4%) with lower FG (β=-0.09±0.01 mmol l(-1), P=3.4 × 10(-12)), T2D risk (OR[95%CI]=0.86[0.76-0.96], P=0.010), early insulin secretion (β=-0.07±0.035 pmolinsulin mmolglucose(-1), P=0.048), but higher 2-h glucose (β=0.16±0.05 mmol l(-1), P=4.3 × 10(-4)). We identify a gene-based association with FG at G6PC2 (pSKAT=6.8 × 10(-6)) driven by four rare protein-coding SNVs (H177Y, Y207S, R283X and S324P). We identify rs651007 (MAF=20%) in the first intron of ABO at the putative promoter of an antisense lncRNA, associating with higher FG (β=0.02±0.004 mmol l(-1), P=1.3 × 10(-8)). Our approach identifies novel coding variant associations and extends the allelic spectrum of variation underlying diabetes-related quantitative traits and T2D susceptibility.

Genome wide association studies (GWAS) for fasting glucose (FG) and insulin (FI) have identified common variant signals which explain 4.8% and 1.2% of trait variance, respectively. It is hypothesized that low-frequency and rare variants could contribute substantially to unexplained genetic variance. To test this, we analyzed exome-array data from up to 33,231 non-diabetic individuals of European ancestry. We found exome-wide significant (P<5×10-7) evidence for two loci not previously highlighted by common variant GWAS: GLP1R (p.Ala316Thr, minor allele frequency (MAF)=1.5%) influencing FG levels, and URB2 (p.Glu594Val, MAF = 0.1%) influencing FI levels. Coding variant associations can highlight potential effector genes at (non-coding) GWAS signals. At the G6PC2/ABCB11 locus, we identified multiple coding variants in G6PC2 (p.Val219Leu, p.His177Tyr, and p.Tyr207Ser) influencing FG levels, conditionally independent of each other and the non-coding GWAS signal. In vitro assays demonstrate that these associated coding alleles result in reduced protein abundance via proteasomal degradation, establishing G6PC2 as an effector gene at this locus. Reconciliation of single-variant associations and functional effects was only possible when haplotype phase was considered. In contrast to earlier reports suggesting that, paradoxically, glucose-raising alleles at this locus are protective against type 2 diabetes (T2D), the p.Val219Leu G6PC2 variant displayed a modest but directionally consistent association with T2D risk. Coding variant associations for glycemic traits in GWAS signals highlight PCSK1, RREB1, and ZHX3 as likely effector transcripts. These coding variant association signals do not have a major impact on the trait variance explained, but they do provide valuable biological insights.

Genome-wide association studies (GWAS) may have reached their limit of detecting common type 2 diabetes (T2D)-associated genetic variation. We evaluated the performance of current polygenic T2D prediction. Using data from the Framingham Offspring (FOS) and the Coronary Artery Risk Development in Young Adults (CARDIA) studies, we tested three hypotheses: 1) a 62-locus genotype risk score (GRSt) improves T2D prediction compared with previous less inclusive GRSt; 2) separate GRS for β-cell (GRSβ) and insulin resistance (GRSIR) independently predict T2D; and 3) the relationships between T2D and GRSt, GRSβ, or GRSIR do not differ between blacks and whites. Among 1,650 young white adults in CARDIA, 820 young black adults in CARDIA, and 3,471 white middle-aged adults in FOS, cumulative T2D incidence was 5.9%, 14.4%, and 12.9%, respectively, over 25 years. The 62-locus GRSt was significantly associated with incident T2D in all three groups. In FOS but not CARDIA, the 62-locus GRSt improved the model C statistic (0.698 and 0.726 for models without and with GRSt, respectively; P < 0.001) but did not materially improve risk reclassification in either study. Results were similar among blacks compared with whites. The GRSβ but not GRSIR predicted incident T2D among FOS and CARDIA whites. At the end of the era of common variant discovery for T2D, polygenic scores can predict T2D in whites and blacks but do not outperform clinical models. Further optimization of polygenic prediction may require novel analytic methods, including less common as well as functional variants., (© 2014 by the American Diabetes Association.)

Background: Common variation at the 11p11.2 locus, encompassing MADD, ACP2, NR1H3, MYBPC3, and SPI1, has been associated in genome-wide association studies with fasting glucose and insulin (FI). In the Cohorts for Heart and Aging Research in Genomic Epidemiology Targeted Sequencing Study, we sequenced 5 gene regions at 11p11.2 to identify rare, potentially functional variants influencing fasting glucose or FI levels., Methods and Results: Sequencing (mean depth, 38×) across 16.1 kb in 3566 individuals without diabetes mellitus identified 653 variants, 79.9% of which were rare (minor allele frequency <1%) and novel. We analyzed rare variants in 5 gene regions with FI or fasting glucose using the sequence kernel association test. At NR1H3, 53 rare variants were jointly associated with FI (P=2.73×10(-3)); of these, 7 were predicted to have regulatory function and showed association with FI (P=1.28×10(-3)). Conditioning on 2 previously associated variants at MADD (rs7944584, rs10838687) did not attenuate this association, suggesting that there are >2 independent signals at 11p11.2. One predicted regulatory variant, chr11:47227430 (hg18; minor allele frequency=0.00068), contributed 20.6% to the overall sequence kernel association test score at NR1H3, lies in intron 2 of NR1H3, and is a predicted binding site for forkhead box A1 (FOXA1), a transcription factor associated with insulin regulation. In human HepG2 hepatoma cells, the rare chr11:47227430 A allele disrupted FOXA1 binding and reduced FOXA1-dependent transcriptional activity., Conclusions: Sequencing at 11p11.2-NR1H3 identified rare variation associated with FI. One variant, chr11:47227430, seems to be functional, with the rare A allele reducing transcription factor FOXA1 binding and FOXA1-dependent transcriptional activity., (© 2014 American Heart Association, Inc.)

Patients with established type 2 diabetes display both β-cell dysfunction and insulin resistance. To define fundamental processes leading to the diabetic state, we examined the relationship between type 2 diabetes risk variants at 37 established susceptibility loci, and indices of proinsulin processing, insulin secretion, and insulin sensitivity. We included data from up to 58,614 nondiabetic subjects with basal measures and 17,327 with dynamic measures. We used additive genetic models with adjustment for sex, age, and BMI, followed by fixed-effects, inverse-variance meta-analyses. Cluster analyses grouped risk loci into five major categories based on their relationship to these continuous glycemic phenotypes. The first cluster (PPARG, KLF14, IRS1, GCKR) was characterized by primary effects on insulin sensitivity. The second cluster (MTNR1B, GCK) featured risk alleles associated with reduced insulin secretion and fasting hyperglycemia. ARAP1 constituted a third cluster characterized by defects in insulin processing. A fourth cluster (TCF7L2, SLC30A8, HHEX/IDE, CDKAL1, CDKN2A/2B) was defined by loci influencing insulin processing and secretion without a detectable change in fasting glucose levels. The final group contained 20 risk loci with no clear-cut associations to continuous glycemic traits. By assembling extensive data on continuous glycemic traits, we have exposed the diverse mechanisms whereby type 2 diabetes risk variants impact disease predisposition., (© 2014 by the American Diabetes Association.)

To further understanding of the genetic basis of type 2 diabetes (T2D) susceptibility, we aggregated published meta-analyses of genome-wide association studies (GWAS), including 26,488 cases and 83,964 controls of European, east Asian, south Asian and Mexican and Mexican American ancestry. We observed a significant excess in the directional consistency of T2D risk alleles across ancestry groups, even at SNPs demonstrating only weak evidence of association. By following up the strongest signals of association from the trans-ethnic meta-analysis in an additional 21,491 cases and 55,647 controls of European ancestry, we identified seven new T2D susceptibility loci. Furthermore, we observed considerable improvements in the fine-mapping resolution of common variant association signals at several T2D susceptibility loci. These observations highlight the benefits of trans-ethnic GWAS for the discovery and characterization of complex trait loci and emphasize an exciting opportunity to extend insight into the genetic architecture and pathogenesis of human diseases across populations of diverse ancestry.

Adiponectin is strongly inversely associated with insulin resistance and type 2 diabetes, but its causal role remains controversial. We used a Mendelian randomization approach to test the hypothesis that adiponectin causally influences insulin resistance and type 2 diabetes. We used genetic variants at the ADIPOQ gene as instruments to calculate a regression slope between adiponectin levels and metabolic traits (up to 31,000 individuals) and a combination of instrumental variables and summary statistics-based genetic risk scores to test the associations with gold-standard measures of insulin sensitivity (2,969 individuals) and type 2 diabetes (15,960 case subjects and 64,731 control subjects). In conventional regression analyses, a 1-SD decrease in adiponectin levels was correlated with a 0.31-SD (95% CI 0.26-0.35) increase in fasting insulin, a 0.34-SD (0.30-0.38) decrease in insulin sensitivity, and a type 2 diabetes odds ratio (OR) of 1.75 (1.47-2.13). The instrumental variable analysis revealed no evidence of a causal association between genetically lower circulating adiponectin and higher fasting insulin (0.02 SD; 95% CI -0.07 to 0.11; N = 29,771), nominal evidence of a causal relationship with lower insulin sensitivity (-0.20 SD; 95% CI -0.38 to -0.02; N = 1,860), and no evidence of a relationship with type 2 diabetes (OR 0.94; 95% CI 0.75-1.19; N = 2,777 case subjects and 13,011 control subjects). Using the ADIPOQ summary statistics genetic risk scores, we found no evidence of an association between adiponectin-lowering alleles and insulin sensitivity (effect per weighted adiponectin-lowering allele: -0.03 SD; 95% CI -0.07 to 0.01; N = 2,969) or type 2 diabetes (OR per weighted adiponectin-lowering allele: 0.99; 95% CI 0.95-1.04; 15,960 case subjects vs. 64,731 control subjects). These results do not provide any consistent evidence that interventions aimed at increasing adiponectin levels will improve insulin sensitivity or risk of type 2 diabetes.

Objective: To test whether pancreatic beta-cell genetic frailty and hypertension (HTN) interact in their associations with change over time in fasting glucose (ΔFG) or type 2 diabetes mellitus (T2D) risk., Methods and Results: We pooled data from 3471 Framingham Offspring Study participants into six ∼4-year periods (15 852 person-examinations; mean age 52; 54% women). We defined two genetic exposures reflecting beta-cell genetic risk burden: single nucleotide polymorphism (SNP) score counts of fasting glucose-associated and T2D-associated risk alleles at 16 and 33 putative beta-cell loci, respectively; and three HTN exposures: HTN versus no-HTN; treated versus untreated HTN; and five mutually exclusive antihypertensive categories (beta-blockers, thiazides, renin-angiotensin system agents, combinations, others) versus untreated HTN. We tested ∼4-year mean ΔFG or odds of T2D by per-risk allele score change and HTN category, seeking genetic score-by-HTN interaction. Genetic scores increased ∼4-year ΔFG (0.6  mg/dl per-risk allele; P = 8.9 × 10(-16)) and T2D-risk (∼17% per-risk allele; P = 2.1 × 10(-7)). As compared to no-HTN, HTN conferred higher ΔFG (2.6 versus 1.7 mg/dl; P < 0.0001) and T2D-risk [odds ratio (OR) = 2.9, 95% confidence interval (CI) 2.8-3.0; P < 0.0001]. As compared to untreated HTN, treated HTN conferred higher ΔFG (3.4 versus 3.0 mg/dl; P < 0.0001) and T2D-risk (OR = 1.4, 95% CI 1.3-1.5; P = 0.02). Beta-blockers (OR = 1.6, 95% CI 1.1-2.4), combinations (OR = 1.6, 95% CI 1.1-2.5), and others (OR = 2.0, 95% CI 1.4-2.9) increased T2D-risk (all P < 0.02). In joint models including interaction terms, all genetic score-by-HTN interaction terms were P value greater than 0.05. In joint models without interaction, fasting glucose-SNP or T2D-SNP genetic scores (both P < 0.001) and HTN (P < 0.0001) independently increased ΔFG or T2D-risk., Conclusion: HTN, HTN treatment, and common fasting glucose-SNP genetic score/T2D-SNP genetic score independently predicted ΔFG and T2D incidence, but did not modify each other's association with ΔFG or T2D risk.

Type 2 diabetes (T2D) disproportionally affects African Americans (AfA) but, to date, genetic variants identified from genome-wide association studies (GWAS) are primarily from European and Asian populations. We examined the single nucleotide polymorphism (SNP) and locus transferability of 40 reported T2D loci in six AfA GWAS consisting of 2,806 T2D case subjects with or without end-stage renal disease and 4,265 control subjects from the Candidate Gene Association Resource Plus Study. Our results revealed that seven index SNPs at the TCF7L2, KLF14, KCNQ1, ADCY5, CDKAL1, JAZF1, and GCKR loci were significantly associated with T2D (P < 0.05). The strongest association was observed at TCF7L2 rs7903146 (odds ratio [OR] 1.30; P = 6.86 × 10⁻⁸). Locus-wide analysis demonstrated significant associations (P(emp) < 0.05) at regional best SNPs in the TCF7L2, KLF14, and HMGA2 loci as well as suggestive signals in KCNQ1 after correction for the effective number of SNPs at each locus. Of these loci, the regional best SNPs were in differential linkage disequilibrium (LD) with the index and adjacent SNPs. Our findings suggest that some loci discovered in prior reports affect T2D susceptibility in AfA with similar effect sizes. The reduced and differential LD pattern in AfA compared with European and Asian populations may facilitate fine mapping of causal variants at loci shared across populations.

Aims/hypothesis: Hyperglycaemia disproportionately affects African-Americans (AfAs). We tested the transferability of 18 single-nucleotide polymorphisms (SNPs) associated with glycaemic traits identified in European ancestry (EuA) populations in 5,984 non-diabetic AfAs., Methods: We meta-analysed SNP associations with fasting glucose (FG) or insulin (FI) in AfAs from five cohorts in the Candidate Gene Association Resource. We: (1) calculated allele frequency differences, variations in linkage disequilibrium (LD), fixation indices (F(st)s) and integrated haplotype scores (iHSs); (2) tested EuA SNPs in AfAs; and (3) interrogated within ± 250 kb around each EuA SNP in AfAs., Results: Allele frequency differences ranged from 0.6% to 54%. F(st) exceeded 0.15 at 6/16 loci, indicating modest population differentiation. All iHSs were <2, suggesting no recent positive selection. For 18 SNPs, all directions of effect were the same and 95% CIs of association overlapped when comparing EuA with AfA. For 17 of 18 loci, at least one SNP was nominally associated with FG in AfAs. Four loci were significantly associated with FG (GCK, p = 5.8 × 10(-8); MTNR1B, p = 8.5 × 10(-9); and FADS1, p = 2.2 × 10(-4)) or FI (GCKR, p = 5.9 × 10(-4)). At GCK and MTNR1B the EuA and AfA SNPs represented the same signal, while at FADS1, and GCKR, the EuA and best AfA SNPs were weakly correlated (r(2) <0.2), suggesting allelic heterogeneity for association with FG at these loci., Conclusions/interpretation: Few glycaemic SNPs showed strict evidence of transferability from EuA to AfAs. Four loci were significantly associated in both AfAs and those with EuA after accounting for varying LD across ancestral groups, with new signals emerging to aid fine-mapping.

To extend understanding of the genetic architecture and molecular basis of type 2 diabetes (T2D), we conducted a meta-analysis of genetic variants on the Metabochip, including 34,840 cases and 114,981 controls, overwhelmingly of European descent. We identified ten previously unreported T2D susceptibility loci, including two showing sex-differentiated association. Genome-wide analyses of these data are consistent with a long tail of additional common variant loci explaining much of the variation in susceptibility to T2D. Exploration of the enlarged set of susceptibility loci implicates several processes, including CREBBP-related transcription, adipocytokine signaling and cell cycle regulation, in diabetes pathogenesis.