Identification of glioblastoma-specific prognostic biomarkers via an integrative analysis of DNA methylation and gene expression

Glioblastoma (GBM) is the most aggressive and lethal tumor of the central nervous system. The present study set out to identify reliable prognostic and predictive biomarkers for patients with GBM. RNA-sequencing data were obtained from The Cancer Genome Atlas database and DNA methylation data were downloaded using the University of California Santa Cruz-Xena database. The expression and methylation differences between patients with GBM, and survival times <1 and ≥1 year were investigated. A protein-protein interaction network was constructed and functional enrichment analyses of differentially expressed and methylated genes were performed. Hub genes were identified using the Cytoscape plug-in cytoHubba software. Survival analysis was performed using the survminer package, in order to determine the prognostic values of the hub genes. The present study identified 71 genes that were hypomethylated and expressed at high levels, and four genes that were hypermethylated and expressed at low levels in GBM. These genes were predominantly enriched in the ‘JAK-STAT signaling pathway’, ‘transcriptional misregulation in cancer’ and the ‘ECM-receptor interaction’, which are associated with GBM development. Among the 24 hub genes identified, 15 possessed potential prognostic value. An integrative analysis approach was implemented in order to analyze the association of DNA methylation with changes in gene expression and to assess the association of gene expression changes with GBM survival time. The results of the present study suggest that these 15 CpG-based genes may be useful and practical tools in predicting the prognosis of patients with GBM. However, future research on gene methylation and/or expression is required in order to develop personalized treatments for patients with GBM.