Skip to main content

Cancer-associated fibroblasts impact the clinical outcome and treatment response in colorectal cancer via immune system modulation: a comprehensive genome-wide analysis

Abstract

Background

Cancer-associated fibroblasts (CAFs) in the tumour microenvironment are associated with poor prognosis and chemoresistance in multiple solid tumours. However, there is a lack of universal measures of CAFs in colorectal cancer (CRC). The aim of this study was to assess the utility of a fibroblast-related gene signature (FRGS) for predicting patient outcomes and reveal its relevant mechanism.

Methods

The GSE39582 dataset, which includes 316 CRC patients who did not receive adjuvant chemotherapy was used as a discovery cohort to identify the prognostic fibroblast-related genes (FRGs). A total of 1352 CRC patients were divided into one training cohort (GSE39582, n = 461) and two validation cohorts (TCGA, n = 338; meta-validation, n = 553) for the construction of the FRGS and the verification of its prognostic value in stage II/III CRC patients. Functional annotation and analysis were performed to explore the underlying mechanism. The ability of the FRGS to predict immunotherapy response was further tested in a clear cell renal cell carcinoma (ccRCC) cohort.

Results

An 11-gene signature that had prognostic value for stage II/III CRC patients in both validation cohorts was developed (TCGA cohort: HR = 1.90, 95% CI 1.16–3.12, P < 0.01; meta-validation cohort: HR = 1.95, 95% CI 1.39–2.73, P < 0.001). A high level of CAFs was correlated with worse prognosis in CRC patients who did not receive adjuvant chemotherapy (HR = 3.63, 95% CI 2.24–5.88, P < 0.001). Importantly, patients in the low-risk group were found to be benefit from chemotherapy (P < 0.01), but not in the high CAF group (P > 0.05). Similar results were found in the TCGA cohort. Integrated with clinical characteristics, the FRGS was confirmed to be an independent prognostic factor in the multivariate analysis after adjustment for tumour TNM stage (GSE39582 cohort: HR = 3.19, 95% CI 1.88–5.41, P < 0.001; TCGA cohort: HR = 5.00, 95% CI 1.58–15.85, P = 0.007; meta-validation cohort: HR = 2.99, 95% CI 1.44–6.21, P = 0.003). Furthermore, the enrichment analysis found that the antitumour immune response was suppressed and the infiltration of CD4 T cells and M1 macrophages was depressed in the high CAF group. The FRGS was also found to have value in predicting for immunotherapy response in the ccRCC cohort.

Conclusions

The 11-gene FRGS had independent prognostic value for CRC patients, as well as utility in the prediction of benefit from chemotherapy. CAFs in the tumour microenvironment might have an impact on the prognosis of CRC patients via inhibiting immune response.

Introduction

Colorectal cancer (CRC) is the third most common and the second most lethal malignancy globally, representing approximately 10% of overall cancer cases and deaths (Sung et al. 2021). Its biological heterogeneity leads to diversity in the survival of CRC patients. Approximately 25% of patients with stage II and III CRC suffer from tumor recurrence despite receiving radical surgery and adjuvant chemotherapy. Moreover, the incidences of chemotherapy-associated death and significant side-effects are approximately 0.5–1% and 20%, respectively (Andre et al. 2009, 2015; Quasar Collaborative Group 2000; Quasar Collaborative Group et al. 2007). Therefore, a new strategy based on reliable markers to stratify subgroups with different risks of treatment response, tumour recurrence and tumour-specific death based on reliable markers is an urgent need for precise treatment of CRC patients.

Cancer-associated fibroblasts (CAFs) are considered to be one of the critical components of the tumour microenvironment (TME) in CRC because they provide physical support for epithelial cells and functioning as key regulators in tumorigenesis in a context-dependent manner (Kobayashi et al. 2019). The infiltration of CAFs has been reported to correlate with CRC patient prognosis and drug response. Previous studies have revealed that CAFs contribute to an immunosuppressive TME and promote tumour invasion and metastasis via the secretion of different kinds of cytokines and chemokines, such as IL-6 and CCL2 (Kalluri 2016). CAFs are also reported to be mediators of response to immune checkpoint inhibitors (Tauriello et al. 2018). Furthermore, emerging evidence has demonstrated that CAFs might confer substantial resistance to cancer therapeutics by impairing drug delivery and several biochemical signaling pathways (Kalluri 2016). Our previous study showed that CAFs could induce the chemoresistance of CRC through enhancement of the stemness of cancer cells (Tang et al. 2018). Additionally, our study showed that characteristics of the TME might be able to predict patient clinical outcome (Zou et al. 2019). The above studies imply that CRC patients might benefit from improved selection based on individual fibroblast-related characteristics. However, there is no clinical method of selecting patients with high levels of fibroblast-related risk factors.

