Identification of a Metabolic Reaction Network from Time-Series Data of Metabolite Concentrations

Recent development of high-throughput analytical techniques has made it possible to qualitatively identify a number of metabolites simultaneously. Correlation and multivariate analyses such as principal component analysis have been widely used to analyse those data and evaluate correlations among the metabolic profiles. However, these analyses cannot simultaneously carry out identification of metabolic reaction networks and prediction of dynamic behaviour of metabolites in the networks. The present study, therefore, proposes a new approach consisting of a combination of statistical technique and mathematical modelling approach to identify and predict a probable metabolic reaction network from time-series data of metabolite concentrations and simultaneously construct its mathematical model. Firstly, regression functions are fitted to experimental data by the locally estimated scatter plot smoothing method. Secondly, the fitted result is analysed by the bivariate Granger causality test to determine which metabolites cause the change in other metabolite concentrations and remove less related metabolites. Thirdly, S-system equations are formed by using the remaining metabolites within the framework of biochemical systems theory. Finally, parameters including rate constants and kinetic orders are estimated by the Levenberg–Marquardt algorithm. The estimation is iterated by setting insignificant kinetic orders at zero, i.e., removing insignificant metabolites. Consequently, a reaction network structure is identified and its mathematical model is obtained. Our approach is validated using a generic inhibition and activation model and its practical application is tested using a simplified model of the glycolysis of Lactococcus lactis MG1363, for which actual time-series data of metabolite concentrations are available. The results indicate the usefulness of our approach and suggest a probable pathway for the production of lactate and acetate. The results also indicate that the approach pinpoints a probable strong inhibition of lactate on the glycolysis pathway.     

Structural analysis of metabolic networks based on flux centrality

Metabolic reactions are fundamental to living organisms, and a large number of reactions simultaneously occur at a given time in living cells transforming diverse metabolites into each other. There has been an ongoing debate on how to classify metabolites with respect to their importance for metabolic performance, usually based on the analysis of topological properties of genome scale metabolic networks. However, none of these studies have accounted quantitatively for flux in metabolic networks, thus lacking an important component of a cell’s biochemistry.We therefore analyzed a genome scale metabolic network of Escherichia coli by comparing growth under 19 different growth conditions, using flux balance analysis and weighted network centrality investigation. With this novel concept of flux centrality we generated metabolite rankings for each particular growth condition. In contrast to the results of conventional analysis of genome scale metabolic networks, different metabolites were top-ranking dependent on the growth condition. At the same time, several metabolites were consistently among the high ranking ones. Those are associated with pathways that have been described by biochemists as the most central part of metabolism, such as glycolysis, tricarboxylic acid cycle and pentose phosphate pathway. The values for the average path length of the analyzed metabolite networks were between 10.5 and 12.6, supporting recent findings that the metabolic network of E. coli is not a small-world network.     

Flux Imbalance Analysis and the Sensitivity of Cellular Growth to Changes in Metabolite Pools

Stoichiometric models of metabolism, such as flux balance analysis (FBA), are classically applied to predicting steady state rates - or fluxes - of metabolic reactions in genome-scale metabolic networks. Here we revisit the central assumption of FBA, i.e. that intracellular metabolites are at steady state, and show that deviations from flux balance (i.e. flux imbalances) are informative of some features of in vivo metabolite concentrations. Mathematically, the sensitivity of FBA to these flux imbalances is captured by a native feature of linear optimization, the dual problem, and its corresponding variables, known as shadow prices. First, using recently published data on chemostat growth of Saccharomyces cerevisae under different nutrient limitations, we show that shadow prices anticorrelate with experimentally measured degrees of growth limitation of intracellular metabolites. We next hypothesize that metabolites which are limiting for growth (and thus have very negative shadow price) cannot vary dramatically in an uncontrolled way, and must respond rapidly to perturbations. Using a collection of published datasets monitoring the time-dependent metabolomic response of Escherichia coli to carbon and nitrogen perturbations, we test this hypothesis and find that metabolites with negative shadow price indeed show lower temporal variation following a perturbation than metabolites with zero shadow price. Finally, we illustrate the broader applicability of flux imbalance analysis to other constraint-based methods. In particular, we explore the biological significance of shadow prices in a constraint-based method for integrating gene expression data with a stoichiometric model. In this case, shadow prices point to metabolites that should rise or drop in concentration in order to increase consistency between flux predictions and gene expression data. In general, these results suggest that the sensitivity of metabolic optima to violations of the steady state constraints carries biologically significant information on the processes that control intracellular metabolites in the cell.     

