Abstract:Background: Although the diagnosis and treatment of glioblastoma (GBM) is significantly improved with recent progresses, there is still a large heterogeneity in therapeutic effects and overall survival. The aim of this study is to analyze gene expressions of transcription factors (TFs) in GBM so as to discover new tumor markers. Methods: Differentially expressed TFs are identified by data mining using public databases. The GBM transcriptome profile is downloaded from The Cancer Genome Atlas (TCGA). The nonnegative matrix factorization (NMF) method is used to cluster the differentially expressed genes to discover hub genes and signal pathways. The TFs affecting the prognosis of GBM are screened by univariate and multivariate COX regression analysis, and the receiver operating characteristic (ROC) curve is determined. The GBM hazard model and nomogram map are constructed by integrating the clinical data. Finally, the TFs involving potential signaling pathways in GBM are screened by Gene Set Enrichment Analysis (GSEA), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Results: There are 68 differentially expressed TFs in GBM, of which 43 genes are upregulated and 25 genes are downregulated. NMF clustering analysis suggested that GBM patients are divided into three groups: Clusters A, B, and C. LHX2, MEOX2, SNAI2, and ZNF22 are identified from the above differential genes by univariate/multivariate regression analysis. The risk score of those four genes are calculated based on the beta coefficient of each gene, and we found that the predictive ability of the risk score gradually increased with the prolonged predicted termination time by time-dependent ROC curve analysis. The nomogram results have showed that the integration of risk score, age, gender, chemotherapy, radiotherapy, and 1p/19q can further improve predictive ability towards the survival of GBM. The pathways in cancer, phosphoinositide 3-kinases (PI3K)-Akt signaling, Hippo signaling, and proteoglycans, are highly enriched in high-risk groups by GSEA. These genes are mainly involved in cell migration, cell adhesion, epithelial-mesenchymal transition (EMT), cell cycle, and other signaling pathways by GO and KEGG analysis. Conclusion: The four-factor combined scoring model of LHX2, MEOX2, SNAI2, and ZNF22 can precisely predict the prognosis of patients with GBM.

Predicting Glioblastoma Prognosis Networks Using Weighted Gene Co-Expression Network Analysis on TCGA Data

Identification of Glioblastoma Gene Prognosis Modules Based on Weighted Gene Co-Expression Network Analysis

Identification of a Specific Gene Module for Predicting Prognosis in Glioblastoma Patients

Identification of candidate biomarkers for GBM based on WGCNA

Identify Signature Regulatory Network for Glioblastoma Prognosis by Integrative Mrna and Mirna Co-Expression Analysis.

Predicting prognostic markers for glioma using gene co-expression network analysis.

Mining prognostic marker of glioma based on weighted gene co-expression network analysis

Weighted Gene Correlation Network Analysis Identifies Microenvironment-Related Genes Signature As Prognostic Candidate for Grade II/III Glioma.

Weighted Gene Coexpression Network Analysis Identifies an Immunogenic Cell Death Signature for Predicting Prognosis and Therapeutic Responses in Glioblastoma

Identification of Survival‑associated Key Genes and Long Non‑coding RNAs in Glioblastoma Multiforme by Weighted Gene Co‑expression Network Analysis.

Identification of Key Genes Involved in the Recurrence of Glioblastoma Multiforme Using Weighted Gene Co-Expression Network Analysis and Differential Expression Analysis.

A Weighted Gene Co-Expression Network Analysis–Derived Prognostic Model for Predicting Prognosis and Immune Infiltration in Gastric Cancer

A Novel Prognostic Signature of Transcription Factors for the Prediction in Patients with GBM.

Identification of Key Immune and Cell Cycle Modules and Prognostic Genes for Glioma Patients through Transcriptome Analysis

Identification of Key Genes Related to Lung Cancer with Brain Metastasis through Weighted Gene Correlation Network Analysis

Identification of WHO II/III Gliomas by 16 Prognostic-related Gene Signatures using Machine Learning Methods

Bioinformatics Analysis of High-Throughput Data to Validate Potential Novel Biomarkers and Small Molecule Drugs for Glioblastoma Multiforme

Identification of novel prognostic targets in glioblastoma using bioinformatics analysis

Weighted gene correlation network analysis identifies RSAD2, HERC5, and CCL8 as prognostic candidates for breast cancer.

Four Specific Biomarkers Associated with the Progression of Glioblastoma Multiforme in Older Adults Identified Using Weighted Gene Co-Expression Network Analysis

Potential Prognostic Predictors and Molecular Targets for Skin Melanoma Screened by Weighted Gene Co-expression Network Analysis