Inverse Association of Female Hormone Replacement Therapy with Age-Related Macular Degeneration and Interactions with ARMS2 Polymorphisms

Abstract

Purpose.: To investigate whether female reproductive history and hormone replacement therapy (HRT) or birth control pills (BCPs) influence risk for age-related macular degeneration (AMD) and whether genetic factors interact with HRT to modulate AMD risk.

Methods.: Related and unrelated female participants (n = 799) were examined and data were analyzed with generalized estimating equations with adjustment for age and smoking. Individuals with AMD grades 1 to 2 were considered to be unaffected (n = 239) and those with grades 3 to 5 were considered affected (n = 560).

Results.: When comparing all cases with controls, significant inverse associations were observed for HRT (odds ratio [OR] = 0.65, 95% CI 0.48–0.90, P = 0.008) and BCPs (OR = 0.60, 95% CI 0.36–0.10, P = 0.048). When analyses were stratified by AMD severity (early versus geographic atrophy versus neovascular), the inverse association remained significant (HRT OR = 0.45, 95% CI 0.30–0.66, P < 0.0001; BCP OR = 0.55, 95% CI 0.32–0.96, P = 0.036) only when comparing neovascular AMD with the control. All pair-wise HRT-genotype and BCP-genotype interactions were examined, to determine whether HRT or BCP modifies the effect of established genetic risk factors. The strongest interactions were observed for HRT x ARMS2 coding SNP (R73H) rs10490923 (P = 0.007) and HRT x ARMS2 intronic SNP rs17623531 (P = 0.019).

Conclusions.: These findings provide the first evidence suggesting that ARMS2 interacts with HRT to modulate AMD risk and are consistent with previous reports demonstrating a protective relationship between exogenous estrogen use and neovascular AMD. These results highlight the genetic and environmental complexity of the etiologic architecture of AMD; however, further replication is necessary to validate them.

The leading cause of irreversible vision loss among the elderly in the United States today is age-related macular degeneration (AMD). ^1,2 The rate of vision loss among those with AMD is predicted to double by 2020 due to the fast-growing population of those over 65 years of age in the United States. ³ AMD is a complex multifactorial disease consisting of both genetic and environmental risk factors that affect the retina, particularly the macula, and often causes various degrees of central visual loss. ^4–7 In the absence of treatment, severe and irreversible vision loss can occur, highlighting the importance of identifying risk factors and those at risk for AMD. The strongest known risk factors for AMD include age, cigarette smoking, and variants in the complement factor H (CFH), complement factor B (CFB), and age-related macular degeneration 2 (ARMS2) genes. ^8–13 These risk factors, however, explain only a portion of AMD risk, and other risk factors have yet to be identified. 

Results in studies suggest that exposure to estrogen resulting from both hormone replacement therapies (HRTs) and reproductive factors is associated with a lower prevalence of AMD ^14,15 and other retinal diseases. ^16–18 In these studies, an inverse association was observed with current and former use of HRT and/or birth control pill (BCP) use among Caucasian and Latino women with neovascular AMD and pooled analyses including neovascular and non-neovascular AMD. The approximate effect size (OR) was ∼0.50 across studies. ^14,15,19,20 These studies have not consistently replicated the findings despite comparable sample sizes, possibly because analyses were performed with slightly different definitions of HRT and BCP use (e.g., binary ever/never versus years of HRT or BCP use). ^21–23 In many of these studies, the sample size was limited for the variables examined. For example in the Eye Disease Case Control study ¹⁴ fewer than 1% of the women took BCPs. Despite this, an effect size (OR = ∼0.50) comparable to that in studies that have reported statistically significant associations with BCP was observed. The consistency of effects (regardless of statistical significance) suggests that limited power may have been the reason for discrepant results when examining BCP use in AMD. Other studies have observed the associations in univariate but not in multivariate analyses in well age-matched case–control studies, possibly due to loss of power by overfitting the multivariate model. ^24,25 Another explanation of these discrepancies is that these studies did not look at genetic risk factors that may interact with HRT or BCP. Misspecifying the model-influencing risk of AMD by neglecting an interaction term may explain why these studies had inconsistent results. ^26,27 Some studies have made adjustments for potential confounders, including genetic risk factors such as Y402H in the CFH gene; however, none of these studies have looked at gene x environment interactions between HRT or BCP and genetic risk factors. 

