Molecular Vision 2013; 19:1471-1481 <>
Received 22 April 2013 | Accepted 05 July 2013 | Published 12 July 2013

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

Louis R. Pasquale,1,2 Stephanie J. Loomis,1 Robert N. Weinreb,3 Jae H. Kang,2 Brian L. Yaspan,4 Jessica Cooke Bailey,5 Douglas Gaasterland,6 Terry Gaasterland,7 Richard K. Lee,8 William K. Scott,8 Paul R. Lichter,9 Donald L. Budenz,10 Yutao Liu,11,12 Tony Realini,13 David S. Friedman,14 Catherine A. McCarty,15 Sayoko E. Moroi,9 Lana Olson,5 Joel S. Schuman,16 Kuldev Singh,17 Douglas Vollrath,18 Gadi Wollstein,16 Donald J. Zack,14 Murray Brilliant,19 Arthur J. Sit,20 William G. Christen,21 John Fingert,22 Peter Kraft,23 Kang Zhang,3 R. Rand Allingham,11 Margaret A. Pericak-Vance,8 Julia E. Richards,9 Michael A. Hauser,11,12 Jonathan L. Haines,5 Janey L. Wiggs1

1Department of Ophthalmology, Mass Eye & Ear Infirmary, Harvard Medical School, Boston, MA; 2Channing Division of Network Medicine, Harvard Medical School, Boston, MA; 3Department of Ophthalmology and Hamilton Glaucoma Center, University of California, San Diego, SD, CA; 4Genentech Inc, San Francisco, CA; 5Center for Human Genetics Research, Vanderbilt University School of Medicine, Nashville TN; 6Eye Doctors of Washington, Chevy Chase, MD; 7Scripps Genome Center, University of California at San Diego, San Diego, CA;8Bascom Palmer Eye Institute and Hussman Institute for Human Genomics, University of Miami Miller School of Medicine, Miami, FL; 9Department of Ophthalmology and Visual Sciences, University of Michigan, Ann Arbor, MI;10Department of Ophthalmology, University of North Carolina, Chapel Hill, NC; 11Department of Ophthalmology, Duke University Medical Center, Durham, NC; 12Medicine, Duke University Medical Center, Durham, NC; 13Department of Ophthalmology, West Virginia University Eye Institute, Morgantown, WV;14Wilmer Eye Institute, Johns Hopkins University Hospital, Baltimore, MD; 15Essentia Institute of Rural Health, Duluth, MN; 16Department of Ophthalmology, UPMC Eye Center, University of Pittsburgh, Pittsburgh, PA;17Department of Ophthalmology, Stanford University, Palo Alto, CA; 18Department of Genetics, Stanford University, Palo Alto, CA; 19Center for Human Genetics, Marshfield Clinic Research Foundation, Marshfield, WI; 20Department of Ophthalmology, Mayo Clinic, Rochester, MN; 21Department of Medicine, Brigham and Women’s Hospital, Boston, MA; 22Departments of Ophthalmology and Anatomy/Cell Biology, University of Iowa, College of Medicine, Iowa City, IO;23Department of Biostatistics, Harvard School of Public Health, Boston, MA.

Correspondence to: Louis R. Pasquale, Massachusetts Eye and Ear Infirmary, 243 Charles Street, Boston, MA 02114; Phone: (617) 573-3674; FAX: (617) 573-3707; email:


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.


There is accumulating evidence that circulating estrogen levels influence processes relevant to glaucoma in the reproductive and post-reproductive years [1]. For example, elevated estrogen levels in pregnancy were associated with an approximate 10% reduction in intraocular pressure (IOP) [2-4], despite a concomitant increase in central corneal thickness [5]. Estrogen levels throughout the normal menstrual cycle influence the neuroretinal rim area measured with confocal scanning laser ophthalmoscopy [6] and mean sensitivity on short wavelength automated perimetry [7] in healthy women without ocular disease. Postmenopausal women have higher IOP than age-matched women who are still menstruating [2], and postmenopausal hormone (PMH) use is associated with reduced IOP in several studies [8-13]. Furthermore, ocular blood flow is reduced with aging, and PMH use seems to restore ocular blood flow resistance indices to levels achieved during the premenopausal period [9,14]. Early age at menopause was associated with increased risk of open angle glaucoma in the Rotterdam study [15]. In the Nurses Health Study (NHS), among women aged 65 or older, entering menopause at age 54 or later was associated with reduced risk of primary open angle glaucoma (POAG) [16]. Furthermore, the NHS revealed that PMH use was also associated with a reduced risk of high-tension POAG among women 65 years of age or older [16]. Collectively, these data suggest that estrogen levels in women may influence POAG, an age- and IOP-related optic neuropathy. Conversely, little is known about the role, if any, sex hormones play in POAG pathogenesis among men.

To date, genetic studies of estrogen’s role in POAG have been restricted to candidate gene studies. For example, given that estrogen receptor 2 (ESR2) is present in retinal ganglion cells of men and women [17], investigators have studied whether estrogen receptor 1 (ESR1) and ESR2 single nucleotide polymorphisms (SNPs) are associated with POAG with conflicting results [18,19]. No group has assessed whether an entire set of estrogen metabolic pathway SNPs is associated with POAG using a gender-specific approach. We hypothesize that if estrogen levels influence glaucoma risk, then a relationship may exist between the aggregate of estrogen pathway SNPs and POAG. This strategy allows for detection of a signal that arises from multiple individual weak genetic effects that impact the glaucomatous process. Researchers have studied the aggregate effects of genetic variants in the estrogen metabolic pathway for association with breast [20] and uterine cancer [21], and these approaches have yielded encouraging results. We assembled high throughput genotyping SNP data from two large data sets—the Glaucoma Genes and Environment (GLAUGEN) study and the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) consortium—into an estrogen metabolic SNP panel and used the Pathway Analysis by Randomization Incorporating Structure (PARIS) analysis software package [22] to test if this panel is associated with POAG. Similarly, we assessed whether collections of SNPs within genes that comprise the pathway are associated with POAG. Since the association between ESR2 SNPs and POAG [18,19] and other complex diseases [23,24] may differ between men and women, we conducted our analyses with consideration of such possible gender differences.


