Identification of Transcriptional Signatures of Colon Tumor Stroma by a Meta-Analysis

Uddin, Md. Nazim; Li, Mengyuan; Wang, Xiaosheng

doi:https://doi.org/10.1155/2019/8752862

Journal of Oncology

On this page

Abstract Background Methods Results Discussion Conclusions Abbreviations Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

Role of the Microenvironment in Gastrointestinal Tumors

View this Special Issue

Research Article | Open Access

Volume 2019 | Article ID 8752862 | https://doi.org/10.1155/2019/8752862

Identification of Transcriptional Signatures of Colon Tumor Stroma by a Meta-Analysis

Md. Nazim Uddin,^1,2,3Mengyuan Li,^1,2,3and Xiaosheng Wang^1,2,3

Guest Editor: Nathaniel Weygant

Received28 Feb 2019

Accepted31 Mar 2019

Published02 May 2019

Abstract

Background. The tumor stroma plays pivotal roles in influencing tumor growth, invasion, and metastasis. Transcriptional signatures of colon tumor stroma (CTS) are significantly associated with prognosis of colon cancer. Thus, identification of the CTS transcriptional features could be useful for colon cancer diagnosis and therapy. Methods. By a meta-analysis of three CTS gene expression profiles datasets, we identified differentially expressed genes (DEGs) between CTS and colon normal stroma. Furthermore, we identified the pathways, upstream regulators, and protein-protein interaction (PPI) network that were significantly associated with the DEGs. Moreover, we analyzed the enrichment levels of immune signatures in CTS. Finally, we identified CTS-associated gene signatures whose expression was significantly associated with prognosis in colon cancer. Results. We identified numerous significantly upregulated genes (such as CTHRC1, NFE2L3, SULF1, SOX9, ENC1, and CCND1) and significantly downregulated genes (such as MYOT, ASPA, KIAA2022, ARHGEF37, BCL-2, and PPARGC1A) in CTS versus colon normal stroma. Furthermore, we identified significantly upregulated pathways in CTS that were mainly involved in cellular development, immune regulation, and metabolism, as well as significantly downregulated pathways in CTS that were mostly metabolism-related. Moreover, we identified upstream TFs (such as SUZ12, NFE2L2, RUNX1, STAT3, and SOX2), kinases (such as MAPK14, CSNK2A1, CDK1, CDK2, and CDK4), and master metabolic transcriptional regulators (MMTRs) (such as HNF1A, NFKB1, ZBTB7A, GATA2, and GATA5) regulating the DEGs. We found that CD8+ T cells were more enriched in CTS than in colon normal stroma. Interestingly, we found that many of the DEGs and their regulators were prognostic markers for colon cancer, including CEBPB, PPARGC1, STAT3, MTOR, BCL2, JAK2, and CDK1. Conclusions. The identification of CTS-specific transcriptional signatures may provide insights into the tumor microenvironment that mediates the development of colon cancer and has potential clinical implications for colon cancer diagnosis and treatment.

1. Background

The tumor stroma is an important component of the tumor microenvironment (TME) and plays key roles in the tumor development [1]. Stromal cells are composed of many different types of cells, including vascular endothelial cells, pericytes, adipocytes, fibroblasts, osteoblasts, chondrocytes, extracellular matrix (ECM), and bone-marrow mesenchymal stromal cells [2]. The tumor stroma can promote ECM remodeling, cellular migration, neoangiogenesis, invasion, immunosurveillance evasion, and drug resistance of tumors [3]. Colorectal cancer (CRC) is the fourth most common cancer and a leading cause of cancer mortality worldwide [4]. Transcriptional signatures of CRC stromal cells have been associated with poor prognosis in CRC [5]. Isella et al. demonstrated that the gene signatures of CRC stromal cells (cancer-associated fibroblasts, leukocytes, and endothelial cells) were significantly upregulated in the stem/serrated/mesenchymal transcriptional subtype of CRC which had a poor prognosis [6]. Calon et al. showed that the CRC stromal transcriptional signatures correlated with disease relapse [5]. These prior studies exhibited the significant roles of tumor stroma in CRC growth, invasion, and metastasis.

In this study, we performed a meta-analysis of three colon tumor stromal transcriptome datasets using the bioinformatics approach. We identified differentially expressed genes (DEGs) between colon tumor stroma (CTS) and normal stroma. On the basis of these DEGs, we identified their associated pathways, upstream regulators, and protein-protein interaction (PPI) network and certain prognostic markers that were associated with survival of colon cancer patients. We also analyzed the enrichment levels of immune signatures in CTS. This study provides insights into CTS molecular features that could have clinical implications for colon cancer diagnosis and treatment.

2. Methods

2.1. Datasets

