Glioblastoma multiforme (GBM) is the most common and malignant brain tumor of the adult central nervous system and is associated with poor prognosis. The present study aimed to identify the hub genes in GBM in order to improve the current understanding of the underlying mechanism of GBM. The RNA‑seq data were downloaded from The Cancer Genome Atlas database. The edgeR package in R software was used to identify differentially expressed genes (DEGs) between two groups: Glioblastoma samples and normal brain samples. Gene Ontology (GO) functional enrichment analysis and the Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis were performed using Database for Annotation, Visualization and Integrated Discovery software. Additionally, Cytoscape and Search Tool for the Retrieval of Interacting Genes/Proteins tools were used for the protein‑protein interaction network, while the highly connected modules were extracted from this network using the Minimal Common Oncology Data Elements plugin. Next, the prognostic significance of the candidate hub genes was analyzed using UALCAN. In addition, the identified hub genes were verified by reverse transcription‑quantitative (RT‑q) PCR. In total, 1,483 DEGs were identified between GBM and control samples, including 954 upregulated genes and 529 downregulated genes (P<0.01; fold‑change >16) and these genes were involved in different GO terms and signaling pathways. Furthermore, CDK1, BUB1, BUB1B, CENPA and GNG3 were identified as key genes in the GBM samples. The UALCAN tool verified that higher expression level of CENPA was relevant to poorer overall survival rates. In conclusion, CDK1, BUB1, BUB1B, CENPA and GNG3 were found to be potential biomarkers for GBM. Additionally, 'cell cycle' and 'γ‑aminobutyric acid signaling' pathways may serve a significant role in the pathogenesis of GBM.
Keywords: glioblastoma multiforme; differentially expressed genes; functional enrichment analyses; Kyoto encyclopedia of Genes and Genomes pathway; prognosis.