Stability of Metabolic Correlations under Changing Environmental Conditions in Escherichia coli – A Systems Approach

Background

Biological systems adapt to changing environments by reorganizing their cellular and physiological program with metabolites representing one important response level. Different stresses lead to both conserved and specific responses on the metabolite level which should be reflected in the underlying metabolic network.

Methodology/Principal Findings

Starting from experimental data obtained by a GC-MS based high-throughput metabolic profiling technology we here develop an approach that: (1) extracts network representations from metabolic condition-dependent data by using pairwise correlations, (2) determines the sets of stable and condition-dependent correlations based on a combination of statistical significance and homogeneity tests, and (3) can identify metabolites related to the stress response, which goes beyond simple observations about the changes of metabolic concentrations. The approach was tested with Escherichia coli as a model organism observed under four different environmental stress conditions (cold stress, heat stress, oxidative stress, lactose diauxie) and control unperturbed conditions. By constructing the stable network component, which displays a scale free topology and small-world characteristics, we demonstrated that: (1) metabolite hubs in this reconstructed correlation networks are significantly enriched for those contained in biochemical networks such as EcoCyc, (2) particular components of the stable network are enriched for functionally related biochemical pathways, and (3) independently of the response scale, based on their importance in the reorganization of the correlation network a set of metabolites can be identified which represent hypothetical candidates for adjusting to a stress-specific response.

Conclusions/Significance

Network-based tools allowed the identification of stress-dependent and general metabolic correlation networks. This correlation-network-based approach does not rely on major changes in concentration to identify metabolites important for stress adaptation, but rather on the changes in network properties with respect to metabolites. This should represent a useful complementary technique in addition to more classical approaches.     

Genetic Basis of Metabolome Variation in Yeast

Metabolism, the conversion of nutrients into usable energy and biochemical building blocks, is an essential feature of all cells. The genetic factors responsible for inter-individual metabolic variability remain poorly understood. To investigate genetic causes of metabolome variation, we measured the concentrations of 74 metabolites across 100 segregants from a Saccharomyces cerevisiae cross by liquid chromatography-tandem mass spectrometry. We found 52 quantitative trait loci for 34 metabolites. These included linkages due to overt changes in metabolic genes, e.g., linking pyrimidine intermediates to the deletion of ura3. They also included linkages not directly related to metabolic enzymes, such as those for five central carbon metabolites to ira2, a Ras/PKA pathway regulator, and for the metabolites, S-adenosyl-methionine and S-adenosyl-homocysteine to slt2, a MAP kinase involved in cell wall integrity. The variant of ira2 that elevates metabolite levels also increases glucose uptake and ethanol secretion. These results highlight specific examples of genetic variability, including in genes without prior known metabolic regulatory function, that impact yeast metabolism.     

K-FIT: An accelerated kinetic parameterization algorithm using steady-state fluxomic data

Kinetic models predict the metabolic flows by directly linking metabolite concentrations and enzyme levels to reaction fluxes. Robust parameterization of organism-level kinetic models that faithfully reproduce the effect of different genetic or environmental perturbations remains an open challenge due to the intractability of existing algorithms. This paper introduces Kinetics-based Fluxomics Integration Tool (K-FIT), a robust kinetic parameterization workflow that leverages a novel decomposition approach to identify steady-state fluxes in response to genetic perturbations followed by a gradient-based update of kinetic parameters until predictions simultaneously agree with the fluxomic data in all perturbed metabolic networks. The applicability of K-FIT to large-scale models is demonstrated by parameterizing an expanded kinetic model for E. coli (307 reactions and 258 metabolites) using fluxomic data from six mutants. The achieved thousand-fold speed-up afforded by K-FIT over meta-heuristic approaches is transformational enabling follow-up robustness of inference analyses and optimal design of experiments to inform metabolic engineering strategies.     

Quantitative Assessment of Thermodynamic Constraints on the Solution Space of Genome-Scale Metabolic Models

