MetaMapp: mapping and visualizing metabolomic data by integrating information from biochemical pathways and chemical and mass spectral similarity

The promise of metabolic phenotyping in gastroenterology and hepatology

Disease risk and treatment response are determined, at the individual level, by a complex history of genetic and environmental interactions, including those with our endogenous microbiomes. Personalized health care requires a deep understanding of patient biology that can now be measured using a range of '-omics' technologies. Patient stratification involves the identification of genetic and/or phenotypic disease subclasses that require different therapeutic strategies. Stratified medicine approaches to disease diagnosis, prognosis and therapeutic response monitoring herald a new dimension in patient care. Here, we explore the potential value of metabolic profiling as applied to unmet clinical needs in gastroenterology and hepatology. We describe potential applications in a number of diseases, with emphasis on large-scale population studies as well as metabolic profiling on the individual level, using spectrometric and imaging technologies that will leverage the discovery of mechanistic information and deliver novel health care solutions to improve clinical pathway management.

Changes in urinary metabolic profiles of colorectal cancer patients enrolled in a prospective cohort study (ColoCare)

INTRODUCTION

Metabolomics is a valuable tool for biomarker screening of colorectal cancer (CRC). In this study, we profiled the urinary metabolomes of patients enrolled in a prospective patient cohort (ColoCare). We aimed to describe changes in the metabolome in the longer clinical follow-up and describe initial predictors as candidate markers with possibly prognostic significance.

METHODS

In total, 199 urine samples from CRC patients pre-surgery (n=97), 1-8 days post-surgery (n=12) and then after 6 and 12 months (n=52 and 38, respectively) were analyzed using both GC-MS and ¹H-NMR. Both datasets were analyzed separately with built in uni- and multivariate analyses of Metaboanalyst 2.0. Furthermore, adjusted linear mixed effects regression models were constructed.

RESULTS

Many concentrations of the metabolites derived from the gut microbiome were affected by CRC surgery, presumably indicating a tumor-induced shift in bacterial species. Associations of the microbial metabolites with disease stage indicate an important role of the gut microbiome in CRC.We were able to differentiate the metabolite profiles of CRC patients prior to surgery from those at any post-surgery timepoint using a multivariate model containing 20 marker metabolites (AUCROC=0.89; 95% CI:0.84-0.95).

CONCLUSION

To the best of our knowledge, this is one of the first metabolomic studies to follow CRC patients in a prospective setting with repeated urine sampling over time. We were able to confirm markers initially identified in case-control studies and pin point metabolites which may serve as candidates for prognostic biomarkers of CRC.

Statistical-based approach in potential diagnostic application of urinary nucleosides in urogenital tract cancer

Aim: We aimed at evaluation the potential diagnostic role of urinary nucleosides in urogenital tract cancer. Materials & methods: Concentrations of 12 nucleosides determined by LC-MS/MS were subjected to correlation, association and interaction analyses. Results: We identified six pairs of nucleosides differently correlated in the group of patients and controls (p < 0.05). N-2-methylguanosine (odds ratio: 4.82; 95% CI: 1.78–12.93; p = 0.002) and N,N-dimethylguanosine (odds ratio: 5.45; 95% CI: 1.78–16.44; p = 0.003), were significantly associated with the disease risk (p-corrected = 0.004). Interaction between N-2-methylguanosine and adenosine (p-interaction = 0.019) suggested their multiplicative effect on the outcome. Conclusion: Urinary nucleosides, namely N,N-dimethylguanosine and N-2-methylguanosine may have the potential to serve as prognostic biomarkers. Gender-specific differences in urogenital tract cancer are likely to occur.

Systemic alterations in the metabolome of diabetic NOD mice delineate increased oxidative stress accompanied by reduced inflammation and hypertriglyceremia

