- Letter to the Editor
- Open Access
Cancer-associated fibroblast-derived gene signatures determine prognosis in colon cancer patients
Molecular Cancer volume 20, Article number: 73 (2021)
Paracrine communication between tumor and surrounding stroma arbitrates the malignant behavior of cancer progression . Fibroblasts, which are the main cell type within the stroma and are called cancer-associated fibroblasts (CAFs), orchestrate the crosstalk with cancer cells [2, 3] and express several markers associated with prognosis . There is increasing evidence that a stroma-specific signature could be used for risk assessment in colon cancer (CC). According to the Consensus Molecular Subtype classification (CMS) in CC, the mesenchymal or CMS4 group is characterized by stromal invasion, extracellular matrix remodeling and TGF-β signaling activation. It is associated with the worst prognosis rates [5, 6]. Genes correlating with the mesenchymal subtype are mostly expressed by CAFs and other stromal cells, rather than by tumor cells . Accordingly, our group defined a gene expression profile associated with CAFs with high pro-migratory effects on colon tumor cells, which was associated with patients’ poor prognosis. These were mostly advanced-stage patients .
The crosstalk between tumor and stromal cells is conducted in part by exosomes that are involved in many tumorogenic processes. ncRNAs contained in exosomes secreted by colon CAFs enhance proliferation and stemness properties of tumoral cells and are involved in chemoresistance . Recently, our group demonstrated that there was a difference in how CAFs and normal colon mucosa fibroblasts (NFs) distributed ncRNA into the exosomal cargo. There was the same difference in distribution for their potential ncRNA target genes in recipient cells .
Our research into CAF profiles and targets of CAF-derived exosome cargo in recipient cells found novel CAF-derived signatures with prognostic value in colon cancer patients. This showed the importance of CAFs and their derived exosomes in tumor progression (Fig. 1).
CAF versus normal fibroblasts’ gene signature as a marker of risk and survival in colon cancer patients
A “CAF signature”, involving 596 protein-coding genes, was identified when gene expression data from colon tumor CAFs and NFs were compared (Additional files 1 and 2). Functional gene-set analysis identified catabolic process, intracellular transport and dimerization and binding of proteins, among others, as enriched biological terms (Additional file 3). The prognostic value and clinical relevance of the new “CAF signature” were evaluated in a meta-dataset of 1235 colon cancer patients. Figure 2a shows the results of the risk prediction for each of the 1235 colon tumors and the best cut-off to divide patients into lower/higher “signature gene score”, based on expression levels of selected genes. Interestingly, the survival analysis showed that a high “signature gene score” subset of patients had shorter overall survival (OS) than the low “signature gene score” patients (Fig. 2b). In line with these data, the prognostic value of our previously reported “pro-migratory CAF signature”  was also confirmed in this cohort of 1235 patients (data not shown).
Prognostic value of a novel CAF-derived exosomes signature
The ncRNA consists of functional RNA molecules that act by means of multiple mechanisms in the regulation of diverse cell functions . In this study, we analyzed a second signature, “CAF-derived exosomes signature”, containing 67 predicted target genes of over-distributed ncRNA genes in CAF-derived exosomes previously described by our group  (Additional file 4). Gene set functional enrichment analysis identified important pathways related to cancer progression and microenvironment regulation, such as proliferation, positive regulation of signaling, regulation of multicellular organismal process and DNA repair processes, among others (Additional file 5). The risk cut-off point (Fig. 2c) of the “CAF-derived exosomes signature” in colon cancer was also evaluated in the same large meta-dataset (N = 1235), confirming the signature’s prognostic value (Fig. 2d).
Overlap of CAF-related signatures and prognostic value, using their combination
Our first signature is related to CAF activation or status, while the second is associated with the crosstalk between fibroblasts and stromal/cancer cells orchestrated by CAF-derived exosomes. Both CAF-related signatures clearly overlapped in the classification of patients as at high or low risk of death (82.6% of concordance). To achieve a better and more balanced comparison of the behavior of the two CAF-related signatures, we undertook a new analysis of only the top 50 genes of each signature (Additional file 6), which confirmed their prognostic value (Fig. 2e and f). The combination of the two signatures (50 + 50 genes) did not greatly improve the separation of the survival curves of low- and high-“signature gene score” patients (Fig. 2g), indicating that there was no overlapping gene between the two gene signatures. Therefore, they include features that are closely related or measure similar characteristics in the colon tumors, though from two perspectives: CAF gene deregulation and fibroblasts/other cells’ crosstalk.
Prognostic value of CAF-related signatures increases in advanced-stage tumors and in CMS4 tumors
Initial analysis of survival of our cohort of patients confirmed the prognostic value of tumor stage (Additional file 7A). Interestingly, when we calculated the risk of each individual by splitting the population into stages (Fig. 3a), we found that all CAF-related signatures predicted higher risk at advanced stages than at early stages. In parallel, the prognostic value of CAF signatures was stronger in the stage III-IV group of patients than in the group of patients at early tumor stages (Fig. 3b-c and Additional file 8A). Our previous studies corroborated these data and showed the greater impact of microenvironment activation in patients’ outcome at advanced CC stages [8, 12]. Since these patients are usually those with worse outcomes, new therapies targeting microenvironment components would improve their clinical management. Specifically, stage III colon cancer patients are likely to develop recurrence. However, a subgroup of stage III patients would have low risk of recurrence and better outcome, with the advantage of shorter chemotherapy treatments that would avoid oxaliplatin cumulative neurotoxicity . However, in daily clinical practice there is not enough evidence for better risk stratification in high/low-risk stage III patients. Our defined CAF-derived signatures with a clear prognostic value might improve patients’ classification and thus cancer patients’ management.
In addition, we confirmed that CMS4 patients had lower survival than CMS1, 2 or 3 ones (Additional file 7B). The Wilcoxon test revealed a higher “signature gene score” predicted by any of the 3 CAF-related signatures in CMS4-type tumors than in CMS1–3 ones (Fig. 3d). Although the Kaplan-Meier survival analysis showed that our CAF-related signatures marked differences in both CMS1–3 and CMS4 tumors (Fig. 3e-f and Additional file 8B), the comparison of the “signature gene score” distributions clearly indicated that in CMS4 risks were higher. CMS4-like subtypes and treatment associations have been studied in some clinical and preclinical studies. However, the clinical translation of the CMS subtypes into treatment decision-making is not widespread . As our defined CAF-derived signatures had a higher prognostic value in those patients in the CMS4-like subtype, this might provide a new framework for patient classification in clinical trials and for the adjustment of suitable treatments. Therapies with the tumor microenvironment as the main target are currently an area of intensive research .
In summary, clinical decisions are based on the usual histological and clinical parameters that do not accurately predict the biological behavior of histologically equal tumors. The risk of recurrence of CC cancer is a critical factor in correctly setting up therapeutic guidelines, which highlights the importance of studying prognosis and predictive biomarkers. These new biomarkers classify patients on the basis of tumor microenvironment differences and define those pathological events at the molecular and cell level associated with patients’ survival. Our results emphasize ongoing efforts to decipher CAF biomarkers in order to advance prognosis and chemotherapeutic responses’ heterogeneity and to improve patients’ survival by tailoring therapies based on specific biological characteristics of CC tumors.
Availability of data and materials
All the data corresponding to the colorectal cancer series used in this study are available in GEO (https://www.ncbi.nlm.nih.gov/geo), which is a public functional genomics data repository. The source identifier (GEO GSM id) of each sample used is indicated in the corresponding Additional File.
Consensus Molecular Subtypes
Gene Ontology-Biological Process
Gene Ontology-Molecular Function
Next Generation Sequencing
Pietras K, Östman A. Hallmarks of cancer: interactions with the tumor stroma. Exp Cell Res. 2010;316(8):1324–31. https://doi.org/10.1016/j.yexcr.2010.02.045.
Xing F, Saidou J, Watabe K. Cancer associated fibroblasts (CAFs) in tumor microenvironment. Front Biosci. 2010;15(1):166–79. https://doi.org/10.2741/3613.
Cirri P, Chiarugi P. Cancer associated fibroblasts: the dark side of the coin. Am J Cancer Res. 2011;1(4):482–97.
Paulsson J, Micke P. Prognostic relevance of cancer-associated fibroblasts in human cancer. Semin Cancer Biol. 2014;25:61–8. https://doi.org/10.1016/j.semcancer.2014.02.006.
Guinney J, Dienstmann R, Wang X, De Reyniès A, Schlicker A, Soneson C, et al. The consensus molecular subtypes of colorectal cancer. Nat Med. 2015;21(11):1350–6. https://doi.org/10.1038/nm.3967.
Conde Moreno E, Pascual A, Prieto-Cuadra D, Laza VF, Molina-Cerrillo J, Ramos-Muñoz ME, et al. Novel Molecular Characterization of Colorectal Primary Tumors Based on miRNAs. Cancers (Basel). 2019;11:346.
Isella C, Terrasi A, Bellomo SE, Petti C, Galatola G, Muratore A, et al. Stromal contribution to the colorectal cancer transcriptome. Nat Genet. 2015;47(4):312–9. https://doi.org/10.1038/ng.3224.
Herrera M, Islam ABMMK, Herrera A, Martín P, García V, Silva J, et al. Functional heterogeneity of cancer-associated fibroblasts from human colon tumors shows specific prognostic gene expression signature. Clin Cancer Res. 2013;19(21):5914–26. https://doi.org/10.1158/1078-0432.CCR-13-0694.
Herrera M, Galindo-Pumariño C, García-Barberán V, Peña C. A snapshot of the tumor microenvironment in colorectal cancer: the liquid biopsy. Int J Mol Sci. 2019;20(23):6016. https://doi.org/10.3390/ijms20236016.
Herrera M, Llorens C, Rodríguez M, Herrera A, Ramos R, Gil B, et al. Differential distribution and enrichment of non-coding RNAs in exosomes from normal and Cancer-associated fibroblasts in colorectal cancer. Mol Cancer. 2018;17(1):114. https://doi.org/10.1186/s12943-018-0863-4.
Szymański M, Barciszewska MZ, Zywicki M, Barciszewski J. Noncoding RNA transcripts. J Appl Genet. 2003;44(1):1–19.
Herrera M, Herrera A, Domínguez G, Silva J, García V, García JM, et al. Cancer-associated fibroblast and M2 macrophage markers together predict outcome in colorectal cancer patients. Cancer Sci. 2013;104(4):437–44. https://doi.org/10.1111/cas.12096.
Grothey A, Sobrero AF, Shields AF, Yoshino T, Paul J, Taieb J, et al. Duration of adjuvant chemotherapy for stage III colon cancer. N Engl J Med. 2018;378(13):1177–88. https://doi.org/10.1056/NEJMoa1713709.
Sveen A, Kopetz S, Lothe RA. Biomarker-guided therapy for colorectal cancer: strength in complexity. Nat Rev Clin Oncol. 2020;17(1):11–32. https://doi.org/10.1038/s41571-019-0241-1.
Sahai E, Astsaturov I, Cukierman E, DeNardo DG, Egeblad M, Evans RM, et al. A framework for advancing our understanding of cancer-associated fibroblasts. Nat Rev Cancer. 2020;20(3):174–86. https://doi.org/10.1038/s41568-019-0238-1.
M. Eaude helped with the English text. We are grateful to lab members for help and advice throughout this research.
This research is supported by PI17/01847, PI18/01020 and PI20/00602 from the Instituto de Salud Carlos III and co-financed by the European Development Regional Fund (FEDER) “A way to achieve Europe” (ERDF); by “CIBER de Cáncer”, CB16/12/00273, CB16/12/00301 and CB16/12/00446, from the Instituto de Salud Carlos III-FEDER; and by the Fundación Científica AECC (a multifaceted approach to targeting pancreatic cancer). The JDLR group also acknowledges the funding provided by the Instituto de Salud Carlos III (ISCiii, AES) in grants PI18/00591 and PT17/0009/0008, co-financed by the European Development Regional Fund (FEDER).
Ethics approval and consent to participate
Not applicable. All the data corresponding to the colorectal cancer series used in this study are available as public functional genomics data repository.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Material and Methods; Clinical and phenotypic characteristics of the colorectal cancer cohort of 1273 patients analyzed in this study.
List of 596 protein-coding genes defined by the “CAF signature”, including: gene symbol; gene description; ENSEMBL gene id; and statistical parameters provided by the SAM algorithm (i.e. fold change in log2 scale, statistic parameter d value, and p value) in the analysis of differential expression.
Functional enrichment analysis with GO-BP and GO-MF of the genes of the “CAF signature” within the database MSigDB. The table shows the most relevant GO terms (top 20 in GO-BP and in GO-MF) associated with genes from the CAF-derived signature.
List of the 67 genes included in the “CAF-derived exosomes signature”. For each gene the table gives the gene symbol and the ENSEMBL gene id (ENSG).
Functional enrichment analysis with GO-BP and GO-MF of the genes of the “CAF-derived exosomes signature” within the database MSigDB. The table shows the most relevant GO terms (top 20 in GO-BP and in GO-MF) associated with genes from this signature.
Top 50 genes of each signature (CAF-signature in blue and CAF-derived exosomes signature in green) and the beta-values corresponding to each gene provided by Cox. * non-coding genes.
Survival analysis of colon cancer cohort including phenotypic factors: tumor stage and CMS. (A) Distribution of the patients at stages I & II versus III & IV clearly showed different prognoses. (B) CMS4 patients showed shorter survival than CMS1, CMS2 and CMS3 ones.
CAF-related signatures showed stronger poor prognosis association in those patients with advanced colon tumors than in initial stages and in those patients with CMS4 than in CMS1–3. Comparison of the prognosis of the colon tumors (1235 samples) in initial stages I and II (660 samples) and in advanced stages III and IV (575 samples) (A), divided into CMS1–3 (593 samples) and CMS4 (237 samples) subtypes (B). The survival analysis used 2 gene signatures: the top 50 genes of the 67-gene signature (“CAF-derived exosomes signature”) and the top 50 of the 596-gene signature (“CAF signature”).
About this article
Cite this article
Herrera, M., Berral-González, A., López-Cade, I. et al. Cancer-associated fibroblast-derived gene signatures determine prognosis in colon cancer patients. Mol Cancer 20, 73 (2021). https://doi.org/10.1186/s12943-021-01367-x
- Colon Cancer
- Tumor microenvironment
- Cancer-associated fibroblasts
- Noncoding RNAs
- Prognostic signatures
- Liquid biopsy