Abstract:Background: Melanoma is typically regarded as the most dangerous form of skin cancer. Although surgical removal of in situ lesions can be used to effectively treat metastatic disease, this condition is still difficult to cure. Melanoma cells are removed in great part due to the action of natural killer (NK) and T cells on the immune system. Still, not much is known about how the activity of NK cell-related pathways changes in melanoma tissue. Thus, we performed a single-cell multi-omics analysis on human melanoma cells in this study to explore the modulation of NK cell activity. Materials and methods: Cells in which mitochondrial genes comprised > 20% of the total number of expressed genes were removed. Gene ontology (GO), gene set enrichment analysis (GSEA), gene set variation analysis (GSVA), and AUCcell analysis of differentially expressed genes (DEGs) in melanoma subtypes were performed. The CellChat package was used to predict cell-cell contact between NK cell and melanoma cell subtypes. Monocle program analyzed the pseudotime trajectories of melanoma cells. In addition, CytoTRACE was used to determine the recommended time order of melanoma cells. InferCNV was utilized to calculate the CNV level of melanoma cell subtypes. Python package pySCENIC was used to assess the enrichment of transcription factors and the activity of regulons in melanoma cell subtypes. Furthermore, the cell function experiment was used to confirm the function of TBX21 in both A375 and WM-115 melanoma cell lines. Results: Following batch effect correction, 26,161 cells were separated into 28 clusters and designated as melanoma cells, neural cells, fibroblasts, endothelial cells, NK cells, CD4+ T cells, CD8+ T cells, B cells, plasma cells, monocytes and macrophages, and dendritic cells. A total of 10137 melanoma cells were further grouped into seven subtypes, i.e., C0 Melanoma BIRC7, C1 Melanoma CDH19, C2 Melanoma EDNRB, C3 Melanoma BIRC5, C4 Melanoma CORO1A, C5 Melanoma MAGEA4, and C6 Melanoma GJB2. The results of AUCell, GSEA, and GSVA suggested that C4 Melanoma CORO1A may be more sensitive to NK and T cells through positive regulation of NK and T cell-mediated immunity, while other subtypes of melanoma may be more resistant to NK cells. This suggests that the intratumor heterogeneity (ITH) of melanoma-induced activity and the difference in NK cell-mediated cytotoxicity may have caused NK cell defects. Transcription factor enrichment analysis indicated that TBX21 was the most important TF in C4 Melanoma CORO1A and was also associated with M1 modules. In vitro experiments further showed that TBX21 knockdown dramatically decreases melanoma cells' proliferation, invasion, and migration. Conclusion: The differences in NK and T cell-mediated immunity and cytotoxicity between C4 Melanoma CORO1A and other melanoma cell subtypes may offer a new perspective on the ITH of melanoma-induced metastatic activity. In addition, the protective factors of skin melanoma, STAT1, IRF1, and FLI1, may modulate melanoma cell responses to NK or T cells.

Single-cell RNA-sequencing analyses identify heterogeneity of CD8+ T cell subpopulations and novel therapy targets in melanoma

Regulatory T cells subgroups in the tumor microenvironment cannot be overlooked: Their involvement in prognosis and treatment strategy in melanoma

Single-Cell Transcriptomic Analysis Reveals the Crosstalk Propensity Between the Tumor Intermediate State and the CD8+ T Exhausted State to be Associated with Clinical Benefits in Melanoma

Generation, Transcriptomic States, and Clinical Relevance of CX3CR1+ CD8 T Cells in Melanoma

The combination of single-cell and RNA sequencing analysis decodes the melanoma tumor microenvironment and identifies novel T cell-associated signature genes

Advancing immunotherapy for melanoma: the critical role of single-cell analysis in identifying predictive biomarkers

Immune and non-immune cell subtypes identify novel targets for prognostic and therapeutic strategy: A study based on intratumoral heterogenicity analysis of multicenter scRNA-seq datasets in lung adenocarcinoma

Establishment of a CD8+ T cells-related prognostic risk model for acral melanoma based on single-cell and bulk RNA sequencing

Bioinformatics analysis of the tumor microenvironment in melanoma - Constructing a prognostic model based on CD8+ T cell-related genes: An observational study

Integrated analysis reveals the pivotal interactions between immune cells in the melanoma tumor microenvironment

Pivotal factors associated with the immunosuppressive tumor microenvironment and melanoma metastasis

Integrative analysis of bulk and single-cell RNA-seq reveals the molecular characterization of the immune microenvironment and oxidative stress signature in melanoma

Phenotypic and spatial heterogeneity of CD8+ tumour infiltrating lymphocytes

Heterogeneity of Exhausted T Cells in the Tumor Microenvironment is Linked to Patient Survival Following Resection in Hepatocellular Carcinoma

A Burned-Out CD8+ T-cell Subset Expands in the Tumor Microenvironment and Curbs Cancer Immunotherapy

A single-cell analysis reveals tumor heterogeneity and immune environment of acral melanoma

An advanced comprehensive muti-cell-type-specific model for predicting anti-PD-1 therapeutic effect in melanoma

Single-cell transcriptome analysis of the heterogeneous effects of differential expression of tumor PD-L1 on responding TCR-T cells

Abstract A019: The role of tumor transcriptional variability in evading acute CD8+ T cell exposure in melanoma

A single-cell atlas reveals tumor heterogeneity and immune environment of acral melanoma

Innovative breakthroughs facilitated by single-cell multi-omics: manipulating natural killer cell functionality correlates with a novel subcategory of melanoma cells