# Integrative single‐cell analysis uncovers distinct tumour microenvironment ecotypes and immune evasion across skin cancers

## TL;DR

This study uses single-cell analysis to uncover immune changes in skin cancers, identifying new cell types and signaling patterns linked to tumor progression and immune evasion.

## Contribution

The study introduces a pan-skin cancer immune remodeling framework and identifies NARS2, NDUFC2, and SPP1 as potential therapeutic targets.

## Key findings

- A malignant NARS2+NDUFC2+ melanoma cell subpopulation is associated with worse survival and reduced MHC-I expression.
- TAMs originate from two sources and show distinct polarization states with different prognostic roles in SCC and melanoma.
- Immune ecotypes transition from T-cell-dominant to desert-like as disease progresses, linked to SPP1+ macrophage signaling.

## Abstract

Skin cancers, including basal cell carcinoma (BCC), squamous cell carcinoma (SCC), cutaneous melanoma (CM) and acral melanoma (AM), exhibit profound heterogeneity in clinical behaviour and therapeutic response. However, how tumour‐immune ecosystems are remodelled across skin cancer types and disease stages, and how these changes influence immune escape and treatment resistance, remain poorly understood.

Here, we integrate single‐cell transcriptomics data from 102 skin cancer samples (including adjacent normal skin, early‐stage and advanced‐stage tumours), with bulk RNA‐seq prognosis cohorts, immunofluorescence staining and in vitro assays to define clinically relevant immune remodelling patterns.

Our analyses identify a malignant NARS2+NDUFC2+
melanoma cell subpopulation, characterised by reduced MHC‐I expression, enriched in advanced‐stage tumours and associated with worse survival and immunotherapy response. CRISPR screening further showed that NARS2 and NDUFC2 are necessary for the proliferation of melanoma cells, highlighting these genes as potential therapeutic targets. Tumour‐associated macrophages (TAMs) originate from both FCN1
+ monocytes and FOLR2+
tissue‐resident macrophages, displaying two polarisation states with distinct prognostic associations. Specifically, pro‐inflammatory CXCL9+CXCL10+
TAMs are enriched in SCC, while tissue‐remodelling SPP1
+ TAMs are predominant in melanoma. Immunofluorescence staining confirmed that SPP1
+ macrophage accumulation correlates with advanced stage, metastasis and poor prognosis in the melanoma cohort. Immune ecotype analysis reveals a transition from ‘T‐cell‐dominant’ ecotypes to ‘desert’ ecotypes as disease advances in BCC, CM and AM. Cell‒cell communication analysis shows that ‘T‐cell‐dominant’ ecotypes have higher MHC‐I signalling pathways in tumour cells, whereas ‘Desert’ ecotypes have higher SPP1+
macrophage signalling, underlining the role of SPP1 on immune remodelling. Functional assays confirm that melanoma cells could drive M2 polarisation and SPP1 upregulation in macrophages. Knocking down or overexpressing SPP1 correspondingly alters M2 gene expression in macrophages.

This study establishes a pan‐skin cancer immune remodelling framework, providing a foundation for biomarker discovery and the development of new immunotherapy strategies.

Graphical highlights

Tumour subpopulations with high malignancy show reduced interferon signalling and reduced MHC‐I expression.

Tumour‐associated macrophages (TAMs) differentiate into two opposing states: SPP1+ TAMs and CXCL9+ TAMs.

Among immune ecotypes, T‐cell‐enriched ecotypes are linked with high MHC‐I expression in tumour cells while desert ecotypes are accompanied by high SPP1+ TAMs.

SPP1 regulates M2 markers in macrophages validated by functional assays.