Identification and Development of a Novel 4-Gene Immune-Related Signature to Predict Osteosarcoma Prognosis

Osteosarcoma (OS) is a malignant disease that develops rapidly and is associated with poor prognosis. Immunotherapy may provide new insights into clinical treatment strategies for OS. The purpose of this study was to identify immune-related genes that could predict OS prognosis. The gene expression profiles and clinical data of 84 OS patients were obtained from the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database. According to non-negative matrix factorization, two molecular subtypes of immune-related genes, C1 and C2, were acquired, and 597 differentially expressed genes between C1 and C2 were identified. Univariate Cox analysis was performed to get 14 genes associated with survival, and 4 genes (GJA5, APBB1IP, NPC2, and FKBP11) obtained through least absolute shrinkage and selection operator (LASSO)-Cox regression were used to construct a 4-gene signature as a prognostic risk model. The results showed that high FKBP11 expression was correlated with high risk (a risk factor), and that high GJA5, APBB1IP, or NPC2 expression was associated with low risk (protective factors). The testing cohort and entire TARGET cohort were used for internal verification, and the independent GSE21257 cohort was used for external validation. The study suggested that the model we constructed was reliable and performed well in predicting OS risk. The functional enrichment of the signature was studied through gene set enrichment analysis, and it was found that the risk score was related to the immune pathway. In summary, our comprehensive study found that the 4-gene signature could be used to predict OS prognosis, and new biomarkers of great significance for understanding the therapeutic targets of OS were identified.