The protective direction of the associations observed with HRT and AMD is relatively consistent across studies and, to some extent, with exogenous estrogen use (e.g., oral contraceptives). ^21–23 However, research has shown that oral contraceptive has been associated with an increase in risk for some retinal conditions. ²⁸ Further studies are warranted to determine what role female reproductive history and HRT have on AMD risk. 

Because of the inconsistencies observed in association studies of AMD with HRT and BCPs, we examined female reproductive history and HRT for association with AMD and tested whether these factors might modify the effects of established genetic risk factors: complement factor H (CFH), complement factor H-related proteins (CFH-R1, R3), complement factor B (CFB), age-related macular degeneration 2 (ARMS2), complement component 3 (C3), HtrA serine peptidase 1 (HTRA1), and apolipoprotein E (APOE). 

Materials and Methods

Study Population

Caucasian female participants were recruited into the study by their physicians (EAP, AA, JLK, SGS), through advertisements in waiting rooms and newsletters, presentations to senior groups, and AMD project Web pages. Informed consent was obtained for all patients, and information was collected in compliance with both the internal review boards of all participating institutions and the Health Insurance Portability and Accountability Act of 1996. The study protocol adhered to the tenets of the Declaration of Helsinki. 

All AMD patients underwent an eye examination that included color fundus photographs for AMD evaluation. Study ophthalmologists examined the color fundus photographs of all participants and assigned AMD grades, as previously described. ^29–31 The overall grade assigned each individual was based on the more severely affected eye. Subjects were assigned a grade of 1 through 5, according to a modified classification from AREDS (Age-Related Eye Disease Study), according to both the Wisconsin grading system and the International Classification System. ^32–34 Individuals with grades 1 and 2 were classified as “unaffected” controls; those with disease grade 3 or higher were considered “affected” cases. A detailed description of the clinical definition for each AMD grade is provided in Table 1. 

Table 1.

 

Description of 560 AMD Cases and 239 Controls

Clinical Description	AMD Grade	Cases (n)	Controls (n)
No drusen or small (<63 μm), nonextensive drusen without RPE abnormalities	1	0	169
Extensive small drusen or nonextensive intermediate drusen (63 μm, 125 μm) and/or RPE hyperpigmentation or hypopigmentation	2	0	70
Extensive intermediate drusen or any large soft drusen (125 μm), including drusenoid RPE detachment	3	172	0
Geographic atrophy (area of RPE atrophy with sharp margins, usually visible choroidal vessels, at least 175 μm in diameter)	4	75	0
Extensive AMD, including nondrusenoid RPE detachment, choroidal neovascularization, subretinal hemorrhage, or fibrosis or photocoagulation scarring consistent with treatment of AMD	5	313	0
Total		560	239
Age at examination, mean y (SD)		75.5 (8.3)	66.8 (8.3)
Smoking, % yes, ever		50	41

Evaluation of Female Reproductive History

Female participants who filled out the self-administered health and activities questionnaire and answered the questions regarding female reproductive history, HRT, and BCP use were included in the analysis. The questionnaire was formatted in large print to accommodate those individuals with low vision. In situations in which a participant was not able to complete the questionnaire without assistance, a project coordinator helped the patient with the form. Female reproductive history was assessed by a series of five questions: (1) “How old were you when you started having menstrual periods?” (2) “Have you reached menopause?” 2a) “If yes: At what age did your menstrual periods stop?” (2b) “Why did they stop?” (3) “Have you ever used hormone pills, creams, or patches to help stop premenstrual symptoms (PMS) or hot flashes, or to prevent bone loss caused by menopause (the change of life)?” (3a) “If yes, what year did you start?” (3b) “Are you currently taking hormones?” (3c) “If No: What year did you stop?” (4) “You may have started and stopped taking hormones several times. How long overall have you taken hormones?” (5) “Many women take birth control pills, shots, or implants to prevent pregnancy or to regulate their menstrual periods. Have you ever taken birth control pills, shots, or implants for birth control or for any other reason?” (5a) “If yes, what year did you start?” (5b) “Are you currently taking birth control?” (5c) “If no, what year did you stop?” (6) “You may have started and stopped taking birth control several times. How long overall have you taken birth control?” 