Constraint-based methods provide powerful computational techniques to allow understanding and prediction of cellular behavior. These methods rely on physiochemical constraints to eliminate infeasible behaviors from the space of available behaviors. One such constraint is thermodynamic feasibility, the requirement that intracellular flux distributions obey the laws of thermodynamics. The past decade has seen several constraint-based methods that interpret this constraint in different ways, including those that are limited to small networks, rely on predefined reaction directions, and/or neglect the relationship between reaction free energies and metabolite concentrations. In this work, we utilize one such approach, thermodynamics-based metabolic flux analysis (TMFA), to make genome-scale, quantitative predictions about metabolite concentrations and reaction free energies in the absence of prior knowledge of reaction directions, while accounting for uncertainties in thermodynamic estimates. We applied TMFA to a genome-scale network reconstruction of Escherichia coli and examined the effect of thermodynamic constraints on the flux space. We also assessed the predictive performance of TMFA against gene essentiality and quantitative metabolomics data, under both aerobic and anaerobic, and optimal and suboptimal growth conditions. Based on these results, we propose that TMFA is a useful tool for validating phenotypes and generating hypotheses, and that additional types of data and constraints can improve predictions of metabolite concentrations.     

Impact of a Transposon Insertion in phzF2 on the Specialized Metabolite Production and Interkingdom Interactions of Pseudomonas aeruginosa

In microbiology, gene disruption and subsequent experiments often center on phenotypic changes caused by one class of specialized metabolites (quorum sensors, virulence factors, or natural products), disregarding global downstream metabolic effects. With the recent development of mass spectrometry-based methods and technologies for microbial metabolomics investigations, it is now possible to visualize global production of diverse classes of microbial specialized metabolites simultaneously. Using imaging mass spectrometry (IMS) applied to the analysis of microbiology experiments, we can observe the effects of mutations, knockouts, insertions, and complementation on the interactive metabolome. In this study, a combination of IMS and liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to visualize the impact on specialized metabolite production of a transposon insertion into a Pseudomonas aeruginosa phenazine biosynthetic gene, phzF2. The disruption of phenazine biosynthesis led to broad changes in specialized metabolite production, including loss of pyoverdine production. This shift in specialized metabolite production significantly alters the metabolic outcome of an interaction with Aspergillus fumigatus by influencing triacetylfusarinine production.     

Integration of Genome-Wide SNP Data and Gene-Expression Profiles Reveals Six Novel Loci and Regulatory Mechanisms for Amino Acids and Acylcarnitines in Whole Blood

Profiling amino acids and acylcarnitines in whole blood spots is a powerful tool in the laboratory diagnosis of several inborn errors of metabolism. Emerging data suggests that altered blood levels of amino acids and acylcarnitines are also associated with common metabolic diseases in adults. Thus, the identification of common genetic determinants for blood metabolites might shed light on pathways contributing to human physiology and common diseases. We applied a targeted mass-spectrometry-based method to analyze whole blood concentrations of 96 amino acids, acylcarnitines and pathway associated metabolite ratios in a Central European cohort of 2,107 adults and performed genome-wide association (GWA) to identify genetic modifiers of metabolite concentrations. We discovered and replicated six novel loci associated with blood levels of total acylcarnitine, arginine (both on chromosome 6; rs12210538, rs17657775), propionylcarnitine (chromosome 10; rs12779637), 2-hydroxyisovalerylcarnitine (chromosome 21; rs1571700), stearoylcarnitine (chromosome 1; rs3811444), and aspartic acid traits (chromosome 8; rs750472). Based on an integrative analysis of expression quantitative trait loci in blood mononuclear cells and correlations between gene expressions and metabolite levels, we provide evidence for putative causative genes: SLC22A16 for total acylcarnitines, ARG1 for arginine, HLCS for 2-hydroxyisovalerylcarnitine, JAM3 for stearoylcarnitine via a trans-effect at chromosome 1, and PPP1R16A for aspartic acid traits. Further, we report replication and provide additional functional evidence for ten loci that have previously been published for metabolites measured in plasma, serum or urine.In conclusion, our integrative analysis of SNP, gene-expression and metabolite data points to novel genetic factors that may be involved in the regulation of human metabolism. At several loci, we provide evidence for metabolite regulation via gene-expression and observed overlaps with GWAS loci for common diseases. These results form a strong rationale for subsequent functional and disease-related studies.     

Organic acid profile of isovaleric acidemia: a comprehensive metabolomics approach

