Adipose Tissue Gene Expression Associations Reveal Hundreds of Candidate Genes for Cardiometabolic Traits

Genome-wide association studies (GWASs) have identified thousands of genetic loci associated with cardiometabolic traits including type 2 diabetes (T2D), lipid levels, body fat distribution, and adiposity, although most causal genes remain unknown. We used subcutaneous adipose tissue RNA-seq data from 434 Finnish men from the METSIM study to identify 9,687 primary and 2,785 secondary cis-expression quantitative trait loci (eQTL; <1 Mb from TSS, FDR < 1%). Compared to primary eQTL signals, secondary eQTL signals were located further from transcription start sites, had smaller effect sizes, and were less enriched in adipose tissue regulatory elements compared to primary signals. Among 2,843 cardiometabolic GWAS signals, 262 colocalized by LD and conditional analysis with 318 transcripts as primary and conditionally distinct secondary cis-eQTLs, including some across ancestries. Of cardiometabolic traits examined for adipose tissue eQTL colocalizations, waist-hip ratio (WHR) and circulating lipid traits had the highest percentage of colocalized eQTLs (15% and 14%, respectively). Among alleles associated with increased cardiometabolic GWAS risk, approximately half (53%) were associated with decreased gene expression level. Mediation analyses of colocalized genes and cardiometabolic traits within the 434 individuals provided further evidence that gene expression influences variant-trait associations. These results identify hundreds of candidate genes that may act in adipose tissue to influence cardiometabolic traits. © 2019 American Society of Human Genetic

Clinical profiles of post-load glucose subgroups and their association with glycaemic traits over time: an IMI-DIRECT study

Aim To examine the hypothesis that, based on their glucose curves during a seven-point oral glucose tolerance test, people at elevated type 2 diabetes risk can be divided into subgroups with different clinical profiles at baseline and different degrees of subsequent glycaemic deterioration.Methods We included 2126 participants at elevated type 2 diabetes risk from the Diabetes Research on Patient Stratification (IMI-DIRECT) study. Latent class trajectory analysis was used to identify subgroups from a seven-point oral glucose tolerance test at baseline and follow-up. Linear models quantified the associations between the subgroups with glycaemic traits at baseline and 18 months.Results At baseline, we identified four glucose curve subgroups, labelled in order of increasing peak levels as 1-4. Participants in Subgroups 2-4, were more likely to have higher insulin resistance (homeostatic model assessment) and a lower Matsuda index, than those in Subgroup 1. Overall, participants in Subgroups 3 and 4, had higher glycaemic trait values, with the exception of the Matsuda and insulinogenic indices. At 18 months, change in homeostatic model assessment of insulin resistance was higher in Subgroup 4 (beta = 0.36, 95% CI 0.13-0.58), Subgroup 3 (beta = 0.30; 95% CI 0.10-0.50) and Subgroup 2 (beta = 0.18; 95% CI 0.04-0.32), compared to Subgroup 1. The same was observed for C-peptide and insulin. Five subgroups were identified at follow-up, and the majority of participants remained in the same subgroup or progressed to higher peak subgroups after 18 months.Conclusions Using data from a frequently sampled oral glucose tolerance test, glucose curve patterns associated with different clinical characteristics and different rates of subsequent glycaemic deterioration can be identified.Molecular Epidemiolog

Identification and Functional Characterization of G6PC2 Coding Variants Influencing Glycemic Traits Define an Effector Transcript at the G6PC2-ABCB11 Locus

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&lt;5&times;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

The role of physical activity in metabolic homeostasis before and after the onset of type 2 diabetes: an IMI DIRECT study.

