metaP-Server: A Web-Based Metabolomics Data Analysis Tool

Metabolomics is an emerging field that is based on the quantitative measurement of as many small organic molecules occurring in a biological sample as possible. Due to recent technical advances, metabolomics can now be used widely as an analytical high-throughput technology in drug testing and epidemiological metabolome and genome wide association studies. Analogous to chip-based gene expression analyses, the enormous amount of data produced by modern kit-based metabolomics experiments poses new challenges regarding their biological interpretation in the context of various sample phenotypes. We developed metaP-server to facilitate data interpretation. metaP-server provides automated and standardized data analysis for quantitative metabolomics data, covering the following steps from data acquisition to biological interpretation: (i) data quality checks, (ii) estimation of reproducibility and batch effects, (iii) hypothesis tests for multiple categorical phenotypes, (iv) correlation tests for metric phenotypes, (v) optionally including all possible pairs of metabolite concentration ratios, (vi) principal component analysis (PCA), and (vii) mapping of metabolites onto colored KEGG pathway maps. Graphical output is clickable and cross-linked to sample and metabolite identifiers. Interactive coloring of PCA and bar plots by phenotype facilitates on-line data exploration. For users of commercial metabolomics kits, cross-references to the HMDB, LipidMaps, KEGG, PubChem, and CAS databases are provided. metaP-server is freely accessible at http://metabolomics.helmholtz-muenchen.de/metap2/

Metabolomics approach reveals effects of antihypertensives and lipid-lowering drugs on the human metabolism

The mechanism of antihypertensive and lipid-lowering drugs on the human organism is still not fully understood. New insights on the drugs' action can be provided by a metabolomics-driven approach, which offers a detailed view of the physiological state of an organism. Here, we report a metabolome-wide association study with 295 metabolites in human serum from 1,762 participants of the KORA F4 (Cooperative Health Research in the Region of Augsburg) study population. Our intent was to find variations of metabolite concentrations related to the intake of various drug classes and-based on the associations found-to generate new hypotheses about on-target as well as off-target effects of these drugs. In total, we found 41 significant associations for the drug classes investigated: For beta-blockers (11 associations), angiotensin-converting enzyme (ACE) inhibitors (four assoc.), diuretics (seven assoc.), statins (ten assoc.), and fibrates (nine assoc.) the top hits were pyroglutamine, phenylalanylphenylalanine, pseudouridine, 1-arachidonoylglycerophosphocholine, and 2-hydroxyisobutyrate, respectively. For beta-blockers we observed significant associations with metabolite concentrations that are indicative of drug side-effects, such as increased serotonin and decreased free fatty acid levels. Intake of ACE inhibitors and statins associated with metabolites that provide insight into the action of the drug itself on its target, such as an association of ACE inhibitors with des-Arg(9)-bradykinin and aspartylphenylalanine, a substrate and a product of the drug-inhibited ACE. The intake of statins which reduce blood cholesterol levels, resulted in changes in the concentration of metabolites of the biosynthesis as well as of the degradation of cholesterol. Fibrates showed the strongest association with 2-hydroxyisobutyrate which might be a breakdown product of fenofibrate and, thus, a possible marker for the degradation of this drug in the human organism. The analysis of diuretics showed a heterogeneous picture that is difficult to interpret. Taken together, our results provide a basis for a deeper functional understanding of the action and side-effects of antihypertensive and lipid-lowering drugs in the general population

Body Fat Free Mass Is Associated with the Serum Metabolite Profile in a Population-Based Study