Therefore, we analyzed the expression of fibroblast-related genes (FRGs) in CRC transcriptional data in this study, and combined multiple fibroblastic genes to construct a prognosis-related signature. Furthermore, the prognostic prediction value of this fibroblast-related gene signature (FRGS) was validated systematically. This study will aid in therapeutic decision making for CRC patients.

Materials and methods

Patient cohorts

Six public cohorts with gene expression data derived from fresh-frozen CRC samples were evaluated retrospectively, including The Cancer Genome Atlas (TCGA) CRC cohort and five datasets from the Gene Expression Omnibus (GEO) database (GSE39582, GSE14333, GSE17536, GSE37892 and GSE33113). The GSE39582 dataset including 309 CRC patients who did not receive adjuvant chemotherapy was used as the training cohort, while the TCGA CRC cohort and a meta-validation cohort containing the other four datasets were used for independent validation. The gene expression profiles of TCGA CRC cohort were obtained from Broad GDAC Firehose (http://gdac.broadinstitute.org/), as transcripts per million (TPM) values of level-three RNA-seq data on a log2 scale were used in our previous study (Zou et al. 2019). The other five datasets were directly obtained from the GEO database with the Bioconductor package ‘GEOquery’. A total of 1656 CRC patients were ultimately included in this study. The batch effects were corrected with the ‘combat’ algorithm implemented in the R package ‘sva’, and the z-scores for each gene were used during the following analyses. Data were collected from Sep 27 to Dec 26, 2018.

Construction and validation of the CAF-related gene signature

To construct a prognostic fibroblast-related gene signature (FRGS), a list of fibroblast-related genes was first identified by including all gene sets using the keyword ‘fibroblast’ in MSigDB (version 6.2), and genes that were measured in each platform were selected. Prognostic genes were further selected using the log-rank test with 1000 randomizations (95% portion of samples each time) to assess the correlation between each gene and patients disease-free survival (DFS) in GSE39582 dataset. CRC patients in all stages were included in this stage. The genes that repeatedly showed significance were selected as the candidates for the FRGS. To minimize the risk of overfitting, a Cox proportional hazards regression model combined with the least absolute shrinkage and selection operator (LASSO) method (glmnet, version 2.0-16) was employed. The penalty parameter was estimated by tenfold cross-validation in the training dataset at one SE beyond the minimum partial likelihood deviance.

To separate patients into low- and high-risk groups, the optimal cut-off of each gene was determined with the time-dependent receiver operating characteristic (ROC) curve (survival ROC, version 1.0.3) for 5-year DFS in the training dataset. The Kaplan–Meier method was used to estimate survival. Then, the prognostic value of the FRGS was assessed in stage II/III CRC patients in the training and independent validation cohorts with univariate and multivariate analyses. A P value less than 0.05 was considered statistically significant.

Functional annotation and analysis

To explore the biological characteristics of the FRGS, enrichment analysis of the genome-wide differentially expressed genes between different CAF risk groups was performed with the R package ‘gProfileR’ with the TCGA CRC dataset. To clarify the potential biological pathways, gene set enrichment analysis (GSEA) was conducted with the Bioconductor package ‘HTSanalyzeR’ (Wang et al. 2011). The Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm was further used to estimate the proportions of stromal and immune cells. The levels of different infiltrating immune cells, such as monocytes, lymphocytes and neutrophils, were calculated with CIBERSORT.

Evaluation of the ability of the FRGS to predict the response to PD-1 treatment

To investigate the value of the FRGS for predicting the response to immune checkpoint inhibitors, a cohort with patients with advanced clear cell renal cell carcinoma who were mostly treated with PD-1 targeted treatment was used, and the patients were divided into a high-risk group and a low-risk group. The progression-free survival (PFS) and overall survival (OS) of the different groups were compared using the Kaplan–Meier method.

Statistical analysis

Statistical analyses were performed with SPSS (version 22.0.0, IBM SPSS statistics, IBM Corporation, Armonk, NY) and R software (version 3.5.1; http://www.Rproject.org). Descriptive statistics, including means and standard deviations (SD) or medians and interquartile ranges (IQR) for continuous factors and frequencies for categorical factors, were computed for all variables. Continuous values were compared between different groups using Student’s t tests. Univariate analysis of the association of the FRGS and other clinicopathological factors with survival was performed using the log-rank test. Factors that were significantly associated with DFS and OS in univariate analyses were included in the Cox proportional hazards regression multivariate analysis. The C-index was calculated with the ‘survcomp’ package (version 1.32.0). A P value less than 0.05 was considered as statistically significant in all tests.

Results

Establishment and training of the FRGS in GSE39582 cohort

The GSE39582 dataset included a total of 309 eligible CRC patients who were not treated with adjuvant chemotherapy, and these patients were enrolled in the analysis as the discovery cohort (Fig. 1A). Among 1531 FRGs from MSigDB v6.2, 1472 FRGs were measured on all platforms and thus were included in this study. After 1000 rounds of Cox univariate regression, 80% repeatable genes were chosen, and 76 FRGs were found to be related to patient DFS. After LASSO Cox regression analysis of stage II/III patients, 11 prognostic FRGs were ultimately confirmed and used to construct the FRGS (Fig. 1B). A satisfactory 5-year DFS cut-off identified via in time-dependent ROC curve analysis was used to train the FRGS for stratification of the high- and low-risk groups (Fig. 1C).

Fig. 1
figure1

Establishment and verification of FRGS. A Schematic flow chart of study design. B A total of 11 fibroblast-related genes selected in LASSO Cox regression. C The optimal cut-off obtained at 5-year in time-dependent ROC curve analysis

The correlation coefficients between the genes in the FRGS and prognostic indexed were obtained (Table 1), and a risk score calculation model was used (Additional file 1: Fig. S1). Stage II/III CRC patients in the GSE39582 dataset (n = 461) were then used as the training cohort. More patients with tumour recurrence were found in the high-risk group than in the low-risk group in both the training and the two validation cohorts (Fig. 2A–C).

Table 1 The list of 11-Gene fibroblastic signatures
Fig. 2
figure2

The outcome of different CAF risk in stage II/III CRC patients. AC The DFS of patients with different CAF risk group in training cohort (A), TCGA cohort (B) and meta-validation cohort (C). DF Kaplan–Meier curves comparing patients with different risk in training cohort (D), TCGA cohort (E) and meta-validation cohort (F)

Validation of the FRGS in the TCGA and meta-validation cohorts

The other five CRC transcriptional datasets with prognostic data, including the TCGA dataset (n = 338) and the meta-validation cohort (n = 553) (combination of GSE17536, GSE33113, GSE37892, and GSE14333), were used as validation cohorts for the evaluation of the prognostic value of the FRGS. The values of the FRGS for patients in all stages were shown in Data S1 (Additional file 2: Data S1). No significant difference was observed among the three cohorts regarding the clinical and pathologic factors (date not shown).

Focusing on the stage II/III patients, the low- and high-risk groups stratified by the same risk score calculation model showed a significant difference in DFS in the training cohort (HR = 2.52, 95% CI 1.77–3.57, P < 0.001) (Fig. 2D), TCGA validation cohort (HR = 1.90, 95% CI 1.16–3.12, P < 0.01) (Fig. 2E) and meta-validation cohort (HR = 1.95, 95% CI 1.39–2.73, P < 0.001) (Fig. 2F). Furthermore, the FRGS showed satisfactory prognostic value in all CRC cohorts regardless of stage (CSE39582 cohort, HR = 1.91, 95% CI 1.41–2.58, P < 0.001; TCGA cohort, HR = 2.06, 95% CI 1.37–3.09, P < 0.001; meta-validation cohort, HR = 1.95, 95% CI 1.42–2.66, P < 0.001).

The FRGS predicts benefit from chemotherapy in stage II/III CRC patients

Univariate and multivariate analyses were employed to investigate the association between sex, age, tumour location, tumour stage, gene aberration status, and FRGS score and CRC patient outcomes. In the univariate analysis, high risk identified by the FRGS was related to significantly worse prognosis in the training cohort (HR = 3.33, 95% CI 1.95–5.68, P < 0.001), TCGA cohort (HR = 5.00, 95% CI 1.58–15.85, P < 0.01) and meta-validation cohort (HR = 2.99, 95% CI 1.44–6.21, P < 0.01) (Table 2). Similarly, FRGS score remained an independent prognostic factor in the multivariate analysis (CSE39582 cohort, HR = 3.19, 95% CI 1.88–5.41, P < 0.001; TCGA cohort, HR = 5.00, 95% CI 1.58–15.85, P < 0.01; meta-validation cohort, HR = 2.99, 95% CI 1.44–6.21, P < 0.01) (Table 2).

Table 2 Univariate and multivariate analysis of FRGS, clinical and pathologic factors with DFS of stage II/III patients in training and validation cohorts

To further clarify whether the FRGS could predict benefit from adjuvant chemotherapy in CRC patients, we focused on the clinical outcome in the nonchemotherapy and chemotherapy groups. In CRC patients who did not receive adjuvant chemotherapy, the DFS of high-risk group was worse than that of the low-risk group in both the training (HR = 3.74, 95% CI 2.22–6.28, P < 0.001, Fig. 3A) and TCGA cohorts (HR = 2.61, 95% CI 1.22–5.59, P = 0.01, Fig. 3B), but this pattern was not observed in patients treated with adjuvant chemotherapy (Fig. 3C, D, P > 0.05). Importantly, patients in the low-risk group were found to be benefit from chemotherapy (Fig. 3E, F, P < 0.01), but this was not seen in the high-risk group (Fig. 3G, H, P > 0.05).

Fig. 3
figure3

Kaplan–Meier plots for validations of FRGS in non-chemotherapy and chemotherapy cohorts. A, B The DFS of patients with different CAF risk group in stage II/III patients without adjuvant chemotherapy: GSE39582 cohort (A) and TCGA cohort (B). C, D The DFS of patients with different CAF risk group in stage II/III patients with adjuvant chemotherapy: GSE39582 cohort (C) and TCGA cohort (D). E, F The DFS of patients with or without adjuvant chemotherapy in stage II/III patients with low CAF risk: GSE39582 cohort (E) and TCGA cohort (F). G, H The DFS of patients with or without adjuvant chemotherapy in stage II/III patients with high CAF risk: GSE39582 cohort (G) and TCGA cohort (H)

Functional assessment of the FRGS

To explore the possible mechanism underlying the FRGS in high-risk patients, gene ontology (GO) enrichment analysis was performed to identify the biological processes of the differentially expressed genes. Interestingly, the most enriched biological processes other than the cell cycle were found to be associated with the immune response (Fig. 4A). To further evaluate the role of the FRGS, GSEA was performed in the TCGA CRC cohort, and the results showed that the CAF-induced microenvironment was significantly related to the inflammatory response, TNF-α, IFN-α, IFN-γ, IL-6 and IL-2 (Fig. 4B), which are related to the immune response network.

Fig. 4
figure4

Functional annotation of FRGS. A Enrichment analysis of the differentially expressed genes between different groups in GO. B GSEA analysis showed inflammatory response, TNF-α, IFN-α, IFN-γ, IL-6 and IL-2 were depressed in high-risk patients. C Infiltrations of immune cells are assessed based on TCGA cohort. D CD4+ T cells and M1 macrophages were depressed in the high fibroblastic tumor

The ESTIMATE algorithm was further used to investigate the immune response in patients in different groups (Fig. 4C). Although no difference was found in the stromal score, immune score or ESTIMATE score between the two groups (data not shown), lower percentages of CD4+ T cells and M1 macrophages were detected in the FRGS-identified high-risk group (Fig. 4D).

Assessment of the ability of the FRGS to predict the response to PD-1 treatment

To evaluate whether the FRGS could be used to predict the response to immune checkpoint inhibitors, a cohort of advanced clear cell renal cell carcinoma (ccRCC) patients was used (Braun et al. 2020). With a cut-off value of − 0.026 (identified in the ROC curve analysis), patients were divided into a high-risk group and a low-risk group. The PFS and OS in the high-risk group were significantly shorter than those in the low-risk group (PFS: HR = 1.3, 95% CI 1.02–1.65, P = 0.0318, OS: HR = 1.39, 95% CI 1.07–1.81, P = 0.0136, Fig. 5A, B). In terms of immunotherapy benefit, patients in the low-risk group were found to have improved PFS (P < 0.0001) and OS (P < 0.0001) compared to those in the high-risk group. After stratification according to the drug response, clinical benefit (CB) and intermediate clinical benefit (ICB) patients in the FRGS-identified low-risk group had a longer DFS time (Fig. 5C), while no clinical benefit (NCB) patients in the FRGS-identified low-risk group had a longer OS time (Fig. 5D).

Fig. 5
figure5

Kaplan–Meier plots for validations of FRGS in a cohort of ccRCC patients who received immunotherapy. A, B The outcome of patients with different CAF risk group in ccRCC patients with immunotherapy: PFS (A) and OS (B). C, D The outcome of patients with different drug response in ccRCC patients with different CAF risk group: PFS (A) and OS (B)

Discussion

Based on our previous genome-wide analysis (Zou et al. 2019), we developed a fibroblast-related gene signature (FRGS) for the prediction of outcome in patients who have undergone resection for stage II/III CRC. A discovery cohort of stage II/III patients who did not receive adjuvant chemotherapy was first used to identify the genes, and the 11 fibroblast-related gene signature was then evaluated in multiple cohorts of stage II/III patients, including one training cohort and two independent validation cohorts. In these cohorts, patients in the high- and low-risk groups differed in HR for DFS by approximately 3 to 5 times in the univariate and multivariate analyses. Moreover, patients in the high-risk group did not benefit from chemotherapy. Some studies have indicated that chemotherapy-treated patients have a less than 5% improvement in the 5-year survival rate (Glynne-Jones et al. 2017; Lee et al. 2018), and the treatment might not reduce the risk of tumour recurrence (Draht et al. 2014; Gunderson et al. 2010; Nazemalhosseini Mojarad et al. 2016). Thus, patients in the high-risk group might be suggested to abandon adjuvant chemotherapy, as there is no improvement in survival, indicating that this new signature might be a useful supplement to the current TNM staging system and modify risk stratification for CRC patients.

Most of the genes in the FRGS have been previously found to be correlated with antitumour immune responses and tumour progression. In addition, Gas6 has been implicated in the promotion of tumour cell proliferation, survival, migration, invasion, angiogenesis, and immune evasion, which it accomplishes via Axl signalling (Tanaka and Siemann 2020). IRF1 has been identified as an important transcription factor for M1 macrophage polarization (Chu et al. 2021), which is consistent with our result that the infiltration of M1 macrophages was lower in the high-risk group (which had lower expression of IRF1). BCL2L1 (Bcl-xL) was found to be a driver in colorectal tumorigenesis and cancer progression (Scherr et al. 2016), which is in agreement with its protumour effect in the FRGS.

An increased proportion of fibroblasts in the tumour microenvironment has been reported to be correlated with poor prognosis in various cancer patients (Nishina et al. 2021; Song et al. 2021; Zhou et al. 2018). In CRC, a recent study showed that the expression of genes enriched in IL-11+ fibroblasts was correlated with reduced recurrence-free survival in patients, as IL-11+ fibroblasts activated both tumour cells and fibroblasts via secretion of IL-11; this phenomenon constituted a feed-forward loop between tumour cells and fibroblasts (Nishina et al. 2021). Although increasing laboratory evidence has revealed that CAFs are associated with tumorigenesis and chemoresistance, there is no tool that can distinguish high/low fibroblast-related risk and predict the prognostic response to chemotherapy in CRC, which hinders the use of tools related to fibroblasts in the clinic. In this study, we selected various FRGs to create a fibroblast-related gene signature (FRGS) for CRC patients. Our results indicated that the FRGS could stratify stage II/III CRC patients into different 5-year DFS rate groups, and the FRGS score also showed a better C-index than the prognostic tool Oncotype DX (Additional file 3: Table S1). These results indicated that the FRGS might be an effective prognostic tool because it can mirror the fibroblast status in CRC tumours.

To develop a novel therapeutic strategy against patients with high risk according to fibroblast-related parameters, it is important to understand the mechanism underlying fibroblast-related poor prognosis. To date, a mechanistic analysis of genome-wide data to determine how the interaction between CAFs and the tumour microenvironment induces a poor chemotherapy response is lacking. Many therapeutic methods may have effects on both tumour cells and CAFs. However, the underlying mechanisms by which chemotherapy affects the population and function of CAFs and how CAFs affect the treatment response remain unclear. Recent studies have indicated that chemotherapy can regulate the activation and function of CAFs. In addition, the cytokines and chemokines secreted by CAFs can also lead to chemoresistance by leading to degradation of extracellular matrix (ECM) and remodelling of blood vessels. Additionally, while tumours are exposed to a chemotherapy-induced toxic environment, the cross-talk between CAFs and tumour cells can also induce chemoresistance (Tang et al. 2018). In this study, functional analysis of the FRGS identified significantly lower immune cell infiltrations in the high-risk group. CAFs and immune cells act as major components and the most important cells in the TME. Recent studies have shown that CAFs are a major source of the immunosuppressive activity in the TME that may dictate the curative effect of immunotherapies (Barrett and Pure 2020). An early study indicated that CAFs contributed to the immunosuppressive phenotype based on a correlation between stromal markers and immunosuppressive cell types such as tumour-associated macrophages (TAMs) and myeloid-derived suppressor cells (MDSCs) (Zhou et al. 2018). CAFs can regulate TAMs by inducing LIF-mediated paracrine signalling, eliciting CCL2 expression for myeloid cell attraction to the TME, and silencing CXCL9 expression to inhibit the recruitment of cytotoxic T-cells (Pascual-Garcia et al. 2019). Most studies have shown the synergistic effect of CAFs and M2 macrophages on tumour progression, CAFs could also inhibit the proinflammatory features of M1 macrophages by reducing proinflammatory cytokine levels, nitric oxide production, migration, and the expression of M1 surface markers (Berzaghi et al. 2019). Patient-derived CAFs can induce the adhesion and enrichment of monocytes by upregulating the expression of VCAM1 and the secretion of IL-8 and promote the recruitment of M2 macrophages into the tumour microenvironment, which together inhibit the function of NK cells (Zhang et al. 2019). Additionally, CAFs also promoted the expression of immune checkpoints, such as TIM-3, PD-1, CTLA-4 and LAG-3, in CD4+ and CD8+ T-cells through the secretion of prostaglandin E2 (PGE2), contributing to a diminished immune function (Gorchs et al. 2019). In line with these data, we also found fewer CD4+ T cells and M1 macrophages in tumours with high levels of fibroblasts. These data suggest that the impact of fibroblasts on the tumour microenvironment is probably associated with immune suppression, which further enhances tumour progression. A microenvironment with high levels of fibroblast-related features could provide tumour cells with resistance to immunotherapy (Ai et al. 2015; Scharping et al. 2017), which has emerged as an effective therapy for CRC, and the FRGS may have the potential to predict the outcome of immunotherapy. Due to the lack of a CRC cohort that received immunotherapy, the FRGS was tested in a cohort of ccRCC patients and was found to have value in predicting the benefit of immunotherapy as well as clinical outcome. Further studies in CRC cohorts will be required to verify this hypothesis.

The limitations of this study are its retrospective design and the sample bias due to intratumour genetic heterogeneity (Mimori et al. 2018). Even though as many available datasets as possible were used to validate the FRGS, challenges might still exist owing to the lack of testing in a prospective cohort. Furthermore, the gene expression signature might be largely affected by sampling bias, and not all the batch effects could be addressed, although constant ordering and exclusion of FRGs were used to reduce the batch effects between different studies. Thus, the FRGS should be validated in a prospective cohort before its clinical application.

Conclusion

In conclusion, the prognostic FRGS is a novel system that can be used to evaluate the DFS for stage II/III CRC patients and to stratify patients based on potential benefit from the adjuvant chemotherapy. CAFs in tumour microenvironment may promote chemoresistance in CRC patients by inhibiting the in situ antitumour immune response.

Availability of data and materials

The datasets generated and analyzed during the current study are available in the GSE39582 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE39582), TCGA (https://www.cancer.gov/about-nci/organization/ccg/research/structural-genomics/tcga), GSE14333 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE14333), GSE17536 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE17536), GSE37892 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE37892) and GSE33113 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE33113).

Abbreviations

CAFs:

Cancer-associated fibroblasts

CB:

Clinical benefit

ccRCC:

Clear cell renal cell carcinoma

CRC:

Colorectal cancer

DFS:

Disease-free survival

ECM:

Extracellular matrix

ESTIMATE:

Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data

FRGs:

Fibroblast-related genes

FRGS:

Fibroblast-related gene signature

GEO:

Gene Expression Omnibus

GO:

Gene ontology

GSEA:

Gene set enrichment analysis

ICB:

Intermediate clinical benefit

IQR:

Interquartile ranges

LASSO:

Least absolute shrinkage and selection operator

MDSCs:

Myeloid-derived suppressor cells

NCB:

No clinical benefit

OS:

Overall survival

PFS:

Progression-free survival

PGE2:

Prostaglandin E2

ROC:

Receiver operating characteristic

SD:

Standard deviations

TAMs:

Tumour-associated macrophages

TCGA:

The Cancer Genome Atlas

TME:

Tumour microenvironment

TPM:

Transcripts per million

References

  1. Ai M, Budhani P, Sheng J, Balasubramanyam S, Bartkowiak T, Jaiswal AR, Ager CR, Haria DD, Curran MA. Tumor hypoxia drives immune suppression and immunotherapy resistance. J Immunother Cancer. 2015;3:P392.

    Article  Google Scholar 

  2. Andre T, Boni C, Navarro M, Tabernero J, Hickish T, Topham C, Bonetti A, Clingan P, Bridgewater J, Rivera F, et al. Improved overall survival with oxaliplatin, fluorouracil, and leucovorin as adjuvant treatment in stage II or III colon cancer in the MOSAIC trial. J Clin Oncol. 2009;27:3109–16.

    CAS  Article  Google Scholar 

  3. Andre T, de Gramont A, Vernerey D, Chibaudel B, Bonnetain F, Tijeras-Raballand A, Scriva A, Hickish T, Tabernero J, Van Laethem JL, et al. Adjuvant fluorouracil, leucovorin, and oxaliplatin in stage II to III colon cancer: updated 10-year survival and outcomes according to BRAF mutation and mismatch repair status of the MOSAIC study. J Clin Oncol. 2015;33:4176–87.

    CAS  Article  Google Scholar 

  4. Barrett RL, Pure E. Cancer-associated fibroblasts and their influence on tumor immunity and immunotherapy. Elife. 2020;9:e57243.

    CAS  Article  Google Scholar 

  5. Berzaghi R, Ahktar MA, Islam A, Pedersen BD, Hellevik T, Martinez-Zubiaurre I. Fibroblast-mediated immunoregulation of macrophage function is maintained after irradiation. Cancers (basel). 2019;11:689.

    CAS  Article  Google Scholar 

  6. Braun DA, Hou Y, Bakouny Z, Ficial M, Sant’ Angelo M, Forman J, Ross-Macdonald P, Berger AC, Jegede OA, Elagina L, et al. Interplay of somatic alterations and immune infiltration modulates response to PD-1 blockade in advanced clear cell renal cell carcinoma. Nat Med. 2020;26:909–18.

    CAS  Article  Google Scholar 

  7. Chu YB, Li J, Jia P, Cui J, Zhang R, Kang X, Lv M, Zhang S. Irf1- and Egr1-activated transcription plays a key role in macrophage polarization: a multiomics sequencing study with partial validation. Int Immunopharmacol. 2021;99:108072.

    CAS  Article  Google Scholar 

  8. Draht MX, Smits KM, Tournier B, Jooste V, Chapusot C, Carvalho B, Cleven AH, Derks S, Wouters KA, Belt EJ, et al. Promoter CpG island methylation of RET predicts poor prognosis in stage II colorectal cancer patients. Mol Oncol. 2014;8:679–88.

    CAS  Article  Google Scholar 

  9. Glynne-Jones R, Wyrwicz L, Tiret E, Brown G, Rodel C, Cervantes A, Arnold D, Committee EG. Rectal cancer: ESMO clinical practice guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2017;28:vi22–40.

    Article  Google Scholar 

  10. Gorchs L, Fernandez Moro C, Bankhead P, Kern KP, Sadeak I, Meng Q, Rangelova E, Kaipe H. Human pancreatic carcinoma-associated fibroblasts promote expression of co-inhibitory markers on CD4(+) and CD8(+) T-cells. Front Immunol. 2019;10:847.

    CAS  Article  Google Scholar 

  11. Gunderson LL, Jessup JM, Sargent DJ, Greene FL, Stewart AK. Revised TN categorization for colon cancer based on national survival outcomes data. J Clin Oncol. 2010;28:264–71.

    Article  Google Scholar 

  12. Kalluri R. The biology and function of fibroblasts in cancer. Nat Rev Cancer. 2016;16:582–98.

    CAS  Article  Google Scholar 

  13. Kobayashi H, Enomoto A, Woods SL, Burt AD, Takahashi M, Worthley DL. Cancer-associated fibroblasts in gastrointestinal cancer. Nat Rev Gastroenterol Hepatol. 2019;16:282–95.

    Article  Google Scholar 

  14. Lee MM, MacKinlay A, Semira C, Schieber C, Jimeno Yepes AJ, Lee B, Wong R, Hettiarachchige CKH, Gunn N, Tie J, et al. Stage-based variation in the effect of primary tumor side on all stages of colorectal cancer recurrence and survival. Clin Colorectal Cancer. 2018;17:e569–77.

    Article  Google Scholar 

  15. Mimori K, Saito T, Niida A, Miyano S. Cancer evolution and heterogeneity. Ann Gastroenterol Surg. 2018;2:332–8.

    Article  Google Scholar 

  16. Nazemalhosseini Mojarad E, Kashfi SM, Mirtalebi H, Taleghani MY, Azimzadeh P, Savabkar S, Pourhoseingholi MA, Jalaeikhoo H, Asadzadeh Aghdaei H, Kuppen PJ, et al. Low level of microsatellite instability correlates with poor clinical prognosis in stage II colorectal cancer patients. J Oncol. 2016;2016:2196703.

    PubMed  Google Scholar 

  17. Nishina T, Deguchi Y, Ohshima D, Takeda W, Ohtsuka M, Shichino S, Ueha S, Yamazaki S, Kawauchi M, Nakamura E, et al. Interleukin-11-expressing fibroblasts have a unique gene signature correlated with poor prognosis of colorectal cancer. Nat Commun. 2021;12:2281.

    CAS  Article  Google Scholar 

  18. Pascual-Garcia M, Bonfill-Teixidor E, Planas-Rigol E, Rubio-Perez C, Iurlaro R, Arias A, Cuartas I, Sala-Hojman A, Escudero L, Martinez-Ricarte F, et al. LIF regulates CXCL9 in tumor-associated macrophages and prevents CD8(+) T cell tumor-infiltration impairing anti-PD1 therapy. Nat Commun. 2019;10:2416.

    Article  Google Scholar 

  19. Quasar Collaborative Group. Comparison of fluorouracil with additional levamisole, higher-dose folinic acid, or both, as adjuvant chemotherapy for colorectal cancer: a randomised trial. Lancet. 2000;355:1588–96.

    Article  Google Scholar 

  20. Quasar Collaborative Group, Gray R, Barnwell J, McConkey C, Hills RK, Williams NS, Kerr DJ. Adjuvant chemotherapy versus observation in patients with colorectal cancer: a randomised study. Lancet. 2007;370:2020–9.

    Article  Google Scholar 

  21. Scharping NE, Menk AV, Whetstone RD, Zeng X, Delgoffe GM. Efficacy of PD-1 blockade is potentiated by metformin-induced reduction of tumor hypoxia. Cancer Immunol Res. 2017;5:9–16.

    CAS  Article  Google Scholar 

  22. Scherr AL, Gdynia G, Salou M, Radhakrishnan P, Duglova K, Heller A, Keim S, Kautz N, Jassowicz A, Elssner C, et al. Bcl-xL is an oncogenic driver in colorectal cancer. Cell Death Dis. 2016;7:e2342.

    CAS  Article  Google Scholar 

  23. Song M, He J, Pan QZ, Yang J, Zhao J, Zhang YJ, Huang Y, Tang Y, Wang Q, He J, et al. Cancer-associated fibroblast-mediated cellular crosstalk supports hepatocellular carcinoma progression. Hepatology. 2021;73:1717–35.

    CAS  Article  Google Scholar 

  24. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71:209–49.

    Article  Google Scholar 

  25. Tanaka M, Siemann DW. Gas6/Axl signaling pathway in the tumor immune microenvironment. Cancers (basel). 2020;12:1850.

    CAS  Article  Google Scholar 

  26. Tang YA, Chen YF, Bao Y, Mahara S, Yatim S, Oguz G, Lee PL, Feng M, Cai Y, Tan EY, et al. Hypoxic tumor microenvironment activates GLI2 via HIF-1alpha and TGF-beta2 to promote chemoresistance in colorectal cancer. Proc Natl Acad Sci USA. 2018;115:E5990–9.

    CAS  Article  Google Scholar 

  27. Tauriello DVF, Palomo-Ponce S, Stork D, Berenguer-Llergo A, Badia-Ramentol J, Iglesias M, Sevillano M, Ibiza S, Canellas A, Hernando-Momblona X, et al. TGFbeta drives immune evasion in genetically reconstituted colon cancer metastasis. Nature. 2018;554:538–43.

    CAS  Article  Google Scholar 

  28. Wang X, Terfve C, Rose JC, Markowetz F. HTSanalyzeR: an R/Bioconductor package for integrated network analysis of high-throughput screens. Bioinformatics. 2011;27:879–80.

    CAS  Article  Google Scholar 

  29. Zhang R, Qi F, Zhao F, Li G, Shao S, Zhang X, Yuan L, Feng Y. Cancer-associated fibroblasts enhance tumor-associated macrophages enrichment and suppress NK cells function in colorectal cancer. Cell Death Dis. 2019;10:273.

    Article  Google Scholar 

  30. Zhou J, Wang XH, Zhao YX, Chen C, Xu XY, Sun Q, Wu HY, Chen M, Sang JF, Su L, et al. Cancer-associated fibroblasts correlate with tumor-associated macrophages infiltration and lymphatic metastasis in triple negative breast cancer patients. J Cancer. 2018;9:4635–41.

    CAS  Article  Google Scholar 

  31. Zou YF, Rong YM, Tan YX, Xiao J, Yu ZL, Chen YF, Ke J, Li CH, Chen X, Wu XJ, et al. A signature of hypoxia-related factors reveals functional dysregulation and robustly predicts clinical outcomes in stage I/II colorectal cancer patients. Cancer Cell Int. 2019;19:243.

    Article  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Yidu Cloud Technology Co., Ltd for the assistance in data processing and other members of The Sixth Affiliated Hospital, Sun Yat-sen University for the support and useful discussions.

Funding

This work was supported by National Key R&D Program of China (No. 2017YFC1308800), National Key Clinical Discipline, National Natural Science Foundation of China (No. 82003197), Science and Technology Planning Project of Guangdong Province, China (No. 2021A0505030028), Science and Technology Planning Project of Guangzhou City (No. 202102020186).

Author information

Affiliations

Authors

Contributions

YC, XW and FG contributed to study concept and design, acquisition, analysis, interpretation of data and drafting of the manuscript. YC, ZY, ML, ZC and YZ contributed to data collections and manuscript review. YC and ZY contributed to study concept and design, analysis and interpretation of data and critical revision of the manuscript for important intellectual content. PL, XW and FG supervised the study. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yu-feng Chen, Xiao-jian Wu or Feng Gao.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Board (IRB) of The Sixth Affiliated Hospital of Sun Yat-sen University. This is a retrospective trial from public datasets with demonstrated minimal risk and we petition BPG for waiver of ethics consent.

Consent for publication

We have obtained consents to publish this paper from all the participants of this study.

Competing interests

The authors declared no financial conflict of interest. The authors declare that there are no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Figure S1. Risk score calculation model of FRGS

Additional file 2:

Data S1. The risk of each patient in the six cohorts.

Additional file 3:

Table S1. Comparation of FRGS and Oncotype in each dataset. 

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chen, Yf., Yu, Zl., Lv, My. et al. Cancer-associated fibroblasts impact the clinical outcome and treatment response in colorectal cancer via immune system modulation: a comprehensive genome-wide analysis. Mol Med 27, 139 (2021). https://doi.org/10.1186/s10020-021-00402-3

Download citation

Keywords

  • Fibroblast-related gene signature
  • Colorectal cancer
  • Prognosis
  • Chemotherapy