Glaucoma is a multi-factorial blinding disease in which genetic factors play an important role. Elevated intraocular pressure is a highly heritable risk factor for primary open angle glaucoma and currently the only target for glaucoma therapy. Our study helps to better understand underlying genetic and molecular mechanisms that regulate intraocular pressure, and identifies a new candidate gene, Cacna2d1, that modulates intraocular pressure and a promising therapeutic, pregabalin, which binds to CACNA2D1 protein and lowers intraocular pressure significantly. Because our study utilizes a genetically diverse population of mice with known sequence variants, we are able to determine that the intraocular pressure-lowering effect of pregabalin is dependent on the Cacna2d1 haplotype. Using human genome-wide association study (GWAS) data, evidence for association of a CACNA2D1 single-nucleotide polymorphism and primary open angle glaucoma is found. Importantly, these results demonstrate that our systems genetics approach represents an efficient method to identify genetic variation that can guide the selection of therapeutic targets.

A new avenue of mining published genome-wide association studies includes the joint analysis of related traits. The power of this approach depends on the genetic correlation of traits, which reflects the number of pleiotropic loci, i.e. genetic loci influencing multiple traits. Here, we applied new meta-analyses of optic nerve head (ONH) related traits implicated in primary open-angle glaucoma (POAG); intraocular pressure and central corneal thickness using Haplotype reference consortium imputations. We performed a multi-trait analysis of ONH parameters cup area, disc area and vertical cup-disc ratio. We uncover new variants; rs11158547 in PPP1R36-PLEKHG3 and rs1028727 near SERPINE3 at genome-wide significance that replicate in independent Asian cohorts imputed to 1000 Genomes. At this point, validation of these variants in POAG cohorts is hampered by the high degree of heterogeneity. Our results show that multi-trait analysis is a valid approach to identify novel pleiotropic variants for ONH.

IMPORTANCE Genetic variants associated with primary open-angle glaucoma (POAG) are known to influence disease risk. However, the clinical effect of associated variants individually or in aggregate is not known. Genetic risk scores (GRS) examine the cumulative genetic load by combining individual genetic variants into a single measure, which is assumed to have a larger effect and increased power to detect relevant disease-related associations. OBJECTIVE To investigate if a GRS that comprised 12 POAG genetic risk variants is associated with age at disease diagnosis. DESIGN, SETTING, AND PARTICIPANTS A cross-sectional study included individuals with POAG and controls from the Glaucoma Genes and Environment (GLAUGEN) study and the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) study. A GRS was formulated using 12 variants known to be associated with POAG, and the alleles associated with increasing risk of POAG were aligned in the case-control sets. In case-only analyses, the association of the GRS with age at diagnosis was analyzed as an estimate of disease onset. Results from cohort-specific analyses were combined with meta-analysis. Data collection started in August 2012 for the NEIGHBOR cohort and in July 2008 for the GLAUGEN cohort and were analyzed starting in March 2018. MAIN OUTCOMES AND MEASURES Association of a 12 single-nucleotide polymorphism POAG GRS with age at diagnosis in individuals with POAG using linear regression. RESULTS The GLAUGEN study included 976 individuals with POAG and 1140 controls. The NEIGHBOR study included 2132 individuals with POAG and 2290 controls. For individuals with POAG, the mean (SD) age at diagnosis was 63.6 (9.8) years in the GLAUGEN cohort and 66.0 (13.7) years in the NEIGHBOR cohort. For controls, the mean (SD) age at enrollment was 65.5 (9.2) years in the GLAUGEN cohort and 68.9 (11.4) years in the NEIGHBOR cohort. All study participants were European white. The GRS was strongly associated with POAG risk in case-control analysis (odds ratio per 1-point increase in score = 1.24; 95% CI, 1.21-1.27; P = 3.4 × 10 −66). In case-only analyses, each higher GRS unit was associated with a 0.36-year earlier age at diagnosis (β = −0.36; 95% CI, −0.56 to −0.16; P = 4.0 × 10 −4). Individuals in the top 5% of the GRS had a mean (SD) age at diagnosis of 5.2 (12.8) years earlier than those in the bottom 5% GRS (61.4 [12.7] vs 66.6 [12.9] years; P = 5.0 × 10 −4). CONCLUSIONS AND RELEVANCE A higher dose of POAG risk alleles was associated with an earlier age at glaucoma diagnosis. On average, individuals with POAG with the highest GRS had 5.2-year earlier age at diagnosis of disease. These results suggest that a GRS that comprised genetic variants associated with POAG could help identify patients with risk of earlier disease onset impacting screening and therapeutic strategies.

PurposeTo identify genetic risk factors contributing to central corneal thickness (CCT) in individuals from South India, a population with a high prevalence of ocular disorders.MethodsOne hundred ninety-five individuals from 15 large South Indian pedigrees were genotyped using the Omni2.5 bead array. Family-based association for CCT was conducted using the score test in MERLIN.ResultsGenome-wide association study (GWAS) identified strongest association for single nucleotide polymorphisms (SNPs) in the first intron of WNT7B and CCT (top SNP rs9330813; β = −0.57, 95% confidence interval [CI]: −0.78 to −0.36; P = 1.7 × 10−7). We further investigated rs9330813 in a Latino cohort and four independent European cohorts. A meta-analysis of these data sets demonstrated statistically significant association between rs9330813 and CCT (β = −3.94, 95% CI: −5.23 to −2.66; P = 1.7 × 10−9). WNT7B SNPs located in the same genomic region that includes rs9330813 have previously been associated with CCT in Latinos but with other ocular quantitative traits related to myopia (corneal curvature and axial length) in a Japanese population (rs10453441 and rs200329677). To evaluate the specificity of the observed WNT7B association with CCT in the South Indian families, we completed an ocular phenome-wide association study (PheWAS) for the top WNT7B SNPs using 45 ocular traits measured in these same families including corneal curvature and axial length. The ocular PheWAS results indicate that in the South Indian families WNT7B SNPs are primarily associated with CCT.ConclusionsThe results indicate robust evidence for association between WNT7B SNPs and CCT in South Indian pedigrees, and suggest that WNT7B SNPs can have population-specific effects on ocular quantitative traits.