Whole exome sequencing implicates eye development, the unfolded protein response and plasma membrane homeostasis in primary open-angle glaucoma

Purpose To identify biological processes associated with POAG and its subtypes, high-tension (HTG) and normal-tension glaucoma (NTG), by analyzing rare potentially damaging genetic variants. Methods A total of 122 and 65 unrelated HTG and NTG participants, respectively, with early onset advanced POAG, 103 non-glaucoma controls and 993 unscreened ethnicity-matched controls were included in this study. Study participants without myocilin disease-causing variants and non-glaucoma controls were subjected to whole exome sequencing on an Illumina HiSeq2000. Exomes of participants were sequenced on an Illumina HiSeq2000. Qualifying variants were rare in the general population (MAF < 0.001) and potentially functionally damaging (nonsense, frameshift, splice or predicted pathogenic using SIFT or Polyphen2 software). Genes showing enrichment of qualifying variants in cases were selected for pathway and network analysis using InnateDB. POAG cases showed enrichment of rare variants in camera-type eye development genes (p = 1.40×10-7, corrected p = 3.28×10-4). Implicated eye development genes were related to neuronal or retinal development. HTG cases were significantly enriched for key regulators in the unfolded protein response (UPR) (p = 7.72×10-5, corrected p = 0.013). The UPR is known to be involved in myocilin-related glaucoma; our results suggest the UPR has a role in non-myocilin causes of HTG. NTG cases showed enrichment in ion channel transport processes (p = 1.05×10-4, corrected p = 0.027) including calcium, chloride and phospholipid transporters involved in plasma membrane homeostasis. Network analysis also revealed enrichment of the MHC Class I antigen presentation pathway in HTG, and the EGFR1 and cell-cycle pathways in both HTG and NTG.