KNApSAcK family databases: integrated metabolite-plant species databases for multifaceted plant research

A database (DB) describing the relationships between species and their metabolites would be useful for metabolomics research, because it targets systematic analysis of enormous numbers of organic compounds with known or unknown structures in metabolomics. We constructed an extensive species-metabolite DB for plants, the KNApSAcK Core DB, which contains 101,500 species-metabolite relationships encompassing 20,741 species and 50,048 metabolites. We also developed a search engine within the KNApSAcK Core DB for use in metabolomics research, making it possible to search for metabolites based on an accurate mass, molecular formula, metabolite name or mass spectra in several ionization modes. We also have developed databases for retrieving metabolites related to plants used for a range of purposes. In our multifaceted plant usage DB, medicinal/edible plants are related to the geographic zones (GZs) where the plants are used, their biological activities, and formulae of Japanese and Indonesian traditional medicines (Kampo and Jamu, respectively). These data are connected to the species-metabolites relationship DB within the KNApSAcK Core DB, keyed via the species names. All databases can be accessed via the website http://kanaya.naist.jp/KNApSAcK_Family/. KNApSAcK WorldMap DB comprises 41,548 GZ-plant pair entries, including 222 GZs and 15,240 medicinal/edible plants. The KAMPO DB consists of 336 formulae encompassing 278 medicinal plants; the JAMU DB consists of 5,310 formulae encompassing 550 medicinal plants. The Biological Activity DB consists of 2,418 biological activities and 33,706 pairwise relationships between medicinal plants and their biological activities. Current statistics of the binary relationships between individual databases were characterized by the degree distribution analysis, leading to a prediction of at least 1,060,000 metabolites within all plants. In the future, the study of metabolomics will need to take this huge number of metabolites into consideration.     

FT-ICR/MS and GC-EI/MS Metabolomics Networking Unravels Global Potato Sprout's Responses to Rhizoctonia solani Infection

The complexity of plant-pathogen interactions makes their dissection a challenging task for metabolomics studies. Here we are reporting on an integrated metabolomics networking approach combining gas chromatography/mass spectrometry (GC/MS) with Fourier transform ion cyclotron resonance/mass spectrometry (FT-ICR/MS) and bioinformatics analyses for the study of interactions in the potato sprout-Rhizoctonia solani pathosystem and the fluctuations in the global metabolome of sprouts. The developed bioanalytical and bioinformatics protocols provided a snapshot of the sprout's global metabolic network and its perturbations as a result of pathogen invasion. Mevalonic acid and deoxy-xylulose pathways were substantially up-regulated leading to the biosynthesis of sesquiterpene alkaloids such as the phytoalexins phytuberin, rishitin, and solavetivone, and steroidal alkaloids having solasodine and solanidine as their common aglycons. Additionally, the perturbation of the sprout's metabolism was depicted in fluctuations of the content of their amino acids pool and that of carboxylic and fatty acids. Components of the systemic acquired resistance (SAR) and hypersensitive reaction (HR) such as azelaic and oxalic acids were detected in increased levels in infected sprouts and strategies of the pathogen to overcome plant defense were proposed. Our metabolic approach has not only greatly expanded the multitude of metabolites previously reported in potato in response to pathogen invasion, but also enabled the identification of bioactive plant-derived metabolites providing valuable information that could be exploited in biotechnology, biomarker-assisted plant breeding, and crop protection for the development of new crop protection agents.     

Optimization and use of peptide mass measurement accuracy in shotgun proteomics

