Abstract
BACKGROUND: Osteosarcoma (OS) patients have a poor response to immunotherapy due to the sheer complexity of the immune system and the nuances of the tumor-immune microenvironment. Methodology. To gain insights into the immune heterogeneity of OS, we identified robust clusters of patients based on the immune gene expression profiles of OS patients in the TARGET database and assessed their reproducibility in an independent cohort collected from the GEO database. The association of comprehensive molecular characterization with reproducible immune subtypes was accessed with ANOVA. Furthermore, we visualized the distribution of individual patients in a tree structure by the graph structure learning-based dimensionality reduction algorithm. RESULTS: We found that 87 OS samples can be divided into 5 immune subtypes, and each of them was associated with distinct clinical outcomes. The immune subtypes also demonstrated widely different patterns in tumor genetic aberrations, tumor-infiltrating, immune cell composition, and cytokine profiles. The immune landscape of OS uncovered the significant intracluster heterogeneity within each subtype and depicted a continuous immune spectrum across patients. CONCLUSION: The established five immune subtypes in our study suggested immune heterogeneity in OS patients and may provide optimal individual immunotherapy for patients exhibiting OS.