Abstract
Immune checkpoint inhibitors currently serve an important role in prolonging patients' overall survival. However, the prognostic signatures of immune checkpoint inhibitors in colorectal cancer (CRC) remain uncertain and more knowledge on the genetic characteristics of colorectal cancer is needed. Patients with CRC from The Cancer Genome Atlas were classified into high-immunity group and low-immunity group based on median scores from single-sample gene set enrichment analysis using the GSVA package. We explored immune status by immune scores, stromal scores and tumor purity scores in ESTIMATE package and surveyed the difference of immune cells distribution with CIBERSORT package. Eighteen genes were selected using the LASSO Cox regression method and a prognostic risk model was constructed. Compared with patients in the low-risk group, those in the high-risk group had a significantly shorter survival time. For assessment of the prognostic validity of the risk model, receiver operating characteristic curves with areas under the curve of 0.769, 0.774 and 0.771 for 1, 3 and 5 years respectively. Differences in molecular mechanisms between high- and low-risk groups were analyzed using the clusterProfiler package. Tumor Immune Dysfunction and Exclusion data were downloaded and analyzed. The top 5 enriched pathways in the high-risk group involved 'calcium signaling', 'dilated cardiomyopathy', 'extracellular matrix receptor interaction', 'hypertrophic cardiomyopathy' and 'neuroactive ligand receptor interaction'. HAMP was identified as a hub gene, which was highly expressed in tumor samples. The results of the present study indicate that the prognostic model based on both immune-related genes and HAMP has the potential to support personalized treatment.