From these questions, the following variables were constructed for analysis: age at menstrual period start (continuous), has menopause been reached (yes, coded 1/no, coded 0), menstrual period stop age (continuous), why did menstrual period stop (categories included, natural menopause, hysterectomy without oophorectomy, partial hysterectomy, radiation chemotherapy, total hysterectomy), HRT (ever, coded 1/never, coded 0), HRT currently (yes, coded 1/no, coded 0), total years of HRT (continuous), BCP use (ever, coded 1/never, coded 0), and total years of BCP use (continuous). For the variable, why did menstrual period stop, each category was tested in a separate analysis treating each as a risk variable (coded 1) comparing it to the baseline (natural menopause). 

Description of the Available Clinical Dataset

As described in Table 1, Caucasian female patients from 492 multiplex and singleton families (560 cases and 239 controls) of non-Hispanic descent who returned the self-administered health and activities questionnaire were included in this analysis. Among these individuals, 392 were unrelated (363 cases and 29 controls), and the remaining individuals made up 100 families, with 66 families containing more than 1 case individual. Of the 239 controls, 169 were classified grade 1 and 70 grade 2. Of the 560 cases, 172 were classified grade 3, 75 grade 4, and 313 grade 5. The mean age at examination for cases (75.5 ± 8.3 years) was higher than the mean age for controls (66.8 ± 8.3 years). Among these individuals, 50% of cases and 41% of controls were smokers (yes, ever). For a more detailed clinical description of the individuals examined in this cohort please refer to Shuler et al. ³⁵  

Genotyping

Genomic DNA was extracted from whole blood (Puregene; Gentra Systems, Minneapolis, MN). DNA was available on 710 (89%) of the 799 Caucasian women with clinical information. Single nucleotide polymorphisms (SNPs) were examined in CFH (11 SNPs), CFH-R1 (2 SNPs), CFH-R3 (1 SNP), CFB (6 SNPs), ARMS2 (9 SNPs), C3 (2 SNPs), HTRA1 (1 SNP), and APOE (2 SNPs defining the E2, E3, and E4 alleles; Supplementary Table S1,). The details regarding SNP selection, genotyping, and quality control measures are described elsewhere. ^{11,30,36–40} A genotyping platform (TaqMan; Applied Biosystems, Inc., [ABI], Foster City, CA) was used to genotype the SNPs (model 7900, ABI using Assays on Demand or Assays by Design when necessary, for the majority of SNPs examined. Assays by Design were used when no predesigned assay was available). Probe and primer sequences for the designed assays are available on request. DNA samples from The Foundation Jean Dausset-Centre d'Etude du Polymorphisme Humain (CEPH) families were duplicated between and across plates for use as a quality control, and laboratory personnel were blinded to the affected status of the individuals being genotyped. 

Statistical Analysis

Analysis of Female Reproductive History and HRT.

To assess the influence of female reproductive history, HRT, and BCP use on risk of AMD, we implemented population averaged generalized estimating equations (GEEs; SAS, Proc GENMOD, ver. 9.1; SAS Institute, Cary, NC). GEE was used to account for the correlations among related individuals, as the data consisted of a combination of related and unrelated female cases and controls. GEE is a valid test of gene x gene and gene x environment interactions in family data. ⁴¹ Population-averaged GEEs with an independence correlation matrix to adjust for correlations among related individuals using a robust variance estimator for each term were used for analyses. Each exposure of interest was included in multivariable analyses adjusting for known AMD risk factors, age at examination (continuous), and smoking status (ever smoked at least 100 cigarettes in lifetime, coded 1, vs. never smoked 100 cigarettes in lifetime, coded 0). Odds ratios (ORs) and 95% confidence intervals (CIs) were estimated from GEE analyses. 

These analyses were performed comparing all AMD cases (moderate AMD, ³ geographic atrophy, ⁴ and neovascular AMD ⁵ ) with controls. However, to investigate whether neovascular AMD significantly differed from early AMD and geographic atrophy AMD, for statistically significant (P < 0.05) associations, we followed up with analyses stratified by AMD grade. 

Genetic Analysis.

Tests for deviations from Hardy-Weinberg equilibrium (HWE), calculations of minor allele frequency, and pair-wise linkage disequilibrium (LD) between SNPs were calculated by using Haploview statistical software. ⁴² For these calculations, one case and one control were randomly selected from each family. Haplotype blocks were defined according the Gabriel et al. ⁴³ algorithm. 

