Non-negative matrix factorization and differential expression analyses identify hub genes linked to progression and prognosis of glioblastoma multiforme

[Display omitted] •GUCA1A, RFC2, GNG11, MMP19, and NRG1 were screened as the real hub genes for the upcoming molecular studies in GBM.•Mitogen-activated protein kinase (MAPK) signaling pathway was found to be the most significant pathway.•The novel application of the NMF bioinformatics pipeline was an effective method to elucidate commonalities and discrepancies between samples of the dataset. One of the most prevailing primary brain tumors in adult human male is glioblastoma multiforme (GBM), which is categorized by rapid cellular growth. Even though the combination therapy comprises surgery, chemotherapy, and adjuvant therapies, the survival rate, on average, is 14.6 months. Glioma stem cells (GSCs) have key roles in tumorigenesis, progression, and defiance against chemotherapy and radiotherapy. In our study, firstly, the gene expression dataset GSE124145 was retrieved; the non-negative matrix factorization (NMF) method was applied on GBM dataset, and differentially expressed genes analysis (DEGs) was performed. After which, overlapping genes between metagenes and DEGs were detected to examine the Gene Ontology (GO) categories in the biological process (BP) in the stemness of GBM. The common hub genes were used to construct protein–protein interaction (PPI) network and further GO, while Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway was utilized to pinpoint the real hub genes. The analysis of hub genes particular for the same GO categories demonstrated that specific hub genes triggered distinct features of the same biological processes. After utilizing GSE124145 and The Cancer Genome Atlas (TCGA) dataset for survival analysis, we screened five real hub genes: GUCA1A, RFC2, GNG11, MMP19, and NRG1, which are strongly associated with the progression and prognosis of GBM. The DEGs analysis revealed that all real hub genes were overexpressed in GBM and TCGA datasets, which further validates our results. The constructed study of PPI, GO, and KEGG pathway on common hub genes was performed. Finally, the KEGG pathways performed on the top 15 candidate hub genes (including six real hub genes) of the PPI network in the GBM gene expression dataset study found mitogen-activated protein kinase (Mapk) signaling pathway to be the most significant pathway. The rest of the hub genes reviewed throughout the analysis might be favorable targets for diagnosing and treating GBM and lower-grade gliomas.