# Annotation-free prediction of immunotherapy response in melanoma using single-cell transcriptomic data

**Authors:** Da Eun Oh, Gaeun Kee, Ji-Hye Oh, Wonkyung Kim, Young Gwang Kang, Chae Won Park, Tae Joon Jun, Chang Ohk Sung

PMC · DOI: 10.1371/journal.pone.0343633 · PLOS One · 2026-02-27

## TL;DR

This study uses AI on single-cell RNA data to predict which melanoma patients will respond to immunotherapy, without needing cell type labels.

## Contribution

The novel contribution is an annotation-free AI model for predicting immunotherapy response using single-cell transcriptomic data.

## Key findings

- An AI model achieved 0.87 AUC in predicting ICI response using unannotated scRNA-seq data.
- CCR7 was linked to better ICI response and survival, while MTRNR2L2 was associated with poor outcomes.
- CCR7 was mainly expressed in B and memory T cells, and was tied to immune subtype and NF-κB activation.

## Abstract

Immune checkpoint inhibitors (ICIs) have transformed the advanced melanoma treatment landscape; however, a subset of patients achieve durable responses. Current biomarkers, such as PD-L1 expression and tumor mutational burden, offer limited predictive power due to the profound heterogeneity of melanoma. Accordingly, we developed an artificial intelligence (AI)-based model to predict ICI responsiveness using single-cell RNA sequencing (scRNA-seq) data without requiring cell type annotation. scRNA-seq data profiled using Smart-seq2 platform were downloaded from a public repository (GEO: GSE120575). From these data, we analyzed 16,290 tumor-infiltrating cells from melanoma scRNA-seq dataset. Various AI-based models, including Extreme Gradient Boosting, Random Forest, Logistic Regression, Support Vector Machine, Feedforward Neural Network, and Convolutional Neural Network were constructed, with the best-performing model achieving an area under the curve of 0.87. This AI-driven approach identified 29 key predictive biomarkers, including CCR7 and MTRNR2L2. Validation using three independent bulk RNA-seq datasets (cBioPortal: DFCI melanoma; ENA: PRJEB23709; GEO: GSE91061) suggested that CCR7 was associated with favorable ICI response and improved survival, whereas MTRNR2L2 showed a tendency toward enrichment in non-responders and poorer outcomes. Cell-type-specific expression analysis revealed that CCR7 was primarily expressed in B cells and memory T cells from responders, whereas MTRNR2L2 was elevated in exhausted and cytotoxic T cells in non-responders. CCR7-positive B cells exhibited activation of the NF-κB pathway and demonstrated prognostic significance independent of the melanoma primary site or histologic subtype. However, among the three molecular subtypes, including immune, keratin, and microphthalmia-associated transcription factor (MITF)-low, CCR7 expression was significantly associated with the immune subtype. Additionally, pathway-level deep learning models reinforced these findings, highlighting immune activation in responders and cell cycle-related signals in non-responders. Our study demonstrates that predictive modeling based on unannotated scRNA-seq data enables clinically relevant biomarker identification, offering a robust approach for patients with stratifying melanoma and guiding personalized immunotherapy.

## Linked entities

- **Genes:** CCR7 (C-C motif chemokine receptor 7) [NCBI Gene 1236], MTRNR2L2 (MT-RNR2 like 2 (pseudogene)) [NCBI Gene 100462981]
- **Diseases:** melanoma (MONDO:0005105)

## Full-text entities

