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
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
Ghorasaini, M.; Mohammed, Y.; Adamski, J.; Bettcher, L.; Bowden, J.A.; Cabruja, M.; ... ; Giera, M. 2021
Modern biomarker and translational research as well as personalized health care studies rely heavily on powerful omics' technologies, including metabolomics and lipidomics. However, to translate... Show moreModern biomarker and translational research as well as personalized health care studies rely heavily on powerful omics' technologies, including metabolomics and lipidomics. However, to translate metabolomics and lipidomics discoveries into a high-throughput clinical setting, standardization is of utmost importance. Here, we compared and benchmarked a quantitative lipidomics platform. The employed Lipidyzer platform is based on lipid class separation by means of differential mobility spectrometry with subsequent multiple reaction monitoring. Quantitation is achieved by the use of 54 deuterated internal standards and an automated informatics approach. We investigated the platform performance across nine laboratories using NIST SRM 1950-Metabolites in Frozen Human Plasma, and three NIST Candidate Reference Materials 8231-Frozen Human Plasma Suite for Metabolomics (high triglyceride, diabetic, and African-American plasma). In addition, we comparatively analyzed 59 plasma samples from individuals with familial hypercholesterolemia from a clinical cohort study. We provide evidence that the more practical methyl-tert-butyl ether extraction outperforms the classic Bligh and Dyer approach and compare our results with two previously published ring trials. In summary, we present standardized lipidomics protocols, allowing for the highly reproducible analysis of several hundred human plasma lipids, and present detailed molecular information for potentially disease relevant and ethnicity-related materials. 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
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
Chouvarine, P.; Giera, M.; Kastenmuller, G.; Artati, A.; Adamski, J.; Bertram, H.; Hansmann, G. 2020
Objective While metabolic dysfunction occurs in several pulmonary arterial hypertension (PAH) animal models, its role in the human hypertensive right ventricle (RV) and lung is not well... Show moreObjective While metabolic dysfunction occurs in several pulmonary arterial hypertension (PAH) animal models, its role in the human hypertensive right ventricle (RV) and lung is not well characterised. We investigated whether circulating metabolite concentrations differ across the hypertensive RV and/or the pulmonary circulation, and correlate with invasive haemodynamic/echocardiographic variables in patients with PAH.Methods Prospective EDTA blood collection during cardiac catheterisation from the superior vena cava (SVC), pulmonary artery (PA) and ascending aorta (AAO) in children with PAH (no shunt) and non-PAH controls (Con), followed by unbiased screens of 427 metabolites and 836 lipid species and fatty acids (FAs) in blood plasma (Metabolon and Lipidyzer platforms). Metabolite concentrations were correlated with echocardiographic and invasive haemodynamic variables.Results Metabolomics/lipidomics analysis of differential concentrations (false discovery rate<0.15) revealed several metabolite gradients in the trans-RV (PA vs SVC) setting. Notably, dicarboxylic acids (eg, octadecanedioate: fold change (FC)_Control=0.77, FC_PAH=1.09, p value=0.044) and acylcarnitines (eg, stearoylcarnitine: FC_Control=0.74, FC_PAH=1.21, p value=0.058). Differentially regulated metabolites were also found in the transpulmonary (AAO vs PA) setting and between-group comparisons, that is, in the SVC (PAH-S VC vs Con-S VC), PA and AAO. Importantly, the differential PAH-metabolite concentrations correlated with numerous outcome-relevant variables (e.g., tricuspid annular plane systolic excursion, pulmonary vascular resistance).Conclusions In PAH, trans-RV and transpulmonary metabolite gradients exist and correlate with haemodynamic determinants of clinical outcome. The most pronounced differential trans-R V gradients are known to be involved in lipid metabolism/lipotoxicity, that is, accumulation of long chain FAs. The identified accumulation of dicarboxylic acids and acylcarnitines likely indicates impaired beta-oxidation in the hypertensive RV and represents emerging biomarkers and therapeutic targets in PAH. 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
Quell, J.D.; Romisch-Margl, W.; Haid, M.; Krumsiek, J.; Skurk, T.; Halama, A.; ... ; Kastenmuller, G. 2019
Kit-based assays, such as AbsoluteIDQ(TM) p150, are widely used in large cohort studies and provide a standardized method to quantify blood concentrations of phosphatidylcholines (PCs). Many... Show moreKit-based assays, such as AbsoluteIDQ(TM) p150, are widely used in large cohort studies and provide a standardized method to quantify blood concentrations of phosphatidylcholines (PCs). Many disease-relevant associations of PCs were reported using this method. However, their interpretation is hampered by lack of functionally-relevant information on the detailed fatty acid side-chain compositions as only the total number of carbon atoms and double bonds is identified by the kit. To enable more substantiated interpretations, we characterized these PC sums using the side-chain resolving Lipidyzer(TM) platform, analyzing 223 samples in parallel to the AbsoluteIDQ(TM). Combining these datasets, we estimated the quantitative composition of PC sums and subsequently tested their replication in an independent cohort. We identified major constituents of 28 PC sums, revealing also various unexpected compositions. As an example, PC 16:0_22:5 accounted for more than 50% of the PC sum with in total 38 carbon atoms and 5 double bonds (PC aa 38:5). For 13 PC sums, we found relatively high abundances of odd-chain fatty acids. In conclusion, our study provides insights in PC compositions in human plasma, facilitating interpretation of existing epidemiological data sets and potentially enabling imputation of PC compositions for future meta-analyses of lipidomics data. Show less
Kidney fibrosis is the common pathophysiological mechanism in end-stage renal disease characterized by excessive accumulation of myofibroblast-derived extracellular matrix. Natriuretic peptides... Show moreKidney fibrosis is the common pathophysiological mechanism in end-stage renal disease characterized by excessive accumulation of myofibroblast-derived extracellular matrix. Natriuretic peptides have been demonstrated to have cyclic guanosine monophosphate (cGMP)-dependent anti-fibrotic properties likely due to interference with pro-fibrotic tissue growth factor beta (TGF-beta) signaling. However, in vivo, natriuretic peptides are rapidly degraded by neutral endopeptidases (NEP). In a unilateral ureteral obstruction (UUO) mouse model for kidney fibrosis we assessed the anti-fibrotic effects of SOL1, an orally active compound that inhibits NEP and endothelin-converting enzyme (ECE). Mice (n= 10 per group) subjected to UUO were treated for 1 week with either solvent, NEP-/ECE-inhibitor SOL1 (two doses), reference NEP-inhibitor candoxatril or the angiotensin II receptor type 1 (AT1)-antagonist losartan. While NEP-inhibitors had no significant effect on blood pressure, they did increase urinary cGMP levels as well as endothelin-1 (ET-1) levels. Immunohistochemical staining revealed a marked decrease in renal collagen (similar to 55% reduction, P< 0.05) and alpha-smooth muscle actin (alpha-SMA; similar to 40% reduction, P< 0.05). Moreover, the number of alpha-SMA positive cells in the kidneys of SOL1-treated groups inversely correlated with cGMP levels consistent with a NEP-dependent anti-fibrotic effect. To dissect the molecular mechanisms associated with the anti-fibrotic effects of NEP inhibition, we performed a 'deep serial analysis of gene expression (Deep SAGE)' transcriptome and targeted metabolomics analysis of total kidneys of all treatment groups. Pathway analyses linked increased cGMP and ET-1 levels with decreased nuclear receptor signaling (peroxisome proliferator-activated receptor [PPAR] and liver X receptor/retinoid X receptor [LXR/RXR] signaling) and actin cytoskeleton organization. Taken together, although our transcriptome and metabolome data indicate metabolic dysregulation, our data support the therapeutic potential of NEP inhibition in the treatment of kidney fibrosis via cGMP elevation and reduced myofibroblast formation. Show less
Introduction Metabolic changes have been frequently associated with Huntington's disease (HD). At the same time peripheral blood represents aminimally invasive sampling avenue with little distress... Show moreIntroduction Metabolic changes have been frequently associated with Huntington's disease (HD). At the same time peripheral blood represents aminimally invasive sampling avenue with little distress to Huntington's disease patients especially when brain or other tissue samples are difficult to collect.Objectives We investigated the levels of 163 metabolites in HD patient and control serum samples in order to identify disease related changes. Additionally, we integrated the metabolomics data with our previously published next generation sequencing-based gene expression data from the same patients in order to interconnect the metabolomics changes with transcriptional alterations. Methods This analysis was performed using targeted metabolomics and flow injection electrospray ionization tandem mass spectrometry in 133 serum samples from 97 Huntington's disease patients (29 pre-symptomatic and 68 symptomatic) and 36 controls.Results By comparing HD mutation carriers with controls we identified 3 metabolites significantly changed in HD (serine and threonine and one phosphatidylcholine-PC ae C36:0) and an additional 8 phosphatidylcholines (PC aa C38:6, PC aa C36:0, PC ae C38:0, PC aa C38:0, PC ae C38:6, PC ae C42:0, PC aa C36:5 and PC ae C36:0) that exhibited a significant association with disease severity. Using workflow based exploitation of pathway databases and by integrating our metabolomics data with our gene expression data from the same patients we identified 4 deregulated phosphatidylcholine metabolism related genes (ALDH1B1, MBOAT1, MTRR and PLB1) that showed significant association with the changes in metabolite concentrations.Conclusion Our results support the notion that phosphatidylcholine metabolism is deregulated in HD blood and that these metabolite alterations are associated with specific gene expression changes. Show less
Mastrokolias, A.; Pool, R.; Mina, E.; Hettne, K.M.; Duijn, E. van; Mast, R.C. van der; ... ; Roon-Mom, W. van 2016
IntroductionMetabolic changes have been frequently associated with Huntington’s disease (HD). At the same time peripheral blood represents a minimally invasive sampling avenue with little distress... Show moreIntroductionMetabolic changes have been frequently associated with Huntington’s disease (HD). At the same time peripheral blood represents a minimally invasive sampling avenue with little distress to Huntington’s disease patients especially when brain or other tissue samples are difficult to collect.ObjectivesWe investigated the levels of 163 metabolites in HD patient and control serum samples in order to identify disease related changes. Additionally, we integrated the metabolomics data with our previously published next generation sequencing-based gene expression data from the same patients in order to interconnect the metabolomics changes with transcriptional alterations.MethodsThis analysis was performed using targeted metabolomics and flow injection electrospray ionization tandem mass spectrometry in 133 serum samples from 97 Huntington’s disease patients (29 pre-symptomatic and 68 symptomatic) and 36 controls.ResultsBy comparing HD mutation carriers with controls we identified 3 metabolites significantly changed in HD (serine and threonine and one phosphatidylcholine—PC ae C36:0) and an additional 8 phosphatidylcholines (PC aa C38:6, PC aa C36:0, PC ae C38:0, PC aa C38:0, PC ae C38:6, PC ae C42:0, PC aa C36:5 and PC ae C36:0) that exhibited a significant association with disease severity. Using workflow based exploitation of pathway databases and by integrating our metabolomics data with our gene expression data from the same patients we identified 4 deregulated phosphatidylcholine metabolism related genes (ALDH1B1, MBOAT1, MTRR and PLB1) that showed significant association with the changes in metabolite concentrations.ConclusionOur results support the notion that phosphatidylcholine metabolism is deregulated in HD blood and that these metabolite alterations are associated with specific gene expression changes. Show less
Berg, R. van den; Mook-Kanamori, D.O.; Donga, E.; Dijk, M. van; Dijk, J.G. van; Lammers, G.J.; ... ; Biermasz, N.R. 2016