Single-marker tests of association were performed with GEEs, with additive genotypic tests of association used that modeled the minor allele as the risk allele (0, homozygous major allele, vs. 1, heterozygous, vs. 2, homozygous minor allele), with adjustment for age at examination and smoking status. Alleles E2, E3, and E4 at APOE were determined from genotypes at two SNPs and coded into two indicator variables for the effects of E2 and E4. Two analyses were performed for APOE. For the first analysis, genotypes were coded E3/E3 = 0, E2/E3 = 1, and E2/E2 = 2 and for the second analysis genotypes were coded E3/E3 = 0, E4/E3 = 1, and E4/E4 = 2. 

Univariate analyses were followed up, and the joint effects of statistically significant variables (from female reproductive history analyses) and established AMD risk genes, including all genotyped SNPs within those genes, were examined for association with AMD susceptibility. When statistically significant interactions were observed, stratified analyses of SNP genotypes and variable subgroups were performed. All SNPs examined are listed in Supplementary Table S1. We performed interaction analyses using additive genotypic models with the minor allele as the risk allele, with the exception of CFH SNP Y402H (rs1061170), ARMS2 SNP rs1040924, and CFB SNP rs641153, which had all been associated with AMD in a dominant model using the major allele homozygote as the referent exposure. All tests were considered significant at a type I error rate of α = 0.05 for the interaction term. These analyses were performed with pooled affected AMD grades and stratified by grade if a significant association was observed in pooled analyses (early AMD, ³ geographic atrophy, ⁴ and neovascular AMD ⁵ ). 

Results

Female Reproductive History and Hormone Use Analyses

GEE multivariable analyses of all cases and controls (Table 2) found significant inverse associations of AMD and both HRT (P = 0.008) and BCPs (P = 0.048), adjusting for the effects of age at enrollment and smoking. However, when analyses were stratified by “affected” AMD grades, the inverse associations increased in significance for both HRT and BCP when comparing neovascular AMD (grade 5) with controls (grades 1 and 2; HRT OR = 0.45, 95% CI 0.30–0.66, P < 0.0001; BCP OR = 0.55, 95% CI 0.32–0.96, P = 0.036). BCPs consistently had inverse associations across AMD grades, however, for HRT, the inverse association was only observed for grade 5 and all grades pooled. 

Table 2.

 

GEE Multivariate Analyses of Female Reproductive History and HRT Use Adjusted for Age at Examination and Smoking Status

Variable	Cases (n)	Controls (n)	OR	95% CI	GEE P
Menstrual period start, y	478	218	0.99	0.93–1.06	0.846
Menopause, yes/no	527/3	232/3	2.43	0.47–12.53	0.290
Menstrual period stop, y	458	210	0.99	0.96–1.01	0.194
Why did menstrual period stop
Normal menopause	301	133	1.00	—	—
Hysterectomy without oophorectomy	78	28	1.66	0.98–2.83	0.061
Partial hysterectomy	26	9	1.29	0.55–3.02	0.564
Radiation chemotherapy	2	3	0.38	0.11–1.27	0.115
Total hysterectomy	105	45	1.02	0.71–1.47	0.900
HRT (ever/never)
All grades	288/229	173/62	0.65	0.48–0.90	0.008
Grade 5	135/154	173/62	0.45	0.30–0.66	<0.0001
Grade 4	37/27	173/62	1.18	0.60–2.29	0.632
Grade 3	116/48	173/62	1.13	0.74–1.73	0.571
HRT currently, yes/no	105/179	64/107	1.44	0.95–2.18	0.086
HRT total, y	199	137	1.00	0.98–1.02	0.921
BCP use, ever/never
All Grades	125/381	130/102	0.60	0.36–1.00	0.048
Grade 5	56/225	130/102	0.55	0.32–0.96	0.036
Grade 4	12/51	130/102	0.56	0.23–1.35	0.195
Grade 3	57/104	130/102	0.61	0.36–1.05	0.077
BCP use total, y	78	87	1.01	0.97–1.05	0.677

 Statistically significant (P < 0.05) results are in bold.

Genetic Analysis

We followed-up the significant associations we observed (HRT and BCP) with interaction analyses examining the joint effects of both these variables and established AMD genetic risk factors. A complete list of the SNPs and genes examined is in Supplementary Table S1. Demographic information, single locus tests of association, and characterization of LD structure on all our AMD samples have been published ^11,30,37 ; however, in this report, we present the analyses of a subset of those samples (female participants who had been genotyped and completed questions regarding reproductive history, HRT, and BCP use). 