Nonobese diabetic (NOD) mice are a commonly used model of type 1 diabetes (T1D). However, not all animals will develop overt diabetes despite undergoing similar autoimmune insult. In this study, a comprehensive metabolomic approach, consisting of gas chromatography time-of-flight (GC-TOF) mass spectrometry (MS), ultra-high-performance liquid chromatography-accurate mass quadruple time-of-flight (UHPLC-qTOF) MS and targeted UHPLC-tandem mass spectrometry-based methodologies, was used to capture metabolic alterations in the metabolome and lipidome of plasma from NOD mice progressing or not progressing to T1D. Using this multi-platform approach, we identified >1,000 circulating lipids and metabolites in male and female progressor and nonprogressor animals (n = 71). Statistical and multivariate analyses were used to identify age- and sex-independent metabolic markers, which best differentiated metabolic profiles of progressors and nonprogressors. Key T1D-associated perturbations were related with 1) increases in oxidation products glucono-δ-lactone and galactonic acid and reductions in cysteine, methionine and threonic acid, suggesting increased oxidative stress; 2) reductions in circulating polyunsaturated fatty acids and lipid signaling mediators, most notably arachidonic acid (AA) and AA-derived eicosanoids, implying impaired states of systemic inflammation; 3) elevations in circulating triacylglyercides reflective of hypertriglyceridemia; and 4) reductions in major structural lipids, most notably lysophosphatidylcholines and phosphatidylcholines. Taken together, our results highlight the systemic perturbations that accompany a loss of glycemic control and development of overt T1D.

Temporal metabolomic responses of cultured HepG2 liver cells to high fructose and high glucose exposures

High fructose consumption has been implicated with deleterious effects on human health, including hyperlipidemia elicited through de novo lipogenesis. However, more global effects of fructose on cellular metabolism have not been elucidated. In order to explore the metabolic impact of fructose-containing nutrients, we applied both GC-TOF and HILIC-QTOF mass spectrometry metabolomic strategies using extracts from cultured HepG2 cells exposed to fructose, glucose, or fructose + glucose. Cellular responses were analyzed in a time-dependent manner, incubated in media containing 5.5 mM glucose + 5.0 mM fructose in comparison to controls incubated in media containing either 5.5 mM glucose or 10.5 mM glucose. Mass spectrometry identified 156 unique known metabolites and a large number of unknown compounds, which revealed metabolite changes due to both utilization of fructose and high-carbohydrate loads independent of hexose structure. Fructose was shown to be partially converted to sorbitol, and generated higher levels of fructose-1-phosphate as a precursor for glycolytic intermediates. Differentially regulated ratios of 3-phosphoglycerate to serine pathway intermediates in high fructose media indicated a diversion of carbon backbones away from energy metabolism. Additionally, high fructose conditions changed levels of complex lipids toward phosphatidylethanolamines. Patterns of acylcarnitines in response to high hexose exposure (10.5 mM glucose or glucose/fructose combination) suggested a reduction in mitochondrial beta-oxidation.

Systemic Metabolomic Changes in Blood Samples of Lung Cancer Patients Identified by Gas Chromatography Time-of-Flight Mass Spectrometry

Lung cancer is a leading cause of cancer deaths worldwide. Metabolic alterations in tumor cells coupled with systemic indicators of the host response to tumor development have the potential to yield blood profiles with clinical utility for diagnosis and monitoring of treatment. We report results from two separate studies using gas chromatography time-of-flight mass spectrometry (GC-TOF MS) to profile metabolites in human blood samples that significantly differ from non-small cell lung cancer (NSCLC) adenocarcinoma and other lung cancer cases. Metabolomic analysis of blood samples from the two studies yielded a total of 437 metabolites, of which 148 were identified as known compounds and 289 identified as unknown compounds. Differential analysis identified 15 known metabolites in one study and 18 in a second study that were statistically different (p-values <0.05). Levels of maltose, palmitic acid, glycerol, ethanolamine, glutamic acid, and lactic acid were increased in cancer samples while amino acids tryptophan, lysine and histidine decreased. Many of the metabolites were found to be significantly different in both studies, suggesting that metabolomics appears to be robust enough to find systemic changes from lung cancer, thus showing the potential of this type of analysis for lung cancer detection.

MetaMapR: pathway independent metabolomic network analysis incorporating unknowns

Metabolic network mapping is a widely used approach for integration of metabolomic experimental results with biological domain knowledge. However, current approaches can be limited by biochemical domain or pathway knowledge which results in sparse disconnected graphs for real world metabolomic experiments. MetaMapR integrates enzymatic transformations with metabolite structural similarity, mass spectral similarity and empirical associations to generate richly connected metabolic networks. This open source, web-based or desktop software, written in the R programming language, leverages KEGG and PubChem databases to derive associations between metabolites even in cases where biochemical domain or molecular annotations are unknown. Network calculation is enhanced through an interface to the Chemical Translation System, which allows metabolite identifier translation between >200 common biochemical databases. Analysis results are presented as interactive visualizations or can be exported as high-quality graphics and numerical tables which can be imported into common network analysis and visualization tools.

