Unraveling intratumoral complexity in metastatic dermatofibrosarcoma protuberans through single-cell RNA sequencing analysis

Dermatofibrosarcoma protuberans (DFSP) stands as a rare and locally aggressive soft tissue tumor, characterized by intricated molecular alterations. The imperative to unravel the complexities of intratumor heterogeneity underscores effective clinical management. Herein, we harnessed single-cell RNA sequencing (scRNA-seq) to conduct a comprehensive analysis encompassing samples from primary sites, satellite foci, and lymph node metastases. Rigorous preprocessing of raw scRNA-seq data ensued, and employing t-distributed stochastic neighbor embedding (tSNE) analysis, we unveiled seven major cell populations and fifteen distinct subpopulations. Malignant cell subpopulations were delineated using infercnv for copy number variation calculations. Functional and metabolic variations of diverse malignant cell populations across samples were deciphered utilizing GSVA and the scMetabolism R packages. Additionally, the exploration of differentiation trajectories within diverse fibroblast subpopulations was orchestrated through pseudotime trajectory analyses employing CytoTRACE and Monocle2, and further bolstered by GO analyses to elucidate the functional disparities across distinct differentiation states. In parallel, we segmented the cellular components of the immune microenvironment and verified the presence of SPP1 + macrophage, which constituted the major constituent in lymph node metastases. Remarkably, the CellChat facilitated a comprehensive intercellular communication analysis. This study culminates in an all-encompassing single-cell transcriptome atlas, propounding novel insights into the multifaceted nature of intratumor heterogeneity and fundamental molecular mechanisms propelling metastatic DFSP.