Isovaleric acidemia (IVA, MIM 248600) can be a severe and potentially life-threatening disease in affected neonates, but with a positive prognosis on treatment for some phenotypes. This study presents the first application of metabolomics to evaluate the metabolite profiles derived from urine samples of untreated and treated IVA patients as well as of obligate heterozygotes. All IVA patients carried the same homozygous c.367 G > A nucleotide change in exon 4 of the IVD gene but manifested phenotypic diversity. Concurrent class analysis (CONCA) was used to compare all the metabolites from the original complete data set obtained from the three case and two control groups used in this investigation. This application of CONCA has not been reported previously, and is used here to compare four different modes of scaling of all metabolites. The variables important in discrimination from the CONCA thus enabled the recognition of different metabolic patterns encapsulated within the data sets that would not have been revealed by using only one mode of scaling. Application of multivariate and univariate analyses disclosed 11 important metabolites that distinguished untreated IVA from controls. These included well-established diagnostic biomarkers of IVA, endogenous detoxification markers, and 3-hydroxycaproic acid, an indicator of ketosis, but not reported previously for this disease. Nine metabolites were identified that reflected the effect of treatment of IVA. They included detoxification products and indicators related to the high carbohydrate and low protein diet which formed the hallmark of the treatment. This investigation also provides the first comparative metabolite profile for heterozygotes of this inherited metabolic disorder. The detection of informative metabolites in even very low concentrations in all three experimental groups highlights the potential advantage of the holistic mode of analysis of inherited metabolic diseases in a metabolomics investigation.     

Metabolite-centric approaches for the discovery of antibacterials using genome-scale metabolic networks

Development of genome-scale metabolic models and various constraints-based flux analyses have enabled more sophisticated examination of metabolism. Recently reported metabolite essentiality studies are also based on the constraints-based modeling, but approaches metabolism from a metabolite-centric perspective, providing synthetic lethal combination of reactions and clues for the rational discovery of antibacterials. In this study, metabolite essentiality analysis was applied to the genome-scale metabolic models of four microorganisms: Escherichia coli, Helicobacter pylori, Mycobacterium tuberculosis and Staphylococcus aureus. Furthermore, chokepoints, metabolites surrounded by enzymes that uniquely consume and/or produce them, were also calculated based on the network properties of the above organisms. A systematic drug targeting strategy was developed by combining information from these two methods. Final drug target metabolites are presented and examined with knowledge from the literature.     

Graph-Based Analysis of the Metabolic Exchanges between Two Co-Resident Intracellular Symbionts,Baumannia cicadellinicola and Sulcia muelleri,with Their Insect Host,Homalodisca coagulata

Endosymbiotic bacteria from different species can live inside cells of the same eukaryotic organism. Metabolic exchanges occur between host and bacteria but also between different endocytobionts. Since a complete genome annotation is available for both, we built the metabolic network of two endosymbiotic bacteria, Sulcia muelleri and Baumannia cicadellinicola, that live inside specific cells of the sharpshooter Homalodisca coagulata and studied the metabolic exchanges involving transfers of carbon atoms between the three. We automatically determined the set of metabolites potentially exogenously acquired (seeds) for both metabolic networks. We show that the number of seeds needed by both bacteria in the carbon metabolism is extremely reduced. Moreover, only three seeds are common to both metabolic networks, indicating that the complementarity of the two metabolisms is not only manifested in the metabolic capabilities of each bacterium, but also by their different use of the same environment. Furthermore, our results show that the carbon metabolism of S. muelleri may be completely independent of the metabolic network of B. cicadellinicola. On the contrary, the carbon metabolism of the latter appears dependent on the metabolism of S. muelleri, at least for two essential amino acids, threonine and lysine. Next, in order to define which subsets of seeds (precursor sets) are sufficient to produce the metabolites involved in a symbiotic function, we used a graph-based method, PITUFO, that we recently developed. Our results highly refine our knowledge about the complementarity between the metabolisms of the two bacteria and their host. We thus indicate seeds that appear obligatory in the synthesis of metabolites are involved in the symbiotic function. Our results suggest both B. cicadellinicola and S. muelleri may be completely independent of the metabolites provided by the co-resident endocytobiont to produce the carbon backbone of the metabolites provided to the symbiotic system (., thr and lys are only exploited by B. cicadellinicola to produce its proteins).     

Systematic NMR Analysis of Stable Isotope Labeled Metabolite Mixtures in Plant and Animal Systems: Coarse Grained Views of Metabolic Pathways