AVAILABILITY AND IMPLEMENTATION

Freely available at http://dgrapov.github.io/MetaMapR/. Requires R and a modern web browser. Installation instructions, tutorials and application examples are available at http://dgrapov.github.io/MetaMapR/.

CONTACT

ofiehn@ucdavis.edu.

Diabetes Associated Metabolomic Perturbations in NOD Mice

Non-obese diabetic (NOD) mice are a widely-used model oftype1 diabetes (T1D). However, not all animals develop overt diabetes. This study examined the circulating metabolomic profiles of NOD mice progressing or not progressing to T1D. Total beta-cell mass was quantified in the intact pancreas using transgenic NOD mice expressinggreen fluorescent protein under the control of mouse insulin I promoter.While both progressor and non-progressor animals displayed lymphocyte infiltration and endoplasmic reticulum stress in the pancreas tissue;overt T1D did not develop until animals lost ~70% of the total beta-cell mass.Gas chromatography time of flight mass spectrometry (GC-TOF) was used to measure >470 circulating metabolites in male and female progressor and non-progressor animals (n=76) across a wide range of ages (neonates to >40-wk).Statistical and multivariate analyses were used to identify age and sex independent metabolic markers which best differentiated progressor and non-progressor animals' metabolic profiles. Key T1D-associated perturbations were related with: (1) increased plasma glucose and reduced 1,5-anhydroglucitol markers of glycemic control; (2) increased allantoin, gluconic acid and nitric oxide-derived saccharic acid markers of oxidative stress; (3) reduced lysine, an insulin secretagogue; (4) increased branched-chain amino acids, isoleucine and valine; (5) reduced unsaturated fatty acids including arachidonic acid; and (6)perturbations in urea cycle intermediates suggesting increased arginine-dependent NO synthesis. Together these findings highlight the strength of the unique approach of comparing progressor and non-progressor NOD mice to identify metabolic perturbations involved in T1D progression.

Large-scale chemical similarity networks for target profiling of compounds identified in cell-based chemical screens

Target identification is one of the most critical steps following cell-based phenotypic chemical screens aimed at identifying compounds with potential uses in cell biology and for developing novel disease therapies. Current in silico target identification methods, including chemical similarity database searches, are limited to single or sequential ligand analysis that have limited capabilities for accurate deconvolution of a large number of compounds with diverse chemical structures. Here, we present CSNAP (Chemical Similarity Network Analysis Pulldown), a new computational target identification method that utilizes chemical similarity networks for large-scale chemotype (consensus chemical pattern) recognition and drug target profiling. Our benchmark study showed that CSNAP can achieve an overall higher accuracy (>80%) of target prediction with respect to representative chemotypes in large (>200) compound sets, in comparison to the SEA approach (60–70%). Additionally, CSNAP is capable of integrating with biological knowledge-based databases (Uniprot, GO) and high-throughput biology platforms (proteomic, genetic, etc) for system-wise drug target validation. To demonstrate the utility of the CSNAP approach, we combined CSNAP's target prediction with experimental ligand evaluation to identify the major mitotic targets of hit compounds from a cell-based chemical screen and we highlight novel compounds targeting microtubules, an important cancer therapeutic target. The CSNAP method is freely available and can be accessed from the CSNAP web server (http://services.mbi.ucla.edu/CSNAP/).

Determining the targets of compounds identified in cell-based high-throughput chemical screens is a critical step for downstream drug development and understanding of compound mechanism of action. However, current computational target prediction approaches like chemical similarity database searches are limited to single or sequential ligand analyses, which limits their ability to accurately deconvolve a large number of compounds that often have chemically diverse structures. Here, we have developed a new computational drug target prediction method, called CSNAP that is based on chemical similarity networks. By clustering diverse chemical structures into distinct sub-networks corresponding to chemotypes, we show that CSNAP improves target prediction accuracy and consistency over a board range of drug classes. We further coupled CSNAP to a mitotic database and successfully determined the major mitotic drug targets of a diverse compound set identified in a cell-based chemical screen. We demonstrate that CSNAP can easily integrate with diverse knowledge-based databases for on/off target prediction and post-target validation, thus broadening its applicability for identifying the targets of bioactive compounds from a wide range of chemical screens.

The human toxome project

