We evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple... Show moreWe evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant associates with multiple molecular phenotypes over multiple genomic regions. The highest proportion of share genetic regulation is detected between gene expression and proteins (66.6%), with a further median shared genetic associations across 49 different tissues of 78.3% and 62.4% between plasma proteins and gene expression. We represent the genetic and molecular associations in networks including 2828 known GWAS variants, showing that GWAS variants are more often connected to gene expression in trans than other molecular phenotypes in the network. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants using different molecular phenotypes in an accessible tissue. Show less
Background In the treatment of type 2 diabetes, GLP-1 receptor agonists lower blood glucose concentrations, body weight, and have cardiovascular benefits. The efficacy and side effects of GLP-1... Show moreBackground In the treatment of type 2 diabetes, GLP-1 receptor agonists lower blood glucose concentrations, body weight, and have cardiovascular benefits. The efficacy and side effects of GLP-1 receptor agonists vary between people. Human pharmacogenomic studies of this inter-individual variation can provide both biological insight into drug action and provide biomarkers to inform clinical decision making. We therefore aimed to identify genetic variants associated with glycaemic response to GLP-1 receptor agonist treatment. Methods In this genome-wide analysis we included adults (aged >= 18 years) with type 2 diabetes treated with GLP-1 receptor agonists with baseline HbA1c of 7% or more (53 mmol/mol) from four prospective observational cohorts (DIRECT, PRIBA, PROMASTER, and GoDARTS) and two randomised clinical trials (HARMONY phase 3 and AWARD). The primary endpoint was HbA1c reduction at 6 months after starting GLP-1 receptor agonists. We evaluated variants in GLP1R, then did a genome-wide association study and gene-based burden tests. Findings 4571 adults were included in our analysis, of these, 3339 (73%) were White European, 449 (10%) Hispanic, 312 (7%) American Indian or Alaskan Native, and 471 (10%) were other, and around 2140 (47%) of the participants were women. Variation in HbA1c reduction with GLP-1 receptor agonists treatment was associated with rs6923761G -> A (Gly168Ser) in the GLP1R (0.08% [95% CI 0.04-0.12] or 0.9 mmol/mol lower reduction in HbA1c per serine, p=6.0 x 10(-5)) and low frequency variants in ARRB1 (optimal sequence kernel association test p=6.7 x 10(-8)), largely driven by rs140226575G -> A (Thr370Met; 0.25% [SE 0.06] or 2.7 mmol/mol [SE 0.7] greater HbA1c reduction per methionine, p=5.2 x 10(-6)). A similar effect size for the ARRB1 Thr370Met was seen in Hispanic and American Indian or Alaska Native populations who have a higher frequency of this variant (6-11%) than in White European populations. Combining these two genes identified 4% of the population who had a 30% greater reduction in HbA1c than the 9% of the population with the worse response. Interpretation This genome-wide pharmacogenomic study of GLP-1 receptor agonists provides novel biological and clinical insights. Clinically, when genotype is routinely available at the point of prescribing, individuals with ARRB1 variants might benefit from earlier initiation of GLP-1 receptor agonists. Copyright (C) 2022 The Author(s). Published by Elsevier Ltd. Show less
OBJECTIVEWe investigated the processes underlying glycemic deterioration in type 2 diabetes (T2D).RESEARCH DESIGN AND METHODSA total of 732 recently diagnosed patients with T2D from the Innovative... Show moreOBJECTIVEWe investigated the processes underlying glycemic deterioration in type 2 diabetes (T2D).RESEARCH DESIGN AND METHODSA total of 732 recently diagnosed patients with T2D from the Innovative Medicines Initiative Diabetes Research on Patient Stratification (IMI DIRECT) study were extensively phenotyped over 3 years, including measures of insulin sensitivity (OGIS), beta-cell glucose sensitivity (GS), and insulin clearance (CLIm) from mixed meal tests, liver enzymes, lipid profiles, and baseline regional fat from MRI. The associations between the longitudinal metabolic patterns and HbA(1c) deterioration, adjusted for changes in BMI and in diabetes medications, were assessed via stepwise multivariable linear and logistic regression.RESULTSFaster HbA(1c) progression was independently associated with faster deterioration of OGIS and GS and increasing CLIm; visceral or liver fat, HDL-cholesterol, and triglycerides had further independent, though weaker, roles (R-2 = 0.38). A subgroup of patients with a markedly higher progression rate (fast progressors) was clearly distinguishable considering these variables only (discrimination capacity from area under the receiver operating characteristic = 0.94). The proportion of fast progressors was reduced from 56% to 8-10% in subgroups in which only one trait among OGIS, GS, and CLIm was relatively stable (odds ratios 0.07-0.09). T2D polygenic risk score and baseline pancreatic fat, glucagon-like peptide 1, glucagon, diet, and physical activity did not show an independent role.CONCLUSIONSDeteriorating insulin sensitivity and beta-cell function, increasing insulin clearance, high visceral or liver fat, and worsening of the lipid profile are the crucial factors mediating glycemic deterioration of patients with T2D in the initial phase of the disease. Stabilization of a single trait among insulin sensitivity, beta-cell function, and insulin clearance may be relevant to prevent progression. Show less
Background The rising prevalence of type 2 diabetes (T2D) poses a major global challenge. It remains unresolved to what extent transcriptomic signatures of metabolic dysregulation and T2D can be... Show moreBackground The rising prevalence of type 2 diabetes (T2D) poses a major global challenge. It remains unresolved to what extent transcriptomic signatures of metabolic dysregulation and T2D can be observed in easily accessible tissues such as blood. Additionally, large-scale human studies are required to further our understanding of the putative inflammatory component of insulin resistance and T2D. Here we used transcriptomics data from individuals with (n = 789) and without (n = 2127) T2D from the IMI-DIRECT cohorts to describe the co-expression structure of whole blood that mainly reflects processes and cell types of the immune system, and how it relates to metabolically relevant clinical traits and T2D. Methods Clusters of co-expressed genes were identified in the non-diabetic IMI-DIRECT cohort and evaluated with regard to stability, as well as preservation and rewiring in the cohort of individuals with T2D. We performed functional and immune cell signature enrichment analyses, and a genome-wide association study to describe the genetic regulation of the modules. Phenotypic and trans-omics associations of the transcriptomic modules were investigated across both IMI-DIRECT cohorts. Results We identified 55 whole blood co-expression modules, some of which clustered in larger super-modules. We identified a large number of associations between these transcriptomic modules and measures of insulin action and glucose tolerance. Some of the metabolically linked modules reflect neutrophil-lymphocyte ratio in blood while others are independent of white blood cell estimates, including a module of genes encoding neutrophil granule proteins with antibacterial properties for which the strongest associations with clinical traits and T2D status were observed. Through the integration of genetic and multi-omics data, we provide a holistic view of the regulation and molecular context of whole blood transcriptomic modules. We furthermore identified an overlap between genetic signals for T2D and co-expression modules involved in type II interferon signaling. Conclusions Our results offer a large-scale map of whole blood transcriptomic modules in the context of metabolic disease and point to novel biological candidates for future studies related to T2D. Show less
Atabaki-Pasdar, N.; Ohlsson, M.; Vinuela, A.; Frau, F.; Pomares-Millan, H.; Haid, M.; ... ; Franks, P.W. 2020
BackgroundNon-alcoholic fatty liver disease (NAFLD) is highly prevalent and causes serious health complications in individuals with and without type 2 diabetes (T2D). Early diagnosis of NAFLD is... Show moreBackgroundNon-alcoholic fatty liver disease (NAFLD) is highly prevalent and causes serious health complications in individuals with and without type 2 diabetes (T2D). Early diagnosis of NAFLD is important, as this can help prevent irreversible damage to the liver and, ultimately, hepatocellular carcinomas. We sought to expand etiological understanding and develop a diagnostic tool for NAFLD using machine learning.Methods and findingsWe utilized the baseline data from IMI DIRECT, a multicenter prospective cohort study of 3,029 European-ancestry adults recently diagnosed with T2D (n= 795) or at high risk of developing the disease (n= 2,234). Multi-omics (genetic, transcriptomic, proteomic, and metabolomic) and clinical (liver enzymes and other serological biomarkers, anthropometry, measures of beta-cell function, insulin sensitivity, and lifestyle) data comprised the key input variables. The models were trained on MRI-image-derived liver fat content (<5% or >= 5%) available for 1,514 participants. We applied LASSO (least absolute shrinkage and selection operator) to select features from the different layers of omics data and random forest analysis to develop the models. The prediction models included clinical and omics variables separately or in combination. A model including all omics and clinical variables yielded a cross-validated receiver operating characteristic area under the curve (ROCAUC) of 0.84 (95% CI 0.82, 0.86;p <0.001), which compared with a ROCAUC of 0.82 (95% CI 0.81, 0.83;p <0.001) for a model including 9 clinically accessible variables. The IMI DIRECT prediction models outperformed existing noninvasive NAFLD prediction tools. One limitation is that these analyses were performed in adults of European ancestry residing in northern Europe, and it is unknown how well these findings will translate to people of other ancestries and exposed to environmental risk factors that differ from those of the present cohort. Another key limitation of this study is that the prediction was done on a binary outcome of liver fat quantity (<5% or >= 5%) rather than a continuous one.ConclusionsIn this study, we developed several models with different combinations of clinical and omics data and identified biological features that appear to be associated with liver fat accumulation. In general, the clinical variables showed better prediction ability than the complex omics variables. However, the combination of omics and clinical variables yielded the highest accuracy. We have incorporated the developed clinical models into a web interface (see:) and made it available to the community. Show less
Dawed, A.Y.; Zhou, K.X.; Leeuwen, N. van; Mahajan, A.; Robertson, N.; Koivula, R.; ... ; IMI DIRECT Consortium 2019
OBJECTIVEGastrointestinal adverse effects occur in 20-30% of patients with metformin-treated type 2 diabetes, leading to premature discontinuation in 5-10% of the cases. Gastrointestinal... Show moreOBJECTIVEGastrointestinal adverse effects occur in 20-30% of patients with metformin-treated type 2 diabetes, leading to premature discontinuation in 5-10% of the cases. Gastrointestinal intolerance may reflect localized high concentrations of metformin in the gut. We hypothesized that reduced transport of metformin via the plasma membrane monoamine transporter (PMAT) and organic cation transporter 1 (OCT1) could increase the risk of severe gastrointestinal adverse effects.RESEARCH DESIGN AND METHODSThe study included 286 severe metformin-intolerant and 1,128 metformin-tolerant individuals from the IMI DIRECT (Innovative Medicines Initiative: DIabetes REsearCh on patient straTification) consortium. We assessed the association of patient characteristics, concomitant medication, and the burden of mutations in the SLC29A4 and SLC22A1 genes on odds of intolerance.RESULTSWomen (P < 0.001) and older people (P < 0.001) were more likely to develop metformin intolerance. Concomitant use of transporter-inhibiting drugs increased the odds of intolerance (odds ratio [OR] 1.72, P < 0.001). In an adjusted logistic regression model, the G allele at rs3889348 (SLC29A4) was associated with gastrointestinal intolerance (OR 1.34, P = 0.005). rs3889348 is the top cis-expression quantitative trait locus for SLC29A4 in gut tissue where carriers of the G allele had reduced expression. Homozygous carriers of the G allele treated with transporter-inhibiting drugs had more than three times higher odds of intolerance compared with carriers of no G allele and not treated with inhibiting drugs (OR 3.23, P < 0.001). Use of a genetic risk score derived from rs3889348 and SLC22A1 variants found that the odds of intolerance were more than twice as high in individuals who carry three or more risk alleles compared with those carrying none (OR 2.15, P = 0.01).CONCLUSIONSThese results suggest that intestinal metformin transporters and concomitant medications play an important role in the gastrointestinal adverse effects of metformin. Show less
OBJECTIVEProgression to insulin therapy in clinically diagnosed type 2 diabetes is highly variable. GAD65 autoantibodies (GADA) are associated with faster progression, but their predictive value is... Show moreOBJECTIVEProgression to insulin therapy in clinically diagnosed type 2 diabetes is highly variable. GAD65 autoantibodies (GADA) are associated with faster progression, but their predictive value is limited. We aimed to determine if a type 1 diabetes genetic risk score (T1D GRS) could predict rapid progression to insulin treatment over and above GADA testing.RESEARCH DESIGN AND METHODSWe examined the relationship between T1D GRS, GADA (negative or positive), and rapid insulin requirement (within 5 years) using Kaplan-Meier survival analysis and Cox regression in 8,608 participants with clinical type 2 diabetes (onset >35 years and treated without insulin for 6 months). T1D GRS was both analyzed continuously (as standardized scores) and categorized based on previously reported centiles of a population with type 1 diabetes (<5th [low], 5th-50th [medium], and >50th [high]).RESULTSIn GADA-positive participants (3.3%), those with higher T1D GRS progressed to insulin more quickly: probability of insulin requirement at 5 years (95% CI): 47.9% (35.0%, 62.78%) (high T1D GRS) vs. 27.6% (20.5%, 36.5%) (medium T1D GRS) vs. 17.6% (11.2%, 27.2%) (low T1D GRS); P = 0.001. In contrast, T1D GRS did not predict rapid insulin requirement in GADA-negative participants (P = 0.4). In Cox regression analysis with adjustment for age of diagnosis, BMI, and cohort, T1D GRS was independently associated with time to insulin only in the presence of GADA: hazard ratio per SD increase was 1.48 (1.15, 1.90); P = 0.002.CONCLUSIONSA T1D GRS alters the clinical implications of a positive GADA test in patients with clinical type 2 diabetes and is independent of and additive to clinical features. Show less