Identification of prognostic genes and construction of a novel gene signature in the skin melanoma based on the tumor microenvironment

Abstract Skin melanoma remains a highly prevalent and yet deadly form of cancer, with the exact degree of melanoma-associated mortality being strongly dependent upon the local tumor microenvironment. The exact composition of stromal and immune cells within this microenvironmental region has the potential to profoundly impact melanoma progression and prognosis. As such, the present study was designed with the goal of clarifying the predictive relevance of stromal and immune cell-related genes in melanoma patients through comprehensive bioinformatics analyses. We therefore analyzed melanoma sample gene expression within The Cancer Genome Atlas database and employed the ESTIMATE algorithm as a means of calculating both stromal and immune scores that were in turn used for identifying differentially expressed genes (DEGs). Subsequently, univariate analyses were used to detect DEGs associated with melanoma patient survival, and through additional functional enrichment analyses, we determined that these survival-related DEGs are largely related to inflammatory and immune responses. A prognostic signature comprised of 10 genes (IL15, CCL8, CLIC2, SAMD9L, TLR2, HLA.DQB1, IGHV1–18, RARRES3, GBP4, APOBEC3G) was generated. This 10-gene signature effectively separated melanoma patients into low- and high-risk groups based upon their survival. These low- and high-risk groups also exhibited distinct immune statuses and differing degrees of immune cell infiltration. In conclusion, our results offer novel insights into a number of microenvironment-associated genes that impact survival outcomes in melanoma patients, potentially highlighting these genes as viable therapeutic targets.