AIM: To explore hub genes for glaucoma based on bioinformatics analysis and an experimental model verification. METHODS: In the Gene Expression Omnibus (GEO) database, the GSE25812 and GSE26299 datasets were selected to analyze differentially expressed genes (DEGs) by the GEO2R tool. Through bioinformatics analysis, 9 hub genes were identified. Receiver operating characteristic (ROC) curves and principal component analysis (PCA) were performed to verify whether the hub gene can distinguish glaucoma from normal eyes. The mouse model of glaucoma was constructed, and the real-time reverse transcriptase-polymerase chain reaction (RT-qPCR) assay was performed to detect the expression levels of hub genes in glaucoma. RESULTS: There were 128 overlapping DEGs in the GSE25812 and GSE26299 datasets, mainly involved in intracellular signalling, cell adhesion molecules and the Ras signalling pathway. A total of 9 hub genes were screened out, including GNAL, BGN, ETS2, FCGP4, MAPK10, MMP15, STAT1, TSPAN8, and VCAM1. The area under the curve (AUC) values of 9 hub genes were greater than 0.8. The PC1 axle could provide a 70.5% interpretation rate to distinguish glaucoma from normal eyes. In the ocular tissues of glaucoma in the mice model, the expression of BGN, ETS2, FCGR4, STAT1, TSPAN8, and VCAM1 was increased, while the expression of GNAL, MAPK10, and MMP15 was decreased. CONCLUSION: Nine hub genes in glaucoma are identified, which may provide new biomarkers and therapeutic targets for glaucoma.