Abstract:Glaucoma is a neurodegenerative disease characterized by progressive optic nerve degeneration and retinal ganglion cell (RGC) loss. In early glaucoma, before obvious axon loss, highly organized pathological processes in RGCs occur sequentially, involving axons, dendrites and synaptic terminals. The optic nerve head (ONH) is the critical structure of early glaucomatous neurodegeneration. Taking advantage of high-throughput data from the ONH and the weighted gene coexpression network analysis (WGCNA) method, the current study aims to gain insight into the full scope of pathological events in early glaucoma and define their chronological sequence. The expression profiles of GSE26299, GSE110019, and GSE139605, which measure ONH gene expression in different glaucoma models, were downloaded from the Gene Expression Omnibus (GEO) database. In GSE26299, which uses 10.5-month-old DBA/2 J mice, WGCNA was utilized to construct a gene coexpression network, and the most significant modules of early (NOE), moderate (MOD) and severe (SEV) glaucoma were identified. The differentially expressed genes (DEGs) of GSE110019 and GSE139605 significantly overlapped with the correlated module of the MOD group, so the 3 gene sets were analyzed together. Pathway enrichment analysis via the GO, KEGG, and Reactome pathways was subsequently performed, followed by protein‒protein interaction (PPI) analysis to screen key genes associated with each stage. Several hub gene expression patterns were identified in a glucocorticoid-induced glaucoma (GIG) model via quantitative PCR and immunostaining. The pink module was positively correlated with the NOE group (r = 0.48, p = 4e-04) and negatively correlated with the glaucoma stage (r = -0.88, p = 3e-17). The genes in the pink module were enriched in the synaptic transmission and axonal transport pathways. The tan module was negatively correlated with the NOE group (r = -0.43, p = 0.002) and positively correlated with the glaucoma stage (r = 0.77, p = 7e-11). The genes in the tan module were associated with pathways such as tight junctions, retinol metabolism, and linoleic acid metabolism. The purple module was positively correlated with the MOD group (r = 0.64, p = 5e-07). The common genes among the purple module and the DEGs of the two other datasets were enriched in pathways related to mitotic cell division, cytokine activity, and the extracellular matrix (ECM). The hub genes identified by PPI included Nrn1 , Cplx1 , Timp1 , and Cdk1 . Quantitative PCR and immunostaining confirmed that Limk1 expression was increased in the ONH of GIG mice. In early glaucomatous neuropathy, intrinsic changes in RGCs precede the activation of glial cells and ECM remodeling. These latter events are common pathological changes observed in the ONH in both cats and mice. Our study may provide new targets for the early detection and treatment of glaucoma.

P16INK4a Upregulation Mediated by SIX6 Defines Retinal Ganglion Cell Pathogenesis in Glaucoma

Integrating Genetic Regulation and Single-Cell Expression with GWAS Prioritizes Causal Genes and Cell Types for Glaucoma

Modeling Glaucoma: Retinal Ganglion Cells Generated from Induced Pluripotent Stem Cells of Patients with SIX6 Risk Allele Show Developmental Abnormalities

High throughput functional profiling of genes at intraocular pressure loci reveals distinct networks for glaucoma

Osteopontin Activates Retinal Microglia Causing Retinal Ganglion Cells Loss Via P38 MAPK Signaling Pathway in Glaucoma.

IOP-induced blood-retinal barrier compromise contributes to RGC death in glaucoma

Trpc6: An Underlying Target For Human Glaucoma

Identification and Validation of Key Biomarkers and Potential Therapeutic Targets for Primary Open-Angle Glaucoma.

Retinal ganglion cell-specific genetic regulation in primary open-angle glaucoma

Identification and validation of oxidative stress-related genes in primary open-angle glaucoma by weighted gene co-expression network analysis and machine learning

Gut-licensed β7 + CD4 + T cells contribute to progressive retinal ganglion cell damage in glaucoma

POU6F2, a risk factor for glaucoma, myopia and dyslexia, labels specific populations of retinal ganglion cells

Integrative transcriptomic and proteomic analysis reveals CD9/ITGA4/PI3K‐Akt axis mediates trabecular meshwork cell apoptosis in human glaucoma

Increased expression of the transient receptor potential cation channel 6 gene in patients with primary open-angle glaucoma.

Mutant RAMP2 Causes Primary Open-Angle Glaucoma Via the CRLR-cAMP Axis

Integrated Aqueous Humor Cerna and Mirna-Tf-mrna Network Analysis Reveals Potential Molecular Mechanisms Governing Primary Open-Angle Glaucoma Pathogenesis

Microglial cGAS-STING signaling underlies glaucoma pathogenesis

Identification of Sequential Molecular Mechanisms and Key Biomarkers in Early Glaucoma by Integrated Bioinformatics Analysis

Elevation of P-Nr2as1232 by Cdk5/p35 Contributes to Retinal Ganglion Cell Apoptosis in a Rat Experimental Glaucoma Model

A Missense Mutation in Pitx2 Leads to Early-Onset Glaucoma Via NRF2-YAP1 Axis