Exploration and Validation of a Novel Inflammatory Response-Associated Gene Signature to Predict Osteosarcoma Prognosis and Immune Infiltration

Prognostic and predictive value of super-enhancer-derived signatures for survival and lung metastasis in osteosarcoma

BACKGROUND

Risk stratification and personalized care are crucial in managing osteosarcoma due to its complexity and heterogeneity. However, current prognostic prediction using clinical variables has limited accuracy. Thus, this study aimed to explore potential molecular biomarkers to improve prognostic assessment.

METHODS

High-throughput inhibitor screening of 150 compounds with broad targeting properties was performed and indicated a direction towards super-enhancers (SEs). Bulk RNA-seq, scRNA-seq, and immunohistochemistry (IHC) were used to investigate SE-associated gene expression profiles in osteosarcoma cells and patient tissue specimens. Data of 212 osteosarcoma patients who received standard treatment were collected and randomized into training and validation groups for retrospective analysis. Prognostic signatures and nomograms for overall survival (OS) and lung metastasis-free survival (LMFS) were developed using Cox regression analyses. The discriminatory power, calibration, and clinical value of nomograms were evaluated.

RESULTS

High-throughput inhibitor screening showed that SEs significantly contribute to the oncogenic transcriptional output in osteosarcoma. Based on this finding, focus was given to 10 SE-associated genes with distinct characteristics and potential oncogenic function. With multi-omics approaches, the hyperexpression of these genes was observed in tumor cell subclusters of patient specimens, which were consistently correlated with poor outcomes and rapid metastasis, and the majority of these identified SE-associated genes were confirmed as independent risk factors for poor outcomes. Two molecular signatures were then developed to predict survival and occurrence of lung metastasis: the SE-derived OS-signature (comprising LACTB, CEP55, SRSF3, TCF7L2, and FOXP1) and the SE-derived LMFS-signature (comprising SRSF3, TCF7L2, FOXP1, and APOLD1). Both signatures significantly improved prognostic accuracy beyond conventional clinical factors.

CONCLUSIONS

Oncogenic transcription driven by SEs exhibit strong associations with osteosarcoma outcomes. The SE-derived signatures developed in this study hold promise as prognostic biomarkers for predicting OS and LMFS in patients undergoing standard treatments. Integrative prognostic models that combine conventional clinical factors with these SE-derived signatures demonstrate substantially improved accuracy, and have the potential to facilitate patient counseling and individualized management.

New gene signature from the dominant infiltration immune cell type in osteosarcoma predicts overall survival

The immune microenvironment of osteosarcoma (OS) has been reported to play an important role in disease progression and prognosis. However, owing to tumor heterogeneity, it is not ideal to predict OS prognosis by examining only infiltrating immune cells. This work aimed to build a prognostic gene signature based on similarities in the immune microenvironments of OS patients. Public datasets were used to examine the correlated genes, and the most consistent dominant infiltrating immune cell type was identified. The LASSO Cox regression model was used to establish a multiple-gene risk prediction signature. A nine-gene prognostic signature was generated from the correlated genes for M0 macrophages and then proven to be effective and reliable in validation cohorts. Signature comparison indicated the priority of the signature. Multivariate Cox regression models indicated that the signature risk score is an independent prognostic factor for OS patients regardless of the Huvos grade in all datasets. In addition, the results of the association between the signature risk score and chemotherapy sensitivity also showed that there was no significant difference in the sensitivity of any drugs between the low- and high-risk groups. A GSEA of GO and KEGG pathways found that antigen processing- and presentation-related biological functions and olfactory transduction receptor signaling pathways have important roles in signature functioning. Our findings showed that M0 macrophages were the dominant infiltrating immune cell type in OS and that the new gene signature is a promising prognostic model for OS patients.

A tumor microenvironment-based prognostic index for osteosarcoma

BACKGROUND

The tumor microenvironment (TME) has a central role in the oncogenesis of osteosarcomas. The composition of the TME is essential for the interaction between tumor and immune cells. The aim of this study was to establish a prognostic index (TMEindex) for osteosarcoma based on the TME, from which estimates about patient survival and individual response to immune checkpoint inhibitor (ICI) therapy can be deduced.

METHODS

Based on osteosarcoma samples from the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database, the ESTIMATE algorithm was used to estimate ImmuneScore and StromalScore. Combined differentially expressed gene analysis, weighted gene co-expression network analyses, the Least Absolute Shrinkage and Selection Operator regression and stepwise regression to construct the TMEindex. The prognostic role of TMEindex was validated in three independent datasets. The molecular and immune characteristics of TMEindex and the impact on immunotherapy were then comprehensively investigated. The expression of TMEindex genes in different cell types and its effects on osteosarcoma cells were explored by scRNA-Seq analysis and molecular biology experiments.

RESULTS

Fundamental is the expression of MYC, P4HA1, RAMP1 and TAC4. Patients with high TMEindex had worse overall survival, recurrence-free survival, and metastasis-free survival. TMEindex is an independent prognostic factor in osteosarcoma. TMEindex genes were mainly expressed in malignant cells. The knockdown of MYC and P4HA1 significantly inhibited the proliferation, invasion and migration of osteosarcoma cells. A high TME index is related to the MYC, mTOR, and DNA replication-related pathways. In contrast, a low TME index is related to immune-related signaling pathways such as the inflammatory response. The TMEindex was negatively correlated with ImmuneScore, StromalScore, immune cell infiltration, and various immune-related signature scores. Patients with a higher TMEindex had an immune-cold TME and higher invasiveness. Patients with a low TME index were more likely to respond to ICI therapy and achieve clinical benefit. In addition, the TME index correlated with response to 29 oncologic drugs.

