Probing the Potential of Defense Response-Associated Genes for Predicting the Progression, Prognosis, and Immune Microenvironment of Osteosarcoma

Abstract

Background: The defense response is a type of self-protective response of the body that protects it from damage by pathogenic factors. Although these reactions make important contributions to the occurrence and development of tumors, the role they play in osteosarcoma (OS), particularly in the immune microenvironment, remains unpredictable. Methods: This study included the clinical information and transcriptomic data of 84 osteosarcoma samples and the microarray data of 12 mesenchymal stem cell samples and 84 osteosarcoma samples. We obtained 129 differentially expressed genes related to the defense response (DRGs) by taking the intersection of differentially expressed genes with genes involved in the defense response pathway, and prognostic genes were screened using univariate Cox regression. Least absolute shrinkage and selection operator (LASSO) penalized Cox regression and multivariate Cox regression were then used to establish a DRG prognostic signature (DGPS) via the stepwise method. DGPS performance was examined using independent prognostic analysis, survival curves, and receiver operating characteristic (ROC) curves. In addition, the molecular and immune mechanisms of adverse prognosis in high-risk populations identified by DGPS were elucidated. The results were well verified by experiments. Result: BNIP3, PTGIS, and ZYX were identified as the most important DRGs for OS progression (hazard ratios of 2.044, 1.485, and 0.189, respectively). DGPS demonstrated outstanding performance in the prediction of OS prognosis (area under the curve (AUC) values of 0.842 and 0.787 in the training and test sets, respectively, adj-p < 0.05 in the survival curve). DGPS also performed better than a recent clinical prognostic approach with an AUC value of only 0.674 [metastasis], which was certified in the subsequent experimental results. These three genes regulate several key biological processes, including immune receptor activity and T cell activation, and they also reduce the infiltration of some immune cells, such as B cells, CD8+ T cells, and macrophages. Encouragingly, we found that DGPS was associated with sensitivity to chemotherapeutic drugs including JNK Inhibitor VIII, TGX221, MP470, and SB52334. Finally, we verified the effect of BNIP3 on apoptosis, proliferation, and migration of osteosarcoma cells through experiments. Conclusions: This study elucidated the role and mechanism of BNIP3, PTGIS, and ZYX in OS progression and was well verified by the experimental results, enabling reliable prognostic means and treatment strategies to be proposed for OS patients.

Integrating single‐cell and multi‐omic approaches reveals Euphorbiae Humifusae Herba‐dependent mitochondrial dysfunction in non‐small‐cell lung cancer

“…Different EHH concentrations were applied to the cell cultures of both A549 and NCI‐H1299. Following a 48‐h period of incubation, each well received 10 μL of CCK‐8 reagent (sourced from Jindo, Japan) and was further incubated for 2 h. 32 , 33 …”

Comprehensive characterization of stemness-related lncRNAs in triple-negative breast cancer identified a novel prognostic signature related to treatment outcomes, immune landscape analysis and therapeutic guidance: a silico analysis with in vivo experiments

“…To explore immune microenvironment in TNBC patients, differences in immune cell infiltration between the two subgroups were evaluated using the CIBERSORT algorithm [ 36 , 37 ]. The CIBERSORT algorithm is a computational tool designed for characterizing the composition of immune cell populations within complex tissues based on gene expression profiles.…”

Revealing the role of regulatory T cells in the tumor microenvironment of lung adenocarcinoma: a novel prognostic and immunotherapeutic signature

“…Gene set variation analysis (GSVA) and gene set enrichment analysis (GSEA) were employed to evaluate the biological characteristics ( 25 , 26 ). The ssGSEA approach was employed to quantify the enrichment scores of 29 immune signatures ( 27 , 28 ).…”