To reveal the potential molecular mechanism of glioblastoma multiforme (GBM) and provide the candidate biomarkers for GBM gene therapy. Microarray dataset GSE50161 was obtained from GEO database. The differentially expressed genes (DEGs) were identified between GBM samples and control samples, followed by the module partition analysis based on WGCNA. Then, the pathway and functional enrichment analyses of DEGs were performed. The hub genes were further investigated, followed by the survival analysis and data validation. A total of 1913 DEGs were investigated between two groups, followed by analysis of 5 modules using WGCNA. These DEGs were mainly enriched in functions like inflammatory response. The hub genes including upregulated N-Myc and STAT Interactor (NMI), Capping Actin Protein-Gelsolin Like (CAPG), and Proteasome Subunit Beta 8 (PSMB8) were revealed as potential liquid biopsy molecules for GBM diagnose. Moreover, Nucleolar and Spindle Associated Protein 1 (NUSAP1) and G Protein-Coupled Receptor 65 (GPR65) were outstanding genes in survival analysis. Our results suggested that CPNE6, HAPLN2, CMTM3, NMI, CAPG, and PSMB8 might be used as potential molecules for liquid biopsy of GBM. NUSAP1 and GPR65 might be novel prognostic targets for GBM gene therapy. Furthermore, the upregulated NMI might play an important role in GBM progression via inflammatory response.