Impact Factor 6.684 | CiteScore 2.7
More on impact ›


Front. Cell Dev. Biol., 16 July 2021 |

Classification of Osteosarcoma Based on Immunogenomic Profiling

Xinwen Wang1,2†, Liangming Wang3†, Weifeng Xu4†, Xinwu Wang5, Dianshan Ke6, Jinluan Lin7, Wanzun Lin8* and Xiaochun Bai1,2*
  • 1The Third Affiliated Hospital of Southern Medical University, Guangzhou, China
  • 2Guangdong Provincial Key Laboratory of Bone and Joint Degeneration Diseases, Guangzhou, China
  • 3Department of Orthopedics, The Second Affiliated Hospital of Fujian Medical University, Quanzhou, China
  • 4Department of Medical Oncology, The Affiliated Cancer Hospital of Zhengzhou University, Zhengzhou, China
  • 5Department of Orthopedics, The First Hospital of Putian City, Putian, China
  • 6Department of Orthopedics, Jiangmen People’s Hospital, Jiangmen, China
  • 7Department of Orthopedics, The First Affiliated Hospital of Fujian Medical University, Fuzhou, China
  • 8Department of Radiation Oncology, Shanghai Proton and Heavy Ion Center, Fudan University Cancer Hospital, Shanghai, China

Accumulating evidence has supported that osteosarcoma is heterogeneous, and several subtypes have been identified based on genomic profiling. Immunotherapy is revolutionizing cancer treatment and is a promising therapeutic strategy. In contrast, few studies have identified osteosarcoma classification based on immune biosignatures, which offer the optimal stratification of individuals befitting immunotherapy. Here, we classified osteosarcoma into two clusters: immunity high and immunity low using the single-sample gene-set enrichment analysis and unsupervised hierarchical clustering. Immunity_H subtype was associated with high immune cells infiltration, a favorable prognosis, benefit to immunotherapy, high human leukocyte antigen gene expression, and activated immune signal pathway indicating an immune-hot phenotype. On the contrary, the Immunity_L subtype was correlated with low immune cell infiltration, poor prognosis, and cancer-related pathway, indicating an immune-cold phenotype. We also identified TYROBP as a key immunoregulatory gene associated with CD8+ T cell infiltration by multiplex immunohistochemistry. Finally, we established an immune-related prognostic model that predicted the survival time of osteosarcoma. In conclusion, we established a new classification system of osteosarcoma based on immune signatures and identified TYROBP as a key immunoregulatory gene. This stratification had significant clinical outcomes for estimating prognosis, as well as the immunotherapy of osteosarcoma patients.


Osteosarcoma, originating from mesenchymal stem cells, is the most common primary malignancy of bone (ElKordy et al., 2018). With the advancement of surgical resection, radiotherapy, and neoadjuvant chemotherapy, the survival rate of osteosarcoma patients has increased up to 60–70% (Kansara et al., 2014; Simpson and Brown, 2018). However, the prognosis of patients with recurrent, metastatic, or unresectable osteosarcomas is extremely poor. The long-term survival rate for patients with primary osteosarcoma is approximately 65%, whereas it is less than 20% for patients with metastatic osteosarcomas (Durfee et al., 2016; Whelan and Davis, 2018). Hence, a novel therapeutic strategy for advanced sarcomas should be developed to improve the outcome of osteosarcomas. Recently, cancer immunotherapy has received tremendous attention for the treatment of numerous refractory malignancies (Baxevanis et al., 2009; Yang, 2015; Emens et al., 2017; O’Donnell et al., 2019; Lin et al., 2020). Hence, it is worth considering immunotherapy for osteosarcoma, as the treatment options for this disease are still limited.

Accumulating evidence has suggested that osteosarcoma tends to be susceptible to immunotherapy (Chen et al., 2021). Osteosarcoma tumors have a high proportion of CD8+ invading lymphocytes relative to other sarcoma subtypes, and the number of infiltrating immune cells associates positively with overall survival (OS) (Wang et al., 2016). Besides, some osteosarcomas have an elevated programmed cell death protein-1 ligand (PD-L1) expression, with a high level of genomic instability, indicating prospective sensitivity to programmed cell death protein-1 (PD-1)/PD-L1 cascade suppressors (Shen et al., 2014). Although some patients benefit from different immunotherapeutic interventions, the majority fail to experience benefit. Thus, considering the low response rate, novel and accurate classification methods should be developed to identify patients befitting immunotherapy.

