Abstract
Glioblastoma (GBM) is the most frequent malignant brain tumor in adults. Our study focused on uncovering differentially expressed genes (DEGs) and their methylation in order to identify novel diagnostic biomarkers and potential treatment targets. Using GBM RNA-sequencing data from The Cancer Genome Atlas (TCGA) database, DEGs between GBM samples and paracancer tissue samples were analyzed. Enrichment analysis for DEGs and transcription factors (TFs) was performed. A total of 1029 upregulated genes and 1542 downregulated genes were identified, which were associated mainly with multiple tumor-related and immune-related pathways such as cell cycle, mitogen-activated protein kinase signaling pathway, leukocyte transendothelial migration, and autoimmune thyroid disease. These DEGs were enriched for 174 TFs, and six TFs were differentially expressed and identified as key TFs in GBM: HOXA3, EN1, ZIC1, and FOXD3 were upregulated, while HLF and EGR3 were downregulated. A total of 1978 DEGs were involved in the regulatory networks of the six key differentially expressed TFs. High expression of EN1 was associated with shorter overall survival, while high expression of EGR3 was associated with shorter recurrence-free survival. The six TFs were differentially methylated in GBM samples compared with paracancer tissues. Our study identifies numerous DEGs and their associated pathways as potential contributors to GBM, particularly the TFs EN1, EGR3, HOXA3, ZIC1, FOXD3, and HLF. The differential expression of these TFs may be unlikely driven by aberrant methylation. These TFs may be useful as diagnostic markers and treatment targets in GBM, and EN1 and EGR3 may have predictive prognostic value.