The immune cell infiltration-associated molecular subtypes and gene signature predict prognosis for osteosarcoma patients.

Host immune dysregulation involves in the initiation and development of osteosarcoma (OS). However, the exact role of immune cells in OS remains unknown. We aimed to distinguish the molecular subtypes and establish a prognostic model in OS patients based on immunocyte infiltration. The gene expression profile and corresponding clinical feature of OS patients were obtained from TARGET and GSE21257 datasets. MCP-counter and univariate Cox regression analyses were applied to identify immune cell infiltration-related molecular subgroups. Functional enrichment analysis and immunocyte infiltration analysis were performed between two subgroups. Furthermore, Cox regression and LASSO analyses were performed to establish the prognostic model for the prediction of prognosis and metastasis in OS patients. The subgroup with low infiltration of monocytic lineage (ML) was related to bad prognosis in OS patients. 435 DEGs were screened between the two subgroups. Functional enrichment analysis revealed these DEGs were involved in immune- and inflammation-related pathways. Three important genes (including TERT, CCDC26, and IL2RA) were identified to establish the prognostic model. The risk model had good prognostic performance for the prediction of metastasis and overall survival in OS patients. A novel stratification system was established based on ML-related signature. The risk model could predict the metastasis and prognosis in OS patients. Our findings offered a novel sight for the prognosis and development of OS.