Exploring the impact of PDGFD in osteosarcoma metastasis through single-cell sequencing analysis
Authors: Huang Yujing , Cao Dongyan , Zhang Manxue , Yang Yue , Niu Gengming , Tang Lina , Shen Zan , Zhang Zhichang , Bai Yueqing , Min Daliu , He Aina 
Published: April 2024

Abstract
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.