A new statistical workflow (R-packages based) to investigate associations between one variable of interest and the metabolome.

INTRODUCTION: In metabolomics, the investigation of associations between the metabolome and one trait of interest is a key research question. However, statistical analyses of such associations are often challenging. Statistical tools enabling resilient verification and clear presentation are therefore highly desired. OBJECTIVES: Our aim is to provide a contribution for statistical analysis of metabolomics data, offering a widely applicable open-source statistical workflow, which considers the intrinsic complexity of metabolomics data. METHODS: We combined selected R packages tailored for all properties of heterogeneous metabolomics datasets, where metabolite parameters typically (i) are analyzed in different matrices, (ii) are measured on different analytical platforms with different precision, (iii) are analyzed by targeted as well as non-targeted methods, (iv) are scaled variously, (v) reveal heterogeneous variances, (vi) may be correlated, (vii) may have only few values or values below a detection limit, or (viii) may be incomplete. RESULTS: The code is shared entirely and freely available. The workflow output is a table of metabolites associated with a trait of interest and a compact plot for high-quality results visualization. The workflow output and its utility are presented by applying it to two previously published datasets: one dataset from our own lab and another dataset taken from the repository MetaboLights. CONCLUSION: Robustness and benefits of the statistical workflow were clearly demonstrated, and everyone can directly re-use it for analysis of own data.

Transferring entropy to the realm of GxG interactions.

Genome-wide association studies are moving to genome-wide interaction studies, as the genetic background of many diseases appears to be more complex than previously supposed. Thus, many statistical approaches have been proposed to detect gene-gene (GxG) interactions, among them numerous information theory-based methods, inspired by the concept of entropy. These are suggested as particularly powerful and, because of their nonlinearity, as better able to capture nonlinear relationships between genetic variants and/or variables. However, the introduced entropy-based estimators differ to a surprising extent in their construction and even with respect to the basic definition of interactions. Also, not every entropy-based measure for interaction is accompanied by a proper statistical test. To shed light on this, a systematic review of the literature is presented answering the following questions: (1) How are GxG interactions defined within the framework of information theory? (2) Which entropy-based test statistics are available? (3) Which underlying distribution do the test statistics follow? (4) What are the given strengths and limitations of these test statistics?

Coding Variation in ANGPTL4, LPL, and SVEP1 and the Risk of Coronary Disease.

BACKGROUND: The discovery of low-frequency coding variants affecting the risk of coronary artery disease has facilitated the identification of therapeutic targets. METHODS: Through DNA genotyping, we tested 54,003 coding-sequence variants covering 13,715 human genes in up to 72,868 patients with coronary artery disease and 120,770 controls who did not have coronary artery disease. Through DNA sequencing, we studied the effects of loss-of-function mutations in selected genes. RESULTS: We confirmed previously observed significant associations between coronary artery disease and low-frequency missense variants in the genes LPA and PCSK9. We also found significant associations between coronary artery disease and low-frequency missense variants in the genes SVEP1 (p.D2702G; minor-allele frequency, 3.60%; odds ratio for disease, 1.14; P=4.2×10(-10)) and ANGPTL4 (p.E40K; minor-allele frequency, 2.01%; odds ratio, 0.86; P=4.0×10(-8)), which encodes angiopoietin-like 4. Through sequencing of ANGPTL4, we identified 9 carriers of loss-of-function mutations among 6924 patients with myocardial infarction, as compared with 19 carriers among 6834 controls (odds ratio, 0.47; P=0.04); carriers of ANGPTL4 loss-of-function alleles had triglyceride levels that were 35% lower than the levels among persons who did not carry a loss-of-function allele (P=0.003). ANGPTL4 inhibits lipoprotein lipase; we therefore searched for mutations in LPL and identified a loss-of-function variant that was associated with an increased risk of coronary artery disease (p.D36N; minor-allele frequency, 1.9%; odds ratio, 1.13; P=2.0×10(-4)) and a gain-of-function variant that was associated with protection from coronary artery disease (p.S447*; minor-allele frequency, 9.9%; odds ratio, 0.94; P=2.5×10(-7)). CONCLUSIONS: We found that carriers of loss-of-function mutations in ANGPTL4 had triglyceride levels that were lower than those among noncarriers; these mutations were also associated with protection from coronary artery disease. (Funded by the National Institutes of Health and others.).