CONCLUSIONS

The TMEindex is a promising biomarker to predict the prognosis of patients with osteosarcoma and their response to ICI therapy, and to distinguish the molecular and immune characteristics.

Strategies to Overcome Resistance to Immune-Based Therapies in Osteosarcoma

Improving the prognosis and cure rate of HGOSs (high-grade osteosarcomas) is an absolute need. Immune-based treatment approaches have been increasingly taken into consideration, in particular for metastatic, relapsed and refractory HGOS patients, to ameliorate the clinical results currently achieved. This review is intended to give an overview on the immunotherapeutic treatments targeting, counteracting or exploiting the different immune cell compartments that are present in the HGOS tumor microenvironment. The principle at the basis of these strategies and the possible mechanisms that HGOS cells may use to escape these treatments are presented and discussed. Finally, a list of the currently ongoing immune-based trials in HGOS is provided, together with the results that have been obtained in recently completed clinical studies. The different strategies that are presently under investigation, which are generally aimed at abrogating the immune evasion of HGOS cells, will hopefully help to indicate new treatment protocols, leading to an improvement in the prognosis of patients with this tumor.

Osteosarcoma subtypes based on platelet-related genes and tumor microenvironment characteristics

Background

Osteosarcoma is a common metastatic tumor in children and adolescents. Because of its easy metastasis, patients often show a poor prognosis. Recently, researchers have found that platelets are closely related to metastasis of a variety of malignant tumors, but the role of platelets related characteristics in osteosarcoma is still unknown. The purpose of this study is to explore the characteristics of platelet-related subtypes and cell infiltration in tumor microenvironment.

Methods

We collected osteosarcoma cohorts from TCGA and GEO databases, and explored the molecular subtypes mediated by platelet-related genes and the related TME cell infiltration according to the expression of platelet-related genes in osteosarcoma. In addition, we also explored the differentially expressed genes (DEGs) among different molecular subtypes and established a protein-protein interaction network (PPI). Then we constructed a platelet scoring model by Univariate cox regression and least absolute shrinkage and selection operator (Lasso) cox regression model to quantify the characteristics of platelet in a single tumor. RT-PCR was used to investigate the expression of six candidate genes in osteosarcoma cell lines and normal osteoblast lines. Finally, we also predicted potential drugs with therapeutic effects on platelet-related subtypes.

Results

We found that platelet-related genes (PRGs) can distinguish osteosarcoma into two different platelet-related subtypes, C1 and C2. And the prognosis of the C2 subtype was significantly worse than that of C1 subtype. The results of ESTIMATE analysis and GO/KEGG enrichment showed that the differences between different subtypes were mainly concentrated in immune response pathways, and the immune response of C2 was inhibited relative to C1. We further studied the relationship between platelet-related subtypes and immune cell infiltration. We found that the distribution of most immune cells in C1 subtype was higher than that in C2 subtype, and there was a correlation between C1 subtype and more immune cells. Finally, we screened the PRGs related to the prognosis of osteosarcoma through Univariate Cox regression, established independent prognostic platelet characteristics consisting of six genes to predict the prognosis of patients with OS, and predicted the drugs that may be used in the treatment of osteosarcoma. RT-PCR was used to verify the expression of candidate genes in osteosarcoma cells.

Conclusion

Platelet scoring model is a significant biomarker, which is of great significance to determine the prognosis, molecular subtypes, characteristics of TME cell infiltration and therapy in patients with OS.

Identification of the Transcription Co-Factor–Related Gene Signature and Risk Score Model for Osteosarcoma

Osteosarcoma is a malignant tumor with a poor prognosis. Nowadays, there is a lack of good methods to assess the prognosis of osteosarcoma patients. Transcription co-factors (TcoFs) play crucial roles in transcriptional regulation through the interaction with transcription factors (TFs). Many studies have revealed that TcoFs are related to many diseases, especially cancer. However, few studies have been reported about prognostic prediction models of osteosarcoma by using TcoF-related genes. In order to construct a prognostic risk model with TcoF-related genes, the mRNA expression data and matched clinical information of osteosarcoma were downloaded from the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database and the Gene Expression Omnibus (GEO) database. TARGET was used as a training set and GSE21257 from GEO was used as a validation set. Univariate Cox regression was performed to select 13 TcoF-related candidate genes, of which five genes (LMO2, MAML3, MTF2, RBPMS, and SIRT1) were finally used to construct the prognostic risk model by LASSO Cox regression analysis. The Kaplan–Meier (K-M) survival curves showed an obvious difference between high- and low-risk groups. The receiver operating characteristic (ROC) curves based on TARGET demonstrated that this risk model was credible (1-year AUC: 0.607; 3-years AUC: 0.713; 5-years AUC: 0.736). Meanwhile, the risk model was associated with immune cells and immune-related functions. By combining the risk score and clinical factors, the nomogram of osteosarcoma was assessed with a C-index of 0.738 to further support the reliability of this 5-gene prognostic risk model. Finally, the expression of TcoF-related genes was validated in different cell lines by quantitative real-time PCR (qRT-PCR) and also in different tissue samples by immunohistochemistry (IHC). In conclusion, the model can predict the prognosis of osteosarcoma patients and may provide novel targets for the treatment of osteosarcoma patients.