Ferroptosis Transcriptional Regulation and Prognostic Impact in Medulloblastoma Subtypes Revealed by RNA-Seq

Medulloblastoma (MB) is the most common malignant brain tumor in children, typically arising during infancy and childhood. Despite multimodal therapies achieving a response rate of 70% in children older than 3 years, treatment remains challenging. Ferroptosis, a form of regulated cell death, can be induced in medulloblastoma cells in vitro using erastin or RSL3. Using two independent medulloblastoma RNA-sequencing cohorts (MB-PBTA and MTAB-10767), we investigated the expression of ferroptosis-related molecules through multiple approaches, including Weighted Gene Co-Expression Network Analysis (WGCNA), molecular subtype stratification, protein-protein interaction (PPI) networks, and univariable and multivariable overall survival analyses. A prognostic expression score was computed based on a cross-validated ferroptosis signature. In training and validation cohorts, the regulation of the ferroptosis transcriptional program distinguished the four molecular subtypes of medulloblastoma. WGCNA identified nine gene modules in the MB tumor transcriptome; five correlated with molecular subtypes, implicating pathways related to oxidative stress, hypoxia, and trans-synaptic signaling. One module, associated with disease recurrence, included epigenetic regulators and nucleosome organizers. Univariable survival analyses identified a 45-gene ferroptosis prognostic signature associated with nutrient sensing, cysteine and methionine metabolism, and trans-sulfuration within a one-carbon metabolism. The top ten unfavorable ferroptosis genes included CCT3, SNX5, SQOR, G3BP1, CARS1, SLC39A14, FAM98A, FXR1, TFAP2C, and ATF4. Patients with a high ferroptosis score showed a worse prognosis, particularly in the G3 and SHH subtypes. The PPI network highlighted IL6 and CBS as unfavorable hub genes. In a multivariable overall survival model, which included gender, age, and the molecular subtype classification, the ferroptosis expression score was validated as an independent adverse prognostic marker (hazard ratio: 5.8; p-value = 1.04 × 10^-9). This study demonstrates that the regulation of the ferroptosis transcriptional program is linked to medulloblastoma molecular subtypes and patient prognosis. A cross-validated ferroptosis signature was identified in two independent RNA-sequencing cohorts, and the ferroptosis score was confirmed as an independent and adverse prognostic factor in medulloblastoma.