The Human Toxome Project, funded as an NIH Transformative Research grant 2011-2016, is focused on developing the concepts and the means for deducing, validating and sharing molecular pathways of toxicity (PoT). Using the test case of estrogenic endocrine disruption, the responses of MCF-7 human breast cancer cells are being phenotyped by transcriptomics and mass-spectroscopy-based metabolomics. The bioinformatics tools for PoT deduction represent a core deliverable. A number of challenges for quality and standardization of cell systems, omics technologies and bioinformatics are being addressed. In parallel, concepts for annotation, validation and sharing of PoT information, as well as their link to adverse outcomes, are being developed. A reasonably comprehensive public database of PoT, the Human Toxome Knowledge-base, could become a point of reference for toxicological research and regulatory test strategies.

Metabolomics and diabetes: analytical and computational approaches

Diabetes is characterized by altered metabolism of key molecules and regulatory pathways. The phenotypic expression of diabetes and associated complications encompasses complex interactions between genetic, environmental, and tissue-specific factors that require an integrated understanding of perturbations in the network of genes, proteins, and metabolites. Metabolomics attempts to systematically identify and quantitate small molecule metabolites from biological systems. The recent rapid development of a variety of analytical platforms based on mass spectrometry and nuclear magnetic resonance have enabled identification of complex metabolic phenotypes. Continued development of bioinformatics and analytical strategies has facilitated the discovery of causal links in understanding the pathophysiology of diabetes and its complications. Here, we summarize the metabolomics workflow, including analytical, statistical, and computational tools, highlight recent applications of metabolomics in diabetes research, and discuss the challenges in the field.

Metabolomic markers of altered nucleotide metabolism in early stage adenocarcinoma

Adenocarcinoma, a type of non-small cell lung cancer, is the most frequently diagnosed lung cancer and the leading cause of lung cancer mortality in the United States. It is well documented that biochemical changes occur early in the transition from normal to cancer cells, but the extent to which these alterations affect tumorigenesis in adenocarcinoma remains largely unknown. Herein, we describe the application of mass spectrometry and multivariate statistical analysis in one of the largest biomarker research studies to date aimed at distinguishing metabolic differences between malignant and nonmalignant lung tissue. Gas chromatography time-of-flight mass spectrometry was used to measure 462 metabolites in 39 malignant and nonmalignant lung tissue pairs from current or former smokers with early stage (stage IA-IB) adenocarcinoma. Statistical mixed effects models, orthogonal partial least squares discriminant analysis and network integration, were used to identify key cancer-associated metabolic perturbations in adenocarcinoma compared with nonmalignant tissue. Cancer-associated biochemical alterations were characterized by (i) decreased glucose levels, consistent with the Warburg effect, (ii) changes in cellular redox status highlighted by elevations in cysteine and antioxidants, alpha- and gamma-tocopherol, (iii) elevations in nucleotide metabolites 5,6-dihydrouracil and xanthine suggestive of increased dihydropyrimidine dehydrogenase and xanthine oxidoreductase activity, (iv) increased 5'-deoxy-5'-methylthioadenosine levels indicative of reduced purine salvage and increased de novo purine synthesis, and (v) coordinated elevations in glutamate and UDP-N-acetylglucosamine suggesting increased protein glycosylation. The present study revealed distinct metabolic perturbations associated with early stage lung adenocarcinoma, which may provide candidate molecular targets for personalizing therapeutic interventions and treatment efficacy monitoring.

Alteration of the fecal microbiota and serum metabolite profiles in dogs with idiopathic inflammatory bowel disease

Idiopathic inflammatory bowel disease (IBD) is a common cause of chronic gastrointestinal (GI) disease in dogs. The combination of an underlying host genetic susceptibility, an intestinal dysbiosis, and dietary/environmental factors are suspected as main contributing factors in the pathogenesis of canine IBD. However, actual mechanisms of the host-microbe interactions remain elusive. The aim of this study was to compare the fecal microbiota and serum metabolite profiles between healthy dogs (n = 10) and dogs with IBD before and after 3 weeks of medical therapy (n = 12). Fecal microbiota and metabolite profiles were characterized by 454-pyrosequencing of 16 S rRNA genes and by an untargeted metabolomics approach, respectively. Significantly lower bacterial diversity and distinct microbial communities were observed in dogs with IBD compared to the healthy control dogs. While Gammaproteobacteria were overrepresented, Erysipelotrichia, Clostridia, and Bacteroidia were underrepresented in dogs with IBD. The functional gene content was predicted from the 16 S rRNA gene data using PICRUSt, and revealed overrepresented bacterial secretion system and transcription factors, and underrepresented amino acid metabolism in dogs with IBD. The serum metabolites 3-hydroxybutyrate, hexuronic acid, ribose, and gluconic acid lactone were significantly more abundant in dogs with IBD. Although a clinical improvement was observed after medical therapy in all dogs with IBD, this was not accompanied by significant changes in the fecal microbiota or in serum metabolite profiles. These results suggest the presence of oxidative stress and a functional alteration of the GI microbiota in dogs with IBD, which persisted even in the face of a clinical response to medical therapy.