Description of the study populations

The GLAUGEN study is part of a network of collaborative studies [25] that has contributed to an improved understanding of genotype-phenotype associations [26] and gene–environment interactions for various complex traits [27], including POAG [28,29]. The GLAUGEN study contains POAG cases and controls derived from the Genetic Etiologies of Primary Open-Angle Glaucoma (GEP) study, the Nurses’ Health Study (NHS), and the Health Professionals Follow-up Study (HPFS). The GEP is a clinical case-control study from Massachusetts Eye and Ear Infirmary while the NHS and HPFS contain cases and controls nested within population-based cohorts at risk for POAG. Details regarding the inclusion and exclusion criteria for the GLAUGEN POAG case-control group have been described at the database of Genotypes and Phenotypes (dbGaP), accession number phs000308.v1.p1, and by Wiggs et al. [30]. The institutional review boards of Massachusetts Eye and Ear Infirmary, Brigham and Women’s Hospital, and Harvard School of Public Health approved this study.

The NEIGHBOR consortium represents the collaborative efforts of 12 institutions where POAG cases and controls were collected using harmonized criteria. Information pertaining to the study sites, design features, inclusion criteria, and clinical variables collected in the NEIGHBOR consortium has been described [31]. Data from the NEIGHBOR consortium have contributed to the understanding of the genetic structure of POAG [29] and central corneal thickness [32,33]. The institutional review boards of the participating institutions also approved this study.

Ophthalmic characteristics of study participants

For the patients, slit-lamp biomicroscopy did not reveal exfoliation syndrome, pigment dispersion syndrome, or other findings that could produce elevated IOP. The iridocorneal filtration angle was judged to be open in both eyes. At least one eye had a cup-disc ratio (CDR) >0.7 and visual field (VF) loss or the VF loss was replicated on a subsequent test regardless of CDR. Eyes with VF loss had defects localized to the nerve fiber layer on reliable tests (fixation loss was ≤33%, false positive rate was ≤20%, and false negative rate was ≤20%). The type of VF test used was not specified. IOP at diagnosis was collected and used to categorize the POAG cases into high-pressure glaucoma (HPG with IOP ≥22 mmHg) and normal pressure glaucoma (NPG with IOP <22 mmHg) subtypes. Controls underwent an eye exam that revealed an IOP <22 mmHg and CDR≤0.6 in both eyes. For members of the NHS and HPFS, the eye exam was performed by the participants’ local eye care provider, and that exam reportedly did not reveal signs of glaucoma.

Genotyping data

The methods used to collect, extract, and plate DNA for GLAUGEN [30] and NEIGHBOR [29] have been previously outlined. Genotyping for GLAUGEN and NEIGHBOR was performed at the Broad Institute (Cambridge, MA) and the Center for Inherited Disease Research (Baltimore, MD), respectively. We used the Illumina Human660W-Quad-v1 array (Illumina; San Diego, CA) to genotype subjects for both studies. The algorithm for calling genotypes, the SNP quality control (QC) filters used, and the preliminary analyses undertaken to discover the predictors of genotype calling rates have been discussed in prior publications [29,30].

Gene association analyses

We used PLINK v1.07 [34] to perform gene association analysis for POAG in GLAUGEN and NEIGHBOR for SNPs that passed the QC process (484,419 SNPs in GLAUGEN and 521,683 SNPs in NEIGHBOR) as previously described [29,30]. Briefly, these analyses were conducted among self-reported Caucasians and adjusted for age, gender, study site, and population structure in both study groups. In GLAUGEN, additional adjustments were made for the DNA extraction method (DNAzol, Beverly, MA; Qiagen or GENTRA, Germatown, MD) and the DNA source (blood or cheek cell sample) because these parameters influenced the genotyping call rates [29,30]. Using the METAL software package [35], we then performed a meta-analysis of the GLAUGEN and NEIGHBOR data sets and obtained p values and odds ratios for the 480,335 overlapping SNPs in the combined data set that passed QC. In this paper, we used the p values and odds ratio for the estrogen metabolism SNPs in the pathway analysis discussed below.

Defining the estrogen pathway and Pathway Analysis by Randomization Incorporating Structure pathway analysis

