Expression profiling in glaucomatous human lamina cribrosa cells based on graph-clustering approach

In primary open angle glaucoma (POAG) patients, elevated intraocular pressure usually leads to extracellular matrix remodeling and astrocytes activation. Thus, lamina cribrosa (LC) cells may play an important role in POAG progression. The objective of this study was to comprehensively explore gene expression profiles in LC cells of POAG patients.

Using the GSE13534 microarray datasets downloaded from Gene Expression Omnibus database, the differentially expressed genes (DEGs) between LC cells from POAG patients and controls were firstly screened based on the classical t-test and false discovery rate <0.05 as a significant threshold. Subsequently, these DEGs were grouped into gene sets using a graph-clustering approach. The underlying molecular mechanisms were investigated by the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis.

A total of 57 DEGs were identified and 478 co-expression relationships were constructed among these DEGs. Among them, cytochrome p450 family 1 subfamily B (CYP1B1), brain-derived neurotrophic factor (BDNF) and myelin basic protein (MBP) showed high-degree relationships and they could interact with several genes. CYP1B1 is an important genetic gene involved in POAG and BDNF is an effective growth neurotrophic factor to weak POAG damage. MBP, versican (VCAN), integrin, alpha 4 (ITGA4) and N-cadherin (CDH2) may be involved in extracellular matrix remodeling in LC cells. FZD2 and FZD7 were enriched in basal cell carcinoma pathway.

The results demonstrate that the genes above may be associated with the pathogenesis of POAG.