Here, we performed the single-sample gene-set enrichment analysis (ssGSEA) of the transcriptome data of 29 immune gene sets and established a novel classification of osteosarcoma. We also tried to characterize the subtype-distinct molecular features, including immune signaling cascades, and key genes and to construct and validate the immune-associated prognostic biosignature.

Materials and Methods


The RNA-seq transcriptome data, as well as the clinicopathological information of osteosarcoma patients, were obtained from TARGET as a training set. Similarly, patients from Gene Expression Omnibus (GSE21257) were abstracted as a validation set.

Single-Sample Gene-Set Enrichment Analysis

For each osteosarcoma dataset, the levels of enrichment of the 29 immune biosignatures were quantified in each osteosarcoma sample using ssGSEA (Barbie et al., 2009). The ssGSEA enrichment score was obtained by “GSVA” and “GSEABase” packages in R software (version 4.0.0). Gene signatures for each immune cell type and immune pathway were obtained from previously published data, representing the overall immune activity of the tumor microenvironment (He et al., 2018).

Clustering of Osteosarcoma Samples

Then, we carried out a hierarchical grouping of osteosarcoma based on the 29 immune biosignature enrichment levels (ssGSEA scores). ESTIMATE was carried out to explore the immune score, as well as the tumor purity for each osteosarcoma sample.

Principal Component Analysis

Principal component analysis (PCA) was performed to compare the transcriptional profiles between the different immune subtypes. The gene names with corresponding expression value and sample information were loaded, and the analysis was conducted by “limma” package using princomp function and visualized by “ggplot2” package in R software (version 4.0.0).

Identification of Differentially Expressed Immune-Related Genes

The differential expression of immune-related genes was identified using the “limma” R package. Adj. p < 0.05, and | fold change| > 2 were defined as the threshold values. The differentially expressed immune-related genes were visualized by the “ggplot2” package in R software (version 4.0.0).

Multiplex Immunohistochemistry

Tissue microarray of osteosarcoma specimens (L714901) was obtained from Biotech Company (Xian, China) and used to explore further the association of TYROBP expression with CD8+ T cell in the osteosarcoma tumor microenvironment.

Multiplex immunohistochemistry was conducted using the sequential staining cycles as follows. In brief, formalin-fixed, paraffin-embedded osteosarcoma tissue sections were deparaffinized and then underwent microwave treatment in citrate for antigen retrieval. Then, they were blocked with 10% normal goat serum and incubated overnight with primary antibodies: mouse anti-CD8 antibody (1:100, ab17147, Abcam) and rabbit anti-TYROBP antibody (1:200, ab124834, Abcam). Next, sections were incubated with the corresponding horseradish peroxidase-conjugated second antibodies (Abcam, CN) for 30 min at room temperature. The antigenic binding sites were visualized using the tyramide signal amplification dye. Fluorescein isothiocyanate-tyramide (1:1,000, G1235, Servicebio) and Cy3-tyramide (1:1,000, G1235, Servicebio) were applied to each antibody.

Constitution of a Risk Model

Critical immune-related genes (IRGs) that presented protein–protein interaction (PPI) network were utilized in Lasso Cox regression to calculate the coefficients; the risk-score formula was formulated as:

Risk score = i =  1 N ( E x p i × C o e i )

where N = 2, Expi was the expression value of every five hypoxia genes, whereas the Coei denoted the corresponding multivariable Cox regression coefficient.

Survival Analyses

Kaplan–Meier curves were used to visualize and compare the OS of osteosarcoma patients in different immunogenomic subtypes via survival and survminer R package. A receiver operating characteristic curve was constructed to verify our risk model accuracy via the survivalROC R package.

Visualization of Protein–Protein Interaction Network

We constricted the PPI network using the STRING database Cytoscape software. Cytoscape1 constitutes a publicly accessible software web source used in visualizing complex networks, as well as for integrating these with any kind of attribute data. The Cytoscape web portal was used to build a protein cross talk association network and investigated the interaction association of the core IRGs.


Immunogenomic Profiling Identifies Two Osteosarcoma Subtypes

