Title : Survival-concordant biomarker discovery in glioblastoma using a diversity-preserving multi-tier machine-learning framework Background: Robust transcriptomic biomarker discovery in glioblastoma multiforme (GBM) is limited by cohort heteroge...
Title : Survival-concordant biomarker discovery in glioblastoma using a diversity-preserving multi-tier machine-learning framework Background: Robust transcriptomic biomarker discovery in glioblastoma multiforme (GBM) is limited by cohort heterogeneity, small sample sizes, and algorithm-dependent feature selection. Conventional machine-learning pipelines often prioritize classification accuracy without considering whether gene-expression changes are concordant with survival risk. Therefore, we developed a survival-concordant, diversity-preserving multi-tier machine-learning framework for biomarker discovery across heterogeneous GBM transcriptomic cohorts. Results: RNA-seq data from heterogeneous public cohorts (TCGA, CGGA, and GEO) were integrated to evaluate three continuous disease-state transitions: Normal-to-Primary (NvP), Normal-to-Recurrenct (NvR), and Primary-to-Recurrenct (PvR). The framework combined a 223-model first-tier machine-learning tournament, diversity-preserving ensemble selection, an expression–prognosis discordance filter, and bootstrap-LASSO-weighted secondary meta-ensemble scoring. The finalized locked-model pipeline achieved robust discriminative performance with AUC values of 0.993, 0.976, and 0.681 for the internal NvP, NvR, and PvR tracks, respectively. Projection onto independent external GEO cohorts (GSE121720 and GSE263588) demonstrated outstanding generalizability (AUCs of 0.966 and 1.000) alongside a 100% feature preservation rate. Integrated multi-cohort prioritization (TCGA and CGGA) successfully identified mechanistically distinct diagnostic modules led by SCN3B (NvP), ORC1 (NvR), and IGFBP2 (PvR). Notably, IGFBP2 showed progressive, monotonic upregulation from normal tissue to primary and recurrent GBM, correlating with poorer survival. Single-cell and immune analyses revealed that the recurrent-specific IGFBP2 and collagen module (COL3A1, COL4A1) reflects malignant, vascular, and myeloid-associated microenvironmental remodeling. Conclusion: This study developed a survival-concordant, diversity-preserving multi-tier machine-learning framework for transcriptomic biomarker discovery across heterogeneous public cohorts. By integrating multi-cohort RNA-seq analysis, ensemble model selection, expression–prognosis discordance filtering, and bootstrap-LASSO-weighted secondary scoring, the framework prioritized biomarkers with both discriminatory value and survival-concordant directionality. Application to glioblastoma successfully decoded three mechanistically distinct disease-state programs: the loss of normal neural identity (SCN3B), the acquisition of unrestrained proliferation (ORC1), and extracellular matrix (ECM) remodeling during recurrence (IGFBP2 and the associated collagen module). Notably, IGFBP2 emerged as a critical bridge marker across the disease continuum, reflecting malignant-cell, vascular, and immune-associated features within the recurrent microenvironment. These findings support the utility of the proposed framework for prioritizing biologically interpretable biomarkers from public transcriptomic datasets and provide a basis for future experimental validation and extension to other cancer types or multi-omic data. Key word: Glioblastoma, machine learning, multi-tier framework