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.

[1]  Thomas Nägele,et al.  Identification of a metabolic bottleneck for cold acclimation in Arabidopsis thaliana. , 2012, The Plant journal : for cell and molecular biology.

[2]  Joost T. van Dongen,et al.  A rapid approach for phenotype‐screening and database independent detection of cSNP/protein polymorphism using mass accuracy precursor alignment , 2008, Proteomics.

[3]  B. Usadel,et al.  Multilevel genomic analysis of the response of transcripts, enzyme activities and metabolites in Arabidopsis rosettes to a progressive decrease of temperature in the non-freezing range. , 2008, Plant, cell & environment.

[4]  Asaph Aharoni,et al.  Expression of a bacterial feedback-insensitive 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase of the shikimate pathway in Arabidopsis elucidates potential metabolic bottlenecks between primary and secondary metabolism. , 2012, The New phytologist.

[5]  Kazuo Shinozaki,et al.  MS/MS spectral tag-based annotation of non-targeted profile of plant secondary metabolites , 2008, The Plant journal : for cell and molecular biology.

[6]  J. Kopka,et al.  Impact of soluble sugar concentrations on the acquisition of freezing tolerance in accessions of Arabidopsis thaliana with contrasting cold adaptation - evidence for a role of raffinose in cold acclimation , 2004 .

[7]  C. Granger Investigating causal relations by econometric models and cross-spectral methods , 1969 .

[8]  J. Morgan,et al.  Mathematical modeling of plant metabolic pathways. , 2002, Metabolic engineering.

[9]  B. M. Lange,et al.  Experimental and mathematical approaches to modeling plant metabolic networks. , 2007, Phytochemistry.

[10]  Wolfram Weckwerth,et al.  Mathematical Modeling of Plant Metabolism―From Reconstruction to Prediction , 2012, Metabolites.

[11]  O. Fiehn,et al.  Process for the integrated extraction, identification and quantification of metabolites, proteins and RNA to reveal their co‐regulation in biochemical networks , 2004, Proteomics.

[12]  W. Weckwerth Metabolomics in systems biology. , 2003, Annual review of plant biology.

[13]  P. Durek,et al.  Metabolic pathway relationships revealed by an integrative analysis of the transcriptional and metabolic temperature stress-response dynamics in yeast. , 2010, Omics : a journal of integrative biology.

[14]  C Giersch Mathematical modelling of metabolism. , 2000, Current opinion in plant biology.

[15]  B. Winkel-Shirley,et al.  Biosynthesis of flavonoids and effects of stress. , 2002, Current opinion in plant biology.

[16]  K. Kloppstech,et al.  The protective functions of carotenoid and flavonoid pigments against excess visible radiation at chilling temperature investigated in Arabidopsisnpq and tt mutants , 2001, Planta.

[17]  Wolfram Weckwerth,et al.  Unpredictability of metabolism—the key role of metabolomics science in combination with next-generation genome sequencing , 2011, Analytical and bioanalytical chemistry.

[18]  Oliver Sawodny,et al.  Mathematical Modeling of the Central Carbohydrate Metabolism in Arabidopsis Reveals a Substantial Regulatory Influence of Vacuolar Invertase on Whole Plant Carbon Metabolism1[W] , 2010, Plant Physiology.

[19]  Gad Galili,et al.  New insights into the shikimate and aromatic amino acids biosynthesis pathways in plants. , 2010, Molecular plant.

[20]  Wolfram Weckwerth,et al.  COVAIN: a toolbox for uni- and multivariate statistics, time-series and correlation network analysis and inverse estimation of the differential Jacobian from metabolomics covariance data , 2012, Metabolomics.

[21]  W. Weckwerth,et al.  Metabolomics Unravel Contrasting Effects of Biodiversity on the Performance of Individual Plant Species , 2010, PloS one.

[22]  S. Abrahams,et al.  The structure of the major anthocyanin in Arabidopsis thaliana. , 2002, Phytochemistry.

[23]  W. Weckwerth Green systems biology - From single genomes, proteomes and metabolomes to ecosystems research and biotechnology. , 2011, Journal of proteomics.

[24]  Yves Gibon,et al.  GMD@CSB.DB: the Golm Metabolome Database , 2005, Bioinform..

[25]  K. Paszkiewicz,et al.  The influence of ascorbate on anthocyanin accumulation during high light acclimation in Arabidopsis thaliana: further evidence for redox control of anthocyanin synthesis. , 2012, Plant, cell & environment.

[26]  D. Xie,et al.  Features of anthocyanin biosynthesis in pap1-D and wild-type Arabidopsis thaliana plants grown in different light intensity and culture media conditions , 2010, Planta.

[27]  Yves Gibon,et al.  Extensive metabolic cross-talk in melon fruit revealed by spatial and developmental combinatorial metabolomics. , 2011, The New phytologist.

[28]  M. Hirai,et al.  Functional genomics by integrated analysis of metabolome and transcriptome of Arabidopsis plants over-expressing an MYB transcription factor. , 2005, The Plant journal : for cell and molecular biology.

[29]  A. Aharoni,et al.  Expression of a bacterial bi-functional chorismate mutase/prephenate dehydratase modulates primary and secondary metabolism associated with aromatic amino acids in Arabidopsis. , 2009, The Plant journal : for cell and molecular biology.

[30]  Jürgen Kurths,et al.  Observing and Interpreting Correlations in Metabolic Networks , 2003, Bioinform..

[31]  G. Coruzzi,et al.  Light represses transcription of asparagine synthetase genes in photosynthetic and nonphotosynthetic organs of plants , 1991, Molecular and cellular biology.

[32]  J. Wolfender,et al.  Evaluation of quadrupole time-of-flight tandem mass spectrometry and ion-trap multiple-stage mass spectrometry for the differentiation of C-glycosidic flavonoid isomers. , 2001, Journal of chromatography. A.