Mass spectrometers that provide high mass accuracy such as FT-ICR instruments are increasingly used in proteomic studies. Although the importance of accurately determined molecular masses for the identification of biomolecules is generally accepted, its role in the analysis of shotgun proteomic data has not been thoroughly studied. To gain insight into this role, we used a hybrid linear quadrupole ion trap/FT-ICR (LTQ FT) mass spectrometer for LC-MS/MS analysis of a highly complex peptide mixture derived from a fraction of the yeast proteome. We applied three data-dependent MS/MS acquisition methods. The FT-ICR part of the hybrid mass spectrometer was either not exploited, used only for survey MS scans, or also used for acquiring selected ion monitoring scans to optimize mass accuracy. MS/MS data were assigned with the SEQUEST algorithm, and peptide identifications were validated by estimating the number of incorrect assignments using the composite target/decoy database search strategy. We developed a simple mass calibration strategy exploiting polydimethylcyclosiloxane background ions as calibrant ions. This strategy allowed us to substantially improve mass accuracy without reducing the number of MS/MS spectra acquired in an LC-MS/MS run. The benefits of high mass accuracy were greatest for assigning MS/MS spectra with low signal-to-noise ratios and for assigning phosphopeptides. Confident peptide identification rates from these data sets could be doubled by the use of mass accuracy information. It was also shown that improving mass accuracy at a cost to the MS/MS acquisition rate substantially lowered the sensitivity of LC-MS/MS analyses. The use of FT-ICR selected ion monitoring scans to maximize mass accuracy reduced the number of protein identifications by 40%.     

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.     

Application of a metabolomic method combining one-dimensional and two-dimensional gas chromatography-time-of-flight/mass spectrometry to metabolic phenotyping of natural variants in rice

We have developed a comprehensive method combining analytical techniques of one-dimensional (1D) and two-dimensional (GC x GC) gas chromatography-time-of-flight (TOF)-mass spectrometry. This method was applied to the metabolic phenotyping of natural variants in rice for the 68 world rice core collection (WRC) and two other varieties. Ten metabolites were selected as metabolite representatives, and the selected ion current of each metabolite peak obtained from both techniques were statistically compared. Our method of combining 1D- and GC x GC-TOF/MS is useful for the metabolic phenotyping of natural variants in rice for further studies in breeding programs.     

Metabolite profiling for plant functional genomics

Multiparallel analyses of mRNA and proteins are central to today's functional genomics initiatives. We describe here the use of metabolite profiling as a new tool for a comparative display of gene function. It has the potential not only to provide deeper insight into complex regulatory processes but also to determine phenotype directly. Using gas chromatography/mass spectrometry (GC/MS), we automatically quantified 326 distinct compounds from Arabidopsis thaliana leaf extracts. It was possible to assign a chemical structure to approximately half of these compounds. Comparison of four Arabidopsis genotypes (two homozygous ecotypes and a mutant of each ecotype) showed that each genotype possesses a distinct metabolic profile. Data mining tools such as principal component analysis enabled the assignment of "metabolic phenotypes" using these large data sets. The metabolic phenotypes of the two ecotypes were more divergent than were the metabolic phenotypes of the single-loci mutant and their parental ecotypes. These results demonstrate the use of metabolite profiling as a tool to significantly extend and enhance the power of existing functional genomics approaches.     

Granger causality in integrated GC–MS and LC–MS metabolomics data reveals the interface of primary and secondary metabolism

Metabolomics has emerged as a key technique of modern life sciences in recent years. Two major techniques for metabolomics in the last 10 years are gas chromatography coupled to mass spectrometry (GC–MS) and liquid chromatography coupled to mass spectrometry (LC–MS). Each platform has a specific performance detecting subsets of metabolites. GC–MS in combination with derivatisation has a preference for small polar metabolites covering primary metabolism. In contrast, reversed phase LC–MS covers large hydrophobic metabolites predominant in secondary metabolism. Here, we present an integrative metabolomics platform providing a mean to reveal the interaction of primary and secondary metabolism in plants and other organisms. The strategy combines GC–MS and LC–MS analysis of the same sample, a novel alignment tool MetMAX and a statistical toolbox COVAIN for data integration and linkage of Granger Causality with metabolic modelling. For metabolic modelling we have implemented the combined GC–LC–MS metabolomics data covariance matrix and a stoichiometric matrix of the underlying biochemical reaction network. The changes in biochemical regulation are expressed as differential Jacobian matrices. Applying the Granger causality, a subset of secondary metabolites was detected with significant correlations to primary metabolites such as sugars and amino acids. These metabolic subsets were compiled into a stoichiometric matrix N. Using N the inverse calculation of a differential Jacobian J from metabolomics data was possible. Key points of regulation at the interface of primary and secondary metabolism were identified.     

