Abstract:Background: The immune microenvironment and oxidative stress of melanoma show significant heterogeneity, which affects tumor growth, invasion and treatment response. Single-cell and bulk RNA-seq data were used to explore the heterogeneity of the immune microenvironment and oxidative stress of melanoma.Methods: The R package Seurat facilitated the analysis of the single-cell dataset, while Harmony, another R package, was employed for batch effect correction. Cell types were classified using Uniform Manifold Approximation and Projection (UMAP). The Secreted Signaling algorithm from CellChatDB.human was applied to elucidate cell-to-cell communication patterns within the single-cell data. Consensus clustering analysis for the skin cutaneous melanoma (SKCM) samples was executed with the R package ConsensusClusterPlus. To quantify immune infiltrating cells, we utilized CIBERSORT, ESTIMATE, and TIMERxCell algorithms provided by the R package Immuno-Oncology Biological Research (IOBR). Single nucleotide variant (SNV) analysis was conducted using Maftools, an R package specifically designed for this purpose. Subsequently, the expression levels of PXDN and PAPSS2 genes were assessed in melanoma tissues compared to adjacent normal tissues. Furthermore, in vitro experiments were conducted to evaluate the proliferation and reactive oxygen species expression in melanoma cells following transfection with siRNA targeting PXDN and PAPSS2.Results: Malignant tumor cell populations were reclassified based on a comprehensive single-cell dataset analysis, which yielded six distinct tumor subsets. The specific marker genes identified for these subgroups were then used to interrogate the Cancer Genome Atlas Skin Cutaneous Melanoma (TCGA-SKCM) cohort, derived from bulk RNA sequencing data, resulting in the delineation of two immune molecular subtypes. Notably, patients within the cluster2 (C2) subtype exhibited a significantly more favorable prognosis compared to those in the cluster1 (C1) subtype. An alignment of immune characteristics was observed between the C2 subtype and unique immune functional tumor cell subsets. Genes differentially expressed across these subtypes were subsequently leveraged to construct a predictive risk model. In vitro investigations further revealed elevated expression levels of PXDN and PAPSS2 in melanoma tissue samples. Functional assays indicated that modulation of PXDN and PAPSS2 expression could influence the production of reactive oxygen species (ROS) and the proliferative capacity of melanoma cells.Conclusion: The constructed six-gene signature can be used as an immune response and an oxidative stress marker to guide the clinical diagnosis and treatment of melanoma.

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

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

Single-cell Characterization of the Cellular Landscape of Acral Melanoma Identifies Novel Targets for Immunotherapy

Supplementary Figure from Single-cell Characterization of the Cellular Landscape of Acral Melanoma Identifies Novel Targets for Immunotherapy

Whole Genome and RNA Sequencing Reveal the Distinct Genomic Landscape of Acral Melanoma.

Mapping the single-cell landscape of acral melanoma and analysis of the molecular regulatory network of the tumor microenvironments

Single-cell RNA sequencing unveils tumor heterogeneity and immune microenvironment between subungual and plantar melanoma

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

Integrative Genomic Profiling Uncovers Therapeutic Targets of Acral Melanoma in Asian Populations

Identification of two molecularly and prognostically distinct subtypes in acral melanoma using network prediction method

Single-Cell RNA Sequencing Reveals Cellular Heterogeneity in an Acral Amelanotic Melanoma after Immunotherapy Treatment.

Network-based Molecular Subtyping of Acral Melanoma

Single-cell profiling of acral melanoma infiltrating lymphocytes reveals a suppressive tumor microenvironment

Integrative molecular and clinical profiling of acral melanoma links focal amplification of 22q11.21 to metastasis

Genomic Heterogeneity and Branched Evolution of Early Stage Primary Acral Melanoma Shown by Multiregional Microdissection Sequencing.

The genetic evolution of acral melanoma

Prognostic Value of Genetic Aberrations and Tumor Immune Microenvironment in Primary Acral Melanoma

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

Comparative transcriptomic analysis reveals differences in gene expression and regulatory pathways between nonacral and acral melanoma in Asian individuals

Delineating the Early Dissemination Mechanisms of Acral Melanoma by Integrating Single-Cell and Spatial Transcriptomic Analyses

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