Although we did not observe a statistically significant interaction (P < 0.05) between HRT or BCPs and the specific SNPs previously associated with AMD, we did observe some significant associations at several other markers in those genes. ARMS2 showed the strongest evidence of association, with a statistically significant interaction with HRT at two SNPs (coding SNP R73H rs10490923 x HRT, P = 0.007; intronic SNP rs17623531 x HRT, P = 0.019; Table 3B). These SNPs did not have statistically significant single-marker effects (Table 3A). We did observe, however, five statistically significant associations (rs10490924, A124205188G, rs3750848, rs3750846, and rs10664316; P < 0.01 for all) at other SNPs in ARMS2. None of these SNPs deviated from HWE. In this subset of women, LD analyses showed that rs10490923 and rs17623531 were in perfect LD, with an r ² = 1.0 (Fig. 1), suggesting that these two results probably represent a single interaction with HRT. Although the two SNPs are in complete LD, they do not have identical test results due to a small number of randomly missing genotypes. When analyses were stratified by AMD grade, the interaction with rs17623531 was not significant, but remained for rs10490923 AMD grades 3 (P = 0.046) and 5 (P = 0.044). 

Table 3.

 

GEE Single Locus Association and Genotype × HRT Interaction Analysis Results for ARMS2

A. Single-Marker Significance
Gene Name	SNP	Allele Frequency (Allele)	GEE P *
Age-related macular degeneration 2 (ARMS2) (Chromosome 10)	rs10490923	0.12(A)/0.88(G)	0.116
	rs2736911	0.11(T)/0.89(C)	0.641
	rs10490924	0.32(T)/0.68(G)	<0.0001
	rs17623531	0.12(A)/.088(G)	0.079
	A124205188G	0.31(G)/0.69(A)	<0.0001
	rs3750848	0.32(G)/0.68(T)	<0.0001
	rs3750846	0.30(G)/0.70(A)	<0.0001
	rs2736912	0.10(A)/0.90(G)	0.943
	rs10664316	0.34(C)/0.66(A)	0.001

B. SNP Interaction
Gene	SNP Function	AMD Grade	Interaction	OR	95% CI	GEE p
Age-related macular degeneration 2 (ARMS2) (Chromosome 10)	Coding (R73H)	All grades	rs10490923 × HRT	0.23	0.08–0.68	0.007
	Coding (R73H)	Grade 3	rs10490923 × HRT	0.26	0.07–0.98	0.046
	Coding (R73H)	Grade 4	rs10490923 × HRT	0.37	0.06–2.52	0.313
	Coding (R73H)	Grade 5	rs10490923 × HRT	0.25	0.06–0.96	0.044
	Intron	All grades	rs17623531 × HRT	0.28	0.09–0.81	0.019
	Intron	Grade 3	rs17623531 × HRT	0.27	0.07–1.01	0.052
	Intron	Grade 4	rs17623531 × HRT	0.44	0.07–2.89	0.391
	Intron	Grade 5	rs17623531 × HRT	0.3	0.08–1.18	0.084

C. SNP Significance in Stratified Analysis
Subgroup	Variable	OR	95% CI	GEE P
No HRT	rs10490923	1.47	0.59–3.65	0.407
	rs17623531	1.21	0.48–3.09	0.686
HRT	rs10490923	0.39	0.22–0.67	0.001
	rs17623531	0.38	0.22–0.66	0.001
rs10490923
GG	HRT	0.91	0.54–1.53	0.710
AG&AA	HRT	0.50	0.30–0.84	0.009
rs17623531
GG	HRT	0.82	0.48–1.40	0.469
AG&GG	HRT	0.54	0.34–0.87	0.011

*  GEE analyses were performed with-adjustment for age-at examiantion and smoking status. Statistically significant results are in bold (P < 0.05).

Figure 1.

 

ARMS2 r ² LD plot for unaffected control subjects. LD plots were generated in Haploview and are presented for one unaffected subject selected from each family. Within each triangle is the pair-wise correlation coefficient (r ²). Standard color coding was used for the Haploview LD plots, for D′ LD plots: white (r ² = 0), shaded gray, (0 < r ² < 1), and black (r ² = 1). Black squares without numbers indicate complete LD (r ² = 1.00).

