Development of an Immune-Related Gene Signature for Prognosis in Melanoma

Abstract

Melanoma remains a potentially deadly malignant tumor. The incidence of melanoma continues to rise. Immunotherapy has become a new treatment method and is widely used in a variety of tumors. Original melanoma data were downloaded from TCGA. ssGSEA was performed to classify them. GSVA software and the "hclust" package were used to analyze the data. The ESTIMATE algorithm screened DEGs. The edgeR package and Venn diagram identified valid immune-related genes. Univariate, LASSO and multivariate analyses were used to explore the hub genes. The "rms" package established the nomogram and calibrated the curve. Immune infiltration data were obtained from the TIMER database. Compared with that of samples in the high immune cell infiltration cluster, we found that the tumor purity of samples in the low immune cell infiltration cluster was higher. The immune score, ESTIMATE score and stromal score in the low immune cell infiltration cluster were lower. In the high immune cell infiltration cluster, the immune components were more abundant, while the tumor purity was lower. The expression levels of TIGIT, PDCD1, LAG3, HAVCR2, CTLA4 and the HLA family were also higher in the high immune cell infiltration cluster. Survival analysis showed that patients in the high immune cell infiltration cluster had shorter OS than patients in the low immune cell infiltration cluster. IGHV1-18, CXCL11, LTF, and HLA-DQB1 were identified as immune cell infiltration-related DEGs. The prognosis of melanoma was significantly negatively correlated with the infiltration of CD4+ T cells, CD8+ T cells, dendritic cells, neutrophils and macrophages. In this study, we identified immune-related melanoma core genes and relevant immune cell subtypes, which may be used in targeted therapy and immunotherapy of melanoma.