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
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
Winters-Van Eekelen, E.; Velde, J.H.P.M. van der; Boone, S.C.; Westgate, K.; Brage, S.; Lamb, H.J.; ... ; Mutsert, R. de 2021
Purpose: It remains unclear to what extent habitual physical activity and sedentary time (ST) are associated with visceral fat and liver fat. We studied the substitution of ST with time spent... Show morePurpose: It remains unclear to what extent habitual physical activity and sedentary time (ST) are associated with visceral fat and liver fat. We studied the substitution of ST with time spent physically active and total body fat (TBF), visceral adipose tissue (VAT), and hepatic triglyceride content (HTGC) in middle-age men and women. Design: In this cross-sectional analysis of the Netherlands Epidemiology of Obesity study, physical activity was assessed in 228 participants using a combined accelerometer and heart rate monitor. TBF was assessed by the Tanita bioelectrical impedance, VAT by magnetic resonance imaging, and HTGC by proton-MR spectroscopy. Behavioral intensity distribution was categorized as ST, time spent in light physical activity (LPA), and moderate to vigorous physical activity (MVPA). To estimate the effect of replacing 30 min.d(-1) of ST with 30 min.d(-1) LPA or MVPA, we performed isotemporal substitution analyses, adjusted for sex, age, ethnicity, education, the Dutch Healthy Diet index, and smoking. Results: Included participants (41% men) had a mean +/- SD age of 56 +/- 6 yr and spent 88 +/- 56 min in MVPA and 9.0 +/- 2.1 h of ST. Replacing 30 min.d(-1) of ST with 30 min of MVPA was associated with 1.3% less TBF (95% confidence interval = -2.0 to -0.7), 7.8 cm(2) less VAT (-11.6 to -4.0), and 0.89 times HTGC (0.82-0.97). Replacement with LPA was not associated with TBF (-0.03%; -0.5 to 0.4), VAT (-1.7 cm(2); -4.4 to 0.9), or HTGC (0.98 times; 0.92-1.04). Conclusions: Reallocation of time spent sedentarywith time spent inMVPA, but not LPA, was associated with less TBF, visceral fat, and liver fat. These findings contribute to the development of more specified guidelines on ST and physical activity. 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
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
