The application of multiple omics technologies in biomedical cohorts has the potential to reveal patient-level disease characteristics and individualized response to treatment. However, the scale... Show moreThe application of multiple omics technologies in biomedical cohorts has the potential to reveal patient-level disease characteristics and individualized response to treatment. However, the scale and heterogeneous nature of multi-modal data makes integration and inference a non-trivial task. We developed a deep-learning-based framework, multi-omics variational autoencoders (MOVE), to integrate such data and applied it to a cohort of 789 people with newly diagnosed type 2 diabetes with deep multi-omics phenotyping from the DIRECT consortium. Using in silico perturbations, we identified drug-omics associations across the multi-modal datasets for the 20 most prevalent drugs given to people with type 2 diabetes with substantially higher sensitivity than univariate statistical tests. From these, we among others, identified novel associations between metformin and the gut microbiota as well as opposite molecular responses for the two statins, simvastatin and atorvastatin. We used the associations to quantify drug-drug similarities, assess the degree of polypharmacy and conclude that drug effects are distributed across the multi-omics modalities. 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
The presentation and underlying pathophysiology of type 2 diabetes (T2D) is complex and heterogeneous. Recent studies attempted to stratify T2D into distinct subgroups using data-driven approaches,... Show moreThe presentation and underlying pathophysiology of type 2 diabetes (T2D) is complex and heterogeneous. Recent studies attempted to stratify T2D into distinct subgroups using data-driven approaches, but their clinical utility may be limited if categorical representations of complex phenotypes are suboptimal. We apply a soft-clustering (archetype) method to characterize newly diagnosed T2D based on 32 clinical variables. We assign quantitative clustering scores for individuals and investigate the associations with glycemic deterioration, genetic risk scores, circulating omics biomarkers, and phenotypic stability over 36 months. Four archetype profiles represent dysfunction patterns across combinations of T2D etiological processes and correlate with multiple circulating biomarkers. One archetype associated with obesity, insulin resistance, dyslipidemia, and impaired 1 beta cell glucose sensitivity corresponds with the fastest disease progression and highest demand for anti-diabetic treatment. We demonstrate that clinical heterogeneity in T2D can be mapped to heterogeneity in individual etiological processes, providing a potential route to personalized treatments. 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
AimSubclasses of different glycaemic disturbances could explain the variation in characteristics of individuals with type 2 diabetes (T2D). We aimed to examine the association between subgroups... Show moreAimSubclasses of different glycaemic disturbances could explain the variation in characteristics of individuals with type 2 diabetes (T2D). We aimed to examine the association between subgroups based on their glucose curves during a five-point mixed-meal tolerance test (MMT) and metabolic traits at baseline and glycaemic deterioration in individuals with T2D.MethodsThe study included 787 individuals with newly diagnosed T2D from the Diabetes Research on Patient Stratification (IMI-DIRECT) Study. Latent class trajectory analysis (LCTA) was used to identify distinct glucose curve subgroups during a five-point MMT. Using general linear models, these subgroups were associated with metabolic traits at baseline and after 18 months of follow up, adjusted for potential confounders.ResultsAt baseline, we identified three glucose curve subgroups, labelled in order of increasing glucose peak levels as subgroup 1-3. Individuals in subgroup 2 and 3 were more likely to have higher levels of HbA1c, triglycerides and BMI at baseline, compared to those in subgroup 1. At 18 months (n = 651), the beta coefficients (95% CI) for change in HbA1c (mmol/mol) increased across subgroups with 0.37 (-0.18-1.92) for subgroup 2 and 1.88 (-0.08-3.85) for subgroup 3, relative to subgroup 1. The same trend was observed for change in levels of triglycerides and fasting glucose.ConclusionsDifferent glycaemic profiles with different metabolic traits and different degrees of subsequent glycaemic deterioration can be identified using data from a frequently sampled mixed-meal tolerance test in individuals with T2D. Subgroups with the highest peaks had greater metabolic risk. Show less
Obura, M.; Beulens, J.W.J.; Slieker, R.; Koopman, A.D.M.; Hoekstra, T.; Nijpels, G.; ... ; Rutters, F. 2020
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... Show moreAim 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. Show less