We generated a custom list of 903 SNPs in 23 genes across 15 chromosomes comprising the estrogen metabolic pathway using the Kyoto Encyclopedia of Genes and Genomes (KEGG) [36] online database and other academic sources [37-40] (see Figure 1). In analyses in which men and women were considered together or men were considered alone, X chromosome variants involved in estrogen metabolism were excluded, yielding a total of 874 SNPs available for study. We uploaded the meta-analyzed p values for genotyped SNPs within these genes ± a 50kB genomic window to PARIS. The statistical approach used by PARIS has been previously described [22]. Conceptually, PARIS first analyzes the genomic features of the pathway, namely, the number of linkage disequilibrium (LD) blocks, the number of SNPs per block, and the number of SNPs not in any LD block for the gene biomarker assembly. The estrogen pathway contains 122 complex features (number of LD blocks with two or more types SNPs) and 123 simple features (number of SNPs not in any LD block). After this step, PARIS compares the overall association between the estrogen pathway and POAG versus the associations between 1,000 randomly generated pathways with similar genetic architecture and POAG. PARIS assesses pathways of interest using a permuted p value threshold of <0.05, which is considered statistically significant. For example, with the estrogen pathway among women, the results indicated 39 significant features; specifically, 15 of 123 simple features (or SNPs) had p value<0.05 for association with POAG, and 24 of 122 complex features were statistically significant in relation to POAG (a complex feature was significant if any SNP in the LD block had a p value<0.05). The PARIS-generated permuted p value=0.006 for the estrogen pathway among women indicated that six of 1,000 random pathways with similar genetic architecture had a higher significant feature count (>39 significant features with p<0.05) in relation to POAG. These analyses were performed for POAG in a gender-specific manner. Subsequently, we considered gender-specific HPG and NPG as outcomes. Finally, we looked at these outcomes when men and women were considered together. To determine which of the genes or SNPs in the estrogen metabolism pathway contributed to any significant signal in the pathway overall, we used the “-I” (Investigate) option within PARIS, which reports the p values of genes and SNPs within the pathway.


This study sample included 3,108 cases (976 from GLAUGEN and 2,132 from NEIGHBOR) and 3,430 controls (1,140 from GLAUGEN and 2,290 from NEIGHBOR; Table 1). 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) among women. Among men, the relation between the estrogen SNP pathway and POAG (permuted p>0.99), HPG (permuted p=0.92), and NPG (permuted p=0.59) was consistently null across the IOP spectrum. The estrogen pathway was not significantly associated with POAG overall (permuted p=0.72), HPG (permuted p=0.12), and NPG (permuted p=0.06) when men and women were considered together.

Analogous to the pathway approach, we assessed the aggregate of SNPs within specific estrogen pathway genes for association with POAG. We performed this analysis to gain insight into which aspects of the pathway contributed to the association we observed among women. ESR1, aldo-keto reductase family 1, member D1 (AKR1D1), cytochrome P450, family 3, subfamily A, polypeptide 4 (CYP3A4), catechol-O-methyltransferase (COMT), and steroid sulfatase (microsomal), isozyme S (STS) were significantly (permuted gene p values of <0.001, 0.05, 0.01, <0.001, and <0.001, respectively) associated with POAG overall among women (Table 2). With the exception of AKR1D1, these genes were also associated with HPG among women (all with permuted gene p<0.001), along with aldo-keto reductase family 1, member C4 (AKR1C4; permuted gene p=0.002) and hydroxysteroid (17-beta) dehydrogenase 1 (HSD17B1; permuted p<0.001). Only COMT and STS were associated with HPG (permuted p<0.001 for both genes) and NPG (permuted gene p=0.01 and p<0.001, respectively); otherwise, there was little overlap between the genes significantly associated with HPG and NPG. Other genes associated with NPG among women included AKR1D1 (permuted gene p=0.046), CYP17A1 (permuted gene p=0.04), ESR2 (permuted gene p<0.001), CYP1A1 (permuted gene p<0.001), and CYP1A2 (permuted gene p=0.004).

Although the COMT SNP set showed a significant association with HPG and NPG among women, the SNP set was not associated with the POAG or POAG subtypes among men (permuted gene p≥0.66). Among men, gene-based analysis revealed one gene associated with POAG overall and with HPG: CYP1B1 (permuted gene p≤0.02). The only other gene associated with HPG among men was CYP11A1 (permuted gene p=0.02). Furthermore, the following estrogen metabolism genes were associated with NPG among men: HSD11B1 (permuted gene p=0.003), CYP1A1 (permuted gene p=0.001), and CYP1A2 (permuted gene p<0.001).

To refine the source of significant association between estrogen metabolic pathways and POAG even further among women, we report the nominally significant SNPs that were in any relevant genes in POAG overall, HPG or NPG (Appendix 1). The p values for all SNPs in the estrogen metabolic pathway are provided in Appendix 2. As with the genes overall, there was some overlap between the significant SNPs in POAG overall and HPG, but little overlap between the HPG and NPG subsets. For example, in COMT, one of two genes that were significant in HPG and NPG, only one significant SNP was shared between the HPG and NPG subsets (rs3804047). Although one COMT SNP (rs2531697) showed particularly strong p values for association with HPG in women (p=4.10E-05), the SNP would not remain significant if one accounted for the multiple comparisons made in this analysis (903 SNPs stratified by gender and IOP yielding a corrected significance level of 1.85E-05); however, our pathway analysis points to a large number of COMT SNPs significant for POAG and its subtypes among women.


The relationship between environmental factors such as PMH use and POAG among women has been previously described [15,16,41]. Candidate gene studies assessing selected estrogen metabolic SNPs [19,42] as well as candidate gene–environment interactions [43,44] in POAG among women have also been reported. This is the first paper to assess a possible association between the estrogen metabolic SNP pathway as a whole and POAG using genome-wide association data. Our results suggest that the estrogen metabolic pathway is associated with POAG among women supporting the notion that multiple genetic signals of modest effect contribute to the glaucomatous process. Furthermore, the association between the estrogen metabolic pathway and POAG among women was more significant in those classified as having HPG. The results highlight a distinct sexual dimorphism in that the pathway is definitively not associated with POAG among men. This latter result also contributes to a null result when the relation between the estrogen pathway and POAG is considered in men and women jointly.