Background

Metabolic phenotyping has become an important ‘bird''s-eye-view’ technology which can be applied to higher organisms, such as model plant and animal systems in the post-genomics and proteomics era. Although genotyping technology has expanded greatly over the past decade, metabolic phenotyping has languished due to the difficulty of ‘top-down’ chemical analyses. Here, we describe a systematic NMR methodology for stable isotope-labeling and analysis of metabolite mixtures in plant and animal systems.

Methodology/Principal Findings

The analysis method includes a stable isotope labeling technique for use in living organisms; a systematic method for simultaneously identifying a large number of metabolites by using a newly developed HSQC-based metabolite chemical shift database combined with heteronuclear multidimensional NMR spectroscopy; Principal Components Analysis; and a visualization method using a coarse-grained overview of the metabolic system. The database contains more than 1000 ¹H and ¹³C chemical shifts corresponding to 142 metabolites measured under identical physicochemical conditions. Using the stable isotope labeling technique in Arabidopsis T87 cultured cells and Bombyx mori, we systematically detected >450 HSQC peaks in each ¹³C-HSQC spectrum derived from model plant, Arabidopsis T87 cultured cells and the invertebrate animal model Bombyx mori. Furthermore, for the first time, efficient ¹³C labeling has allowed reliable signal assignment using analytical separation techniques such as 3D HCCH-COSY spectra in higher organism extracts.

Conclusions/Significance

Overall physiological changes could be detected and categorized in relation to a critical developmental phase change in B. mori by coarse-grained representations in which the organization of metabolic pathways related to a specific developmental phase was visualized on the basis of constituent changes of 56 identified metabolites. Based on the observed intensities of ¹³C atoms of given metabolites on development-dependent changes in the 56 identified ¹³C-HSQC signals, we have determined the changes in metabolic networks that are associated with energy and nitrogen metabolism.     

Metabolic profiles of sunflower genotypes with contrasting response to Sclerotinia sclerotiorum infection

We report a comprehensive primary metabolite profiling of sunflower (Helianthus annuus) genotypes displaying contrasting behavior to Sclerotinia sclerotiorum infection. Applying a GC-MS-based metabolite profiling approach, we were able to identify differential patterns involving a total of 63 metabolites including major and minor sugars and sugar alcohols, organic acids, amino acids, fatty acids and few soluble secondary metabolites in the sunflower capitulum, the main target organ of pathogen attack. Metabolic changes and disease incidence of the two contrasting genotypes were determined throughout the main infection period (R5.2-R6). Both point-by-point and non-parametric statistical analyses showed metabolic differences between genotypes as well as interaction effects between genotype and time after inoculation. Network correlation analyses suggested that these metabolic changes were synchronized in a time-dependent manner in response to the pathogen. Concerted differential metabolic changes were detected to a higher extent in the susceptible, rather than the resistant genotype, thereby allowing differentiation of modules composed by intermediates of the same pathway which are highly interconnected in the susceptible line but not in the resistant one. Evaluation of these data also demonstrated a genotype specific regulation of distinct metabolic pathways, suggesting the importance of detection of metabolic patterns rather than specific metabolite changes when looking for metabolic markers differentially responding to pathogen infection. In summary, the GC-MS strategy developed in this study was suitable for detection of differences in carbon primary metabolism in sunflower capitulum, a tissue which is the main entry point for this and other pathogens which cause great detrimental impact on crop yield.     

Thermodynamically based profiling of drug metabolism and drug-drug metabolic interactions: a case study of acetaminophen and ethanol toxic interaction

Drug-drug metabolic interactions can result in unwanted side effects, including reduced drug efficacy and formation of toxic metabolic intermediates. In this work, thermodynamic constraints on non-equilibrium metabolite concentrations are used to reveal the biochemical interactions between the metabolic pathways of ethanol and acetaminophen (N-acetyl-p-aminophenol), two drugs known to interact unfavorably. It is known that many reactions of these pathways are coupled to the central energy metabolic reactions through a number of metabolites and the cellular redox potential. Based on these observations, a metabolic network model has been constructed and a database of thermodynamic properties for all participating metabolites and reactions has been compiled. Constraint-based computational analysis of the feasible metabolite concentrations reveals that the non-toxic pathways for APAP metabolism and the pathway for detoxifying N-acetyl-p-benzoquinoneimine (NAPQI) are inhibited by network interactions with ethanol metabolism. These results point to the potential utility of thermodynamically based profiling of metabolic network interactions in screening of drug candidates and analysis of potential toxicity.     