Genetic variants associated with celiac disease and the risk for coronary artery disease.

Epidemiological evidence suggests that patients with celiac disease are at increased risk for coronary artery disease (CAD). Genetic-epidemiological analyses identified many single nucleotide polymorphisms (SNPs) associated with celiac disease. If there is a causal relation between celiac disease and CAD, one might expect that risk alleles primarily associated with celiac disease also increase the risk of CAD. In this study we identified from literature 41 SNPs that have been previously described to be genome-wide associated with celiac disease (p < 5 × 10(-08)). These SNPs were evaluated for their association with CAD in the Coronary ARtery DIsease Genome-wide Replication and Meta-analysis (CARDIoGRAM) dataset, a meta-analysis comprising genome-wide SNP association data from 22,233 CAD cases and 64,762 controls. 24 out of 41 (58.5 %) risk alleles for celiac disease displayed a positive association with CAD (CAD-OR range 1.001-1.081). The remaining risk alleles for celiac disease (n = 16) revealed CAD-ORs of ≤1.0 (range 0.951-1.0). The proportion of CAD associated alleles was greater but did not differ significantly from the proportion of 50 % expected by chance (p = 0.069). One SNP (rs653178 at the SH2B3/ATXN2 locus) displayed study-wise statistically significant association with CAD with directionality consistent effects on celiac disease and CAD. However, the effect of this locus is most likely driven by pleiotropic effects on multiple other diseases. In conclusion, this genetically based approach provided no convincing evidence that SNPs associated with celiac disease contribute to the risk of CAD. Hence, common non-genetic factors may play a more important role explaining the coincidence of these two complex disease conditions.

Gut Microbiome Analysis for Personalized Nutrition: The State of Science.

Whereas most concepts of personalized nutrition (PN) in the past, included genotyping, recent years have brought new approaches that include microbiome analysis to optimize recommendations for diet and lifestyle changes. The new approach, offered by companies, that microbiome analysis provides a real benefit to either more concise recommendations or for increased compliance to PN, is largely lacking scientific validation. Although the microbiome field shows enormous proliferation, it has some major flaws that make its use in the public health domain currently critical. Starting with the quality and representative character of the stool samples, its processing and analysis as well as assembly of metagenome data and the interpretation. Moreover, there is still no consensus of what constitutes a "normal/healthy" microbiome, nor what features characterize a dysbiotic microbiome. And, based on hundreds of individual parameters and environmental factors, the intestinal microbiome shows a huge variability and consequently changing one factor-such as food intake-is likely to have a limited impact in achieving optimized health. The present review intends to summarize the state of consolidated knowledge on human gut microbiome in the context of diet and disease, its key features, and its influencing factors as well as its "add-on" quality for PN offers.

Acute effects of moderate vs. vigorous endurance exercise on urinary metabolites in healthy, young, physically active men-A multi-platform metabolomics approach.

Introduction: Endurance exercise alters whole-body as well as skeletal muscle metabolism and physiology, leading to improvements in performance and health. However, biological mechanisms underlying the body's adaptations to different endurance exercise protocols are not entirely understood. Methods: We applied a multi-platform metabolomics approach to identify urinary metabolites and associated metabolic pathways that distinguish the acute metabolic response to two endurance exercise interventions at distinct intensities. In our randomized crossover study, 16 healthy, young, and physically active men performed 30 min of continuous moderate exercise (CME) and continuous vigorous exercise (CVE). Urine was collected during three post-exercise sampling phases (U01/U02/U03: until 45/105/195 min post-exercise), providing detailed temporal information on the response of the urinary metabolome to CME and CVE. Also, fasting spot urine samples were collected pre-exercise (U00) and on the following day (U04). While untargeted two-dimensional gas chromatography-mass spectrometry (GC×GC-MS) led to the detection of 608 spectral features, 44 metabolites were identified and quantified by targeted nuclear magnetic resonance (NMR) spectroscopy or liquid chromatography-mass spectrometry (LC-MS). Results: 104 urinary metabolites showed at least one significant difference for selected comparisons of sampling time points within or between exercise trials as well as a relevant median fold change >1.5 or <0. 6 ¯ (NMR, LC-MS) or >2.0 or <0.5 (GC×GC-MS), being classified as either exercise-responsive or intensity-dependent. Our findings indicate that CVE induced more profound alterations in the urinary metabolome than CME, especially at U01, returning to baseline within 24 h after U00. Most differences between exercise trials are likely to reflect higher energy requirements during CVE, as demonstrated by greater shifts in metabolites related to glycolysis (e.g., lactate, pyruvate), tricarboxylic acid cycle (e.g., cis-aconitate, malate), purine nucleotide breakdown (e.g., hypoxanthine), and amino acid mobilization (e.g., alanine) or degradation (e.g., 4-hydroxyphenylacetate). Discussion: To conclude, this study provided first evidence of specific urinary metabolites as potential metabolic markers of endurance exercise intensity. Future studies are needed to validate our results and to examine whether acute metabolite changes in urine might also be partly reflective of mechanisms underlying the health- or performance-enhancing effects of endurance exercise, particularly if performed at high intensities.

