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.

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

Single-Cell Characterization of the Immune Microenvironment of Melanoma Brain and Leptomeningeal Metastases.

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

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

Creating a multifaceted prognostic model for cutaneous melanoma: the convergence of single-cell and bulk sequencing with machine learning

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

Bioinformatics Analysis to Identify the Critical Genes, Micrornas and Long Noncoding RNAs in Melanoma

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

Multiplexed Imaging Mass Cytometry Reveals Distinct Tumor-Immune Microenvironments Linked to Immunotherapy Responses in Melanoma

Tumor Immune Microenvironment Characterization Identifies Prognosis and Immunotherapy-Related Gene Signatures in Melanoma

Bioinformatic Analysis Identifies Potential Key Genes in the Pathogenesis of Melanoma

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

Decoding Immune Heterogeneity of Melanoma and identifying immune-prognostic hub genes

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

A Novel Gene Signature Based on Immunogenic Cell Death-Related Genes Predicts the Prognosis and Immune Infiltration Status of Melanoma Patients

Unveiling the link between inflammasomes and skin cutaneous melanoma: Insights into expression patterns and immunotherapy response prediction.

Single-Cell Profiling Reveals Global Immune Responses during the Progression of Murine Epidermal Neoplasms

Pivotal factors associated with the immunosuppressive tumor microenvironment and melanoma metastasis

Exploring Tumor Immune Microenvironment and Its Associations With Molecular Characteristics in Melanoma

Exploring a Specialized Programmed-Cell Death Patterns to Predict the Prognosis and Sensitivity of Immunotherapy in Cutaneous Melanoma Via Machine Learning.

Detailed Analysis of Molecular Mechanisms in Primary and Metastatic Melanoma.