Metabolomic correlation-network modules in Arabidopsis based on a graph-clustering approach
- Atsushi Fukushima^{1},
- Miyako Kusano^{1, 2},
- Henning Redestig^{1},
- Masanori Arita^{1, 3, 4} and
- Kazuki Saito^{1, 5}Email author
https://doi.org/10.1186/1752-0509-5-1
© Fukushima et al; licensee BioMed Central Ltd. 2011
Received: 25 June 2010
Accepted: 1 January 2011
Published: 1 January 2011
Abstract
Background
Deciphering the metabolome is essential for a better understanding of the cellular metabolism as a system. Typical metabolomics data show a few but significant correlations among metabolite levels when data sampling is repeated across individuals grown under strictly controlled conditions. Although several studies have assessed topologies in metabolomic correlation networks, it remains unclear whether highly connected metabolites in these networks have specific functions in known tissue- and/or genotype-dependent biochemical pathways.
Results
In our study of metabolite profiles we subjected root tissues to gas chromatography-time-of-flight/mass spectrometry (GC-TOF/MS) and used published information on the aerial parts of 3 Arabidopsis genotypes, Col-0 wild-type, methionine over-accumulation 1 (mto1), and transparent testa4 (tt4) to compare systematically the metabolomic correlations in samples of roots and aerial parts. We then applied graph clustering to the constructed correlation networks to extract densely connected metabolites and evaluated the clusters by biochemical-pathway enrichment analysis. We found that the number of significant correlations varied by tissue and genotype and that the obtained clusters were significantly enriched for metabolites included in biochemical pathways.
Conclusions
We demonstrate that the graph-clustering approach identifies tissue- and/or genotype-dependent metabolomic clusters related to the biochemical pathway. Metabolomic correlations complement information about changes in mean metabolite levels and may help to elucidate the organization of metabolically functional modules.
Keywords
Background
Combining and integrating different 'omics' data such as transcript-, protein-, and metabolite levels and enzyme activities is essential for a full understanding of the nature of the cellular metabolism as a system [1–4]. With respect to transcript levels, a large amount of microarray data is publicly available for Arabidopsis thaliana, a model plant. Such large datasets facilitate the construction of gene co-expression databases [5] and the survey of transcriptome organization [6–8]. Integrating transcript- and metabolite data by, for example, studying the correlation relationships among profiled data, facilitates the characterization of unknown gene functions, and furthers our understanding of plant cellular systems [9–11].
The correlation between variables (e.g. genes and metabolites) is also important for multivariate statistical analyses such as principal component analysis (PCA) and hierarchical cluster analysis. Typical metabolite-profiling data show a few, but significant correlations among metabolite levels when data sampling is repeated across individuals grown under strictly controlled conditions [12]. The metabolomic correlation as well as gene co-expression are not always in agreement with known biochemical pathways. Metabolomic correlation approaches have highlighted some properties (e.g. modularity and scale-freeness) in several species including plants [13–16]. Steuer et al. [17], who provided a relationship between the structure of a metabolomic-correlation network and a metabolic reaction network using a Jacobian matrix, found that the relationship is not simple. They pointed out that small fluctuations such as glucose availability can result in a certain correlation pattern and persist through metabolic pathways. Using metabolic control analysis (MCA) and correlation analysis based on metabolomic data, Camacho et al. [18] suggested that metabolites are strongly correlated when they respond in the same directions to all perturbations (fluctuations) in enzyme levels. For example, mass conservation and chemical equilibrium were suggested as one origin of a high correlation. Muller-Linow et al. [19] applied network similarity, a graph-theoretic parameter, to compare metabolomic correlation networks with biochemical reactions derived from the KEGG database [20]. They reported that these networks were in disagreement and that closeness in metabolomic correlation is not an indicator of closeness in biochemical networks. Studies on the effect of changes in environmental conditions and temporal- and spatial assessments of the topology of metabolomic correlation networks have been reported [19, 21, 22]. Further investigation of the properties of metabolomic correlation networks may discover whether highly connected metabolites, the so-called 'modules', in the correlation network reflect known biochemical pathways.
We investigated similarities and dissimilarities in metabolomic correlations in the aerial parts of 3 Arabidopsis genotypes, Col-0 wild-type (WT), methionine-over accumulation 1 (mto1) [23], and transparent testa4 (tt4) [24]. Elsewhere [25] we reported that the mutation in cystathionine γ-synthase (CGS) and/or the over-accumulation of methionine (Met) strongly affect the correlation networks in aerial parts of mto1. In the present study, using gas chromatography-time-of-flight/mass spectrometry (GC-TOF/MS), we measured the relative metabolite levels in root samples of the 3 Arabidopsis genotypes to assess tissue- and/or genotype-dependent changes in their metabolite levels. We systematically compared the metabolomic correlations observed in 2 different datasets, the roots and the aerial parts. Multivariate statistical analyses showed the distinct metabolome of these plants and tissues. We then constructed correlation networks by pair-wise correlation between the metabolites and performed graph clustering using the DPClus algorithm [26] that efficiently extracts densely connected metabolites in a large-scale network. We then evaluated the obtained clusters with KEGG [20] enrichment analysis. Our results demonstrate that changes in each network topology are tissue- and/or genotype-dependent and that they reflect, at least partially, known biochemical pathways in Arabidopsis.
Results
Metabolic phenotypes of the roots of 3 Arabidopsis genotypes
The number of significant correlations depends on the tissue/genotype and their differential correlations
We first examined the number of significant correlations between tissues. We found that the number of significant correlations in WT and mto1 roots was larger than in their aerial parts (Figure 3A); it was also much larger in the aerial parts of tt4 than its roots. Then, we compared differences between correlations in roots and aerial parts using Fisher Z-transformation (see Methods). Similarly, we also compared differential correlations between genotypes. The number of statistical differences between correlations in roots and aerial parts in WT, mto1, and tt4 was 7, 8, and 21, although there were no significant changes in the correlations between genotypes except for 6 differential correlations between WT and tt4 in roots (fdr < 0.05) (Figure 3B and Additional file 4). Our findings imply that in both tissue-types, the mto1 mutation strongly affects the metabolomic correlation network, and that there are many correlations intensified by a lack of the chalcone synthase (CHS) gene. Our observations also suggest that the tendency for changes in the number of significant correlations is similar within the studied tissues of the 3 genotypes.
Construction of a metabolomic correlation network and selection of the correlation threshold
For visualization and to gain insights into the metabolomic correlations for each tissue and genotype we constructed correlation networks. The selection of an appropriate threshold is important for the construction of such a correlation network, if it is too high the generated network will be sparse, if too low, only a large connected component will be produced. To assess threshold-dependent changes in the network topology we computed 6 graph-theoretic measures (Additional files 5 and 6), i.e., the graph density, clustering coefficient, average degree, average path length, number of connected components, and the number of edges (for example, see [30]). We also compared the statistics of the real data and 100 randomized data (dots in Additional files 5 and 6). We observed that the number of connected components showed a transition from small to large at a correlation threshold of 0.5 (arrows in Additional file 5). Therefore, we set the threshold at r ≥ 0.5. Although a threshold does not guarantee explicit biological significance, we found that examining various statistics among the constructed networks was useful for selection.
Graph clustering identifies modules significantly enriched for metabolites contained in biochemical pathways
Two-by-two contingency table
within DPClus cluster | without DPClus cluster | total | |
---|---|---|---|
KEGG pathway | a | b | a+b |
non-KEGG pathway | c | d | c+d |
totals | a+c | b+d | n |
List of enriched KEGG pathways in each cluster detected by DPClus
p-value_{min} | Metabolism | |
---|---|---|
WT, roots | ||
Clus1 | 0.0002 | Aminoacyl-tRNA biosynthesis |
Clus2 | n.s. | - |
Clus3 | 0.0061 | Steroid biosynthesis |
Clus4 | 0.0186 | Galactose metabolism |
Clus5 | 0.0041 | Carbon fixation in photosynthetic organisms |
mto1, roots | ||
Clu1 | 0.0021 | Alanine, aspartate and glutamate metabolism |
Clu2 | 0.0311 | Ascorbate and aldarate metabolism |
Clu3 | 0.0025 | Steroid biosynthesis |
Clu4 | 0.0009 | Biosynthesis of unsaturated fatty acids |
Clu5 | n.s. | - |
tt4, roots | ||
Clus1 | 0.0224 | Biosynthesis of alkaloids derived from shikimate pathway |
Clus2 | 0.0102 | Steroid biosynthesis |
Clus3 | 0.0012 | Nitrogen metabolism |
Clus4 | 0.0139 | Purine metabolism |
WT, aerial parts | ||
Clus1 | 0.0017 | Galactose metabolism |
Clus2 | 0.0422 | Biosynthesis of plant hormones |
Clus3 | 0.0041 | Glucosinolate biosynthesis |
Clus4 | 0.0220 | Glycerophospholipid metabolism |
Clus5 | 0.0005 | Biosynthesis of unsaturated fatty acids |
mto1, aerial parts | ||
Clus1 | 0.0205 | Alanine, aspartate and glutamate metabolism |
Clus2 | 0.0019 | Galactose metabolism |
Clus3 | n.s. | - |
Clus4 | 0.0146 | Tropane, piperidine and pyridine alkaloid biosynthesis |
Clus5 | n.s. | - |
Clus6 | n.s. | - |
Clus7 | n.s. | - |
tt4, aerial parts | ||
Clus1 | 0.0197 | Aminoacyl-tRNA biosynthesis |
Clus2 | 0.0038 | Steroid biosynthesis |
Clus3 | 0.0000 | Galactose metabolism |
Clus4 | 0.0045 | Biosynthesis of unsaturated fatty acids |
Clus5 | n.s. | - |
Predominant pathways in clusters reflect tissue- and/or genotype-dependent metabolomic regulation
Except for the clusters found in the aerial parts of WT and mto1, the detected clusters tended to include steroid biosynthetic pathways (Table 2). The strong correlation between sterols confirmed our earlier observations for aerial parts [25]. Metabolites belonging to fatty-acid biosynthesis tended to cluster together. The clusters we commonly found in the aerial parts of the 3 genotypes were galactose metabolism. The largest cluster for the aerial parts and roots of mto1 were 'alanine, aspartate, and glutamate metabolism.' This suggests that the mutation resulting in Met over-accumulation can affect the metabolic pathway of the aspartate family. Metabolite enrichment analysis using other mto mutants also extracted the same pathways [31]. The largest cluster in the roots of tt4 was 'biosynthesis of alkaloids derived from the shikimate pathway.' For both tissues, tt4, which lacks the gene encoding CHS [24], showed less pronounced changes in the metabolite levels belonging to this biosynthetic pathway (see Discussion and Additional file 2).
Discussion
We focused on the changes in the topology of metabolomic correlation networks of 3 Arabidopsis genotypes (WT, mto1, and tt4). When comparing the WT, the increases and decreases in the number of significant correlations for mto1 and tt4 were similar in their aerial parts and roots. Our comparative studies showed strongly altered metabolomic correlations in the roots of mto1. As did Kusano et al. [25] in aerial parts, we observed a marked increase in the number of significant correlations in the roots of tt4 compared to the roots of WT. As our study imposed no environmental perturbations on the 3 genotypes, we suggest that the observed changes in the topology of metabolomic correlation networks are attributable to underlying tissue- and/or genotype-dependent biochemical regulations.
Graph clustering using DPClus [26] yielded densely connected metabolites on the metabolomic correlation networks. KEGG enrichment analysis to assess the statistical significance of the detected clusters demonstrated specific differences in the clusters in the enriched pathways. We posit that the assigned KEGG pathways for each cluster reflect differences in underlying genetic properties and in tissue- and/or genotype-dependent biochemical regulations.
The largest cluster in mto1 roots was assigned as the'alanine, aspartate, and glutamate metabolism'. Since mto1 mutants are characterized by mutations involving the aspartate pathway, the approach based on graph-clustering proved itself useful for characterizing genotypes. Although there were no significant enriched KEGG pathways in cluster 5 in Figure 4, this cluster may represent the extensive coordination among biosynthetic pathways involved in glycolysis, the tricarboxylic acid (TCA)-cycle, and stress-responsive metabolites in mto1 roots. Our approach may yield new insights into the organization of metabolites in the functional pathways of a given organism. In tt4 roots, the largest cluster was 'biosynthesis of alkaloids derived from the shikimate pathway'. This cluster contained metabolites associated with the biosynthetic pathways of the TCA cycle and with glycolysis because this KEGG pathway covered extensive pathways from glycolysis to alkaloids via shikimate for the biosynthesis of aromatic amino acids including phenylalanine, tyrosine, and tryptophan. The mutation in the CHS gene appears not to affect the biosynthesis of alkaloids in Arabidopsis. We posit that the model plant Arabidopsis lacks the production of complex alkaloids, although homologues for many enzymes associated with alkaloid biosynthesis have been detected in its genome [32]. Our multivariate statistical analysis suggests that the mutation of CHS may influence nitrogen assimilation (Additional file 3). Furthermore, the third cluster (Clus3) in tt4 roots was significantly enriched for metabolites contained in 'nitrogen metabolism' (Table 2), suggesting that our approach may reflect underlying changes in the metabolism by lacking of CHS in Arabidopsis. Taken together, our observations demonstrate that variations in the topology of correlation networks are tissue and/or genotype-dependent and reflect at least partially known biochemical pathways in Arabidopsis.
Earlier studies indicated that the interpretation of metabolomic correlations requires careful evaluation [12, 17, 18]. Our work demonstrates that graph clustering can be used to gather metabolites belonging to the metabolic pathways that change in response to different regulations that in turn are dependent on the tissue and/or genotype. Our analyses differ from those of [19] in that we examined the densely connected metabolites in the correlation network by statistical KEGG enrichment analysis instead of using network similarity or proximity. Graph clustering approaches have been effectively applied to gene co-expression networks for extracting functional, densely connected genes [6, 33–35]. We showed that the approach is also effective for metabolomic correlations.
Although we only evaluated 3 genotypes, the approach based on graph clustering seems to be useful and applicable to the monitoring of changes in the topology of correlation networks. As in the gene co-expression analysis that involves thousands of microarrays, there were few negative correlations in this study. For instance, the number of significant negative correlations (r < 0 and fdr < 0.05) in the roots and aerial parts of WT was 17 and 6, respectively. Therefore, we deliberately did not use negative correlations. Other similarity measures such as mutual information [36] and partial correlation [37], instead of Spearman correlation coefficients, may be applied in this approach. Further studies are necessary to assess to what extent these measures affect the graph clustering results of metabolomic correlation networks.
There are few systematic comparisons of metabolomic correlations among time series data [19]. Temporal and spatial assessment of the topology in the metabolomic correlation network may contribute to the development of novel data-mining methods and to the discovery of biomarkers [38]. There are also few published systematic metabolomic comparisons under stress conditions. Szymanski et al. [21] demonstrated that the correlation network approach can help to identify stress-dependent regulatory mechanisms in Escherichia coli. Our current report is a logical extension of our continued effort to understand the regulation of the primary (especially Met biosynthesis [25]) and the secondary (especially flavonoid biosynthesis [39]) metabolism in Arabidopsis. Our approach contributes to the generation of new testable hypotheses for further experiments and may expand our fundamental understanding of the metabolic behaviors affected by genetic and/or environmental perturbations. Such information will contribute to the characterization of unknown gene function(s) and help in the high-throughput screening of metabolic phenotypes ranging from experimental models to crops [40] and disease studies [41].
Conclusions
This study demonstrated that (1) differences in mutations and/or tissues affect changes in the topology of metabolomic correlation networks; (2) graph clustering yields tissue- and/or genotype-dependent metabolomic clusters related to respective biochemical pathways; and (3) metabolomic correlations complement information on changes in mean metabolite levels and the approach based on metabolomic correlations yields insights into the organization and regulation of metabolic modules.
Methods
Plant materials, growth, and harvest
For metabolomic correlation analysis we sampled tissues from the roots of Arabidopsis thaliana Col-0 wild-type (WT), methionine over-accumulation 1 (mto1) [23], and transparent testa4 (tt4) [24]. Data on the aerial part were from [25]. The number of biological replicates in roots was 53 (17 × WT, 16 × mto1, and 20 × tt4). We pooled 3 root samples as a batch for metabolite profiling.
Metabolite profiling
As described in [25], each sample was extracted, derivatized, and analyzed by GC-TOF/MS. Briefly, we pre-processed all raw data using custom MATLAB software (version 6.5; Mathworks, Natick, MA, USA) for hyphenated data analysis (HDA) [42]; it performs baseline correction, peak alignment, and peak deconvolution. For metabolite identification we used the Golm Metabolome Database (GMD) [27, 43] and our in-house mass spectral libraries.
Statistical data analysis
Metabolomic data were log2 transformed and then statistically analyzed using the rank product method [44] to identify differentially changed metabolites with the Bioconductor 'RankProd' package. Significantly changed metabolites showed a false discovery rate (FDR) < 0.05. The FDR value in the rank product was obtained with 1,000 random permutations. PCA was performed with the Bioconductor 'pcaMethods' package [45].
To reduce outlier-related artifacts, the Spearman rank correlation (r), which is relatively robust to outliers and is based only on monotonicity, was used for comparative metabolomic correlations. The statistical significance of the correlation coefficient (r) was tested against the hypothesis of no correlation (r = 0) according to t-statistics: , where n shows the sample size.
where r_{1} and r_{2} are the correlation coefficient for each of the two tissues (or two genotypes), and n_{1} and n_{2} indicate the number of replicates for each of the two tissues (or genotypes) for each metabolite-metabolite pair. The local false-discovery rate (fdr) was determined with the 'fdrtool' package [29]. The significance level was set at fdr < 0.05. All statistical tests were performed with the R program http://www.r-project.org/.
Network analyses and graph clustering
To identify co-accumulated metabolite groups we used DPClus [26], a graph clustering algorithm that can extract densely connected nodes as a cluster. It is based on density- and periphery tracking of clusters. DPClus is freely available from http://kanaya.naist.jp/DPClus/. In this study, we used the overlapping-mode with the DPClus settings because we are confident that extracting clusters with the overlapping-mode is consistent with the overlapping of many of the metabolic pathways and protein complexes. We set the parameter settings of cluster property c_{ p } ; density values were set to 0.5. The resulting clusters are listed in Additional file 7. All network statistics such as graph density were calculated in R with the 'igraph' package [46].
Generation of randomized clusters and assessment of clustered modules in metabolomic correlation networks
To assess the statistical significance of clusters obtained by graph clustering, we subjected datasets of the aerial parts and roots of the 3 Arabidopsis genotypes and 100 randomly-generated sets of clusters to KEGG enrichment. Each randomized cluster was created by permutation of the metabolite names without changing the cluster size. The assessment was done using the S-value [6]: , where n shows the number of clusters and i a cluster. This value is based on the best p-value, p_{ min } , for KEGG pathway enrichment in each cluster. The best p-values were averaged over all clusters to provide the S-value.
The significance of the KEGG pathway in clusters is represented by the p-value, which shows Fisher's exact probabilities based on two-by-two contingency tables: Table 1 is a 2 × 2 table ("KEGG pathway", "non-KEGG pathway", "within DPClus cluster", "without DPClus cluster") used by Fisher's exact test, where significance is defined as a p-value (obtained from the test for DPClus clusters) that is less than or equal to a pre-set level of significance (p < 0.05). The probability of obtaining any set of values was given by the hypergeometric distribution. A full list of the KEGG pathways used in this study is shown in Additional file 9.
Declarations
Acknowledgements
The authors thank Dr. Satoshi Naito, Hokkaido University, and Dr. Satoshi Kitamura, Japan Atomic Energy Research Institute, for providing mto1 and tt4 seeds, respectively, Mr. Makoto Kobayashi and Ms. Naomi Hayashi, RIKEN Plant Science Center, for technical assistance, Drs. Keiko Yonekura-Sakakibara and Takayuki Tohge, RIKEN Plant Science Center, for discussions, Ms. Ursula Petralia for editing this manuscript, Mr. Tetsuya Sakurai and Mr. Kenji Akiyama, RIKEN Plant Science Center, for computational assistance, and two anonymous referees for their critical reading of the manuscript and their insightful comments.
Authors’ Affiliations
References
- Bino RJ, Hall RD, Fiehn O, Kopka J, Saito K, Draper J, Nikolau BJ, Mendes P, Roessner-Tunali U, Beale MH, et al.: Potential of metabolomics as a functional genomics tool. Trends Plant Sci. 2004, 9 (9): 418-425. 10.1016/j.tplants.2004.07.004View ArticlePubMedGoogle Scholar
- Hall RD: Plant metabolomics: from holistic hope, to hype, to hot topic. New Phytol. 2006, 169 (3): 453-468. 10.1111/j.1469-8137.2005.01632.xView ArticlePubMedGoogle Scholar
- Stitt M, Sulpice R, Keurentjes J: Metabolic networks: how to identify key components in the regulation of metabolism and growth. Plant Physiol. 2010, 152 (2): 428-444. 10.1104/pp.109.150821PubMed CentralView ArticlePubMedGoogle Scholar
- Saito K, Matsuda F: Metabolomics for functional genomics, systems biology, and biotechnology. Annu Rev Plant Biol. 2010, 61: 463-489. 10.1146/annurev.arplant.043008.092035View ArticlePubMedGoogle Scholar
- Saito K, Hirai MY, Yonekura-Sakakibara K: Decoding genes with coexpression networks and metabolomics - 'majority report by precogs'. Trends Plant Sci. 2008, 13 (1): 36-43. 10.1016/j.tplants.2007.10.006View ArticlePubMedGoogle Scholar
- Mentzen WI, Wurtele ES: Regulon organization of Arabidopsis. BMC Plant Biol. 2008, 8: 99- 10.1186/1471-2229-8-99PubMed CentralView ArticlePubMedGoogle Scholar
- Fukushima A, Wada M, Kanaya S, Arita M: SVD-based anatomy of gene expressions for correlation analysis in Arabidopsis thaliana. DNA Res. 2008, 15 (6): 367-374. 10.1093/dnares/dsn025PubMed CentralView ArticlePubMedGoogle Scholar
- Mao L, Van Hemert JL, Dash S, Dickerson JA: Arabidopsis gene co-expression network and its functional modules. BMC Bioinformatics. 2009, 10: 346- 10.1186/1471-2105-10-346PubMed CentralView ArticlePubMedGoogle Scholar
- Oksman-Caldentey KM, Saito K: Integrating genomics and metabolomics for engineering plant metabolic pathways. Curr Opin Biotechnol. 2005, 16 (2): 174-179. 10.1016/j.copbio.2005.02.007View ArticlePubMedGoogle Scholar
- Yuan JS, Galbraith DW, Dai SY, Griffin P, Stewart CN: Plant systems biology comes of age. Trends Plant Sci. 2008, 13 (4): 165-171. 10.1016/j.tplants.2008.02.003View ArticlePubMedGoogle Scholar
- Fukushima A, Kusano M, Redestig H, Arita M, Saito K: Integrated omics approaches in plant systems biology. Curr Opin Chem Biol. 2009, 13 (5-6): 532-538. 10.1016/j.cbpa.2009.09.022View ArticlePubMedGoogle Scholar
- Steuer R: Review: on the analysis and interpretation of correlations in metabolomic data. Brief Bioinform. 2006, 7 (2): 151-158. 10.1093/bib/bbl009View ArticlePubMedGoogle Scholar
- Roessner U, Luedemann A, Brust D, Fiehn O, Linke T, Willmitzer L, Fernie A: Metabolic profiling allows comprehensive phenotyping of genetically or environmentally modified plant systems. Plant Cell. 2001, 13 (1): 11-29. 10.1105/tpc.13.1.11PubMed CentralView ArticlePubMedGoogle Scholar
- Kose F, Weckwerth W, Linke T, Fiehn O: Visualizing plant metabolomic correlation networks using clique-metabolite matrices. Bioinformatics. 2001, 17 (12): 1198-1208. 10.1093/bioinformatics/17.12.1198View ArticlePubMedGoogle Scholar
- Martins AM, Camacho D, Shuman J, Sha W, Mendes P, Shulaev V: A Systems Biology Study of Two Distinct Growth Phases of Saccharomyces cerevisiae Cultures. Curr Genomics. 2004, 5: 649-663. 10.2174/1389202043348643.View ArticleGoogle Scholar
- Weckwerth W, Loureiro ME, Wenzel K, Fiehn O: Differential metabolic networks unravel the effects of silent plant phenotypes. Proc Natl Acad Sci USA. 2004, 101 (20): 7809-7814. 10.1073/pnas.0303415101PubMed CentralView ArticlePubMedGoogle Scholar
- Steuer R, Kurths J, Fiehn O, Weckwerth W: Observing and interpreting correlations in metabolomic networks. Bioinformatics. 2003, 19 (8): 1019-1026. 10.1093/bioinformatics/btg120View ArticlePubMedGoogle Scholar
- Camacho D, de la Fuente A, Mendes P: The origin of correlations in metabolomics data. Metabolomics. 2005, 1 (1): 53-63. 10.1007/s11306-005-1107-3.View ArticleGoogle Scholar
- Muller-Linow M, Weckwerth W, Hutt MT: Consistency analysis of metabolic correlation networks. BMC Syst Biol. 2007, 1: 44- 10.1186/1752-0509-1-44PubMed CentralView ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M: KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010, D355-360. 38 DatabaseView ArticleGoogle Scholar
- Szymanski J, Jozefczuk S, Nikoloski Z, Selbig J, Nikiforova V, Catchpole G, Willmitzer L: Stability of metabolic correlations under changing environmental conditions in Escherichia coli--a systems approach. PLoS ONE. 2009, 4 (10): e7441- 10.1371/journal.pone.0007441PubMed CentralView ArticlePubMedGoogle Scholar
- Urano K, Maruyama K, Ogata Y, Morishita Y, Takeda M, Sakurai N, Suzuki H, Saito K, Shibata D, Kobayashi M, et al.: Characterization of the ABA-regulated global responses to dehydration in Arabidopsis by metabolomics. Plant J. 2009, 57 (6): 1065-1078. 10.1111/j.1365-313X.2008.03748.xView ArticlePubMedGoogle Scholar
- Inaba K, Fujiwara T, Hayashi H, Chino M, Komeda Y, Naito S: Isolation of an Arabidopsis thaliana Mutant, mto1, That Overaccumulates Soluble Methionine (Temporal and Spatial Patterns of Soluble Methionine Accumulation). Plant Physiol. 1994, 104 (3): 881-887.PubMed CentralPubMedGoogle Scholar
- Shikazono N, Yokota Y, Kitamura S, Suzuki C, Watanabe H, Tano S, Tanaka A: Mutation rate and novel tt mutants of Arabidopsis thaliana induced by carbon ions. Genetics. 2003, 163 (4): 1449-1455.PubMed CentralPubMedGoogle Scholar
- Kusano M, Fukushima A, Arita M, Jonsson P, Moritz T, Kobayashi M, Hayashi N, Tohge T, Saito K: Unbiased characterization of genotype-dependent metabolic regulations by metabolomic approach in Arabidopsis thaliana. BMC Syst Biol. 2007, 1: 53- 10.1186/1752-0509-1-53PubMed CentralView ArticlePubMedGoogle Scholar
- Altaf-Ul-Amin M, Shinbo Y, Mihara K, Kurokawa K, Kanaya S: Development and implementation of an algorithm for detection of protein complexes in large interaction networks. BMC Bioinformatics. 2006, 7: 207- 10.1186/1471-2105-7-207PubMed CentralView ArticlePubMedGoogle Scholar
- Schauer N, Steinhauser D, Strelkov S, Schomburg D, Allison G, Moritz T, Lundgren K, Roessner-Tunali U, Forbes MG, Willmitzer L, et al.: GC-MS libraries for the rapid identification of metabolites in complex biological samples. FEBS Lett. 2005, 579 (6): 1332-1337. 10.1016/j.febslet.2005.01.029View ArticlePubMedGoogle Scholar
- Redestig H, Kusano M, Fukushima A, Matsuda F, Saito K, Arita M: Consolidating metabolite identifiers to enable contextual and multi-platform metabolomics data analysis. BMC Bioinformatics. 2010, 11 (1): 214- 10.1186/1471-2105-11-214PubMed CentralView ArticlePubMedGoogle Scholar
- Strimmer K: fdrtool: a versatile R package for estimating local and tail area-based false discovery rates. Bioinformatics. 2008, 24 (12): 1461-1462. 10.1093/bioinformatics/btn209View ArticlePubMedGoogle Scholar
- Barabasi AL, Oltvai ZN: Network biology: understanding the cell's functional organization. Nat Rev Genet. 2004, 5 (2): 101-113. 10.1038/nrg1272View ArticlePubMedGoogle Scholar
- Kusano M, Fukushima A, Redestig H, Kobayashi M, Otsuki H, Onouchi H, Naito S, Hirai MY, Saito K: Comparative metabolomics charts the impact of genotype-dependent methionine accumulation in Arabidopsis thaliana. Amino Acids. 2010, 39 (4): 1013-1021. 10.1007/s00726-010-0562-yView ArticlePubMedGoogle Scholar
- Facchini PJ, Bird DA, St-Pierre B: Can Arabidopsis make complex alkaloids?. Trends Plant Sci. 2004, 9 (3): 116-122. 10.1016/j.tplants.2004.01.004View ArticlePubMedGoogle Scholar
- Freeman TC, Goldovsky L, Brosch M, van Dongen S, Maziere P, Grocock RJ, Freilich S, Thornton J, Enright AJ: Construction, visualisation, and clustering of transcription networks from microarray expression data. PLoS Comput Biol. 2007, 3 (10): 2032-2042. 10.1371/journal.pcbi.0030206View ArticlePubMedGoogle Scholar
- Ma S, Gong Q, Bohnert HJ: An Arabidopsis gene network based on the graphical Gaussian model. Genome Res. 2007, 17 (11): 1614-1625. 10.1101/gr.6911207PubMed CentralView ArticlePubMedGoogle Scholar
- Fukushima A, Kanaya S, Arita M: Characterizing gene coexpression modules in Oryza sativa based on a graph-clustering approach. Plant Biotechnol. 2009, 26: 485-493.View ArticleGoogle Scholar
- Steuer R, Kurths J, Daub CO, Weise J, Selbig J: The mutual information: detecting and evaluating dependencies between variables. Bioinformatics. 2002, 18 (Suppl 2): S231-240.View ArticlePubMedGoogle Scholar
- de la Fuente A, Bing N, Hoeschele I, Mendes P: Discovery of meaningful associations in genomic data using partial correlation coefficients. Bioinformatics. 2004, 20 (18): 3565-3574. 10.1093/bioinformatics/bth445View ArticlePubMedGoogle Scholar
- Tarpley L, Duran AL, Kebrom TH, Sumner LW: Biomarker metabolites capturing the metabolite variance present in a rice plant developmental period. BMC Plant Biol. 2005, 5: 8- 10.1186/1471-2229-5-8PubMed CentralView ArticlePubMedGoogle Scholar
- Yonekura-Sakakibara K, Tohge T, Matsuda F, Nakabayashi R, Takayama H, Niida R, Watanabe-Takahashi A, Inoue E, Saito K: Comprehensive flavonol profiling and transcriptome coexpression analysis leading to decoding gene-metabolite correlations in Arabidopsis. Plant Cell. 2008, 20 (8): 2160-2176. 10.1105/tpc.108.058040PubMed CentralView ArticlePubMedGoogle Scholar
- Albinsky D, Kusano M, Higuchi M, Hayashi N, Kobayashi M, Fukushima A, Mori M, Ichikawa T, Matsui K, Kuroda H, et al.: Metabolomic Screening Applied to Rice FOX Arabidopsis Lines Leads to the Identification of a Gene-Changing Nitrogen Metabolism. Mol Plant. 2010, 3 (1): 125-142. 10.1093/mp/ssp069View ArticlePubMedGoogle Scholar
- de la Fuente A: From 'differential expression' to 'differential networking' - identification of dysfunctional regulatory networks in diseases. Trends Genet. 2010, 26 (7): 326-333. 10.1016/j.tig.2010.05.001View ArticlePubMedGoogle Scholar
- Jonsson P, Johansson ES, Wuolikainen A, Lindberg J, Schuppe-Koistinen I, Kusano M, Sjostrom M, Trygg J, Moritz T, Antti H: Predictive metabolite profiling applying hierarchical multivariate curve resolution to GC-MS data--a potential tool for multi-parametric diagnosis. J Proteome Res. 2006, 5 (6): 1407-1414. 10.1021/pr0600071View ArticlePubMedGoogle Scholar
- Kopka J, Schauer N, Krueger S, Birkemeyer C, Usadel B, Bergmuller E, Dormann P, Weckwerth W, Gibon Y, Stitt M, et al.: GMD@CSB.DB: the Golm Metabolome Database. Bioinformatics. 2005, 21 (8): 1635-1638. 10.1093/bioinformatics/bti236View ArticlePubMedGoogle Scholar
- Breitling R, Armengaud P, Amtmann A, Herzyk P: Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Lett. 2004, 573 (1-3): 83-92. 10.1016/j.febslet.2004.07.055View ArticlePubMedGoogle Scholar
- Stacklies W, Redestig H, Scholz M, Walther D, Selbig J: pcaMethods--a bioconductor package providing PCA methods for incomplete data. Bioinformatics. 2007, 23 (9): 1164-1167. 10.1093/bioinformatics/btm069View ArticlePubMedGoogle Scholar
- Csardi G, Nepusz T: The igraph software package for complex network research. Int J Complex Syst. 2006, http://necsi.org/events/iccs6/viewabstract.php?id=88Google Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.