ESR1 codes for estrogen receptor alpha and is expressed in the retina of both sexes [17]. In a gene association study of 87 incident POAG cases and 3,616 controls, a haplotype of two ESR1 SNPs was not associated with open angle glaucoma in women or men in the Rotterdam study [19]. Our study with 3,108 cases and 3,430 controls contains genotypes on 127 ESR1 SNPs. Collectively, ESR1 SNPs contributed to the association between the estrogen pathways and POAG overall and in HPG in women (see Table 2; gene permuted p<0.001). The signal from ESR1 was modest as the ESR1 gene permuted p value for association with POAG among women was 0.31 in NEIGHBOR alone (2132 cases) and 0.26 in GLAUGEN alone (976 cases), only reaching significance in the combined data set. Furthermore, ESR1 was modestly associated with POAG overall when men and women were considered together (gene permuted p=0.023).

Exogenous estrogen is protective in various animal models of neurodegeneration [45,46] including experimental models of glaucoma [47,48]. This raises questions regarding whether polymorphisms in ESR2 (which codes for the estrogen receptor found on retinal ganglion cells [17]) may alter affinity for estrogen and affect retinal ganglion cell viability in glaucoma. A Japanese group found an association between ESR2rs1256031 and rs4986938 variants and HPG in a study including 212 female cases and 191 controls [42]. Our study included 859 women with HPG and 1,937 female controls and contained genotypes on 38 ESR2 SNPs including several SNPs in high LD with rs1256031 and rs4986938. Collectively, we found that ESR2 SNPs were associated with NPG (gene permuted p<0.001) but not HPG (gene permuted p=0.10) among women. Furthermore, the rs1256031 and rs4986938 variants were specifically not associated with HPG among the women in our data set. The Rotterdam study of European Caucasians did not find associations between these two ESR2 SNPs and open angle glaucoma in women but found an association in men. We did not find an association between the ESR2 collection of SNPs and POAG in men (gene permuted p>0.99) or for the rs1256031 and rs4986938 SNPs in particular (p≥0.45). The conflicting results may be due to differences in sample size and population stratification. Assessing large data sets and evaluating the collective effects of genes in a pathway and the pathway as a whole in relation to glaucoma provide greater power to find existing associations. Overall, our data set suggests that estrogen receptor polymorphisms are important in POAG among women, with ESR1 SNPs playing an important role in HPG and ESR2 SNPS playing a role in NPG.

COMT codes for catechol-O-methyltransferase that catalyzes the transfer of a methyl group from S-adenosyl methionine to a hydroxyl group on estrogen, dopamine, and epinephrine. The enzyme COMT is involved in methylating two derivatives of estradiol—2-hydroxyestradiol and 4-hydroxyestradiol—and perhaps reducing estrogen bioavailability. Variants in COMT were associated with HPG and NPG in women but not men. Although no other group has reported on associations between COMT variants and open angle glaucoma, we saw some evidence of internal replication within our data sets. Among women, the permuted p value for association between COMT and HPG was 0.01 in NEIGHBOR and 0.02 in GLAUGEN. Interestingly, there is robust evidence that gene variants in COMT differentially impact brain function in men and women [20,49-55].

Cytochrome P450 enzymes have an important role in estrogen metabolism, and several genes coding for these liver enzymes contributed to the association between the estrogen pathway and POAG in women. CYP3A4 showed an association with HTG in women while CYP17A1, CYPIA1, and CYPIA2 showed an association with NPG among women. Selected CYP1B1 SNPs have been studied in association with POAG. A prior meta-analysis of the CYP1BI SNPs in association with POAG consisting of 1,953 cases and 1,341 controls that found no significant association [56] is largely consistent with our larger study (3,108 cases and 3,430 controls), which found no association between the CYP1B1 aggregate of SNPs and POAG among men and women together (permuted gene p=0.098). The role of cytochrome P450 metabolism in glaucoma is unknown, but it is intriguing to speculate that metabolism of intracellular estrogen is involved. Members of the cytochrome P450 family could alter estrogen metabolism over time [57] and impact the development of open angle glaucoma.

This study has limitations. First, none of the individual SNPs in the estrogen metabolic pathway were significantly associated with POAG either overall or in either gender if a strict Bonferroni correction is imposed. Nonetheless, existing genome-wide studies of POAG are underpowered to find biologically meaningful SNPs of modest effects due to the need to curb the false discovery rate. The purpose of this analysis is to use a hypothesis-driven approach that involves SNP subsets to find potentially important associations that provide insight into POAG pathogenesis. Second, no other current data set of comparable size can be used to replicate the findings we present here; however, we demonstrated internal replication between NEIGHBOR and GLAUGEN for our one novel finding, the association between COMT and HPG in women. Therefore, even though it is intriguing that common polymorphisms in a network of estrogen-metabolizing genes may be relevant to POAG development among women, caution is warranted until these findings can be replicated in other large data sets. Finally, although the estrogen pathway SNPs were associated with POAG overall in women, it is not possible to identify functional networks with discrete effect size on the glaucomatous process. Such analysis is the focus of future research.

This is the largest study (6,538 subjects) assessing the relation between aggregate effects of 903 estrogen-metabolizing SNPs in relation to POAG. We find that this SNP collection is associated with POAG in women but not men. These data may shed light on the many epidemiological data and basic science research implicating estrogen levels in POAG pathogenesis among women. The work may open new avenues of research into the role of estrogen in the development of POAG and may lead to new gender-specific prevention strategies that are tailored to an individual’s genetic makeup.

Appendix 1.

Appendix 2.