To characterise the influence of the fat free mass on the metabolite profile in serum samples from participants of the population-based KORA (Cooperative Health Research in the Region of Augsburg) S4 study. Analyses were based on metabolite profile from 965 participants of the S4 and 890 weight-stable subjects of its seven-year follow-up study (KORA F4). 190 different serum metabolites were quantified in a targeted approach including amino acids, acylcarnitines, phosphatidylcholines (PCs), sphingomyelins and hexose. Associations between metabolite concentrations and the fat free mass index (FFMI) were analysed using adjusted linear regression models. To draw conclusions on enzymatic reactions, intra-metabolite class ratios were explored. Pairwise relationships among metabolites were investigated and illustrated by means of Gaussian graphical models (GGMs). We found 339 significant associations between FFMI and various metabolites in KORA S4. Among the most prominent associations (p-values 4.75 × 10(-16)-8.95 × 10(-06)) with higher FFMI were increasing concentrations of the branched chained amino acids (BCAAs), ratios of BCAAs to glucogenic amino acids, and carnitine concentrations. For various PCs, a decrease in chain length or in saturation of the fatty acid moieties could be observed with increasing FFMI, as well as an overall shift from acyl-alkyl PCs to diacyl PCs. These findings were reproduced in KORA F4. The established GGMs supported the regression results and provided a comprehensive picture of the relationships between metabolites. In a sub-analysis, most of the discovered associations did not exist in obese subjects in contrast to non-obese subjects, possibly indicating derangements in skeletal muscle metabolism. A set of serum metabolites strongly associated with FFMI was identified and a network explaining the relationships among metabolites was established. These results offer a novel and more complete picture of the FFMI effects on serum metabolites in a data-driven network

Genetics Meets Metabolomics: A Genome-Wide Association Study of Metabolite Profiles in Human Serum

The rapidly evolving field of metabolomics aims at a comprehensive measurement of ideally all endogenous metabolites in a cell or body fluid. It thereby provides a functional readout of the physiological state of the human body. Genetic variants that associate with changes in the homeostasis of key lipids, carbohydrates, or amino acids are not only expected to display much larger effect sizes due to their direct involvement in metabolite conversion modification, but should also provide access to the biochemical context of such variations, in particular when enzyme coding genes are concerned. To test this hypothesis, we conducted what is, to the best of our knowledge, the first GWA study with metabolomics based on the quantitative measurement of 363 metabolites in serum of 284 male participants of the KORA study. We found associations of frequent single nucleotide polymorphisms (SNPs) with considerable differences in the metabolic homeostasis of the human body, explaining up to 12% of the observed variance. Using ratios of certain metabolite concentrations as a proxy for enzymatic activity, up to 28% of the variance can be explained (p-values 10−16 to 10−21). We identified four genetic variants in genes coding for enzymes (FADS1, LIPC, SCAD, MCAD) where the corresponding metabolic phenotype (metabotype) clearly matches the biochemical pathways in which these enzymes are active. Our results suggest that common genetic polymorphisms induce major differentiations in the metabolic make-up of the human population. This may lead to a novel approach to personalized health care based on a combination of genotyping and metabolic characterization. These genetically determined metabotypes may subscribe the risk for a certain medical phenotype, the response to a given drug treatment, or the reaction to a nutritional intervention or environmental challenge

Metabolites of milk intake: a metabolomic approach in UK twins with findings replicated in two European cohorts

Purpose: Milk provides a significant source of calcium, protein, vitamins and other minerals to Western populations throughout life. Due to its widespread use, the metabolic and health impact of milk consumption warrants further investigation and biomarkers would aid epidemiological studies.  Methods: Milk intake assessed by a validated food frequency questionnaire was analyzed against fasting blood metabolomic profiles from two metabolomic platforms in females from the TwinsUK cohort (n = 3559). The top metabolites were then replicated in two independent populations (EGCUT, n = 1109 and KORA, n = 1593), and the results from all cohorts were meta-analyzed.  Results: Four metabolites were significantly associated with milk intake in the TwinsUK cohort after adjustment for multiple testing (P < 8.08 × 10−5) and covariates (BMI, age, batch effects, family relatedness and dietary covariates) and replicated in the independent cohorts. Among the metabolites identified, the carnitine metabolite trimethyl-N-aminovalerate (β = 0.012, SE = 0.002, P = 2.98 × 10−12) and the nucleotide uridine (β = 0.004, SE = 0.001, P = 9.86 × 10−6) were the strongest novel predictive biomarkers from the non-targeted platform. Notably, the association between trimethyl-N-aminovalerate and milk intake was significant in a group of MZ twins discordant for milk intake (β = 0.050, SE = 0.015, P = 7.53 × 10−4) and validated in the urine of 236 UK twins (β = 0.091, SE = 0.032, P = 0.004). Two metabolites from the targeted platform, hydroxysphingomyelin C14:1 (β = 0.034, SE = 0.005, P = 9.75 × 10−14) and diacylphosphatidylcholine C28:1 (β = 0.034, SE = 0.004, P = 4.53 × 10−16), were also replicated.  Conclusions: We identified and replicated in independent populations four novel biomarkers of milk intake: trimethyl-N-aminovalerate, uridine, hydroxysphingomyelin C14:1 and diacylphosphatidylcholine C28:1. Together, these metabolites have potential to objectively examine and refine milk-disease associations

