Identification of key modules and hub genes in glioblastoma multiforme based on co-expression network analysis

Glioblastoma multiforme (GBM) is the most malignant primary tumour in the central nervous system, but the molecular mechanisms underlying its pathogenesis remain unclear. In this study, data set GSE50161 was used to construct a co-expression network for weighted gene co-expression network analysis. Two modules (dubbed brown and turquoise) were found to have the strongest correlation with GBM. Functional enrichment analysis indicated that the brown module was involved in the cell cycle, DNA replication, and pyrimidine metabolism. The turquoise module was primarily related to circadian rhythm entrainment, glutamatergic synapses, and axonal guidance. Hub genes were screened by survival analysis using The Cancer Genome Atlas and Human Protein Atlas databases and further tested using the GSE4290 and Gene Expression Profiling Interactive Analysis databases. The eight hub genes (NUSAP1, SHCBP1, KNL1, SULT4A1, SLC12A5, NUF2, NAPB, and GARNL3) were verified at both the transcriptional and translational levels, and these gene expression levels were significant based on the World Health Organization classification system. These hub genes may be potential biomarkers and therapeutic targets for the accurate diagnosis and management of GBM.