Comprehensive bioinformatics analysis of trabecular meshwork gene expression data to unravel the molecular pathogenesis of primary open-angle glaucoma

Purpose Performing bioinformatics analyses using trabecular meshwork (TM) gene expression data in order to further elucidate the molecular pathogenesis of primary open-angle glaucoma (POAG), and to identify candidate target genes. Methods A systematic search in Gene Expression Omnibus and ArrayExpress was conducted, and quality control and preprocessing of the data was performed with ArrayAnalysis.org. Molecular pathway overrepresentation analysis was performed with PathVisio using pathway content from three pathway databases: WikiPathways, KEGG and Reactome. In addition, Gene Ontology (GO) analysis was performed on the gene expression data. The significantly changed pathways were clustered into functional categories which were combined into a network of connected genes. Results Ninety-two significantly changed pathways were clustered into five functional categories: extracellular matrix (ECM), inflammation, complement activation, senescence and Rho GTPase signalling. ECM included pathways involved in collagen, actin and cell-matrix interactions. Inflammation included pathways entailing NF-kappa B and arachidonic acid. The network analysis showed that several genes overlap between the inflammation cluster on the one hand, and the ECM, complement activation and senescence clusters on the other hand. GO analysis, identified additional clusters, related to development and corticosteroids. Conclusion This study provides an overview of the processes involved in the molecular pathogenesis of POAG in the TM. The results show good face validity and confirm findings from histological, biochemical, genome-wide association and transcriptomics studies. The identification of known points of action for drugs, such as Rho GTPase, arachidonic acid, NF-kappa B, prostaglandins and corticosteroid clusters, supports the value of this approach to identify potential drug targets.