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.

Cell type identification from single-cell transcriptomes in melanoma

Identification of Transcriptional Heterogeneity and Construction of a Prognostic Model for Melanoma Based on Single-Cell and Bulk Transcriptome Analysis

P02.05 Comprehensive Profiling of Tumor Heterogeneity and Its Microenvironment in Advanced Non-Small Cell Lung Cancer at Single Cell Resolution

Identification of Molecular Subtypes for Integrated Multi-Omics Analysis for Use in Guiding Precision Medicine in Hepatocellular Carcinoma.

Unravelling tumour cell diversity and prognostic signatures in cutaneous melanoma through machine learning analysis

High-Dimensional Single-Cell Transcriptomics in Melanoma and Cancer Immunotherapy

Bioinformatics data combined with single-cell analysis reveals patterns of immunoinflammatory infiltration and cell death in melanoma

Integration of Bulk RNA Sequencing and Single-Cell RNA Sequencing to Reveal Uveal Melanoma Tumor Heterogeneity and Cells Related to Survival

Identifying Cell-Type Specific Genes and Expression Rules Based on Single-Cell Transcriptomic Atlas Data

Single-cell RNA sequencing reveals melanoma cell state-dependent heterogeneity of response to MAPK inhibitors

Discovery and Validation of a Metastasis-Related Prognostic and Diagnostic Biomarker for Melanoma Based on Single Cell and Gene Expression Datasets.

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

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

A fine-grained cell type deconvolution of bulk transcriptome from single-cell RNA-seq data reveals myeloid cells heterogeneity in lung adenocarcinomas

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

Identification of cell-type-specific genes in multimodal single-cell data using deep neural network algorithm

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

Identification of feature genes and molecular mechanisms involved in cell communication in uveal melanoma through analysis of single‐cell sequencing data

Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq

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

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