CAF signature predicts the prognosis of colorectal cancer patients: A retrospective study based on bulk RNA sequencing and single-cell RNA sequencing data

The incidence rate and mortality rate of colorectal cancer (CRC) ranks third and second globally. Cancer-associated fibroblasts (CAFs) are the major constituent of the stromal cells in the tumor microenvironment (TME) and are closely associated with patients’ prognoses. Our study intended to establish a prognostic model for CRC using hallmark genes of CAFs. The expression values of genes and clinicopathological characteristics of patients were enrolled from the cancer genome atlas database as well as the gene expression omnibus database. The single-cell RNA sequencing data were collected and analyzed in the deeply integrated human single-cell omics database and cancer single-cell expression map databases. The ESTIMATE algorithm was applied to access the infiltration levels of immune and stromal cells. The prognostic genes were selected by the Cox regression analysis and the prognostic signature was constructed by the least absolute shrinkage and selection operator algorithm. gene set enrichment analysis was used to explore the enriched gene sets. In this study, based on bulk RNA sequencing and single-cell RNA sequencing data, and we found that more CAFs were infiltrated in the tumor microenvironment and consisted of 3 subtypes. Then we constructed a prognostic signature for CRC using hallmark genes of CAFs and proved that this signature exhibited high values to predict the overall survival of CRC patients in independent training and validating cohorts. Besides, function enrichment analysis revealed that our prognostic model was significantly associated with immune regulation. Further analysis showed that the infiltrated levels of tumor-suppressing immune cells and the expression of higher immune checkpoint genes in CRC tissues were higher in patients with high-risk scores. Furthermore, immunohistochemistry analysis exhibited that these genes in our prognostic signature were markedly upregulated in CRC tissues. We first constructed a signature based on CAFs hallmark genes to predict the survival of CRC patients and further revealed that the tumor-suppressing microenvironment and dysregulated immune checkpoint genes in CRC tissues were partly responsible for the poor prognosis of patients.