AIMS/HYPOTHESIS: It is well established that physical activity, abdominal ectopic fat and glycaemic regulation are related but the underlying structure of these relationships is unclear. The previously proposed twin-cycle hypothesis (TC) provides a mechanistic basis for impairment in glycaemic control through the interactions of substrate availability, substrate metabolism and abdominal ectopic fat accumulation. Here, we hypothesise that the effect of physical activity in glucose regulation is mediated by the twin-cycle. We aimed to examine this notion in the Innovative Medicines Initiative Diabetes Research on Patient Stratification (IMI DIRECT) Consortium cohorts comprised of participants with normal or impaired glucose regulation (cohort 1: N ≤ 920) or with recently diagnosed type 2 diabetes (cohort 2: N ≤ 435). METHODS: We defined a structural equation model that describes the TC and fitted this within the IMI DIRECT dataset. A second model, twin-cycle plus physical activity (TC-PA), to assess the extent to which the effects of physical activity in glycaemic regulation are mediated by components in the twin-cycle, was also fitted. Beta cell function, insulin sensitivity and glycaemic control were modelled from frequently sampled 75 g OGTTs (fsOGTTs) and mixed-meal tolerance tests (MMTTs) in participants without and with diabetes, respectively. Abdominal fat distribution was assessed using MRI, and physical activity through wrist-worn triaxial accelerometry. Results are presented as standardised beta coefficients, SE and p values, respectively. RESULTS: The TC and TC-PA models showed better fit than null models (TC: χ2 = 242, p = 0.004 and χ2 = 63, p = 0.001 in cohort 1 and 2, respectively; TC-PA: χ2 = 180, p = 0.041 and χ2 = 60, p = 0.008 in cohort 1 and 2, respectively). The association of physical activity with glycaemic control was primarily mediated by variables in the liver fat cycle. CONCLUSIONS/INTERPRETATION: These analyses partially support the mechanisms proposed in the twin-cycle model and highlight mechanistic pathways through which insulin sensitivity and liver fat mediate the association between physical activity and glycaemic control.S.Bra. was funded by the UK Medical Research Council [MC_UU_12015/3]

Parental origin of sequence variants associated with complex diseases

To access publisher full text version of this article. Please click on the hyperlink in Additional Links fieldEffects of susceptibility variants may depend on from which parent they are inherited. Although many associations between sequence variants and human traits have been discovered through genome-wide associations, the impact of parental origin has largely been ignored. Here we show that for 38,167 Icelanders genotyped using single nucleotide polymorphism (SNP) chips, the parental origin of most alleles can be determined. For this we used a combination of genealogy and long-range phasing. We then focused on SNPs that associate with diseases and are within 500 kilobases of known imprinted genes. Seven independent SNP associations were examined. Five-one with breast cancer, one with basal-cell carcinoma and three with type 2 diabetes-have parental-origin-specific associations. These variants are located in two genomic regions, 11p15 and 7q32, each harbouring a cluster of imprinted genes. Furthermore, we observed a novel association between the SNP rs2334499 at 11p15 and type 2 diabetes. Here the allele that confers risk when paternally inherited is protective when maternally transmitted. We identified a differentially methylated CTCF-binding site at 11p15 and demonstrated correlation of rs2334499 with decreased methylation of that site.info:eu-repo/grantAgreement/EC/FP7/21807

A new strategy for enhancing imputation quality of rare variants from next-generation sequencing data via combining SNP and exome chip data