We assessed 29 immune-linked gene sets that represented multiple immune cell types, cascades, and functions. The enrichment levels or the activity of immune cells and signaling cascades were quantified in the osteosarcoma samples using ssGSEA analysis. Totally, we hierarchically categorized osteosarcoma in two cohorts (TARGET and GSE21257). Noteworthy, all cohorts presented similar clustering findings, with two clusters clearly being identified. We termed the classes as immunity high (Immunity_H) and immunity low (Immunity_L). Overall, the Immunity_H subtype was correlated with high immune cell infiltration and activated immune pathways indicating an immune-hot phenotype, whereas the Immunity_L subtype was associated with low immune cell infiltration indicating an immune-cold phenotype (Figures 1A,B). Moreover, the immune scores were remarkably elevated in Immunity_H than that in Immunity_L in TARGET and GSE21257 datasets (Figure 1C). Besides, opposite trends were reported with regards to tumor purity between the two subtypes, with tumor purity escalating from Immunity_L to Immunity_H (Figure 1D).


Figure 1. Hierarchical clustering of osteosarcoma identified two immune subtypes in TARGET (A) and GSE21257 datasets (B). Correlation between immunogenomic subtypes and immune cell invasion level in immune score (C) and tumor purity (D). Violin plots particularly indicated differences in two subtypes. ****p < 0.0001.

Altogether, these data indicate that Immunity_H has the greatest number of immune cells, which may benefit from immune checkpoint blockade therapy.

Immunity_H Subtype Is Associated With Immune-Related Signaling Pathways, a Favorable Prognosis and a Potential Benefit to Immunotherapy

Principal component analysis was then performed to compare the transcriptional profiles between the Immunity_H and Immunity_L subtypes, which displayed a clear distinction. In detail, PCA showed that the samples from the two clusters were well separated from each other (Figure 2A). To further identify relevant signaling cascades activated in the Immunity_H and the Immunity_L groups, we performed a GSEA analysis. Gene sets were differentially enriched in the Immunity_H groups of the TARGET web portal, as they were associated with processes that trigger immunity, such as dendritic cell antigen processing and presentation, interleukin-10 secretion, interleukin-12 secretion, and natural killer cell chemotaxis (Figure 2B). On the contrary, the cancer-linked signaling cascades were hyperactivated in the Immunity_L group, including the positive modulation of G1/S transition of the mitotic cell cycle and vascular endothelial cell proliferation. A similar trend was observed in the GES21257 dataset (Figure 2C). This finding verified the increased immune activity in Immunity_H.


