DNA damage repair molecular subtype derived immune signature applicable for the prognosis and immunotherapy response prediction in colon cancer

BACKGROUND: The DNA damage repair (DDR) pathway is one of the pathways of tumor pathogenesis, but its relationship with the immunophenotype has not been clarified in colon cancer (CC). METHODS: We identified the differentially expressed immune-related genes (DEIRGs) between two DDR molecular subtypes, namely, C1 and C2, and used univariate Cox analysis and least absolute shrinkage and selection operator (LASSO) penalized Cox regression analysis to construct the risk score in the training cohort [n=1,009, a combination of The Cancer Genome Atlas (TCGA) and GSE39582]. Regarding the median risk score as the unified cutoff to classify the patients into high- and low-risk groups. Two independent cohorts (GSE17538, n=232; GSE38832, n=122) were used for external validation of the prognostic value of the risk score. The IMvigor210 cohort (n=348) was used to test the predictive value of the risk score for immunotherapy response. Gene set variation analysis (GSVA) and gene set enrichment analysis (GSEA) were performed to discover the underlying mechanism. Immune cell infiltration was quantified by the single sample gene set enrichment analysis (ssGSEA) algorithm. RESULTS: The high-risk group showed significantly reduced overall survival (OS), disease-specific survival (DSS), disease-free survival (DFS), progression-free survival (PFS), and relapse-free survival (RFS) compared to the low-risk group, and the two groups differed significantly in lymphatic invasion, American Joint Committee on Cancer (AJCC) TNM stage, preoperative carcinoembryonic antigen (CEA) level, etc. The enrichment levels of pathways related to colorectal cancer, epithelial-mesenchymal transition (EMT), angiogenesis, hypoxia, P53, TGF-β, KRAS signaling, etc., were upregulated in the high-risk group, but DDR-related pathways were defective in the high-risk group. The immunophenotypes of the high-risk group tended to be desert and excluded, and the risk score of patients who responded to immunotherapy was significantly lower than that of patients who did not respond to immunotherapy. The higher the infiltration levels of gamma delta T cells (γδ T cells), immature dendritic cells, and T follicular helper (Tfh) cells, the more significant adverse impact on the prognosis of CC patients was exhibited and an obviously positive correlation with the risk score was showed. CONCLUSIONS: An immune gene risk score associated with the DDR molecular subtype was built and verified herein; that is applicable to the prognosis and immunotherapy response prediction in CC.