Multitrait analysis of glaucoma identifies new risk loci and enables polygenic prediction of disease susceptibility and progression

Primary open angle glaucoma (POAG) is a leading cause of blindness world-wide. To identify new susceptibility loci, we meta-analyzed GWAS results from 8 independent studies from the United States (3,853 cases and 33,480 controls) and investigated the most significant SNPs in two Australian studies (1,252 cases and 2,592 controls), 3 European studies (875 cases and 4,107 controls) and a Singaporean Chinese study (1,037 cases and 2,543 controls). A meta-analysis of top SNPs identified three novel loci: rs35934224[T] within TXNRD2 (odds ratio (OR) = 0.78, P = 4.05×10−11 encoding a mitochondrial protein required for redox homeostasis; rs7137828[T] within ATXN2 (OR = 1.17, P = 8.73×10−10), and rs2745572[A] upstream of FOXC1 (OR = 1.17, P = 1.76×10−10). Using RT-PCR and immunohistochemistry, we show TXNRD2 and ATXN2 expression in retinal ganglion cells and the optic nerve head. These results identify new pathways underlying POAG susceptibility and suggest novel targets for preventative therapies.

Intraocular pressure (IOP) is currently the sole modifiable risk factor for primary open-angle glaucoma (POAG), one of the leading causes of blindness worldwide. Both IOP and POAG are highly heritable. We report a combined analysis of participants from the UK Biobank (n = 103,914) and previously published data from the International Glaucoma Genetic Consortium (n = 29,578) that identified 101 statistically independent genome-wide-significant SNPs for IOP, 85 of which have not been previously reported. We examined these SNPs in 11,018 glaucoma cases and 126,069 controls, and 53 SNPs showed evidence of association. Gene-based tests implicated an additional 22 independent genes associated with IOP. We derived an allele score based on the IOP loci and loci influencing optic nerve head morphology. In 1,734 people with advanced glaucoma and 2,938 controls, participants in the top decile of the allele score were at increased risk (odds ratio (OR) = 5.6; 95% confidence interval (CI): 4.1-7.6) of glaucoma relative to the bottom decile.

Primary open-angle glaucoma (POAG), the most common optic neuropathy, is a heritable disease. Siblings of POAG cases have a ten-fold increased risk of developing the disease. Intraocular pressure (IOP) and optic nerve head characteristics are used clinically to predict POAG risk. We conducted a genome-wide association meta-analysis of IOP and optic disc parameters and validated our findings in multiple sets of POAG cases and controls. Using imputation to the 1000 genomes (1000G) reference set, we identified 9 new genomic regions associated with vertical cup-disc ratio (VCDR) and 1 new region associated with IOP. Additionally, we found 5 novel loci for optic nerve cup area and 6 for disc area. Previously it was assumed that genetic variation influenced POAG either through IOP or via changes to the optic nerve head; here we present evidence that some genomic regions affect both IOP and the disc parameters. We characterized the effect of the novel loci through pathway analysis and found that pathways involved are not entirely distinct as assumed so far. Further, we identified a novel association between CDKN1A and POAG. Using a zebrafish model we show that six6b (associated with POAG and optic nerve head variation) alters the expression of cdkn1a. In summary, we have identified several novel genes influencing the major clinical risk predictors of POAG and showed that genetic variation in CDKN1A is important in POAG risk.

Primary open-angle glaucoma (POAG) is a major cause of irreversible blindness worldwide. We performed a genome-wide association study in an Australian discovery cohort comprising 1,155 advanced POAG cases and 1,992 controls. Association of the top SNPs from the discovery stage was investigated in two Australian replication cohorts (total 932 cases, 6,862 controls) and two US replication cohorts (total 2,616 cases, 2,634 controls). Meta-analysis of all cohorts revealed three novel loci associated with development of POAG. These loci are located upstream of ABCA1 (rs2472493 [G] OR=1.31, P= 2.1 × 10-19), within AFAP1 (rs4619890 [G] OR=1.20, P= 7.0 × 10-10) and within GMDS (rs11969985 [G] OR=1.31, and P= 7.7 × 10-10). Using RT-PCR and immunolabelling, we also showed that these genes are expressed within human retina, optic nerve and trabecular meshwork and that ABCA1 and AFAP1 are also expressed in retinal ganglion cells.