Background: Rare variants have gathered increasing attention as a possible alternative source of missing heritability. Since next generation sequencing technology is not yet cost-effective for large-scale genomic studies, a widely used alternative approach is imputation. However, the imputation approach may be limited by the low accuracy of the imputed rare variants. To improve imputation accuracy of rare variants, various approaches have been suggested, including increasing the sample size of the reference panel, using sequencing data from study-specific samples (i.e., specific populations), and using local reference panels by genotyping or sequencing a subset of study samples. While these approaches mainly utilize reference panels, imputation accuracy of rare variants can also be increased by using exome chips containing rare variants. The exome chip contains 250 K rare variants selected from the discovered variants of about 12,000 sequenced samples. If exome chip data are available for previously genotyped samples, the combined approach using a genotype panel of merged data, including exome chips and SNP chips, should increase the imputation accuracy of rare variants. Results: In this study, we describe a combined imputation which uses both exome chip and SNP chip data simultaneously as a genotype panel. The effectiveness and performance of the combined approach was demonstrated using a reference panel of 848 samples constructed using exome sequencing data from the T2D-GENES consortium and 5,349 sample genotype panels consisting of an exome chip and SNP chip. As a result, the combined approach increased imputation quality up to 11 %, and genomic coverage for rare variants up to 117.7 % (MAF < 1 %), compared to imputation using the SNP chip alone. Also, we investigated the systematic effect of reference panels on imputation quality using five reference panels and three genotype panels. The best performing approach was the combination of the study specific reference panel and the genotype panel of combined data. Conclusions: Our study demonstrates that combined datasets, including SNP chips and exome chips, enhances both the imputation quality and genomic coverage of rare variants

Meta-analysis of type 2 Diabetes in African Americans Consortium

Type 2 diabetes (T2D) is more prevalent in African Americans than in Europeans. However, little is known about the genetic risk in African Americans despite the recent identification of more than 70 T2D loci primarily by genome-wide association studies (GWAS) in individuals of European ancestry. In order to investigate the genetic architecture of T2D in African Americans, the MEta-analysis of type 2 DIabetes in African Americans (MEDIA) Consortium examined 17 GWAS on T2D comprising 8,284 cases and 15,543 controls in African Americans in stage 1 analysis. Single nucleotide polymorphisms (SNPs) association analysis was conducted in each study under the additive model after adjustment for age, sex, study site, and principal components. Meta-analysis of approximately 2.6 million genotyped and imputed SNPs in all studies was conducted using an inverse variance-weighted fixed effect model. Replications were performed to follow up 21 loci in up to 6,061 cases and 5,483 controls in African Americans, and 8,130 cases and 38,987 controls of European ancestry. We identified three known loci (TCF7L2, HMGA2 and KCNQ1) and two novel loci (HLA-B and INS-IGF2) at genome-wide significance (4.15 × 10(-94)<P<5 × 10(-8), odds ratio (OR)  = 1.09 to 1.36). Fine-mapping revealed that 88 of 158 previously identified T2D or glucose homeostasis loci demonstrated nominal to highly significant association (2.2 × 10(-23) < locus-wide P<0.05). These novel and previously identified loci yielded a sibling relative risk of 1.19, explaining 17.5% of the phenotypic variance of T2D on the liability scale in African Americans. Overall, this study identified two novel susceptibility loci for T2D in African Americans. A substantial number of previously reported loci are transferable to African Americans after accounting for linkage disequilibrium, enabling fine mapping of causal variants in trans-ethnic meta-analysis studies.Peer reviewe

New insights into the genetic etiology of Alzheimer's disease and related dementias

Characterization of the genetic landscape of Alzheimer's disease (AD) and related dementias (ADD) provides a unique opportunity for a better understanding of the associated pathophysiological processes. We performed a two-stage genome-wide association study totaling 111,326 clinically diagnosed/'proxy' AD cases and 677,663 controls. We found 75 risk loci, of which 42 were new at the time of analysis. Pathway enrichment analyses confirmed the involvement of amyloid/tau pathways and highlighted microglia implication. Gene prioritization in the new loci identified 31 genes that were suggestive of new genetically associated processes, including the tumor necrosis factor alpha pathway through the linear ubiquitin chain assembly complex. We also built a new genetic risk score associated with the risk of future AD/dementia or progression from mild cognitive impairment to AD/dementia. The improvement in prediction led to a 1.6- to 1.9-fold increase in AD risk from the lowest to the highest decile, in addition to effects of age and the APOE ε4 allele

Determinants of penetrance and variable expressivity in monogenic metabolic conditions across 77,184 exomes

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