Identification of immune-related genes with prognostic significance in the microenvironment of cutaneous melanoma

Cutaneous melanoma is one of the most aggressive cancers characterized by increasing incidence and mortality. In recent years, the emergence of immunotherapy has greatly raised the survival rate of patients suffering from cutaneous melanoma, yet some sufferers remain to have poor outcomes after treatment mainly due to the tumor microenvironment (TME). In this study, cutaneous melanoma-associated TME was systematically analyzed using the ESTIMATE algorithm based on the gene transcriptome data obtained from the TCGA database. Totally, 471 patients were included and 553 TME-related genes were screened. Afterwards, a 3-gene signature-based model (CLEC4A, GBP4, KIR2DL4) was constructed via univariate Cox, LASSO, and multivariate Cox regression analyses. To validate the validity of this model, ROC analysis was conducted, and the model was further validated to be an independent prognostic biomarker through univariate and multivariate regression analyses. Finally, the three genes in the model were studied by GSEA and GSVA for their biological significance. We found that the three genes could promote cancer immune response predominantly through affecting immune-related pathways such as antigen processing and presentation, and they may help tumor cells in escaping from surveillance of the immune system when their expression levels were decreased. Additionally, we as well discovered that the expression of the three genes was significantly and positively correlated with the infiltration of related immune cells, but negatively associated with tumor purity. Overall, this study comprehensively analyzed the TME of cutaneous melanoma, identified related biomarkers, and discovered their association with immune system.