Estrogen pathway polymorphisms in relation to primary open angle glaucoma: an analysis accounting for gender from the United States.

PURPOSE: Circulating estrogen levels are relevant in glaucoma phenotypic traits. We assessed the association between an estrogen metabolism single nucleotide polymorphism (SNP) panel in relation to primary open angle glaucoma (POAG), accounting for gender. METHODS: We included 3,108 POAG cases and 3,430 controls of both genders from the Glaucoma Genes and Environment (GLAUGEN) study and the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) consortium genotyped on the Illumina 660W-Quad platform. We assessed the relation between the SNP panels representative of estrogen metabolism and POAG using pathway- and gene-based approaches with the Pathway Analysis by Randomization Incorporating Structure (PARIS) software. PARIS executes a permutation algorithm to assess statistical significance relative to the pathways and genes of comparable genetic architecture. These analyses were performed using the meta-analyzed results from the GLAUGEN and NEIGHBOR data sets. We evaluated POAG overall as well as two subtypes of POAG defined as intraocular pressure (IOP) ≥22 mmHg (high-pressure glaucoma [HPG]) or IOP <22 mmHg (normal pressure glaucoma [NPG]) at diagnosis. We conducted these analyses for each gender separately and then jointly in men and women. RESULTS: Among women, the estrogen SNP pathway was associated with POAG overall (permuted p=0.006) and HPG (permuted p<0.001) but not NPG (permuted p=0.09). Interestingly, there was no relation between the estrogen SNP pathway and POAG when men were considered alone (permuted p>0.99). Among women, gene-based analyses revealed that the catechol-O-methyltransferase gene showed strong associations with HTG (permuted gene p≤0.001) and NPG (permuted gene p=0.01). CONCLUSIONS: The estrogen SNP pathway was associated with POAG among women.

Association of a Primary Open-Angle Glaucoma Genetic Risk Score With Earlier Age at Diagnosis.

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.

A Common Variant in MIR182 Is Associated With Primary Open-Angle Glaucoma in the NEIGHBORHOOD Consortium.

PURPOSE: Noncoding 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). METHODS: Using 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. RESULTS: Only 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. CONCLUSIONS: Our integrative study is the first to associate rs76481776 with POAG via elevated miR-182 expression.

DNA copy number variants of known glaucoma genes in relation to primary open-angle glaucoma.

PURPOSE: We examined the role of DNA copy number variants (CNVs) of known glaucoma genes in relation to primary open angle glaucoma (POAG). METHODS: Our study included DNA samples from two studies (NEIGHBOR and GLAUGEN). All the samples were genotyped with the Illumina Human660W_Quad_v1 BeadChip. After removing non-blood-derived and amplified DNA samples, we applied quality control steps based on the mean Log R Ratio and the mean B allele frequency. Subsequently, data from 3057 DNA samples (1599 cases and 1458 controls) were analyzed with PennCNV software. We defined CNVs as those ≥5 kilobases (kb) in size and interrogated by ≥5 consecutive probes. We further limited our investigation to CNVs in known POAG-related genes, including CDKN2B-AS1, TMCO1, SIX1/SIX6, CAV1/CAV2, the LRP12-ZFPM2 region, GAS7, ATOH7, FNDC3B, CYP1B1, MYOC, OPTN, WDR36, SRBD1, TBK1, and GALC. RESULTS: Genomic duplications of CDKN2B-AS1 and TMCO1 were each found in a single case. Two cases carried duplications in the GAS7 region. Genomic deletions of SIX6 and ATOH7 were each identified in one case. One case carried a TBK1 deletion and another case carried a TBK1 duplication. No controls had duplications or deletions in these six genes. A single control had a duplication in the MYOC region. Deletions of GALC were observed in five cases and two controls. CONCLUSIONS: The CNV analysis of a large set of cases and controls revealed the presence of rare CNVs in known POAG susceptibility genes. Our data suggest that these rare CNVs may contribute to POAG pathogenesis and merit functional evaluation.