When stratified analyses were performed (stratifying by both genotypes and HRT) we observed that significance at these SNPs was present only in individuals who indicated use of HRT (rs10490923 OR = 0.39, 95% CI 0.22–0.67, P = 0.0007; rs17623531 OR = 0.38, CI 0.22–0.66, P = 0.0006; Table 3C). We also observed that HRT was most significant with the rs10490923 AG&GG (OR = 0.50, 95% CI 0.30–0.84, P = 0.009) and rs17623531 AG&GG (OR = 0.54, 95% CI 0.34–0.87, P = 0.011) genotypes, for which risk was based on having at least one copy of the minor allele (Table 3C). All gene x environment interaction analyses were performed with adjustment for age at examination and smoking. A significant test statistic for the interaction term would indicate that both were synergistic, inverse associations with AMD. 

Discussion

In the present study, we observed significant inverse associations with HRT and BCP use and AMD risk and significant interactions between ARMS2 and HRT. Our strongest interaction was observed with two variants in ARMS2, including one intronic and one coding variant, and HRT adjusted for age at enrollment and smoking. This inverse association was strongest in pooled affected AMD grades 3 to 5. These interactions did not include the specific variants previously associated with AMD risk, although these associations were also in close proximity to the coding variant previously associated with AMD, ARMS2 rs10490924. These results are consistent with some studies that have observed protective effects of HRT and oral contraceptive use in Caucasian and Latino AMD patients ^14,15,20 and in contrast to others in which these effects have not been observed. ^21–23 The inconsistent replication of the associations with BCPs and HRT are not attributable to study design differences, as both studies that have replicated and those that have not have consisted primarily of a cohort study design. Our strongest associations with HRT and BCPs were observed in those individuals with grade 5 or neovascular AMD. This effect, specific to those individuals with the neovascular AMD, has also been observed in previous studies (OR = ∼0.50). ¹⁵  

HRT has been associated with reduced risk of cardiovascular disease in postmenopausal women, with evidence coming from both epidemiologic and clinical studies. ^44,45 However, conflicting results of randomized clinical trials performed by the Women's Health Initiative (WHI) have shown no benefit of oral BCP or HRT on cardiovascular risk ⁴⁶ in some studies and protective effects in others. ⁴⁷ The protective effect that has been observed has been attributed to the lowering of low-density lipoprotein (LDL) cholesterol levels and the elevation of high-density lipoprotein (HDL) cholesterol levels, as well as favorable changes in fibrinogen and PAI-1, enhanced blood flow, and antioxidant properties. ^48–50 Estrogen has also been associated with lowered risk of mortality among women. ¹⁵  

Recently, the ARMS2 A69S polymorphism (rs10490924), which has been associated with AMD, ^12,13 has also been associated with coronary artery disease. ⁵¹ This association with cardiovascular disease is not surprising, given that AMD is a risk factor for cardiovascular disease and has been positively associated with increased mortality among women with AMD compared with women without AMD. ⁵² ARMS2 is a gene of unknown function on chromosome 10q25. It is located in close proximity to PLEKHA, a phosphoinositol-binding protein, and PRSS11, a protein containing both insulin-like growth factor–binding and serine protease domains. ARMS2 A69S is located between 200 to 500 bp from the ARMS2 variants we observed to interact with HRT, although they are not in high r ² with this variant. Presently, the biology underlying the interaction between ARMS2 and HRT is unknown; however, these results suggest that the impact of menopause on blood lipids and lipoproteins puts women at higher risk of several cardiovascular diseases, and women who have specific ARMS2 variants and who take HRT are at decreased risk of these diseases. ⁵³  

A biological pathway to AMD that features HRT may include estrogen receptor α (ESR1) and estrogen receptor β (ESR2). Both ESR1 and ESR2 proteins have been observed in the human retina, suggesting that estrogen plays a role in the pathogenesis of AMD. ^54,55 Estrogen controls expression of chitinase 3-like-1 protein (YKL-40) found in choroidal neovascular membranes, if levels of estrogen are reduced then this can lead to the upregulation of YKL-40 and neovascular AMD. ⁵⁶ Estrogen provides an antioxidant effect by inhibiting lipid peroxidation, which provides protection against oxidative damage in the retina caused by the aging process. ^57,58 A haplotype (PvuII_XbaI haplotype 1) in ESR1 has recently been identified as being associated with lowered serum estradiol levels. ⁵⁹ Studies by Boekhoorn et al., ⁶⁰ using the Rotterdam cohort study population, have shown that PvuII_XbaI haplotype 1 in ESR1 is associated with risk for late AMD, with a stronger effect observed for neovascular AMD (hazard ratio = 4.29, 95% CI 1.47–12.49, with adjustment for age, sex, smoking, and CFH). ⁶⁰ They did not, however, test for a direct interaction between this ESR1 haplotype and HRT or oral contraceptive use. 

