Integrated aqueous humor ceRNA and miRNA-TF-mRNA network analysis reveals potential molecular mechanisms governing primary open-angle glaucoma pathogenesis

Primary open-angle glaucoma (POAG) is the leading cause of blindness globally, which develops through complex and poorly understood biological mechanisms. Herein, we conducted an integrated bioinformatics analysis of extant aqueous humor (AH) gene expression datasets in order to identify key genes and regulatory mechanisms governing POAG progression. We downloaded AH gene expression datasets (GSE101727 and GSE105269) corresponding to healthy controls and POAG patients from the Gene Expression Omnibus. We then identified mRNAs, microRNAs (miRNAs), and long non-coding RNAs (lncRNAs) that were differentially expressed (DE) between control and POAG patients. DEmRNAs and DElncRNAs were then subjected to pathway enrichment analyses, after which a protein-protein interaction (PPI) network was generated. This network was then expanded to establish lncRNA-miRNA-mRNA and miRNA-transcription factor(TF)-mRNA networks. In total, the GSE101727 dataset was used to identify 2746 DElncRNAs and 2208 DEmRNAs, while the GSE105269 dataset was used to identify 45 DEmiRNAs. We ultimately constructed a competing endogenous RNA (ceRNA) network incorporating 37, 5, and 14 of these lncRNAs, miRNAs and mRNAs, respectively. The proteins encoded by these 14 hub mRNAs were found to be significantly enriched for activities that may be linked to POAG pathogenesis. In addition, we generated a miRNA-TF-mRNA regulatory network containing 2 miRNAs (miR-135a-5p and miR-139-5p), 5 TFs (TGIF2, TBX5, HNF1A, TCF3, and FOS) and 5 mRNAs (SHISA7, ST6GAC2, TXNIP, FOS, and DCBLD2). The SHISA7, ST6GAC2, TXNIP, FOS, and DCBLD2 genes that may be viable therapeutic targets for the prevention or treatment of POAG, and regulated by the TFs (TGIF2, HNF1A, TCF3, and FOS).