The functional cancer map: A systems-level synopsis of genetic deregulation in cancer
© Krupp et al; licensee BioMed Central Ltd. 2011
Received: 12 November 2010
Accepted: 30 June 2011
Published: 30 June 2011
Cancer cells are characterized by massive dysegulation of physiological cell functions with considerable disruption of transcriptional regulation. Genome-wide transcriptome profiling can be utilized for early detection and molecular classification of cancers. Accurate discrimination of functionally different tumor types may help to guide selection of targeted therapy in translational research. Concise grouping of tumor types in cancer maps according to their molecular profile may further be helpful for the development of new therapeutic modalities or open new avenues for already established therapies.
Complete available human tumor data of the Stanford Microarray Database was downloaded and filtered for relevance, adequacy and reliability. A total of 649 tumor samples from more than 1400 experiments and 58 different tissues were analyzed. Next, a method to score deregulation of KEGG pathway maps in different tumor entities was established, which was then used to convert hundreds of gene expression profiles into corresponding tumor-specific pathway activity profiles. Based on the latter, we defined a measure for functional similarity between tumor entities, which yielded to phylogeny of tumors.
We provide a comprehensive, easy-to-interpret functional cancer map that characterizes tumor types with respect to their biological and functional behavior. Consistently, multiple pathways commonly associated with tumor progression were revealed as common features in the majority of the tumors. However, several pathways previously not linked to carcinogenesis were identified in multiple cancers suggesting an essential role of these pathways in cancer biology. Among these pathways were 'ECM-receptor interaction', 'Complement and Coagulation cascades', and 'PPAR signaling pathway'.
The functional cancer map provides a systematic view on molecular similarities across different cancers by comparing tumors on the level of pathway activity. This work resulted in identification of novel superimposed functional pathways potentially linked to cancer biology. Therefore, our work may serve as a starting point for rationalizing combination of tumor therapeutics as well as for expanding the application of well-established targeted tumor therapies.
Keywordscancer systems biology prognostic marker microarray bioinformatics
Cancer is one of the leading causes of death . If segregated by age, cancer, has even surpassed the mortality of heart diseases and has become the deadliest disease in the world, with 12.7 million new cases and 7.6 million deaths worldwide . Furthermore, due to the steadily increase in aging and expansion of the world's population an ongoing increase also in the burden of cancer is predicted for the next decades . Nevertheless, since the year 1990, cancer death rate continuously declined despite increasing numbers of new cancer cases , indicating the improvements in diagnostics and development of novel more effective therapies. Therein gene expression profiling made a considerable contribution. This technology can be used to classify diverse types of tumors [3, 4] and predict outcome of patients [5–7] as well as response to chemotherapies [8, 9]. Several studies utilized global transcriptome and genome profiling to identify phenotypically similar cancers [10–12]. Segal et al. successfully introduced a module map for cancers , consisting of gene sets that act in concert to carry out specific function. More recently, a meta-analyses tool for different cancer tissues was established . All of these studies revealed associations between molecular profiles and biological behavior of different tumors, without a detailed integration of regulated genes and signaling pathways. Utilization of these functional associations to place global transcriptome data in a biological context might significantly improve the limited knowledge about common molecular changes across different tumor types.
Here, we employed an extensive functional genomics approach to analyze gene expression profiles of the complete human cancer related microarray datasets from the Stanford Microarray Database . Our aim was to develop a functional cancer map focusing on signaling and metabolic events across multiple cancer types to provide a deeper understanding of common and distinct molecular mechanisms leading to carcinogenesis. This novel approach provides a systematic and unbiased view on cancer development by investigating signaling pathways and functional networks in an easy-to-interpret way. In addition, we demonstrate tumor phylogeny based on semantic aggregated profiles by identifying tumor mechanisms of individual tumor entities from publically available gene expression data. Our analysis revealed important signaling pathways commonly dysregulated in several tumor entities. These results could serve as a starting point for the efficient combination of tumor therapeutics as well as for expanding the application area of well-established tumor therapies.
Gene expression data
Gene expression data were obtained from microarray experiments, downloaded from the Stanford Microarray Database . Ignoring cell line and time course experiments, we collected a dataset of 1402 human two-color cDNA microarrays corresponding to 58 discrete tissues types and about 15000 individual genes. Screening for tumor relevant information in our dataset resulted in 836 tumor arrays, which could be assigned to 33 tumor tissue types defined as tumor entities in our analysis. Additionally 16 tumor classes were allocated to tumor entities representing the tumor differentiation and developmental origin of the tissue, e.g. tumor entities 'Prostate Tumor, Grade T2b', 'Prostate Tumor, Grade T3a' and 'Prostate Tumor, Grade 3b' were assigned to the tumor class 'Prostate'. Downloaded microarray experiments had been analyzed against diverse human controls . In order to avoid divergence due to diverse control data sets, we filtered all tumor classes for common applied controls used in the respective microarray experiments. We always used those data with the largest number of individual datasets with the same one control per tumor class for further evaluation. All other data were discarded. Finally our dataset was reduced to 649 microarrays belonging to 28 tumor entities respectively 16 tumor classes, which were analyzed with the R Project statistical software  to identify differentially expressed genes as follows:
Initially the raw data were loaded into R, invalid cDNA probes removed (e.g. flag 'contaminated' or 'failed') and a background correction was applied. After the calculation of the log2 fold changes relative to the control RNA the quantile normalization was applied to each experiment. In all these initial steps we made use of the 'limma' R package . Next, we calculated the variance as well as mean for each gene within the corresponding tumor entities. Moreover, the pooled variance for each gene was calculated to identify global outliers. Pooled variance was defined as weighted mean of the variances within each tumor entity and weighted by the relative sizes of the respective entity. Genes (426 in total) with a pooled variance of all genes greater than the third interquartile range of the pooled variance (3iQR = 0.038) were defined as outliers and excluded from further analysis, to ensure the comparability and correct technical bias due to heterogeneity of the dataset. Additionally, tumor entity specific outliers were eliminated by considering those genes whose variance were greater than the 3iQR of the corresponding tumor entity; depending on tumor entity between 364 and 1053 genes were excluded from the dataset.
Metabolic reaction and signaling event information were acquired from the PATHWAY database of the Kyoto Encyclopedia of Genes and Genomes (KEGG) . These pathway maps and pathway modules are categorized into metabolism, genetic information processing, environmental information processing, cellular process, organismal systems and human diseases in several organisms. Our analysis was based on the KEGG v47 and comprises 210 human pathway maps in XML file format corresponding to 6 categories and 32 subcategories. The data are publically available and can be accessed from the FTP ftp://ftp.genome.jp/pub/. The KEGG pathway data were parsed by a Perl  script using the XML::simple module, resulting in 10001 gene-pathway relations belonging to 208 pathway maps and 4186 unique genes (figure 1b).
Integration of gene expression and functional data
For each tumor entity - KEGG pathway pair, we tested whether the set of genes that are differentially expressed in the given tumor entity is enriched in the set of genes belonging to the given KEGG pathway. Testing was performed using an one-sided exact Fisher test. This resulted in a 28 tumor entities × 208 KEGG pathways matrix containing the p-values of each test.
Using the web-based tool PubMatrix  for literature mining of PubMed , we analyzed our significantly enriched KEGG pathway maps against an assortment of tumor associated terms to identify pathways previously mentioned in context to cancer. Tumor associated keywords used in the PubMatrix searches included "Neoplasm", "Tumors", "Tumor", "Benign Neoplasms", "Neoplasms, Benign", "Benign Neoplasm", "Neoplasm, Benign", "Cancer", "Cancers" and "neoplasms[mesh]". PubMatrix is limited to 10000 results for each search and modifier terms. All search results can be found in the additional files (Additional files 5).
Calculation of tumor phylogeny
In a second step, we used the binarized functional cancer map (1 respectively 0 representing significant resp. non-significant entries) to calculate a similarity score between tumor entities. The similarity between two tumor entities T1 and T2, was quantified by calculating the Cohen's Kappa Coefficient  of corresponding columns in the binarized functional cancer map. This coefficient is a statistical measure of agreement for qualitative items.
Gene expression data
The tumor sample distribution of the 649 microarrays with corresponding controls among the tumor classes and tumor entities used in our analysis
Number of Samples
Acute Lymphoblastic Leukemia Peripheral Blood, Grade Pre-B
Acute Lymphoblastic Leukemia Bone Marrow, Grade Common
Acute Lymphoblastic Leukemia Bone Marrow, Grade Pre-B
Extraskeletal Myxoid Chondrosarcoma
Breast Carcinoma Invasive Ductal Carcinoma, Grade 2
Breast Carcinoma Invasive Ductal Carcinoma, Grade 3
Breast Carcinoma Invasive Lobular Carcinoma, Grade 2
Human Hepatocellular Carcinoma
Solitary Fibrous Tumor
Adenocarcinoma, Pancreas, Grade II
Adenocarcinoma, Pancreas, Grade III
Prostate Tumor, Grade T2b
Prostate Tumor, Grade T3a
Prostate Tumor, Grade T3b
Conventional Renal Cell Carcinoma, Grade 1
Conventional Renal Cell Carcinoma, Grade 2
Conventional Renal Cell Carcinoma, Grade 3
Conventional Renal Cell Carcinoma, Grade 4
Stomach Cancer Before Treatment
Stomach Cancer, Placebo After 1 Year
Exposing functional mechanisms
Focusing on deregulated genetic signaling events across multiple tumors as well as metabolic events in tumor progression, we analyzed the differentially expressed genes in context of the functional categories as well as signaling and metabolic pathway maps provided by the Kyoto Encyclopedia of Genes and Genomes (KEGG) . After assigning the KEGG data to gene expression data, a total of 141 KEGG pathway maps were significantly enriched and revealed four well conserved subcategories among the tumor entities: 'Cell Growth and Death', 'Signaling Molecules and Interactions', 'Immune System' and 'Immune System Diseases'. The most highly conserved KEGG pathway was 'Cell adhesion molecules (CAMs)' (p value ranging between 5.58e-13 and 0.01) which is significantly enriched in 26 tumor entities, followed by the KEGG pathway map 'Asthma' (1.37e-12 - 0.01), 'Cell cycle' (1.42e-09 - 0.02) and 'Hematopoietic cell lineage' (1.07e-7 - 0.04) in 25 tumor entities. 'Antigen processing and presentation' (2.59e10 - 0.02), 'ECM-receptor interaction' (7.11e-9 - 0.04), 'Complement and coagulation cascades' (1.66e-15 - 0.03) and 'Graft-versus-host disease' (9.37e-10 - 0.04) KEGG pathway maps were found to be significantly enriched in 24 tumor entities. As a result, we generated a functional cancer map to visualize this large amount of multiple functional tumor data into a clearly structured compendium (figure 2).
To evaluate the significantly enriched KEGG pathway maps in context of published literature to tumor progression already described and mostly unknown we utilized the web-based tool PubMatrix . The 'Cell cycle' pathway was the most cited pathway with focus on tumor progression displaying more than 1000000 search results per modifier terms in total. Furthermore, 'Melanoma' and 'Acute myeloid leukemia' with 718282 respectively 493765 tumor associated publications were ranked at the top. Taking into count highly enriched KEGG pathway maps, among the KEGG pathway map 'Cell cycle' in KEGG class 'Cellular Processes', the 'Nucleotide excision repair' (174594) in 'Genetic Information Processing' as well as 'Wnt signaling pathway' (21699) in 'Environmental Information Processing' were highly linked to tumor publications. In contrast, highly enriched KEGG pathway maps 'ECM-receptor interaction' (4), the 'Complement and Coagulation cascades' (35) and the 'PPAR signaling pathway' (238) were mostly undescribed. The complete search results can be found in the additional files (Additional file 5).
Exposing tumor phylogeny
Further use of the functional cancer map was made to calculate tumor phylogeny by laying down new starting points for the efficient combination of tumor therapeutics as well as for expanding the application area of well-established tumor therapies. We postulated that tumor entities with comparable significantly enriched KEGG pathway maps demonstrate comparable activation of biochemical processes and, thus, biochemical behavior. As a measure of closer or distant relationship between two tumor entities we applied the Cohen's Kappa Coefficient  to the binarized functional cancer map data. Analyzing phylogeny between diverse tumor entities, highest concordance was found between tumor entities within the same tissue class. The tumor tissue 'Conventional Renal Cell Carcinoma, Grade 2' and 'Conventional Renal Cell Carcinoma, Grade 3' showed the highest concordance κ = 0.77, followed by 'Breast Carcinoma' and 'Breast Carcinoma Invasive Ductal Carcinoma, Grade 2' (κ = 0.74) and 'Breast Carcinoma Invasive Ductal Carcinoma, Grade 2' and 'Breast Carcinoma Invasive Ductal Carcinoma, Grade 2' (κ = 0.70). Examining the tumor entities belonging to different tumor classes, 'Breast Carcinoma Invasive Ductal Carcinoma, Grade 2' and 'Adenocarcinoma, Pancreas, Grade II' (κ = 0.69) displayed the most functional similarity. Also, high similarity in biological behavior was found between 'Conventional Renal Cell Carcinoma, Grade 2' and 'Adenocarcinoma, Pancreas, Grade II' (κ = 0.63) as well as 'Acute Lymphoblastic Leukemia Peripheral Blood, Grade Pre-B' and 'Acute Lymphoblastic Leukemia Bone Marrow, Grade Common' (κ = 0.61) (figure 3).
Despite comprehensive gene expression sets available on multiple human tumor entities, only few studies have attempted to assess and compare biological behavior of the tumors on the basis of functional microarray analysis and evaluation of commonly deregulated and conserved genetic signaling and metabolic events. Moreover, previous analyses focused on individual mechanisms or did not compose tumor-specific gene expression information into detailed, structured maps, impeding a systems view on cancer. However, in particular this kind of information could potentially guide medical research in the context of new applications and novel combinations of existing chemotherapeutics. Therefore, we here present gene expression information of tumors according to pathway activity profiles deemed to better reflect biological behavior. Our analysis of multiple tumor entities using these activity profiles resulted in a clearly structured functional cancer map (figure 2). Note: Due to experimental variation, wrong annotation, and simplifying model assumptions, our pathway signature almost surely fails to identify all active cellular processes in the respective experiment/cancer type . But this is not the goal of the present work, our goal is to identify and distinguish relevant features of multiple tumor tissues; we demonstrate that the proposed analysis is sensitive enough for this purpose.
Exposing functional mechanisms
As a proof of principle and supporting our approach, we found well documented biological events such as cell cycle and apoptosis to be highly conserved among the diverse tumor entities. These findings confirmed the reliability of the generated functional cancer map. Deregulated balance between proliferation and cell death represents a pro-tumorigenic mechansim in human carcinogenesis . Both activation of proliferative signals and an inhibition of death process, leading to survival and consecutively to proliferation, may contribute essential steps in the progression to cancer in affected cells. Among highly enriched pathways in multiple tumors were 'Cell cycle', 'Signaling Molecules and Interactions', and 'Cell Growth and Death' summarizing genes involved in cell proliferation and apoptosis, demonstrating that key tumorigenic features of tumor development were identified by our analysis. Furthermore multiple functional pathways within the categories, 'Immune System' and 'Immune System Diseases' summarizing genes involved in immune response and immunological changes were also highly enriched in almost all tumor entities. Changes in these processes have been extensively documented to be involved in cancer development and growth [34–38]. Whenever normal body cells turn into tumor cells, various antigens on their surface change and were constantly secreted into the circulatory system. The so called tumor antigens were recognized by the immune system and cytotoxic T cells, natural killer cells and macrophages were advised to catch and eliminated the cell that undergo malignant transformation. A disturbed immune surveillance may lead to an increased tumor growth [34, 38].
However, besides extensively studied and repeatedly documented signaling pathways highlighted in our analysis several pathways were found to be commonly enriched among differentially regulated genes which have not been clearly attributed to tumorigenesis. Thus, our work identified key roles of additional signaling pathways in cancer biology. Deregulations or mutations in the molecules involved in the 'Cell Cycle' and 'ECM-receptor interaction' pathway, preserving the physical link between extracellular matrix and actin cytoskeleton, are also closely connected to tumor progression . Besides the 'Cell Cycle' pathway, the 'ECM-receptor interaction' pathway interacts with the 'Cell adhesion molecules (CAMs)'. This balance between cell adhesion and extra cellular molecules is essential for normal cell survival, imbalance between those pathways leads to detaching cells from the extra cellular matrix and therefore enhance metastasis [40, 41]. However, the ECM-receptor interaction and mechanisms of focal adhesion were not in the focus of biomedical research on cancer mechanisms. Over the past decades, the majority of cancer related genetic studies were focussed on genes involved in transcription and cell cycle control. This is reflected in a tenfold higher number of abstracts containing the search terms 'cancer and cell cycle' than 'cancer and extracellular matrix or microenvironment' (122 594 vs. 20 287 articles) when entered in a PubMed query. This may be due to the fact that the plethora of biological functions of ECM proteins has been underestimated for a long time. For a long time functional properties of the extracellular microenvironment and matrix was believed to be restricted to form the "glue" between cells, responsible for simply to maintain shape and coherence of tissues and organs as well as a reservoir of body fluids. In the last decade the microenvironment gained increasing attention for many other processes including tumorigenesis  e.g. the breast microenvironment plays an important and complex role in hepatocarciongenesis. Therefore, cancers are currently regarded as heterogeneous multicellular entities containing cells of multiple lineages who interact with one another, the extracellular matrix (ECM), and soluble molecules in their vicinity are dynamic. Overall, these interactions favour cancer cell proliferation, movement, differentiation, and ECM metabolism, while simultaneously restricting cell death, stationary polarized growth, and ECM stability . However, the microenvironment is also essential to key processes such as angiogenesis. Another important factor is the appreciation that cell migration, matrix and tissue remodeling are not unique properties of cancer cell growths but instead are tightly regulated programs normally utilized during development and in adult tissues responding to acute tissue injury . Reciprocal interactions between responding "normal" cells, their mediators, structural components of the ECM, and genetically altered neoplastic cells regulate all aspects of tumourigenicity. An additional role was suggested for complex interaction of the extra cellular tumour stroma with (mesenchymal) stem cells involved in cancer growth . Thus, demonstrating that functional pathways involved in ECM formation and loosening of focal adhesions are differentially regulated in almost all cancer entities, our results clearly advocate for enhanced efforts in unravelling molecular mechanisms of cancer development influenced by these pathways.
Almost no studies were published on the functional pathway of 'Complement and Coagulation cascades'. Few studies demonstrated a potential role of individual genes within this pathway in cancer development. Some members of the complement system can be modulated by the antifibrinolytic protein, activated thrombin-activatable fibrinolysis inhibitor (TAFIa), by inactivating the anaphylatoxins C3a, C4a and C5a and dysregulation of those genes has been implicated in tumor growth [45, 46]. In addition a link between the haemostatic system and both tumor stroma and metastasis has been noticed in animal and in vitro models. The relationship between cancer and blood coagulation has been described on a diverse pathophysiological level. Cancer development and growth was shown to go along with hypercoagulable state and predisposes to activation of platelets and thrombosis, blood coagulation and fibrinolysis interfering with tumor cell biology, tumor growth, angiogenesis and metastatic process. However, the detailed molecular mechanisms of anti-coagulants such as warfarin on cancer remain elusive [47, 48]. Given the strong enrichment of differentially regulated genes linked to the functional pathway of complement and coagulation, our data suggested that it may be worth to further elucidate exact underlying molecular mechanisms as they may hold key regulatory mechanisms in cancer development and growth.
Another promising candidate for further evaluation was the highly conserved PPAR signaling pathway. It is known that the peroxisome proliferation-activated receptor (PPAR) signaling pathway regulates a multitude of genes important for diverse cellular functions including cell proliferation, cell differentiation, immune response and apoptosis. Once activated, PPARgamma will preferentially bind with retinoid X receptor alpha and signal antiproliferative, antiangiogenic, and prodifferentiation pathways in several tissue types. Since ligands and other agents influencing this pathway have revealed anticancer effects in a variety of human cancers [49, 50] efforts to target this signaling pathway should certainly be intensified given the broad range of tumors with deregulated PPAR signaling.
Finally, we noticed an enrichment of genes invovlved in the functional pathways "Diabetes type I" and "gluycolysis/gluconeogesis". Besides our data several clinical reports support a synergism between Diabetes mellitus and cancer development. Konishi I et al. reported a correlation between a diabetic pattern of the 75 g oral glucose tolerance test OGTT and hepatocarcinogenesis and considered it a significant risk factor for hepatocellular carcinoma (HCC) in patients with hepatitis C virus . Also Huo et al. found that Diabetes mellitus independently predicts decreased survival in HCC patients . Finally, El-Serag et al. reported that Diabetes mellitus is associated with an increased risk for HCC . Interestingly the PPARgamma activators are utilized for the treatment of diabetes mellitus pointing towards a congruent signaling in of PPAR pathway and diabetes signatures .
In summary, there is accumulating evidence for a synergy between diabetes mellitus and an increased risk for and decreased survival of patients with HCC. Our data suggests that important interactions may be identified by analyzing the enriched functional pathways of "Diabetes type I" and "gluycolysis/gluconeogesis" previously not extensively studied in context of cancer.
Exposing tumor phylogeny
We introduced a measure of profile similarity by which tumor entities were structured based on semantic aggregated profiles, which means significant, individual tumor mechanisms calculated out of gene expression data, into a dendrogram (figure 2). In order to get a measurement for genetic distance and thus to estimate the affinity between the diverse tumor entities, we applied the Cohen's Kappa Coefficient  to our binarized functional cancer map and clustered the distant relation matrix (figure 3). Predominant clustering of the tumor entities to their corresponding tumor class is indication for the robustness of the tumor phylogeny.
Notably, examining downregulated genes and clustering of the tumor entities by means of KEGG pathways, we found our sample to almost perfectly cluster according to the germ layer of the tissue of tumor origin. However, this finding needs to be further analyzed and applied to robust statistical testing [54, 55].
Multiple tumor entities have been analyzed for signaling and metabolic events in tumor progression and composed into a functional cancer map. The resulting detailed, structured map should lead to a more systematic view on common and individual tumor mechanisms. This work lead to the identification of novel superimposed functional pathways closely linked to cancer biology. Among these pathways were 'ECM-receptor interaction', 'Complement and Coagulation cascades', and 'PPAR signaling pathway'. Furthermore tumor phylogeny was calculated with respect to their biological behavior. Both, the functional cancer map paired with the indication for tumor phylogeny could serve as a starting point for the efficient combination of targeted tumor therapeutics as well as for expanding the application area of well-established tumor therapies.
- Jemal A, Siegel R, Ward E, Hao Y, Xu J, Thun MJ: Cancer statistics, 2009. CA Cancer J Clin. 2009, 59 (4): 225-249. 10.3322/caac.20006.View ArticlePubMedGoogle Scholar
- Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin DM: Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int J Cancer. 2010Google Scholar
- Lee JS, Heo J, Libbrecht L, Chu IS, Kaposi-Novak P, Calvisi DF, Mikaelyan A, Roberts LR, Demetris AJ, Sun Z, et al: A novel prognostic subtype of human hepatocellular carcinoma derived from hepatic progenitor cells. Nat Med. 2006, 12 (4): 410-416. 10.1038/nm1377.View ArticlePubMedGoogle Scholar
- Alizadeh AA, Eisen MB, Davis RE, Ma C, Lossos IS, Rosenwald A, Boldrick JC, Sabet H, Tran T, Yu X, et al: Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling. Nature. 2000, 403 (6769): 503-511. 10.1038/35000501.View ArticlePubMedGoogle 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, et al: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415 (6871): 530-536. 10.1038/415530a.View ArticleGoogle Scholar
- van Galen JC, Hoefnagel JJ, Vermeer MH, Willemze R, Dijkman R, Tensen CP, de Boer WP, Meijer CJ, Oudejans JJ: Profiling of apoptosis genes identifies distinct types of primary cutaneous large B cell lymphoma. J Pathol. 2008, 215 (3): 340-346. 10.1002/path.2374.View ArticlePubMedGoogle Scholar
- Conley AP, Trent J, Zhang W: Recent progress in the genomics of soft tissue sarcomas. Curr Opin Oncol. 2008, 20 (4): 395-399. 10.1097/CCO.0b013e328302edc0.View ArticlePubMedGoogle Scholar
- Villanueva A, Toffanin S, Llovet JM: Linking molecular classification of hepatocellular carcinoma and personalized medicine: preliminary steps. Curr Opin Oncol. 2008, 20 (4): 444-453. 10.1097/CCO.0b013e328302c9e9.View ArticlePubMedPubMed CentralGoogle Scholar
- Straver ME, Glas AM, Hannemann J, Wesseling J, van de Vijver MJ, Rutgers EJ, Vrancken Peeters MJ, van Tinteren H, Van't Veer LJ, Rodenhuis S: The 70-gene signature as a response predictor for neoadjuvant chemotherapy in breast cancer. Breast Cancer Res Treat. 2010, 119 (3): 551-558. 10.1007/s10549-009-0333-1.View ArticlePubMedGoogle Scholar
- Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, Barrette TR, Anstet MJ, Kincead-Beal C, Kulkarni P, et al: Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia. 2007, 9 (2): 166-180. 10.1593/neo.07112.View ArticlePubMedPubMed CentralGoogle Scholar
- Segal E, Friedman N, Koller D, Regev A: A module map showing conditional activity of expression modules in cancer. Nat Genet. 2004, 36 (10): 1090-1098. 10.1038/ng1434.View ArticlePubMedGoogle Scholar
- Xu M, Kao MC, Nunez-Iglesias J, Nevins JR, West M, Zhou XJ: An integrative approach to characterize disease-specific pathways and their coordination: a case study in cancer. BMC Genomics. 2008, 9 (Suppl 1): S12-10.1186/1471-2164-9-S1-S12.View ArticlePubMedPubMed CentralGoogle Scholar
- Hubble J, Demeter J, Jin H, Mao M, Nitzberg M, Reddy TBK, Wymore F, Zachariah ZK, Sherlock G, Ball CA: Implementation of GenePattern within the Stanford Microarray Database. Nucleic Acids Res. 2009, 37 (Database): D898-901. 10.1093/nar/gkn786.View ArticlePubMedGoogle Scholar
- R Development Core Team: R: A Language and Environment for Statistical Computing. 2008, R Foundation for Statistical Computing, Vienna, Austria, [http://www.R-project.org]Google Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article 3Google Scholar
- Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, et al: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36 (Database): D480-484.View ArticlePubMedGoogle Scholar
- Perl: [http://www.perl.org/]
- Python: [http://www.python.org]
- Becker KG, Hosack DA, Dennis G, Lempicki RA, Bright TJ, Cheadle C, Engel J: PubMatrix: a tool for multiplex literature mining. BMC Bioinformatics. 2003, 4: 61-10.1186/1471-2105-4-61.View ArticlePubMedPubMed CentralGoogle Scholar
- Sayers EW, Barrett T, Benson DA, Bolton E, Bryant SH, Canese K, Chetvernin V, Church DM, Dicuccio M, Federhen S, et al: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2010, 38 (Database): D5-16. 10.1093/nar/gkp967.View ArticlePubMedGoogle Scholar
- Gwet KL: Handbook of Inter-Rater Reliability. 2010, ISBN 978-0970806222, 2Google Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.View ArticlePubMedPubMed CentralGoogle Scholar
- Romesburg HC: Cluster Analysis for Researchers. 2004, 340-ISBN 1-4116-0617-5Google Scholar
- Ishii T, Yasuchika K, Suemori H, Nakatsuji N, Ikai I, Uemoto S: Alpha-fetoprotein producing cells act as cancer progenitor cells in human cholangiocarcinoma. Cancer Lett. 2010, 294 (1): 25-34. 10.1016/j.canlet.2010.01.019.View ArticlePubMedGoogle Scholar
- Isonishi S, Ogura A, Kiyokawa T, Suzuki M, Kunito S, Hirama M, Tachibana T, Ochiai K, Tanaka T: Alpha-fetoprotein (AFP)-producing ovarian tumor in an elderly woman. Int J Clin Oncol. 2009, 14 (1): 70-73. 10.1007/s10147-008-0800-4.View ArticlePubMedGoogle Scholar
- Kesik V, Ozcan A, Sari E, Guven A, Kismet E, Koseoglu V: Alpha-fetoprotein producing tumor cells in children with Wilms' tumor. Fetal Pediatr Pathol. 2010, 29 (3): 127-132. 10.3109/15513811003777284.View ArticlePubMedGoogle Scholar
- Veeck J, Breuer E, Rose M, Chorovicer M, Naami A, Bektas N, Alkaya S, von Serenyi S, Horn F, Hartmann A, et al: [Novel prognostic marker in invasive breast cancer. ITIH5 expression is abrogated by aberrant promoter methylation]. Pathologe. 2008, 29 (Suppl 2): 338-346.View ArticlePubMedGoogle Scholar
- Veeck J, Chorovicer M, Naami A, Breuer E, Zafrakas M, Bektas N, Durst M, Kristiansen G, Wild PJ, Hartmann A, et al: The extracellular matrix protein ITIH5 is a novel prognostic marker in invasive node-negative breast cancer and its aberrant expression is caused by promoter hypermethylation. Oncogene. 2008, 27 (6): 865-876. 10.1038/sj.onc.1210669.View ArticlePubMedGoogle Scholar
- Himmelfarb M, Klopocki E, Grube S, Staub E, Klaman I, Hinzmann B, Kristiansen G, Rosenthal A, Durst M, Dahl E: ITIH5, a novel member of the inter-alpha-trypsin inhibitor heavy chain family is downregulated in breast cancer. Cancer Lett. 2004, 204 (1): 69-77. 10.1016/j.canlet.2003.09.011.View ArticlePubMedGoogle Scholar
- Mizutani K, Onda M, Asaka S, Akaishi J, Miyamoto S, Yoshida A, Nagahama M, Ito K, Emi M: Overexpressed in anaplastic thyroid carcinoma-1 (OEATC-1) as a novel gene responsible for anaplastic thyroid carcinoma. Cancer. 2005, 103 (9): 1785-1790. 10.1002/cncr.20988.View ArticlePubMedGoogle Scholar
- Yuan RH, Jeng YM, Pan HW, Hu FC, Lai PL, Lee PH, Hsu HC: Overexpression of KIAA0101 predicts high stage, early tumor recurrence, and poor prognosis of hepatocellular carcinoma. Clin Cancer Res. 2007, 13 (18 Pt 1): 5368-5376.View ArticlePubMedGoogle Scholar
- Montaner D, Minguez P, Al-Shahrour F, Dopazo J: Gene set internal coherence in the context of functional profiling. BMC Genomics. 2009, 10: 197-10.1186/1471-2164-10-197.View ArticlePubMedPubMed CentralGoogle Scholar
- Kanzler S, Galle PR: Apoptosis and the liver. Semin Cancer Biol. 2000, 10 (3): 173-184. 10.1006/scbi.2000.0318.View ArticlePubMedGoogle Scholar
- Banchereau J, Palucka AK: Dendritic cells as therapeutic vaccines against cancer. Nat Rev Immunol. 2005, 5 (4): 296-306. 10.1038/nri1592.View ArticlePubMedGoogle Scholar
- Crowe NY, Coquet JM, Berzins SP, Kyparissoudis K, Keating R, Pellicci DG, Hayakawa Y, Godfrey DI, Smyth MJ: Differential antitumor immunity mediated by NKT cell subsets in vivo. J Exp Med. 2005, 202 (9): 1279-1288. 10.1084/jem.20050953.View ArticlePubMedPubMed CentralGoogle Scholar
- Dunn GP, Old LJ, Schreiber RD: The three Es of cancer immunoediting. Annu Rev Immunol. 2004, 22: 329-360. 10.1146/annurev.immunol.22.012703.104803.View ArticlePubMedGoogle Scholar
- Fernandez NC, Levraud JP, Haddada H, Perricaudet M, Kourilsky P: High frequency of specific CD8+ T cells in the tumor and blood is associated with efficient local IL-12 gene therapy of cancer. J Immunol. 1999, 162 (1): 609-617.PubMedGoogle Scholar
- Fernandez NC, Lozier A, Flament C, Ricciardi-Castagnoli P, Bellet D, Suter M, Perricaudet M, Tursz T, Maraskovsky E, Zitvogel L: Dendritic cells directly trigger NK cell functions: cross-talk relevant in innate anti-tumor immune responses in vivo. Nat Med. 1999, 5 (4): 405-411. 10.1038/7403.View ArticlePubMedGoogle Scholar
- Kuphal S, Bauer R, Bosserhoff AK: Integrin signaling in malignant melanoma. Cancer Metastasis Rev. 2005, 24 (2): 195-222. 10.1007/s10555-005-1572-1.View ArticlePubMedGoogle Scholar
- Frisch SM, Ruoslahti E: Integrins and anoikis. Curr Opin Cell Biol. 1997, 9 (5): 701-706. 10.1016/S0955-0674(97)80124-X.View ArticlePubMedGoogle Scholar
- Meredith JE, Fazeli B, Schwartz MA: The extracellular matrix as a cell survival factor. Mol Biol Cell. 1993, 4 (9): 953-961.View ArticlePubMedPubMed CentralGoogle Scholar
- Bissell MJ, Radisky DC, Rizki A, Weaver VM, Petersen OW: The organizing principle: microenvironmental influences in the normal and malignant breast. Differentiation. 2002, 70 (9-10): 537-546. 10.1046/j.1432-0436.2002.700907.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Schwock J, Dhani N, Hedley DW: Targeting focal adhesion kinase signaling in tumor growth and metastasis. Expert Opin Ther Targets. 2010, 14 (1): 77-94. 10.1517/14728220903460340.View ArticlePubMedGoogle Scholar
- Benjamin LC, Allan JS, Madsen JC: Cytokines in immunity and allograft rejection. Crit Rev Immunol. 2002, 22 (4): 269-279.PubMedGoogle Scholar
- An Y, Bekesova S, Edwards N, Goldman R: Peptides in low molecular weight fraction of serum associated with hepatocellular carcinoma. Dis Markers. 2010, 29 (1): 11-20.View ArticlePubMedPubMed CentralGoogle Scholar
- Beppu T, Gil-Bernabe P, Boveda-Ruiz D, D'Alessandro-Gabazza C, Matsuda Y, Toda M, Miyake Y, Shiraki K, Murata M, Murata T, et al: High incidence of tumors in diabetic thrombin activatable fibrinolysis inhibitor and apolipoprotein E double deficient mice. J Thromb Haemost. 2010Google Scholar
- Gerotziafas GT, Papageorgiou C, Hatmi M, Samama MM, Elalamy I: Clinical studies with anticoagulants to improve survival in cancer patients. Pathophysiol Haemost Thromb. 2008, 36 (3-4): 204-211.View ArticlePubMedGoogle Scholar
- Pengo V, Denas G, Jose SP, Pengo MF: Cancer prevention and vitamin K antagonists: an overview. Thromb Res. 2010, 125 (Suppl 2): S103-105.View ArticlePubMedGoogle Scholar
- Ondrey F: Peroxisome proliferator-activated receptor gamma pathway targeting in carcinogenesis: implications for chemoprevention. Clin Cancer Res. 2009, 15 (1): 2-8. 10.1158/1078-0432.CCR-08-0326.View ArticlePubMedGoogle Scholar
- Sertznig P, Seifert M, Tilgen W, Reichrath J: Peroxisome proliferator-activated receptor (PPAR) and vitamin D receptor (VDR) signaling pathways in melanoma cells: promising new therapeutic targets?. J Steroid Biochem Mol Biol. 2010, 121 (1-2): 383-386. 10.1016/j.jsbmb.2010.03.003.View ArticlePubMedGoogle Scholar
- Konishi I, Hiasa Y, Shigematsu S, Hirooka M, Furukawa S, Abe M, Matsuura B, Michitaka K, Horiike N, Onji M: Diabetes pattern on the 75 g oral glucose tolerance test is a risk factor for hepatocellular carcinoma in patients with hepatitis C virus. Liver Int. 2009, 29 (8): 1194-1201. 10.1111/j.1478-3231.2009.02043.x.View ArticlePubMedGoogle Scholar
- Huo TI, Hsu CY, Huang YH, Hsia CY, Lin HC, Lee PC, Loong CC, Chiang JH, Chiou YY, Lee SD: Diabetes mellitus as an independent prognostic predictor and its association with renal dysfunction in patients with hepatocellular carcinoma. Liver Int. 2010, 30 (2): 198-207. 10.1111/j.1478-3231.2009.02143.x.View ArticlePubMedGoogle Scholar
- El-Serag HB, Hampel H, Javadi F: The association between diabetes and hepatocellular carcinoma: a systematic review of epidemiologic evidence. Clin Gastroenterol Hepatol. 2006, 4 (3): 369-380. 10.1016/j.cgh.2005.12.007.View ArticlePubMedGoogle Scholar
- Reya T, Morrison SJ, Clarke MF, Weissman IL: Stem cells, cancer, and cancer stem cells. Nature. 2001, 414 (6859): 105-111. 10.1038/35102167.View ArticlePubMedGoogle Scholar
- Eaves CJ: Cancer stem cells: Here, there, everywhere?. Nature. 2008, 456 (7222): 581-582. 10.1038/456581a.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/4/53/prepub