Single-cell and bulk transcriptomic integration identifies a four-gene prognostic score linked to immunometabolic features in osteosarcoma
- 1epartment of Neurosurgery and Oncology, College of Medicine, The Affiliated Hospital of Southwest Jiaotong University, The Third People's Hospital of Chengdu, Chengdu, Sichuan, China
Abstract
BackgroundRecurrent and metastatic osteosarcoma calls for prognostic markers that connect tumor-cell states with the metabolic and microenvironmental context of disease progression. We investigated whether recurrence-derived single-cell signals could identify an outcome-associated transcriptional program and functionally informative candidate genes.MethodsSingle-cell data from GSE152048 were integrated with the 84-patient TARGET-OS cohort to derive a four-gene RiskScore and characterize molecular subtypes, pathway activity, and immune/stromal associations. COL13A1- and KIF25-directed perturbations were examined using transcript and protein measurements, colony formation, and wound-closure assays in 143B and HOS cells.ResultsOxidative phosphorylation enrichment characterized both recurrence-derived osteoblast candidates and an outcome-associated molecular subtype. The CGREF1–COL13A1–KIF25–RPL35AP7 score achieved an apparent C-index of 0.814 and 1-, 3-, and 5-year areas under the receiver operating characteristic curve of 0.855, 0.852, and 0.845. RiskScore remained associated with overall survival after adjustment for age, sex, and metastatic status (hazard ratio 2.413, 95% confidence interval 1.769–3.290). Higher scores tracked lower B-lineage, neutrophil, and fibroblast signals. Gene-directed treatments reduced wound closure in both cell lines. COL13A1-directed treatment reduced colony formation in both cell lines, and KIF25-directed treatment reduced colony formation in HOS cells (Dunnett-adjusted P < 0.05).ConclusionThis four-gene signature connects osteosarcoma outcome stratification with an immunometabolic and stromal context and identifies COL13A1 and KIF25 as candidates for focused functional investigation.
Related
- Clinical, single-cell, and spatial transcriptomic integration identifies an RCD-high inflammatory myeloid program for early neurological deterioration and outcome stratification after acute ischemic stroke
- Multi-Omic Integration of the miRNome and Transcriptome in Dermal Fibroblasts from Hypermobile Ehlers–Danlos Syndrome and Hypermobility Spectrum Disorders Identifies Convergent Candidate miRNA-mRNA Regulatory Networks Linking Extracellular Matrix Remodeling and Inflammation
- Integrated Spatial and Bulk Transcriptomic Analysis Identifies a Multicellular Ecosystem Involved in Progression of Young Breast Cancer.
- Integrated analysis of transcriptomic and proteomic profiles identifies divergent regulatory networks for intramuscular fat in Chinese versus Western Pigs
- Single-Cell Transcriptomic Analysis Identifies a Translation/RBP-Related Fibroblast Signature for Drug Repurposing in Idiopathic Pulmonary Fibrosis
BibTeX
@article{li2026singlecell,
title = {{Single-cell and bulk transcriptomic integration identifies a four-gene prognostic score linked to immunometabolic features in osteosarcoma}},
author = {Jingyue Li and Yijin Wang and Dan Wen and Shanshan Zhao and Qiang Li and Xinhao Peng and Jun Zhang and Sunfu Zhang},
journal = {Frontiers in Cell and Developmental Biology},
year = {2026},
month = nov,
doi = {10.3389/fcell.2026.2010244},
url = {https://doi.org/10.3389/fcell.2026.2010244},
note = {CC BY},
}
RIS
TY - JOUR
TI - Single-cell and bulk transcriptomic integration identifies a four-gene prognostic score linked to immunometabolic features in osteosarcoma
AU - Jingyue Li
AU - Yijin Wang
AU - Dan Wen
AU - Shanshan Zhao
AU - Qiang Li
AU - Xinhao Peng
AU - Jun Zhang
AU - Sunfu Zhang
PY - 2026
DA - 2026/11/13
JO - Frontiers in Cell and Developmental Biology
DO - 10.3389/fcell.2026.2010244
UR - https://doi.org/10.3389/fcell.2026.2010244
AB - BackgroundRecurrent and metastatic osteosarcoma calls for prognostic markers that connect tumor-cell states with the metabolic and microenvironmental context of disease progression. We investigated whether recurrence-derived single-cell signals could identify an outcome-associated transcriptional program and functionally informative candidate genes.MethodsSingle-cell data from GSE152048 were integrated with the 84-patient TARGET-OS cohort to derive a four-gene RiskScore and characterize molecular subtypes, pathway activity, and immune/stromal associations. COL13A1- and KIF25-directed perturbations were examined using transcript and protein measurements, colony formation, and wound-closure assays in 143B and HOS cells.ResultsOxidative phosphorylation enrichment characterized both recurrence-derived osteoblast candidates and an outcome-associated molecular subtype. The CGREF1–COL13A1–KIF25–RPL35AP7 score achieved an apparent C-index of 0.814 and 1-, 3-, and 5-year areas under the receiver operating characteristic curve of 0.855, 0.852, and 0.845. RiskScore remained associated with overall survival after adjustment for age, sex, and metastatic status (hazard ratio 2.413, 95% confidence interval 1.769–3.290). Higher scores tracked lower B-lineage, neutrophil, and fibroblast signals. Gene-directed treatments reduced wound closure in both cell lines. COL13A1-directed treatment reduced colony formation in both cell lines, and KIF25-directed treatment reduced colony formation in HOS cells (Dunnett-adjusted P < 0.05).ConclusionThis four-gene signature connects osteosarcoma outcome stratification with an immunometabolic and stromal context and identifies COL13A1 and KIF25 as candidates for focused functional investigation.
SN - 2296-634X
ER -