Figure 2. (A) Principal component analysis plots. Red dots represent Immunity_H samples, and blue dots represent Immunity_L samples; (B,C) GSEA identified differential immune-related signaling cascades enriched in Immunity_H, as well as Immunity_L subtypes in TAGRT and GSE21257 datasets. Size of dots represented -log10 (p-value); (D) Kaplan–Meier curves of OS in Immunity_H and Immunity_L subtypes; (E) response rate to immunotherapy predicted by TIDE (a computational framework developed to evaluate potential of tumor immune escape from gene expression profiles of cancer samples:; (F) half-maximal inhibitory concentration scores of gemcitabine and mitomycin C in immunity_H and immunity_L subtypes. *p < 0.05, ***p < 0.001, and ****p < 0.0001.

Previous studies showed that osteosarcoma with elevated immune activity and high immune cell infiltration was associated with more favorable clinical outcomes. In line with these pieces of evidence, survival analyses showed that these immune subtypes had distinct clinical outcomes. The Immunity_H subtype had a better survival prognosis than the Immunity_L subtypes (Figure 2D).

We then used TIDE (a computational framework developed to evaluate the potential of tumor immune escape from the gene expression profiles of cancer samples:2) to evaluate the potential clinical efficacy of immunotherapy in these immune subtypes. In our results, the Immunity_H subtype had a higher response rate than the Immunity_L subgroup, implying that Immunity_H patients could benefit more from immune checkpoint inhibitor (ICI) therapy (Figure 2E). Besides, chemosensitivity of these subtypes was assessed by the Genomics of Drug Sensitivity in Cancer database, and the results showed that half-maximal inhibitory concentration scores of gemcitabine and mitomycin C in immunity_H were significantly lower than immunity_L subtype, indicating sensitivity to chemotherapy (Figure 2F).

Differences in Human Leukocyte Antigen Genes and Immune Checkpoint Expression Related to Immune Phenotypes

Tumor-related antigen presentation via MHC class I complexes is a prerequisite for immune surveillance and is instrumental for the clinical response of immunotherapies targeting immune checkpoints. Hence, we assessed transcript levels of 24 human leukocyte antigen (HLA) genes in TARGET and GSE21257 validation cohorts. Most immune HLA genes exhibited markedly lower expression levels in Immunity_L, indicating impaired antigen presentation on tumor cells as an escape mechanism from immune surveillance (Figures 3A,B). Expression of immune checkpoints has been exploited as predictive biomarkers for checkpoint inhibitor-based immunotherapy; we further compared the differences in immune checkpoints (PD-L1, PD-1, CTLA4, TIM3, LAG3, and TIGIT) related to immunity high and immunity low subtypes. Interestingly, high immune checkpoint expression was observed in the Immunity_H subtype, indicating a potential benefit from checkpoint inhibitor-based immunotherapy (Figures 3C,D).


Figure 3. Box plots presented differential expression of HLA genes (A,B) and multiple immune checkpoints (C,D) between Immune_H and Immune_L subtypes. **p < 0.01, ***p < 0.001, and ****p < 0.0001.

Identification of Key Genes Between Immunity_H and Immunity_L Subtypes

We further identified the differentially expressed genes (DEGs) between Immunity_H and Immunity_L groups and revealed the critical immune-linked genes. A total of 303 DEGs in the TARGET cohort and 173 genes in the GES21257 cohort were identified with an intersection of 83 genes in the two datasets, and most of the genes were upregulated in the Immunity_H group (Figures 4A,B). PPI network was then performed to reveal the key genes associated with the immune subtypes. The top 10 immune-related genes with the highest degrees of interaction were determined, including TYROBP, ITGB2, LCP2, C1QB, C1QC, CD74, HLA-DRA, CXCL10, CCL5, and CXCL9, suggesting their pivotal role in modulating tumor immunity (Figure 4C).


Figure 4. (A) Volcano plot presented distribution of differentially expressed genes (DEGs) quantified among Immune_H and Immune_L subtypes with threshold of | log2 Fold change| > 1 and P < 0.05 in TARGET and GSE21257 cohort; (B) Venn plot showed intersection of differentially expressed genes among Immune_H and Immune_L subtypes; (C) PPI network analysis showed interactions among 83 DEGs.

TYROBP Expression Is Associated With CD8+ T Cell Infiltration

According to PPI analysis, TYROBP showed the highest degrees of interaction with other DEGs, suggesting its important role in driving the immunophenotype. We further investigated the association of TYROBP expression with immune cells in the osteosarcoma tumor microenvironment. Using the CIBERSORT method in combination with the LM22 signature matrix, we calculated differences in the immune infiltration of 22 immune cell types between low- and high-TYROBP groups. Figure 5A summarized the 22 immune cell landscapes among osteosarcoma patients in the TARGET dataset. The results showed that patients with high TYROBP expression had significantly higher proportions of immune cells (e.g., CD8 T cell, CD4 naive T cell, CD4 memory activated T cell, macrophages M0, macrophages M1, and macrophages M2) (Figure 5B). Among the immune cells, CD8+ T cells, also known as cytotoxic T cells, play a vital role in antitumor immunity by killing cancer cells. Our results demonstrated that TYROBP expression showed a positive correlation with CD8+ cell infiltration (Figure 5C), and a high TYROBP expression group was infiltrated with a high CD8+ T cell level (Figure 5D). To further validate the association between TYROBP and CD8+ T cells in osteosarcoma, we performed multiplex immunohistochemistry. In line with results from bioinformatics analysis, multiplexed immunohistochemistry showed that high TYROBP expression was associated with high CD8+ T cell levels in the tumor microenvironment (Figure 5E).


Figure 5. (A) Relative proportion of immune infiltration in Immune_H and Immune_L groups; (B) box plots visualizing significantly different immune cells between Immune_H and Immune_L groups; (C) correlation of TYROBP and CD8+ T cell infiltration; (D) box plot showed CD8+ T cell infiltration level in TYROBP high and low groups; (E) multiplex immunohistochemistry validated differentially infiltrated CD8+ T cell in TYROBP high and low groups. **p < 0.01, ***p < 0.001, and ****p < 0.0001.

Construction and Confirmation of the Immune-Related Prognostic Biosignature

To improve the accuracy and reliability of the predictive model, we treated TARGET as the training set and the GES21257 as validation sets. Lasso Cox regression analysis was subsequently carried out to construct an immune-associated predictive model using 34 critical IRGs presented in the PPI network (Figures 6A,B). A risk-score formula was established as follows: risk score = −0.267 × FPR1–0.001 × FCER1G. As shown in Figures 6A,B, the heatmap showed that the expression level of FPR1 and FCER1G was decreased, accompanied by the higher risk scores. Additionally, we examined the correlation between risk score and survival status (Figure 6C). Our data showed that the number of alive statuses in the low-risk group was markedly higher relative to the high-risk group (Figure 6D).


Figure 6. (A,B) Lasso Cox analysis uncovered two genes most associated with OS in TARGET dataset. (C,D) Kaplan–Meier curves of OS for LGG patients based on IPS in TARGET cohort and GSE21257 cohort. (E,F) Risk scores distribution, survival status of each patient, and heatmaps of prognostic four-gene biosignature in TARGET and GSE21257 datasets.

To explore the prognostic significance of the risk model in osteosarcoma, we performed a Kaplan–Meier analysis. As indicated in Figure 6, the high-risk score was correlated with dismal OS in the TARGET dataset (cutoff values = −0.44), which was further confirmed by the GSE21257 dataset (cutoff values = −2.4) (Figure 6E).


In this study, we focused on uncovering immune-associated osteosarcoma subtypes through immunogenomic profiling. Our study demonstrates that osteosarcoma could be categorized into immunity low and immunity high subtypes. We established that this categorization was reproducible, as well as predictable. The immunity high osteosarcoma subtype was enriched with high immune cell infiltration and immunoactive signaling pathway indicating an immune hot phenotype, whereas immunity low was associated with low immune cell infiltration and immunosuppressive signaling pathway indicating an immune cold phenotype. Moreover, we developed and validated an immune-related prognostic signature.

Immunity_H had a higher level of immune cell invasion, as well as antitumor immune activities, such as high levels of immunoactivated cell infiltration, and activated immune signaling pathway. This analysis validated the increased antitumor immune activity in the Immunity_H subtype. Indeed, accumulating evidence has documented that the density and level of tumor-infiltrating lymphocytes are positively correlated with survival prognosis in multiple cancers (Hendry et al., 2017; Jiang et al., 2017; Wang et al., 2020; Yao et al., 2017).

The HLA genes code for MHC I, as well as MHC II molecules, which deliver pathogen-originated short peptides to T cells and trigger an adaptive immune response (Apanius et al., 2017; Danchin et al., 2004). Previous researches have reported that the expression of HLA may regulate ICI response in metastatic melanoma (Chowell et al., 2018; Garrido et al., 2016; O’Donnell et al., 2019). The study demonstrated that elevated tumor-distinct expression of MHC-I was pivotal for the response to therapy with anti-CTLA-4. In the meantime, tumor-specific MHC-II expression was critical for the response to the therapy with anti-PD-1 (Johnson et al., 2016). Additionally, recent studies have identified predictive biomarkers (e.g., PD-L1, PD-1, and CTLA-4) that can be exploited to predict response to ICI therapy (Gibney et al., 2016; Lin et al., 2021). In our study, in the Immunity_H subtype, the HLA genes and immune checkpoint expression levels were also the highest compared with the Immunity_L subtype. In contrast, we revealed that Immunity_High exhibited strong immune activity given the high immune invasion level, percentage of distinct tumor-infiltrating lymphocytes, HLA richness, and level of expression and exhibited relatively improved clinical outcomes.

We also identified several key genes with the highest interaction degrees between Immune_H and Immunity_L subtypes, including TYROBP, ITGB2, LCP2, C1QB, C1QC, CD74, HLA-DRA, CXCL10, CCL5, and CXCL9. Previous studies have demonstrated that TYROBP, also known as DAP12, is overexpressed in various cancers. Functionally, TYROBP encodes a transmembrane signaling polypeptide on the surface of a variety of immune cells that contains an immunoreceptor tyrosine-based activation motif in its cytoplasmic domain and mediates signaling transductions (Dietrich et al., 2000; Shabo et al., 2013; Töpfer et al., 2015). It has been reported that TYROBP expression is associated with CD8 T cell infiltration in gastric cancer and clear cell renal cell carcinoma (Jiang et al., 2020; Wu et al., 2020). In line with these pieces of evidence, our result showed that TYROBP expression is positively correlated with CD8 T cell infiltration in the osteosarcoma tumor microenvironment. CD74 plays an important function in the processing of MHC class II antigens through the stabilization of peptide-free class II alpha/beta heterodimers in a complex shortly following their synthesis and directs the transportation of the complex to the endosomal/lysosomal system from the endoplasmic reticulum (Badve et al., 2002; Denzin et al., 1996). High CD74 expression in Immune_H indicated a high antigen presentation and is important for the clinical response of immunotherapies. Pro-inflammatory cytokine participates in a wide variety of processes such as differentiation, chemotaxis, activation of peripheral immune cells, and regulation of immune cell growth (Morrell et al., 2014). CXCL10 is an important pro-inflammatory cytokine and has been reported to mediate the mobilization of tumor-inhibitory CXCR3+ T cells, as well as natural killer cells, into solid cancers (Tokunaga et al., 2018). Respectively, a high CXCL10 concentration is correlated with a higher immune cell infiltration and better survival in several malignancies (Bronger et al., 2016).

In our research, we established a novel classification of osteosarcoma based on immunogenomic profiling. We also construct and validate the immune-associated prognostic signature. This classification of osteosarcoma based on immunogenomic profiling may provide valuable information for immunotherapy strategies in osteosarcoma patients. Besides, we also identified TYROBP as a key immunoregulatory gene. Nevertheless, it should be noted that our findings require further validation in vitro or in vivo and a larger sample cohort. Our findings should be interpreted with this limitation in mind.

In conclusion, we established a novel osteosarcoma classification based on the differences in the transcriptome of 29 immune biosignatures in tumor samples. This classification had significant clinical outcomes for estimating the prognosis, as well as the responsiveness of immunotherapy.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Review Committee of the Third Affiliated Hospital of Southern Medical University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

XeW, LW, and WX: conceptualization, methodology, and writing. XuW: software. DK: validation. JL: data curation. WL and XB: project administration and funding acquisition. All authors have read and agreed to the published version of the manuscript.


This work was supported by Jiangmen Medical and Health Technology Plan Project (Project No.2019B003).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


  1. ^
  2. ^


Apanius, V., Penn, D., Slev, P. R., Ruff, L. R., and Potts, W. K. (2017). The nature of selection on the major histocompatibility complex. Crit. Rev. Immunol. 37, 75–120. doi: 10.1615/CritRevImmunol.v37.i2-6.10

PubMed Abstract | CrossRef Full Text | Google Scholar

Badve, S., Deshpande, C., Hua, Z., and Lögdberg, L. (2002). Expression of invariant chain (CD 74) and major histocompatibility complex (MHC) class II antigens in the human fetus. J. Histochem. Cytochem. 50, 473–482. doi: 10.1177/002215540205000404

PubMed Abstract | CrossRef Full Text | Google Scholar

Barbie, D. A., Tamayo, P., Boehm, J. S., Kim, S. Y., Moody, S. E., Dunn, I. F., et al. (2009). Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature 462, 108–112. doi: 10.1038/nature08460

PubMed Abstract | CrossRef Full Text | Google Scholar

Baxevanis, C. N., Perez, S. A., and Papamichail, M. (2009). Cancer immunotherapy. Crit. Rev. Clin. Lab. Sci. 46, 167–189. doi: 10.1080/10408360902937809

PubMed Abstract | CrossRef Full Text | Google Scholar

Bronger, H., Singer, J., Windmüller, C., Reuning, U., Zech, D., Delbridge, C., et al. (2016). CXCL9 and CXCL10 predict survival and are regulated by cyclooxygenase inhibition in advanced serous ovarian cancer. Br. J. Cancer 115, 553–563. doi: 10.1038/bjc.2016.172

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C., Xie, L., Ren, T., Huang, Y., Xu, J., and Guo, W. (2021). Immunotherapy for osteosarcoma: fundamental mechanism, rationale, and recent breakthroughs. Cancer Lett. 500, 1–10. doi: 10.1016/j.canlet.2020.12.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Chowell, D., Morris, L. G. T., Grigg, C. M., Weber, J. K., Samstein, R. M., Makarov, V., et al. (2018). Patient HLA class I genotype influences cancer response to checkpoint blockade immunotherapy. Science 359, 582–587. doi: 10.1126/science.aao4572

PubMed Abstract | CrossRef Full Text | Google Scholar

Danchin, E., Vitiello, V., Vienne, A., Richard, O., Gouret, P., McDermott, M. F., et al. (2004). The major histocompatibility complex origin. Immunol. Rev. 198, 216–232. doi: 10.1111/j.0105-2896.2004.00132.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Denzin, L. K., Hammond, C., and Cresswell, P. (1996). HLA-DM interactions with intermediates in HLA-DR maturation and a role for HLA-DM in stabilizing empty HLA-DR molecules. J. Exp. Med. 184, 2153–2165. doi: 10.1084/jem.184.6.2153

PubMed Abstract | CrossRef Full Text | Google Scholar

Dietrich, J., Cella, M., Seiffert, M., Bühring, H. J., and Colonna, M. (2000). Cutting edge: signal-regulatory protein beta 1 is a DAP12-associated activating receptor expressed in myeloid cells. J. Immunol. 164, 9–12. doi: 10.4049/jimmunol.164.1.9

PubMed Abstract | CrossRef Full Text | Google Scholar

Durfee, R. A., Mohammed, M., and Luu, H. H. (2016). Review of osteosarcoma and current management. Rheumatol. Ther. 3, 221–243. doi: 10.1007/s40744-016-0046-y

PubMed Abstract | CrossRef Full Text | Google Scholar

ElKordy, M. A., ElBaradie, T. S., ElSebai, H. I., KhairAlla, S. M., and Amin, A. A. E. (2018). Osteosarcoma of the jaw: challenges in the diagnosis and treatment. J. Egypt Natl. Canc. Inst. 30, 7–11. doi: 10.1016/j.jnci.2018.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Emens, L. A., Ascierto, P. A., Darcy, P. K., Demaria, S., Eggermont, A. M. M., Redmond, W. L., et al. (2017). Cancer immunotherapy: opportunities and challenges in the rapidly evolving clinical landscape. Eur. J. Cancer 81, 116–129. doi: 10.1016/j.ejca.2017.01.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Garrido, F., Aptsiauri, N., Doorduijn, E. M., Garcia Lora, A. M., and van Hall, T. (2016). The urgent need to recover MHC class I in cancers for effective immunotherapy. Curr. Opin. Immunol. 39, 44–51. doi: 10.1016/j.coi.2015.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Gibney, G. T., Weiner, L. M., and Atkins, M. B. (2016). Predictive biomarkers for checkpoint inhibitor-based immunotherapy. Lancet Oncol. 17, e542–e551. doi: 10.1016/s1470-2045(16)30406-5

CrossRef Full Text | Google Scholar

He, Y., Jiang, Z., Chen, C., and Wang, X. (2018). Classification of triple-negative breast cancers based on immunogenomic profiling. J. Exp. Clin. Cancer Res. 37:327. doi: 10.1186/s13046-018-1002-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Hendry, S., Salgado, R., Gevaert, T., Russell, P. A., John, T., Thapa, B., et al. (2017). Assessing tumor-infiltrating lymphocytes in solid tumors: a practical review for pathologists and proposal for a standardized method from the international immuno-oncology biomarkers working group: part 2: TILs in melanoma, gastrointestinal tract carcinomas, non-small cell lung carcinoma and mesothelioma, endometrial and ovarian carcinomas, squamous cell carcinoma of the head and neck, genitourinary carcinomas, and primary brain tumors. Adv. Anat. Pathol. 24, 311–335. doi: 10.1097/pap.0000000000000161

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, D., Liu, Y., Wang, H., Wang, H., Song, Q., Sujie, A., et al. (2017). Tumour infiltrating lymphocytes correlate with improved survival in patients with esophageal squamous cell carcinoma. Sci. Rep. 7:44823. doi: 10.1038/srep44823

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, J., Ding, Y., Wu, M., Lyu, X., Wang, H., Chen, Y., et al. (2020). Identification of TYROBP and C1QB as two novel key genes with prognostic value in gastric cancer by network analysis. Front. Oncol. 10:1765. doi: 10.3389/fonc.2020.01765

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, D. B., Estrada, M. V., Salgado, R., Sanchez, V., Doxie, D. B., Opalenik, S. R., et al. (2016). Melanoma-specific MHC-II expression represents a tumour-autonomous phenotype and predicts response to anti-PD-1/PD-L1 therapy. Nat. Commun. 7:10582. doi: 10.1038/ncomms10582

PubMed Abstract | CrossRef Full Text | Google Scholar

Kansara, M., Teng, M. W., Smyth, M. J., and Thomas, D. M. (2014). Translational biology of osteosarcoma. Nat. Rev. Cancer 14, 722–735. doi: 10.1038/nrc3838

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, W., Qiu, X., Sun, P., Ye, Y., Huang, Q., Kong, L., et al. (2021). Association of IDH mutation and 1p19q co-deletion with tumor immune microenvironment in lower-grade glioma. Mol. Ther. Oncolytics 21, 288–302. doi: 10.1016/j.omto.2021.04.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, W., Xu, Y., Chen, X., Liu, J., Weng, Y., Zhuang, Q., et al. (2020). Radiation-induced small extracellular vesicles as “carriages” promote tumor antigen release and trigger antitumor immunity. Theranostics 10, 4871–4884. doi: 10.7150/thno.43539

PubMed Abstract | CrossRef Full Text | Google Scholar

Morrell, C. N., Aggrey, A. A., Chapman, L. M., and Modjeski, K. L. (2014). Emerging roles for platelets as immune and inflammatory cells. Blood 123, 2759–2767. doi: 10.1182/blood-2013-11-462432

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Donnell, J. S., Teng, M. W. L., and Smyth, M. J. (2019). Cancer immunoediting and resistance to T cell-based immunotherapy. Nat. Rev. Clin. Oncol. 16, 151–167. doi: 10.1038/s41571-018-0142-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Shabo, I., Olsson, H., Stål, O., and Svanvik, J. (2013). Breast cancer expression of DAP12 is associated with skeletal and liver metastases and poor survival. Clin. Breast Cancer 13, 371–377. doi: 10.1016/j.clbc.2013.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, J. K., Cote, G. M., Choy, E., Yang, P., Harmon, D., Schwab, J., et al. (2014). Programmed cell death ligand 1 expression in osteosarcoma. Cancer Immunol. Res. 2, 690–698. doi: 10.1158/2326-6066.Cir-13-0224

PubMed Abstract | CrossRef Full Text | Google Scholar

Simpson, E., and Brown, H. L. (2018). Understanding osteosarcomas. JAAPA 31, 15–19. doi: 10.1097/01.JAA.0000541477.24116.8d

CrossRef Full Text | Google Scholar

Tokunaga, R., Zhang, W., Naseem, M., Puccini, A., Berger, M. D., Soni, S., et al. (2018). CXCL9, CXCL10, CXCL11/CXCR3 axis for immune activation–a target for novel cancer therapy. Cancer Treat. Rev. 63, 40–47. doi: 10.1016/j.ctrv.2017.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Töpfer, K., Cartellieri, M., Michen, S., Wiedemuth, R., Müller, N., Lindemann, D., et al. (2015). DAP12-based activating chimeric antigen receptor for NK cell tumor immunotherapy. J. Immunol. 194, 3201–3212. doi: 10.4049/jimmunol.1400330

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Tian, S., Sun, J., Zhang, J., Lin, L., and Hu, C. (2020). The presence of tumour-infiltrating lymphocytes (TILs) and the ratios between different subsets serve as prognostic factors in advanced hypopharyngeal squamous cell carcinoma. BMC Cancer 20:731. doi: 10.1186/s12885-020-07234-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Li, B., Ren, Y., and Ye, Z. (2016). T-cell-based immunotherapy for osteosarcoma: challenges and opportunities. Front. Immunol. 7:353. doi: 10.3389/fimmu.2016.00353

PubMed Abstract | CrossRef Full Text | Google Scholar

Whelan, J. S., and Davis, L. E. (2018). Osteosarcoma, chondrosarcoma, and chordoma. J. Clin. Oncol. 36, 188–193. doi: 10.1200/jco.2017.75.1743

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, P., Xiang, T., Wang, J., Lv, R., and Wu, G. (2020). TYROBP is a potential prognostic biomarker of clear cell renal cell carcinoma. FEBS Open Bio. 10, 2588–2604. doi: 10.1002/2211-5463.12993

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y. (2015). Cancer immunotherapy: harnessing the immune system to battle cancer. J. Clin. Invest. 125, 3335–3337. doi: 10.1172/jci83871

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, W., He, J. C., Yang, Y., Wang, J. M., Qian, Y. W., Yang, T., et al. (2017). The prognostic value of tumor-infiltrating lymphocytes in hepatocellular carcinoma: a systematic review and meta-analysis. Sci. Rep. 7:7525. doi: 10.1038/s41598-017-08128-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: osteosarcoma, immune subtype, immune checkpoint inhibitors, TYROBP, prognosis

Citation: Wang X, Wang L, Xu W, Wang X, Ke D, Lin J, Lin W and Bai X (2021) Classification of Osteosarcoma Based on Immunogenomic Profiling. Front. Cell Dev. Biol. 9:696878. doi: 10.3389/fcell.2021.696878

Received: 18 April 2021; Accepted: 14 June 2021;
Published: 16 July 2021.

Edited by:

Andrea Del Fattore, Bambino Gesù Children Hospital (IRCCS), Italy

Reviewed by:

Zhan Wang, Zhejiang University, China
Yongbiao Huang, Huazhong University of Science and Technology, China

Copyright © 2021 Wang, Wang, Xu, Wang, Ke, Lin, Lin and Bai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Wanzun Lin,; Xiaochun Bai,

These authors have contributed equally to this work