Comprehensive analysis of angiogenesis pattern and related immune landscape for individual treatment in osteosarcoma

Postoperative recurrence and metastasis are the main reasons for the poor prognosis of osteosarcoma (OS). Currently, an ideal predictor for not only prognosis but also drug sensitivity and immunotherapy responses in OS patients is urgently needed. Angiogenesis plays a crucial role in tumour progression, which suggests its immense potential for predicting prognosis and responses to immunotherapy for OS. Angiogenesis patterns in OS were explored in depth in this study to construct a prognostic model called ANGscore and clarify the underlying mechanism involved in the immune microenvironment. The efficacy and robustness of the model were validated in multiple datasets, including bulk RNA-seq datasets (TARGET-OS, GSE21257), a single-cell RNA-seq dataset (GSE152048) and immunotherapy-related datasets (GSE91061, GSE173839). OS patients with a high ANGscore had a worse prognosis, accompanied by the immune desert phenotype. Pseudotime and cellular communication analyses in scRNA-seq data revealed that as the ANGscore increased, the malignant degree of cells increased, and IFN-γ signalling was involved in tumour progression and regulation of the tumour immune microenvironment. Furthermore, the ANGscore was associated with immune cell infiltration and the response rate to immunotherapy. OS patients with high ANGscore might be resistant to uprosertib, and be sensitive to VE821, AZD6738 and BMS.345541. In conclusion, we established a novel ANGscore system by comprehensively analysing the expression pattern of angiogenesis genes, which can accurately differentiate the prognosis and immune characteristics of OS populations. Additionally, the ANGscore can be used for patient stratification during immunotherapy, and guide individualized treatment strategies.