Integrative transcriptomic, spatial and functional-genomic analysis identifies a UFMylation-related vascular-stromal program and prioritizes WWTR1 in glioblastoma.

Abstract

Glioblastoma (GBM) contains spatially organized stress-adaptive and vascular niches. Because transcript abundance does not measure UFM1 conjugation, we asked whether a UFMylation-related transcriptional axis identifies a reproducible tissue program and alters candidate prioritization. In 518 unique primary TCGA-GBM tumors profiled on the Affymetrix HT Human Genome U133A array, weighted gene co-expression network analysis of 8,000 variable genes yielded 12 modules. The 278-gene green module ranked first across nine prespecified traits (mean |r|=0.637). Direct overlap comprised 1/3 measurable UFMylation-core, 5/19 ER-stress/UPR, and 2/15 proteostasis genes; after excluding overlapping genes, correlations with the green eigengene remained significant (r = 0.373, 0.831, 0.639, and 0.699 for UFMylation-core, ER-stress/UPR, proteostasis, and composite scores, respectively). The green score was associated with overall survival per standard-deviation increase (HR 1.17, 95% CI 1.07-1.28), although clinical adjustment attenuated the estimate. In a 10-sample single-cell dataset, sample-level scores were higher in pericytes and endothelial cells than in malignant cells. Donor-aware IvyGAP analysis supported regional organization, whereas one Visium section showed stronger concordance with ER-stress/UPR and mesenchymal scores than with the UFMylation-core score. CellChat indicated pathway-selective rather than global remodeling of inferred vascular communication. Layer ablation moved WWTR1 from rank 48 using WGCNA alone to rank 4 overall and rank 1 among non-common-essential genes after cross-platform integration. These findings define an ER-stress/mesenchymal-weighted, UFMylation-related vascular-stromal transcriptional association and nominate WWTR1 for experimental testing.