Principles and practice of lipidomics

The technical advances in mass spectrometry, particularly the development of (ultra)-high-resolution/mass accuracy measurement capabilities in combination with refinement of soft ionization techniques, have increased the application and success of lipidomics to answer biological questions in relation to lipid metabolism. Together with other omics technologies, lipidomics has become an important tool to practice systems biology as lipids comprise a very significant part of the metabolome and play pleiotropic roles in cellular functions. As an increasing number of disorders are linked to lipid metabolism, lipidomics is used to search for biomarkers, understand disease mechanism and follow the efficacy of therapeutic options. This review provides a first introduction to the major methodological strategies currently used for mass spectrometry-based lipidomics and associated data pre-processing and analysis.

Metabolomics in psoriatic disease: pilot study reveals metabolite differences in psoriasis and psoriatic arthritis

Importance: While “omics” studies have advanced our understanding of inflammatory skin diseases, metabolomics is mostly an unexplored field in dermatology.

Objective: We sought to elucidate the pathogenesis of psoriatic diseases by determining the differences in metabolomic profiles among psoriasis patients with or without psoriatic arthritis and healthy controls.

Design: We employed a global metabolomics approach to compare circulating metabolites from patients with psoriasis, psoriasis and psoriatic arthritis, and healthy controls.

Setting: Study participants were recruited from the general community and from the Psoriasis Clinic at the University of California Davis in United States.

Participants: We examined metabolomic profiles using blood serum samples from 30 patients age and gender matched into three groups: 10 patients with psoriasis, 10 patients with psoriasis and psoriatic arthritis and 10 control participants.

Main outcome(s) and measures(s): Metabolite levels were measured calculating the mean peak intensities from gas chromatography time-of-flight mass spectrometry.

Results: Multivariate analyses of metabolomics profiles revealed altered serum metabolites among the study population. Compared to control patients, psoriasis patients had a higher level of alpha ketoglutaric acid (Pso: 288 ± 88; Control: 209 ± 69; p=0.03), a lower level of asparagine (Pso: 5460 ± 980; Control: 7260 ± 2100; p=0.02), and a lower level of glutamine (Pso: 86000 ± 20000; Control: 111000 ± 27000; p=0.02). Compared to control patients, patients with psoriasis and psoriatic arthritis had increased levels of glucuronic acid (Pso + PsA: 638 ± 250; Control: 347 ± 61; p=0.001). Compared to patients with psoriasis alone, patients with both psoriasis and psoriatic arthritis had a decreased level of alpha ketoglutaric acid (Pso + PsA: 186 ± 80; Pso: 288 ± 88; p=0.02) and an increased level of lignoceric acid (Pso + PsA: 442 ± 280; Pso: 214 ± 64; p=0.02).

Conclusions and relevance: The metabolite differences help elucidate the pathogenesis of psoriasis and psoriatic arthritis and they may provide insights for therapeutic development.

The blood exposome and its role in discovering causes of disease

BACKGROUND

Since 2001, researchers have examined the human genome (G) mainly to discover causes of disease, despite evidence that G explains relatively little risk. We posit that unexplained disease risks are caused by the exposome (E; representing all exposures) and G × E interactions. Thus, etiologic research has been hampered by scientists' continuing reliance on low-tech methods to characterize E compared with high-tech omics for characterizing G.

OBJECTIVES

Because exposures are inherently chemical in nature and arise from both endogenous and exogenous sources, blood specimens can be used to characterize exposomes. To explore the "blood exposome" and its connection to disease, we sought human blood concentrations of many chemicals, along with their sources, evidence of chronic-disease risks, and numbers of metabolic pathways.

METHODS

From the literature we obtained human blood concentrations of 1,561 small molecules and metals derived from foods, drugs, pollutants, and endogenous processes. We mapped chemical similarities after weighting by blood concentrations, disease-risk citations, and numbers of human metabolic pathways.

