mzMatch-ISO: an R tool for the annotation and relative quantification of isotope-labelled mass spectrometry data

Systems-Level Annotation of a Metabolomics Data Set Reduces 25 000 Features to Fewer than 1000 Unique Metabolites

When using liquid chromatography/mass spectrometry (LC/MS) to perform untargeted metabolomics, it is now routine to detect tens of thousands of features from biological samples. Poor understanding of the data, however, has complicated interpretation and masked the number of unique metabolites actually being measured in an experiment. Here we place an upper bound on the number of unique metabolites detected in Escherichia coli samples analyzed with one untargeted metabolomics method. We first group multiple features arising from the same analyte, which we call "degenerate features", using a context-driven annotation approach. Surprisingly, this analysis revealed thousands of previously unreported degeneracies that reduced the number of unique analytes to ∼2961. We then applied an orthogonal approach to remove nonbiological features from the data using the 13C-based credentialing technology. This further reduced the number of unique analytes to less than 1000. Our 90% reduction in data is 5-fold greater than previously published studies. On the basis of the results, we propose an alternative approach to untargeted metabolomics that relies on thoroughly annotated reference data sets. To this end, we introduce the creDBle database ( http://creDBle.wustl.edu ), which contains accurate mass, retention time, and MS/MS fragmentation data as well as annotations of all credentialed features.

MetExtract II: A Software Suite for Stable Isotope-Assisted Untargeted Metabolomics

Stable isotope labeling (SIL) techniques have the potential to enhance different aspects of liquid chromatography–high-resolution mass spectrometry (LC-HRMS)-based untargeted metabolomics methods including metabolite detection, annotation of unknown metabolites, and comparative quantification. In this work, we present MetExtract II, a software toolbox for detection of biologically derived compounds. It exploits SIL-specific isotope patterns and elution profiles in LC-HRMS(/MS) data. The toolbox consists of three complementary modules: M1 (AllExtract) uses mixtures of uniformly highly isotope-enriched and native biological samples for selective detection of the entire accessible metabolome. M2 (TracExtract) is particularly suited to probe the metabolism of endogenous or exogenous secondary metabolites and facilitates the untargeted screening of tracer derivatives from concurrently metabolized native and uniformly labeled tracer substances. With M3 (FragExtract), tandem mass spectrometry (MS/MS) fragments of corresponding native and uniformly labeled ions are evaluated and automatically assigned with putative sum formulas. Generated results can be graphically illustrated and exported as a comprehensive data matrix that contains all detected pairs of native and labeled metabolite ions that can be used for database queries, metabolome-wide internal standardization, and statistical analysis. The software, associated documentation, and sample data sets are freely available for noncommercial use at http://metabolomics-ifa.boku.ac.at/metextractII.

Characterization of the Role of the Malate Dehydrogenases to Lung Tumor Cell Survival

Cellular compartmentalization of biochemical processes in eukaryotic cells is critical for many functions including shuttling of reducing equivalents across membranes. Although coordination of metabolic flux between different organelles is vital for cell physiology, its impact on tumor cell survival is not well understood. By using an integrative approach, we have dissected the role of the key metabolic enzymes Malate dehydrogenases (MDH1 and MDH2) to the survival of Non-small Cell Lung Carcinomas. Here, we report that while both the MDH1 (cytosolic) and the MDH2 (mitochondrial) enzymes display elevated levels in patients compared to normal counterparts, only high expression of MDH1 is associated with poor prognosis. We further show that the MDH1 enzymatic activity is significantly higher in NSCLC cells than that of MDH2. Accordingly, genetic depletion of MDH1 leads to significantly higher toxicity than depletion of MDH2. These findings provide molecular insights into the metabolic characteristics of the malate isoenzymes and mark MDH1 as a potential therapeutic target in these tumors.

From chromatogram to analyte to metabolite. How to pick horses for courses from the massive web resources for mass spectral plant metabolomics

The grand challenge currently facing metabolomics is the expansion of the coverage of the metabolome from a minor percentage of the metabolic complement of the cell toward the level of coverage afforded by other post-genomic technologies such as transcriptomics and proteomics. In plants, this problem is exacerbated by the sheer diversity of chemicals that constitute the metabolome, with the number of metabolites in the plant kingdom generally considered to be in excess of 200 000. In this review, we focus on web resources that can be exploited in order to improve analyte and ultimately metabolite identification and quantification. There is a wide range of available software that not only aids in this but also in the related area of peak alignment; however, for the uninitiated, choosing which program to use is a daunting task. For this reason, we provide an overview of the pros and cons of the software as well as comments regarding the level of programing skills required to effectively exploit their basic functions. In addition, the torrent of available genome and transcriptome sequences that followed the advent of next-generation sequencing has opened up further valuable resources for metabolite identification. All things considered, we posit that only via a continued communal sharing of information such as that deposited in the databases described within the article are we likely to be able to make significant headway toward improving our coverage of the plant metabolome.

Analysis of stable isotope assisted metabolomics data acquired by GC-MS

Stable isotope assisted metabolomics (SIAM) measures the abundance levels of metabolites in a particular pathway using stable isotope tracers (e.g., ¹³C, ¹⁸O and/or ¹⁵N). We report a method termed signature ion approach for analysis of SIAM data acquired on a GC-MS system equipped with an electron ionization (EI) ion source. The signature ion is a fragment ion in EI mass spectrum of a derivatized metabolite that contains all atoms of the underivatized metabolite, except the hydrogen atoms lost during derivatization. In this approach, GC-MS data of metabolite standards were used to recognize the signature ion from the EI mass spectra acquired from stable isotope labeled samples, and a linear regression model was used to deconvolute the intensity of overlapping isotopologues. A mixture score function was also employed for cross-sample chromatographic peak list alignment to recognize the chromatographic peaks generated by the same metabolite in different samples, by simultaneously evaluating the similarity of retention time and EI mass spectrum of two chromatographic peaks. Analysis of a mixture of 16 ¹³C-labeled and 16 unlabeled amino acids showed that the signature ion approach accurately identified and quantified all isotopologues. Analysis of polar metabolite extracts from cells respectively fed with uniform ¹³C-glucose and ¹³C-glutamine further demonstrated that this method can also be used to analyze the complex data acquired from biological samples.

Analysis of Stable Isotope Assisted Metabolomics Data Acquired by High Resolution Mass Spectrometry

Stable isotope assisted metabolomics (SIAM) uses stable isotope tracers to support studies of biochemical mechanisms. We report a suite of data analysis algorithms for automatic analysis of SIAM data acquired on a high resolution mass spectrometer. To increase the accuracy of isotopologue assignment, metabolites detected in the unlabeled samples were used as reference metabolites to generate possible isotopologue candidates for analysis of peaks detected in the labeled samples. An iterative linear regression model was developed to deconvolute the overlapping isotopic peaks of isotopologues present in a full MS spectrum, where the threshold for the weight factor was determined by a simulation study assuming different levels of Gaussian white noise contamination. A normalization method enabling isotope ratio-based normalization was implemented to study the difference of isotopologue abundance distribution between sample groups. The developed method can analyze SIAM data acquired by direct infusion MS and LC-MS, and can handle metabolite tracers containing different tracer elements. Analysis of SIAM data acquired from mixtures of known compounds showed that the developed algorithms accurately identify metabolites and quantify stable isotope enrichment. Application of SIAM data acquired from a biological study further demonstrated the effectiveness and accuracy of the developed method for analysis of complex samples.

Navigating freely-available software tools for metabolomics analysis

INTRODUCTION

The field of metabolomics has expanded greatly over the past two decades, both as an experimental science with applications in many areas, as well as in regards to data standards and bioinformatics software tools. The diversity of experimental designs and instrumental technologies used for metabolomics has led to the need for distinct data analysis methods and the development of many software tools.

OBJECTIVES

To compile a comprehensive list of the most widely used freely available software and tools that are used primarily in metabolomics.

METHODS

The most widely used tools were selected for inclusion in the review by either ≥ 50 citations on Web of Science (as of 08/09/16) or the use of the tool being reported in the recent Metabolomics Society survey. Tools were then categorised by the type of instrumental data (i.e. LC-MS, GC-MS or NMR) and the functionality (i.e. pre- and post-processing, statistical analysis, workflow and other functions) they are designed for.

RESULTS

A comprehensive list of the most used tools was compiled. Each tool is discussed within the context of its application domain and in relation to comparable tools of the same domain. An extended list including additional tools is available at https://github.com/RASpicer/MetabolomicsTools which is classified and searchable via a simple controlled vocabulary.

CONCLUSION

This review presents the most widely used tools for metabolomics analysis, categorised based on their main functionality. As future work, we suggest a direct comparison of tools' abilities to perform specific data analysis tasks e.g. peak picking.

Optimized Method for Untargeted Metabolomics Analysis of MDA-MB-231 Breast Cancer Cells

Cancer cells often have dysregulated metabolism, which is largely characterized by the Warburg effect-an increase in glycolytic activity at the expense of oxidative phosphorylation-and increased glutamine utilization. Modern metabolomics tools offer an efficient means to investigate metabolism in cancer cells. Currently, a number of protocols have been described for harvesting adherent cells for metabolomics analysis, but the techniques vary greatly and they lack specificity to particular cancer cell lines with diverse metabolic and structural features. Here we present an optimized method for untargeted metabolomics characterization of MDA-MB-231 triple negative breast cancer cells, which are commonly used to study metastatic breast cancer. We found that an approach that extracted all metabolites in a single step within the culture dish optimally detected both polar and non-polar metabolite classes with higher relative abundance than methods that involved removal of cells from the dish. We show that this method is highly suited to diverse applications, including the characterization of central metabolic flux by stable isotope labelling and differential analysis of cells subjected to specific pharmacological interventions.

MIA: non-targeted mass isotopolome analysis

MIA detects and visualizes isotopic enrichment in gas chromatography electron ionization mass spectrometry (GC-EI-MS) datasets in a non-targeted manner. It provides an easy-to-use graphical user interface that allows for visual mass isotopomer distribution analysis across multiple datasets. MIA helps to reveal changes in metabolic fluxes, visualizes metabolic proximity of isotopically enriched compounds and shows the fate of the applied stable isotope labeled tracer.

AVAILABILITY AND IMPLEMENTATION

Linux and Windows binaries, documentation, and sample data are freely available for download at http://massisotopolomeanalyzer.lu MIA is a stand-alone application implemented in C ++  and based on Qt5, NTFD and the MetaboliteDetector framework.

CONTACT

karsten.hiller@uni.lu.

Stable Isotope-Assisted Evaluation of Different Extraction Solvents for Untargeted Metabolomics of Plants

The evaluation of extraction protocols for untargeted metabolomics approaches is still difficult. We have applied a novel stable isotope-assisted workflow for untargeted LC-HRMS-based plant metabolomics , which allows for the first time every detected feature to be considered for method evaluation. The efficiency and complementarity of commonly used extraction solvents, namely 1 + 3 (v/v) mixtures of water and selected organic solvents (methanol, acetonitrile or methanol/acetonitrile 1 + 1 (v/v)), with and without the addition of 0.1% (v/v) formic acid were compared. Four different wheat organs were sampled, extracted and analysed by LC-HRMS. Data evaluation was performed with the in-house-developed MetExtract II software and R. With all tested solvents a total of 871 metabolites were extracted in ear, 785 in stem, 733 in leaf and 517 in root samples, respectively. Between 48% (stem) and 57% (ear) of the metabolites detected in a particular organ were found with all extraction mixtures, and 127 of 996 metabolites were consistently shared between all extraction agent/organ combinations. In aqueous methanol, acidification with formic acid led to pronounced pH dependency regarding the precision of metabolite abundance and the number of detectable metabolites, whereas extracts of acetonitrile-containing mixtures were less affected. Moreover, methanol and acetonitrile have been found to be complementary with respect to extraction efficiency. Interestingly, the beneficial properties of both solvents can be combined by the use of a water-methanol-acetonitrile mixture for global metabolite extraction instead of aqueous methanol or aqueous acetonitrile alone.

Unsupervised Identification of Isotope-Labeled Peptides

In vivo isotopic labeling coupled with high-resolution proteomics is used to investigate primary metabolism in techniques such as stable isotope probing (protein-SIP) and peptide-based metabolic flux analysis (PMFA). Isotopic enrichment of carbon substrates and intracellular metabolism determine the distribution of isotopes within amino acids. The resulting amino acid mass distributions (AMDs) are convoluted into peptide mass distributions (PMDs) during protein synthesis. With no a priori knowledge on metabolic fluxes, the PMDs are therefore unknown. This complicates labeled peptide identification because prior knowledge on PMDs is used in all available peptide identification software. An automated framework for the identification and quantification of PMDs for nonuniformly labeled samples is therefore lacking. To unlock the potential of peptide labeling experiments for high-throughput flux analysis and other complex labeling experiments, an unsupervised peptide identification and quantification method was developed that uses discrete deconvolution of mass distributions of identified peptides to inform on the mass distributions of otherwise unidentifiable peptides. Uniformly (13)C-labeled Escherichia coli protein was used to test the developed feature reconstruction and deconvolution algorithms. The peptide identification was validated by comparing MS(2)-identified peptides to peptides identified from PMDs using unlabeled E. coli protein. Nonuniformly labeled Glycine max protein was used to demonstrate the technology on a representative sample suitable for flux analysis. Overall, automatic peptide identification and quantification were comparable or superior to manual extraction, enabling proteomics-based technology for high-throughput flux analysis studies.

Bridging the gap between non-targeted stable isotope labeling and metabolic flux analysis

BACKGROUND

Metabolism gained increasing interest for the understanding of diseases and to pinpoint therapeutic intervention points. However, classical metabolomics techniques only provide a very static view on metabolism. Metabolic flux analysis methods, on the other hand, are highly targeted and require detailed knowledge on metabolism beforehand.

RESULTS

We present a novel workflow to analyze non-targeted metabolome-wide stable isotope labeling data to detect metabolic flux changes in a non-targeted manner. Furthermore, we show how similarity-analysis of isotopic enrichment patterns can be used for pathway contextualization of unidentified compounds. We illustrate our approach with the analysis of changes in cellular metabolism of human adenocarcinoma cells in response to decreased oxygen availability. Starting without a priori knowledge, we detect metabolic flux changes, leading to an increased glutamine contribution to acetyl-CoA production, reveal biosynthesis of N-acetylaspartate by N-acetyltransferase 8-like (NAT8L) in lung cancer cells and show that NAT8L silencing inhibits proliferation of A549, JHH-4, PH5CH8, and BEAS-2B cells.

CONCLUSIONS

Differential stable isotope labeling analysis provides qualitative metabolic flux information in a non-targeted manner. Furthermore, similarity analysis of enrichment patterns provides information on metabolically closely related compounds. N-acetylaspartate and NAT8L are important players in cancer cell metabolism, a context in which they have not received much attention yet.

geoRge: A Computational Tool To Detect the Presence of Stable Isotope Labeling in LC/MS-Based Untargeted Metabolomics

Studying the flow of chemical moieties through the complex set of metabolic reactions that happen in the cell is essential to understanding the alterations in homeostasis that occur in disease. Recently, LC/MS-based untargeted metabolomics and isotopically labeled metabolites have been used to facilitate the unbiased mapping of labeled moieties through metabolic pathways. However, due to the complexity of the resulting experimental data sets few computational tools are available for data analysis. Here we introduce geoRge, a novel computational approach capable of analyzing untargeted LC/MS data from stable isotope-labeling experiments. geoRge is written in the open language R and runs on the output structure of the XCMS package, which is in widespread use. As opposed to the few existing tools, which use labeled samples to track stable isotopes by iterating over all MS signals using the theoretical mass difference between the light and heavy isotopes, geoRge uses unlabeled and labeled biologically equivalent samples to compare isotopic distributions in the mass spectra. Isotopically enriched compounds change their isotopic distribution as compared to unlabeled compounds. This is directly reflected in a number of new m/z peaks and higher intensity peaks in the mass spectra of labeled samples relative to the unlabeled equivalents. The automated untargeted isotope annotation and relative quantification capabilities of geoRge are demonstrated by the analysis of LC/MS data from a human retinal pigment epithelium cell line (ARPE-19) grown on normal and high glucose concentrations mimicking diabetic retinopathy conditions in vitro. In addition, we compared the results of geoRge with the outcome of X(13)CMS, since both approaches rely entirely on XCMS parameters for feature selection, namely m/z and retention time values. To ensure data traceability and reproducibility, and enabling for comparison with other existing and future approaches, raw LC/MS files have been deposited in MetaboLights (MTBLS213) and geoRge is available as an R script at https://github.com/jcapelladesto/geoRge.

Non-targeted Tracer Fate Detection

Stable isotopes have been used to trace atoms through metabolism and quantify metabolic fluxes for several decades. Only recently non-targeted stable isotope labeling approaches have emerged as a powerful tool to gain insights into metabolism. However, the manual detection of isotopic enrichment for a non-targeted analysis is tedious and time consuming. To overcome this limitation, the non-targeted tracer fate detection (NTFD) algorithm for the automated metabolome-wide detection of isotopic enrichment has been developed. NTFD detects and quantifies isotopic enrichment in the form of mass isotopomer distributions (MIDs) in an automated manner, providing the means to trace functional groups, determine MIDs for metabolic flux analysis, or detect tracer-derived molecules in general. Here, we describe the algorithmic background of NTFD, discuss practical considerations for the freely available NTFD software package, and present potential applications of non-targeted stable isotope labeling analysis.

Probing the metabolic network in bloodstream-form Trypanosoma brucei using untargeted metabolomics with stable isotope labelled glucose

Metabolomics coupled with heavy-atom isotope-labelled glucose has been used to probe the metabolic pathways active in cultured bloodstream form trypomastigotes of Trypanosoma brucei, a parasite responsible for human African trypanosomiasis. Glucose enters many branches of metabolism beyond glycolysis, which has been widely held to be the sole route of glucose metabolism. Whilst pyruvate is the major end-product of glucose catabolism, its transamination product, alanine, is also produced in significant quantities. The oxidative branch of the pentose phosphate pathway is operative, although the non-oxidative branch is not. Ribose 5-phosphate generated through this pathway distributes widely into nucleotide synthesis and other branches of metabolism. Acetate, derived from glucose, is found associated with a range of acetylated amino acids and, to a lesser extent, fatty acids; while labelled glycerol is found in many glycerophospholipids. Glucose also enters inositol and several sugar nucleotides that serve as precursors to macromolecule biosynthesis. Although a Krebs cycle is not operative, malate, fumarate and succinate, primarily labelled in three carbons, were present, indicating an origin from phosphoenolpyruvate via oxaloacetate. Interestingly, the enzyme responsible for conversion of phosphoenolpyruvate to oxaloacetate, phosphoenolpyruvate carboxykinase, was shown to be essential to the bloodstream form trypanosomes, as demonstrated by the lethal phenotype induced by RNAi-mediated downregulation of its expression. In addition, glucose derivatives enter pyrimidine biosynthesis via oxaloacetate as a precursor to aspartate and orotate.

In this work we have followed the distribution of carbon derived from glucose in bloodstream form trypanosomes, the causative agent of African trypanosomiasis, revealing it to enter a diverse range of metabolites. The work involved using ¹³C-labelled glucose and following the fate of the labelled carbon with an LC-MS based metabolomics platform. Beyond glycolysis and the oxidative branch of the pentose phosphate pathway the label entered lipid biosynthesis both through glycerol 3-phosphate and also acetate. Glucose derived carbon also entered nucleotide synthesis through ribose and pyrimidine synthesis through oxaloacetate-derived aspartate. Appreciable quantities of the carboxylic acids succinate and malate were identified, although labelling patterns indicate they are not TCA cycle derived. Amino sugars and sugar nucleotides were also labelled as was inositol used in protein modification but not in inositol phospholipid headgroup production. We confirm active and essential oxaloacetate production in bloodstream form trypanosomes and show that phosphoenolpyruvate carboxykinase is essential to these parasites using RNA interference. The amount of glucose entering these metabolites is minor compared to the quantity that enters pyruvate excreted from the cell, but the observation that enzymes contributing to the metabolism of glucose beyond glycolysis can be essential offers potential new targets for chemotherapy against trypanosomiasis.

LC-MS-based absolute metabolite quantification: application to metabolic flux measurement in trypanosomes

Human African trypanosomiasis is a neglected tropical disease caused by the protozoan parasite, Trypanosoma brucei. In the mammalian bloodstream, the trypanosome's metabolism differs significantly from that of its host. For example, the parasite relies exclusively on glycolysis for energy source. Recently, computational and mathematical models of trypanosome metabolism have been generated to assist in understanding the parasite metabolism with the aim of facilitating drug development. Optimisation of these models requires quantitative information, including metabolite concentrations and/or metabolic fluxes that have been hitherto unavailable on a large scale. Here, we have implemented an LC-MS-based method that allows large scale quantification of metabolite levels by using U-¹³C-labelled E.coli extracts as internal standards. Known amounts of labelled E. coli extract were added into the parasite samples, as well as calibration standards, and used to obtain calibration curves enabling us to convert intensities into concentrations. This method allowed us to reliably quantify the changes of 43 intracellular metabolites and 32 extracellular metabolites in the medium over time. Based on the absolute quantification, we were able to compute consumption and production fluxes. These quantitative data can now be used to optimise computational models of parasite metabolism.

Integrating proteomics profiling data sets: a network perspective

Understanding disease mechanisms often requires complex and accurate integration of cellular pathways and molecular networks. Systems biology offers the possibility to provide a comprehensive map of the cell's intricate wiring network, which can ultimately lead to decipher the disease phenotype. Here, we describe what biological pathways are, how they function in normal and abnormal cellular systems, limitations faced by databases for integrating data, and highlight how network models are emerging as a powerful integrative framework to understand and interpret the roles of proteins and peptides in diseases.

Untargeted profiling of tracer-derived metabolites using stable isotopic labeling and fast polarity-switching LC-ESI-HRMS

An untargeted metabolomics workflow for the detection of metabolites derived from endogenous or exogenous tracer substances is presented. To this end, a recently developed stable isotope-assisted LC-HRMS-based metabolomics workflow for the global annotation of biological samples has been further developed and extended. For untargeted detection of metabolites arising from labeled tracer substances, isotope pattern recognition has been adjusted to account for nonlabeled moieties conjugated to the native and labeled tracer molecules. Furthermore, the workflow has been extended by (i) an optional ion intensity ratio check, (ii) the automated combination of positive and negative ionization mode mass spectra derived from fast polarity switching, and (iii) metabolic feature annotation. These extensions enable the automated, unbiased, and global detection of tracer-derived metabolites in complex biological samples. The workflow is demonstrated with the metabolism of (13)C9-phenylalanine in wheat cell suspension cultures in the presence of the mycotoxin deoxynivalenol (DON). In total, 341 metabolic features (150 in positive and 191 in negative ionization mode) corresponding to 139 metabolites were detected. The benefit of fast polarity switching was evident, with 32 and 58 of these metabolites having exclusively been detected in the positive and negative modes, respectively. Moreover, for 19 of the remaining 49 phenylalanine-derived metabolites, the assignment of ion species and, thus, molecular weight was possible only by the use of complementary features of the two ion polarity modes. Statistical evaluation showed that treatment with DON increased or decreased the abundances of many detected metabolites.

ALLocator: an interactive web platform for the analysis of metabolomic LC-ESI-MS datasets, enabling semi-automated, user-revised compound annotation and mass isotopomer ratio analysis

Adduct formation, fragmentation events and matrix effects impose special challenges to the identification and quantitation of metabolites in LC-ESI-MS datasets. An important step in compound identification is the deconvolution of mass signals. During this processing step, peaks representing adducts, fragments, and isotopologues of the same analyte are allocated to a distinct group, in order to separate peaks from coeluting compounds. From these peak groups, neutral masses and pseudo spectra are derived and used for metabolite identification via mass decomposition and database matching. Quantitation of metabolites is hampered by matrix effects and nonlinear responses in LC-ESI-MS measurements. A common approach to correct for these effects is the addition of a U-¹³C-labeled internal standard and the calculation of mass isotopomer ratios for each metabolite. Here we present a new web-platform for the analysis of LC-ESI-MS experiments. ALLocator covers the workflow from raw data processing to metabolite identification and mass isotopomer ratio analysis. The integrated processing pipeline for spectra deconvolution “ALLocatorSD” generates pseudo spectra and automatically identifies peaks emerging from the U-¹³C-labeled internal standard. Information from the latter improves mass decomposition and annotation of neutral losses. ALLocator provides an interactive and dynamic interface to explore and enhance the results in depth. Pseudo spectra of identified metabolites can be stored in user- and method-specific reference lists that can be applied on succeeding datasets. The potential of the software is exemplified in an experiment, in which abundance fold-changes of metabolites of the l-arginine biosynthesis in C. glutamicum type strain ATCC 13032 and l-arginine producing strain ATCC 21831 are compared. Furthermore, the capability for detection and annotation of uncommon large neutral losses is shown by the identification of (γ-)glutamyl dipeptides in the same strains. ALLocator is available online at: https://allocator.cebitec.uni-bielefeld.de. A login is required, but freely available.

Stable isotope-labeling studies in metabolomics: new insights into structure and dynamics of metabolic networks

The rapid emergence of metabolomics has enabled system-wide measurements of metabolites in various organisms. However, advances in the mechanistic understanding of metabolic networks remain limited, as most metabolomics studies cannot routinely provide accurate metabolite identification, absolute quantification and flux measurement. Stable isotope labeling offers opportunities to overcome these limitations. Here we describe some current approaches to stable isotope-labeled metabolomics and provide examples of the significant impact that these studies have had on our understanding of cellular metabolism. Furthermore, we discuss recently developed software solutions for the analysis of stable isotope-labeled metabolomics data and propose the bioinformatics solutions that will pave the way for the broader application and optimal interpretation of system-scale labeling studies in metabolomics.

Credentialing features: a platform to benchmark and optimize untargeted metabolomic methods

The aim of untargeted metabolomics is to profile as many metabolites as possible, yet a major challenge is comparing experimental method performance on the basis of metabolome coverage. To date, most published approaches have compared experimental methods by counting the total number of features detected. Due to artifactual interference, however, this number is highly variable and therefore is a poor metric for comparing metabolomic methods. Here we introduce an alternative approach to benchmarking metabolome coverage which relies on mixed Escherichia coli extracts from cells cultured in regular and ¹³C-enriched media. After mass spectrometry-based metabolomic analysis of these extracts, we “credential” features arising from E. coli metabolites on the basis of isotope spacing and intensity. This credentialing platform enables us to accurately compare the number of nonartifactual features yielded by different experimental approaches. We highlight the value of our platform by reoptimizing a published untargeted metabolomic method for XCMS data processing. Compared to the published parameters, the new XCMS parameters decrease the total number of features by 15% (a reduction in noise features) while increasing the number of true metabolites detected and grouped by 20%. Our credentialing platform relies on easily generated E. coli samples and a simple software algorithm that is freely available on our laboratory Web site (http://pattilab.wustl.edu/software/credential/). We have validated the credentialing platform with reversed-phase and hydrophilic interaction liquid chromatography as well as Agilent, Thermo Scientific, AB SCIEX, and LECO mass spectrometers. Thus, the credentialing platform can readily be applied by any laboratory to optimize their untargeted metabolomic pipeline for metabolite extraction, chromatographic separation, mass spectrometric detection, and bioinformatic processing.

Novel strategy for non-targeted isotope-assisted metabolomics by means of metabolic turnover and multivariate analysis

Isotope-labeling is a useful technique for understanding cellular metabolism. Recent advances in metabolomics have extended the capability of isotope-assisted studies to reveal global metabolism. For instance, isotope-assisted metabolomics technology has enabled the mapping of a global metabolic network, estimation of flux at branch points of metabolic pathways, and assignment of elemental formulas to unknown metabolites. Furthermore, some data processing tools have been developed to apply these techniques to a non-targeted approach, which plays an important role in revealing unknown or unexpected metabolism. However, data collection and integration strategies for non-targeted isotope-assisted metabolomics have not been established. Therefore, a systematic approach is proposed to elucidate metabolic dynamics without targeting pathways by means of time-resolved isotope tracking, i.e., "metabolic turnover analysis", as well as multivariate analysis. We applied this approach to study the metabolic dynamics in amino acid perturbation of Saccharomyces cerevisiae. In metabolic turnover analysis, 69 peaks including 35 unidentified peaks were investigated. Multivariate analysis of metabolic turnover successfully detected a pathway known to be inhibited by amino acid perturbation. In addition, our strategy enabled identification of unknown peaks putatively related to the perturbation.

Phosphoenolpyruvate carboxylase identified as a key enzyme in erythrocytic Plasmodium falciparum carbon metabolism

Phospoenolpyruvate carboxylase (PEPC) is absent from humans but encoded in the Plasmodium falciparum genome, suggesting that PEPC has a parasite-specific function. To investigate its importance in P. falciparum, we generated a pepc null mutant (D10^Δpepc), which was only achievable when malate, a reduction product of oxaloacetate, was added to the growth medium. D10^Δpepc had a severe growth defect in vitro, which was partially reversed by addition of malate or fumarate, suggesting that pepc may be essential in vivo. Targeted metabolomics using ¹³C-U-D-glucose and ¹³C-bicarbonate showed that the conversion of glycolytically-derived PEP into malate, fumarate, aspartate and citrate was abolished in D10^Δpepc and that pentose phosphate pathway metabolites and glycerol 3-phosphate were present at increased levels. In contrast, metabolism of the carbon skeleton of ¹³C,¹⁵N-U-glutamine was similar in both parasite lines, although the flux was lower in D10^Δpepc; it also confirmed the operation of a complete forward TCA cycle in the wild type parasite. Overall, these data confirm the CO₂ fixing activity of PEPC and suggest that it provides metabolites essential for TCA cycle anaplerosis and the maintenance of cytosolic and mitochondrial redox balance. Moreover, these findings imply that PEPC may be an exploitable target for future drug discovery.

The genome of the human malaria parasite Plasmodium falciparum encodes a protein called phosphoenolpyruvate carboxylase (PEPC) absent from the human host. PEPC is known to fix CO₂ to generate metabolites used for energy metabolism in plants and bacteria, but its function in malaria parasites remained an enigma. Our study aimed to elucidate the role and importance of PEPC in P. falciparum in its host red blood cell by generating a gene deletion mutant in P. falciparum. This was only achievable in the presence of high concentrations of malate were added to the culture medium. The mutant generated (D10^Δpepc) had a severe growth defect, which was rescued partially by malate or fumarate (but not any other downstream metabolites), suggesting that they feed into the same metabolic pathway. Using heavy isotope labelled ¹³C-U-D-glucose and ¹³C-bicarbonate we showed that PECP has an important role in intermediary carbon metabolism and is vital for the maintenance of cytosolic and mitochondrial redox balance. Together these findings imply that PEPC may be an exploitable target for future drug discovery.

A novel stable isotope labelling assisted workflow for improved untargeted LC-HRMS based metabolomics research

Many untargeted LC-ESI-HRMS based metabolomics studies are still hampered by the large proportion of non-biological sample derived signals included in the generated raw data. Here, a novel, powerful stable isotope labelling (SIL)-based metabolomics workflow is presented, which facilitates global metabolome extraction, improved metabolite annotation and metabolome wide internal standardisation (IS). The general concept is exemplified with two different cultivation variants, (1) co-cultivation of the plant pathogenic fungi Fusarium graminearum on non-labelled and highly ¹³C enriched culture medium and (2) experimental cultivation under native conditions and use of globally U-¹³C labelled biological reference samples as exemplified with maize and wheat. Subsequent to LC-HRMS analysis of mixtures of labelled and non-labelled samples, two-dimensional data filtering of SIL specific isotopic patterns is performed to better extract truly biological derived signals together with the corresponding number of carbon atoms of each metabolite ion. Finally, feature pairs are convoluted to feature groups each representing a single metabolite. Moreover, the correction of unequal matrix effects in different sample types and the improvement of relative metabolite quantification with metabolome wide IS are demonstrated for the F. graminearum experiment. Data processing employing the presented workflow revealed about 300 SIL derived feature pairs corresponding to 87-135 metabolites in F. graminearum samples and around 800 feature pairs corresponding to roughly 350 metabolites in wheat samples. SIL assisted IS, by the use of globally U-¹³C labelled biological samples, reduced the median CV value from 7.1 to 3.6 % for technical replicates and from 15.1 to 10.8 % for biological replicates in the respective F. graminearum samples.