Analysis of Immune-Stromal Score-Based Gene Signature and Molecular Subtypes in Osteosarcoma: Implications for Prognosis and Tumor Immune Microenvironment

Objective: Infiltrating immune and stromal cells are essential for osteosarcoma progression. This study set out to analyze immune-stromal score-based gene signature and molecular subtypes in osteosarcoma. Methods: The immune and stromal scores of osteosarcoma specimens from the TARGET cohort were determined by the ESTIMATE algorithm. Then, immune-stromal score-based differentially expressed genes (DEGs) were screened, followed by univariate Cox regression analysis. A LASSO regression analysis was applied for establishing a prognostic model. The predictive efficacy was verified in the GSE21257 dataset. Associations between the risk scores and chemotherapy drug sensitivity, immune/stromal scores, PD-1/PD-L1 expression, immune cell infiltrations were assessed in the TARGET cohort. NMF clustering analysis was employed for characterizing distinct molecular subtypes based on immune-stromal score-based DEGs. Results: High immune/stromal scores exhibited the prolonged survival duration of osteosarcoma patients. Based on 85 prognosis-related stromal-immune score-based DEGs, a nine-gene signature was established. High-risk scores indicated undesirable prognosis of osteosarcoma patients. The AUCs of overall survival were 0.881 and 0.849 in the TARGET cohort and GSE21257 dataset, confirming the well predictive performance of this signature. High-risk patients were more sensitive to doxorubicin and low-risk patients exhibited higher immune/stromal scores, PD-L1 expression, and immune cell infiltrations. Three molecular subtypes were characterized, with distinct clinical outcomes and tumor immune microenvironment. Conclusion: This study developed a robust prognostic gene signature as a risk stratification tool and characterized three distinct molecular subtypes for osteosarcoma patients based on immune-stromal score-based DEGs, which may assist decision-making concerning individualized therapy and follow-up project.

Automated summary

High immune/stromal scores were associated with prolonged survival duration in osteosarcoma patients. A nine-gene signature was established based on 85 prognosis-related stromal-immune score-based DEGs, with high-risk scores indicating unfavorable prognosis for osteosarcoma patients. The predictive performance of this signature was confirmed by AUC values of 0.881 and 0.849 in the TARGET cohort and GSE21257 dataset.