RESULTS

Blood concentrations spanned 11 orders of magnitude and were indistinguishable for endogenous and food chemicals and drugs, whereas those of pollutants were 1,000 times lower. Chemical similarities mapped by disease risks were equally distributed by source categories, but those mapped by metabolic pathways were dominated by endogenous molecules and essential nutrients.

CONCLUSIONS

For studies of disease etiology, the complexity of human exposures motivates characterization of the blood exposome, which includes all biologically active chemicals. Because most small molecules in blood are not human metabolites, investigations of causal pathways should expand beyond the endogenous metabolome.

MetDisease--connecting metabolites to diseases via literature

MOTIVATION

In recent years, metabolomics has emerged as an approach to perform large-scale characterization of small molecules in biological systems. Metabolomics posed a number of bioinformatics challenges associated in data analysis and interpretation. Genome-based metabolic reconstructions have established a powerful framework for connecting metabolites to genes through metabolic reactions and enzymes that catalyze them. Pathway databases and bioinformatics tools that use this framework have proven to be useful for annotating experimental metabolomics data. This framework can be used to infer connections between metabolites and diseases through annotated disease genes. However, only about half of experimentally detected metabolites can be mapped to canonical metabolic pathways. We present a new Cytoscape 3 plug-in, MetDisease, which uses an alternative approach to link metabolites to disease information. MetDisease uses Medical Subject Headings (MeSH) disease terms mapped to PubChem compounds through literature to annotate compound networks.

AVAILABILITY AND IMPLEMENTATION

MetDisease can be downloaded from http://apps.cytoscape.org/apps/metdisease or installed via the Cytoscape app manager. Further information about MetDisease can be found at http://metdisease.ncibi.org

CONTACT

akarnovs@med.umich.edu

SUPPLEMENTARY INFORMATION

Supplementary Data are available at Bioinformatics online.

Serum metabolomics identifies citrulline as a predictor of adverse outcomes in an equine model of gut-derived sepsis

Acute laminitis is an inflammatory disease of the equine foot that often occurs secondarily to sepsis or systemic inflammation associated with gastrointestinal disease. It has been suggested that laminitis is similar to multiple organ dysfunction syndrome in humans, although in horses the weight-bearing laminar epithelium of the foot appears to be the tissue most sensitive to insult and the first "organ" to fail. Metabolomics performed on serum samples collected before (Con) and after (Lmn) experimental induction of gastrointestinal-associated sepsis in six horses detected 1,177 metabolites of both mammalian and bacterial origin in equine serum. Network and correlation analyses suggested a dysregulation of fatty acid metabolism in the Lmn group, as well as an accumulation of organic acids such as lactate. Furthermore, concentrations of the amino acid citrulline were decreased in Lmn samples from all study animals, suggesting that citrulline might be useful as a biomarker to identify critically ill animals that are at risk of developing laminitis. We therefore established normal ranges of plasma citrulline concentrations in a separate group of horses (n = 36) and tested the ability of citrulline to predict adverse outcomes (laminitis or death) in critically ill horses (n = 23). Plasma citrulline was significantly lower in critically ill horses that went on to experience adverse outcomes (n = 6). Further study is required to accurately determine a diagnostic cutoff, but the present data are suggestive of the predictive value of citrulline as a biomarker for laminar failure in equine sepsis.

Inactivation of metabolic genes causes short- and long-range dys-regulation in Escherichia coli metabolic network

The metabolic network in E. coli can be severely affected by the inactivation of metabolic genes that are required to catabolize a nutrient (D-galactose). We hypothesized that the resulting accumulation of small molecules can yield local as well as systemic effects on the metabolic network. Analysis of metabolomics data in wild-type and D-galactose non-utilizing mutants, galT, galU and galE, reveal the large metabolic differences between the wild-type and the mutants when the strains were grown in D-galactose. Network mapping suggested that the enzymatic defects affected the metabolic modules located both at short- and long-ranges from the D-galactose metabolic module. These modules suggested alterations in glutathione, energy, nucleotide and lipid metabolism and disturbed carbon to nitrogen ratio in mutant strains. The altered modules are required for normal cell growth for the wild-type strain, explaining why the cell growth is inhibited in the mutants in the presence of D-galactose. Identification of these distance-based dys-regulations would enhance the systems level understanding of metabolic networks of microorganisms having importance in biomedical and biotechnological research.

