Exploring the Impact of PDGFD in Osteosarcoma Metastasis Through Single-cell Sequencing Analysis

Purpose: The overall survival rate for metastatic osteosarcoma hovers around 20%. Responses to second-line chemotherapy, targeted therapies, and immunotherapies have demonstrated limited efficacy in metastatic osteosarcoma. Our objective is to validate differentially expressed genes and signaling pathways between non-metastatic and metastatic osteosarcoma, employing single-cell RNA sequencing (scRNA-seq) and additional functional investigations. We aim to enhance comprehension of metastatic mechanisms and potentially unveil a therapeutic target.

Methods: scRNA-seq was performed on two primary osteosarcoma lesions (1 non-metastatic and 1 metastatic). Seurat package facilitated dimensionality reduction and cluster identification. Copy number variation (CNV) was predicted using InferCNV. CellChat characterized ligand-receptor-based intercellular communication networks. Differentially expressed genes underwent GO function enrichment analysis and GSEA. Validation was achieved through the GSE152048 dataset, which identified PDGFD-PDGFRB as a common ligand-receptor pair with significant contribution. Immunohistochemistry assessed PDGFD and PDGFRB expression, while multicolor immunofluorescence and flow cytometry provided insight into spatial relationships and the tumor immune microenvironment. Kaplan-Meier survival analysis compared metastasis-free survival and overall survival between high and low levels of PDGFD and PDGFRB. Manipulation of PDGFD expression in primary osteosarcoma cells examined invasion abilities and related markers.

Results: Ten clusters encompassing osteoblasts, osteoclasts, osteocytes, fibroblasts, pericytes, endothelial cells, myeloid cells, T cells, B cells, and proliferating cells were identified. Osteoblasts, osteoclasts, and osteocytes exhibited heightened CNV levels. Ligand-receptor-based communication networks exposed significant fibroblast crosstalk with other cell types, and the PDGF signaling pathway was activated in non-metastatic osteosarcoma primary lesion. These results were corroborated by the GSE152048 dataset, confirming the prominence of PDGFD-PDGFRB as a common ligand-receptor pair. Immunohistochemistry demonstrated considerably greater PDGFD expression in non-metastatic osteosarcoma tissues and organoids, correlating with extended metastasis-free and overall survival. PDGFRB expression showed no significant variation between non-metastatic and metastatic osteosarcoma, nor strong correlations with survival times. Multicolor immunofluorescence suggested co-localization of PDGFD with PDGFRB. Flow cytometry unveiled a highly immunosuppressive microenvironment in metastatic osteosarcoma. Manipulating PDGFD expression demonstrated altered invasive abilities and marker expressions in primary osteosarcoma cells from both non-metastatic and metastatic lesions.

Conclusions: scRNA-seq illuminated the activation of the PDGF signaling pathway in primary lesion of non-metastatic osteosarcoma. PDGFD displayed an inhibitory effect on osteosarcoma metastasis, likely through the suppression of the EMT signaling pathway.