Funding/Support: A Horizon Grant to MEEI from Allergan (Irvine, CA) supported the collection of some glaucoma feature data. The Harvard Glaucoma Center of Excellence and Margolis fund (Boston, MA) support LRP and JLW. LRP, JER and JLW are also supported by Research to Prevent Blindness, Inc. (New York, NY). The Arthur Ashley Foundation also supports Dr. Pasquale. The Glaucoma Research Foundation (San Francisco, CA), American Health Assistance Foundation (Clarksburg, MD), and the Glaucoma Foundation (New York, NY) support YL. The following National Institutes of Health grants support the maintenance of the Nurses Health Study and Health Professionals Follow-up, allowing these health professionals to contribute to this analysis: CA87969, CA49449, UM1 CA167552, and HL35464. The following grants from the National Human Genome Research Institute (Bethesda, MD) supported GLAUGEN: HG004728 (LRP), HG004424 (Broad Institute to support genotyping), HG004446 (C. Laurie, U. Washington, to support genotype data cleaning and analysis). Genotyping services for the NEIGHBOR study were provided by the Center for Inherited Disease Research (CIDR) and were supported by the National Eye Institute through grant HG005259-01 (JLW). Additionally CIDR is funded through a federal contract from the National Institutes of Health to The Johns Hopkins University, contract number HHSN268200782096C. The National Eye Institute (Bethesda, MD) through ARRA grants EY015872 (JLW) and EY019126 (MAH) supported the collection and processing of samples for the NEIGHBOR dataset. Funding for the collection of cases and controls was provided by National Institutes of Health (Bethesda, MD) grants: EY015543 (RRA), EY006827 (DG), HL073389 (Hauser, E); EY13315 (MAH); EY09611 (Hankinson, S), and EY015473 (LRP). EY009149 (PRL), HG004608 (CAM), EY008208 (Medeiros, P.), EY015473 (LRP), EY012118 (MAP-V), EY015682 (TR), EY011671 (JER), EY09580 (JER), EY013178 (JSS), EY015872 (JLW), EY010886 (JLW), EY009847 (JLW), EY011008 (Zangwill, L), EY144428 (KZ), EY144448 (KZ), EY18660 (KZ). None of the authors have any commercial interests in the subject of the manuscript or in entities discussed in the manuscript. This work was presented in part at the American Society of Human Genetics Meeting, November 2012, San Francisco, CA and at the ARVO Meeting May 2013, Seattle Washington.


  1. Vajaranant TS, Pasquale LR. Estrogen deficiency accelerates aging of the optic nerve. Menopause. 2012; 19:942-7. [PMID: 22415565]
  2. Qureshi IA. Intraocular pressure: association with menstrual cycle, pregnancy and menopause in apparently healthy women. Chin J Physiol. 1995; 38:229-34. [PMID: 8925675]
  3. Phillips CI, Gore SM. Ocular hypotensive effect of late pregnancy with and without high blood pressure. Br J Ophthalmol. 1985; 69:117-9. [PMID: 3966998]
  4. Green K, Phillips CI, Cheeks L, Siagle T. Aqueous humor flow rate and intraocular pressure during and after pregnancy. Ophthalmic Res. 1988; 20:353-7. [PMID: 3237393]
  5. Weinreb RN, Lu A, Beeson C. Maternal corneal thickness during pregnancy. Am J Ophthalmol. 1988; 105:258-60. [PMID: 3344782]
  6. Akar ME, Taskin O, Yucel I, Akar Y. The effect of the menstrual cycle on optic nerve head analysis in healthy women. Acta Ophthalmol Scand. 2004; 82:741-5. [PMID: 15606474]
  7. Yucel I, Akar ME, Dora B, Akar Y, Taskin O, Ozer HO. Effect of the menstrual cycle on standard achromatic and blue-on-yellow visual field analysis of women with migraine. Can J Ophthalmol. 2005; 40:51-7. [PMID: 15825530]
  8. Affinito P, Di Spiezio Sardo A, Di Carlo C, Sammartino A, Tommaselli GA, Bifulco G, Loffredo A, Loffredo M, Nappi C. Effects of hormone replacement therapy on ocular function in postmenopause. Menopause. 2003; 10:482-7. [PMID: 14501611]
  9. Altintaş O, Caglar Y, Yuksel N, Demirci A, Karabas L. The effects of menopause and hormone replacement therapy on quality and quantity of tear, intraocular pressure and ocular blood flow. Ophthalmologica. 2004; 218:120-9. [PMID: 15004502]
  10. Sator MO, Joura EA, Frigo P, Kurz C, Metka M, Hommer A, Huber JC. Hormone replacement therapy and intraocular pressure. Maturitas. 1997; 28:55-8. [PMID: 9391995]
  11. Sorrentino C, Affinito P, Mattace Raso F, Loffredo M, Merlino P, Loffredo A, Palomba S, Nappi C. Effect of hormone replacement therapy on postmenopausal ocular function. Minerva Ginecol. 1998; 50:19-24. [PMID: 9577151]
  12. Tint NL, Alexander P, Tint KM, Vasileiadis GT, Yeung AM, Azuara-Blanco A. Hormone therapy and intraocular pressure in nonglaucomatous eyes. Menopause. 2010; 17:157-60. [PMID: 19770781]
  13. Uncu G, Avci R, Uncu Y, Kaymaz C, Develioglu O. The effects of different hormone replacement therapy regimens on tear function, intraocular pressure and lens opacity. Gynecol Endocrinol. 2006; 22:501-5. [PMID: 17071534]
  14. Harris-Yitzhak M, Harris A, Ben-Refael Z, Zarfati D, Garzozi HJ, Martin BJ. Estrogen-replacement therapy: effects on retrobulbar hemodynamics. Am J Ophthalmol. 2000; 129:623-8. [PMID: 10844054]
  15. Hulsman CA, Westendorp IC, Ramrattan RS, Wolfs RC, Witteman JC, Vingerling JR, Hofman A, de Jong PT. Is open-angle glaucoma associated with early menopause? The Rotterdam Study. Am J Epidemiol. 2001; 154:138-44. [PMID: 11447046]
  16. Pasquale LR, Rosner BA, Hankinson SE, Kang JH. Attributes of female reproductive aging and their relation to primary open-angle glaucoma: a prospective study. J Glaucoma. 2007; 16:598-605. [PMID: 18091177]
  17. Munaut C, Lambert V, Nöel A, Frankenne F, Deprez M, Foidart JM, Rakic JM. Presence of oestrogen receptor type beta in human retina. Br J Ophthalmol. 2001; 85:877-82. [PMID: 11423466]
  18. Mabuchi F, Sakurada Y, Kashiwagi K, Yamagata Z. IIjuima H, Tsukahara S. Estrogen receptor beta gene polymorphisms and intraocular pressure elevation in female patients with primary open-angle glaucoma. Am J Ophthalmol. 2010; 149:826-30. [PMID: 20399928]
  19. de Voogd S, Wolfs RC, Jansonius NM, Uitterlinden AG, Pols HA, Hofman A, de Jong PT. Estrogen receptors alpha and beta and the risk of open-angle glaucoma: the Rotterdam Study. Arch Ophthalmol. 2008; 126:110-4. [PMID: 18195227]
  20. Cerne JZ, Pohar-Perme M, Novakovic S, Frkovic-Grazio S, Stegel V, Gersak K. Combined effect of CYP1B1, COMT, GSTP1, and MnSOD genotypes and risk of postmenopausal breast cancer. J Gynecol Oncol. 2011; 22:110-9. [PMID: 21860737]
  21. Yang HP, Gonzalez Bosquet J, Li Q, Platz EA, Brinton LA, Sherman ME, Lacey JV, , Jr Gaudet MM, Burdettte LA, Figueroa JD, Ciampa JG, Lissowska J, Peplonska B, Chanock SJ, Garcia-Closas M. Common genetic variation in the sex hormone metabolic pathway and endometrial cancer risk: pathway-based evaluation of candidate genes. Carcinogenesis. 2010; 31:827-33. [PMID: 20053928]
  22. Yaspan BL, Bush WS, Torstenson ES, Ma D, Pericak-Vance MA, Ritchie MD, Sutcliffe JS, Haines JL. Genetic analysis of biological pathway data through genomic randomization. Hum Genet. 2011; 129:563-71. [PMID: 21279722]
  23. Rexrode KM, Ridker PM, Hegener HH, Buring JE, Manson JE, Zee RY. Polymorphisms and haplotypes of the estrogen receptor-beta gene (ESR2) and cardiovascular disease in men and women. Clin Chem. 2007; 53:1749-56. [PMID: 17702854]
  24. Pirskanen M, Hiltunen M, Mannermaa A, Helisalmi S, Lehtovirta M, Hänninen T, Soinien H. Estrogen receptor beta gene variants are associated with increased risk of Alzheimer's disease in women. Eur J Hum Genet. 2005; 13:1000-6. [PMID: 15944651]
  25. Cornelis MC, Agrawal A, Cole JW, Hansel NN, Barnes KC, Beaty TH, Bennett SN, Bierut LJ, Boerwinkle E, Doheny KF, Feenstra B, Feingold E, Fornage M, Haiman CA, Harris EL, Hayes MG, Heit JA, Hu FB, Kang JH, Laurie CC, Ling H, Manolio TA, Marazita ML, Mathias RA, Mirel DB, Paschall J, Pasquale LR, Pugh EW, Rice JP, Udren J, van Dam RM, Wang X, Wiggs JL, Williams K, Yu K, GENEVA consortium. The Gene, Environment Association Studies consortium (GENEVA): maximizing the knowledge obtained from GWAS by collaboration across studies of multiple conditions. Genet Epidemiol. 2010; 34:364-72. [PMID: 20091798]
  26. Chu AY, Workalemahu T, Paynter NP, Rose LM, Giulianni F, Tanaka T, Ngwa JS, CHARGE Nutrition Working Group. QiQ, Curhan GC, Rimm EB, Hunter DJ, Pasquale LR, Ridker PM, Hu FB, Chasma DI, Qi L, DietGen Consortium. Novel locus including FGF21 is associated with dietary macronutrient intake. Hum Mol Genet. 2013; 22:1895-1902. [PMID: 23372041]
  27. Qi Q, Chu AY, Kang JH, Jensen MK, Curhan GC, Pasquale LR, Riker PM, Hunter DJ, Willett WC, Rimm EB, Chasman DI, Hu FB, Qi L. Sugar-sweetened beverages and genetic risk of obesity. N Engl J Med. 2012; 367:1387-96. [PMID: 22998338]
  28. Pasquale LR, Loomis SJ, Kang JH, Yaspan BL, Abdrabou W, Budenz DL, Chen TC, Delbono E, Friedman DS, Gaasterland D, Gaasterland T, Grosskreutz CL, Lee RK, Lichter PR, Liu Y, McCarty CA, Moroi SE, Olson LM, Realini T, Rhee DJ, Schuman JS, Singh K, Vollrath D, Wollstein G, Zack DJ, Allingham RR, Pericak-Vance MA, Weinreb RN, Zhang K, Hauser MA, Richards JE, Haines JL, Wiggs JL. CDKN2B-AS1 genotype-glaucoma feature correlations in primary open-angle glaucoma patients from the United States. Am J Ophthalmol. 2013; 155:342-53. [PMID: 23111177]
  29. Wiggs JL, Yaspan B, Hauser M, Kang JH, Allingham RR, Olson LM, Abdrabou W, Fan BJ, Wang DY, Brodeur W, Budenz DL, Caprioli J, Crenshaw A, Crooks K, Delbono E, Doheny KF, Friedman DS, Gaasterland D, Gaasterland T, Laurie C, Lee RK, Lichter PR, Loomis S, Liu Y, Medeiros FA, McCarty C, Mirel D, Moroi SE, Musch DC, Realini A, Rezsa FW, Schuman JS, Scott K, Singh K, Stein JD, Trager EH, Vanveldhuisen P, Vollrath D, Wollstein G, Yoneyama S, Zhang K, Weinreb RN, Ernst J, Kellis M, Masuda T, Zack D, Richards JE, Pericak-Vance M, Pasquale LR, Haines JL. Common variants in the CDKN2BAS region on 9p21 and a novel region on 8q23 are associated with the normal tension form of open-angle glaucoma with increased susceptibility to optic nerve degeneration. PLoS Genet. 2012; 8:e1002654 [PMID: 22570617]
  30. Wiggs JL, Kang JH, Yaspan BL, Mirel DB, Laurie C, Crenshaw A, Brodeur W, Gogarten S, Olson LM, Abdrabou W, Delbono E, Loomis S, Haines JL, Pasquale LR. GENEVA Consortium. Common variants near CAV1 and CAV2 are associated with primary open-angle glaucoma in Caucasians from the USA. Hum Mol Genet. 2011; 20:4707-13. [PMID: 21873608]
  31. Wiggs JL, Hauser M, Abdrabou W, Allingham RR, Budenz DL, Delbono E, Friedman DS, Kang JH, Gaasterland D, Gaasterland T, Lee RK, Lichter PR, Loomis S, Liu Y, McCarty C, Medeiros FA, Moroi SE, Olson LM, Realini A, Richards JE, Rozsa FW, Schuman JS, Singh K, Stein JD, Vollrath D, Weinreb RN, Wollstein G, Yaspan BL, Yoneyama S, Zack D, Zhang K, Pericak-Vance M, Pasquale LR, Haines JL. The NEIGHBOR Consortium Primary Open Angle Glaucoma Genome-wide Association Study: Rationale, study design and clinical variables. J Glaucoma. 2012;
    Epub ahead of print
    [PMID: 22828004]
  32. Lu Y, Vitart V, Burdon KP, Khor CC, Bykhovskaya Y, Mirshahi A, Hewitt AW, Koehn D, Hysi PG, Ramdas WD, Zeller T, Vithana EN, Cornes BK, Tay WT, Tai ES, Cheng CY, Liu J, Foo JN, Saw SM, Thorleifsson G, Stefansson K, Dimasi DP, Mills RA, Mountain J, Ang W, Hoehn R, Verhoeven VJ, Grus F, Wolfs R, Castagne R, Lackner KJ, Springelkamp H, Yang J, Jonasson F, Leung DY, Chen LJ, Tham CC, Rudan I, Vatavuk Z, Hayward C, Gibson J, Cree AJ, MacLeod A, Ennis S, Polasek O, Campbell H, Wilson JF, Ulmer M, Li J, Yaspan BL, Ozel AB, Richards JE, Moroi SE, Haines JL, Kang JH, Pasquale LR, Allingham RR, Ashley-Koch A, NEIGHBOR Consortium. Mitchell P, Wang JJ, Wright AF, Pennell C, Spector TD, Young TL, Klaver CC, Martin NG, Montgomery GW, Anderson MG, Aung T, Willoughby CE, Wiggs JL, Pang CP, Thorsteinsdottir U, Lotery AJ, Hammond CJ, van Diujn CM, Hauser MA, Rabinowitz YS, Pfeiffer N, Mackey DA, Criag JE, Macgregor S, Wong TY. Genome-wide association analyses identify multiple loci associated with central corneal thickness and keratoconus. Nat Genet. 2013; 45:155-63. [PMID: 23291589]
  33. Ulmer M, Li J, Yaspan BL, Richards JE, Moroi SE, Hawthrone F, Budenz DL, Friedman DS, Gaasterland D, Haines J, Kang JH, Lee R, Lichter P, Liu Y, Pasquale LR, Pericak-Vance M, Realini A, Schuman JS, Singh K, Vollrath D, Weinreb R, Wollstein G, Zack DJ, Zhang K, Young T, Allingham RR, Wiggs JL, Ashley-Koch A, Hauser MA. Genome-wide analysis of central corneal thickness in primary open-angle glaucoma cases in the NEIGHBOR and GLAUGEN consortia. Invest Ophthalmol Vis Sci. 2012; 53:4468-74. [PMID: 22661486]
  34. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007; 81:559-75. [PMID: 17701901]
  35. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010; 26:2190-1. [PMID: 20616382]
  36. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000; 28:27-30. [PMID: 10592173]
  37. Hefler LA, Grimm C, Heinze G, Schneeberger C, Mueller MW, Muendlein A, Huber JC, Leodolter S, Tempfer CB. Estrogen-metabolizing gene polymorphisms and age at natural menopause in Caucasian women. Hum Reprod. 2005; 20:1422-7. [PMID: 15774541]
  38. Cooper GS, Miller FW, Pandey JP. The role of genetic factors in autoimmune disease: implications for environmental research. Environ Health Perspect. 1999; 107Suppl 5:693-700. [PMID: 10502533]
  39. Lakhani NJ, Sarkar MA, Venitz J, Figg WD. 2-Methoxyestradiol, a promising anticancer agent. Pharmacotherapy. 2003; 23:165-72. [PMID: 12587805]
  40. Wishart DS, Jewison T, Guo AC, Wilson M, Knox C, Liu Y, Djoumbou Y, Mandel R, Aziat F, Dong E, Bouatra S, Sinelnikov I, Arndt D, Xia J, Liu P, Yallou F, Bjorndahl T, Perez-Pineiro R, Eisner R, Allen F, Neveu V, Greiner R. Scalbert. HMDB 3.0 - - The Human Metabolome Database in 2013. Nucleic Acids Res. 2013; 41Database issue:D801-7. [PMID: 23161693]
  41. Lee AJ, Mitchell P, Rochtchina E, Healey PR. Female reproductive factors and open angle glaucoma: the Blue Mountains Eye Study. Br J Ophthalmol. 2003; 87:1324-8. [PMID: 14609824]
  42. Mabuchi F, Sakurada Y, Kashiwagi K, Yamagata Z, Iijima H, Tsukahara S. Estrogen receptor beta gene polymorphism and intraocular pressure elevation in female patients with primary open-angle glaucoma. Am J Ophthalmol. 2010; 149:826-30. [PMID: 20399928]
  43. Kang JH, Wiggs JL, Rosner BA, Hankinson SE, Abdrabou W, Fan BJ, Haines J, Pasquale LR. Endothelial nitric oxide synthase gene variants and primary open-angle glaucoma: interactions with sex and postmenopausal hormone use. Invest Ophthalmol Vis Sci. 2010; 51:971-9. [PMID: 19815736]
  44. Kang JH, Wiggs JL, Abdrabou W, Haines J, Pasquale LR. Reproductive Factors and NOS3 Variant Interactions in Primary Open-Angle Glaucoma. Mol Vis. 2011; 17:2544-51. [PMID: 22025889]
  45. Morissette M, Al Sweidi S, Callier S, DiPaolo T. Estrogen and SERM neuroprotection in animal models of Parkinson's disease. Mol Cell Endocrinol. 2008; 290:60-9. [PMID: 18515001]
  46. Carroll JC, Rosario ER, Villamagna A, Pike CJ. Continuous and cyclic progesterone differentially interact with estradiol in the regulation of Alzheimer-like pathology in female 3xTransgenic-Alzheimer's disease mice. Endocrinology. 2010; 151:2713-22. [PMID: 20410196]
  47. Zhou X, Li F, Ge J, Sarkisian SR, , Jr Tomita H, Zaharia A, Chodosh J, Cao W. Retinal ganglion cell protection by 17-beta-estradiol in a mouse model of inherited glaucoma. Dev Neurobiol. 2007; 67:603-16. [PMID: 17443811]
  48. Russo R, Cavaliere F, Watanabe C, Nucci C, Bagetta G, Corasaniti MT, Sakurada S, Morrone LA. 17Beta-estradiol prevents retinal ganglion cell loss induced by acute rise of intraocular pressure in rat. Prog Brain Res. 2008; 173:583-90. [PMID: 18929136]
  49. Domschke K, Baune BT, Havlik L, Stuhrmann A, Suslow T, Kugel H, Zwanzger P, Grotegerd D, Sehlmeyer C, Arolt V, Dannlowski U. Catechol-O-methyltransferase gene variation: impact on amygdala response to aversive stimuli. Neuroimage. 2012; 60:2222-9. [PMID: 22387174]
  50. Domschke K, Deckert J, O'Donovan MC, Glatt SJ. Meta-analysis of COMT val158met in panic disorder: ethnic heterogeneity and gender specificity. Am J Med Genet B Neuropsychiatr Genet. 2007; 144B:667-73. [PMID: 17357147]
  51. Beuten J, Payne TJ, Ma JZ, Li MD. Significant association of catechol-O-methyltransferase (COMT) haplotypes with nicotine dependence in male and female smokers of two ethnic populations. Neuropsychopharmacology. 2006; 31:675-84. [PMID: 16395295]
  52. Chen C, Chen C, Moyzis R, Dong Q, He Q, Zhu B, Li J, Li H, Li J, Lessard J. Sex modulates the associations between the COMT gene and personality traits. Neuropsychopharmacology. 2011; 36:1593-8. [PMID: 21471954]
  53. Enoch MA, Waheed JF, Harris CR, Albaugh B, Goldman D. Sex differences in the influence of COMT Val158Met on alcoholism and smoking in plains American Indians. Alcohol Clin Exp Res. 2006; 30:399-406. [PMID: 16499480]
  54. Holtzer R, Ozelius L, Xue X, Wang T, Lipton RB, Verghese J. Differential effects of COMT on gait and executive control in aging. Neurobiol Aging. 2010; 31:523-31. [PMID: 18547681]
  55. Serý O, Didden W, Mikes V, Pitelova R, Znojil V, Zvolsky P. The association between high-activity COMT allele and alcoholism. Neuroendocrinol Lett. 2006; 27:231-5. [PMID: 16648777]
  56. Dong S, Yang J, Yu W, Kota P, Xia X, Xu H. No association of genetic polymorphisms in CYP1B1 with primary open-angle glaucoma: a meta- and gene-based analysis. Mol Vis. 2012; 18:786-96. [PMID: 22509109]
  57. Hanna IH, Dawling S, Roodi N, Guengerich FP, Parl FF. Cytochrome P450 1B1 (CYP1B1) pharmacogenetics: association of polymorphisms with functional differences in estrogen hydroxylation activity. Cancer Res. 2000; 60:3440-4. [PMID: 10910054]