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.
Methods
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.
Results
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'.
Conclusion
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.
Cancer is one of the leading causes of death [1]. 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 [2]. 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 [2]. Nevertheless, since the year 1990, cancer death rate continuously declined despite increasing numbers of new cancer cases [1], 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 [11], 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 [10]. 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 [13]. 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.
Methods
Gene expression data
Gene expression data were obtained from microarray experiments, downloaded from the Stanford Microarray Database [13]. 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 [13]. 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 [14] 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 [15]. 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.
A gene was considered upregulated (respectively downregulated) if it had a M-value > 1.0 (resp. < -1.0) and a p-value smaller than 0.01 in a one sided one-sample t-test. This procedure defined between 391 and 1433 up- and between 424 and 2150 downregulated tumor associated genes (figure 1a). Note: Although the p-values produced this way are merely used as a score and are not interpreted as type-1 error rates, it might be interesting to assess the significance of the individual tests. We therefore performed an FDR correction using "p.adjust" and provided the list of multiple testing corrected p-values in the additional files (Additional file 1).
Figure 1
Schematic drawing of the workflow - a) gene expression data, b) functional data and c) the functional cancer map.
Functional data
Metabolic reaction and signaling event information were acquired from the PATHWAY database of the Kyoto Encyclopedia of Genes and Genomes (KEGG) [16]. 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 [17] 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.
To get a better overview, we implemented a Python [18] script, using the PIL module, to visualize significant enriched KEGG pathway maps (p < 0.05) within the tumor entities into a clear structured functional cancer map (figure 1c) (figure 2). The same procedure as described here for the set of all differentially expressed genes can be applied to the set of upregulated resp. downregulated genes separately. The corresponding functional cancer maps are given in the additional files (Additional file 2, Additional file 3). Note: A separately generated functional cancer map without prior filtering for outlier genes (pooled variance and tissue specific variance) was added to the additional files (Additional file 4).
Figure 2
Functional cancer map: Functional expression profile of significant enriched KEGG pathway maps across 28 tumor entities assigned to 16 tumor classes. The significance of a tumor entity to a KEGG pathway map is symbolized by the diameter of the spot, the greater the diameter the more significant enriched is the corresponding tumor entity - KEGG pathway map relation. The color of the spot gives information on the tumor entity and its belonging tumor class. Information on the KEGG classification schema is supported by the color gradient on the left side.
Literature Mining
Using the web-based tool PubMatrix [19] for literature mining of PubMed [20], 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 [21] of corresponding columns in the binarized functional cancer map. This coefficient is a statistical measure of agreement for qualitative items.
We calculated the quadratic agreement matrix among the tumor entities. The result was visualized with the R function heatmap [22]. Euclidean distance average linkage hierarchical clustering [23] was chosen for the dendrogram included in the plot (figure 3). Visualizations of the tumor similarity matrices with respect to upregulated and downregulated genes are given in the additional files (Additional file 6, Additional file 7).
Figure 3
Tumor phylogeny: Calculation of the tumor phylogeny was done by applying the Cohen's Kappa Coefficient to the binarized functional cancer map, tumor entities - KEGG pathway map relations respectfully.
Results
Gene expression data
In order to investigate individual and conserved signaling and metabolic events in tumor progression, we analyzed a dataset of 649 human two-color cDNA microarray tumor samples obtained from the Stanford Microarray Database [13]. These samples covered 28 different tumor entities which were further categorized into 16 tumor classes (table1). Analyzing these gene expression profiles revealed a total of 10385 genes differentially expressed in at least one of the tumor entities. Among highly conserved differentially expressed genes were apolipoprotein (APOH), previously not mentioned in context with cancer, T cell receptor beta constant 1 (TRBC1) a widely undescribed gene and also not associated in context with cancer and alpha-fetoprotein (AFP), a major plasma protein reported to be associated with several carcinomas [24–26]. Furthermore Orosomucoid 2 (ORM2) another plasma protein whose specific function has not yet been determined and not linked to carcinogenesis was found to be regulated in 25 tumor tissues. Furthermore, two genes were identified to be regulated in 24 tumor tissues, inter-alpha (globulin) inhibitor H2 (ITIH2), a gene involved in extracellular matrix stabilization and related to prevention of tumor metastasis [27–29] and the functionally undefined gene KIAA0101 which is associated with hepatocellular and thyroid carcinoma [30, 31].
Table 1
The tumor sample distribution of the 649 microarrays with corresponding controls among the tumor classes and tumor entities used in our analysis
Acute Lymphoblastic Leukemia Bone Marrow, Grade Common
Stratagene
35
Bone marrow
Acute Lymphoblastic Leukemia Bone Marrow, Grade Pre-B
Stratagene
33
Bone
Extraskeletal Myxoid Chondrosarcoma
CRH
10
Breast
Breast Carcinoma
Stratagene
94
Breast
Breast Carcinoma Invasive Ductal Carcinoma, Grade 2
Stratagene
31
Breast
Breast Carcinoma Invasive Ductal Carcinoma, Grade 3
Stratagene
19
Breast
Breast Carcinoma Invasive Lobular Carcinoma, Grade 2
Stratagene
19
Liver
Human Hepatocellular Carcinoma
Stratagene
19
Lymhpocytes
Follicular Lymphoma
CRG
24
Mixture
Solitary Fibrous Tumor
CRG
7
Muscle
Leiomyosarcoma
CRG
7
Ovarian
Ovarian Carcinoma
CRD
39
Pancreas
Adenocarcinoma, Pancreas
CRH
10
Pancreas
Adenocarcinoma, Pancreas, Grade II
CRH
5
Pancreas
Adenocarcinoma, Pancreas, Grade III
CRH
6
Prostate
Prostate Tumor, Grade T2b
Stratagene
23
Prostate
Prostate Tumor, Grade T3a
Stratagene
19
Prostate
Prostate Tumor, Grade T3b
Stratagene
9
Renal
Conventional Renal Cell Carcinoma, Grade 1
Stratagene
9
Renal
Conventional Renal Cell Carcinoma, Grade 2
Stratagene
34
Renal
Conventional Renal Cell Carcinoma, Grade 3
Stratagene
94
Renal
Conventional Renal Cell Carcinoma, Grade 4
Stratagene
40
Skin
Dermatofibrosarcoma Protuberans
CRG
6
Soft Tissue
Synovial Sarcoma
CRG
5
Stomach
Stomach Cancer Before Treatment
CRF
14
Stomach
Stomach Cancer, Placebo After 1 Year
CRF
13
Testis
Testis Seminoma
CRG
16
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) [16]. 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).
Literature Mining
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 [19]. 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 [21] 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).
Discussion
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 [32]. 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 [33]. 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 [39]. 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 [42] 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 [43]. 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 [43]. 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 [44]. 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 [51]. Also Huo et al. found that Diabetes mellitus independently predicts decreased survival in HCC patients [52]. Finally, El-Serag et al. reported that Diabetes mellitus is associated with an increased risk for HCC [53]. Interestingly the PPARgamma activators are utilized for the treatment of diabetes mellitus pointing towards a congruent signaling in of PPAR pathway and diabetes signatures [50].
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 [21] 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].
Conclusion
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.
Declarations
Authors’ Affiliations
(1)
Department of Medicine I, Johannes Gutenberg University
(2)
and Molecular Biotechnology, Bioquant, Institut of Pharmacy, University of Heidelberg, INF 267
(3)
and Department of Cytopathology, Institute for Pathology, Johannes Gutenberg University, Heinrich Heine University
(4)
Gene Center Munich, Department of Chemistry and Biochemistry, Ludwig-Maximilians-University
References
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
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
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
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
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.
