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.

Correction: Clinical-mediated discovery of pyroptosis in CD8 + T cell and NK cell reveals melanoma heterogeneity by single-cell and bulk sequence

Author Correction: PD-L1-mediated gasdermin C expression switches apoptosis to pyroptosis in cancer cells and facilitates tumour necrosis

Multiomics characterization of pyroptosis in the tumor microenvironment and therapeutic relevance in metastatic melanoma

Publisher Correction: Global characterization of T cells in non-small-cell lung cancer by single-cell sequencing

Correction to: The phosphorylation of CD147 by Fyn plays a critical role for melanoma cells growth and metastasis

Publisher Correction: Molecular analysis of primary melanoma T cells identifies patients at risk for metastatic recurrence

Correction To: Apoptosis of Murine Melanoma B16-BL6 Cells Induced by Quercetin Targeting Mitochondria, Inhibiting Expression of PKC-α and Translocating PKC-δ

Exploring the role of pyroptosis in shaping the tumor microenvironment of colorectal cancer by bulk and single-cell RNA sequencing

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

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

Author Correction: Neoadjuvant Immune Checkpoint Blockade in High-Risk Resectable Melanoma

Pan-cancer analyses of pyroptosis with functional implications for prognosis and immunotherapy in cancer

A new pyroptosis model can predict the immunotherapy response and immune microenvironment characteristics and prognosis of patients with cutaneous melanoma based on TCGA and GEO databases

Corrigendum: Identification of Breast Cancer Immune Subtypes by Analyzing Bulk Tumor and Single Cell Transcriptomes

Pyroptosis-related gene mediated modification patterns and immune cell infiltration landscapes in cutaneous melanoma to aid immunotherapy

Integrated bulk and single-cell transcriptomes reveal pyroptotic signature in prognosis and therapeutic options of hepatocellular carcinoma by combining deep learning

Author Correction: Cutaneous T cell lymphoma atlas reveals malignant T H 2 cells supported by a B cell-rich tumor microenvironment

Correction: Corrigendum: CRISPR-Cas9 mediated efficient PD-1 disruption on human primary T cells from cancer patients

Identification and Validation of Two Heterogeneous Molecular Subtypes and a Prognosis Predictive Model for Hepatocellular Carcinoma Based on Pyroptosis

Correction To: Predictive Biomarkers and Mechanisms Underlying Resistance to PD1/PD-L1 Blockade Cancer Immunotherapy.

The expression pattern of pyroptosis-related genes predicts the prognosis and drug response of melanoma