# Prognostic value and immune infiltration of anoikis-related genes in osteosarcoma

**Authors:** Jianming Mo, Hening Li, Hui Wang, Xu Fang, Jie Ma, Kaiwei Chen

PMC · DOI: 10.3389/fmed.2025.1669470 · 2025-11-06

## TL;DR

This study builds a model using genes related to anoikis to predict osteosarcoma prognosis and immune response, validated with single-cell RNA sequencing and experiments.

## Contribution

The novel contribution is a prognostic model for osteosarcoma based on anoikis-related genes and single-cell RNA sequencing data.

## Key findings

- A four-gene model (MYC, BNIP3, IGFBP5, SPP1) accurately predicts osteosarcoma prognosis with AUC values of 0.836 for 1-, 3-, and 5-year survival.
- The model shows significant differences in immune cell infiltration between risk groups and correlates with immune checkpoint expression.
- qRT-PCR and immunohistochemistry confirmed elevated gene expression in osteosarcoma cells and poor prognosis patients.

## Abstract

Currently, there is no research on building osteosarcoma (OS) prognostic models based on single-cell RNA sequencing (scRNA-seq) and anoikis-related genes (ARGs).

Differential genes between osteoblasts cells and osteosarcoma cells were identified using scRNA-seq, and ARGs were determined by Genecard database. Lasso regression was employed to investigate hub genes and construct the model based on TARGET. Kaplan-Meier survival analysis was applied to compare the survival differences. ROC curves were used to evaluate the predictive performance of the model. CIBERSORT and ESTIMATE algorithms were conducted to calculate immune cell infiltration abundance. Finally, qRT-PCR and immunohistochemistry experiments were conducted to validate the results.

A predictive model containing four modeling genes (MYC, BNIP3, IGFBP5, and SPP1) was successfully constructed, with AUC values of 0.836, 0.837, and 0.836 for 1-, 3-, and 5-year patient prognosis, respectively. Importantly, the model also showed good predictive value in two validation set. The infiltration of immune cells in different risk groups showed significant differences. The modeling genes were associated with the expression of various immune checkpoints and the response to immune therapy. qRT-PCR showed MYC, BNIP3, IGFBP5, and SPP1 substantially exhibited a trend of high expression in osteosarcoma cells. Immunohistochemistry suggested that in osteosarcoma patients with poorer prognosis, the expression of these four hub genes was significantly elevated.

We have developed an effective model for predicting osteosarcoma prognosis and immune response, which may provide valuable insights for osteosarcoma prognostic evaluation and immune therapy strategies.

## Linked entities

- **Genes:** MYC (MYC proto-oncogene, bHLH transcription factor) [NCBI Gene 4609], BNIP3 (BCL2 interacting protein 3) [NCBI Gene 664], IGFBP5 (insulin like growth factor binding protein 5) [NCBI Gene 3488], SPP1 (secreted phosphoprotein 1) [NCBI Gene 6696]
- **Diseases:** osteosarcoma (MONDO:0002623)

## Full-text entities

- **Genes:** MYC (MYC proto-oncogene, bHLH transcription factor) [NCBI Gene 4609] {aka MRTL, MYCC, bHLHe39, c-Myc}, BNIP3 (BCL2 interacting protein 3) [NCBI Gene 664] {aka HABON, NIP3}, IGFBP5 (insulin like growth factor binding protein 5) [NCBI Gene 3488] {aka IBP5}, SPP1 (secreted phosphoprotein 1) [NCBI Gene 6696] {aka BNSP, BSPI, ETA-1, OPN}
- **Diseases:** OS (MESH:D012516)
- **Species:** Homo sapiens (human, species) [taxon 9606]

## Figures

8 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12631436/full.md

---
Source: https://tomesphere.com/paper/PMC12631436