Abstract
Purpose The mechanisms contributing to recurrence of glioblastoma (GBM), an aggressive neuroepithelial brain tumor, remain unknown. We have recently shown that nuclear respiratory factor 1 (NRF1) is an oncogenic transcription factor and its transcriptional activity is associated with the progression and prognosis of GBM. Herein, we extend our efforts to (1) identify influential NRF1-driven gene and microRNA (miRNA) expression for the aggressiveness of mesenchymal GBM; and (2) understand the molecular basis for its poor response to therapy. Methods Clinical data and RNA-Seq from four independent GBM cohorts were analyzed by Bayesian Network Inference with Java Objects (BANJO) and Markov chain Monte Carlo (MCMC)-based gene order to identify molecular drivers of mesenchymal GBM as well as prognostic indicators of poor response to radiation and chemotherapy. Results We are the first to report sex-specific NRF1 motif enriched gene signatures showing increased susceptibility to GBM. Risk estimates for GBM were increased by greater than 100-fold with the joint effect of NRF1-driven gene signatures-CDK4, DUSP6, MSH2, NRF1, and PARK7 in female GBM patients and CDK4, CASP2, H6PD, and NRF1 in male GBM patients. NRF1-driven causal Bayesian network genes were predictive of poor survival and resistance to chemoradiation in IDH1 wild-type mesenchymal GBM patients. NRF1-regulatable miRNAs were also associated with poor response to chemoradiation therapy in female IDH1 wild-type mesenchymal GBM. Stable overexpression of NRF1 reprogramed human astrocytes into neural stem cell-like cells expressing SOX2 and nestin. These cells differentiated into neurons and form tumorospheroids. Conclusions In summary, our novel discovery shows that NRF1-driven causal genes and miRNAs involved in cancer cell stemness and mesenchymal features contribute to cancer aggressiveness and recurrence of aggressive therapy-resistant glioblastoma.