Glycogenomics as a mass spectrometry-guided genome-mining method for microbial glycosylated molecules

Glycosyl groups are an essential mediator of molecular interactions in cells and on cellular surfaces. There are very few methods that directly relate sugar-containing molecules to their biosynthetic machineries. Here, we introduce glycogenomics as an experiment-guided genome-mining approach for fast characterization of glycosylated natural products (GNPs) and their biosynthetic pathways from genome-sequenced microbes by targeting glycosyl groups in microbial metabolomes. Microbial GNPs consist of aglycone and glycosyl structure groups in which the sugar unit(s) are often critical for the GNP's bioactivity, e.g., by promoting binding to a target biomolecule. GNPs are a structurally diverse class of molecules with important pharmaceutical and agrochemical applications. Herein, O- and N-glycosyl groups are characterized in their sugar monomers by tandem mass spectrometry (MS) and matched to corresponding glycosylation genes in secondary metabolic pathways by a MS-glycogenetic code. The associated aglycone biosynthetic genes of the GNP genotype then classify the natural product to further guide structure elucidation. We highlight the glycogenomic strategy by the characterization of several bioactive glycosylated molecules and their gene clusters, including the anticancer agent cinerubin B from Streptomyces sp. SPB74 and an antibiotic, arenimycin B, from Salinispora arenicola CNB-527.

Human metabolomics: strategies to understand biology

Metabolomics provides a direct functional read-out of the physiological status of an organism and is in principle ideally suited to describe someone's health status. Whereas only a limited number of small metabolites are used in the clinics, in inborn errors of metabolism an extensive repertoire of metabolites are used as biomarkers. We discuss that the proper clinical phenotyping is crucial to find biomarkers and obtain biological insights for multifactorial diseases. This requires to study the phenotype dynamics including the concepts of homeostasis and allostasis, that is, the ability to adapt and cope with a challenge. We also elaborate that biology-driven metabolomics platforms (i.e. development of metabolomics technology driven by the need of studying and answering important biomedical questions) addressing clinically relevant pathways and at the same time providing absolute concentrations are key to allow discovery and validation of biomarkers across studies and labs. Following individuals over years will require high throughput metabolomics approaches, which are emerging for nuclear magnetic resonance spectroscopy and direct-infusion mass spectrometry, but should also include the biochemical networks needed for personalized health monitoring.

Pharmacometabolomics reveals racial differences in response to atenolol treatment

Antihypertensive drugs are among the most commonly prescribed drugs for chronic disease worldwide. The response to antihypertensive drugs varies substantially between individuals and important factors such as race that contribute to this heterogeneity are poorly understood. In this study we use metabolomics, a global biochemical approach to investigate biochemical changes induced by the beta-adrenergic receptor blocker atenolol in Caucasians and African Americans. Plasma from individuals treated with atenolol was collected at baseline (untreated) and after a 9 week treatment period and analyzed using a GC-TOF metabolomics platform. The metabolomic signature of atenolol exposure included saturated (palmitic), monounsaturated (oleic, palmitoleic) and polyunsaturated (arachidonic, linoleic) free fatty acids, which decreased in Caucasians after treatment but were not different in African Americans (p<0.0005, q<0.03). Similarly, the ketone body 3-hydroxybutyrate was significantly decreased in Caucasians by 33% (p<0.0001, q<0.0001) but was unchanged in African Americans. The contribution of genetic variation in genes that encode lipases to the racial differences in atenolol-induced changes in fatty acids was examined. SNP rs9652472 in LIPC was found to be associated with the change in oleic acid in Caucasians (p<0.0005) but not African Americans, whereas the PLA2G4C SNP rs7250148 associated with oleic acid change in African Americans (p<0.0001) but not Caucasians. Together, these data indicate that atenolol-induced changes in the metabolome are dependent on race and genotype. This study represents a first step of a pharmacometabolomic approach to phenotype patients with hypertension and gain mechanistic insights into racial variability in changes that occur with atenolol treatment, which may influence response to the drug.

Computational tools for the secondary analysis of metabolomics experiments

