Metabolite profiling of human colon carcinoma – deregulation of TCA cycle and amino acid turnover
- Carsten Denkert†1Email author,
- Jan Budczies†1, 2,
- Wilko Weichert†1,
- Gert Wohlgemuth3,
- Martin Scholz3,
- Tobias Kind3,
- Silvia Niesporek1,
- Aurelia Noske1,
- Anna Buckendahl1,
- Manfred Dietel1 and
- Oliver Fiehn3
© Denkert et al; licensee BioMed Central Ltd. 2008
Received: 19 December 2007
Accepted: 18 September 2008
Published: 18 September 2008
Apart from genetic alterations, development and progression of colorectal cancer has been linked to influences from nutritional intake, hyperalimentation, and cellular metabolic changes that may be the basis for new diagnostic and therapeutic approaches. However, in contrast to genomics and proteomics, comprehensive metabolomic investigations of alterations in malignant tumors have rarely been conducted.
In this study we investigated a set of paired samples of normal colon tissue and colorectal cancer tissue with gas-chromatography time-of-flight mass-spectrometry, which resulted in robust detection of a total of 206 metabolites. Metabolic phenotypes of colon cancer and normal tissues were different at a Bonferroni corrected significance level of p = 0.00170 and p = 0.00005 for the first two components of an unsupervised PCA analysis. Subsequent supervised analysis found 82 metabolites to be significantly different at p < 0.01. Metabolites were connected to abnormalities in metabolic pathways by a new approach that calculates the distance of each pair of metabolites in the KEGG database interaction lattice. Intermediates of the TCA cycle and lipids were found down-regulated in cancer, whereas urea cycle metabolites, purines, pyrimidines and amino acids were generally found at higher levels compared to normal colon mucosa.
This study demonstrates that metabolic profiling facilitates biochemical phenotyping of normal and neoplastic colon tissue at high significance levels and points to GC-TOF-based metabolomics as a new method for molecular pathology investigations.
Colorectal cancer is the third most common cancer for both sexes, with an estimated number of 153760 new cases and 52180 estimated deaths in the United States in 2007 . Despite advances in diagnosis and surgical treatment of colorectal cancer, incidence as well as mortality of colorectal cancer has decreased only slightly in the last twenty years. While the disease is curable in early stages, the risk of recurrence and metastasis is substantially higher for advanced tumors.
In addition to genetic alterations, the development of colorectal cancer has been linked to a dysregulation of energy homeostasis which is induced by hyperalimentation and adipositas as well as to influences from nutritional intake [2–6]. Recently it has been suggested to use this metabolic deregulation as the basis for new therapeutic interventions .
While genetic alterations have been extensively characterized in colon cancer, the metabolic remodeling that occurs downstream from genomic and proteomic alterations have not been analyzed to a great extent, so far. This is in part due to the fact that low molecular metabolites are comparably difficult to analyze in a large-scale -omics approach. Comparably new technologies allow comprehensive and quantitative investigation of a multitude of different metabolites, which is called "metabolomics" in analogy to the terms "transcriptomics" and "proteomics" [8–10] currently, the analysis by gas chromatography (GC) coupled with time-of flight (TOF) mass spectrometry is used as a standard method for studying primary metabolism [11–13]. While these methods are state-of-the-art in the analysis of plant metabolism-, the translation to cancer research for the characterization of tumor types in molecular pathology is still a challenge .
For cancer research, the ultimate aim is to define metabolic profiles in normal, precancerous and cancerous tissues and to link the dysregulation of tumor metabolism to clinical outcome and treatment response. To reach this aim, three steps are necessary. First, it is necessary to demonstrate that metabolite levels can be determined in human tissue samples that were collected during routine surgical procedures. Second, it is necessary to identify as many metabolites as possible, to go beyond a simple pattern analysis towards an understanding of the functional alterations of metabolism in tumor tissue. Third, we have to develop tools that allow us to easily interpret metabolic data by connection with the existing knowledge on metabolic pathways that is described in databases such as the Kyoto encyclopedia of genes and genomes (KEGG) database .
In human samples, there may be general differences in metabolism that are related to the genetic background or – more likely – to the nutritional intake of the individuals studied. It is therefore not clear, if the metabolic differences between tumor tissue and normal tissue are strong enough to be detected despite the different genetic, metabolic and nutritional background of individual patients.
Therefore, it was one aim of this study to investigate if a set of paired samples of normal colon tissue and colorectal cancer tissue from individual patients can be used for metabolic profiling with GC-MS to detect and interpret molecular changes in tumor tissue and to detect metabolic patterns associated with different biological entities.
We evaluated different bioinformatical strategies to detect those metabolites that are differentially expressed in normal tissue and cancer tissue. To facilitate the interpretation of the results with respect to changes in metabolic pathways, we have developed a new method that projects the metabolite interactions from the multidimensional KEGG interaction lattice to a one-dimensional axis (PROFILE). This method uses the relational information from all metabolic pathways described in the KEGG database but focuses on those metabolites and reactions that can be observed in our investigation and thus builds the bridge to functional interpretation of the metabolomic changes in colon cancer.
Study Population and histopathological examination
For GC-TOF analysis 45 colon samples were examined at the Institute of Pathology, Charité Hospital, Berlin, Germany. The tissue specimens included 27 primary colon carcinomas and 18 normal mucosa samples. For 15 cases paired samples of cancer tissue and normal tissue were available from the same patient. The tissues were dissected by a senior pathologist in the operating room and was immediately frozen in liquid nitrogen and stored at -80°C. Additional H&E sections were performed for histopathological evaluation.
Fresh-frozen biopsy tissues (approximately 5 mg fresh weight) were prepared by grinding in 2 ml Eppendorf tubes for 30 s at 25 s-1 using 3 mm i.d. metal balls in a MM300 ball mill (Retsch, Germany). Subsequent extraction was carried using 1 ml of a one phase mixture of chloroform:methanol:water (2:5:2, v/v/v) at -20°C for 5 min. Tubes were centrifuged for 30 s at 14,000 g and the supernatant was collected and concentrated to complete dryness. Samples were derivatized for GC-TOF analysis as previously published . In order to avoid cross-contamination of samples, 1.5 μl of the derivatized solution was injected in the split less mode into a thermodesorption unit (DTD, ATAS GL, Zoetermeer, Netherlands) equipped with automatic exchange of liners and micro inserts. The sample was introduced at 40°C using a programmable temperature vaporization OPTIC3 injector (ATAS GL, Zoetermeer, Netherlands) and heated to 290°C using a 4°C/min ramp. Mass spectrometry analysis was carried out using an Agilent 6890 gas chromatography oven (Hewlett-Packard, Atlanta, GA, USA) which was coupled to a Pegasus III time of flight (TOF) mass spectrometer from Leco (St Joseph, MI, USA). A MDN-35 fused silica capillary column of 30 m length, 0.32 mm I.D. and 0.25 μm film thickness was used for separation using a start temperature of 85°C which was ramped by 15°C/min to 360°C. Mass spectra were acquired for a scan range of 83–500 m/z and an acquisition rate of 20 spectra per second. The ionization mode was electron impact at 70 eV. The temperature for the ion source was set to 250°C.
As an additional quality control we used stable isotope labeled cholesterol as internal standard for lipophilic and high boiling compounds for every sample as further control on analytical variation. The subsequent analysis of quality control charts of this deuterated cholesterol standard did not show trends over the duration of data acquisition.
Data processing and normalization
Two stages of raw data processing were employed. In the initial step, automated peak detection and mass spectral deconvolution was performed by the Leco ChromaTOF software (v2.32). For each sample, around 700 spectra were exported with absolute intensities (peak heights). These spectra were further processed by the in-house programmed database BinBase . All known artifact peaks such as internal standards, column bleed, plasticizers or reagent peaks were excluded from the result sheets. The algorithm underlying BinBase effectively removes inconsistent signals and noise peaks, yielding a total number of 206 metabolic signals that were reliably determined across the entire sample set. Of these peaks, 107 were structurally identified as known metabolites by comparison to mass spectra and retention indices of customized reference mass spectral libraries that were acquired with authentic standard compounds under identical data acquisition parameters. For metabolites that were below detection limits at a given sample, or for which mass spectra did not match the quality criteria underlying the BinBase algorithm, the data set resulted in missing values. In order to minimize the number of such missing values, only compounds were taken into account that were consistently detected in at least 85% of samples. The metabolite data were normalized relative to the sum of the 107 known metabolites in each sample and transformed to the log scale keeping the missing values in the data set. These missing values were further kept in the dataset for metabolite-wise analyses (t-test, fold change calculation). For collective analyses (PCA, clustering, classification) we have replaced missing metabolite measurements by the corresponding arithmetic means over all samples. The entire data analysis was performed within programming and visualization environment R . In additional to the normalization strategy described above, additional analysis were performed using the raw (unnormalized) data as input.
Principal component analysis (PCA)
Principal component analysis was performed using the function prcomp() in the R package stats. In principal component analysis the original set of metabolites is reduced to a new set of principal components that retain the variance-covariance structure of the data, but use lesser dimensions of data space. Out of the 45 principal components, the fist (PC1) as well as the second (PC2) turned out be significantly different between colon carcinomas and normal mucosa (Welch's t-test, Bonferroni corrected p-values). The values of the first two PCs were plotted with designation of the cases as colon carcinoma or normal mucosa. Further box plots served as visualisation of the discriminative power of the first and the second PC.
Detection of metabolic changes
Alterations between colon carcinomas and normal mucosa were evaluated by thresholds on the fold change and Welch's t-test p-values. The results of three different selection procedures (p < 0.05, p < 0.01, p < 0.00024) were validated by repeated (n = 1000) random permutations of the samples. The false discovery rates (FDR) for the corresponding metabolite lists were estimated as the ratio (nexp/nobs) between the number of observed significant metabolites between colon carcinomas and mucosa tissues (nobs) and the number of metabolites that were expected to be significant by chance from the permutation distribution (nexp).
Prediction of the tissue type
Predictive models were derived in a two-step procedure consisting of a feature selection step followed by the proper construction of the classifier. We have applied the nearest centroid classification (NCC) and compared it to more complex methods like linear discriminant analysis (LDA) and linear support vector machines (SVM). We used an own implementation of the nearest centroid rule, the implementation of LDA in the R package MAS and the implementation of SVMs in the R package e1071. Feature selection was performed by ranking according to the result of Welch's t-test on the training set. The top 2, 3, ..., 206 metabolites were used as input for the construction of the classifier and classification results were monitored in dependence of the number of features. The classification results were obtained in a leave-one-out approach and reported separately the carcinomas and normal tissues. Beyond that we have studied the stability of the classification results under different choices of the training data. To this end we have employed a protocol similar to that in  and trained the classifier on balanced data sets of different size.
Functional analysis of gene signatures by PROFILE clustering
Pro jection f rom i nteraction l attice is a new method for the interpretation of metabolomic changes by the integration of pathway information. Based on the KEGG REACTION data base [cf. ref. ] a metabolism network was build by joining two compounds A and B with an edge whenever there was a reaction that converted A into B. Such reaction was assumed to exist if and only if A and B where annotated as reaction pair with the qualifier "main" in KEGG REACTION. The final network consisted of 4206 metabolites that were annotated as "main" in at least one reaction pair. Based on these concepts, the bio-chemical distance between two metabolites was defined as the minimal length of a joining path within the network. The function johnson.all.pairs.sp() from the R package RBGL was used to calculate the bio-chemical distance between all pairs of the 4206 metabolites. The RBGL package is an interface to the Boost C++ library for graph algorithms.
This result was useful for the interpretation of metabolomic changes in colon tissue. Our study included 206 measurements of which 107 could be mapped to chemical compounds and metabolite names. Out of these compounds, 84 were registered in the KEGG database and 71 were a main reaction partner in at least one of the reactions annotated in KEGG REACTION. After projection to these 71 metabolites, the network turned out to be composed out of a single connection component, i.e. pair of metabolites could be joined by a path inside the network. Next we performed an agglomerative hierarchical clustering of these compounds with respect to the bio-chemical distance to project the complexity of the metabolomic network to an one-dimensional axis. In doing so the average linkage method was used to calculate distances between clusters. Finally, the fold changes between cancer and normal tissues were plotted against the functional axis.
Metabolic profiling of tissue samples by GC-TOF
Clinicopathological characteristics of the study cohort
Number of cases (%)
Normal colon tissue
Cases with paired normal and carcinoma tissue
Characteristics of colon carcinomas
Unsupervised evaluation of metabolite signatures using principal component analysis (PCA)
Identification of differential metabolites by supervised analysis
Estimation of false discovery rate and the number of false positives for three different p-value thresholds. The estimates result from a permutation procedure that is detailed in the material and methods section.
number of detected metabolites
False discovery rate (FDR)
Estimated number of false positive metabolites
As an additional validation step, we used a permutation analysis to estimate the number of false positives detected by different p-value thresholds. After 1000 permutations of the data set, the number of detected metabolites in all of the shuffled data sets was lower than the number of regulated metabolites between cancer and normal tissues. This was true for all investigated thresholds (p = 0.05, p = 0.01, and p = 0.00024). Finally, the false discovery rate (FDR) and the number of false positive detections were estimated from the permutation analysis, see table 2. Only two false positives are expected within in the list of 82 metabolites that was generated by the threshold p = 0.01.
Classification models for prediction of tissue-type
Sensitivity and specificity for the detection of colon cancer with the nearest centroid classifier based on the levels of 10, 25, 50, 100 or 206 metabolites.
number of metabolites in the classifier
Interpretation of metabolic changes in tumor tissue using pathway information from the KEGG database (PROFILE)
Metabolic profiling analyses have the advantage that reliable and standardized information on metabolic pathways exists in the KEGG database. However, with the current techniques in data acquisition and deconvolution it is not possible to measure all the metabolites corresponding to each pathway. Many compounds involved in pathways are intermediates with too low concentrations to be reliably detected. Others are involved in more than one pathway, depending on the localization and distribution between the cellular compartments. KEGG can be interpreted as a multidimensional metabolic network that can be visualized by the KEGG pathway diagrams. However, these diagrams are only partially suitable for visualization and interpretation of metabolic results, since the detected metabolites usually come from different pathways and not all metabolites from a single pathway can be detected. Furthermore, the various cross-links between different pathways that may be relevant as well are very difficult to visualize with the conventional approach of using full KEGG pathway maps. Efforts have been carried out to reduce complexity by constructing static pathway maps ; however, such approaches are not flexible enough to integrate novel information and are too superficial to depict the different roles and connections of metabolites in metabolic networks. In previous reports  we have proposed to exploit differences in metabolic network graphs. On the one hand, such graphs are powerful because both identified and unidentified metabolic signals can be taken into account. On the other hand, the biological relevance and pathway connection of unknown metabolites cannot be known. Therefore, we constrain analyses here to only structurally identified compounds that can be mapped to pathways.
We have developed a new approach to visualize metabolomics data in context of the KEGG database. PROFILE is based on projection of the metabolite interactions from the multidimensional KEGG interaction lattice to a one-dimensional axis. PROFILE starts with the calculation of the metabolic distance of each pair of metabolites in the KEGG interaction lattice. This first step is performed using all metabolites in the KEGG database, without a restriction to those that can be actually measured. Distances between all metabolites are described as the number of main biochemical reactions that are needed to connect these metabolites. The result of this first step is a 4206 × 4206 data matrix that represents metabolite distance across the complete KEGG REACTION database. In a subsequent step, this matrix is reduced to visualize relationships between those metabolites that can actually be measured by GC-TOF profiling of colon tissues.
As shown in figure 7, the PROFILE method clusters many compounds according to major metabolic pathways such as urea cycle, purines or TCA cycle intermediates. Within these pathways, compounds are usually deregulated in the same direction, i.e. urea cycle intermediates, purines and most amino acids up regulated in carcinoma tissues and TCA cycle intermediates and free fatty acids down regulated. Consequently, the clustering of results according to metabolic distances helps annotating significant differences in regulation to pathways, in contrary to the alternative display of metabolic differences by p-values or by fold-changes. Interestingly, the most significant metabolic difference was found for up regulation of beta-alanine in colon carcinoma.
Our study demonstrates that large scale metabolic profiling using GC-TOF mass spectrometry and database annotation yields numerous significant differences between colon carcinoma and normal colon mucosa. We have utilized both unsupervised and supervised approaches to investigate these metabolic differences. The metabolite signatures are capable of predicting the status (normal tissue or colon carcinoma) of a previously unknown test tumor at sensitivity and specificity around 95%. Importantly, we could show that the classification results are robust against different choices of the classificator and the training set (figure 6). Regarding the extent of changes detected, it is important to note that – from a tumorbiological point of view – the comparison of normal tissue and carcinoma tissue means that two completely different entities are compared. Therefore we would expect that comparison of those two tissue types leads to a large set of completely different biomarkers. Similar results have been reported in gene expression analysis. For example in the study by Hlubek et al.  39% of transcripts were differentially expressed between colon tumor center and colon normal tissue. Since metabolites are regarded as an amplified output of a biological system, the expected metabolite changes could be even more prominent compared to genomics, as shown in our analysis.
Many of these metabolic events can be ascribed to known metabolic dysregulation in cancer thus validating the method itself. Metabolites involved in the citric acid cycle were generally found at lower amounts in cancer tissues compared to normal colon samples, in accordance with results published earlier . Purines were detected at increased levels in malignant tissues as indicator for higher capacity for DNA synthetic capacity. Similarly, we found almost all amino acids to be up regulated in carcinoma tissues, which may be interpreted as reflecting cellular needs for higher turnover of structural proteins. This finding is in agreement with earlier publications for select amino acids, notably glutamate and aspartate . Similarly, the high GABA contents had been described in colon cancer tissues in a previous study . Certain amino acids are synthesized by mammalian metabolic routes, often using TCA intermediates as precursor such as alpha-ketoglutarate for glutamate and its derived amino acids, and oxaloacetate for aspartate-derived amino acids. With higher needs in amino acids but lower use of the TCA cycle, an alternative route is needed to deliver carbon backbones for such TCA-derived intermediates. Such higher import may be accomplished by up regulation of amino acid transporter, facilitating higher cellular needs for energy metabolism as well as delivering carbon backbones for biosynthesis of cellular molecules. This interpretation is supported by our finding of increased levels in urea cycle intermediates in colon carcinoma tissues, indicating higher turnover of amino acids. Interestingly, beta-alanine was found as the most upregulated (f.c. = 4.9) metabolite in carcinoma tissues with very high statistical significance (p = 5.8e-13). In humans, beta-alanine is a unidirectional catabolic product from aspartate in a decarboxylation reaction (EC:18.104.22.168) or by catabolic routes from pyrimidine metabolism (EC:22.214.171.124). However, the eventual fate of beta-alanine in humans is yet unclear, since no enzymes are known that would transfer its backbone into acetyl-CoA or towards pantothenate metabolism, as it occurs in other species. We therefore suggest that beta-alanine might be important for metabolic alterations in colon cancer. In addition, we found that not all amino acids were up regulated in the same manner. In fact, the glutamate/glutamine ratio was greatly altered in comparison to normal colon tissue, indicating a lesser role of aminotransferase reactions utilizing glutamine or less need for transport of nitrogen across cells.
A limitation of this study was found in the need of normalizing the raw data to the total sum of known metabolites. The normalization strategy was developed in analogy to gene expression studies and was chosen because frozen tissue sections (as detailed in the methods section) should not be weighed on fine balances in order to preserve the cold chain and to prevent reactivation of metabolism prior to extraction. We found that the raw data for carcinoma tissues were significantly higher (p = 0.03) than those for normal tissue, relating to roughly a 33% increase in overall metabolic levels. However, this might be due to either a higher number of tumor cells per area of tissue or for generally enhanced metabolism. More detailed studies would be needed to address this question. In an additional validation using the raw (unnormalized) data as input the major metabolic differences between both tissue types could be detected, as well, suggesting that the major metabolite differences are not dependent on the normalization strategy.
The total time between surgery and freezing tissues was kept as minimal as possible due to the fact that the frozen section pathology laboratory was directly adjacent to the operating room. Nevertheless, clinical and pathological workflows do not allow for exact measures and timing of tissue dissection parameters for samples collected during routine surgical interventions. In addition, depending on the surgical technique there is a variable amount of intraoperative tissue ischemia due to surgical ligation of blood vessels. This fact may account for the inability to quantify glycolytic intermediates which have such a high turnover in non-frozen tissues that these are found to be depleted if metabolism is not immediately quenched after disruption of blood flow. To minimize unrelated technical noise related to surgical procedures we have chosen to compare tumor tissue and normal tissue collected during the same surgery.
We have used a metabolomic approach by GC-TOF mass spectrometry in order to gain a broad overview over primary metabolism at limited costs but at high sensitivity and selectivity. In total more than 100 compounds could be identified by chemical structure from as little as 5 mg fresh tissue, which compares favorably to reports using one dimensional 1H-NMR data acquisition. Specifically, we here demonstrate for the first time the efficacy of an automated annotation using a customized database approach. On the one hand, the BinBase database unambiguously identifies chemically or biochemically known compounds that are utilized for pathway mapping. On the other hand, the database also facilitates adding novel and potentially unique metabolic signals that yet are to be structurally identified but that nevertheless were often found to be differentially regulated at high significance levels. These compounds are stored in the database by unique identifiers combining mass spectra and retention index information that enable re-using these database entries for later studies aimed at validating initial biomarkers or at structural identification of these metabolic signals.
In the study presented here, we have focused on using information from identified compounds by developing and applying a new biochemical mapping method, PROFILE. This method maximizes the interpretability of results by facilitating physiological and biochemical understanding of metabolic alterations in carcinoma. Specifically, PROFILE leads to simplified output of results than mapping on single pathway maps from KEGG which would focus on a small number of select metabolites rather than taking into account the relative distances of metabolites across the metabolic network.
As a conclusion, our results show that metabolic signatures as well as individual metabolites can be detected from fresh-frozen tumor tissue of colon cancer and that these alterations can be linked to relevant biochemical pathways. Based on our results, we suggest that metabolomics is a promising approach complementary to transcriptomics and proteomics for analyses of changes in the malignant phenotype. As metabolites constitute the amplified output of a biological system, their quantitative and qualitative analysis will be relevant for tumor biology in different types of investigations. Databases such as the one presented here will enable comparisons of findings across studies and laboratories. Metabolomics can be used for biochemical classification of different tumor types and for comparison of malignant tumors with their corresponding normal tissue. Recently, it has been suggested that therapeutic approaches directed against metabolic abnormalities may be useful in the treatment of malignant tumors [27, 28]. In this context, the metabolic profiling approach described here may be useful to monitor the complex changes in tumor metabolism that may occur under these treatments. Furthermore, analysis of metabolic alterations may be used as a new method for molecular pathology to develop classifiers for therapy response prediction, which may ultimately lead to the identification of new prognostic markers. With the combination of advanced instrumentation, standardized database algorithms and the development of tools for interpretation of data, our study provides a methodological basis for these further investigations.
We would like to thank Ms. Ines Koch, Ms. Lisa Glanz and Ms. Petra Wachs for their excellent technical assistance. Ms. Martina Eickmann is gratefully acknowledged for editorial assistance.
- Jemal A, Siegel R, Ward E, Murray T, Xu J, Thun MJ: Cancer statistics, 2007. CA Cancer J Clin. 2007, 57: 43-66.View ArticlePubMed
- Cousins RJ: Nutritional regulation of gene expression. Am J Med. 1999, 106: 20S-23S.View ArticlePubMed
- Jump DB: Fatty acid regulation of gene transcription. Crit Rev Clin Lab Sci. 2004, 41: 41-78.View ArticlePubMed
- Mariadason JM, Corner GA, Augenlicht LH: Genetic reprogramming in pathways of colonic cell maturation induced by short chain fatty acids: comparison with trichostatin A, sulindac, and curcumin and implications for chemoprevention of colon cancer. Cancer Res. 2000, 60: 4561-72.PubMed
- Gunter MJ, Leitzmann MF: Obesity and colorectal cancer: epidemiology, mechanisms and candidate genes. J Nutr Biochem. 2006, 17: 145-56.View ArticlePubMed
- Mason JB: Diet, folate, and colon cancer. Curr Opin Gastroenterol. 2002, 18: 229-34.View ArticlePubMed
- Bonnet S, Archer SL, Allalunis-Turner J, Haromy A, Beaulieu C, Thompson R, Lee CT, Lopaschuk GD, Puttagunta L, Bonnet S, Harry G, Hashimoto K, Porter CJ, Andrade MA, Thebaud B, Michelakis ED: A mitochondria-K+ channel axis is suppressed in cancer and its normalization promotes apoptosis and inhibits cancer growth. Cancer Cell. 2007, 11: 37-51.View ArticlePubMed
- Fiehn O: Metabolomics – the link between genotypes and phenotypes. Plant Mol Biol. 2002, 48: 155-71.View ArticlePubMed
- Bino RJ, Hall RD, Fiehn O, Kopka J, Saito K, Draper J, Nikolau BJ, Mendes P, Roessner-Tunali U, Beale MH: Potential of metabolomics as a functional genomics tool. Trends Plant Sci. 2004, 9: 418-25.View ArticlePubMed
- Oliver SG: Functional genomics: lessons from yeast. Phil Trans R Soc Lond B. 2002, 357: 17-23.View Article
- Denkert C, Budczies J, Kind T, Weichert W, Tablack P, Sehouli J, Niesporek S, Konsgen D, Dietel M, Fiehn O: Mass spectrometry-based metabolic profiling reveals different metabolite patterns in invasive ovarian carcinomas and ovarian borderline tumors. Cancer Res. 2006, 66: 10795-804.View ArticlePubMed
- Catchpole GS, Beckmann M, Enot DP: Hierarchical metabolomics demonstrates substantial compositional similarity between genetically modified and conventional potato crops. Proc Natl Acad Sci USA. 2005, 102: 14458-62.PubMed CentralView ArticlePubMed
- Fiehn O, Kopka J, Dormann P, Altmann T, Trethewey RN, Willmitzer L: Metabolite profiling for plant functional genomics. Nat Biotechnol. 2000, 18: 1157-61.View ArticlePubMed
- Griffin JL, Shockcor JP: Metabolic profiles of cancer cells. Nature reviews cancer. 2004, 4: 551-561.View ArticlePubMed
- Kyoto Encyclopedia of Genes and Genomes (KEGG). http://www.genome.ad.jp/kegg/
- Scholz M, Fiehn O: SetupX – A Public Study Design Database for Metabolomic Projects. Pacific Symp Biocomp. 2007, 12: 169-180.
- R Development Core Team: R: A language and environment for statistical computing. 2005, ISBN 3-900051-07-0, (R Foundation for Statistical Computing, Vienna, Austria), http://www.R-project.org
- Fan RE, Chen PH, Lin CJ: Working set selection using the second order information for training SVM. Journal of Machine Learning Research. 2005, 6: 1889-1918.
- Venables WN, Ripley BD: Modern Applied Statistics with S. 2002, Springer, Fourth.View Article
- Wessels LF, Reinders MJ, Hart AA, Veenman CJ, Dai H, He YD, van't Veer LJ: protocol for building and evaluating predictors of disease state based on microarray data. Bioinformatics. 2005, 21: 3755-62.View ArticlePubMed
- Usadel B, Nagel A, Thimm O, Redestig H, Blaesing OE, Palacios-Rojas N, Selbig J, Hannemann J, Piques MC, Steinhauser D, Scheible WR, Gibon Y, Morcuende R, Weicht D, Meyer S, Stitt M: Extension of the visualization tool MapMan to allow statistical analysis of arrays, display of corresponding genes, and comparison with known responses. Plant Physio. 2005, 138: 1195-204.View Article
- Kose F, Weckwerth W, Linke T, Fiehn O: Visualizing plant metabolomic correlation networks using clique-metabolite matrices. Bioinformatics. 2001, 17: 1198-208.View ArticlePubMed
- Hlubek F, Brabletz T, Budczies J, Pfeiffer S, Jung A, Kirchner T: Heterogeneous expression of Wnt/beta-catenin target genes within colorectal cancer. Int J Cancer. 121 (9): 1941-8. 2007 Nov 1.
- Warburg O: Über den Stoffwechsel der Tumoren. 1930, London: Constable.
- Okada A, Takehara H, Yoshida K, Nishi M, Miyake H, Kita Y, Komi N: Increased aspartate and glutamate levels in both gastric and colon cancer tissues. Tokushima J Exp Med. 1993, 40: 19-25.PubMed
- Kleinrok Z, Matuszek M, Jesipowicz J, Matuszek B, Opolski A, Radzikowski C: GABA content and GAD activity in colon tumors taken from patients with colon cancer or from xenografted human colon cancer cells growing as s.c. tumors in athymic nu/nu mice. J Physiol Pharmacol. 1998, 49: 303-10.PubMed
- Pelicano H, Martin DS, Xu RH, Huang P: Glycolysis inhibition for anticancer treatment. Oncogene. 2006, 25: 4633-46.View ArticlePubMed
- Pan JG, Mak TW: Metabolic targeting as an anticancer strategy: dawn of a new era?. Sci STKE. 2007, 381: 14-
- Michiels S, Koscielny S, Hill C: Prediction of cancer outcome with microarrays: a multiple random validation strategy. Lancet. 2005, 365: 488-92.View ArticlePubMed
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.