Search

Primary open-angle glaucoma (POAG), is a heritable common cause of blindness world-wide. To identify risk loci, we conduct a large multi-ethnic meta-analysis of genome-wide association studies on a total of 34,179 cases and 349,321 controls, identifying 44 previously unreported risk loci and confirming 83 loci that were previously known. The majority of loci have broadly consistent effects across European, Asian and African ancestries. Cross-ancestry data improve fine-mapping of causal variants for several loci. Integration of multiple lines of genetic evidence support the functional relevance of the identified POAG risk loci and highlight potential contributions of several genes to POAG pathogenesis, including SVEP1, RERE, VCAM1, ZNF638, CLIC5, SLC2A12, YAP1, MXRA5, and SMAD6. Several drug compounds targeting POAG risk genes may be potential glaucoma therapeutic candidates.

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.

ImportancePrimary open-angle glaucoma presents with increased prevalence and a higher degree of clinical severity in populations of African ancestry compared with European or Asian ancestry. Despite this, individuals of African ancestry remain understudied in genomic research for blinding disorders.ObjectivesTo perform a genome-wide association study (GWAS) of African ancestry populations and evaluate potential mechanisms of pathogenesis for loci associated with primary open-angle glaucoma.Design, settings, and participantsA 2-stage GWAS with a discovery data set of 2320 individuals with primary open-angle glaucoma and 2121 control individuals without primary open-angle glaucoma. The validation stage included an additional 6937 affected individuals and 14 917 unaffected individuals using multicenter clinic- and population-based participant recruitment approaches. Study participants were recruited from Ghana, Nigeria, South Africa, the United States, Tanzania, Britain, Cameroon, Saudi Arabia, Brazil, the Democratic Republic of the Congo, Morocco, Peru, and Mali from 2003 to 2018. Individuals with primary open-angle glaucoma had open iridocorneal angles and displayed glaucomatous optic neuropathy with visual field defects. Elevated intraocular pressure was not included in the case definition. Control individuals had no elevated intraocular pressure and no signs of glaucoma.ExposuresGenetic variants associated with primary open-angle glaucoma.Main outcomes and measuresPresence of primary open-angle glaucoma. Genome-wide significance was defined as P C) with primary open-angle glaucoma (odds ratio [OR], 1.32 [95% CI, 1.20-1.46]; P = 2 × 10-8). The association was validated in an analysis of an additional 6937 affected individuals and 14 917 unaffected individuals (OR, 1.15 [95% CI, 1.09-1.21]; P

Importance:Primary open-angle glaucoma presents with increased prevalence and a higher degree of clinical severity in populations of African ancestry compared with European or Asian ancestry. Despite this, individuals of African ancestry remain understudied in genomic research for blinding disorders. Objectives:To perform a genome-wide association study (GWAS) of African ancestry populations and evaluate potential mechanisms of pathogenesis for loci associated with primary open-angle glaucoma. Design, Settings, and Participants:A 2-stage GWAS with a discovery data set of 2320 individuals with primary open-angle glaucoma and 2121 control individuals without primary open-angle glaucoma. The validation stage included an additional 6937 affected individuals and 14 917 unaffected individuals using multicenter clinic- and population-based participant recruitment approaches. Study participants were recruited from Ghana, Nigeria, South Africa, the United States, Tanzania, Britain, Cameroon, Saudi Arabia, Brazil, the Democratic Republic of the Congo, Morocco, Peru, and Mali from 2003 to 2018. Individuals with primary open-angle glaucoma had open iridocorneal angles and displayed glaucomatous optic neuropathy with visual field defects. Elevated intraocular pressure was not included in the case definition. Control individuals had no elevated intraocular pressure and no signs of glaucoma. Exposures:Genetic variants associated with primary open-angle glaucoma. Main Outcomes and Measures:Presence of primary open-angle glaucoma. Genome-wide significance was defined as P C) with primary open-angle glaucoma (odds ratio [OR], 1.32 [95% CI, 1.20-1.46]; P = 2 × 10-8). The association was validated in an analysis of an additional 6937 affected individuals and 14 917 unaffected individuals (OR, 1.15 [95% CI, 1.09-1.21]; P

ImportanceGenetic 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.ObjectiveTo investigate if a GRS that comprised 12 POAG genetic risk variants is associated with age at disease diagnosis.Design, setting, and participantsA 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 measuresAssociation of a 12 single-nucleotide polymorphism POAG GRS with age at diagnosis in individuals with POAG using linear regression.ResultsThe 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 relevanceA 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.

Importance:Mutations in the myocilin (MYOC) gene are the most common molecularly defined cause of primary open-angle glaucoma that typically occurs in patients with high intraocular pressures (IOP). One MYOC mutation, p.Gln368Ter, has been associated with as many as 1.6% of primary open-angle glaucoma cases that had a mean maximum recorded IOP of 30 mm Hg. However, to our knowledge, the role of the p.Gln368Ter mutation in patients with normal-tension glaucoma (NTG) with an IOP of 21 mm Hg or lower has not been investigated. Objective:To evaluate the role of the p.Gln368Ter MYOC mutation in patients with NTG. Design, Setting, and Participants:In this case-control study of the prevalence of the p.Gln368Ter mutation in patients with NTG, cohort 1 was composed of 772 patients with NTG and 2152 controls from the United States (Iowa, Minnesota, and New York) and England and cohort 2 was composed of 561 patients with NTG and 2606 controls from the Massachusetts Eye and Ear Infirmary and the NEIGHBORHOOD consortium. Genotyping was conducted using real-time polymerase chain reaction that was confirmed with Sanger sequencing, the imputation of genome-wide association study data, or an analysis of whole-exome sequence data. Data analysis occurred between April 2007 and April 2018. Main Outcomes and Measures:Comparison of the frequency of the p.Gln368Ter MYOC mutation between NTG cases and controls with the Fisher exact test. Results:Of 6091 total participants, 3346 (54.9%) were women and 5799 (95.2%) were white. We detected the p.Gln368Ter mutation in 7 of 772 patients with NTG (0.91%) and 7 of 2152 controls (0.33%) in cohort 1 (P = .03). In cohort 2, we detected the p.Gln368Ter mutation in 4 of 561 patients with NTG (0.71%) and 10 of 2606 controls (0.38%; P = .15). When the cohorts were analyzed as a group, the p.Gln368Ter mutation was associated with NTG (odds ratio, 2.3; 95% CI, 0.98-5.3; P = .04). Conclusions and Relevance:In cohorts 1 and 2, the p.Gln368Ter mutation in MYOC was found in patients with IOPs that were 21 mm Hg or lower (NTG), although at a frequency that is lower than previously detected in patients with higher IOP. These data suggest that the p.Gln368Ter mutation may be associated with glaucoma in patients with normal IOPs as well as in patients with IOPs that are greater than 21 mm Hg.

Emmanuelle Souzeau, who contributed to analysis of data, was inadvertently omitted from the author list in the originally published version of this Article. This has now been corrected in both the PDF and HTML versions of the Article.

Primary open angle glaucoma (POAG) is a complex disease with a major genetic contribution. Its prevalence varies greatly among ethnic groups, and is up to five times more frequent in black African populations compared to Europeans. So far, worldwide efforts to elucidate the genetic complexity of POAG in African populations has been limited. We conducted a genome-wide association study in 1113 POAG cases and 1826 controls from Tanzanian, South African and African American study samples. Apart from confirming evidence of association at TXNRD2 (rs16984299; OR[T] 1.20; P = 0.003), we found that a genetic risk score combining the effects of the 15 previously reported POAG loci was significantly associated with POAG in our samples (OR 1.56; 95% CI 1.26-1.93; P = 4.79 × 10-5). By genome-wide association testing we identified a novel candidate locus, rs141186647, harboring EXOC4 (OR[A] 0.48; P = 3.75 × 10-8), a gene transcribing a component of the exocyst complex involved in vesicle transport. The low frequency and high degree of genetic heterogeneity at this region hampered validation of this finding in predominantly West-African replication sets. Our results suggest that established genetic risk factors play a role in African POAG, however, they do not explain the higher disease load. The high heterogeneity within Africans remains a challenge to identify the genetic commonalities for POAG in this ethnicity, and demands studies of extremely large size.

Central corneal thickness (CCT) is a highly heritable trait associated with complex eye diseases such as keratoconus and glaucoma. We perform a genome-wide association meta-analysis of CCT and identify 19 novel regions. In addition to adding support for known connective tissue-related pathways, pathway analyses uncover previously unreported gene sets. Remarkably, >20% of the CCT-loci are near or within Mendelian disorder genes. These included FBN1, ADAMTS2 and TGFB2 which associate with connective tissue disorders (Marfan, Ehlers-Danlos and Loeys-Dietz syndromes), and the LUM-DCN-KERA gene complex involved in myopia, corneal dystrophies and cornea plana. Using index CCT-increasing variants, we find a significant inverse correlation in effect sizes between CCT and keratoconus (r = -0.62, P = 5.30 × 10-5) but not between CCT and primary open-angle glaucoma (r = -0.17, P = 0.2). Our findings provide evidence for shared genetic influences between CCT and keratoconus, and implicate candidate genes acting in collagen and extracellular matrix regulation.

Purpose:Sex hormones may be associated with primary open-angle glaucoma (POAG), although the mechanisms are unclear. We previously observed that gene variants involved with estrogen metabolism were collectively associated with POAG in women but not men; here we assessed gene variants related to testosterone metabolism collectively and POAG risk. Methods:We used two datasets: one from the United States (3853 cases and 33,480 controls) and another from Australia (1155 cases and 1992 controls). Both datasets contained densely called genotypes imputed to the 1000 Genomes reference panel. We used pathway- and gene-based approaches with Pathway Analysis by Randomization Incorporating Structure (PARIS) software to assess the overall association between a panel of single nucleotide polymorphisms (SNPs) in testosterone metabolism genes and POAG. In sex-stratified analyses, we evaluated POAG overall and POAG subtypes defined by maximum IOP (high-tension [HTG] or normal tension glaucoma [NTG]). Results:In the US dataset, the SNP panel was not associated with POAG (permuted P = 0.77), although there was an association in the Australian sample (permuted P = 0.018). In both datasets, the SNP panel was associated with POAG in men (permuted P ≤ 0.033) and not women (permuted P ≥ 0.42), but in gene-based analyses, there was no consistency on the main genes responsible for these findings. In both datasets, the testosterone pathway association with HTG was significant (permuted P ≤ 0.011), but again, gene-based analyses showed no consistent driver gene associations. Conclusions:Collectively, testosterone metabolism pathway SNPs were consistently associated with the high-tension subtype of POAG in two datasets.

Purpose:Sex hormones may be associated with primary open-angle glaucoma (POAG), although the mechanisms are unclear. We previously observed that gene variants involved with estrogen metabolism were collectively associated with POAG in women but not men; here we assessed gene variants related to testosterone metabolism collectively and POAG risk. Methods:We used two datasets: one from the United States (3853 cases and 33,480 controls) and another from Australia (1155 cases and 1992 controls). Both datasets contained densely called genotypes imputed to the 1000 Genomes reference panel. We used pathway- and gene-based approaches with Pathway Analysis by Randomization Incorporating Structure (PARIS) software to assess the overall association between a panel of single nucleotide polymorphisms (SNPs) in testosterone metabolism genes and POAG. In sex-stratified analyses, we evaluated POAG overall and POAG subtypes defined by maximum IOP (high-tension [HTG] or normal tension glaucoma [NTG]). Results:In the US dataset, the SNP panel was not associated with POAG (permuted P = 0.77), although there was an association in the Australian sample (permuted P = 0.018). In both datasets, the SNP panel was associated with POAG in men (permuted P ≤ 0.033) and not women (permuted P ≥ 0.42), but in gene-based analyses, there was no consistency on the main genes responsible for these findings. In both datasets, the testosterone pathway association with HTG was significant (permuted P ≤ 0.011), but again, gene-based analyses showed no consistent driver gene associations. Conclusions:Collectively, testosterone metabolism pathway SNPs were consistently associated with the high-tension subtype of POAG in two datasets.

Central corneal thickness (CCT) is one of the most heritable ocular traits and it is also a phenotypic risk factor for primary open angle glaucoma (POAG). The present study uses the BXD Recombinant Inbred (RI) strains to identify novel quantitative trait loci (QTLs) modulating CCT in the mouse with the potential of identifying a molecular link between CCT and risk of developing POAG. The BXD RI strain set was used to define mammalian genomic loci modulating CCT, with a total of 818 corneas measured from 61 BXD RI strains (between 60-100 days of age). The mice were anesthetized and the eyes were positioned in front of the lens of the Phoenix Micron IV Image-Guided OCT system or the Bioptigen OCT system. CCT data for each strain was averaged and used to QTLs modulating this phenotype using the bioinformatics tools on GeneNetwork (www.genenetwork.org). The candidate genes and genomic loci identified in the mouse were then directly compared with the summary data from a human POAG genome wide association study (NEIGHBORHOOD) to determine if any genomic elements modulating mouse CCT are also risk factors for POAG.This analysis revealed one significant QTL on Chr 13 and a suggestive QTL on Chr 7. The significant locus on Chr 13 (13 to 19 Mb) was examined further to define candidate genes modulating this eye phenotype. For the Chr 13 QTL in the mouse, only one gene in the region (Pou6f2) contained nonsynonymous SNPs. Of these five nonsynonymous SNPs in Pou6f2, two resulted in changes in the amino acid proline which could result in altered secondary structure affecting protein function. The 7 Mb region under the mouse Chr 13 peak distributes over 2 chromosomes in the human: Chr 1 and Chr 7. These genomic loci were examined in the NEIGHBORHOOD database to determine if they are potential risk factors for human glaucoma identified using meta-data from human GWAS. The top 50 hits all resided within one gene (POU6F2), with the highest significance level of p = 10-6 for SNP rs76319873. POU6F2 is found in retinal ganglion cells and in corneal limbal stem cells. To test the effect of POU6F2 on CCT we examined the corneas of a Pou6f2-null mice and the corneas were thinner than those of wild-type littermates. In addition, these POU6F2 RGCs die early in the DBA/2J model of glaucoma than most RGCs. Using a mouse genetic reference panel, we identified a transcription factor, Pou6f2, that modulates CCT in the mouse. POU6F2 is also found in a subset of retinal ganglion cells and these RGCs are sensitive to injury.

Primary open-angle glaucoma (POAG) is the most common chronic optic neuropathy worldwide. Epidemiological studies show a robust positive relation between intraocular pressure (IOP) and POAG and modest positive association between IOP and blood pressure (BP), while the relation between BP and POAG is controversial. The International Glaucoma Genetics Consortium (n=27 558), the International Consortium on Blood Pressure (n=69 395), and the National Eye Institute Glaucoma Human Genetics Collaboration Heritable Overall Operational Database (n=37 333), represent genome-wide data sets for IOP, BP traits and POAG, respectively. We formed genome-wide significant variant panels for IOP and diastolic BP and found a strong relation with POAG (odds ratio and 95% confidence interval: 1.18 (1.14-1.21), P=1.8 × 10-27) for the former trait but no association for the latter (P=0.93). Next, we used linkage disequilibrium (LD) score regression, to provide genome-wide estimates of correlation between traits without the need for additional phenotyping. We also compared our genome-wide estimate of heritability between IOP and BP to an estimate based solely on direct measures of these traits in the Erasmus Rucphen Family (ERF; n=2519) study using Sequential Oligogenic Linkage Analysis Routines (SOLAR). LD score regression revealed high genetic correlation between IOP and POAG (48.5%, P=2.1 × 10-5); however, genetic correlation between IOP and diastolic BP (P=0.86) and between diastolic BP and POAG (P=0.42) were negligible. Using SOLAR in the ERF study, we confirmed the minimal heritability between IOP and diastolic BP (P=0.63). Overall, IOP shares genetic basis with POAG, whereas BP has limited shared genetic correlation with IOP or POAG.

ObjectiveSeveral attributes of female reproductive history, including age at natural menopause (ANM), have been related to primary open-angle glaucoma (POAG). We assembled 18 previously reported common genetic variants that predict ANM to determine their association with ANM or POAG.MethodsUsing data from the Nurses' Health Study (7,143 women), we validated the ANM weighted genetic risk score in relation to self-reported ANM. Subsequently, to assess the relation with POAG, we used data from 2,160 female POAG cases and 29,110 controls in the National Eye Institute Glaucoma Human Genetics Collaboration Heritable Overall Operational Database (NEIGHBORHOOD), which consists of 8 datasets with imputed genotypes to 5.6+ million markers. Associations with POAG were assessed in each dataset, and site-specific results were meta-analyzed using the inverse weighted variance method.ResultsThe genetic risk score was associated with self-reported ANM (P = 2.2 × 10) and predicted 4.8% of the variance in ANM. The ANM genetic risk score was not associated with POAG (Odds Ratio (OR) = 1.002; 95% Confidence Interval (CI): 0.998, 1.007; P = 0.28). No single genetic variant in the panel achieved nominal association with POAG (P ≥0.20). Compared to the middle 80 percent, there was also no association with the lowest 10 percentile or highest 90 percentile of genetic risk score with POAG (OR = 0.75; 95% CI: 0.47, 1.21; P = 0.23 and OR = 1.10; 95% CI: 0.72, 1.69; P = 0.65, respectively).ConclusionsA genetic risk score predicting 4.8% of ANM variation was not related to POAG; thus, genetic determinants of ANM are unlikely to explain the previously reported association between the two phenotypes.

The cornea is the outermost layer of the eye and is a vital component of focusing incoming light on the retina. Central corneal thickness (CCT) is now recognized to have a significant role in ocular health and is a risk factor for various ocular diseases, such as keratoconus and primary open angle glaucoma. Most previous genetic studies utilized European and Asian subjects to identify genetic loci associated with CCT. Minority populations, such as Latinos, may aid in identifying additional loci and improve our understanding of the genetic architecture of CCT. In this study, we conducted a genome-wide association study (GWAS) in Latinos, a traditionally understudied population in genetic research, to further identify loci contributing to CCT. Study participants were genotyped using either the Illumina OmniExpress BeadChip (∼730K markers) or the Illumina Hispanic/SOL BeadChip (∼2.5 million markers). All study participants were 40 years of age and older. We assessed the association between individual single nucleotide polymorphisms (SNPs) and CCT using linear regression, adjusting for age, gender and principal components of genetic ancestry. To expand genomic coverage and to interrogate additional SNPs, we imputed SNPs from the 1000 Genomes Project reference panels. We identified a novel SNP, rs10453441 (P = 6.01E-09), in an intron of WNT7B that is associated with CCT. Furthermore, WNT7B is expressed in the human cornea. We also replicated 11 previously reported loci, including IBTK, RXRA-COL5A1, COL5A1, FOXO1, LRRK1 and ZNF469 (P

Glaucoma is a highly heritable disease that is a leading cause of blindness worldwide. Here, we identified heterozygous thrombospondin 1 (THBS1) missense alleles altering p.Arg1034, a highly evolutionarily conserved amino acid, in 3 unrelated and ethnically diverse families affected by congenital glaucoma, a severe form of glaucoma affecting children. [Thbs1.sup.R1034C]-mutant mice had elevated intraocular pressure (IOP), reduced ocular fluid outflow, and retinal ganglion cell loss. Histology revealed an abundant, abnormal extracellular accumulation of THBS1 with abnormal morphology of juxtacanalicular trabecular meshwork (TM), an ocular tissue critical for aqueous fluid outflow. Functional characterization showed that the THBS1 missense alleles found in affected individuals destabilized the THBS1 C-terminus, causing protein misfolding and extracellular aggregation. Analysis using a range of amino acid substitutions at position R1034 showed that the extent of aggregation was correlated with the change in protein-folding free energy caused by variations in amino acid structure. Extracellular matrix (ECM) proteins, especially fibronectin, which bind to THBS1, also accumulated within THBS1 deposits. These results show that missense variants altering THBS1 p.Arg1034 can cause elevated IOP through a mechanism involving impaired TM fluid outflow in association with accumulation of aggregated THBS1 in the ECM of juxtacanalicular meshwork with altered morphology., Introduction Glaucoma is a progressive blinding disease that causes loss of retinal ganglion cells (RGCs) and irreversible degeneration of the optic nerve. Elevated intraocular pressure (IOP), currently the only modifiable [...]

PurposeRecent studies indicate that mitochondrial proteins may contribute to the pathogenesis of primary open-angle glaucoma (POAG). In this study, we examined the association between POAG and common variations in gene-encoding mitochondrial proteins.MethodsWe examined genetic data from 3430 POAG cases and 3108 controls derived from the combination of the GLAUGEN and NEIGHBOR studies. We constructed biological-system coherent mitochondrial nuclear-encoded protein gene-sets by intersecting the MitoCarta database with the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. We examined the mitochondrial gene-sets for association with POAG and with normal-tension glaucoma (NTG) and high-tension glaucoma (HTG) subsets using Pathway Analysis by Randomization Incorporating Structure.ResultsWe identified 22 KEGG pathways with significant mitochondrial protein-encoding gene enrichment, belonging to six general biological classes. Among the pathway classes, mitochondrial lipid metabolism was associated with POAG overall (P = 0.013) and with NTG (P = 0.0006), and mitochondrial carbohydrate metabolism was associated with NTG (P = 0.030). Examining the individual KEGG pathway mitochondrial gene-sets, fatty acid elongation and synthesis and degradation of ketone bodies, both lipid metabolism pathways, were significantly associated with POAG (P = 0.005 and P = 0.002, respectively) and NTG (P = 0.0004 and P < 0.0001, respectively). Butanoate metabolism, a carbohydrate metabolism pathway, was significantly associated with POAG (P = 0.004), NTG (P = 0.001), and HTG (P = 0.010).ConclusionsWe present an effective approach for assessing the contributions of mitochondrial genetic variation to open-angle glaucoma. Our findings support a role for mitochondria in POAG pathogenesis and specifically point to lipid and carbohydrate metabolism pathways as being important.

PurposeNoncoding microRNAs (miRNAs) have been implicated in the pathogenesis of glaucoma. We aimed to identify common variants in miRNA coding genes (MIR) associated with primary open-angle glaucoma (POAG).MethodsUsing the NEIGHBORHOOD data set (3853 cases/33,480 controls with European ancestry), we first assessed the relation between 85 variants in 76 MIR genes and overall POAG. Subtype-specific analyses were performed in high-tension glaucoma (HTG) and normal-tension glaucoma subsets. Second, we examined the expression of miR-182, which was associated with POAG, in postmortem human ocular tissues (ciliary body, cornea, retina, and trabecular meshwork [TM]), using miRNA sequencing (miRNA-Seq) and droplet digital PCR (ddPCR). Third, miR-182 expression was also examined in human aqueous humor (AH) by using miRNA-Seq. Fourth, exosomes secreted from primary human TM cells were examined for miR-182 expression by using miRNA-Seq. Fifth, using ddPCR we compared miR-182 expression in AH between five HTG cases and five controls.ResultsOnly rs76481776 in MIR182 gene was associated with POAG after adjustment for multiple comparisons (odds ratio [OR] = 1.23, 95% confidence interval [CI]: 1.11-1.42, P = 0.0002). Subtype analysis indicated that the association was primarily in the HTG subset (OR = 1.26, 95% CI: 1.08-1.47, P = 0.004). The risk allele T has been associated with elevated miR-182 expression in vitro. Data from ddPCR and miRNA-Seq confirmed miR-182 expression in all examined ocular tissues and TM-derived exosomes. Interestingly, miR-182 expression in AH was 2-fold higher in HTG patients than nonglaucoma controls (P = 0.03) without controlling for medication treatment.ConclusionsOur integrative study is the first to associate rs76481776 with POAG via elevated miR-182 expression.

Primary open-angle glaucoma (POAG) is a leading cause of blindness worldwide. To identify new susceptibility loci, we performed meta-analysis on genome-wide association study (GWAS) results from eight independent studies from the United States (3,853 cases and 33,480 controls) and investigated the most significantly associated SNPs in two Australian studies (1,252 cases and 2,592 controls), three European studies (875 cases and 4,107 controls) and a Singaporean Chinese study (1,037 cases and 2,543 controls). A meta-analysis of the top SNPs identified three new associated loci: rs35934224[T] in TXNRD2 (odds ratio (OR) = 0.78, P = 4.05 × 10(-11)) encoding a mitochondrial protein required for redox homeostasis; rs7137828[T] in 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 new targets for preventative therapies.