Identification of subtypes correlated with tumor immunity and immunotherapy in cutaneous melanoma

Because immune checkpoint inhibitors (ICIs) are effective for a subset of melanoma patients, identification of melanoma subtypes responsive to ICIs is crucial. We performed clustering analyses to identify immune subtypes of melanoma based on the enrichment levels of 28 immune cells using transcriptome datasets for six melanoma cohorts, including four cohorts not treated with ICIs and two cohorts treated with ICIs. We identified three immune subtypes (Im-H, Im-M, and Im-L), reproducible in these cohorts. Im-H displayed strong immune signatures, low stemness and proliferation potential, genomic stability, high immunotherapy response rate, and favorable prognosis. Im-L showed weak immune signatures, high stemness and proliferation potential, genomic instability, low immunotherapy response rate, and unfavorable prognosis. The pathways highly enriched in Im-H included immune, MAPK, apoptosis, calcium, VEGF, cell adhesion molecules, focal adhesion, gap junction, and PPAR. The pathways highly enriched in Im-L included Hippo, cell cycle, and ErbB. Copy number alterations correlated inversely with immune signatures in melanoma, while tumor mutation burden showed no significant correlation. The molecular features correlated with favorable immunotherapy response included immune-promoting signatures and pathways of PPAR, MAPK, VEGF, calcium, and glycolysis/gluconeogenesis. Our data recapture the immunological heterogeneity in melanoma and provide clinical implications for the immunotherapy of melanoma.