Search

In cross-platform analyses of 174 metabolites, we identify 499 associations (P

Funder: Medical Research Council, Understanding the genetic architecture of host proteins interacting with SARS-CoV-2 or mediating the maladaptive host response to COVID-19 can help to identify new or repurpose existing drugs targeting those proteins. We present a genetic discovery study of 179 such host proteins among 10,708 individuals using an aptamer-based technique. We identify 220 host DNA sequence variants acting in cis (MAF 0.01-49.9%) and explaining 0.3-70.9% of the variance of 97 of these proteins, including 45 with no previously known protein quantitative trait loci (pQTL) and 38 encoding current drug targets. Systematic characterization of pQTLs across the phenome identified protein-drug-disease links and evidence that putative viral interaction partners such as MARK3 affect immune response. Our results accelerate the evaluation and prioritization of new drug development programmes and repurposing of trials to prevent, treat or reduce adverse outcomes. Rapid sharing and detailed interrogation of results is facilitated through an interactive webserver (https://omicscience.org/apps/covidpgwas/).

Late-onset Alzheimer’s disease (AD) can, in part, be considered a metabolic disease. Besides age, female sex and APOE ε4 genotype represent strong risk factors for AD that also give rise to large metabolic differences. We systematically investigated group-specific metabolic alterations by conducting stratified association analyses of 139 serum metabolites in 1,517 individuals from the AD Neuroimaging Initiative with AD biomarkers. We observed substantial sex differences in effects of 15 metabolites with partially overlapping differences for APOE ε4 status groups. Several group-specific metabolic alterations were not observed in unstratified analyses using sex and APOE ε4 as covariates. Combined stratification revealed further subgroup-specific metabolic effects limited to APOE ε4+ females. The observed metabolic alterations suggest that females experience greater impairment of mitochondrial energy production than males. Dissecting metabolic heterogeneity in AD pathogenesis can therefore enable grading the biomedical relevance for specific pathways within specific subgroups, guiding the way to personalized medicine., Sex and the APOE ε4 genotype are important risk factors for late-onset Alzheimer’s disease. In the current study, the authors investigate how sex and APOE ε4 genotype modify the association between Alzheimer’s disease biomarkers and metabolites in serum.

Aim. Cardiac arrhythmias comprise a major health and economic burden and are associated with significant morbidity and mortality, including cardiac failure, stroke and sudden cardiac death (SCD). Development of efficient preventive and therapeutic strategies is hampered by incomplete knowledge of disease mechanisms and pathways. Our aim is to identify novel mechanisms underlying cardiac arrhythmia and SCD using an unbiased approach. Methods and Results. We employed a phenotype-driven N-ethyl-N-nitrosourea (ENU) mutagenesis screen and identified a mouse line with a high incidence of sudden death at young age (6-9 weeks) in the absence of prior symptoms. Affected mice were found to be homozygous for the nonsense mutation Bcat2p.Q300*/p.Q300* in the Bcat2 gene encoding branched chain amino acid transaminase 2. At the age of 4-5 weeks, Bcat2p.Q300*/p.Q300* mice displayed drastic increase of plasma levels of branch chain amino acids (BCAAs – leucine, isoleucine, valine) due to the incomplete catabolism of BCAAs, in addition to inducible arrhythmias ex vivo as well as cardiac conduction and repolarization disturbances. In line with these findings, plasma BCAA levels were positively correlated to ECG indices of conduction and repolarization in the German community-based KORA F4 Study. Isolated cardiomyocytes from Bcat2p.Q300*/p.Q300* mice revealed action potential (AP) prolongation, pro-arrhythmic events (early and late afterdepolarizations, triggered APs) and dysregulated calcium homeostasis. Incubation of human pluripotent stem cell-derived cardiomyocytes with elevated concentration of BCAAs induced similar calcium dysregulation and pro-arrhythmic events which were prevented by rapamycin, demonstrating the crucial involvement of mTOR pathway activation. Conclusions. Our findings identify for the first time a causative link between elevated BCAAs and arrhythmia, which has implications for arrhythmogenesis in conditions associated with BCAA metabolism dysregulation such as diabetes, metabolic syndrome and heart failure. Translational perspectives. Development of efficient anti-arrhythmic strategies is hampered by incomplete knowledge of disease mechanisms. Using an unbiased approach, we here identified for the first time a pro-arrhythmic effect of increased levels of branched chain amino acids (BCAAs). This is of particular relevance for conditions associated with BCAA dysregulation and increased arrhythmia risk, including heart failure, obesity and diabetes, as well as for athletes supplementing their diet with BCAAs. Such metabolic dysregulation is potentially modifiable through dietary interventions, paving the way for novel preventive strategies. Our findings furthermore identify mTOR inhibition as a potential anti-arrhythmic strategy in patients with metabolic syndrome.

Multimorbidity, the simultaneous presence of multiple chronic conditions, is an increasing global health problem and research into its determinants is of high priority. We used baseline untargeted plasma metabolomics profiling covering >1,000 metabolites as a comprehensive readout of human physiology to characterize pathways associated with and across 27 incident noncommunicable diseases (NCDs) assessed using electronic health record hospitalization and cancer registry data from over 11,000 participants (219,415 person years). We identified 420 metabolites shared between at least 2 NCDs, representing 65.5% of all 640 significant metabolite-disease associations. We integrated baseline data on over 50 diverse clinical risk factors and characteristics to identify actionable shared pathways represented by those metabolites. Our study highlights liver and kidney function, lipid and glucose metabolism, low-grade inflammation, surrogates of gut microbial diversity and specific health-related behaviors as antecedents of common NCD multimorbidity with potential for early prevention. We integrated results into an open-access webserver ( https://omicscience.org/apps/mwasdisease/ ) to facilitate future research and meta-analyses., We are grateful to all the participants who have been part of the project and to the many members of the study teams at the University of Cambridge who have enabled this research. The EPIC-Norfolk study (https://doi.org/10.22025/2019.10.105.00004) has received funding from the Medical Research Council (MR/N003284/1 and MC-UU_12015/1) and Cancer Research UK (C864/A14136). Metabolite measurements in the EPIC-Norfolk study were supported by the MRC Cambridge Initiative in Metabolic Science (MR/L00002/1) and the Innovative Medicines Initiative Joint Undertaking under EMIF grant agreement no. 115372. MP was supported by a fellowship of the German Research Foundation (DFG 1446/2-1). JR is supported by the German Federal Ministry of Education and Research (BMBF) within the framework of the e:Med research and funding concept (grant no. 01ZX1912D).

Biological exploration of early biomarkers for chronic kidney disease (CKD) in (pre)diabetic individuals is crucial for personalized management of diabetes. Here, we evaluated two candidate biomarkers of incident CKD (sphingomyelin (SM) C18:1 and phosphatidylcholine diacyl (PC aa) C38:0) concerning kidney function in hyperglycemic participants of the Cooperative Health Research in the Region of Augsburg (KORA) cohort, and in two biofluids and six organs of leptin receptor-deficient (db/db) mice and wild type controls. Higher serum concentrations of SM C18:1 and PC aa C38:0 in hyperglycemic individuals were found to be associated with lower estimated glomerular filtration rate (eGFR) and higher odds of CKD. In db/db mice, both metabolites had a significantly lower concentration in urine and adipose tissue, but higher in the lungs. Additionally, db/db mice had significantly higher SM C18:1 levels in plasma and liver, and PC aa C38:0 in adrenal glands. This cross-sectional human study confirms that SM C18:1 and PC aa C38:0 associate with kidney dysfunction in pre(diabetic) individuals, and the animal study suggests a potential implication of liver, lungs, adrenal glands, and visceral fat in their systemic regulation. Our results support further validation of the two phospholipids as early biomarkers of renal disease in patients with (pre)diabetes.

Background Periodontitis is among the most common chronic diseases worldwide, and it is one of the main reasons for tooth loss. Comprehensive profiling of the metabolite content of the saliva can enable the identification of novel pathways associated with periodontitis and highlight non-invasive markers to facilitate time and cost-effective screening efforts for the presence of periodontitis and the prediction of tooth loss. Methods We first investigated cross-sectional associations of 13 oral health variables with saliva levels of 562 metabolites, measured by untargeted mass spectrometry among a sub-sample (n = 938) of the Study of Health in Pomerania (SHIP-2) using linear regression models adjusting for common confounders. We took forward any candidate metabolite associated with at least two oral variables, to test for an association with a 5-year tooth loss over and above baseline oral health status using negative binomial regression models. Results We identified 84 saliva metabolites that were associated with at least one oral variable cross-sectionally, for a subset of which we observed robust replication in an independent study. Out of 34 metabolites associated with more than two oral variables, baseline saliva levels of nine metabolites were positively associated with a 5-year tooth loss. Across all analyses, the metabolites 2-pyrrolidineacetic acid and butyrylputrescine were the most consistent candidate metabolites, likely reflecting oral dysbiosis. Other candidate metabolites likely reflected tissue destruction and cell proliferation. Conclusions Untargeted metabolic profiling of saliva replicated metabolic signatures of periodontal status and revealed novel metabolites associated with periodontitis and future tooth loss.

Metabolomics in the Alzheimer’s Disease Neuroimaging Initiative cohort provides a powerful tool for mapping biochemical changes in Alzheimer’s disease, and a unique opportunity to learn about the association between circulating blood metabolites and brain amyloid-β deposition in Alzheimer’s disease. We examined 140 serum metabolites and their associations with brain amyloid-β deposition, cognition and conversion from mild cognitive impairment to Alzheimer’s disease in the Alzheimer’s Disease Neuroimaging Initiative. Processed [18F] Florbetapir PET images were used to perform a voxel-wise statistical analysis of the effect of metabolite levels on amyloid-β accumulation across the whole brain. We performed a multivariable regression analysis using age, sex, body mass index, apolipoprotein E ε4 status and study phase as covariates. We identified nine metabolites as significantly associated with amyloid-β deposition after multiple comparison correction. Higher levels of one acylcarnitine (C3; propionylcarnitine) and one biogenic amine (kynurenine) were associated with decreased amyloid-β accumulation and higher memory scores. However, higher levels of seven phosphatidylcholines (lysoPC a C18:2, PC aa C42:0, PC ae C42:3, PC ae C44:3, PC ae C44:4, PC ae C44:5 and PC ae C44:6) were associated with increased brain amyloid-β deposition. In addition, higher levels of PC ae C44:4 were significantly associated with lower memory and executive function scores and conversion from mild cognitive impairment to Alzheimer’s disease dementia. Our findings suggest that dysregulation of peripheral phosphatidylcholine metabolism is associated with earlier pathological changes noted in Alzheimer’s disease as measured by brain amyloid-β deposition as well as later clinical features including changes in memory and executive functioning. Perturbations in phosphatidylcholine metabolism may point to issues with membrane restructuring leading to the accumulation of amyloid-β in the brain. Additional studies are needed to explore whether these metabolites play a causal role in the pathogenesis of Alzheimer’s disease or if they are biomarkers for systemic changes during preclinical phases of the disease., Nho et al. report that serum-based metabolites are associated with brain amyloid-β deposition, cognition and disease progression from mild cognitive impairment to Alzheimer’s disease in a large cohort of older adults. Results suggest that dysregulation of very long chain peripheral phospholipid metabolism is associated with earlier pathological changes in Alzheimer’s disease., Graphical Abstract Graphical Abstract

Shared metabolomic patterns at delivery have been suggested to underlie the mother-to-child transmission of adverse metabolic health. This study aimed to investigate whether mothers with gestational diabetes mellitus (GDM) and their offspring show similar metabolomic patterns several years postpartum. Targeted metabolomics (including 137 metabolites) was performed in plasma samples obtained during an oral glucose tolerance test from 48 mothers with GDM and their offspring at a cross-sectional study visit 8 years after delivery. Partial Pearson&rsquo, s correlations between the area under the curve (AUC) of maternal and offspring metabolites were calculated, yielding so-called Gaussian graphical models. Spearman&rsquo, s correlations were applied to investigate correlations of body mass index (BMI), Matsuda insulin sensitivity index (ISI-M), dietary intake, and physical activity between generations, and correlations of metabolite AUCs with lifestyle variables. This study revealed that BMI, ISI-M, and the AUC of six metabolites (carnitine, taurine, proline, SM(-OH) C14:1, creatinine, and PC ae C34:3) were significantly correlated between mothers and offspring several years postpartum. Intergenerational metabolite correlations were independent of shared BMI, ISI-M, age, sex, and all other metabolites. Furthermore, creatinine was correlated with physical activity in mothers. This study suggests that there is long-term metabolic programming in the offspring of mothers with GDM and informs us about targets that could be addressed by future intervention studies.

Metabolomics can be a tool to identify dietary biomarkers. However, reported food-metabolite associations have been inconsistent, and there is a need to explore further associations. Our aims were to confirm previously reported food-metabolite associations and to identify novel food-metabolite associations. We conducted a cross-sectional analysis of data from 849 participants (57% men) of the PopGen cohort. Dietary intake was obtained using FFQ and serum metabolites were profiled by an untargeted metabolomics approach. We conducted a systematic literature search to identify previously reported food-metabolite associations and analyzed these associations using linear regression. To identify potential novel food-metabolite associations, datasets were split into training and test datasets and linear regression models were fitted to the training datasets. Significant food-metabolite associations were evaluated in the test datasets. Models were adjusted for covariates. In the literature, we identified 82 food-metabolite associations. Of these, 44 associations were testable in our data and confirmed associations of coffee with 12 metabolites, of fish with five, of chocolate with two, of alcohol with four, and of butter, poultry and wine with one metabolite each. We did not identify novel food-metabolite associations, however, some associations were sex-specific. Potential use of some metabolites as biomarkers should consider sex differences in metabolism.

Early and precise identification of individuals with prediabetes and type 2 diabetes (T2D) at risk for progressing to chronic kidney disease (CKD) is essential to prevent complications of diabetes. Here, we identify and evaluate prospective metabolite biomarkers and the best set of predictors of CKD in the longitudinal, population-based Cooperative Health Research in the Region of Augsburg (KORA) cohort by targeted metabolomics and machine learning approaches. Out of 125 targeted metabolites, sphingomyelin C18:1 and phosphatidylcholine diacyl C38:0 were identified as candidate metabolite biomarkers of incident CKD specifically in hyperglycemic individuals followed during 6.5 years. Sets of predictors for incident CKD developed from 125 metabolites and 14 clinical variables showed highly stable performances in all three machine learning approaches and outperformed the currently established clinical algorithm for CKD. The two metabolites in combination with five clinical variables were identified as the best set of predictors, and their predictive performance yielded a mean area value under the receiver operating characteristic curve of 0.857. The inclusion of metabolite variables in the clinical prediction of future CKD may thus improve the risk prediction in people with prediabetes and T2D. The metabolite link with hyperglycemia-related early kidney dysfunction warrants further investigation.

Saliva, as a biofluid, is inexpensive and non-invasive to obtain, and provides a vital tool to investigate oral health and its interaction with systemic health conditions. There is growing interest in salivary biomarkers for systemic diseases, notably cardiovascular disease. Whereas hundreds of genetic loci have been shown to be involved in the regulation of blood metabolites leading to unprecedented insights into the pathogenesis of complex human diseases, little is known about the impact of host genetics on salivary metabolites. Here we report the first genome-wide association study exploring 476 salivary metabolites in 1,419 subjects of European ancestry from the TwinsUK cohort (discovery phase). A total of 14 salivary metabolites were significantly associated (p−10) with genetic variants that mapped to 11 distinct loci, most of which replicated in the Study of Health in Pomerania (SHIP-2) cohort. Interestingly, while only a limited number of the loci that are known to regulate blood metabolites were also associated with salivary metabolites in our study, we identified several novel saliva-specific locus-metabolite associations, including associations for theAGMAT(with the metabolites 4-guanidinobutanoate and beta-guanidinopropanoate),ATP13A5(with the metabolite creatinine) andDPYS(with the metabolites 3-ureidopropionate and 3-ureidoisobutyrate) loci. Our study suggests that there are biological pathways which are specific to the regulation of the salivary metabolome. In addition, some of our findings may have clinical relevance, such as the utility of the pyrimidine (uracil) degradation metabolites in predicting 5-fluorouracil toxicity and the role of the agmatine pathway metabolites as biomarkers of oral health.

Obesity is one of the major risk factor for cardiovascular and metabolic diseases. A disproportional accumulation of fat at visceral (VAT) compared to subcutaneous sites (SAT) has been suspected as a key detrimental event. We used non-targeted metabolomics profiling to reveal metabolic pathways associated with higher VAT or SAT amount among subjects free of metabolic diseases to identify possible contributing metabolic pathways. The study population comprised 491 subjects [mean (standard deviation): age 44.6 yrs (13.0), body mass index 25.4 kg/m² (3.6), 60.1% females] without diabetes, hypertension, dyslipidemia, the metabolic syndrome or impaired renal function. We associated MRI-derived fat amounts with mass spectrometry-derived metabolites in plasma and urine using linear regression models adjusting for major confounders. We tested for sex-specific effects using interactions terms and performed sensitivity analyses for the influence of insulin resistance on the results. VAT and SAT were significantly associated with 155 (101 urine) and 49 (29 urine) metabolites, respectively, of which 45 (27 urine) were common to both. Major metabolic pathways were branched-chain amino acid metabolism (partially independent of insulin resistance), surrogate markers of oxidative stress and gut microbial diversity, and cortisol metabolism. We observed a novel positive association between VAT and plasma levels of the potential pharmacological agent piperine. Sex-specific effects were only a few, e.g. the female-specific association between VAT and O-methylascorbate. In brief, higher VAT was associated with an unfavorable metabolite profile in a sample of healthy, mostly non-obese individuals from the general population and only few sex-specific associations became apparent.

Correlation networks are frequently used to statistically extract biological interactions between omics markers. Network edge selection is typically based on the statistical significance of the correlation coefficients. This procedure, however, is not guaranteed to capture biological mechanisms. We here propose an alternative approach for network reconstruction: a cutoff selection algorithm that maximizes the overlap of the inferred network with available prior knowledge. We first evaluate the approach on IgG glycomics data, for which the biochemical pathway is known and well-characterized. Importantly, even in the case of incomplete or incorrect prior knowledge, the optimal network is close to the true optimum. We then demonstrate the generalizability of the approach with applications to untargeted metabolomics and transcriptomics data. For the transcriptomics case, we demonstrate that the optimized network is superior to statistical networks in systematically retrieving interactions that were not included in the biological reference used for optimization., Correlation network inference is typically based on the significance of the correlation coefficients, but this procedure is not guaranteed to capture biological mechanisms. Here, the authors develop a cutoff selection algorithm that maximizes the overlap between inferred networks and prior knowledge.

ObjectiveWhile 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.MethodsProspective 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.ResultsMetabolomics/lipidomics analysis of differential concentrations (false discovery rateConclusionsIn PAH, trans-RV and transpulmonary metabolite gradients exist and correlate with haemodynamic determinants of clinical outcome. The most pronounced differential trans-RV 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 β-oxidation in the hypertensive RV and represents emerging biomarkers and therapeutic targets in PAH.

Selenoprotein P (SELENOP) has been previously related to various metabolic traits with partially conflicting results. The identification of SELENOP-associated metabolites, using an untargeted metabolomics approach, may provide novel biological insights relevant to disentangle the role of SELENOP in human health. In this cross-sectional study, 572 serum metabolites were identified by comparing the obtained LC–MS/MS spectra with spectra stored in Metabolon’s spectra library. Serum SELENOP levels were measured in 832 men and women using an ELISA kit. Circulating SELENOP levels were associated with 24 out of 572 metabolites after accounting for the number of independent dimensions in the metabolomics data, including inverse associations with alanine, glutamate, leucine, isoleucine and valine, an unknown compound X-12063, urate and the peptides gamma-glutamyl-leucine, and N-acetylcarnosine. Positive associations were observed between SELENOP and several lipid compounds. Of the identified metabolites, each standard deviation increase in the branched-chain amino acids (isoleucine, leucine, valine), alanine and gamma-glutamyl-leucine was related to higher odds of having T2DM [OR (95% CI): 1.96 (1.41–2.73); 1.62 (1.15–2.28); 1.94 (1.45–2.60), 1.57 (1.17–2.11), and 1.52 (1.13–2.05), respectively]. Higher serum SELENOP levels were associated with an overall healthy metabolomics profile, which may provide further insights into potential mechanisms of SELENOP-associated metabolic disorders.

Objective Thyroid hormones are ubiquitously involved in human metabolism. However, the precise molecular patterns associated with alterations in thyroid hormones levels remain to be explored in detail. A number of recent studies took great advantage of metabolomics profiling to outline the metabolic actions of thyroid hormones in humans. Methods Among 952 participants in the Study of Health in Pomerania, data on serum free thyroxine (FT4) and thyrotropin and comprehensive nontargeted metabolomics data from plasma and urine samples were available. Linear regression analyses were performed to assess the association between FT4 or thyrotropin and metabolite levels. Results and Conclusion After accounting for major confounders, 106 of 613 plasma metabolites were significantly associated with FT4. The associations in urine were minor (12 of 587). Most of the plasma metabolites consisted of lipid species, and subsequent analysis of highly resolved lipoprotein subclasses measured by proton nuclear magnetic resonance spectroscopy revealed a consistent decrease in several of these species (e.g., phospholipids) and large low-density lipoprotein and small high-density lipoprotein particles. The latter was unique to men. Several polyunsaturated and saturated fatty acids displayed an association with FT4 in women only. A random forest-based variable selection approach using phenotypic characteristics revealed higher alcohol intake in men and an adverse thyroid state and menopause in women as the putative mediating factors. In general, our observations have confirmed the lipolytic and lipogenic effect of thyroid hormones even in the physiological range and revealed different phenotypic characteristics (e.g., lifestyle differences) as possible confounders for sex-specific findings.

Metabolites are small molecules that are intermediates or products of metabolism, many of which are freely filtered by the kidneys. In addition, the kidneys have a central role in metabolite anabolism and catabolism, as well as in active metabolite reabsorption and/or secretion during tubular passage. This review article illustrates how the coupling of genomics and metabolomics in genome-wide association analyses of metabolites can be used to illuminate mechanisms underlying human metabolism, with a special focus on insights relevant to nephrology. First, genetic susceptibility loci for reduced kidney function and chronic kidney disease (CKD) were reviewed systematically for their associations with metabolite concentrations in metabolomics studies of blood and urine. Second, kidney function and CKD-associated metabolites reported from observational studies were interrogated for metabolite-associated genetic variants to generate and discuss complementary insights. Finally, insights originating from the simultaneous study of both blood and urine or by modeling intermetabolite relationships are summarized. We also discuss methodologic questions related to the study of metabolite concentrations in urine as well as among CKD patients. In summary, genome-wide association analyses of metabolites using metabolite concentrations quantified from blood and/or urine are a promising avenue of research to illuminate physiological and pathophysiological functions of the kidney.

Using oral contraceptives has been implicated in the aetiology of stress-related disorders like depression. Here, we followed the hypothesis that oral contraceptives deregulate the HPA-axis by elevating circulating cortisol levels. We report for a sample of 233 pre-menopausal women increased circulating cortisol levels in those using oral contraceptives. For women taking oral contraceptives, we observed alterations in circulating phospholipid levels and elevated triglycerides and found evidence for increased glucocorticoid signalling as the transcript levels of the glucocorticoid-regulated genes DDIT4 and FKBP5 were increased in whole blood. The effects were statistically mediated by cortisol. The associations of oral contraceptives with higher FKBP5 mRNA and altered phospholipid levels were modified by rs1360780, a genetic variance implicated in psychiatric diseases. Accordingly, the methylation pattern of FKBP5 intron 7 was altered in women taking oral contraceptives depending on the rs1360780 genotype. Moreover, oral contraceptives modified the association of circulating cortisol with depressive symptoms, potentially explaining conflicting results in the literature. Finally, women taking oral contraceptives displayed smaller hippocampal volumes than non-using women. In conclusion, the integrative analyses of different types of physiological data provided converging evidence indicating that oral contraceptives may cause effects analogous to chronic psychological stressors regarding the regulation of the HPA axis.

Alzheimer’s disease (AD) is the most common cause of dementia. The mechanism of disease development and progression is not well understood, but increasing evidence suggests multifactorial etiology, with a number of genetic, environmental, and aging-related factors. There is a growing body of evidence that metabolic defects may contribute to this complex disease. To interrogate the relationship between system level metabolites and disease susceptibility and progression, the AD Metabolomics Consortium (ADMC) in partnership with AD Neuroimaging Initiative (ADNI) is creating a comprehensive biochemical database for patients in the ADNI1 cohort. We used the Biocrates Bile Acids platform to evaluate the association of metabolic levels with disease risk and progression. We detail the quantitative metabolomics data generated on the baseline samples from ADNI1 and ADNIGO/2 (370 cognitively normal, 887 mild cognitive impairment, and 305 AD). Similar to our previous reports on ADNI1, we present the tools for data quality control and initial analysis. This data descriptor represents the third in a series of comprehensive metabolomics datasets from the ADMC on the ADNI., Measurement(s)bile acid • MedicationTechnology Type(s)ultra high performance liquid chromatography with mass spectrometer • Resource InformaticsFactor Type(s)age • biological sex • cognitive stateSample Characteristic - OrganismHomo sapiens Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.9724652

Alzheimer’s disease (AD) is the leading cause of dementia, with metabolic dysfunction seen years before the emergence of clinical symptoms. Increasing evidence suggests a role for primary and secondary bile acids, the end-product of cholesterol metabolism, influencing pathophysiology in AD. In this study, we analyzed transcriptomes from 2114 post-mortem brain samples from three independent cohorts and identified that the genes involved in alternative bile acid synthesis pathway was expressed in brain compared to the classical pathway. These results were supported by targeted metabolomic analysis of primary and secondary bile acids measured from post-mortem brain samples of 111 individuals. We reconstructed brain region-specific metabolic networks using data from three independent cohorts to assess the role of bile acid metabolism in AD pathophysiology. Our metabolic network analysis suggested that taurine transport, bile acid synthesis and cholesterol metabolism differed in AD and cognitively normal individuals. Using the brain transcriptional regulatory network, we identified putative transcription factors regulating these metabolic genes and influencing altered metabolism in AD. Intriguingly, we find bile acids from the brain metabolomics whose synthesis cannot be explained by enzymes we find in the brain, suggesting they may originate from an external source such as the gut microbiome. These findings motivate further research into bile acid metabolism and transport in AD to elucidate their possible connection to cognitive decline.

Kit-based assays, such as AbsoluteIDQ™ 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™ platform, analyzing 223 samples in parallel to the AbsoluteIDQ™. 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 datasets and potentially enabling imputation of PC compositions for future meta-analyses of lipidomics data.

Recent studies have provided evidence that late-onset Alzheimer’s disease (AD) can in part be considered a metabolic disease. Besides age, female sex andAPOEε4 genotype represent strong risk factors for AD. They also both give rise to large metabolic differences, suggesting that metabolic aspects of AD pathogenesis may differ between males and females and betweenAPOEε4 carriers and non-carriers. We systematically investigated group-specific metabolic alterations by conducting stratified association analyses of 140 metabolites measured in serum samples of 1,517 AD neuroimaging initiative subjects, with AD biomarkers for Aβ and tau pathology and neurodegeneration. We observed substantial sex differences in effects of 15 metabolites on AD biomarkers with partially overlapping differences forAPOEε4 status groups. These metabolites highlighted several group-specific alterations not observed in unstratified analyses using sex andAPOEε4 as covariates. Combined stratification by both variables uncovered further subgroup-specific metabolic effects limited to the group with presumably the highest AD risk:APOEε4+ females. Pathways linked to the observed metabolic alterations suggest that females experience more expressed impairment of mitochondrial energy production in AD than males. These findings indicate that dissecting metabolic heterogeneity in AD pathogenesis may enable grading of the biomedical relevance of specific pathways for specific subgroups. Extending our approach beyond simple one- or two-fold stratification may thus guide the way to personalized medicine.Significance statementResearch provides substantial evidence that late-onset Alzheimer’s disease (AD) is a metabolic disease. Besides age, female sex and APOEε4 genotype represent strong risk factors for AD, and at the same time give rise to large metabolic differences. Our systematic investigation of sex and APOE ε4 genotype differences in the link between metabolism and measures of pre-symptomatic AD using stratified analysis revealed several group-specific metabolic alterations that were not observed without sex and genotype stratification of the same cohort. Pathways linked to the observed metabolic alterations suggest females are more affected by impairment of mitochondrial energy production in AD than males, highlighting the importance of tailored treatment approaches towards a precision medicine approach.

Determination of metabolomic signatures of pulmonary function and chronic obstructive pulmonary disease (COPD) in the general population could aid in identification and understanding of early disease processes. Metabolome measurements were performed on serum from 4742 individuals (2354 African-Americans and 1529 European-Americans from the Atherosclerosis Risk in Communities study and 859 Europeans from the Cooperative Health Research in the Region of Augsburg study). We examined 368 metabolites in relation to cross-sectional measures of forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC), their ratio (FEV1/FVC) and COPD using multivariable regression followed by meta-analysis. At a false discovery rate of 0.05, 95 metabolites were associated with FEV1 and 100 with FVC (73 overlapping), including inverse associations with branched-chain amino acids and positive associations with glutamine. Ten metabolites were associated with FEV1/FVC and seventeen with COPD (393 cases). Enriched pathways of amino acid metabolism were identified. Associations with FEV1 and FVC were not driven by individuals with COPD. We identified novel metabolic signatures of pulmonary function and COPD in African and European ancestry populations. These may allow development of biomarkers in the general population of early disease pathogenesis, before pulmonary function has decreased to levels diagnostic for COPD.

The kidneys integrate information from continuous systemic processes related to the absorption, distribution, metabolism and excretion (ADME) of metabolites. To identify underlying molecular mechanisms, we performed genome-wide association studies of the urinary concentrations of 1,172 metabolites among 1,627 patients with reduced kidney function. The 240 unique metabolite–locus associations (metabolite quantitative trait loci, mQTLs) that were identified and replicated highlight novel candidate substrates for transport proteins. The identified genes are enriched in ADME-relevant tissues and cell types, and they reveal novel candidates for biotransformation and detoxification reactions. Fine mapping of mQTLs and integration with single-cell gene expression permitted the prioritization of causal genes, functional variants and target cell types. The combination of mQTLs with genetic and health information from 450,000 UK Biobank participants illuminated metabolic mediators, and hence, novel urinary biomarkers of disease risk. This comprehensive resource of genetic targets and their substrates is informative for ADME processes in humans and is relevant to basic science, clinical medicine and pharmaceutical research. Genome-wide association analysis of 1,172 urinary metabolites identifies 240 metabolite–locus associations that when combined with UK Biobank data suggest novel metabolic mediators of disease and markers of disease risk.

BACKGROUND: Numerous observational studies have observed associations between vitamin D deficiency and cardiometabolic diseases, but these findings might be confounded by obesity. A characterization of the metabolic profile associated with serum 25-hydroxyvitamin D [25(OH)D] levels, in general and stratified by abdominal obesity, may help to untangle the relationship between vitamin D, obesity and cardiometabolic health. METHODS: Serum metabolomics measurements were obtained from a nuclear magnetic resonance spectroscopy (NMR)- and a mass spectrometry (MS)-based platform. The discovery was conducted in 1726 participants of the population-based KORA-F4 study, in which the associations of the concentrations of 415 metabolites with 25(OH)D levels were assessed in linear models. The results were replicated in 6759 participants (NMR) and 609 (MS) participants, respectively, of the population-based FINRISK 1997 study. RESULTS: Mean [standard deviation (SD)] 25(OH)D levels were 15.2 (7.5) ng/ml in KORA F4 and 13.8 (5.9) ng/ml in FINRISK 1997; 37 metabolites were associated with 25(OH)D in KORA F4 at P&thinsp

This study aims at identifying metabolites that significantly associate with self-reported joint symptoms (diagnostic) and metabolites that can predict the change from a symptom-free status to the development of self-reported joint symptoms after a 7 years period (prognostic). More than 300 metabolites were analyzed for 2246 subjects from the longitudinal study of the KORA (Cooperative Health Research in the Region of Augsburg, Germany), specifically the fourth survey S4 and its 7-year follow-up study F4. Two types of self-reported symptoms, chronic joint inflammation and worn out joints, were used for the analyses. Diagnostic analysis identified dysregulated metabolites in cases with symptoms compared with controls. Prognostic analysis identified metabolites that differentiate subjects in S4 who remained symptom-free after 7 years (F4) from those who developed any combination of symptoms. 48 metabolites were identified as nominally significantly (p0.05) associated with the self-reported symptoms in the diagnostic analysis, among which steroids show Bonferroni significance. 45 metabolites were identified as nominally significantly associated with developing symptoms after 7 years, among which hippurate showed Bonferroni significance. We show that metabolic profiles of self-reported joint symptoms are in line with metabolites known to associate with various forms of arthritis and suggest that future studies may benefit from that by investigating the possible use of self-reporting/questionnaire along with metabolic markers for the early referral of patients for further diagnostic workup and treatment of arthritis.

Background: In numerous studies based predominantly on rodent models, administration of 3,5-diiodo-L-thyronine (3,5-T2), a metabolite of the thyroid hormones (TH) thyroxine (T4) and triiodo-L-thyronine (T3), was reported to cause beneficial health effects, including reversal of steatohepatosis and prevention of insulin resistance, in most instances without adverse thyrotoxic side effects. However, the empirical evidence concerning the physiological relevance of endogenously produced 3,5-T2 in humans is comparatively poor. Therefore, to improve the understanding of 3,5-T2-related metabolic processes, we performed a comprehensive metabolomic study relating serum 3,5-T2 concentrations to plasma and urine metabolite levels within a large general population sample. Methods: Serum 3,5-T2 concentrations were determined for 856 participants of the population-based Study of Health in Pomerania-TREND (SHIP-TREND). Plasma and urine metabolome data were generated using mass spectrometry and nuclear magnetic resonance spectroscopy, allowing quantification of 613 and 578 metabolites in plasma and urine, respectively. To detect thyroid function-independent significant 3,5-T2—metabolite associations, linear regression analyses controlling for major confounders, including thyrotropin and free T4, were performed. The same analyses were carried out using a sample of 16 male healthy volunteers treated for 8 weeks with 250 μg/day levothyroxine to induce thyrotoxicosis. Results: The specific molecular fingerprint of 3,5-T2 comprised 15 and 73 significantly associated metabolites in plasma and urine, respectively. Serum 3,5-T2 concentrations were neither associated with classical thyroid function parameters nor altered during experimental thyrotoxicosis. Strikingly, many metabolites related to coffee metabolism, including caffeine and paraxanthine, formed the clearest positively associated molecular signature. Importantly, these associations were replicated in the experimental human thyrotoxicosis model. Conclusion: The molecular fingerprint of 3,5-T2 demonstrates a clear and strong positive association of the serum levels of this TH metabolite with plasma levels of compounds indicating coffee consumption, therefore pointing to the liver as an organ, the metabolism of which is strongly affected by coffee. Furthermore, 3,5-T2 serum concentrations were found not to be directly TH dependent. Considering the beneficial health effects of 3,5-T2 administration observed in animal models and those of coffee consumption demonstrated in large epidemiological studies, one might speculate that coffee-stimulated hepatic 3,5-T2 production or accumulation represents an important molecular link in this connection.

Background: Genome-wide association studies have identified hundreds of loci that influence a wide variety of complex human traits; however, little is known regarding the biological mechanism of action of these loci. The recent accumulation of functional genomics ("omics"), including metabolomics data, has created new opportunities for studying the functional role of specific changes in the genome. Functional genomic data are characterized by their high dimensionality, the presence of (strong) statistical dependency between traits, and, potentially, complex genetic control. Therefore, the analysis of such data requires specific statistical genetics methods. Results: To facilitate our understanding of the genetic control of omics phenotypes, we propose a trait-centered, network-based conditional genetic association (cGAS) approach for identifying the direct effects of genetic variants on omics-based traits. For each trait of interest, we selected from a biological network a set of other traits to be used as covariates in the cGAS. The network can be reconstructed either from biological pathway databases (a mechanistic approach) or directly from the data, using a Gaussian graphical model applied to the metabolome (a data-driven approach). We derived mathematical expressions that allow comparison of the power of univariate analyses with conditional genetic association analyses. We then tested our approach using data from a population-based Cooperative Health Research in the region of Augsburg (KORA) study (n = 1,784 subjects, 1.7 million single-nucleotide polymorphisms) with measured data for 151 metabolites. Conclusions: We found that compared to single-trait analysis, performing a genetic association analysis that includes biologically relevant covariates can either gain or lose power, depending on specific pleiotropic scenarios, for which we provide empirical examples. In the context of analyzed metabolomics data, the mechanistic network approach had more power compared to the data-driven approach. Nevertheless, we believe that our analysis shows that neither a prior-knowledge-only approach nor a phenotypic-data-only approach is optimal, and we discuss possibilities for improvement.

IntroductionIncreasing evidence suggests a role for the gut microbiome in central nervous system disorders and specific role for the gut-brain axis in neurodegeneration. Bile acids (BA), products of cholesterol metabolism and clearance, are produced in the liver and are further metabolized by gut bacteria. They have major regulatory and signaling functions and seem dysregulated in Alzheimer disease (AD).MethodsSerum levels of 15 primary and secondary BAs and their conjugated forms were measured in 1,464 subjects including 370 cognitively normal older adults (CN), 284 with early mild cognitive impairment (MCI), 505 with late MCI, and 305 AD cases enrolled in the AD Neuroimaging Initiative. We assessed associations of BA profiles including selected ratios with diagnosis, cognition, and AD-related genetic variants, adjusting for cofounders and multiple testing.ResultsIn AD compared to CN, we observed significantly lower serum concentrations of a primary BA (cholic acid CA) and increased levels of the bacterially produced, secondary BA, deoxycholic acid (DCA), and its glycine and taurine conjugated forms. An increased ratio of DCA:CA, which reflects 7α-dehydroxylation of CA by gut bacteria, strongly associated with cognitive decline, a finding replicated in serum and brain samples in the Rush Religious Orders and Memory and Aging Project. Several genetic variants in immune response related genes implicated in AD showed associations with BA profiles.ConclusionWe report for the first time an association between altered BA profile, genetic variants implicated in AD and cognitive changes in disease using a large multicenter study. These findings warrant further investigation of gut dysbiosis and possible role of gut liver brain axis in the pathogenesis of AD.

BACKGROUNDUntargeted mass spectrometry (MS)-based metabolomics data often contain missing values that reduce statistical power and can introduce bias in epidemiological studies. However, a systematic assessment of the various sources of missing values and strategies to handle these data has received little attention. Missing data can occur systematically, e.g. from run day-dependent effects due to limits of detection (LOD); or it can be random as, for instance, a consequence of sample preparation.METHODSWe investigated patterns of missing data in an MS-based metabolomics experiment of serum samples from the German KORA F4 cohort (n = 1750). We then evaluated 31 imputation methods in a simulation framework and biologically validated the results by applying all imputation approaches to real metabolomics data. We examined the ability of each method to reconstruct biochemical pathways from data-driven correlation networks, and the ability of the method to increase statistical power while preserving the strength of established genetically metabolic quantitative trait loci.RESULTSRun day-dependent LOD-based missing data accounts for most missing values in the metabolomics dataset. Although multiple imputation by chained equations (MICE) performed well in many scenarios, it is computationally and statistically challenging. K-nearest neighbors (KNN) imputation on observations with variable pre-selection showed robust performance across all evaluation schemes and is computationally more tractable.CONCLUSIONMissing data in untargeted MS-based metabolomics data occur for various reasons. Based on our results, we recommend thatKNN-based imputation is performed on observations with variable pre-selection since it showed robust results in all evaluation schemes.Key messagesUntargeted MS-based metabolomics data show missing values due to both batch-specific LOD-based and non-LOD-based effects.Statistical evaluation of multiple imputation methods was conducted on both simulated and real datasets.Biological evaluation on real data assessed the ability of imputation methods to preserve statistical inference of biochemical pathways and correctly estimate effects of genetic variants on metabolite levels.KNN-based imputation on observations with variable pre-selection andK= 10 showed robust performance for all data scenarios across all evaluation schemes.

Background The N-terminal prohormone of brain natriuretic peptide (NT-proBNP) is an important biomarker for the diagnosis of heart failure. Apart from this and only recently recognized, NT-proBNP levels associate with higher HDL- and lower LDL-cholesterol levels comprising a favorable blood lipid profile. To further examine this observation, the lipoprotein profile in relation to NT-proBNP was examined in-depth by proton nuclear magnetic resonance spectroscopy (1H-NMR). We complemented this investigation with a state-of-the-art untargeted metabolomics approach. Methods Lipoprotein particles were determined by 1H-NMR spectroscopy in 872 subjects without self-reported diabetes from the population-based Study of Health in Pomerania (SHIP)-TREND with available NT-proBNP measurements. Comprehensive metabolomics data for plasma and urine samples were obtained. Linear regression models were performed to assess the associations between serum concentrations of NT-proBNP and the metabolites/lipoprotein particles measured in plasma or urine. Results An increase in serum NT-proBNP was associated with a benefical lipoprotein profile, including a decrease in VLDL, IDL and LDL-particles along with an increase in large HDL particles. These findings were replicated in a second independent cohort. Serum concentrations of NT-proBNP showed significant inverse associations with seven plasma metabolites while associations with 39 urinary metabolites, mostly comprising amino acids and related intermediates, were identified. Mediation analyses revealed adiponection as mediating factor for the associations observed with lipoproteins particles. Conclusions Most of the metabolic changes associated with NT-proBNP implicate significant influence on the blood lipid profile besides vasodilatory and the diuretic action of BNP signaling. Our data suggest that the more favorable lipoprotein profile as associated with elevated NT-proBNP concentrations in mainly cardiac healthy individuals might relate to adiponectin signaling indicating even indirect cardio-protective effects for NT-proBNP. Electronic supplementary material The online version of this article (10.1186/s12933-018-0765-1) contains supplementary material, which is available to authorized users.

IntroductionBile acids (BAs) are the end products of cholesterol metabolism produced by human and gut microbiome co-metabolism. Recent evidence suggests gut microbiota influence pathological features of Alzheimer’s disease (AD) including neuroinflammation and amyloid-β deposition.MethodSerum levels of 20 primary and secondary BA metabolites from the AD Neuroimaging Initiative (n=1562) were measured using targeted metabolomic profiling. We assessed the association of BAs with the “A/T/N” (Amyloid, Tau and Neurodegeneration) biomarkers for AD: CSF biomarkers, atrophy (MRI), and brain glucose metabolism ([18F]FDG-PET).ResultsOf 23 BA and relevant calculated ratios, three BA signatures were associated with CSF Aβ1-42 (“A”) and three with CSF p-tau181 (“T”) (corrected ppConclusionThis is the first study to show serum-based BA metabolites are associated with “A/T/N” AD biomarkers, providing further support for a role of BA pathways in AD pathophysiology. Prospective clinical observations and validation in model systems are needed to assess causality and specific mechanisms underlying this association.

Exposure to an intrauterine hyperglycaemic environment has been suggested to increase the offspring’s later risk for being overweight or having metabolic abnormalities, but conclusive evidence for pregnancies affected by maternal type 1 diabetes is still lacking. This study aims to analyse the relationship between maternal type 1 diabetes and the offspring’s metabolic health and investigate whether birthweight and/or changes in the offspring’s metabolome are in the potential pathway. We analysed data from 610 and 2169 offspring having a first-degree relative with type 1 diabetes from the TEENDIAB and BABYDIAB/BABYDIET cohorts, respectively. Anthropometric and metabolic outcomes, assessed longitudinally at 0.3–18 years of age, were compared between offspring of mothers with type 1 diabetes and offspring of non-diabetic mothers but with fathers or siblings with type 1 diabetes using mixed regression models. Non-targeted metabolomic measurements were carried out in 500 individuals from TEENDIAB and analysed with maternal type 1 diabetes and offspring overweight status. The offspring of mothers with type 1 diabetes had a higher BMI SD score (SDS) and an increased risk for being overweight than the offspring of non-diabetic mothers (e.g. OR for overweight status in TEENDIAB 2.40 [95% CI 1.41, 4.06]). Further, waist circumference SDS, fasting levels of glucose, insulin and C-peptide, and insulin resistance and abdominal obesity were significantly increased in the offspring of mothers with type 1 diabetes, even when adjusted for potential confounders and birthweight. Metabolite patterns related to androgenic steroids and branched-chain amino acids were found to be associated with offspring’s overweight status, but no significant associations were observed between maternal type 1 diabetes and metabolite concentrations in the offspring. Maternal type 1 diabetes is associated with offspring’s overweight status and metabolic health in later life, but this is unlikely to be caused by alterations in the offspring’s metabolome.

International audience; Background and aims: Preclinical experiments on animal models are essential to understand the mechanisms of cardiovascular disease (CVD). Metabolomics allows access to the metabolic perturbations associated with CVD in heart and vessels. Here we assessed which potential animal CVD model most closely mimics the serum metabolite signature of increased carotid intima-media thickness (cIMT) in humans, a clinical parameter widely accepted as a surrogate of CVD. Methods: A targeted mass spectrometry assay was used to quantify and compare a series of blood me-tabolites between 1362 individuals (KORA F4 cohort) and 5 animal CVD models: ApoE À/À , Ldlr À/À , and klotho-hypomorphic mice (kl/kl) and SHRSP rats with or without salt feeding. The metabolite signatures were obtained using linear regressions adjusted for various co-variates. Results: In human, increased cIMT [quartile Q4 vs. Q1] was associated with 26 metabolites (9 acylcar-nitines, 2 lysophosphatidylcholines, 9 phosphatidylcholines and 6 sphingomyelins). Acylcarnitines correlated preferentially with serum glucose and creatinine. Phospholipids correlated preferentially with cholesterol (total and LDL). The human signature correlated positively and significantly with Ldlr À/À and ApoE À/À mice, while correlation with kl/kl mice and SHRP rats was either negative and non-significant. Human and Ldlr À/À mice shared 11 significant metabolites displaying the same direction of regulation: 5 phosphatidylcholines, 1 lysophosphatidylcholines, 5 sphingomyelins; ApoE À/À mice shared 10. Conclusions: The human cIMT signature was partially mimicked by Ldlr À/À and ApoE À/À mice. These animal models might help better understand the biochemical and molecular mechanisms involved in the vessel metabolic perturbations associated with, and contributing to metabolic disorders in CVD.

An “omics” investigation of chronic widespread pain in 2 population samples reveals a genetic variant in CYP3A5 as underlying steroid hormone abnormalities seen in CWP., Chronic widespread musculoskeletal pain (CWP) is common, having a population prevalence of 10%. This study aimed to define the biological basis of the CWP/body mass association by using a systems biology approach. Adult female twins (n = 2444) from the TwinsUK registry who had extensive clinical, anthropometric, and “omic” data were included. Nontargeted metabolomics screening including 324 metabolites was carried out for CWP and body composition using dual-energy X-ray absorptiometry. The biological basis of these associations was explored through a genome-wide association study and replicated in an independent population sample (Cooperative Health Research in the Region of Augsburg [KORA] study, n = 2483). A causal role for the genetic variants identified was sought in CWP using a Mendelian randomisation study design. Fat mass/height2 was the body composition variable most strongly associated with CWP (TwinsUK: P = 2.4 × 10−15 and KORA: P = 1.59 × 10−10). Of 324 metabolites examined, epiandrosterone sulfate (EAS) was highly associated with both CWP (P = 1.05 × 10−09 in TwinsUK and P = 3.70 × 10−06 in KORA) and fat mass/height2. Genome-wide association study of EAS identified imputed single nucleotide polymorphism rs1581492 at 7q22.1 to be strikingly associated with EAS levels (P ≤ 2.49 × 10−78), and this result was replicated in KORA (P = 2.12 × 10−9). Mendelian randomization by rs1581492 genotype showed that EAS is unlikely to be causally related to CWP. Using an agnostic omics approach to focus on the association of CWP with body mass index, we have confirmed a steroid hormone association and identified a genetic variant upstream of the CYP genes, which likely controls this response. This study suggests that steroid hormone abnormalities result from pain rather than causing it, and EAS may provide a biomarker that identifies subgroups at risk of CWP.

Summary Age is the strongest risk factor for many diseases including neurodegenerative disorders, coronary heart disease, type 2 diabetes and cancer. Due to increasing life expectancy and low birth rates, the incidence of age‐related diseases is increasing in industrialized countries. Therefore, understanding the relationship between diseases and aging and facilitating healthy aging are major goals in medical research. In the last decades, the dimension of biological data has drastically increased with high‐throughput technologies now measuring thousands of (epi) genetic, expression and metabolic variables. The most common and so far successful approach to the analysis of these data is the so‐called reductionist approach. It consists of separately testing each variable for association with the phenotype of interest such as age or age‐related disease. However, a large portion of the observed phenotypic variance remains unexplained and a comprehensive understanding of most complex phenotypes is lacking. Systems biology aims to integrate data from different experiments to gain an understanding of the system as a whole rather than focusing on individual factors. It thus allows deeper insights into the mechanisms of complex traits, which are caused by the joint influence of several, interacting changes in the biological system. In this review, we look at the current progress of applying omics technologies to identify biomarkers of aging. We then survey existing systems biology approaches that allow for an integration of different types of data and highlight the need for further developments in this area to improve epidemiologic investigations.

Genome-wide association studies with metabolomics (mGWAS) identify genetically influenced metabotypes (GIMs), their ensemble defining the heritable part of every human's metabolic individuality. Knowledge of genetic variation in metabolism has many applications of biomedical and pharmaceutical interest, including the functional understanding of genetic associations with clinical endpoints, design of strategies to correct dysregulations in metabolic disorders, and the identification of genetic effect modifiers of metabolic disease biomarkers. Furthermore, it has been shown that GIMs provide testable hypotheses for functional genomics and metabolomics and for the identification of novel gene functions and metabolite identities. mGWAS with growing sample sizes and increasingly complex metabolic trait panels are being conducted, allowing for more comprehensive and systems based downstream analyses. The generated large data sets of genetic associations can now be mined by the biomedical research community and provide valuable resources for hypothesis driven studies. In this review, we provide a brief summary of the key aspects of mGWAS, followed by an update of recently published mGWAS. We then discuss new approaches of integrating and exploring mGWAS results and finish by presenting selected applications of GIMs in recent studies.

Alzheimer’s disease (AD) is the most common neurodegenerative disease presenting major health and economic challenges that continue to grow. Mechanisms of disease are poorly understood but significant data point to metabolic defects that might contribute to disease pathogenesis. The Alzheimer Disease Metabolomics Consortium (ADMC) in partnership with Alzheimer Disease Neuroimaging Initiative (ADNI) is creating a comprehensive biochemical database for AD. Using targeted and non- targeted metabolomics and lipidomics platforms we are mapping metabolic pathway and network failures across the trajectory of disease. In this report we present quantitative metabolomics data generated on serum from 199 control, 356 mild cognitive impairment and 175 AD subjects enrolled in ADNI1 using AbsoluteIDQ-p180 platform, along with the pipeline for data preprocessing and medication classification for confound correction. The dataset presented here is the first of eight metabolomics datasets being generated for broad biochemical investigation of the AD metabolome. We expect that these collective metabolomics datasets will provide valuable resources for researchers to identify novel molecular mechanisms contributing to AD pathogenesis and disease phenotypes.

IntroductionThe Alzheimer's Disease Research Summits of 2012 and 2015 incorporated experts from academia, industry, and nonprofit organizations to develop new research directions to transform our understanding of Alzheimer's disease (AD) and propel the development of critically needed therapies. In response to their recommendations, big data at multiple levels are being generated and integrated to study network failures in disease. We used metabolomics as a global biochemical approach to identify peripheral metabolic changes in AD patients and correlate them to cerebrospinal fluid pathology markers, imaging features, and cognitive performance.MethodsFasting serum samples from the Alzheimer's Disease Neuroimaging Initiative (199 control, 356 mild cognitive impairment, and 175 AD participants) were analyzed using the AbsoluteIDQ-p180 kit. Performance was validated in blinded replicates, and values were medication adjusted.Results Multivariable-adjusted analyses showed that sphingomyelins and ether-containing phosphatidylcholines were altered in preclinical biomarker-defined AD stages, whereas acylcarnitines and several amines, including the branched-chain amino acid valine and α-aminoadipic acid, changed in symptomatic stages. Several of the analytes showed consistent associations in the Rotterdam, Erasmus Rucphen Family, and Indiana Memory and Aging Studies. Partial correlation networks constructed for Aβ1–42, tau, imaging, and cognitive changes provided initial biochemical insights for disease-related processes. Coexpression networks interconnected key metabolic effectors of disease.DiscussionMetabolomics identified key disease-related metabolic changes and disease-progression-related changes. Defining metabolic changes during AD disease trajectory and its relationship to clinical phenotypes provides a powerful roadmap for drug and biomarker discovery.

The identification of phenotype-driven network modules in complex, multifluid metabolomics data poses a considerable challenge for statistical analysis and result interpretation. This is the case for phenotypes with only few associations ('sparse' effects), but, in particular, for phenotypes with a large number of metabolite associations ('dense' effects). Herein, we postulate that examining the data at different layers of resolution, from metabolites to pathways, will facilitate the interpretation of modules for both the sparse and the dense cases. We propose an approach for the phenotype-driven identification of modules on multifluid networks based on untargeted metabolomics data of plasma, urine, and saliva samples from the German Study of Health in Pomerania (SHIP-TREND) study. We generated a hierarchical, multifluid map of metabolism covering both metabolite and pathway associations using Gaussian graphical models. First, this map facilitates a fundamental understanding of metabolism within and across fluids for our study, and can serve as a valuable and downloadable resource. Second, based on this map, we then present an algorithm to identify regulated modules that associate with factors such as gender and insulin-like growth factor I (IGF-I) as examples of traits with dense and sparse associations, respectively. We found IGF-I to associate at the rather fine-grained metabolite level, while gender shows well-interpretable associations at pathway level. Our results confirm that a holistic and interpretable view of metabolic changes associated with a phenotype can only be obtained if different layers of metabolic resolution from multiple body fluids are considered., Metabolomics: phenotype-driven modules in a multifluid metabolic map Metabolism consists of complex interactions across various organs and body fluids, which poses a substantial challenge for the analysis of metabolic data. To address this problem, Jan Krumsiek from Helmholtz Zentrum München and colleagues used metabolomics measurements of plasma, urine, and saliva from 1000 people to statistically reconstruct a map of interactions in human metabolism. Based on this map, a novel approach that identifies highly correlated biochemical modules that are associated with a given phenotype, was tested for gender and insulin-like growth factor I (IGF-I). The identified modules provided insights into the interaction between metabolome and phenotype that reach beyond what can be found by commonly used statistical approaches for metabolomics. The approach is generic and can be readily applied to new datasets by other colleagues from the field.

We appreciate the comments of Jegatheesan et al. (1) on our article (2), as well as the remarks of Irving et al. (3). Jegatheesan et al. (1) also summarize citrulline’s beneficial effects, specifically in nonalcoholic fatty liver disease, and characterize this amino acid as a coin with two sides. They suggest citrulline as a marker of insulin resistance as well as an insulin-sensitizing amino acid (1). Although this is an interesting question, it was not the subject of our study and therefore could not be addressed (2). We would like to emphasize that we did not speculate about a direct association between high citrulline and insulin resistance. In fact, …

Genome-wide association scans with high-throughput metabolic profiling provide unprecedented insights into how genetic variation influences metabolism and complex disease. Here we report the most comprehensive exploration of genetic loci influencing human metabolism thus far, comprising 7,824 adult individuals from 2 European population studies. We report genome-wide significant associations at 145 metabolic loci and their biochemical connectivity with more than 400 metabolites in human blood. We extensively characterize the resulting in vivo blueprint of metabolism in human blood by integrating it with information on gene expression, heritability and overlap with known loci for complex disorders, inborn errors of metabolism and pharmacological targets. We further developed a database and web-based resources for data mining and results visualization. Our findings provide new insights into the role of inherited variation in blood metabolic diversity and identify potential new opportunities for drug development and for understanding disease.

BackgroundGenome-wide association studies (GWAS) have identified hundreds of loci influencing complex human traits, however, their biological mechanism of action remains mostly unknown. Recent accumulation of functional genomics (‘omics’) including metabolomics data opens up opportunities to provide a new insight into the functional role of specific changes in the genome. Functional genomic data are characterized by high dimensionality, presence of (strong) statistical dependencies between traits, and, potentially, complex genetic control. Therefore, analysis of such data asks for development of specific statistical genetic methods.ResultsWe propose a network-based, conditional approach to evaluate the impact of genetic variants on omics phenotypes (conditional GWAS, cGWAS). For each trait of interest, based on biological network, we select a set of other traits to be used as covariates in GWAS. The network could be reconstructed either from biological pathway databases or directly from the data. We evaluated our approach using data from a population-based KORA study (n=1,784, 1.7 M SNPs) with measured metabolomics data (151 metabolites) and demonstrated that our approach allows for identification of up to five additional loci not detected by conventional GWAS. We show that this gain in power is achieved through increased precision of genetic effect estimates, and in presence of specific ‘contra-intuitive’ pleiotropic scenarios (when genetic and environmental sources of covariance are acting in opposite manner). We justify existence of such scenarios, and discuss possible applications of our method beyond metabolomics.ConclusionsWe demonstrate that in context of metabolomics network-based, conditional genome-wide association analysis is able to dramatically increase power of identification of loci with specific ‘contra-intuitive’ pleiotropic architecture. Our method has modest computational costs, can utilize summary level GWAS data, and is applicable to other omics data types. We anticipate that application of our method to new and existing data sets will facilitate progress in understanding genetic bases of control of molecular and complex phenotypes.Short abstractWe propose a network-based, conditional approach for genome-wide analysis of multivariate omics phenotypes. Our methods can incorporate prior biological knowledge about biological pathways from external sources. We evaluated our approach using metabolomics data and demonstrated that our approach has bigger power and allows for identification of additional loci. We show that gain in power is achieved through increased precision of genetic effect estimates, and in presence of specific ‘contra-intuitive’ pleiotropic scenarios (when genetic and environmental sources of covariance are acting in opposite manner). We justify existence of such scenarios, and discuss possible applications of our method beyond metabolomics.

Using a nontargeted metabolomics approach of 447 fasting plasma metabolites, we searched for novel molecular markers that arise before and after hyperglycemia in a large population-based cohort of 2,204 females (115 type 2 diabetic [T2D] case subjects, 192 individuals with impaired fasting glucose [IFG], and 1,897 control subjects) from TwinsUK. Forty-two metabolites from three major fuel sources (carbohydrates, lipids, and proteins) were found to significantly correlate with T2D after adjusting for multiple testing; of these, 22 were previously reported as associated with T2D or insulin resistance. Fourteen metabolites were found to be associated with IFG. Among the metabolites identified, the branched-chain keto-acid metabolite 3-methyl-2-oxovalerate was the strongest predictive biomarker for IFG after glucose (odds ratio [OR] 1.65 [95% CI 1.39–1.95], P = 8.46 × 10−9) and was moderately heritable (h2 = 0.20). The association was replicated in an independent population (n = 720, OR 1.68 [ 1.34–2.11], P = 6.52 × 10−6) and validated in 189 twins with urine metabolomics taken at the same time as plasma (OR 1.87 [1.27–2.75], P = 1 × 10−3). Results confirm an important role for catabolism of branched-chain amino acids in T2D and IFG. In conclusion, this T2D-IFG biomarker study has surveyed the broadest panel of nontargeted metabolites to date, revealing both novel and known associated metabolites and providing potential novel targets for clinical prediction and a deeper understanding of causal mechanisms.

Summary Background Individuals with negative affectivity who are inhibited in social situations are characterized as distressed, or Type D, and have an increased risk of cardiovascular disease (CVD). The underlying biomechanisms that link this psychological affect to a pathological state are not well understood. This study applied a metabolomic approach to explore biochemical pathways that may contribute to the Type D personality. Methods Type D personality was determined by the Type D Scale-14. Small molecule biochemicals were measured using two complementary mass-spectrometry based metabolomics platforms. Metabolic profiles of Type D and non-Type D participants within a population-based study in Southern Germany were compared in cross-sectional regression analyses. The PHQ-9 and GAD-7 instruments were also used to assess symptoms of depression and anxiety, respectively, within this metabolomic study. Results 668 metabolites were identified in the serum of 1502 participants (age 32–77); 386 of these individuals were classified as Type D. While demographic and biomedical characteristics were equally distributed between the groups, a higher level of depression and anxiety was observed in Type D individuals. Significantly lower levels of the tryptophan metabolite kynurenine were associated with Type D (p-value corrected for multiple testing = 0.042), while no significant associations could be found for depression and anxiety. A Gaussian graphical model analysis enabled the identification of four potentially interesting metabolite networks that are enriched in metabolites (androsterone sulfate, tyrosine, indoxyl sulfate or caffeine) that associate nominally with Type D personality. Conclusions This study identified novel biochemical pathways associated with Type D personality and demonstrates that the application of metabolomic approaches in population studies can reveal mechanisms that may contribute to psychological health and disease.

Although association studies have unveiled numerous correlations of biochemical markers with age and age-related diseases, we still lack an understanding of their mutual dependencies. To find molecular pathways that underlie age-related diseases as well as their comorbidities, we integrated aging markers from four different high-throughput omics datasets, namely epigenomics, transcriptomics, glycomics and metabolomics, with a comprehensive set of disease phenotypes from 510 participants of the TwinsUK cohort. We used graphical random forests to assess conditional dependencies between omics markers and phenotypes while eliminating mediated associations. Applying this novel approach for multi-omics data integration yields a model consisting of seven modules that represent distinct aspects of aging. These modules are connected by hubs that potentially trigger comorbidities of age-related diseases. As an example, we identified urate as one of these key players mediating the comorbidity of renal disease with body composition and obesity. Body composition variables are in turn associated with inflammatory IgG markers, mediated by the expression of the hormone oxytocin. Thus, oxytocin potentially contributes to the development of chronic low-grade inflammation, which often accompanies obesity. Our multi-omics graphical model demonstrates the interconnectivity of age-related diseases and highlights molecular markers of the aging process that might drive disease comorbidities.

BACKGROUND/OBJECTIVES: Fatty liver disease (FLD) is an important intermediate trait along the cardiometabolic disease spectrum and strongly associates with type 2 diabetes. Knowledge of biological pathways implicated in FLD is limited. An untargeted metabolomic approach might unravel novel pathways related to FLD. SUBJECTS/METHODS: In a population-based sample (n=555) from Northern Germany, liver fat content was quantified as liver signal intensity using magnetic resonance imaging. Serum metabolites were determined using a non-targeted approach. Partial least squares regression was applied to derive a metabolomic score, explaining variation in serum metabolites and liver signal intensity. Associations of the metabolomic score with liver signal intensity and FLD were investigated in multivariable-adjusted robust linear and logistic regression models, respectively. Metabolites with a variable importance in the projection >1 were entered in in silico overrepresentation and pathway analyses. RESULTS: In univariate analysis, the metabolomics score explained 23.9% variation in liver signal intensity. A 1-unit increment in the metabolomic score was positively associated with FLD (n=219; odds ratio: 1.36; 95% confidence interval: 1.27-1.45) adjusting for age, sex, education, smoking and physical activity. A simplified score based on the 15 metabolites with highest variable importance in the projection statistic showed similar associations. Overrepresentation and pathway analyses highlighted branched-chain amino acids and derived gamma-glutamyl dipeptides as significant correlates of FLD. CONCLUSIONS: A serum metabolomic profile was associated with FLD and liver fat content. We identified a simplified metabolomics score, which should be evaluated in prospective studies.European Journal of Clinical Nutrition advance online publication, 5 April 2017; doi:10.1038/ejcn.2017.43.

Identification of novel biomarkers for type 2 diabetes and their genetic determinants could lead to improved understanding of causal pathways and improve risk prediction.In this study, we used data from non-targeted metabolomics performed using liquid chromatography coupled with tandem mass spectrometry in three Swedish cohorts (Uppsala Longitudinal Study of Adult Men [ULSAM], n = 1138; Prospective Investigation of the Vasculature in Uppsala Seniors [PIVUS], n = 970; TwinGene, n = 1630). Metabolites associated with impaired fasting glucose (IFG) and/or prevalent type 2 diabetes were assessed for associations with incident type 2 diabetes in the three cohorts followed by replication attempts in the Cooperative Health Research in the Region of Augsburg (KORA) S4 cohort (n = 855). Assessment of the association of metabolite-regulating genetic variants with type 2 diabetes was done using data from a meta-analysis of genome-wide association studies.Out of 5961 investigated metabolic features, 1120 were associated with prevalent type 2 diabetes and IFG and 70 were annotated to metabolites and replicated in the three cohorts. Fifteen metabolites were associated with incident type 2 diabetes in the four cohorts combined (358 events) following adjustment for age, sex, BMI, waist circumference and fasting glucose. Novel findings included associations of higher values of the bile acid deoxycholic acid and monoacylglyceride 18:2 and lower concentrations of cortisol with type 2 diabetes risk. However, adding metabolites to an existing risk score improved model fit only marginally. A genetic variant within the CYP7A1 locus, encoding the rate-limiting enzyme in bile acid synthesis, was found to be associated with lower concentrations of deoxycholic acid, higher concentrations of LDL-cholesterol and lower type 2 diabetes risk. Variants in or near SGPP1, GCKR and FADS1/2 were associated with diabetes-associated phospholipids and type 2 diabetes.We found evidence that the metabolism of bile acids and phospholipids shares some common genetic origin with type 2 diabetes.Metabolomics data have been deposited in the Metabolights database, with accession numbers MTBLS93 (TwinGene), MTBLS124 (ULSAM) and MTBLS90 (PIVUS).

Angiotensin-I-converting enzyme (ACE) inhibitors are an important class of antihypertensives whose action on the human organism is still not fully understood. Although it is known that ACE especially cleaves COOH-terminal dipeptides from active polypeptides, the whole range of substrates and products is still unknown. When analyzing the action of ACE inhibitors, effects of genetic variation on metabolism need to be considered since genetic variance in the ACE gene locus was found to be associated with ACE-concentration in blood as well as with changes in the metabolic profiles of a general population. To investigate the interactions between genetic variance at the ACE-locus and the influence of ACE-therapy on the metabolic status we analyzed 517 metabolites in 1,361 participants from the KORA F4 study. We replicated our results in 1,964 individuals from TwinsUK. We observed differences in the concentration of five dipeptides and three ratios of di- and oligopeptides between ACE inhibitor users and non-users that were genotype dependent. Such changes in the concentration affected major homozygotes, and to a lesser extent heterozygotes, while minor homozygotes showed no or only small changes in the metabolite status. Two of these resulting dipeptides, namely aspartylphenylalanine and phenylalanylserine, showed significant associations with blood pressure which qualifies them - and perhaps also the other dipeptides - as readouts of ACE-activity. Since so far ACE activity measurement is substrate specific due to the usage of only one oligopeptide, taking several dipeptides as potential products of ACE into account may provide a broader picture of the ACE activity.