Metabolic Profiling Reveals Distinct Variations Linked to Nicotine Consumption in Humans — First Results from the KORA Study

Exposure to nicotine during smoking causes a multitude of metabolic changes that are poorly understood. We quantified and analyzed 198 metabolites in 283 serum samples from the human cohort KORA (Cooperative Health Research in the Region of Augsburg). Multivariate analysis of metabolic profiles revealed that the group of smokers could be clearly differentiated from the groups of former smokers and non-smokers. Moreover, 23 lipid metabolites were identified as nicotine-dependent biomarkers. The levels of these biomarkers are all up-regulated in smokers compared to those in former and non-smokers, except for three acyl-alkyl-phosphatidylcholines (e.g. plasmalogens). Consistently significant results were further found for the ratios of plasmalogens to diacyl-phosphatidylcolines, which are reduced in smokers and regulated by the enzyme alkylglycerone phosphate synthase (alkyl-DHAP) in both ether lipid and glycerophospholipid pathways. Notably, our metabolite profiles are consistent with the strong down-regulation of the gene for alkyl-DHAP (AGPS) in smokers that has been found in a study analyzing gene expression in human lung tissues. Our data suggest that smoking is associated with plasmalogen-deficiency disorders, caused by reduced or lack of activity of the peroxisomal enzyme alkyl-DHAP. Our findings provide new insight into the pathophysiology of smoking addiction. Activation of the enzyme alkyl-DHAP by small molecules may provide novel routes for therapy

Hundreds of variants clustered in genomic loci and biological pathways affect human height

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

Variation in the SERPINA6SERPINA1 locusalters morning plasma cortisol, hepatic corticosteroid binding globulin expression, gene expressionin peripheral tissues, and risk of cardiovascular disease

The stress hormone cortisol modulates fuel metabolism, cardiovascular homoeostasis, mood, inflammation and cognition. The CORtisol NETwork (CORNET) consortium previously identified a single locus associated with morning plasma cortisol. Identifying additional genetic variants that explain more of the variance in cortisol could provide new insights into cortisol biology and provide statistical power to test the causative role of cortisol in common diseases. The CORNET consortium extended its genome-wide association meta-analysis for morning plasma cortisol from 12,597 to 25,314 subjects and from ~2.2 M to ~7 M SNPs, in 17 population-based cohorts of European ancestries. We confirmed the genetic association with SERPINA6/SERPINA1. This locus contains genes encoding corticosteroid binding globulin (CBG) and α1-antitrypsin. Expression quantitative trait loci (eQTL) analyses undertaken in the STARNET cohort of 600 individuals showed that specific genetic variants within the SERPINA6/SERPINA1 locus influence expression of SERPINA6 rather than SERPINA1 in the liver. Moreover, trans-eQTL analysis demonstrated effects on adipose tissue gene expression, suggesting that variation

Meta-analysis of exome array data identifies six novel genetic loci for lung function [version 1; peer review:1 approved, 1 approved with reservations]

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV1), forced vital capacity (FVC) and the ratio of FEV1 to FVC (FEV1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P&lt;2•8x10 -7 ) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs (SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.</p