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

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.

No abstract

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.

Primary open-angle glaucoma is the most common optic neuropathy and an important cause of irreversible blindness worldwide. The optic nerve head or optic disc is divided in two parts: a central cup (without nerve fibers) surrounded by the neuroretinal rim (containing axons of the retinal ganglion cells). The International Glaucoma Genetics Consortium conducted a meta-analysis of genome-wide association studies consisting of 17,248 individuals of European ancestry and 6,841 individuals of Asian ancestry. The outcomes of the genome-wide association studies were disc area and cup area. These specific measurements describe optic nerve morphology in another way than the vertical cup-disc ratio, which is a clinically used measurement, and may shed light on new glaucoma mechanisms. We identified 10 new loci associated with disc area (CDC42BPA, F5, DIRC3, RARB, ABI3BP, DCAF4L2, ELP4, TMTC2, NR2F2, and HORMAD2) and another 10 new loci associated with cup area (DHRS3, TRIB2, EFEMP1, FLNB, FAM101, DDHD1, ASB7, KPNB1, BCAS3, and TRIOBP). The new genes participate in a number of pathways and future work is likely to identify more functions related to the pathogenesis of glaucoma.

Variation in FOXC1 and PITX2 is associated with Axenfeld-Rieger syndrome, characterised by structural defects of the anterior chamber of the eye and a range of systemic features. Approximately half of all affected individuals will develop glaucoma, but the age at diagnosis and the phenotypic spectrum have not been well defined. As phenotypic heterogeneity is common, we aimed to delineate the age-related penetrance and the full phenotypic spectrum of glaucoma in FOXC1 or PITX2 carriers recruited through a national disease registry. All coding exons of FOXC1 and PITX2 were directly sequenced and multiplex ligation-dependent probe amplification was performed to detect copy number variation. The cohort included 53 individuals from 24 families with disease-associated FOXC1 or PITX2 variants, including one individual diagnosed with primary congenital glaucoma and five with primary open-angle glaucoma. The overall prevalence of glaucoma was 58.5% and was similar for both genes (53.3% for FOXC1 vs 60.9% for PITX2, P=0.59), however, the median age at glaucoma diagnosis was significantly lower in FOXC1 (6.0±13.0 years) compared with PITX2 carriers (18.0±10.6 years, P=0.04). The penetrance at 10 years old was significantly lower in PITX2 than FOXC1 carriers (13.0% vs 42.9%, P=0.03) but became comparable at 25 years old (71.4% vs 57.7%, P=0.38). These findings have important implications for the genetic counselling of families affected by Axenfeld-Rieger syndrome, and also suggest that FOXC1 and PITX2 contribute to the genetic architecture of primary glaucoma subtypes.

PURPOSE To investigate the presence of TBK1 copy number variations in a large, well-characterized Australian cohort of patients with glaucoma comprising both normal-tension glaucoma and high-tension glaucoma cases. DESIGN A retrospective cohort study. METHODS DNA samples from patients with normal-tension glaucoma and high-tension glaucoma and unaffected controls were screened for TBK1 copy number variations using real-time quantitative polymerase chain reaction. Samples with additional copies of the TBK1 gene were further tested using custom comparative genomic hybridization arrays. RESULTS Four out of 334 normal-tension glaucoma cases (1.2%) were found to carry TBK1 copy number variations using quantitative polymerase chain reaction. One extra dose of the TBK1 gene (duplication) was detected in 3 normal-tension glaucoma patients, while 2 extra doses of the gene (triplication) were detected in a fourth normal-tension glaucoma patient. The results were further confirmed by custom comparative genomic hybridization arrays. Further, the TBK1 copy number variation segregated with normal-tension glaucoma in the family members of the probands, showing an autosomal dominant pattern of inheritance. No TBK1 copy number variations were detected in 1045 Australian patients with high-tension glaucoma or in 254 unaffected controls. CONCLUSION We report the presence of TBK1 copy number variations in our Australian normal-tension glaucoma cohort, including the first example of more than 1 extra copy of this gene in glaucoma patients (gene triplication). These results confirm TBK1 to be an important cause of normal-tension glaucoma, but do not suggest common involvement in high-tension glaucoma.