Identification of Novel Molecular Therapeutic Targets and Their Potential Prognostic Biomarkers Based on Cytolytic Activity in Skin Cutaneous Melanoma

Skin cutaneous melanoma (SKCM) attracts attention worldwide for its extremely high malignancy. A novel term cytolytic activity (CYT) has been introduced as a potential immunotherapy biomarker associated with counter-regulatory immune responses and enhanced prognosis in tumors. In this study, we extracted all datasets of SKCM patients, namely, RNA sequencing data and clinical information from The Cancer Genome Atlas (TCGA) database and the Gene Expression Omnibus (GEO) database, conducted differential expression analysis to yield 864 differentially expressed genes (DEGs) characteristic of CYT and used non-negative matrix factorization (NMF) method to classify molecular subtypes of SKCM patients. Among all genes, 14 hub genes closely related to prognosis for SKCM were finally screen out. Based on these genes, we constructed a 14-gene prognostic risk model and its robustness and strong predictive performance were further validated. Subsequently, the underlying mechanisms in tumor pathogenesis and prognosis have been defined from a number of perspectives, namely, tumor mutation burden (TMB), copy number variation (CNV), tumor microenvironment (TME), infiltrating immune cells, gene set enrichment analysis (GSEA) and immune checkpoint inhibitors (ICIs). Furthermore, combined with GTEx database and HPA database, the expression of genes in the model was verified at the transcriptional level and protein level, and the relative importance of genes in the model was described by random forest algorithm. In addition, the model was used to predict the difference in sensitivity of SKCM patients to chemotherapy and immunotherapy. Finally, a nomogram was constructed to better aid clinical diagnosis.