Abstract
Tumor immunity is closely associated with the prognosis of tumors, including osteosarcoma (OS). The aim of the present study was to construct an immune-related prognostic index (PI) to predict the prognosis of OS. Herein, OS expression data were sourced from the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database. We divided the OS patients into nonmetastatic and metastatic groups, allowing differentially immune-related genes (DIRGs) to be selected. After DIRGs were further investigated by enrichment analysis, four keys prognostic IRGs (CD79A, CSF3R, MTNR1B and NPPC) were identified using a Cox proportional hazards model. Then, an immune-related prognostic index was constructed. Finally, gene set enrichment analysis (GSEA) was employed to further explore the underlying mechanisms. The difference in tumor-infiltrating immune cell (TIIC) abundance was also discussed. In our study, eight upregulated genes and 30 downregulated genes were identified. Several Gene Ontology (GO) terms and the most significantly enriched KEGG pathways were immune-associated functions and pathways. Four genes, including CD79A, CSF3R, MTNR1B and NPPC, were used to establish a risk assessment model for evaluating OS prognosis. GSEA revealed that the risk score was related to cytokine receptor interaction and to the chemokine and B cell receptor signaling pathways. Furthermore, high risk markedly related to the infiltration of several immune cell types, including M2 macrophages, naïve CD4 T cells, and CD8 T cells. In sum, we developed a survival model for OS. The underlying molecular mechanisms of the high-risk group may affect immune-related biological processes and TIICs.Abbreviations TARGET: Therapeutically Applicable Research To Generate Effective Treatments; PI: Prognostic index; OS: Osteosarcoma; DIRGs: Differentially immune-related genes; GSEA: Gene set enrichment analysis; TIIC: Tumor-infiltrating immune cell.