Analysis of Intercellular Communication in the Osteosarcoma Microenvironment Based on Single Cell Sequencing Data

Osteosarcoma (OS) is the most common primary bone cancer in children and young adults, patient survival rates have not improved in recent decades. To further understand the interrelationship between different cell types in the tumor microenvironment of osteosarcoma, we comprehensively analyzed single-cell sequencing data from six patients with untreated osteosarcoma. Nine major cell types were identified from a total of 46,046 cells based on unbiased clustering of gene expression profiles and canonical markers. Osteosarcoma from different patients display heterogeneity in cellular composition. Myeloid cells were the most commonly represented cell type, followed by osteoblastic and TILs. Copy number variation (CNV) results identified amplifications and deletions in malignant osteoblastic cells and fibroblasts. Trajectory analysis based on RNA velocity showed that osteoclasts in the OS microenvironment could be differentiated from myeloid cells. Furthermore, we explored the intercellular communications in OS microenvironment and identified multiple ligand-receptor pairs between myeloid cells, osteoblastic cells and their cells, including 21 ligand-receptor pair genes that significantly associated with survival outcomes. Importantly, we found chemotherapy may have an effect on cellular communication in the OS microenvironment by analyzing single-cell sequencing data from seven primary osteosarcoma patients who received chemotherapy. We believe these observations will improve our understanding of potential mechanisms of microenvironment contributions to OS progression and help identify potential targets for new treatment development in the future.