We searched the NCBI Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/) using the keywords “colon cancer,” “stroma,” and “tumor stroma” and identified three CTS gene expression profiles datasets (GSE31279, GSE35602, and GSE46824) [7–9]. In survival analyses, we used the TCGA colon cancer dataset (https://portal.gdc.cancer.gov/) and a SurvExpress (http://bioinformatica.mty.itesm.mx/SurvExpress) built-in dataset (colon metabase) [10]. A summary of these datasets is shown in Supplementary Table S1.

2.2. Identification of DEGs between CTS and Normal Stroma

We used the web tool Network Analyst [11] to identify the DEGs between CTS and normal stroma. The ComBat method [12] in the tool was utilized to remove batch effects from the three CTS datasets (Supplementary Figure S1). Each individual dataset was normalized by base-2 log transformation and quantile normalization, and the R package “limma” was utilized to identify the DEGs between CTS and normal stroma. A meta-analysis of the three datasets was performed using Cochran's combination test [13]. The false discovery rate (FDR), calculated by the Benjamini–Hochberg method [14], was used to adjust for multiple tests. We determined the DEGs with a threshold of absolute combined effect size (ES) >0.82 and FDR<0.05.

2.3. Gene-Set Enrichment Analysis

We performed gene-set enrichment analysis of the DEGs by GSEA [15]. The KEGG pathways significantly associated with the upregulated and the downregulated DEGs were identified (FDR < 0.05), respectively.

2.4. Identification of Transcription Factors (TFs), Kinases, and Master Metabolic Transcriptional Regulators (MMTRs) That Are Significantly Associated with the DEGs

To link gene expression signatures to upstream cell signaling networks, we used eXpression2Kinases [16] to identify the upstream TFs and kinases that regulate the DEGs and utilized iRegulon [17] to identify the MMTRs of the DEGs.

2.5. Identification of PPI Network of the DEGs

We employed Network Analyst [11] to construct a PPI network of the DEGs [11]. Two types of modules (function-first modules and connection-first modules) of the PPI network were extracted. The function-first modules (FFMs) were constructed by pathway enrichment analysis and the connection-first modules (CFMs) were identified by the random walk-based algorithm [18].

2.6. Comparison of the Enrichment Levels of CD8+ T Cells between Two Classes of Samples

The enrichment level of CD8+ T cells in a sample was evaluated by the expression level of CD8A. We compared the enrichment levels of CD8+ T cells between two groups of samples using Student's t-test.

2.7. Identification of DEGs between High-Stroma-Content and Low-Stroma-Content TCGA Colon Cancer Samples

We used ESTIMATE [19] to quantify the intratumoral stromal content (stroma score) of TCGA colon cancer samples. We identified the DEGs between high-stroma-content (stroma score > median) and low-stroma-content (stroma score < median) tumors using Student's t-test.

2.8. Survival Analyses

We compared the overall survival (OS) and the disease-free survival (DFS) of colon cancer patients classified based on gene expression levels (expression levels > median versus expression levels < median). Kaplan-Meier survival curves were used to show the survival differences, and the log-rank test was utilized to evaluate the significance of survival differences. The individual prognostic genes were identified and were fitted in a multivariate Cox regression model. SurvExpress [10] was used for the multivariate survival analysis.

3. Results

3.1. Identification of DEGs between CTS and Normal Stroma

We identified 694 DEGs between CTS and normal stroma by the meta-analysis. These DEGs included 295 downregulated and 399 upregulated genes in CTS (Supplementary Tables S2 and S3). Figure 1 shows the top 25 upregulated and top 25 downregulated genes in CTS ranked on the basis of the combined ES (the detailed results of statistical analysis for the top 10 upregulated and top 10 downregulated genes in CTS are shown in Supplementary Tables S4). CTHRC1, a gene involved in vascular remodeling, bone formation, and developmental morphogenesis, was upregulated in CTS with the highest ES. It has been shown that CTHRC1 could promote human CRC cell proliferation and invasion by activating Wnt/PCP signaling [20]. This gene also plays an important role in promoting ovarian cancer cell adhesion, migration, and metastasis through the activation of integrin β3/FAK signaling [21]. NFE2L3, a gene regulating the cell cycle progression in colon cancer [22], was upregulated in CTS with the second highest ES. Interestingly, both CTHRC1 and NFE2L3 have been indicated as useful biomarker candidates for CRC diagnosis because of their overexpression in adenomas and CRC relative to normal tissue [23]. SULF1, whose expression in tumor stroma is a prognostic marker in advanced pancreatic cancer [24], was upregulated in CTS with the third highest ES. The overexpression of this gene has been associated with a poor prognosis in urothelial carcinoma [25]. SOX9, the gene upregulated in CTS with the fourth highest ES, has been shown to be overexpressed in CRC and its overexpression was an independent adverse prognosticator in CRC [26]. Some other genes upregulated in CTS have been demonstrated to be overexpressed in CRC and their expression was negatively associated with CRC prognosis, such as ENC1, CCND1, VCAN, SEMA5A, and NOS3 [27–31]. Interestingly, both PCDH17 and BCL6B were upregulated in CTS, while they had reduced expression in CRC [32, 33]. It indicates that PCDH17 and BCL6B could be specifically expressed in CTS cells but not in colon cancer cells.

Many of the significantly downregulated genes in CTS have been associated with CRC [34–37]. For example, MYOT, ASPA, and KIAA2022 were downregulated in CRC [34], the downregulation of ARHGEF37 was associated with a poor prognosis in CRC [35], higher expression levels of BCL-2 were correlated with a better survival prognosis in CRC [36], and PPARGC1A was a negative predictor for CRC prognosis [37].

Altogether, a number of the abnormally expressed genes in CTS compared to colon normal stroma identified by the meta-analysis have been associated with CRC pathology and prognosis.

3.2. Identification of Pathways Significantly Associated with the DEGs

GSEA [15] identified 44 KEGG pathways that were significantly associated with the upregulated genes in CTS. These pathways were mainly involved in cellular development (p53 signaling, Wnt signaling, apoptosis, Notch signaling, focal adhesion, endocytosis, ECM-receptor interaction, cell adhesion molecules, adherens junction, tight junction, gap junction, and regulation of actin cytoskeleton), immune regulation (leukocyte transendothelial migration, complement and coagulation cascades, natural killer cell mediated cytotoxicity, Toll-like receptor, chemokine signaling, and cytokine-cytokine receptor interaction), and metabolism (purine metabolism and pyrimidine metabolism) (Figure 2, Supplementary Table S5). Previous studies have shown that some of these pathways were significantly associated with colon cancer [38–41]. For example, the Wnt and Notch pathways were associated with colon cancer development [38, 39]. The cytokine-cytokine receptor interaction pathway was significantly enriched in CRC [34]. The ECM and ECM-associated proteins [39], the glycosaminoglycan metabolism, and chondroitin sulfate/dermatan sulfate metabolism pathways played key roles in mediating tumor microenvironment [40, 41].

In addition, GSEA identified six KEGG pathways that were significantly associated with the downregulated genes in CTS (Supplementary Figure S2). Most of these pathways were metabolism-related, including purine metabolism, histidine metabolism, glycine, serine, and threonine metabolism, and drug metabolism-cytochrome p450. These pathways have been associated with colon and other cancers [42–44]. For example, impaired purine metabolism was associated with the progression of cancer [42]. Histidine metabolism could boost cancer therapy [43]. Cytochrome P450 enzymes were associated with the metabolism of anticancer drugs and their expression was associated with a poor prognosis in CRC patents [44].

3.3. Identification of Upstream TFs, Kinases, and MMTRs Significantly Associated with the DEGs

We identified 11 significant upstream TFs regulating the DEGs, including SUZ12, NFE2L2, RUNX1, ESR1, STAT3, TCF3, FOSL2, SALL4, AR, SMC3, and SOX2, of which the genes encoding RUNX1 and SALL4 were upregulated in CTS (Figure 3(a)). Most of these TFs have been associated with colon cancer [45–49]. For example, SUZ12 was the most significant upstream TF which could contribute to the CRC development [45]. RUNX1 mutations were associated with the CRC risk [46]. TCF3 and FOSL2 were associated with the tumorigenesis of CRC [47, 48]. The overexpression of SOX2 was associated with the progression and a poor prognosis in colon cancer [49].

(a)

(b)

(c)

Moreover, we identified 124 significant protein kinases that regulate the DEGs (Figure 3(b), Supplementary Table S6). These kinases mainly included cell cycle regulation kinases (CDKs), signaling MAP kinases (MAPKs, MAP2Ks, and MAP3Ks), and ribosomal kinases (RPS6KA1, RPS6KA3, and RPS6KA5). MAPK14 was the most significant upstream kinase negatively regulating the formation of colitis-associated colon tumors [50]. Furthermore, we constructed a TF-kinase interaction network of these TFs and kinases (Figure 3(c)). In the network, the most connected TFs included SUZ12, NFE2L2, RUNX1, STAT3, FOSL2, AR, SMC3, ESR1, and TCF3, and the most connected kinases included MAPK14, CDK1, CSNK2A1, CDK2, MAPK3, HIPK2, ERK1, and CDK4. It indicates that the cell cycle regulation may play a pivotal role in CTS.

MMTRs are interesting biomarkers and targets for metabolism-targeted cancer therapy [51]. We identified 9 (HNF1A, NFKB1, ZBTB7A, ATF6, TEAD4, TFAP2B, JAZF1, FNTB, and EP300) and 12 (PKNOX2, GATA2, MAPK10, TEAD1, TOX, MEF2A, GATA5, ELK1, MAZ, NHLH1, ATF1, and RAD21) MMTRs for the upregulated and the downregulated genes in CTS, respectively (Supplementary Table S7), and built the regulatory networks associated with these MMTRs (Figure 4). In the networks, ATF6 (activating transcription factor 6), a TF regulating unfolded protein response during endoplasmic reticulum (ER) stress, targeted 163 upregulated genes, and PKNOX2 (PBX/knotted 1 homeobox 2), which plays key roles in regulating cell proliferation, differentiation, and death, targeted 131 downregulated genes. Interestingly, two members of the GATA family of TFs (GATA2 and GATA5) were the MMTRs that regulated the downregulated genes in CTS (Figure 4(b)).

(a)

(b)

Altogether, the identification of upstream TFs, kinases, and MMTRs significantly associated with the DEGs may provide insights into the TME that mediates the development of colon cancer.

3.4. CD8+ T Cells Are More Enriched in CTS than in Normal Stroma

We compared the enrichment levels of CD8+ T cells between CTS and normal stroma and found that CD8+ T cells showed significantly higher enrichment levels in CTS than in normal stroma (Student's t-test, p=0.016) (Figure 5). This suggests an antitumor immune response activity in the TME of colon cancer.

3.5. Identification of Prognostic Factors in Colon Cancer Based on the DEGs and Their Upstream Regulators

We investigated the association between the transcriptional signatures of CTS and survival prognosis (overall survival (OS) and disease-free survival (DFS)) in the TCGA colon cancer dataset. The transcriptional signatures included the top 10 upregulated and top 10 downregulated genes in CTS on the basis of ES, 45 hub genes (≥3 degrees) from the zero-order PPI network of the DEGs (Supplementary Table S8), and the genes encoding 11 TFs, 124 kinases, and 21 MMTRs regulating the DEGs. We found that the expression of many of these transcriptional signatures was significantly associated with the survival of colon cancer patients. For example, the expression of CEBPB, a gene significantly upregulated in CTS and a hub node in the PPI network, had a significant negative correlation with OS in colon cancer (Figure 6(a)). The negative correlation between CEBPB expression and survival has also been demonstrated in other cancer types, such as high-grade serous ovarian cancer [52]. PPARGC1 was significantly downregulated in CTS and was a hub node in the PPI network, while its expression had a significant positive correlation with OS in colon cancer (Figure 6(a)). PPARGC1A was indicated as a tumor suppressor in colon cancer [53] and ovarian cancer [54], as well as a negative prognostic biomarker for CRC [37]. Our data indicate that the deregulation of these genes in CTS is prognostic for colon cancer patients.

(a)

(b)

(c)

Figure 6

The CTS gene signatures whose expression is associated with prognosis in colon cancer. (a) Kaplan-Meier survival curves show the gene signatures whose expression is significantly associated with overall survival (OS) in colon cancer in the TCGA colon cancer dataset (log-rank test, p<0.05). (b) Multivariate Cox regression analysis shows that the OS-associated CTS gene signatures are prognostic for OS in colon cancer in a SurvExpress built-in dataset (colon metabase) [10]. (c) Multivariate Cox regression analysis shows that the DFS-associated CTS gene signatures are prognostic for DFS in colon cancer in a SurvExpress built-in dataset (colon metabase) [10]. DFS: disease-free survival.

Among the upstream regulators (TFs, kinases, and MMTRs) of the DEGs, the expression of STAT3, RPS6KA5, IKBKE, ERBB2, MTOR, and NFKB1 had a positive correlation with OS in colon cancer, while the expression of CDK1, CDK5, and BRD2 had a negative correlation with OS in colon cancer (Figure 6(a)). The deregulation of these genes has been associated with tumor progression in a wide variety of cancer types [55–60].

In addition, we identified 18 transcriptional signatures of CTS whose expression was significantly associated with DFS in colon cancer individually (Supplementary Figure S3). These genes included CEBPB, BCL2, PAN2, NOS3, FTL, ARHGEF37, SMC3, EP300, JAK2, RPS6KA3, RPS6KA1, PRKACA, HIPK1, HIPK2, MAPK8, GSK3A, CLK2, and CDK3. It indicates that these CTS transcriptional signatures could be biomarkers for colon cancer relapse.

Furthermore, we used the multivariate analysis to validate the association between the prognostic CTS transcriptional signatures and survival using the colon metabase data [10]. For OS analysis, a total of 482 patients were split into two groups: high-risk group (N=241) versus low-risk group (N=241) based on the prognostic index (Supplementary Figure S4A). As expected, the high-risk group had worse OS than the low-risk group (Figure 6(b)). Similarly, for DFS analysis, we divided patients into the high-risk group (N=272) and the low-risk group (N=273) based on the prognostic index (Supplementary Figure S4B) and found that the high-risk group had worse DFS compared to the low-risk group (Figure 6(c)). These results proved the prognostic value of these CTS transcriptional signatures in colon cancer.

4. Discussion

The tumor stroma constitutes an important component of the TME that mediates tumor growth, immune evasion, and metastasis [1]. Thus, it is important to identify molecular features in the tumor stroma. To this end, we performed a meta-analysis of three CTS transcriptome datasets for identifying CTS-associated transcriptional signatures. We identified a number of upregulated and downregulated genes in CTS compared to colon normal stroma. Furthermore, we identified upregulated and downregulated pathways significantly associated with these deregulated genes in CTS. The upregulated pathways were mainly involved in cellular development, immune regulation, and metabolism, and the downregulated pathways were mostly metabolism-related. These results revealed the abnormal alterations of cellular development, immune regulation, and metabolism pathways in CTS. We found that CD8+ T cells were more enriched in CTS than in colon normal stroma, suggesting an immune infiltration microenvironment in CTS. Furthermore, we identified numerous CTS transcriptional signatures whose expression was significantly associated with prognosis in colon cancer, such as CEBPB, PPARGC1, STAT3, MTOR, BCL2, JAK2, and CDK1. These transcriptional signatures are mainly involved in immune regulation (CEBPB, STAT3, and JAK2), metabolism (PPARGC1 and MTOR), cell cycle (CDK1), and apoptosis (BCL2), suggesting that the deregulation of these pathways in CTS may contribute to the altered prognosis in colon cancer.

To verify the association of the identified transcriptional signatures with CTS, we analyzed the TCGA colon cancer dataset. We divided these cancers into high-stroma-content and low-stroma-content groups on the basis of their intratumoral stromal content evaluated by ESTIMATE [19] and found that 153 upregulated genes in CTS had significantly higher expression levels in the high-stroma-content group than in the low-stroma-content group. These genes included 18 hub genes in the PPI network of DEGs and 6 TFs, 40 kinases, and 12 MMTRs encoding genes that regulated the DEGs (Supplementary Figure S5, Table S9). We also found 27 downregulated genes in CTS which had significantly lower expression levels in the high-stroma-content group, including 14 hub genes, and genes encoding 2 TFs, 18 kinases, and 3 MMTRs (Supplementary Figure S5, Table S9). Interestingly, most of the downregulated hub genes in CTS were also downregulated in the high-stroma-content colon cancers (Supplementary Figure S5). These results indicate that many transcriptional signatures of CTS identified by the meta-analysis are tumor stroma-specific. In addition, we found that CD8+ T cells had significantly higher enrichment levels in CTS versus colon normal stroma (Student's t-test, p=0.016), as well as in the high-stroma-content colon cancers versus the low-stroma-content colon cancers (Student's t-test, p=3.310⁻⁸). It indicates that CD8+ T cells tend to have elevated infiltration in the TME of colon cancer. Interestingly, we found that the higher enrichment levels of CD8+ T cells were associated with better DFS in the low-stroma-content colon cancers, but not in the high-stroma-content colon cancers (Figure 7). It suggests that the immune cells exert an antitumor effect only when they have infiltrated into tumor cells and that the immune cells in the tumor stroma may not have such a direct antitumor effect.

This study has identified a number of CTS-associated transcriptional signatures that could be biomarkers for colon cancer diagnosis and prognosis and may provide therapeutic targets for colon cancer. However, to translate these findings into clinical application, further experimental and clinical validation would be necessary.

5. Conclusions

The identification of CTS-specific transcriptional features may provide insights into the mechanism that mediates the development of colon cancer and thus has potential clinical implications for colon cancer diagnosis and treatment.

Abbreviations

CTS:	Colon tumor stroma
DEGs:	Differentially expressed genes
PPI:	Protein-protein interaction
ES:	Effect size
TME:	Tumor microenvironment
ECM:	Extracellular matrix
TFs:	Transcription factors
MMTRs:	Master metabolic transcription factors
OS:	Overall survival
DFS:	Disease-free survival
CRC:	Colorectal cancer
GEO:	Gene Expression Omnibus
TCGA:	The Cancer Genome Atlas
FDR:	False discovery rate
GSEA:	Gene-set enrichment analysis
FFM:	Function-first module
CFM:	Connection-first module.

Data Availability

The datasets (GSE31279, GSE35602, and GSE46824) were downloaded from the NCBI GEO database (https://www.ncbi.nlm.nih.gov/geo/), the TCGA colon cancer dataset was downloaded from the website https://portal.gdc.cancer.gov/, and the colon metabase dataset was from SurvExpress (http://bioinformatica.mty.itesm.mx/SurvExpress).

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Md. Nazim Uddin performed data analyses and helped write and prepare the manuscript. Mengyuan Li performed data analyses and helped prepare the manuscript. Xiaosheng Wang conceived the research, designed analysis strategies, and wrote the manuscript. All the authors read and approved the final manuscript.

Acknowledgments

This work was supported by China Pharmaceutical University (grant numbers 3150120001 and 2632018YX01 to Xiaosheng Wang).

Supplementary Materials

Supplementary Figure S1: illustration of PCA and density plots as validation tools for batch effect removal. Plot of principal components: (A) before batch effect removal and (B) after batch effect removal. Plot of density: (C) before batch effect removal and (D) after batch effect removal. The multidimensional scaling of the datasets demonstrates that, before the batch effect adjustment, each dataset obviously separated from all the others, whereas, after batch effect adjustment, samples from all the datasets are incorporated clearly. Supplementary Figure S2: six KEEG pathways that are significantly associated with the downregulated genes in colon tumor stroma relative to normal stroma. Supplementary Figure S3: the genes (DEGs and their upstream regulators) whose expression is significantly associated with disease-free survival in colon cancer (log-rank test, p<0.05). Supplementary Figure S4: patients divided into the high-risk group and the low-risk group based on the prognostic gene signatures identified. A. Overall survival. B. Disease-free survival. Supplementary Figure S5: numbers of overlapping genes between the DEGs between colon tumor stroma and normal stroma and their upstream regulators and the DEGs between high-stroma-content and low-stroma-content colon cancers. UP TCGA: upregulated differentially expressed genes between high-stroma-content and low-stroma-content TCGA colon cancer samples and DOWN TCGA: downregulated differentially expressed genes between high-stroma-content and low-stroma-content TCGA colon cancer samples. Supplementary Table S1: a summary of the datasets used in this study. Supplementary Table S2: upregulated genes in colon tumor stroma versus colon normal stroma. Supplementary Table S3: downregulated genes in colon tumor stroma versus colon normal stroma. Supplementary Table S4: the top 10 upregulated and top 10 downregulated genes in colon tumor stroma. Supplementary Table S5: 44 KEGG pathways that were significantly associated with the upregulated genes in colon tumor stroma (CTS). Supplementary Table S6: upstream transcription factors and kinases regulating the differentially expressed genes between colon tumor stroma and normal stroma. Supplementary Table S7: master metabolic transcriptional regulators (MMTRs) (iRegulon normalized enrichment score NES > 3.0) regulating the differentially expressed genes between colon tumor stroma and normal stroma. Supplementary Table S8: hub genes in the protein-protein interaction network of the differentially expressed genes between colon tumor stroma and normal stroma. Supplementary Table S9: overlapping genes between the DEGs between colon tumor stroma and normal stroma and their upstream regulators and the DEGs between high-stroma-content and low-stroma-content colon cancers. (Supplementary Materials)

References

K. C. Valkenburg, A. E. de Groot, and K. J. Pienta, “Targeting the tumour stroma to improve cancer therapy,” Nature Reviews Clinical Oncology, vol. 15, no. 6, pp. 366–381, 2018.
View at: Publisher Site | Google Scholar
K. M. Bussard, L. Mutkus, K. Stumpf, C. Gomez-Manzano, and F. C. Marini, “Tumor-associated stromal cells as key contributors to the tumor microenvironment,” Breast Cancer Research, vol. 18, no. 1, p. 84, 2016.
View at: Google Scholar
D. Hanahan and L. M. Coussens, “Accessories to the crime: functions of cells recruited to the tumor microenvironment,” Cancer Cell, vol. 21, no. 3, pp. 309–322, 2012.
View at: Publisher Site | Google Scholar
F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: A Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.
View at: Publisher Site | Google Scholar
A. Calon, E. Lonardo, and A. Berenguer-Llergo, “Stromal gene expression defines poor-prognosis subtypes in colorectal cancer,” Nature Genetics, vol. 47, pp. 320–329, 2015.
View at: Publisher Site | Google Scholar
C. Isella, A. Terrasi, S. E. Bellomo et al., “Stromal contribution to the colorectal cancer transcriptome,” Nature Genetics, vol. 47, no. 4, pp. 312–319, 2015.
View at: Publisher Site | Google Scholar
M. Abba, S. Laufs, M. Aghajany, B. Korn, A. Benner, and H. Allgayer, “Look who's talking: Deregulated signaling in colorectal cancer,” Cancer Genomics & Proteomics, vol. 9, no. 1, pp. 15–25, 2012.
View at: Google Scholar
N. Nishida, M. Nagahara, T. Sato et al., “Microarray analysis of colorectal cancer stromal tissue reveals upregulation of two oncogenic miRNA clusters,” Clinical Cancer Research, vol. 18, no. 11, pp. 3054–3070, 2012.
View at: Publisher Site | Google Scholar
“GEO Accession viewer,” https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE46824.
View at: Google Scholar
R. Aguirre-Gamboa, H. Gomez-Rueda, E. Martínez-Ledesma et al., “SurvExpress: an online biomarker validation tool and database for cancer gene expression data using survival analysis,” PLoS ONE, vol. 8, no. 9, p. e74250, 2013.
View at: Publisher Site | Google Scholar
J. Xia, E. E. Gill, and R. E. W. Hancock, “NetworkAnalyst for statistical, visual and network-based meta-analysis of gene expression data,” Nature Protocols, vol. 10, no. 6, pp. 823–844, 2015.
View at: Publisher Site | Google Scholar
W. E. Johnson, C. Li, and A. Rabinovic, “Adjusting batch effects in microarray expression data using empirical Bayes methods,” Biostatistics, vol. 8, no. 1, pp. 118–127, 2007.
View at: Publisher Site | Google Scholar
W. G. Cochran, “The combination of estimates from different experiments,” Biometrics, vol. 10, no. 1, pp. 101–129, 1954.
View at: Google Scholar
Y. Benjamini and Y. Hochberg, “Controlling the false discovery rate: a practical and powerful approach to multiple testing,” Journal of the Royal Statistical Society B: Methodological, vol. 57, no. 1, pp. 289–300, 1995.
View at: Google Scholar | MathSciNet
A. Subramanian, P. Tamayo, V. K. Mootha et al., “Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 102, no. 43, pp. 15545–15550, 2005.
View at: Publisher Site | Google Scholar
D. J. B. Clarke, M. V. Kuleshov, B. M. Schilder et al., “EXpression2Kinases (X2K) Web: Linking expression signatures to upstream cell signaling networks,” Nucleic Acids Research, vol. 46, no. W1, pp. W171–W179, 2018.
View at: Publisher Site | Google Scholar
R. Janky, A. Verfaillie, H. Imrichová et al., “iRegulon: from a gene list to a gene regulatory network using large motif and track collections,” PLoS Computational Biology, vol. 10, no. 7, Article ID e1003731, 2014.
View at: Publisher Site | Google Scholar
D. Petrochilos, A. Shojaie, J. Gennari, and N. Abernethy, “Using random walks to identify cancer-associated modules in expression data,” BioData Mining, vol. 6, p. 17, 2013.
View at: Google Scholar
K. Yoshihara, M. Shahmoradgoli, E. Martínez et al., “Inferring tumour purity and stromal and immune cell admixture from expression data,” Nature Communications, vol. 4, p. 2612, 2013.
View at: Publisher Site | Google Scholar
X.-M. Yang, H.-Y. You, Q. Li et al., “CTHRC1 promotes human colorectal cancer cell proliferation and invasiveness by activating Wnt/PCP signaling,” International Journal of Clinical and Experimental Pathology, vol. 8, no. 10, pp. 12793–12801, 2015.
View at: Google Scholar
B. Guo, H. Yan, L. Li, K. Yin, F. Ji, and S. Zhang, “Collagen triple helix repeat containing 1 (CTHRC1) activates Integrin β3/FAK signaling and promotes metastasis in ovarian cancer,” Journal of Ovarian Research, vol. 10, no. 1, p. 69, 2017.
View at: Google Scholar
A. M. Chowdhury, H. Katoh, A. Hatanaka et al., “Multiple regulatory mechanisms of the biological function of NRF3 (NFE2L3) control cancer cell proliferation,” Scientific Reports, vol. 7, no. 1, p. 12494, 2017.
View at: Publisher Site | Google Scholar
M. Palma, L. Lopez, M. García et al., “Detection of collagen triple helix repeat containing-1 and nuclear factor (erythroid-derived 2)-like 3 in colorectal cancer,” BMC Clinical Pathology, vol. 12, p. 2, 2012.
View at: Publisher Site | Google Scholar
R. M. S. Gill, A. Michael, L. Westley, H. M. Kocher, J. I. Murphy, and G. K. Dhoot, “SULF1/SULF2 splice variants differentially regulate pancreatic tumour growth progression,” Experimental Cell Research, vol. 324, no. 2, pp. 157–171, 2014.
View at: Publisher Site | Google Scholar
H. Lee, B. Yeh, T. Chan et al., “Sulfatase-1 overexpression indicates poor prognosis in urothelial carcinoma of the urinary bladder and upper tract,” Oncotarget, vol. 8, no. 29, pp. 47216–47229, 2017.
View at: Publisher Site | Google Scholar
B. Lü, Y. Fang, J. Xu et al., “Analysis of SOX9 expression in colorectal cancer,” American Journal of Clinical Pathology, vol. 130, no. 6, pp. 897–904, 2008.
View at: Publisher Site | Google Scholar
S. Arikan, C. Cacina, E. Guler, S. Çulcu, G. Tuna, and I. Yaylim-Eraltan, “The effects of NOS3 Glu298Asp variant on colorectal cancer risk and progression in Turkish population,” Molecular Biology Reports, vol. 39, no. 3, pp. 3245–3249, 2012.
View at: Publisher Site | Google Scholar
S. Chida, H. Okayama, M. Noda et al., “Stromal VCAN expression as a potential prognostic biomarker for disease recurrence in stage II-III colon cancer,” Carcinogenesis, vol. 37, no. 9, pp. 878–887, 2016.
View at: Publisher Site | Google Scholar
S. Demirkol, I. Gomceli, M. Isbilen et al., “A combined ULBP2 and SEMA5A expression signature as a prognostic and predictive biomarker for colon cancer,” Journal of Cancer, vol. 8, no. 7, pp. 1113–1122, 2017.
View at: Publisher Site | Google Scholar
M. Fujita, Y. Furukawa, T. Tsunoda, T. Tanaka, M. Ogawa, and Y. Nakamura, “Up-regulation of the ectodermal-neural cortex 1 (ENC1) gene, a downstream target of the β-catenin/T-cell factor complex, in colorectal carcinomas,” Cancer Research, vol. 61, no. 21, pp. 7722–7726, 2001.
View at: Google Scholar
Y. Li, J. Wei, C. Xu, Z. Zhao, and T. You, “Prognostic significance of cyclin D1 expression in colorectal cancer: a meta-analysis of observational studies,” PLoS ONE, vol. 9, no. 4, Article ID e94508, 2014.
View at: Publisher Site | Google Scholar
X. Hu, X. Sui, L. Li et al., “Protocadherin 17 acts as a tumour suppressor inducing tumour cell apoptosis and autophagy, and is frequently methylated in gastric and colorectal cancers,” The Journal of Pathology, vol. 229, no. 1, pp. 62–73, 2013.
View at: Publisher Site | Google Scholar
Y. Gu, A. Li, H. Sun et al., “BCL6B suppresses proliferation and migration of colorectal carcinoma cells through inhibition of the PI3K/AKT signaling pathway,” International Journal of Molecular Medicine, vol. 41, no. 5, pp. 2660–2668, 2018.
View at: Publisher Site | Google Scholar
Y. Guo, Y. Bao, M. Ma, and W. Yang, “Identification of key candidate genes and pathways in colorectal cancer by integrated bioinformatical analysis,” International Journal of Molecular Sciences, vol. 18, no. 4, 2017.
View at: Publisher Site | Google Scholar
Q. Sun, M. Li, and X. Wang, “The Cancer Omics Atlas: An integrative resource for cancer omics annotations,” BMC Medical Genomics, vol. 11, no. 1, p. 63, 2018.
View at: Google Scholar
Q. Huang, S. Li, P. Cheng et al., “High expression of anti-apoptotic protein Bcl-2 is a good prognostic factor in colorectal cancer: result of a meta-analysis,” World Journal of Gastroenterology, vol. 23, no. 27, pp. 5018–5033, 2017.
View at: Publisher Site | Google Scholar
S.-H. Yun, M.-S. Roh, J.-S. Jeong, and J.-I. Park, “Peroxisome proliferator-activated receptor γ coactivator-1α is a predictor of lymph node metastasis and poor prognosis in human colorectal cancer,” Annals of Diagnostic Pathology, vol. 33, pp. 11–16, 2018.
View at: Publisher Site | Google Scholar
F. E. Bertrand, C. W. Angus, W. J. Partis, and G. Sigounas, “Developmental pathways in colon cancer,” Cell Cycle, vol. 11, no. 23, pp. 4344–4351, 2012.
View at: Publisher Site | Google Scholar
A. M. Socovich and A. Naba, “The cancer matrisome: From comprehensive characterization to biomarker discovery,” Seminars in Cell & Developmental Biology, 2018.
View at: Google Scholar
N. Afratis, C. Gialeli, D. Nikitovic et al., “Glycosaminoglycans: key players in cancer cell biology and treatment,” FEBS Journal, vol. 279, pp. 1177–1197, 2012.
View at: Publisher Site | Google Scholar
Basappa, K. S. Rangappa, and K. Sugahara, “Roles of glycosaminoglycans and glycanmimetics in tumor progression and metastasis,” Glycoconjugate Journal, vol. 31, no. 6-7, pp. 461–467, 2014.
View at: Publisher Site | Google Scholar
J. Yin, W. Ren, X. Huang, J. Deng, T. Li, and Y. Yin, “Potential mechanisms connecting purine metabolism and cancer therapy,” Frontiers in Immunology, vol. 9, 2018.
View at: Publisher Site | Google Scholar
C. Frezza, “Histidine metabolism boosts cancer therapy,” Nature, vol. 559, no. 7715, pp. 484-485, 2018.
View at: Publisher Site | Google Scholar
M. Kumarakulasingham, P. H. Rooney, S. R. Dundas et al., “Cytochrome P450 profile of colorectal cancer: Identification of markers of prognosis,” Clinical Cancer Research, vol. 11, no. 10, pp. 3758–3765, 2005.
View at: Publisher Site | Google Scholar
Y.-L. Liu, X. Gao, Y. Jiang et al., “Expression and clinicopathological significance of EED, SUZ12 and EZH2 mRNA in colorectal cancer,” Journal of Cancer Research and Clinical Oncology, vol. 141, no. 4, pp. 661–669, 2015.
View at: Publisher Site | Google Scholar
M. L. Slattery, A. Lundgreen, J. S. Herrick, B. J. Caan, J. D. Potter, and R. K. Wolff, “Associations between genetic variation in RUNX1, RUNX2, RUNX3, MAPK1 and eIF4E and risk of colon and rectal cancer: additional support for a TGF-β-signaling pathway,” Carcinogenesis, vol. 32, no. 3, pp. 318–326, 2011.
View at: Publisher Site | Google Scholar
K. Taniue, A. Kurimoto, Y. Takeda et al., “ASBEL –TCF3 complex is required for the tumorigenicity of colorectal cancer cells,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 113, no. 45, pp. 12739–12744, 2016.
View at: Publisher Site | Google Scholar
T. Bae, K. Rho, J. W. Choi, K. Horimoto, W. Kim, and S. Kim, “Identification of upstream regulators for prognostic expression signature genes in colorectal cancer,” BMC Systems Biology, vol. 7, no. 1, p. 86, 2013.
View at: Publisher Site | Google Scholar
T. J. Miller, M. J. McCoy, C. Hemmings, M. K. Bulsara, B. Iacopetta, and C. F. Platell, “The prognostic value of cancer stem-like cell markers SOX2 and CD133 in stage III colon cancer is modified by expression of the immune-related markers FoxP3, PD-L1 and CD3,” Pathology, vol. 49, no. 7, pp. 721–730, 2017.
View at: Publisher Site | Google Scholar
J. Gupta, I. delBarcoBarrantes, A. Igea et al., “Dual function of p38α MAPK in colon cancer: suppression of colitis-associated tumor initiation but requirement for cancer cell survival,” Cancer Cell, vol. 25, no. 4, pp. 484–500, 2014.
View at: Publisher Site | Google Scholar
S. R. Rosario, M. D. Long, H. C. Affronti, A. M. Rowsam, K. H. Eng, and D. J. Smiraglia, “Pan-cancer analysis of transcriptional metabolic dysregulation using the cancer genome atlas,” Nature Communications, vol. 9, 2018.
View at: Publisher Site | Google Scholar
D. Liu, X.-X. Zhang, M.-C. Li et al., “C/EBPβ enhances platinum resistance of ovarian cancer cells by reprogramming H3K79 methylation,” Nature Communications, vol. 9, no. 1, p. 1739, 2018.
View at: Google Scholar
J. Feilchenfeldt, M.-A. Bründler, C. Soravia, M. Tötsch, and C. A. Meier, “Peroxisome proliferator-activated receptors (PPARs) and associated transcription factors in colon cancer: reduced expression of PPARgamma-coactivator 1 (PGC-1),” Cancer Letters, vol. 203, no. 1, pp. 25–33, 2014.
View at: Publisher Site | Google Scholar
Y. Zhang, Y. Ba, C. Liu et al., “PGC-1alpha induces apoptosis in human epithelial ovarian cancer cells through a PPARgamma-dependent pathway,” Cell Research, vol. 17, no. 4, pp. 363–373, 2007.
View at: Google Scholar
J. Concetti and C. L. Wilson, “NFKB1 and cancer: friend or foe?” Cells, vol. 7, no. 9, 2018.
View at: Publisher Site | Google Scholar
U. Saini, S. Naidu, A. C. Elnaggar et al., “Elevated STAT3 expression in ovarian cancer ascites promotes invasion and metastasis: A potential therapeutic target,” Oncogene, vol. 36, no. 2, pp. 168–181, 2017.
View at: Publisher Site | Google Scholar
W. J. Yang, Y. Qu, B. X. Tan et al., “Prognostic significance of preoperative IKBKE expression in esophageal squamous cell carcinoma,” OncoTargets and Therapy, vol. 11, pp. 1305–1314, 2018.
View at: Publisher Site | Google Scholar
W. Sung, Y. Lin, P. Wu et al., “High nuclear/cytoplasmic ratio of Cdk1 expression predicts poor prognosis in colorectal cancer patients,” BMC Cancer, vol. 14, p. 951, 2014.
View at: Publisher Site | Google Scholar
K. Zhuang, J. Zhang, M. Xiong et al., “CDK5 functions as a tumor promoter in human colorectal cancer via modulating the ERK5–AP-1 axis,” Cell Death & Disease, vol. 7, no. 10, p. e2415, 2016.
View at: Publisher Site | Google Scholar
Y. Hu, J. Zhou, F. Ye et al., “BRD4 inhibitor inhibits colorectal cancer growth and metastasis,” International Journal of Molecular Sciences, vol. 16, no. 1, pp. 1928–1948, 2015.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Md. Nazim Uddin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

2381

Downloads

1446

Citations