Metabolic responses of Thellungiella halophila/salsuginea to biotic and abiotic stresses: Metabolite profiles and quantitative analyses

The metabolite profiles of the model crucifer Thellungiella salsuginea (salt cress) ecotype Shandong subjected to various biotic and abiotic stresses were analyzed using HPLC-DAD-ESI-MS. Two different cruciferous microbial pathogens, Albugo candida, a biotrophic oomycete, and Leptosphaeria maculans, a necrotrophic fungus, elicited formation of the phytoalexins wasalexins A and B without causing visual damage on inoculated leaves. Analyses of non-polar and polar metabolites led to elucidation of the chemical structures of five metabolites: 4′-O-(E)-sinapoyl-7-methoxyisovitexin-2″-O-β-d-glucopyranoside, 4′-O-(E)-sinapoylisovitexin-2″-O-β-d-glucopyranoside, 4-O-β-d-glucopyranosyl-7-hydroxymatairesinol, 5′-O-β-d-glucopyranosyldihydroneoascorbigen and 3-O-β-d-glucopyranosylthiane. 3-O-β-d-glucopyranosylthiane, an unique metabolite for which we suggest the name glucosalsuginin, is proposed to derive from the glucosinolate glucoberteroin. In addition, the identification of a broad range of polar metabolites identical to those of other crucifers was carried out. Quantification of several metabolites over a period of eight days showed that concentrations of the polar phytoanticipin 4-methoxyglucobrassicin increased substantially in leaves irradiated with UV light (λ_max 254 nm) relative to control leaves, but not in leaves subjected to other stresses.     

Control analysis of metabolic networks. 2. Total differentials and general formulation of the connectivity relations

The mathematical background of the connectivity relations of metabolic control theory is analysed. The connectivity relations are shown to reflect general properties of total differentials of reaction rate vi, flux J, and metabolite concentration Xj. Connectivity relations hold for any metabolic network in which all vi are homogeneous functions of enzyme concentration Ei. This notion allows established algebraic methods to be used for the formulation of connectivity relations for metabolic systems in which numerous constraints are imposed on metabolite concentrations. A general procedure to derive connectivity relations for such metabolic systems is given. To encourage a broader audience to apply control theory to physiological systems, an easy-to-use graphical procedure is derived for formulating connectivity relations for biochemical systems in which no metabolite is involved in more than one constraint.     

Differences between human plasma and serum metabolite profiles

Background

Human plasma and serum are widely used matrices in clinical and biological studies. However, different collecting procedures and the coagulation cascade influence concentrations of both proteins and metabolites in these matrices. The effects on metabolite concentration profiles have not been fully characterized.

Methodology/Principal Findings

We analyzed the concentrations of 163 metabolites in plasma and serum samples collected simultaneously from 377 fasting individuals. To ensure data quality, 41 metabolites with low measurement stability were excluded from further analysis. In addition, plasma and corresponding serum samples from 83 individuals were re-measured in the same plates and mean correlation coefficients (r) of all metabolites between the duplicates were 0.83 and 0.80 in plasma and serum, respectively, indicating significantly better stability of plasma compared to serum (p = 0.01). Metabolite profiles from plasma and serum were clearly distinct with 104 metabolites showing significantly higher concentrations in serum. In particular, 9 metabolites showed relative concentration differences larger than 20%. Despite differences in absolute concentration between the two matrices, for most metabolites the overall correlation was high (mean r = 0.81±0.10), which reflects a proportional change in concentration. Furthermore, when two groups of individuals with different phenotypes were compared with each other using both matrices, more metabolites with significantly different concentrations could be identified in serum than in plasma. For example, when 51 type 2 diabetes (T2D) patients were compared with 326 non-T2D individuals, 15 more significantly different metabolites were found in serum, in addition to the 25 common to both matrices.

Conclusions/Significance

Our study shows that reproducibility was good in both plasma and serum, and better in plasma. Furthermore, as long as the same blood preparation procedure is used, either matrix should generate similar results in clinical and biological studies. The higher metabolite concentrations in serum, however, make it possible to provide more sensitive results in biomarker detection.