This article has Open Peer Review reports available.
Analyzing the regulation of metabolic pathways in human breast cancer
- Gunnar Schramm†1,
- Eva-Maria Surmann†1,
- Stefan Wiesberg2,
- Marcus Oswald2,
- Gerhard Reinelt2,
- Roland Eils1, 3Email author and
- Rainer König1Email author
© Schramm et al; licensee BioMed Central Ltd. 2010
Received: 2 December 2009
Accepted: 10 September 2010
Published: 10 September 2010
Tumor therapy mainly attacks the metabolism to interfere the tumor's anabolism and signaling of proliferative second messengers. However, the metabolic demands of different cancers are very heterogeneous and depend on their origin of tissue, age, gender and other clinical parameters. We investigated tumor specific regulation in the metabolism of breast cancer.
For this, we mapped gene expression data from microarrays onto the corresponding enzymes and their metabolic reaction network. We used Haar Wavelet transforms on optimally arranged grid representations of metabolic pathways as a pattern recognition method to detect orchestrated regulation of neighboring enzymes in the network. Significant combined expression patterns were used to select metabolic pathways showing shifted regulation of the aggressive tumors.
Besides up-regulation for energy production and nucleotide anabolism, we found an interesting cellular switch in the interplay of biosynthesis of steroids and bile acids. The biosynthesis of steroids was up-regulated for estrogen synthesis which is needed for proliferative signaling in breast cancer. In turn, the decomposition of steroid precursors was blocked by down-regulation of the bile acid pathway.
We applied an intelligent pattern recognition method for analyzing the regulation of metabolism and elucidated substantial regulation of human breast cancer at the interplay of cholesterol biosynthesis and bile acid metabolism pointing to specific breast cancer treatment.
Breast cancer is a prevalent disease and a leading cause of cancer death in women . Worldwide, breast cancer is the second most common type of cancer after lung cancer and the fifth most common cause of cancer death. Breast cancer patients with the same stage of disease can have very different treatment responses and overall outcome. Clinical predictive factors like age, tumor size, lymph node status, histological and pathological grade or hormone-receptor status, often fail to accurately predict clinical outcome, distant metastasis and recurrence of the cancer.
Chemotherapy and hormonal therapy reduces the risk of distant metastases by approximately one third. However, 70-80% of the patients would have survived without it (five years of follow-up) [2–4]. A more accurate means of prognosis and selection of therapy would substantially improve disease-free and overall survival of breast cancer patients . Cancer cells acquire their hallmarks of malignancy through the accumulation of advantageous gene activation and inactivation events over long periods of time . Nevertheless, the molecular basis of breast cancer tumorigenesis remains poorly understood. A long-standing strategy for cancer treatment is to attack basic tumor metabolism by inhibiting nucleotide biosynthesis [7, 8] and DNA production .
Besides this, over the past decade there have been exciting developments in analyzing large scale gene expression profiles. This improved the understanding of the tumors' composition and behavior to develop new targets for therapy . Many studies of gene expression identified expression profiles that are prognostic for patients with breast cancer. However, comparisons of the lists of genes derived from these studies showed that they overlap only slightly due to differences in the patient cohorts, microarray platforms, and mathematical methods of analysis . One strategy to tackle this problem is to map lists of differentially expressed genes on groups of genes with related functions according to the information provided by several databases such as Gene Ontology  and KEGG . Finding enrichments of specific gene sets related to certain phenotypes or cell states yields a functional grouping of differentially expressed genes which can be related to their pathogenic behavior and can lead to more robust results in comparison to the analysis of single genes.
To detect the enrichment of gene sets, commonly, a list of significantly differentially expressed genes is identified and statistical tests applied, such as Fisher's exact test and χ2 test. In a different approach, a gene-specific statistics, known as the "local" statistics, measures the strength of association between the gene expression and the phenotype for each gene. A global statistics for a gene set is then constructed as a function of local statistics for each gene in it. The significance is assessed by permutation tests . A global test for gene sets to associate gene expression with clinical outcome was presented by Goeman and co-workers and enabled determining whether the global expression pattern of a group of genes is significantly related to clinical outcome of interest using a linear regression approach . In another approach, expression levels of all genes in the gene sets are combined and presented as gene specific features. These features are then compared between the treatment and the control groups to identify significantly affected gene sets . In general, these methods test the association of all genes in a gene set with the phenotypes, whereas often only genes in a subset of the gene set are associated with the phenotype. Some of the genes may not belong to the set due to incompleteness or errors in the available data. Additionally, even if all genes in the gene set have apparently the same function, or belong to the same process, it is likely that only a few genes are associated with the phenotype.
To overcome these and other gene-associated problems, transcriptional data was analyzed using topology information of cellular networks. Topological information derived from the metabolic network was connected by calculating Z-scores of highly correlated sub-networks . Chuang and co-workers improved classification of breast cancers with expression patterns of small subnets of a signal transduction network . Substantial new genetic mediators for prostate cancer were found using reverse engineered gene networks in combination with gene expression profiles . Rapaport and co-workers found gene expression patterns of neighboring genes in the network yielding good classification of the profiled samples by calculating Fourier transformations and rejecting high frequency signals . Common gene expression levels of neighboring nodes in a metabolic network were calculated by averaging over all neighbors of a gene, and revealed several interesting regulated pathways for the human immune system . However, these approaches were not developed to detect highly contrasting expression of neighboring genes that undergo a switch-like shift of regulation in a tumor cell. Importantly, especially these switches can be highly relevant to identify potential drug targets that specifically attack the tumor at nodes at which it redirects fluxes in the network to establish parasitic advantages.
In our study, gene expression profiles of breast tumors having an "unfavorable" prognosis were compared to breast tumors with a "favorable" prognosis. We wanted to track how the aggressive (unfavorable) tumors have specifically regulated their metabolism to optimize their oncogenetic fitness, and to elucidate ways to severely perturb this process. For this, we used an approach that detects orchestrated regulation of neighboring enzymes in the metabolic network. We mapped gene expression data onto optimally arranged grid representations of pathways of the metabolic network and applied Haar wavelet transforms onto defined pathways of the network to combine gene expression values from neighboring enzymes. These combined features were tested using a non-parametric test (Wilcoxon) if they could separate samples from different treatments. Metabolic pathways were selected that had features with the most discriminative gene expression patterns. We detected a substantially higher number of significant gene expression patterns in comparison to commonly used enrichment tests. We revealed 19 significant metabolic pathways including increased purine and pyrimidine biosynthesis which were needed for increased mitosis cycles. Furthermore, we found pathways for increased energy metabolism (glycolysis, pyruvate metabolism and fructose/mannose metabolism). Interestingly, we observed the regulation of a possible a cellular switch in the pathway for bile acid biosynthesis redirecting the metabolic flux to the synthesis of steroids while preventing degradation into bile acids.
Results and Discussion
Significantly differentially regulated pathways.
Number of reactions
Down-regulated in unfavorable
Up-regulated in unfavorable
Alanine and aspartate metabolism
Valine, leucine and isoleucine degradation
Fatty acid metabolism
Biosynthesis of Steroids
Glycine, Serine and Threonine metabolism
Bile acid biosynthesis
Pentose phosphate pathway
Inositol phosphate metabolism
Fructose and mannose metabolism
Pyrimidine and purine metabolism
Most up-regulated reactions were identified in the pyrimidine (P = 2.25E - 04) and purine (P = 1.21E - 03) metabolism. These pathways were up-regulated to enable enforced nucleotide biosynthesis for increased cell cycle activity of the aggressive tumors. Nearly all enzymes involved in the biosynthetic pathway for nucleotides were up-regulated, such as enzymes converting substrates to dNTPs and polyribonucleotide nucleotidyltransferases (EC 18.104.22.168), incorporating dNTPs into DNA. Enzymes reversing pyrimidine and purine anabolism, such as enzymes degrading dNTPs, were down-regulated (ECs 22.214.171.124, 126.96.36.199 and 188.8.131.52). Reactions involved in RNA synthesis were partially up-regulated to increase protein biosynthesis (EC 184.108.40.206). Furthermore, reactions which were responsible for synthesizing adenosine were up-regulated (EC 220.127.116.11). Adenosine was shown to be angiogenic, cyto-protective and anti-inflammatory in several tissues and contributed to more aggressive behavior and metastasis of cancer cells .
Pathways for energy supply were significantly up-regulated
We detected significant differential expression patterns in glycolysis (P = 3.85E - 03), pyruvate (P = 1.20E - 02) and fructose/mannose (P = 1.55E - 02) metabolism. They were mostly up-regulated to generate sufficient energy and metabolites for fast-growing cancer cells. Genes of the glycolysis pathway have been found to be over-expressed in a set of 24 cancers, while other pathways showed less consistent up-regulation so far . Glycolysis is increased in cancers and this generates ATP by oxidative phosphorylation [23, 24]. Furthermore, elevated glycolytic flux may support the production of important metabolites that contribute to essential cellular processes like fatty acid synthesis, nucleotide synthesis and to the better protection of oxidative stress [25, 26].
However, the extent of up-regulated genes in glycolysis varies for different cancer types. In our study, all except one differentially regulated reactions were up-regulated, including enzymes such as glyceraldehyde-3-phosphate dehydrogenase (EC 18.104.22.168), phosphor glycerate mutase (EC 22.214.171.124), phosphoglycerate kinase (EC 126.96.36.199), triosephosphate isomerase (EC 188.8.131.52) and fructose-bisphosphate aldolase (EC 184.108.40.206). These enzymes were also described as up-regulated by Altenberg and co-workers . Inhibitors of glycolysis, such as the glucose analog 2-deoxyglucose, which binds and suppresses hexokinase I, and arsenate, that causes arsenolysis in glyceraldehyde-3-phosphate dehydrogenase, as well as 3-bromopyruvate, an inhibitor of hexokinase II, have already been developed to target this metabolic abnormality and could effectively kill cancer cells .
Inhibition of glycolysis was effective even in killing cancer cells with a multidrug resistance (MDR) phenotype. It is known that cells expressing MDR proteins require ATP as their energy source to export drugs out of the cell. Thus, to overcome this drug resistance, depletion of cellular ATP causes the excretion to fail and consequently, cancer cells become more sensitive to anti-cancer therapy .
We found a smaller pattern of up-regulated reactions in the pyruvate metabolism. Significantly up-regulated were: pyruvate kinase (EC 220.127.116.11) which is responsible for ATP production in glycolysis, and the conversion of pyruvate to malate via (S)-Malate:NADP+ oxidoreductase (ECs 18.104.22.168, 22.214.171.124, 126.96.36.199). Aldehyde dehydrogenase (EC 188.8.131.52) was down-regulated to prevent degradation of fatty acids. Oxaloacetate was reported to directly induce cell proliferation by increasing DNA synthesis . Surprisingly, our results did not show any increase in lactate dehydrogenase, which is often up-regulated in cancer cells . Highly expressed reactions in fructose and mannose metabolism were mainly involved in the conversion of D-fructose (EC 184.108.40.206), D-fructose-1P (ECs 220.127.116.11, 18.104.22.168, 22.214.171.124) and D-mannose (ECs 126.96.36.199, 188.8.131.52, 184.108.40.206) into glyceraldehyd-3P, a substrate of glycolysis, leading to increased ATP production.
Biosynthesis of steroids and bile acids
Besides this, there have been intensive cancer investigations on bile acids as they have very heterogeneous effects on tumor cells and carcinogenesis. They are known enhancers for invasiveness of colon cancers [31, 32]. In turn, bile acids can induce apoptosis either specifically (receptor-mediated interactions)  or, in high concentrations, non-specifically (as detergents) through mitochondrial destabilization and oxidative stress . Non-toxic doses of deoxycholic acid (DCA), chenodeoxycholic acid (CDCA) and lithocholic acid (LCA) induced differentiation in promyelocytic leukemia cell lines . 6ECDCA is a synthetic bile acid derivative and can act as a selective FXR ligand to promote differentiation of preadipocyte cell lines . In the blood plasma of postmenopausal women with newly diagnosed breast cancer, elevated concentrations of DCA were detected which may have been released from osteoblasts to induce migration of the cancer . This would explain a down-regulation in breast cancer cells as they may use elevated bile acid levels from the blood for detaching and infiltrating while sustaining their endogenous cholesterol synthesis.
Bile acids are normally predominantly produced by hepatocytes which may explain the elevated levels of bile acids in the blood when elevated breast cancer estrogen is decomposed into bile acids in the liver. In conclusion, this regulation may reflect a cellular switch which redirects the metabolic flux from degradation of cholesterol (via bile acids) into the conversion of steroid hormones, suggesting drug targets in the biosynthesis of steroids, such as mevalonate (diphospho) decarboxylase (EC 220.127.116.11) and mevalonate kinase (EC 18.104.22.168).
Analyzing a second dataset
To verify our findings we analyzed an independent gene expression dataset of breast cancer [36, 37] consisting of 414 samples. The data was downloaded from Array Express (GSE6532, http://www.ebi.ac.uk/microarray-as/ae/) and normalized as described in Methods. We divided the dataset into 250 tumors showing a favorable prognosis and 61 tumors with an unfavorable prognosis following our criteria as described for the other dataset (in Methods). 103 tumors were discarded as they did not fit in any of the prognostic groups. The significance level was again set to P = 0.05. Our analysis revealed 12 pathways with significantly differentially regulated patterns (see Additional file 1 Table A1). Although the overlap of survival predicting genes between different breast cancer studies is often rather small , four pathways from the first analyzed study also showed a significant regulation pattern in this second dataset, i.e. histidine metabolism, valine, leucine and isoleucine degradation, bile acid biosynthesis and galactose metabolism. Interestingly also the citrate cycle showed a significant regulation pattern which compares to the pathways for carbon hydrate metabolism we found in the first analyzed dataset, and C21-steroid hormone metabolism which compares to steroid biosynthesis found in the first analyzed study.
Performing a network based gene expression analysis revealed interesting insights into the metabolism and regulation of aggressive breast tumor cells. Expected differentially regulated pathways in cancer could be confirmed, e.g. the aggressive tumors showed significantly up-regulated pathways for purine and pyrimidine synthesis to maintain elevated proliferation, as well as the up-regulation of glycolysis and pyruvate metabolism to support energy supply for the tumor. The analysis revealed insights into differentially regulated metabolic pathways in breast cancer cells, which may support the induction of proliferation by inositol signal transduction cascades (inositol phosphate metabolism) and steroid hormones (biosynthesis of steroids). An interesting view in breast cancer metabolism was observed in the interplay of biosynthesis of steroids and bile acids, the latter of which was down-regulated possibly to convert cholesterol into proliferative acting steroid hormones. Such a cellular switch would be interesting to compare to other tumors and tissues, also in order to define a tumor specific therapy. Performing Fisher's exact tests as a standard enrichment analysis revealed the lysine biosynthesis pathway to be significantly enriched with differentially expressed genes. The GSEA method yielded pyrimidine biosynthesis. In addition to this pathway, our method yielded eighteen further significant regulation patterns. However, not all pathways detected by such enrichment tests and also our pattern analysis may be relevant for the pathogenesis of the analyzed disease. Some pathways in our results, such as valine, leucine and isoleucine degradation, could not be associated with oncogenesis and may need further examination.
Complex regulated pathways which were relevant to breast tumors with unfavorable prognosis were detected and described in a straightforward manner. The global analysis of network patterns offered a good insight into the regulation of metabolism in breast tumors and may support revealing new potential targets for drug design, specifically in the interplay of the biosynthesis of bile acids and steroids.
We wanted to identify pathways and regions of pathways that showed patterns of differentially regulated reactions. To find such patterns we applied wavelet transforms on the mapped gene expression data. For the application of wavelet transforms and the subsequent feature calculations we needed to arrange the metabolic pathways onto a 2-dimensional grid (Figure 2 illustrates this procedure with a small example pathway).
Preparing the microarray data
Normalized gene expression data was taken from a published study  of breast-cancer samples of 295 women (diagnosis between 1984 and 1995) with age ≤ 53 years and no previous history of cancer, except for non-melanoma skin cancer. The data was downloaded from Rosetta Inpharmatics at http://www.rii.com/publications/2002/nejm.html. The gene expression profiles were derived by using oligonucleotide microarrays from Agilent Technologies http://www.agilent.de. Data on relapse-free survival (defined as the time to the first event) and overall survival were available for all patients. Most of the patients had breast cancer of stages one and two. 165 had received local therapy alone, 20 had received tamoxifen only, 20 had received tamoxifen plus chemotherapy, and 90 had received chemotherapy only. 151 patients were diagnosed as lymph-node-negative and 144 as lymph-node-positive. The tumors were primary invasive breast carcinoma that were less than 5 cm in diameter at pathological examination . To differentiate between tumors with "favorable" and "unfavorable" prognosis, samples were separated according to their "time recurrence" (period of time of relapse after the first event of breast cancer, denoted in years) and "event of metastasis" (if metastasis occurs during this period). A time recurrence above 5 years indicated a "favorable" prognosis (196 patients), whereas samples were classified as samples with unfavorable prognosis if they showed a time recurrence of less than three years and the occurrence of metastasis (54 patients). Samples with ambiguous information were discarded (45 samples) and identified either by a time recurrence of less than three years without the event of metastasis or by a time recurrence above three years but below five years. A Remark: A new ethical approval was not needed for our study as all our analyses based on published microarray datasets (which also applies to the second analyzed dataset).
Assembling the metabolic pathways
Pathways were defined according to curated pathway maps of the KEGG database  (version from February 4th, 2009). Two reactions were neighbors if a metabolite existed that was the product of one and the substrate of the other. We defined reactions as the nodes and metabolites as the edges between them. Pathways without any connected reaction were discarded. This resulted in 99 pathways with 1826 different reactions. Each pathway was represented by its adjacency-matrix. An entry at row a and column b was set to one if there existed a metabolite that was produced by reaction a and consumed by reaction b or vice versa. The sizes of the symmetric adjacency-matrices were between 2 × 2 and 92 × 92 reactions.
Ordering the two-dimensional pathway representation with the grid arrangement method
The right-hand sides lb(G') for the different sub-graph motifs needed to be determined only once as they are independent of G. Furthermore, the motifs were defined in a pre-processing step and could therefore be separated quickly during the optimization process. We considered the sub-graph motifs of star graphs, cliques consisting of up to 10 vertices and odd cycles (2k + 1-cycles) for k = 1,2. Moreover, a certain class of trees with maximum vertex degree Δ(T) ≤ 4 decreased computation time and enhanced separation ability. Furthermore, calculation time was reduced by symmetry breaking constraints eliminating all but a few representative embeddings from each equivalence class of symmetrical embeddings. For this, grid symmetries due to translation, rotation and reflection of the embeddings were considered as well as vertex subsets which inner permutations didn't change the value of the objective function. After solving the optimization problem the values of the variables x vij allowed to construct the embedding of the metabolic pathway on a 2-dimensional grid placing reaction v on the grid point with coordinates (i,j). On the grid representations of the pathways we then applied our feature generation method.
Pattern recognition of gene regulation on the metabolic network
Neighboring enzymes on the 2-dimensional grid representation were grouped by combining their gene expression values with wavelet transforms. These transforms yielded combined expression values ("features") of low pass filters to detect similar expression changes and high pass filters to detect contrasting regulation patterns. The wavelet calculation is described in the following. Each cluster-matrix was divided into 2 × 2 pixeled disjoint sub-sections (e.g. a cluster matrix of size 8 × 8 was divided into 16 sub-sections). Clusters with non-fitting sizes (e.g. 3 × 3, 5 × 5,) were extended with rows and columns of zeros to yield matrices that could be divided into 2 × 2 pixels sub-sections. For each sub-section, all combinations of row-wise and column-wise mean and differences, respectively, were calculated. This yielded 4 combined values for each 2 × 2 pixels sub-section: 1st: mean of the mean of the upper and mean of the lower row, 2nd: difference of the mean of the upper and the mean of the lower row, 3rd: mean of the difference of the upper and the difference of the lower row, and, 4th: difference of the difference of the upper and the difference of the lower row. All four combined values for each 2 × 2 pixels sub-section were stored and applied as features for the classifier. This was done for all sub-sections of the matrix. All 1st combined values (mean of means) were taken for a new matrix and were again grouped into 2 × 2 fractions that were combined in the same manner, yielding again 4 new features for every fraction. This procedure was repeated until no further grouping was possible. The discriminative behavior of all non-trivial features was tested using Wilcoxon rank tests. Pathways were ranked according to their best discriminating features. Features were regarded as significant if they had p-values ≤ 0.05 after correction for multiple testing using the Bonferroni method [40, 41]. Only pathways with more than three significantly, differentially regulated reactions and genes were further investigated to focus on the most relevant features. Two reactions that consisted of exactly the same genes were counted as one reaction.
This work was funded within the BMBF-FORSYS consortium Viroquant (#0313923), the Helmholtz Alliance on Systems Biology and the Nationales Genom-Forschungs-Netz (NGFN+) for the neuroblastoma project, ENGINE.
- Oakman C, Bessi S, Zafarana E, Galardi F, Biganzoli L, Di Leo A: New diagnostics and biological predictors of outcome in early breast cancer. Breast Cancer Research. 2009, 205: 1-11.Google Scholar
- van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, Peterse HL, van der Kooy K, Marton MJ, Witteveen AT, Schreiber GJ, Kerkhoven RM, Roberts C, Linsley PS, Bernards R, Friend SH: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415: 530-536. 10.1038/415530a.View ArticlePubMedGoogle Scholar
- Polychemotherapy for early breast cancer: an overview of the randomised trials. Early Breast Cancer Trialists' Collaborative Group. Lancet. 1998, 352: 930-942. 10.1016/S0140-6736(98)03301-7.Google Scholar
- Tamoxifen for early breast cancer: an overview of the randomised trials. Early Breast Cancer Trialists' Collaborative Group. Lancet. 1998, 351: 1451-1467. 10.1016/S0140-6736(97)11423-4.Google Scholar
- Van de Vijver MJ, He YD, van't Veer LJ, Dai H, Hart AAM, Voskuil DW, Schreiber GJ, Peterse JL, Roberts C, Marton MJ, Parrish M, Atsma D, Witteveen A, Glas A, Delahaye L, van der Velde T, Bartelink H, Rodenhuis S, Rutgers ET, Friend SH, Bernards R: A Gene-Expression Signature as a Predictor of Survival in Breast Cancer. The New England Journal of Medicine. 2002, 347: 1999-2009. 10.1056/NEJMoa021967.View ArticlePubMedGoogle Scholar
- Fan C, Oh DS, Wessels L, Weigelt B, Nuyten DSA, Nobel AB, van't Veer LJ, Perou CM: Concordance among Gene-Expression-Based Predictors for Breast Cancer. The New England Journal of Medicine. 2006, 355: 560-569. 10.1056/NEJMoa052933.View ArticlePubMedGoogle Scholar
- Chen L, Pankiewicz KW: Recent development of IMP dehydrogenase inhibitors for the treatment of cancer. Curr Opim Drug Discov Devel. 2007, 10: 403-412.Google Scholar
- Lui MS: Biochemical pharmacology of acivicin in rat hepatoma cells. Biochem Pharamcol. 1982, 31: 3469-3473. 10.1016/0006-2952(82)90628-1.View ArticleGoogle Scholar
- Pedersen-Bjergaard J: Risk of acute nonlymphocytic leukemia and preleukemia in patients treated with cyclophosphamide for non-Hodgkin's lymphomas. Comparison with results obtained in patients treated for Hodgkin's disease and ovarian carcinoma with other alkylating agents. Ann Intern Med. 1985, 103: 195-200.View ArticlePubMedGoogle Scholar
- Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, Richter J, Rubin GM, Blake JA, Bult C, Dolan M, Drabkin H, Eppig JT, Hill DP, Ni L, Ringwald M, Balakrishnan R, Cherry JM, Christie KR, Costanzo MC, Dwight SS, Engel S, Fisk DG, Hirschman JE, Hong EL, Nash RS, Sethuraman A, Theesfeld CL, Botstein D, Dolinski K, Feierbach B, Berardini T, Mundodi S, Rhee SY, Apweiler R, Barrell D, Camon E, Dimmer E, Lee V, Chisholm R, Gaudet P, Kibbe W, Kishore R, Schwarz EM, Sternberg P, Gwinn M, Hannick L, Wortman J, Berriman M, Wood V, de la Cruz N, Tonellato P, Jaiswal P, Seigfried T, White R: The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 2004, 32: D258-261. 10.1093/nar/gkh066.View ArticlePubMedGoogle Scholar
- Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36: D480-484. 10.1093/nar/gkm882.View ArticlePubMedGoogle Scholar
- Pitman EJG: Significance tests which may be applied to samples from any population. Royal Statistical Society Supplement. 1937, 4: 119-130. 10.2307/2984124. 225-232View ArticleGoogle Scholar
- Goeman JJ, van de Geer SA, de Kort F, van Houwelingen HC: A global test for groups of genes: testing association with a clinical outcome. Bioinformatics. 2004, 20: 93-99. 10.1093/bioinformatics/btg382.View ArticlePubMedGoogle Scholar
- Yan X, Sun F: Testing gene set enrichments for subset of genes: Sub-GSE. BMC Bioinformatics. 2008, 9: 1-15. 10.1186/1471-2105-9-1.View ArticleGoogle Scholar
- Patil KR, Nielsen J: Uncovering transcriptional regulation of metabolism by using metabolic network topology. Proc Natl Acad Sci USA. 2005, 102: 2685-2689. 10.1073/pnas.0406811102.View ArticlePubMedPubMed CentralGoogle Scholar
- Chuang HY, Lee E, Liu YT, Lee D, Ideker T: Network-based classification of breast cancer metastasis. Mol Syst Biol. 2007, 3: 140-10.1038/msb4100180.View ArticlePubMedPubMed CentralGoogle Scholar
- Ergun A, Lawrence CA, Kohanski MA, Brennan TA, Collins JJ: A network biology approach to prostate cancer. Mol Syst Biol. 2007, 3: 82-10.1038/msb4100125.View ArticlePubMedPubMed CentralGoogle Scholar
- Rapaport F, Zinovyev A, Dutreix M, Barillot E, Vert JP: Classification of microarray data using gene networks. BMC Bioinformatics. 2007, 8: 35-10.1186/1471-2105-8-35.View ArticlePubMedPubMed CentralGoogle Scholar
- Nacu S, Critchley-Thorne R, Lee P, Holmes S: Gene Expression Network Analysis, and Applications to Immunology. Bioinformatics. 2007Google Scholar
- Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, Houstis N, Daly MJ, Patterson N, Mesirov JP, Golub TR, Tamayo P, Spiegelman B, Lander ES, Hirschhorn JN, Altshuler D, Groop LC: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34: 267-273. 10.1038/ng1180.View ArticlePubMedGoogle Scholar
- Spychala J: Regulation and function of ecto-5'-nucleotidase and adenosine in cancer. 39 th Meeting of the Polish Biochemical society. 2003, 185.Google Scholar
- Altenberg B, Greulich KO: Genes of glycolysis are ubiquitously overexpressed in 24 cancer classes. Genomics. 2004, 84: 1014-1020. 10.1016/j.ygeno.2004.08.010.View ArticlePubMedGoogle Scholar
- Guppy M, Greiner E, Brand K: The role of the Crabtree effect and an endogenous fuel in the energy metabolism of resting and proliferating thymocytes. Eur J Biochem. 1993, 212: 95-99. 10.1111/j.1432-1033.1993.tb17637.x.View ArticlePubMedGoogle Scholar
- Pfeiffer T, Schuster S, Bonhoeffer S: Cooperation and competition in the evolution of ATP-producing pathways. Science. 2001, 292: 504-507. 10.1126/science.1058079.View ArticlePubMedGoogle Scholar
- Feron O: Pyruvate into lactate and back: from the Warburg effect to symbiotic energy fuel exchange in cancer cells. Radiother Oncol. 2009, 92: 329-333. 10.1016/j.radonc.2009.06.025.View ArticlePubMedGoogle Scholar
- Vander Heiden MG, Cantley LC, Thompson CB: Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009, 324: 1029-1033. 10.1126/science.1160809.View ArticlePubMedPubMed CentralGoogle Scholar
- Xu R, Pelicano H, Zhou Y, Carew JS, Feng L, Bhalla KN, Keating MJ, Huang P: Inhibition of glycolysis in cancer cells: a novel strategy to overcome drug resistance associated with mitochondrial respiratory defect and hypoxia. Cancer Research. 2005, 65: 613-621. 10.1158/0008-5472.CAN-04-4313.View ArticlePubMedGoogle Scholar
- Li Y: Oxaloacetate induces DNA synthesis and Mitosis in primary cultured Rat Hepatocytes in the basence of EGF. Biochemical and Biophysical Research communications. 1993, 193: 1339-1346. 10.1006/bbrc.1993.1772.View ArticlePubMedGoogle Scholar
- Suzuki T, Miki Y, Nakamura Y, Moriya T, Ito K, Ohuchi N, Sasano H: Sex steroid-producing enzymes in human breast cancer. Endocrine-Related Cancer. 2005, 12: 701-720. 10.1677/erc.1.00834.View ArticlePubMedGoogle Scholar
- Zimber A, Gespach C: Bile acids and derivatives, their nuclear receptors FXR, PXR and ligands: role in health and disease and their therapeutic potential. Anticancer Agents Med Chem. 2008, 8: 540-563.View ArticlePubMedGoogle Scholar
- Debruyne PR, Bruyneel EA, Li X, Zimber A, Gespach C, Mareel MM: The role of bile acids in carcinogenesis. Mutat Res. 2001, 480-481: 359-369.View ArticlePubMedGoogle Scholar
- Pai R, Tarnawski AS, Tran T: Deoxycholic acid activates beta-catenin signaling pathway and increases colon cell cancer growth and invasiveness. Mol Biol Cell. 2004, 15: 2156-2163. 10.1091/mbc.E03-12-0894.View ArticlePubMedPubMed CentralGoogle Scholar
- Garewal H, Bernstein H, Bernstein C, Sampliner R, Payne C: Reduced bile acid-induced apoptosis in "normal" colorectal mucosa: a potential biological marker for cancer risk. Cancer Res. 1996, 56: 1480-1483.PubMedGoogle Scholar
- Katona BW, Anant S, Covey DF, Stenson WF: Characterization of Enantiomeric Bile Acid-induced Apoptosis in Colon Cancer Cell Lines. The journal of biological chemistry. 2009, 284: 3354-3364. 10.1074/jbc.M805804200.View ArticlePubMedPubMed CentralGoogle Scholar
- Zimber A, Chedeville A, Gespach C, Abita JP: Inhibition of proliferation and induction of monocytic differentiation on HL60 human promyelocytic leukemia cells treated with bile acids in vitro. Int J Cancer. 1994, 59: 71-77. 10.1002/ijc.2910590115.View ArticlePubMedGoogle Scholar
- Loi S, Haibe-Kains B, Desmedt C, Lallemand F, Tutt AM, Gillet C, Ellis P, Harris A, Bergh J, Foekens JA, Klijn JG, Larsimont D, Buyse M, Bontempi G, Delorenzi M, Piccart MJ, Sotiriou C: Definition of clinically distinct molecular subtypes in estrogen receptor-positive breast carcinomas through genomic grade. J Clin Oncol. 2007, 25: 1239-1246. 10.1200/JCO.2006.07.1522.View ArticlePubMedGoogle Scholar
- Loi S, Haibe-Kains B, Desmedt C, Wirapati P, Lallemand F, Tutt AM, Gillet C, Ellis P, Ryder K, Reid JF, Daidone MG, Pierotti MA, Berns EM, Jansen MP, Foekens JA, Delorenzi M, Bontempi G, Piccart MJ, Sotiriou C: Predicting prognosis using molecular profiling in estrogen receptor-positive breast cancer treated with tamoxifen. BMC Genomics. 2008, 9: 239-10.1186/1471-2164-9-239.View ArticlePubMedPubMed CentralGoogle Scholar
- Ein-Dor L, Kela I, Getz G, Givol D, Domany E: Outcome signature genes in breast cancer: is there a unique set?. Bioinformatics. 2005, 21: 171-178. 10.1093/bioinformatics/bth469.View ArticlePubMedGoogle Scholar
- Nemhauser G, Wolsey L: Integer and Combinatorial Optimization. 1999, John Wiley & Sons, HobokenGoogle Scholar
- Bonferroni CE: Il calcolo delle assicurazioni su gruppi di test. Studi in Onore del Professore Salvatore Ortu Carboni. Rome, Italy. 1935, 13-60.Google Scholar
- Gordi T, Khamis H: Simple solution to a common statistical problem: interpreting multiple tests. Clin Ther. 2004, 26: 780-786. 10.1016/S0149-2918(04)90078-1.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/3/39/prepub
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.