Although these results are interesting, there are some limitations to our study. Although our associations with HRT and AMD and the interaction between ARMS2 and HRT remain significant after an FDR correction (q* = 0.20) ⁶¹ for the number of variables and genes examined, further studies are necessary to replicate and extend these findings in independent cohorts. Another factor that may affect our findings is the potential bias caused by our case group's being somewhat younger than our control group. To address this potential bias, we included a term in the multivariate model to adjust for confounding by age. This adjustment, however, may not take into account other unobserved factors such as a different medical awareness of the use of HRT and/or cycles of popularity of BCP by age and case/control status. As a result, we also performed analyses in which our data were divided into quartiles (<60, 60–69, 70–79, 80+; results not shown) for the multivariate analyses of BCP and HRT with adjustment for smoking status. In these analyses, we observed that the protective effect was consistent across age ranges for both HRT and BCP, although the significance of the effect was dependent on the number of samples present for that age quartile. Examination of these factors would be beyond the scope of the present study but they are worth examining in follow-up studies of HRT and BCP, to demonstrate that the protective direction of the effects of HRT and BCP was independent of age quartile. In the past, our group has observed an association between ARMS2 A69S and smoking history. ³⁷ We directly tested all the ARMS2 SNPs for interaction with smoking and did not observe a statistically significant interaction between the ARMS2 SNPs that interacted with HRT and smoking. Therefore, we did not adjust for a smoking interaction in our models. One factor that may explain why we observed the association in variants in the 5′ region of ARMS2 but not at the previously observed ARMS2 variant rs10490924 may be the SNPs tagging variants upstream in the HTRA1 gene located approximately 6.1 kb downstream of ARMS2. HTRA1 has been associated with AMD in a white and Hong Kong Chinese cohort. ^62,63 It may be that we are observing associations with ARMS2 due to a functional variant in HTRA1. We did not observe an interaction with the previously associated HTRA1 SNP rs11200638; however, further studies examining the relationship between these ARMS2 variants and variants in HTRA1 should be conducted. 

Our data indicate that both HRT and oral contraceptives have significant protective effects in women with AMD, particularly those with the neovascular form. We speculate that this effect may be due to the role that estrogen plays in the eye. These findings also lead us to conclude that HRT in women does not seem to have a significant effect on the occurrence of moderate AMD. If these results are validated, they may indicate that treatment with hormones reduces the risk of soft drusen and could be beneficial for prevention of later-stage AMD. Haan et al., ⁶⁴ although with a small sample size, supported this conclusion in a study suggesting that HRT in individuals with neovascular AMD reduces their risk of soft drusen. These findings, although inconsistently replicated in the past, support a role for estrogen levels in AMD susceptibility and suggest that interactions with ARMS2 genetic variants may explain a portion of AMD risk. Further studies replicating these effects in an independent cohort are necessary to validate these findings. Expression studies evaluating the effect of variation in the 5′ region of ARMS2 should be conducted to elucidate the mechanism underlying this gene's effect on AMD risk. 

Supplementary Materials

Supplementary Table S1 - (PDF) 

Footnotes

 Supported by National Eye Institute Grant EY12118 (MP-V, JLH); National Center for Research Resources Grant M01 RR-00095 to Vanderbilt University; and National Eye Institute Center Grant P30-EY014801; and by an unrestricted grant to the University of Miami from Research to Prevent Blindness, New York, NY.

Footnotes

 Disclosure: D.R. Velez Edwards, None; P. Gallins, None; M. Polk, None; J. Ayala-Haedo, None; S.G. Schwartz, Genentech (F), Pfizer (I), P; J.L. Kovach, None; K. Spencer, None; G. Wang, None; A. Agarwal, P; E.A. Postel, P; J.L. Haines, P; M. Pericak-Vance, P; W.K. Scott, P

The authors thank the patients, their families, and the control individuals who participated in the study.