Identification of genes involved in glaucoma pathogenesis using combined network analysis and empirical studies

Glaucoma is a leading cause of blindness. We aimed in this study to identify genes that may make subtle and cumulative contributions to glaucoma pathogenesis. To this end, we identified molecular interactions and pathways that include transcription factors (TFs) FOXC1, PITX2, PAX6 and NFKB1 and various microRNAs including miR-204 known to have relevance to trabecular meshwork (TM) functions and/or glaucoma. TM tissue is involved in glaucoma pathogenesis. In-house microarray transcriptome results and data sources were used to identify target genes of the regulatory molecules. Bioinformatics analyses were done to filter TM and glaucoma relevant genes. These were submitted to network-creating softwares to define interactions, pathways and a network that would include the genes. The network was stringently scrutinized and minimized, then expanded by addition of microarray data and data on TF and microRNA-binding sites. Selected features of the network were confirmed by empirical studies such as dual luciferase assays, real-time PCR and western blot experiments and apoptosis assays. MYOC, WDR36, LTPBP2, RHOA, CYP1B1, OPA1, SPARC, MEIS2, PLEKHG5, RGS5, BBS5, ALDH1A1, NOMO2, CXCL6, FMNL2, ADAMTS5, CLOCK and DKK1 were among the genes included in the final network. Pathways identified included those that affect ECM properties, IOP, ciliary body functions, retinal ganglion cell viability, apoptosis, focal adhesion and oxidative stress response. The identification of many genes potentially involved in glaucoma pathology is consistent with its being a complex disease. The inclusion of several known glaucoma-related genes validates the approach used.