Deciphering the heterogeneity and immunosuppressive function of regulatory T cells in osteosarcoma using single-cell RNA transcriptome

Debin Cheng; Zhao Zhang; Zhenzhou Mi; Weidong Tao; Dong Liu; Jun Fu; Hongbin Fan

doi:10.1016/j.compbiomed.2023.107417

Deciphering the heterogeneity and immunosuppressive function of regulatory T cells in osteosarcoma using single-cell RNA transcriptome

Comput Biol Med. 2023 Oct:165:107417. doi: 10.1016/j.compbiomed.2023.107417. Epub 2023 Sep 1.

Authors

Debin Cheng¹, Zhao Zhang¹, Zhenzhou Mi¹, Weidong Tao¹, Dong Liu¹, Jun Fu¹, Hongbin Fan²

Affiliations

¹ Department of Orthopaedic Surgery, Xi Jing Hospital, The Fourth Military Medical University, Xi'an, China.
² Department of Orthopaedic Surgery, Xi Jing Hospital, The Fourth Military Medical University, Xi'an, China. Electronic address: fanhb@fmmu.edu.cn.

PMID: 37669584
DOI: 10.1016/j.compbiomed.2023.107417

Abstract

Osteosarcoma (OS) is a highly invasive malignant neoplasm with poor prognosis. The tumor microenvironment (TME) plays an essential role in the occurrence and development of OS. Regulatory T cells (Tregs) are known to facilitate immunosuppression, tumor progression, invasion, and metastasis. However, the effect of Tregs in the TME of OS remains unclear. In this study, single-cell RNA sequencing (scRNA-seq) data was used to identify Tregs and various other cell clusters in the TME of OS. Gene set variation analysis (GSVA) was used to investigate the signaling pathways in Tregs from OS and adjacent tissues. The CellChat and iTALK packages were used to analyze cellular communication. In addition, a prognostic model was established based on the Tregs-specific genes using bulk RNA-seq from the TARGET database, and it was verified using a Gene Expression Omnibus dataset. The pRRophetic package was used to predict drug sensitivity. Immunohistochemistry was used to verify the expression of candidate genes in OS. Based on the above methods, we showed that the OS samples were highly infiltrated with Tregs. GSVA revealed that oxidative phosphorylation, angiogenesis and mammalian target of rapamycin complex 1 (mTORC1) were highly activated in Tregs from OS compared with those from adjacent tissues. Using cellular communication analysis, we found that Tregs interacted with osteoblastic, endothelial, and myeloid cells via C-X-C motif chemokine ligand (CXCL) signaling; particularly, they strongly affected the expression of C-X-C motif chemokine receptor 4 (CXCR4) and interacted with other cell clusters through CXCL12/transforming growth factor β1 (TGFB1) to collectively enable tumor growth and progression. Subsequently, two Tregs-specific genes-CD320 and MAF-were screened through univariate, least absolute shrinkage and selection operator regression (LASSO) and multivariate analysis to construct a prognostic model, which showed excellent prognostic accuracy in two independent cohorts. In addition, drug sensitivity analysis demonstrated that OS patients at high Tregs risk were sensitive to sunitinib, sorafenib, and axitinib. We also used immunohistochemistry to validate that CD320 and MAF were significantly upregulated in OS tissues compared with adjacent tissues. Overall, this study reveals the heterogeneity of Tregs in the OS TME, providing new insights into the invasion and treatment of this cancer.

Keywords: Osteosarcoma; Prognosis; Regulatory T cells; Single cell RNA sequencing; Tumor microenvironment.