Metabolomics experiments have become commonplace in a wide variety of disciplines. By identifying and quantifying metabolites researchers can achieve a systems level understanding of metabolism. These studies produce vast swaths of data which are often only lightly interpreted due to the overwhelmingly large amount of variables that are measured. Recently, a number of computational tools have been developed which enable much deeper analysis of metabolomics data. These data have been difficult to interpret as understanding the connections between dozens of altered metabolites has often relied on the biochemical knowledge of researchers and their speculations. Modern biochemical databases provide information about the interconnectivity of metabolism which can be automatically polled using metabolomics secondary analysis tools. Starting with lists of altered metabolites, there are two main types of analysis: enrichment analysis computes which metabolic pathways have been significantly altered whereas metabolite mapping contextualizes the abundances and significances of measured metabolites into network visualizations. Many different tools have been developed for one or both of these applications. In this review the functionality and use of these software is discussed. Together these novel secondary analysis tools will enable metabolomics researchers to plumb the depths of their data and produce farther reaching biological conclusions than ever before.

Prioritization of putative metabolite identifications in LC-MS/MS experiments using a computational pipeline

One of the major bottle-necks in current LC-MS-based metabolomic investigations is metabolite identification. An often-used approach is to first look up metabolites from databases through peak mass, followed by verification of the obtained putative identifications using MS/MS data. However, the mass-based search may provide inappropriate putative identifications when the observed peak is from isotopes, fragments, or adducts. In addition, a large fraction of peaks is often left with multiple putative identifications. To differentiate these putative identifications, manual verification of metabolites through comparison between biological samples and authentic compounds is necessary. However, such experiments are laborious, especially when multiple putative identifications are encountered. It is desirable to use computational approaches to obtain more reliable putative identifications and prioritize them before performing experimental verification of the metabolites. In this article, a computational pipeline is proposed to assist metabolite identification with improved metabolome coverage and prioritization capability. Multiple publicly available software tools and databases, along with in-house developed algorithms, are utilized to fully exploit the information acquired from LC-MS/MS experiments. The pipeline is successfully applied to identify metabolites on the basis of LC-MS as well as MS/MS data. Using accurate masses, retention time values, MS/MS spectra, and metabolic pathways/networks, more appropriate putative identifications are retrieved and prioritized to guide subsequent metabolite verification experiments.

Recent progress in the development of metabolome databases for plant systems biology

Metabolomics has grown greatly as a functional genomics tool, and has become an invaluable diagnostic tool for biochemical phenotyping of biological systems. Over the past decades, a number of databases involving information related to mass spectra, compound names and structures, statistical/mathematical models and metabolic pathways, and metabolite profile data have been developed. Such databases complement each other and support efficient growth in this area, although the data resources remain scattered across the World Wide Web. Here, we review available metabolome databases and summarize the present status of development of related tools, particularly focusing on the plant metabolome. Data sharing discussed here will pave way for the robust interpretation of metabolomic data and advances in plant systems biology.

Induced pluripotent stem cells show metabolomic differences to embryonic stem cells in polyunsaturated phosphatidylcholines and primary metabolism

Induced pluripotent stem cells are different from embryonic stem cells as shown by epigenetic and genomics analyses. Depending on cell types and culture conditions, such genetic alterations can lead to different metabolic phenotypes which may impact replication rates, membrane properties and cell differentiation. We here applied a comprehensive metabolomics strategy incorporating nanoelectrospray ion trap mass spectrometry (MS), gas chromatography-time of flight MS, and hydrophilic interaction- and reversed phase-liquid chromatography-quadrupole time-of-flight MS to examine the metabolome of induced pluripotent stem cells (iPSCs) compared to parental fibroblasts as well as to reference embryonic stem cells (ESCs). With over 250 identified metabolites and a range of structurally unknown compounds, quantitative and statistical metabolome data were mapped onto a metabolite networks describing the metabolic state of iPSCs relative to other cell types. Overall iPSCs exhibited a striking shift metabolically away from parental fibroblasts and toward ESCs, suggestive of near complete metabolic reprogramming. Differences between pluripotent cell types were not observed in carbohydrate or hydroxyl acid metabolism, pentose phosphate pathway metabolites, or free fatty acids. However, significant differences between iPSCs and ESCs were evident in phosphatidylcholine and phosphatidylethanolamine lipid structures, essential and non-essential amino acids, and metabolites involved in polyamine biosynthesis. Together our findings demonstrate that during cellular reprogramming, the metabolome of fibroblasts is also reprogrammed to take on an ESC-like profile, but there are select unique differences apparent in iPSCs. The identified metabolomics signatures of iPSCs and ESCs may have important implications for functional regulation of maintenance and induction of pluripotency.