Recent progress in mass spectrometry for single-cell metabolomics

Metabolites are the end products of cellular vital activities and can reflect the state of cellular to a certain extent. Rapid change of metabolites and the low abundance of signature metabolites cause difficulties in single-cell detection, which is a great challenge in single-cell metabolomics analysis. Mass spectrometry (MS) is a powerful tool that uniquely suited to detect intracellular small-molecule metabolites and has shown good application in single-cell metabolite analysis. In this mini-review, we describe three types of emerging technologies for MS-based single-cell metabolic analysis in recent years, including nano-ESI-MS based single-cell metabolomics analysis, high-throughput analysis via flow cytometry, and cellular metabolic imaging analysis. These techniques provide a large amount of single-cell metabolic data, allowing the potential of MS in single-cell metabolic analysis is gradually being explored and is of great importance in disease and life science research.     

Metabolic source isotopic pair labeling and genome-wide association are complementary tools for the identification of metabolite–gene associations in plants

The optimal extraction of information from untargeted metabolomics analyses is a continuing challenge. Here, we describe an approach that combines stable isotope labeling, liquid chromatography– mass spectrometry (LC–MS), and a computational pipeline to automatically identify metabolites produced from a selected metabolic precursor. We identified the subset of the soluble metabolome generated from phenylalanine (Phe) in Arabidopsis thaliana, which we refer to as the Phe-derived metabolome (FDM) In addition to identifying Phe-derived metabolites present in a single wild-type reference accession, the FDM was established in nine enzymatic and regulatory mutants in the phenylpropanoid pathway. To identify genes associated with variation in Phe-derived metabolites in Arabidopsis, MS features collected by untargeted metabolite profiling of an Arabidopsis diversity panel were retrospectively annotated to the FDM and natural genetic variants responsible for differences in accumulation of FDM features were identified by genome-wide association. Large differences in Phe-derived metabolite accumulation and presence/absence variation of abundant metabolites were observed in the nine mutants as well as between accessions from the diversity panel. Many Phe-derived metabolites that accumulated in mutants also accumulated in non-Col-0 accessions and was associated to genes with known or suspected functions in the phenylpropanoid pathway as well as genes with no known functions. Overall, we show that cataloguing a biochemical pathway’s products through isotopic labeling across genetic variants can substantially contribute to the identification of metabolites and genes associated with their biosynthesis.

An isotopic labeling and LC–MS pipeline to identify metabolites produced from Phe and its integration with genome-wide association identifies genes associated with the phenylpropanoid pathway.     

Assessment of Metabolome Annotation Quality: A Method for Evaluating the False Discovery Rate of Elemental Composition Searches

Background

In metabolomics researches using mass spectrometry (MS), systematic searching of high-resolution mass data against compound databases is often the first step of metabolite annotation to determine elemental compositions possessing similar theoretical mass numbers. However, incorrect hits derived from errors in mass analyses will be included in the results of elemental composition searches. To assess the quality of peak annotation information, a novel methodology for false discovery rates (FDR) evaluation is presented in this study. Based on the FDR analyses, several aspects of an elemental composition search, including setting a threshold, estimating FDR, and the types of elemental composition databases most reliable for searching are discussed.

Methodology/Principal Findings

The FDR can be determined from one measured value (i.e., the hit rate for search queries) and four parameters determined by Monte Carlo simulation. The results indicate that relatively high FDR values (30–50%) were obtained when searching time-of-flight (TOF)/MS data using the KNApSAcK and KEGG databases. In addition, searches against large all-in-one databases (e.g., PubChem) always produced unacceptable results (FDR >70%). The estimated FDRs suggest that the quality of search results can be improved not only by performing more accurate mass analysis but also by modifying the properties of the compound database. A theoretical analysis indicates that FDR could be improved by using compound database with smaller but higher completeness entries.

