Identification and analysis of a CD8+ T cell-related prognostic signature for colorectal cancer based on bulk RNA sequencing and scRNA sequencing data: A STROBE-compliant retrospective study

Colorectal cancer (CRC) is one of the most common malignancies worldwide, leading to a large number of cancer-related mortalities. Aberrant CD8+ T cell infiltration plays a critical role in tumor progression and patient prognosis. This study aimed to identify a prognostic model for CRC based on CD8+ T cell-related genes. The infiltration levels of immune cells in CRC tissues were accessed by the ESTIMATE algorithm. Weighted gene co-expression network analysis (WGCNA) analysis was used to select CD8+ T cell-related genes. Prognostic genes were identified using Cox regression analysis and Kaplan–Meier curves. The least absolute shrinkage and selection operator (LASSO) algorithm was used to construct prognostic models. Gene set enrichment analysis (GSEA) was performed to annotate enriched gene sets. Single-cell RNA (scRNA) sequencing analysis was used to examine gene expression in different cell types. We found that the downregulated infiltration level of CD8+ T cells was an independent prognostic factor for CRC and selected a cluster of differentially expressed genes correlated with CD8+ T cell infiltration (CD8TDEGs). Subsequently, we identified 18 prognostic CD8TDEGs, according to which patients were reclassified into two clusters with distinct overall survival. Seven prognostic CD8TDEGs were selected to calculate the constructed prognostic model’s risk scores. Interestingly, although CRC tissues with higher risk scores had higher infiltration levels of CD8+ T cells, the level of immune checkpoint genes was also high. Moreover, the scRNA-sequencing analysis showed that the expression levels of CD8TDEGs in the prognostic model varied among different types of cells. This study constructed a novel prognostic model for CRC and provided a foundation for targeting CD8+ T cell infiltration to improve the survival of CRC patients.