- **Genes:** CD3D (CD3 delta subunit of T-cell receptor complex) [NCBI Gene 915] {aka CD3-DELTA, CD3DELTA, IMD19, T3D}, HLA-A (major histocompatibility complex, class I, A) [NCBI Gene 3105] {aka HLAA}, IFNG (interferon gamma) [NCBI Gene 3458] {aka IFG, IFI, IMD69}, IFITM3 (interferon induced transmembrane protein 3) [NCBI Gene 10410] {aka 1-8U, DSPA2b, IP15}, CD274 (CD274 molecule) [NCBI Gene 29126] {aka ADMIO5, B7-H, B7H1, PD-L1, PDCD1L1, PDCD1LG1}, EGFR (epidermal growth factor receptor) [NCBI Gene 1956] {aka ERBB, ERBB1, ERRP, HER1, NISBD2, NNCIS}, RPL9P8 (ribosomal protein L9 pseudogene 8) [NCBI Gene 254948] {aka RPL9P9, RPL9_14_1458, RPL9_15_1459}, CGAS (cyclic GMP-AMP synthase) [NCBI Gene 115004] {aka C6orf150, D4, MB21D1, h-cGAS}, MTRNR2L2 (MT-RNR2 like 2 (pseudogene)) [NCBI Gene 100462981] {aka HN2}, CD8A (CD8 subunit alpha) [NCBI Gene 925] {aka CD8, CD8alpha, IMD116, Leu2, p32}, TIGIT (T cell immunoreceptor with Ig and ITIM domains) [NCBI Gene 201633] {aka VSIG9, VSTM3, WUCAM}, LAG3 (lymphocyte activating 3) [NCBI Gene 3902] {aka CD223}, CTLA4 (cytotoxic T-lymphocyte associated protein 4) [NCBI Gene 1493] {aka ALPS5, CD, CD152, CELIAC3, CTLA-4, GRD4}, TET3 (tet methylcytosine dioxygenase 3) [NCBI Gene 200424] {aka BEFAHRS, hCG_40738}, CXCL10 (C-X-C motif chemokine ligand 10) [NCBI Gene 3627] {aka C7, IFI10, INP10, IP-10, SCYB10, crg-2}, CCL5 (C-C motif chemokine ligand 5) [NCBI Gene 6352] {aka D17S136E, RANTES, SCYA5, SIS-delta, SISd, TCP228}, STING1 (stimulator of interferon response cGAMP interactor 1) [NCBI Gene 340061] {aka ERIS, MITA, MPYS, NET23, SAVI, STING}, BRAF (B-Raf proto-oncogene, serine/threonine kinase) [NCBI Gene 673] {aka B-RAF1, B-raf, BRAF-1, BRAF1, NS7, RAFB1}, MTRNR2L8 (MT-RNR2 like 8 (pseudogene)) [NCBI Gene 100463486] {aka HN8}, HSPA1B (heat shock protein family A (Hsp70) member 1B) [NCBI Gene 3304] {aka HSP70-1, HSP70-1B, HSP70-2, HSP70.1, HSP70.2, HSP72}, CXCL9 (C-X-C motif chemokine ligand 9) [NCBI Gene 4283] {aka CMK, Humig, MIG, SCYB9, crg-10}, IFI6 (interferon alpha inducible protein 6) [NCBI Gene 2537] {aka 6-16, FAM14C, G1P3, IFI-6-16, IFI616}, MAP2K7 (mitogen-activated protein kinase kinase 7) [NCBI Gene 5609] {aka JNKK2, MAPKK7, MEK, MEK 7, MKK7, PRKMK7}, CD3E (CD3 epsilon subunit of T-cell receptor complex) [NCBI Gene 916] {aka CD3epsilon, IMD18, T3E, TCRE}, NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790] {aka CVID12, EBP-1, KBF1, NF-kB, NF-kB1, NF-kappa-B1}, HAVCR2 (hepatitis A virus cellular receptor 2) [NCBI Gene 84868] {aka CD366, HAVcr-2, KIM-3, SPTCL, TIM3, TIMD-3}, PDCD1 (programmed cell death 1) [NCBI Gene 5133] {aka ADMIO4, AIMTBS, CD279, PD-1, PD1, SLEB2}, HSPA1A (heat shock protein family A (Hsp70) member 1A) [NCBI Gene 3303] {aka HEL-S-103, HSP70, HSP70-1, HSP70-1A, HSP70-2, HSP70.1}, MITF (melanocyte inducing transcription factor) [NCBI Gene 4286] {aka CMM8, COMMAD, MI, MITF-A, WS2, WS2A}, PTPRC (protein tyrosine phosphatase receptor type C) [NCBI Gene 5788] {aka B220, CD45, CD45R, GP180, IMD105, L-CA}, CCR7 (C-C motif chemokine receptor 7) [NCBI Gene 1236] {aka BLR2, CC-CKR-7, CCR-7, CD197, CDw197, CMKBR7}, MT2A (metallothionein 2A) [NCBI Gene 4502] {aka MT-2, MT-II, MT2}
- **Diseases:** PD (MESH:D018450), cutaneous melanoma (MESH:C562393), autoimmune thyroid disease (MESH:D013967), NM (MESH:C536816), metastasis (MESH:D009362), hypoxia (MESH:D000860), oncogenesis (MESH:D063646), PR (MESH:D004828), Cancer (MESH:D009369), SD (MESH:D060050), inflammatory (MESH:D007249), disease (MESH:D004194), Malignant melanoma (MESH:D008545)
- **Chemicals:** formalin (MESH:D005557), steroid (MESH:D013256), serine (MESH:D012694), Glyphosate (MESH:C010974), glycine (MESH:D005998), pentose-phosphate (MESH:D010428), paraffin (MESH:D010232)
- **Species:** Homo sapiens (human, species) [taxon 9606]

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12948085/full.md

## Figures

7 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12948085/full.md

## References

67 references — full list in the complete paper: https://tomesphere.com/paper/PMC12948085/full.md

---
Source: https://tomesphere.com/paper/PMC12948085