What is the promise of personalised nutrition?

Personalised nutrition (PN) is an emerging field that bears great promise. Several definitions of PN have been proposed and different modelling approaches have been used to claim PN effects. We tentatively propose to group these approaches into two categories, which we term outcome-based and population reference approaches, respectively. Understanding the fundamental differences between these two types of modelling approaches may allow a more realistic appreciation of what to expect from PN interventions presently and may be helpful for designing and planning future studies investigating PN interventions.

Exploring the Diversity of Sugar Compounds in Healthy, Prediabetic, and Diabetic Volunteers.

SCOPE: Diabetes is thought to primarily represent a disturbance of carbohydrate metabolism; however, population studies employing metabolomics have mainly identified plasma amino acids and lipids, or their products, as biomarkers. In this pilot study, the aim is to analyze a wide spectrum of sugar compounds in the fasting state and during an oral glucose tolerance test (OGTT) in healthy, prediabetic, and type 2 diabetic volunteers. METHODS AND RESULTS: The three volunteer groups underwent a standard OGTT. Plasma samples obtained in the fasting state, 30 and 90 min after the OGTT, are subjected to a semitargeted GC-MS (gas chromatography-mass spectrometry) sugar profiling. Overall, 40 sugars are detected in plasma, of which some are yet unknown to change during an OGTT. Several sugars (e.g., trehalose) reveal significant differences between the volunteer groups both in fasting plasma and in distinct time courses after the OGTT. This suggests an endogenous production from orally absorbed glucose and/or an insulin-dependent production/removal from plasma. CONCLUSION: It is demonstrated that more sugars than expected can be found in human plasma. Since some of these show characteristic differences depending on health status, it may be worthwhile to assess their usability as biomarkers for diagnosing early-stage insulin resistance and type 2 diabetes.

Nutrimetabolomics: An Integrative Action for Metabolomic Analyses in Human Nutritional Studies.

The life sciences are currently being transformed by an unprecedented wave of developments in molecular analysis, which include important advances in instrumental analysis as well as biocomputing. In light of the central role played by metabolism in nutrition, metabolomics is rapidly being established as a key analytical tool in human nutritional studies. Consequently, an increasing number of nutritionists integrate metabolomics into their study designs. Within this dynamic landscape, the potential of nutritional metabolomics (nutrimetabolomics) to be translated into a science, which can impact on health policies, still needs to be realized. A key element to reach this goal is the ability of the research community to join, to collectively make the best use of the potential offered by nutritional metabolomics. This article, therefore, provides a methodological description of nutritional metabolomics that reflects on the state-of-the-art techniques used in the laboratories of the Food Biomarker Alliance (funded by the European Joint Programming Initiative "A Healthy Diet for a Healthy Life" (JPI HDHL)) as well as points of reflections to harmonize this field. It is not intended to be exhaustive but rather to present a pragmatic guidance on metabolomic methodologies, providing readers with useful "tips and tricks" along the analytical workflow.

The complex human urinary sugar profile: determinants revealed in the cross-sectional KarMeN study.

