Abstract
Immunogenic cell death (ICD) is a type of regulated cell death that can activate adaptive immune response, and its ability to reshape the tumor microenvironment via multiple mechanisms may contribute to immunotherapy. The treatment options for patients with skin cutaneous melanoma (SKCM) vary based on BRAF V600E statuses. However, all standard treatments include immunotherapy. Therefore, it is critical to identify ICD-associated signatures that can help classify patients according to benefits from ICD immunotherapy. In this study, data on melanoma samples with BRAF V600E mutation (BRAF V600E-mutant melanoma) and melanoma samples with wild-type BRAF V600E alleles (BRAF V600E WT melanoma) were collected from The Cancer Genome Atlas (TCGA) database. The ICD-related (ICD-high and ICD-low) subgroups of patients with BRAF V600E WT melanoma were established via consensus clustering. The analyses of survival, differentially expressed genes (DEGs), functional annotation, and immune landscape were performed in these two subgroups. Results showed that ICD-high subgroup was correlated with a positive overall survival (OS) and active tumor immune landscape. A model comprising seven prognosis ICD-related gene biomarkers was developed. Survival analysis and receiver operating characteristic (ROC) curve evaluation in both cohorts with BRAF V600E WT and BRAF V600E-mutant melanoma showed an accurate prognostic estimation of ICD-related risk signature. There was a correlation between immune cell infiltration and immunotherapy response and risk score. Thus, the ICD risk signature was closely associated with the tumor's immune microenvironment. Our results may provide insights to further individualize and improve precision therapeutic decision-making in BRAF V600E-mutant and WT melanoma.