Conclusions/Significance

High accuracy mass analysis, such as Fourier transform (FT)-MS, is needed for reliable annotation (FDR <10%). In addition, a small, customized compound database is preferable for high-quality annotation of metabolome data.     

Untargeted Plasma Metabolomics Identifies Endogenous Metabolite with Drug-like Properties in Chronic Animal Model of Multiple Sclerosis

We performed untargeted metabolomics in plasma of B6 mice with experimental autoimmune encephalitis (EAE) at the chronic phase of the disease in search of an altered metabolic pathway(s). Of 324 metabolites measured, 100 metabolites that mapped to various pathways (mainly lipids) linked to mitochondrial function, inflammation, and membrane stability were observed to be significantly altered between EAE and control (p < 0.05, false discovery rate <0.10). Bioinformatics analysis revealed six metabolic pathways being impacted and altered in EAE, including α-linolenic acid and linoleic acid metabolism (PUFA). The metabolites of PUFAs, including ω-3 and ω-6 fatty acids, are commonly decreased in mouse models of multiple sclerosis (MS) and in patients with MS. Daily oral administration of resolvin D1, a downstream metabolite of ω-3, decreased disease progression by suppressing autoreactive T cells and inducing an M2 phenotype of monocytes/macrophages and resident brain microglial cells. This study provides a proof of principle for the application of metabolomics to identify an endogenous metabolite(s) possessing drug-like properties, which is assessed for therapy in preclinical mouse models of MS.     

The yeast metabolome addressed by electrospray ionization mass spectrometry: Initiation of a mass spectral library and its applications for metabolic footprinting by direct infusion mass spectrometry

Mass spectrometry (MS) has been a major driver for metabolomics, and gas chromatography (GC)-MS has been one of the primary techniques used for microbial metabolomics. The use of liquid chromatography (LC)-MS has however been limited, but electrospray ionization (ESI) is very well suited for ionization of microbial metabolites without any previous derivatization needed. To address the capabilities of ESI-MS in detecting the metabolome of Saccharomyces cerevisiae, the in silico metabolome of this organism was used as a template to present a theoretical metabolome. This showed that in combination with the specificity of MS up to 84% of the metabolites can be identified in a high mass accuracy ESI-spectrum. A total of 66 metabolites were systematically analyzed by positive and negative ESI-MS/MS with the aim of initiating a spectral library for ESI of microbial metabolites. This systematic analysis gave insight into the ionization and fragmentation characteristics of the different metabolites. With this insight, a small study of metabolic footprinting with ESI-MS demonstrated that biological information can be extracted from footprinting spectra. Statistical analysis of the footprinting data revealed discriminating ions, which could be assigned using the in silico metabolome. By this approach metabolic footprinting can advance from a classification method that is used to derive biological information based on guilt-by-association, to a tool for extraction of metabolic differences, which can guide new targeted biological experiments. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Microbial metabolism of catechin stereoisomers by human faecal microbiota: Comparison of targeted analysis and a non-targeted metabolomics method

Microbial metabolism of the stereoisomers (+)-catechin and (−)-epicatechin was compared by two analytical techniques, GC/MS for quantitative targeted analysis and GC×GC-TOF for global characterization of the metabolome, using human faecal microbiota as an inoculum of converting microbiota. The ring-fission site changed when the inocula originated from two different groups of donors, but dehydroxylation progressed similarly regardless of the inoculum. Whereas GC/MS proved to be an appropriate tool for the study of specific expected metabolites of catechin stereoisomers, GC×GC-TOF-based metabolomics analysis also revealed new metabolites not included in the targeted analyses. Quantitation and verification of identification can also be performed in a metabolomics platform, if authentic standards are available.