Background: Although sugars and sugar derivatives are an important class of metabolites involved in many physiologic processes, there is limited knowledge on their occurrence and pattern in biofluids. Objective: Our aim was to obtain a comprehensive urinary sugar profile of healthy participants and to demonstrate the wide applicability and usefulness of this sugar profiling approach for nutritional as well as clinical studies. Design: In the cross-sectional KarMeN study, the 24-h urine samples of 301 healthy participants on an unrestricted diet, assessed via a 24-h recall, were analyzed by a newly developed semitargeted gas chromatography-mass spectrometry (GC-MS) profiling method that enables the detection of known and unknown sugar compounds. Statistical analyses were performed with respect to associations of sex and diet with the urinary sugar profile. Results: In total, 40 known and 15 unknown sugar compounds were detected in human urine, ranging from mono- and disaccharides, polyols, and sugar acids to currently unknown sugar-like compounds. A number of rarely analyzed sugars were found in urine samples. Maltose was found in statistically higher concentrations in the urine of women compared with men and was also associated with menopausal status. Further, a number of individual sugar compounds associated with the consumption of specific foods, such as avocado, or food groups, such as alcoholic beverages and dairy products, were identified. Conclusions: We here provide data on the complex nature of the sugar profile in human urine, of which some compounds may have the potential to serve as dietary markers or early disease biomarkers. Thus, comprehensive urinary sugar profiling not only has the potential to increase our knowledge of host sugar metabolism, but can also reveal new dietary markers after consumption of individual food items, and may lead to the identification of early disease biomarkers in the future. The KarMeN study was registered at drks.de as DRKS00004890.

Associations of current diet with plasma and urine TMAO in the KarMeN study: direct and indirect contributions.

SCOPE: Knowledge on the influence of current diet on trimethylamine-N-oxide (TMAO) levels in humans is still inconsistent. Thus, we aimed to investigate associations of current diet with urine and plasma TMAO levels and to determine the effect of different foods on TMAO variation. METHODS AND RESULTS: TMAO concentrations of 297 healthy individuals were assessed using 1 H-NMR spectroscopy for 24 h urine collection and spot urine, and LC-MS for plasma. Of 35 assessed food groups, those with a correlation of ρ >|0.15| with plasma or urine TMAO levels were further investigated in multivariate linear regression models showing current fish and (red) meat consumption as plausible dietary sources of TMAO. Overall, explained variance of TMAO levels by current diet and co-variables (age, sex, lean body mass, glomerular filtration rate) was small. Associations with urine and plasma concentrations differed depending on the TMAO source. Fish consumption was associated with urine and plasma TMAO concentrations, whereas meat consumption was only associated with TMAO concentrations in plasma. Furthermore, associations of plasma TMAO concentration with fish consumption were two times stronger than with meat consumption. CONCLUSION: Meat and fish consumption differentially affects TMAO concentrations in body fluids. Only a small fraction of variance is explained by current diet.

Systematic Evaluation of Pleiotropy Identifies 6 Further Loci Associated With Coronary Artery Disease.

BACKGROUND: Genome-wide association studies have so far identified 56 loci associated with risk of coronary artery disease (CAD). Many CAD loci show pleiotropy; that is, they are also associated with other diseases or traits. OBJECTIVES: This study sought to systematically test if genetic variants identified for non-CAD diseases/traits also associate with CAD and to undertake a comprehensive analysis of the extent of pleiotropy of all CAD loci. METHODS: In discovery analyses involving 42,335 CAD cases and 78,240 control subjects we tested the association of 29,383 common (minor allele frequency >5%) single nucleotide polymorphisms available on the exome array, which included a substantial proportion of known or suspected single nucleotide polymorphisms associated with common diseases or traits as of 2011. Suggestive association signals were replicated in an additional 30,533 cases and 42,530 control subjects. To evaluate pleiotropy, we tested CAD loci for association with cardiovascular risk factors (lipid traits, blood pressure phenotypes, body mass index, diabetes, and smoking behavior), as well as with other diseases/traits through interrogation of currently available genome-wide association study catalogs. RESULTS: We identified 6 new loci associated with CAD at genome-wide significance: on 2q37 (KCNJ13-GIGYF2), 6p21 (C2), 11p15 (MRVI1-CTR9), 12q13 (LRP1), 12q24 (SCARB1), and 16q13 (CETP). Risk allele frequencies ranged from 0.15 to 0.86, and odds ratio per copy of the risk allele ranged from 1.04 to 1.09. Of 62 new and known CAD loci, 24 (38.7%) showed statistical association with a traditional cardiovascular risk factor, with some showing multiple associations, and 29 (47%) showed associations at p < 1 × 10-4 with a range of other diseases/traits. CONCLUSIONS: We identified 6 loci associated with CAD at genome-wide significance. Several CAD loci show substantial pleiotropy, which may help us understand the mechanisms by which these loci affect CAD risk.