A novel prognostic model for cutaneous melanoma based on an immune-related gene signature and clinical variables

Abundant evidence has indicated that the prognosis of cutaneous melanoma (CM) patients is highly complicated by the tumour immune microenvironment. We retrieved the clinical data and gene expression data of CM patients in The Cancer Genome Atlas (TCGA) database for modelling and validation analysis. Based on single-sample gene set enrichment analysis (ssGSEA) and consensus clustering analysis, CM patients were classified into three immune level groups, and the differences in the tumour immune microenvironment and clinical characteristics were evaluated. Seven immune-related CM prognostic molecules, including three mRNAs (SUCO, BTN3A1 and TBC1D2), three lncRNAs (HLA-DQB1-AS1, C9orf139 and C22orf34) and one miRNA (hsa-miR-17-5p), were screened by differential expression analysis, ceRNA network analysis, LASSO Cox regression analysis and univariate Cox regression analysis. Their biological functions were mainly concentrated in the phospholipid metabolic process, transcription regulator complex, protein serine/threonine kinase activity and MAPK signalling pathway. We established a novel prognostic model for CM integrating clinical variables and immune molecules that showed promising predictive performance demonstrated by receiver operating characteristic curves (AUC ≥ 0.74), providing a scientific basis for predicting the prognosis and improving the clinical outcomes of CM patients.