Table 3_Longitudinal single-cell analysis identifies a chemotherapy-suppressed proliferative program associated with adverse prognosis in neuroblastoma.xlsx

Prognostic value was evaluated in a primary diagnostic bulk RNA-seq validation cohort and three independent external diagnostic bulk RNA-seq cohorts. Clinical model increment, cross-cohort meta-analysis, clinical feature association, and Hallmark correlation analyses were performed.</p>Results<p>Paired pseudo-bulk analysis showed that PTX malignant-neuroblast profiles had reduced activity of a proliferative-replication program enriched for E2F targets, G2/M checkpoint, MYC targets, DNA repair, and mitotic spindle pathways, together with relative enrichment of adaptive-state features related to hypoxia, inflammatory signaling, apoptosis, and epithelial-mesenchymal transition.

In the primary validation cohort, high PTX-down proliferation score was associated with inferior OS and EFS and remained prognostic after adjustment for age, COG-defined high-risk classification, MYCN status, and INSS stage. Adding the PTX-down score improved C-index and 1-, 3-, and 5-year time-dependent AUC beyond both COG and clinical models. Across four bulk RNA-seq cohorts, higher PTX-down score was consistently associated with adverse OS (random-effects pooled HR = 1.83, 95% CI 1.20–2.80; P = 0.005).

The score was also higher in tumors with MYCN amplification, stage 4 disease, and age ≥ 18 months. Core-gene and Hallmark correlation analyses defined the PTX-down score as an E2F/G2M/MYC-driven proliferative axis.</p>Conclusion<p>Longitudinal single-cell analysis identifies a chemotherapy-suppressed proliferative program in malignant neuroblasts. Projection of this program onto diagnostic bulk tumors captures high-risk biology, improves prognostic discrimination, and predicts adverse outcome in neuroblastoma.</p>