PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2013; 8(1): e53830.
Published online 2013 January 11. doi:  10.1371/journal.pone.0053830
PMCID: PMC3543264

Insights into the Genetic Architecture of Early Stage Age-Related Macular Degeneration: A Genome-Wide Association Study Meta-Analysis

Xiao-Ping Miao, Editor

Abstract

Genetic factors explain a majority of risk variance for age-related macular degeneration (AMD). While genome-wide association studies (GWAS) for late AMD implicate genes in complement, inflammatory and lipid pathways, the genetic architecture of early AMD has been relatively under studied. We conducted a GWAS meta-analysis of early AMD, including 4,089 individuals with prevalent signs of early AMD (soft drusen and/or retinal pigment epithelial changes) and 20,453 individuals without these signs. For various published late AMD risk loci, we also compared effect sizes between early and late AMD using an additional 484 individuals with prevalent late AMD. GWAS meta-analysis confirmed previously reported association of variants at the complement factor H (CFH) (peak P = 1.5×10−31) and age-related maculopathy susceptibility 2 (ARMS2) (P = 4.3×10−24) loci, and suggested Apolipoprotein E (ApoE) polymorphisms (rs2075650; P = 1.1×10−6) associated with early AMD. Other possible loci that did not reach GWAS significance included variants in the zinc finger protein gene GLI3 (rs2049622; P = 8.9×10−6) and upstream of GLI2 (rs6721654; P = 6.5×10−6), encoding retinal Sonic hedgehog signalling regulators, and in the tyrosinase (TYR) gene (rs621313; P = 3.5×10−6), involved in melanin biosynthesis. For a range of published, late AMD risk loci, estimated effect sizes were significantly lower for early than late AMD. This study confirms the involvement of multiple established AMD risk variants in early AMD, but suggests weaker genetic effects on the risk of early AMD relative to late AMD. Several biological processes were suggested to be potentially specific for early AMD, including pathways regulating RPE cell melanin content and signalling pathways potentially involved in retinal regeneration, generating hypotheses for further investigation.

Introduction

Age-related macular degeneration (AMD) is the most common cause of irreversible blindness in the elderly in many countries [1]. Clinically, the disease is believed to develop via a series of progressive stages [2]. Early AMD is characterised by abnormalities in the retinal pigment epithelium (RPE) and the deposition of small extracellular deposits, called drusen, between Bruch’s membrane and the RPE [3]. As the disease progresses, large drusen (>125 µm in diameter) with or without fuzzy edges (soft indistinct or distinct drusen) may slowly disappear and be replaced by regions of retinal depigmentation (decreased pigment). Reticular drusen and large, soft drusen involving a large area of the macula are strong indicators of increased risk of progression to late, vision-impairing forms of this disease [4]. The two main phenotypes of late AMD are exudative AMD and geographic atrophy (GA). Exudative AMD is typified by sub-retinal neovascularisation with sensory retinal and/or RPE detachment, sub-retinal and/or sub-RPE haemorrhage followed by sub-retinal scarring. GA involves gradual degeneration and disappearance of RPE and photoreceptor cells within the macular area [3].

AMD is known to result from a complex interplay of multiple environmental and genetic factors. Disease risk increases strongly with age, with additional risk factors including smoking, possibly high body mass index, hypertension and cardiovascular disease [5], [6]. However, genetic factors have been shown to be strongly associated with AMD risk, with an increased genetic burden associated particularly with late forms of the disease [7]. For example, heritability estimates from twin studies are 0.45 for early AMD [8], but 0.71 for late AMD [9].

Case-control genome-wide association studies (GWAS) of late AMD have provided evidence for individual, large-effect risk variants in the CFH, ARMS2, C2/CFB and C3 genes (estimated per allele odds ratios ranging from 1.7–4.5), and for additional, smaller-effect variants in CFI, FRK/COL10A1, TNFRSF10A, LIPC, CETP, TIMP3, REST and VEGFA (odds ratios ranging from 1.15–1.4) [10], [11], [12], [13], [14], [15]. Together with biological evidence, these findings highlight the importance of immune, inflammatory and lipid metabolic pathways in late AMD pathogenesis, while also implicating angiogenic, apoptotic, extracellular-matrix remodelling and melanosome trafficking processes [10], [11], [12], [13], [14], [15].

The vast majority of published GWAS have been conducted in samples of either exclusively or predominant late AMD cases (particularly late exudative AMD). This likely reflects the tendency for only symptomatic, late AMD cases to attend clinics, thus forming the majority of ascertained AMD samples available. In contrast, persons with early AMD are usually asymptomatic and less likely to be seen except in population-based studies. Therefore the genetic architecture of early AMD has been relatively under researched and is poorly understood [16].

Given the substantial heritability of AMD and the knowledge of modifiable risk factors that may reduce risk of progression to late, vision-impairing forms of this disease, improved understanding of the genetic architecture of early AMD may also be important. With this aim, we conducted a GWAS meta-analysis of early AMD, including approximately 4,000 well-characterised early AMD cases and 20,000 strictly defined controls without any drusen or with hard drusen only. Using a small set of mutually exclusive, late AMD cases, we also compared genetic effect sizes at validated AMD risk loci between early and late stages of the disease, to determine their relative importance for different disease stages. This study amalgamates a number of large population-based cohorts with GWAS data and AMD grading from retinal photographs.

Materials and Methods

Study Populations

The primary GWAS meta-analysis for early AMD was conducted in five European-ancestry cross-sectional cohorts (Table 1). These were recruited in the USA, Europe and Australia and contributed a total of 3,772 prevalent cases of early AMD and 16,033 contemporaneous controls from the Age, Gene/Environment Susceptibility-Reykjavik Study (AGES) [17], the Atherosclerosis Risk in Communities (ARIC) study [18], the Cardiovascular Health Study (CHS) [19], the Blue Mountains Eye Study (BMES) [20] and three distinct cohorts from the Rotterdam Study (RS) [21]: RS-I, RS-II and RS-III. In addition, two Asian-ancestry cross-sectional cohorts were included in secondary analyses; these included a total 264 prevalent early AMD cases and 3,926 controls from the Singapore Indian Eye Study (SINDI) [22] and the Singapore Malay Eye Study (SiMES) [23]. Subsequent candidate SNP meta-analyses for late AMD were performed in the European-ancestry cohorts including a total 498 prevalent cases of late AMD and 16,033 controls.

Table 1
Population characteristics for the individual cohorts.

Each cohort obtained approval from relevant institutional review boards, and all participants provided written informed consent in accordance with the Declaration of Helsinki. All participating studies approved guidelines for collaboration, and a working group agreed, in advance, on phenotype definition, covariate selection and analytic plans for within-study analyses and meta-analyses of results. Details of each participating study are described below and in Tables S1, S2 in File S2.

Phenotype Definitions

The same AMD phenotype definitions based on photographic grading were used for all the cohorts. Early AMD was defined as the presence of soft drusen (>63 µm) alone, retinal pigment epithelium (RPE) depigmentation alone or a combination of soft drusen with increased retinal pigment and/or depigmentation in the absence of late AMD. Late AMD was defined as the presence of exudative AMD or GA, as described in the International AMD classification [24]. Controls had no soft (distinct or indistinct) drusen or retinal pigment abnormalities (either depigmentation or increased pigment), and no signs of early or late AMD; controls were permitted to have hard drusen. Presence and severity of AMD lesions were assigned following the Wisconsin Age-Related Maculopathy grading system [25], based on masked assessment of fundus photographs. For the AGES, BMES, RS and Singapore-based cohorts, photographs were examined for both eyes using a retinal camera with pharmacological mydriasis, with cases satisfying AMD diagnostic criteria for at least one eye. For the ARIC and CHS cohorts, case and control diagnoses were based on examination of one, randomly selected eye using a retinal camera without mydriasis. Additional details are provided in Table S1 in File S2.

Genotyping

The AGES and CHS samples were genotyped using the Illumina Human370 CNV quad array. The ARIC sample was genotyped using the Affymetrix SNP 6.0 GeneChip. The RS-I sample was genotyped using the Illumina Infinium II HumanHap550 chip. The RS-II, RS-III, SiMES and SINDI samples were genotyped using the Illumina Human610-Quad array. The BMES sample was genotyped using the Illumina Human670-Quad v1 custom array. All cohorts applied similar quality control (QC) procedures to genotype data prior to analysis. Briefly, this involved excluding SNPs with low call rate, pronounced deviation from Hardy-Weinberg equilibrium or low minor allele frequency, and individuals with low call rate, discrepant clinical and genotypic gender, evidence of cryptic relatedness based on IBS sharing (one member of each pair excluded) or outlying continental ancestry based on principal components analysis. The particular procedures and thresholds used for each study are detailed in Table S2 in File S2.

Following genotype quality control, all cohorts were imputed to approximately 2.5 million HapMap Phase II SNPs. European ancestry cohorts used the Caucasian (CEU) reference panel, and Asian ancestry cohorts used the combined Caucasian, Chinese, Japanese and African reference (CEU+CHB+JPT+YRI) panel. This provided data for a large set of overlapping SNPs for GWAS meta-analyses.

Statistical Analysis

Primary meta-analysis of genetic associations with early AMD compared to controls

Individual European-ancestry cohorts conducted GWAS of early AMD. For each SNP, logistic regression was performed using a one-degree of freedom trend test assuming an additive effect of SNP allele dosage; logistic models included age, gender and the first two ancestry principal components as covariates. Regression coefficients and their standard errors were estimated to represent the change in log odds of AMD resulting from each additional copy of the test allele, after adjusting for specified covariates. Analyses were performed using the ProbABEL program [26] for AGES and ARIC cohorts, mach2dat [27] for BMES, R software for CHS and GRIMP software for the RS cohorts [28].

SNPs were included in meta-analyses if they exhibited an imputation quality score (defined as the ratio of observed to expected dosage variance) >0.3, minor allele frequency >0.01 and valid data in at least two cohorts. Strand data was available for each cohort and all results were synchronised to the forward strand. Fixed effects, inverse variance-weighted meta-analyses were performed using METAL software. The pre-specified threshold for genome-wide significant association was set at 5×10−8. SNPs not meeting this threshold but achieving P<1×10−5 were considered as suggestively associated; that is, potentially true signals, whose association may become more significant with increased sample size. Between-cohort heterogeneity of estimated SNP association with disease was assessed using I2 and Cochran’s Q metrics. Quantile-Quantile (Q-Q) plots and Manhattan plots summarising the meta-analysis results were produced using R software. Regional association results were plotted using LocusZoom software [29]. Meta-analysis genomic control inflation factors (λGC) were calculated as previously described [30], as were standardised values for a sample of 1000 cases and 1000 controls (λ1000) [30].

Comparison of effect sizes between early and advanced AMD for established AMD risk variants

In addition to the primary GWAS of early AMD, in European-ancestry cohorts, we performed association analyses for advanced AMD using a set of confirmed AMD-associated SNPs. Asian-ancestry cohorts had insufficient advanced AMD cases to perform these analyses. Confirmed AMD-associated SNPs were required to demonstrate genome-wide significant association (P<5×10−8) with advanced AMD in a published report, and were identified using the NIH National Human Genome Research Institute (NHGRI) Catalog of Published Genome-Wide Association Studies with search terms “Age-related macular degeneration” and “Age-related macular degeneration (wet)”. The identity and association evidence for SNPs returned by the NHGRI search were corroborated with original published reports in each case. In turn, the original reports were checked for appropriate AMD phenotype definition and searched to ensure no SNPs meeting eligibility criteria had been omitted from the NHGRI catalog. Statistical analyses of individual cohorts and meta-analyses of estimated SNP effect sizes were performed as described above.

For the selected set of confirmed AMD-associated SNPs, estimated SNP effect sizes for late and early AMD were statistically compared using summary effect estimates from the European-ancestry meta-analyses. The difference between the regression coefficients was assessed using a two-tailed Z test [31] formulated as

equation image

where An external file that holds a picture, illustration, etc.
Object name is pone.0053830.e002.jpg, An external file that holds a picture, illustration, etc.
Object name is pone.0053830.e003.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0053830.e004.jpg, An external file that holds a picture, illustration, etc.
Object name is pone.0053830.e005.jpg represent the meta-analysis summary effect estimates and their standard errors for late and early AMD, respectively. A total of 21 SNPs were assessed for effect size differences. Each test was performed at a pre-specified significance level of 0.00238, corresponding to a family-wise type I error rate of 0.05 after Bonferroni adjustment. We note that this threshold is somewhat conservative, as multiple, correlated SNPs were included from the CFH and ARMS2 loci.

Secondary, Trans-ethnic Meta-analysis of Genetic Associations with Early AMD

To extend on the results of the European-ancestry analysis, a secondary meta-analysis was conducted including all European-ancestry and Singapore-based, Asian-ancestry cohorts. Analyses in individual cohorts were performed via logistic regression as described for the primary European analysis. Fixed effects, inverse variance-weighted meta-analyses including all European (AGES, ARIC, BMES, CHS and RS) and Singapore-based (SiMES and SINDI) cohorts were performed using METAL. For a set of selected, validated AMD-associated variants (described above), we also tabulated ancestry-specific effect sizes separately for European- and Asian-ancestry cohorts.

Results

Primary Meta-analysis of Genetic Associations with Early AMD

For each of the individual cohorts, summary demographic characteristics are shown by AMD status in Table 1. A Q-Q plot summarising the distribution of all P-values showed excellent agreement with the expected null distribution throughout all but the extreme tail of the distribution (Figures S2 and S5a in File S1). The overall meta-analysis genomic control inflation factor (λGC = 1.021) indicated an absence of confounding by population stratification or other artefacts, as did the scaled, standardised value (λ1000 = 1.003).

The primary meta-analysis detected highly significant association of numerous SNPs within CFH and ARMS2/HTRA1 loci with risk of early AMD (Table 2, Table S4 in File S3, Figure S1 in File S1). The peak results at these loci were detected for the CFH SNP rs1329424 (OR[T] = 1.41, 95% CI: 1.33–1.50, P-value = 1.5×10−31) and the ARMS2/HTRA1 SNP rs3793917 (OR[G] = 1.43, 95% CI: 1.34–1.54, P-value = 4.3×10−24).

Table 2
Results for SNPs showing suggestive evidence of association (P<1×10−5) in the primary (European-ancestry) meta-analysis of early AMD.

Aside from SNPs at the CFH and ARMS/HTRA1 loci, no additional variants reached genome-wide significance. However, a number of confirmed and plausible candidate genes contained SNPs demonstrating possible association with early AMD (5×10−8<P<1×10−5) (Table 2 and Table S5 in File S3). These included SNPs in the ApoE/TOMM40 gene cluster on chromosome 19, the GLI3 gene on chromosome 7 and the tyrosinase gene (TYR) on chromosome 11 (see Figures 1, ,2,2, ,3),3), and SNPs in the vicinity of the complement-related genes CD46 and MBL2.

Figure 1
Regional association results for chromosome 7 SNPs in the GLI-Kruppel family member GLI3 gene.
Figure 2
Regional association results for chromosome 11 SNPs in the tyrosinase precursor (TYR) gene.
Figure 3
Regional association results for chromosome 19 SNPs in the PVRL2/APOE/TOMM40 gene cluster.

Comparison of Effect Sizes between Early and Late AMD for Established AMD Risk Variants

SNPs in genes C2/CFB, C3, CFI, CETP, TIMP3, TNFRSF10A, FRK/COL10A1, REST and LIPC (upstream) that have previously demonstrated genome-wide significance with advanced AMD showed only modest association with early AMD (Table 3). Of these published associated variants, the strongest associations with early AMD were observed for rs429608 in C2/CFB (OR[G] = 1.18, 95% CI: 1.08–1.27, P-value = 9.6×10−5), rs2230199 in C3 (OR[C] = 1.18, 95% CI: 1.08–1.29, P-value = 2.5×10−4) and rs13278062 in TNFRSF10A (OR[T] = 1.11, 95% CI: 1.05–1.18, P-value = 5×10−4). For many published SNPs we detected some evidence of association with both late and early AMD (at α = 0.05), observing similar odds ratios for late AMD as those previously reported (Table 3). Notably, for most of these SNPs the estimated effect size for late AMD was several-fold higher than that for early AMD. The absolute fold-change was ~3-4-fold for risk variants in CFH, ~3-fold for variants in ARMS2/HTRA1, ~3-6-fold for variants in C2/CFB, ~4-fold for a SNP in CETP and ~2-fold for variants in C3 and FRK/COL10A1. For SNPs in CFH, ARMS2/HTRA1, and C2/CFB, the heterogeneity of effect size between early and late AMD was statistically significant after multiple testing adjustment (P<0.05). The only SNP showing some evidence for association but a weaker effect for late AMD was rs10033900 in CFI, which showed very marginal association (P = 0.038) with early AMD but no association with late AMD, and had a slightly higher estimated effect for early than late AMD (1.8-fold greater); however the difference in effect size was not statistical significant (P>0.05). Previously published AMD risk variants in REST, LIPC upstream and TIMP3 showed no evidence of association (at α = 0.05) with either late or early AMD in our study.

Table 3
Comparison of estimated effect sizes for early versus advanced AMD for published SNPs showing genome-wide significant association with AMD.

Due to the limited number of late AMD cases, our direct, internal comparison provided modest power to detect effect size differences between early and late AMD. Hence, we repeated the effect size comparisons for validated AMD-associated SNPs, using published effect estimates from well-powered GWAS of late AMD. Table 4 shows analogous comparisons as shown Table 3, with late AMD effect estimates replaced with published effect estimates. This analysis produced largely consistent results as the internal comparison. However, in addition to SNPs in CFH, ARMS2/HTRA1 and C2/CFB numerous additional SNPs showed significant effect size differences between early and late AMD – including SNPs in REST, TNFRSF10A, C3, TIMP3 and upstream of LIPC.

Table 4
Comparison of effect sizes for early AMD from this study versus published effect estimates for late AMD.

Secondary, Trans-ethnic Meta-analysis of Genetic Associations with Early AMD

Secondary meta-analysis including all European samples plus the two Asian-ancestry samples produced very similar results to the primary analysis. There was no evidence of systematic bias of test statistics (Figure S4 in File S1; λGC = 1.022, λ1000 = 1.003) and SNPs in the CFH and ARMS2/HTRA1 loci reached genome-wide significance (Figure S3 in File S1). SNPs in the vicinity of ApoE achieved slightly greater significance (peak P = 9.7×10−7 at rs6857), while the peak evidence for TYR (P = 9×10−6 at rs10830228) and GLI3 (P = 1×10−5 at rs2049622) was slightly reduced. Several additional intergenic regions were also identified at P<10−5 (Table S6 in File S3).

For the set of selected, validated AMD-associated variants (described above), ancestry-specific effect sizes are shown separately for European and Asian ancestry cohorts in Table 5. Considering the low power of the Singapore-based analysis – with <300 early AMD cases represented across the two constituent samples – the effect sizes in the European and Asian-ancestry cohorts were largely consistent for the candidate SNPs, with confidence intervals for effect sizes overlapping for 8 of the 16 assessed SNPs. The principal exception was for the 6 SNPs located within the CFH gene, which showed strong association in the European analysis, but no association in our small Singapore-based analysis.

Table 5
Tabulation of estimated effect sizes for European-ancestry versus Asian ancestry analyses of early AMD, for the validated AMD-associated SNPs shown in Table 3.

Discussion

We report herein findings from GWAS meta-analysis restricted to individuals with early AMD, prior to the development of late, vision-impairing disease. As expected, SNPs at the CFH and ARMS2/HTRA1 loci, which were discovered for late AMD, comprised strong and significant risk factors for the development of early AMD as well. However, these variants were shown to impart a 3–4 times significantly weaker risk for early than late AMD. Similar effect size ratios for early versus late AMD were observed for additional risk variants in C2, CFB, C3 and CETP, with the difference reaching significance for a SNP in C2/CFB. These findings are consistent with those of studies showing that risk variants in CFH, ARMS2, C2, CFB and C3 influence not only risk of late AMD (based on comparison with normal subjects), but also risk of progression from early to late stages of disease (based on comparison with early AMD cases that did not progress) [32], [33], [34], [35]. A possible explanation for these findings is that the estimated genetic effect size for late AMD comprises cumulative effects of an initial development of early AMD signs and progression from early to late stages of the disease.

It is noteworthy that our large study offered similar sample size to assess genetic association with early AMD to a number of previous GWAS that assessed the associations with late AMD [10], [11], [13], [14], [15], but in contrast to other studies, we detected no significant variants besides CFH and ARMS2 variants at genome-wide significance. This may suggest weaker genetic effects for early AMD, and a possible absence of large effect loci unique to early AMD.

We speculate that the genetic component of early AMD may include risk variants both shared with late AMD, unique to early AMD, and perhaps unique to specific lesions classified here altogether as early AMD. Many confirmed AMD risk variants appear to confer both modest risk of early AMD and late AMD; the presence of these variants may encourage both early disease incidence and steady disease progression throughout time. Additionally, we found suggestive evidence for possible early AMD variants in several novel and biologically plausible genes not previously associated with late AMD. These variants (e.g. within TYR, GLI3 and upstream of GLI2) may influence only the initial development of early pathologic features of AMD such as RPE abnormalities and large soft drusen, with alternative pathways mediating disease progression once early stage changes have accumulated.

The GLI3 and GLI2 transcription factors are critical mediators of Sonic Hedgehog (Shh) protein signalling, with species-conserved roles in a range of developmental and adult processes [36]. Notably, Shh signalling is required for the development of Muller glial stem cells [37], a retinal progenitor cell (RPC) with both proliferative and neurogenic ability [38]. Although retinal regeneration has not yet been demonstrated in the adult human, in fish and amphibians acute retinal damage stimulates RPCs of the ciliary marginal zone (CMZ) to proliferate and fully regenerate the multilayered retina throughout life [39](24). Muller glia isolated from the adult human retinal marginal region can also re-enter the cell cycle and demonstrate proliferative potential in vitro [38]. The presence of a Muller cell-specific regulatory region upstream of the HTRA1 promoter [40], and the upregulation of glial fibrillary acidic protein (GFAP) immunoreactivity in the Müller cells of donor retina with drusen [41] may be considered supporting evidence for the involvement of Muller glia in early stages of AMD. It is also interesting to note that GLI3 and GLI2 are zinc finger proteins [42], as zinc supplements have been shown reduce the risk of progression from early to late AMD [43].

We also detected possible association of multiple variants in the TYR gene encoding tyrosinase, the enzyme catalysing the first step of melanin biosynthesis. Melanin is produced within pigmented cells including skin melanocytes and retinal pigment epithelium (RPE) cells [44] and deficiencies of melanin in RPE cells have been associated with AMD [44], [45]. Variants in TYR, including SNPs correlated with our lead SNP, have previously shown strong association with skin, hair and eye colour [46], [47], tanning ability [48], vitiligo [49] and melanoma risk [50], [51] via GWAS. In addition to protecting from sunlight-induced damage [52], melanin is also an efficient antioxidant, reducing the oxidative stress resulting from lipid peroxidation and reactive oxygen species (ROS) generation [53]. Specifically, elevated melanin content in RPE cells has been shown to reduce accretion of the auto-oxidant lipofuscin [54], a key waste product of photoreceptor outer segment phagocytosis implicated in extracellular drusen formation [3]. This antioxidant function of melanin is consistent with evidence supporting supplemental lutein and zeaxanthin in modifying the course of AMD [55].

The detection of suggestively associated SNPs in the vicinity of the Apolipoprotein E gene (ApoE) is not novel, but is consistent with early implication of ApoE polymorphisms in AMD [56], [57], [58], [59], [60]. Our most strongly associated SNPs have also previously shown unequivocal association with Alzheimer’s disease diagnosis, pathologic features and biomarkers [61], [62], [63], [64], [65], as well as lipid levels and cardiovascular traits [66], [67] in GWAS.

A potential limitation of this study is possible measurement error in early AMD ascertainment. The milder nature of early AMD signs could result in misclassification of some controls, particularly in the two studies (ARIC and CHS) that examined only non-stereoscopic color retinal photographs of one eye per subject, using non-mydriatic cameras. An estimated misclassification rate for cases and controls could be around 10–20% if using data from one eye per subject. Misclassification would be expected to bias genetic effect estimates for early AMD towards the null, and thus increase the apparent effect differences between early and late AMD. We note though, that the majority of participating cohorts in our study examined stereoscopic retinal photographs of both eyes per subject. Inter-center grading reliability of early AMD has been assessed in a different study across three grading centres (the Wisconsin, BMES and RS) based on stereoscopic images, and showed, for example, 70.2% exact agreement for a 5 step AMD severity scale, or 90.4% agreement if allowing 1 step difference between Wisconsin and BMES graders (personal communication with R.Klein). Furthermore, in sensitivity analyses excluding the ARIC and CHS cohorts that used non-stereoscopic photographs of one eye per subject, very little difference in effect estimates was observed for the selected candidate SNPs. In fact, effect estimates tended to be slightly smaller than in the full European meta-analysis (the mean difference in regression coefficients was -0.012, with a standard deviation of only 0.025). If measurement error was a major contributor to the observed effect size differences between early and late AMD, we would have expected the opposite result, i.e. larger effect sizes after reducing measurement error. In addition, a recent Age-relate Eye Disease Study (AREDS) report involving longitudinally assessment of stereoscopic images of AMD patients also documented differential effects of different AMD candidate SNPs on various stages of AMD (from normal to early and then late AMD) [35]. Taken together, this suggests that misclassification error is not an important contributor to the observed effect size differences between early and late AMD detected in this study.

We also acknowledge the potential influence of increased genetic and phenotypic heterogeneity of early AMD, compared to late AMD. If the phenotypic signs of early AMD have a more complex genetic basis, including a larger number of contributing genetic variants, it is possible that a particular AMD risk variant may act in a smaller proportion of subjects than in late AMD. Under this scenario, the effect size in the relevant early AMD subset may be similar to the effect upon late AMD, but the sample-wide estimate will be reduced. Genetic heterogeneity may partially result from phenotypic heterogeneity, if alternative biological processes lead to a clinical presentation that appears as early AMD, but does not progress to late AMD. Although very difficult to quantify or control, such genetic heterogeneity may have contributed to our observed effect estimate differences.

Finally, our analyses of late AMD had limited power, resulting in wide confidence intervals for effect estimates and reducing power to detect effect size differences between early and late AMD. We note however, that analogous comparisons using published estimates from well-powered GWAS of late AMD produced strikingly similar results. As expected, numerous additional SNPs showed significant effect size differences between early and late AMD, but the general pattern of several-fold increased effect sizes for late AMD was clearly consistent with results of the direct internal comparison conducted in our study.

In conclusion, the results of our study confirm the involvement of several established late AMD risk variants in early AMD and provide additional, suggestive evidence for possible risk variants and biological pathways specific to early AMD. These include TYR SNPs previously associated with skin and eye pigmentation, and variants in and upstream of GLI3 and GLI2, respectively, potentially influencing retinal regeneration following injury. Our study also demonstrates that many established late AMD genetic risk variants showed reduced effects on early AMD compared to late AMD. Further research should seek to clarify the underlying biological processes involved in early AMD, potentially uncovering novel preventative therapies to prevent the progression of early to late, vision-threatening stages of AMD.

Web Resources

ProbABEL program: http://mga.bionet.nsc.ru/~yurii/ABEL/

mach2dat software: http://www.sph.umich.edu/csg/yli/mach/index.html

METAL software: http://www.sph.umich.edu/csg/abecasis/Metal/index.html

R software: http://www.r-project.org

(NHGRI) Catalog of Published Genome-Wide Association Studies: http://www.genome.gov/gwastudies/

Supporting Information

File S1

Supplementary Figures S1– S5.

(DOC)

File S2

Supplementary Tables S1– S3.

(DOC)

File S3

Supplementary Tables S4– S6.

(XLS)

Funding Statement

Age, Gene/Environment Susceptibility–Reykjavik Study has been funded by National Institutes of Health (NIH) contract N01-AG-12100, the NIA and NEI Intramural Research Programs at the National Institutes of Health (ZIAAG007380 and ZIAEY000401), Hjartavernd (the Icelandic Heart Association), and the Althingi (the Icelandic Parliament). The Atherosclerosis Risk in Communities Study is carried out as a collaborative study supported by National Heart, Lung, and Blood Institute contracts HHSN268201100005C, HHSN268201100006C, HHSN268201100007C, HHSN268201100008C, HHSN268201100009C, HHSN268201100010C, HHSN268201100011C, HHSN268201100012C, R01HL087641, R01HL59367 and R01HL086694; National Human Genome Research Institute contract U01HG004402; and National Institutes of Health contract HHSN268200625226C. Infrastructure was partly supported by Grant Number UL1RR025005, a component of the National Institutes of Health and NIH Roadmap for Medical Research. The Blue Mountains Eye Study was supported by the Australian National Health & Medical Research Council (NHMRC) project grants (IDs 974159, 991407, 211069 and 457349) and Centre for Clinical Research Excellence (CCRE) in Translational Clinical Research in Eye Diseases, CCRE in TCR-Eye (ID 529923). The Blue Mountains Eye Study GWAS and genotyping costs were supported by Australian NHMRC project grants (IDs 512423, 475604,529912 and 590204), and the Wellcome Trust, United Kingdom, as part of Wellcome Trust Case Control Consortium 2 (grant IDs 085475/B/08/Z and 085475/08/Z). EGH (631096), PNB (APP1028444) and JJW (358702 and 632909) are supported by the NHMRC fellowship scheme. This CHS research was supported by NHLBI contracts N01-HC-85239, N01-HC-85079 through N01-HC-85086; N01-HC-35129, N01 HC-15103, N01 HC-55222, N01-HC-75150, N01-HC-45133, HHSN268201200036C and NHLBI grants HL080295, HL087652, HL105756 with additional contribution from NINDS. Additional support was provided through AG-023629, AG-15928, AG-20098, and AG-027058 from the NIA. See also http://www.chs-nhlbi.org/pi.htm. DNA handling and genotyping was supported in part by National Center of Advancing Translational Technologies CTSI grant UL1TR000124 and National Institute of Diabetes and Digestive and Kidney Diseases grant DK063491 to the Southern California Diabetes Endocrinology Research Center and Cedars-Sinai Board of Governors' Chair in Medical Genetics (JIR). The GWA database of the Rotterdam Study was funded through the Netherlands Organization of Scientific Research NWO (nr. 175.010.2005.011). The Rotterdam Study is supported by the Erasmus Medical Center and Erasmus University, Rotterdam; the Netherlands Organization for Scientific Research (NWO), the Netherlands Organization for Health Research and Development (ZonMw), the Research Institute for Diseases in the Elderly (RIDE), the Ministry of Education, Culture and Science, the Ministry for Health, Welfare and Sports, the European Commission (DG XII), and the Municipality of Rotterdam. The ophthalmologic part of the Rotterdam Study was supported by Lijf en Leven, Krimpen a/d Lek; MD Fonds, Utrecht. Oogfonds Nederland, Utrecht; Stichting Nederlands Oogheelkundig Onderzoek, Nijmegen/Rotterdam; Swart van Essen, Rotterdam; Netherlands Organisation for Scientific Research (NWO); Bevordering van Volkskracht, Rotterdam; Blindenhulp, The Hague; Rotterdamse Vereniging Blindenbelangen, Rotterdam; OOG, The Hague; Algemene Nederlandse Vereniging ter Voorkoming van Blindheid, Doorn; Blinden-Penning, Amsterdam; Blindenhulp, ’s Gravenzande; Henkes Stichting, Rotterdam; Topcon Europe BV, Capelle aan de IJssel; Medical Workshop BV, Groningen; all in the Netherlands; Heidelberg Engineering, Dossenheim, Germany. The Singapore Indian Eye Study (SINDI) was funded by grants from Biomedical Research Council of Singapore (BMRC 08/1/35/19/550), National Medical Research Council of Singapore (NMRC/STaR/0003/2008 and CG/SERI/2010), and American Health Assistance Foundation, USA. The Singapore BioBank and the Genome Institute of Singapore, Agency for Science, Technology and Research, Singapore provided services for tissue archival and genotyping, respectively. The Singapore Malay Eye Study (SiMES) is funded by National Medical Research Council (grants 0796/2003, IRG07nov013, IRG09nov014, STaR/0003/2008 and CG/SERI/2010) and Biomedical Research Council of Singapore (grants 09/1/35/19/616), Singapore and the American Health Assistance Foundation, USA. The Singapore Tissue Network and the Genome Institute of Singapore, Agency for Science, Technology and Research, Singapore provided services for tissue archival and genotyping, respectively. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Lim LS, Mitchell P, Seddon JM, Holz FG, Wong TY (2012) Age-related macular degeneration. Lancet 379: 1728–1738. [PubMed]
2. Sallo FB, Peto T, Leung I, Xing W, Bunce C, et al. (2009) The International Classification system and the progression of age-related macular degeneration. Curr Eye Res 34: 238–240. [PubMed]
3. Kinnunen K, Petrovski G, Moe MC, Berta A, Kaarniranta K (2011) Molecular mechanisms of retinal pigment epithelium damage and development of age-related macular degeneration. Acta Ophthalmol. [PubMed]
4. Abdelsalam A, Del Priore L, Zarbin MA (1999) Drusen in age-related macular degeneration: pathogenesis, natural course, and laser photocoagulation-induced regression. Surv Ophthalmol 44: 1–29. [PubMed]
5. Chakravarthy U, Wong TY, Fletcher A, Piault E, Evans C, et al. (2010) Clinical risk factors for age-related macular degeneration: a systematic review and meta-analysis. BMC Ophthalmol 10: 31. [PMC free article] [PubMed]
6. Thornton J, Edwards R, Mitchell P, Harrison RA, Buchan I, et al. (2005) Smoking and age-related macular degeneration: a review of association. Eye (Lond) 19: 935–944. [PubMed]
7. Heiba IM, Elston RC, Klein BE, Klein R (1994) Sibling correlations and segregation analysis of age-related maculopathy: the Beaver Dam Eye Study. Genet Epidemiol 11: 51–67. [PubMed]
8. Hammond CJ, Webster AR, Snieder H, Bird AC, Gilbert CE, et al. (2002) Genetic influence on early age-related maculopathy: a twin study. Ophthalmology 109: 730–736. [PubMed]
9. Seddon JM, Cote J, Page WF, Aggen SH, Neale MC (2005) The US twin study of age-related macular degeneration: relative roles of genetic and environmental influences. Arch Ophthalmol 123: 321–327. [PubMed]
10. Arakawa S, Takahashi A, Ashikawa K, Hosono N, Aoi T, et al. (2011) Genome-wide association study identifies two susceptibility loci for exudative age-related macular degeneration in the Japanese population. Nat Genet 43: 1001–1004. [PubMed]
11. Chen W, Stambolian D, Edwards AO, Branham KE, Othman M, et al. (2010) Genetic variants near TIMP3 and high-density lipoprotein-associated loci influence susceptibility to age-related macular degeneration. Proc Natl Acad Sci U S A 107: 7401–7406. [PubMed]
12. Klein RJ, Zeiss C, Chew EY, Tsai JY, Sackler RS, et al. (2005) Complement factor H polymorphism in age-related macular degeneration. Science 308: 385–389. [PMC free article] [PubMed]
13. Kopplin LJ, Igo RP Jr, Wang Y, Sivakumaran TA, Hagstrom SA, et al. (2010) Genome-wide association identifies SKIV2L and MYRIP as protective factors for age-related macular degeneration. Genes Immun 11: 609–621. [PMC free article] [PubMed]
14. Neale BM, Fagerness J, Reynolds R, Sobrin L, Parker M, et al. (2010) Genome-wide association study of advanced age-related macular degeneration identifies a role of the hepatic lipase gene (LIPC). Proc Natl Acad Sci U S A 107: 7395–7400. [PubMed]
15. Yu Y, Bhangale TR, Fagerness J, Ripke S, Thorleifsson G, et al. (2011) Common variants near FRK/COL10A1 and VEGFA are associated with advanced age-related macular degeneration. Hum Mol Genet 20: 3699–3709. [PMC free article] [PubMed]
16. Francis PJ, Klein ML (2011) Update on the role of genetics in the onset of age-related macular degeneration. Clin Ophthalmol 5: 1127–1133. [PMC free article] [PubMed]
17. Harris TB, Launer LJ, Eiriksdottir G, Kjartansson O, Jonsson PV, et al. (2007) Age, Gene/Environment Susceptibility-Reykjavik Study: multidisciplinary applied phenomics. Am J Epidemiol 165: 1076–1087. [PMC free article] [PubMed]
18. The ARIC Investigators (1989) The Atherosclerosis Risk in Communities (ARIC) Study: design and objectives. Am J Epidemiol 129: 687–702. [PubMed]
19. Fried LP, Borhani NO, Enright P, Furberg CD, Gardin JM, et al. (1991) The Cardiovascular Health Study: design and rationale. Ann Epidemiol 1: 263–276. [PubMed]
20. Mitchell P, Smith W, Attebo K, Wang JJ (1995) Prevalence of age-related maculopathy in Australia. The Blue Mountains Eye Study. Ophthalmology 102: 1450–1460. [PubMed]
21. Hofman A, van Duijn CM, Franco OH, Ikram MA, Janssen HL, et al. (2011) The Rotterdam Study: 2012 objectives and design update. Eur J Epidemiol 26: 657–686. [PMC free article] [PubMed]
22. Zheng Y, Lavanya R, Wu R, Wong WL, Wang JJ, et al. (2011) Prevalence and causes of visual impairment and blindness in an urban Indian population: the Singapore Indian Eye Study. Ophthalmology 118: 1798–1804. [PubMed]
23. Foong AW, Saw SM, Loo JL, Shen S, Loon SC, et al. (2007) Rationale and methodology for a population-based study of eye diseases in Malay people: The Singapore Malay eye study (SiMES). Ophthalmic Epidemiol 14: 25–35. [PubMed]
24. Bird AC, Bressler NM, Bressler SB, Chisholm IH, Coscas G, et al. (1995) An international classification and grading system for age-related maculopathy and age-related macular degeneration. The International ARM Epidemiological Study Group. Surv Ophthalmol 39: 367–374. [PubMed]
25. Klein R, Davis MD, Magli YL, Segal P, Klein BE, et al. (1991) The Wisconsin age-related maculopathy grading system. Ophthalmology 98: 1128–1134. [PubMed]
26. Aulchenko YS, Ripke S, Isaacs A, van Duijn CM (2007) GenABEL: an R library for genome-wide association analysis. Bioinformatics 23: 1294–1296. [PubMed]
27. Li Y, Willer C, Sanna S, Abecasis G (2009) Genotype imputation. Annu Rev Genomics Hum Genet 10: 387–406. [PMC free article] [PubMed]
28. Estrada K, Abuseiris A, Grosveld FG, Uitterlinden AG, Knoch TA, et al. (2009) GRIMP: a web- and grid-based tool for high-speed analysis of large-scale genome-wide association using imputed data. Bioinformatics 25: 2750–2752. [PMC free article] [PubMed]
29. Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, et al. (2010) LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 26: 2336–2337. [PMC free article] [PubMed]
30. de Bakker PI, Ferreira MA, Jia X, Neale BM, Raychaudhuri S, et al. (2008) Practical aspects of imputation-driven meta-analysis of genome-wide association studies. Hum Mol Genet 17: R122–128. [PMC free article] [PubMed]
31. Clogg CC, Petkova E, Haritou A (1995) Statistical methods for comparing regression coefficients between models. American Journal of Sociology 100: 1261–1293.
32. Francis PJ, Hamon SC, Ott J, Weleber RG, Klein ML (2009) Polymorphisms in C2, CFB and C3 are associated with progression to advanced age related macular degeneration associated with visual loss. J Med Genet 46: 300–307. [PubMed]
33. Seddon JM, Francis PJ, George S, Schultz DW, Rosner B, et al. (2007) Association of CFH Y402H and LOC387715 A69S with progression of age-related macular degeneration. JAMA 297: 1793–1800. [PubMed]
34. Seddon JM, Reynolds R, Yu Y, Daly MJ, Rosner B (2011) Risk models for progression to advanced age-related macular degeneration using demographic, environmental, genetic, and ocular factors. Ophthalmology 118: 2203–2211. [PubMed]
35. Yu Y, Reynolds R, Rosner B, Daly MJ, Seddon JM (2012) Prospective assessment of genetic effects on progression to different stages of age-related macular degeneration using multistate Markov models. Invest Ophthalmol Vis Sci 53: 1548–1556. [PMC free article] [PubMed]
36. Riobo NA, Manning DR (2007) Pathways of signal transduction employed by vertebrate Hedgehogs. Biochem J 403: 369–379. [PubMed]
37. Muto A, Iida A, Satoh S, Watanabe S (2009) The group E Sox genes Sox8 and Sox9 are regulated by Notch signaling and are required for Muller glial cell development in mouse retina. Exp Eye Res 89: 549–558. [PubMed]
38. Bhatia B, Singhal S, Lawrence JM, Khaw PT, Limb GA (2009) Distribution of Muller stem cells within the neural retina: evidence for the existence of a ciliary margin-like zone in the adult human eye. Exp Eye Res 89: 373–382. [PubMed]
39. Raymond PA, Hitchcock PF (1997) Retinal regeneration: common principles but a diversity of mechanisms. Adv Neurol 72: 171–184. [PubMed]
40. Friedrich U, Myers CA, Fritsche LG, Milenkovich A, Wolf A, et al. (2011) Risk- and non-risk-associated variants at the 10q26 AMD locus influence ARMS2 mRNA expression but exclude pathogenic effects due to protein deficiency. Hum Mol Genet 20: 1387–1399. [PMC free article] [PubMed]
41. Wu KH, Madigan MC, Billson FA, Penfold PL (2003) Differential expression of GFAP in early v late AMD: a quantitative analysis. Br J Ophthalmol 87: 1159–1166. [PMC free article] [PubMed]
42. Ruiz i Altaba A (1999) Gli proteins encode context-dependent positive and negative functions: implications for development and disease. Development 126: 3205–3216. [PubMed]
43. Evans J (2008) Antioxidant supplements to prevent or slow down the progression of AMD: a systematic review and meta-analysis. Eye (Lond) 22: 751–760. [PubMed]
44. Schraermeyer U, Heimann K (1999) Current understanding on the role of retinal pigment epithelium and its pigmentation. Pigment Cell Res 12: 219–236. [PubMed]
45. Weiter JJ, Delori FC, Wing GL, Fitch KA (1985) Relationship of senile macular degeneration to ocular pigmentation. Am J Ophthalmol 99: 185–187. [PubMed]
46. Stokowski RP, Pant PV, Dadd T, Fereday A, Hinds DA, et al. (2007) A genomewide association study of skin pigmentation in a South Asian population. Am J Hum Genet 81: 1119–1132. [PubMed]
47. Sulem P, Gudbjartsson DF, Stacey SN, Helgason A, Rafnar T, et al. (2007) Genetic determinants of hair, eye and skin pigmentation in Europeans. Nat Genet 39: 1443–1452. [PubMed]
48. Nan H, Kraft P, Qureshi AA, Guo Q, Chen C, et al. (2009) Genome-wide association study of tanning phenotype in a population of European ancestry. J Invest Dermatol 129: 2250–2257. [PMC free article] [PubMed]
49. Jin Y, Birlea SA, Fain PR, Gowan K, Riccardi SL, et al. (2010) Variant of TYR and autoimmunity susceptibility loci in generalized vitiligo. N Engl J Med 362: 1686–1697. [PMC free article] [PubMed]
50. Barrett JH, Iles MM, Harland M, Taylor JC, Aitken JF, et al. (2011) Genome-wide association study identifies three new melanoma susceptibility loci. Nat Genet 43: 1108–1113. [PMC free article] [PubMed]
51. Bishop DT, Demenais F, Iles MM, Harland M, Taylor JC, et al. (2009) Genome-wide association study identifies three loci associated with melanoma risk. Nat Genet 41: 920–925. [PMC free article] [PubMed]
52. Hu DN, Simon JD, Sarna T (2008) Role of ocular melanin in ophthalmic physiology and pathology. Photochem Photobiol 84: 639–644. [PubMed]
53. Memoli S, Napolitano A, d'Ischia M, Misuraca G, Palumbo A, et al. (1997) Diffusible melanin-related metabolites are potent inhibitors of lipid peroxidation. Biochim Biophys Acta 1346: 61–68. [PubMed]
54. Sundelin SP, Nilsson SE, Brunk UT (2001) Lipofuscin-formation in cultured retinal pigment epithelial cells is related to their melanin content. Free Radic Biol Med 30: 74–81. [PubMed]
55. Ma L, Dou HL, Wu YQ, Huang YM, Huang YB, et al. (2012) Lutein and zeaxanthin intake and the risk of age-related macular degeneration: a systematic review and meta-analysis. Br J Nutr 107: 350–359. [PubMed]
56. Baird PN, Guida E, Chu DT, Vu HT, Guymer RH (2004) The epsilon2 and epsilon4 alleles of the apolipoprotein gene are associated with age-related macular degeneration. Invest Ophthalmol Vis Sci 45: 1311–1315. [PubMed]
57. Klaver CC, Kliffen M, van Duijn CM, Hofman A, Cruts M, et al. (1998) Genetic association of apolipoprotein E with age-related macular degeneration. Am J Hum Genet 63: 200–206. [PubMed]
58. Schmidt S, Klaver C, Saunders A, Postel E, De La Paz M, et al. (2002) A pooled case-control study of the apolipoprotein E (APOE) gene in age-related maculopathy. Ophthalmic Genet 23: 209–223. [PubMed]
59. Thakkinstian A, Bowe S, McEvoy M, Smith W, Attia J (2006) Association between apolipoprotein E polymorphisms and age-related macular degeneration: A HuGE review and meta-analysis. Am J Epidemiol 164: 813–822. [PubMed]
60. Zareparsi S, Reddick AC, Branham KE, Moore KB, Jessup L, et al. (2004) Association of apolipoprotein E alleles with susceptibility to age-related macular degeneration in a large cohort from a single center. Invest Ophthalmol Vis Sci 45: 1306–1310. [PubMed]
61. Harold D, Abraham R, Hollingworth P, Sims R, Gerrish A, et al. (2009) Genome-wide association study identifies variants at CLU and PICALM associated with Alzheimer's disease. Nat Genet 41: 1088–1093. [PMC free article] [PubMed]
62. Heinzen EL, Need AC, Hayden KM, Chiba-Falek O, Roses AD, et al. (2010) Genome-wide scan of copy number variation in late-onset Alzheimer's disease. J Alzheimers Dis 19: 69–77. [PMC free article] [PubMed]
63. Kim S, Swaminathan S, Shen L, Risacher SL, Nho K, et al. (2011) Genome-wide association study of CSF biomarkers Abeta1–42, t-tau, and p-tau181p in the ADNI cohort. Neurology 76: 69–79. [PMC free article] [PubMed]
64. Naj AC, Beecham GW, Martin ER, Gallins PJ, Powell EH, et al. (2010) Dementia revealed: novel chromosome 6 locus for late-onset Alzheimer disease provides genetic evidence for folate-pathway abnormalities. PLoS Genet 6. [PMC free article] [PubMed]
65. Seshadri S, Fitzpatrick AL, Ikram MA, DeStefano AL, Gudnason V, et al. (2010) Genome-wide analysis of genetic loci associated with Alzheimer disease. JAMA 303: 1832–1840. [PMC free article] [PubMed]
66. Aulchenko YS, Ripatti S, Lindqvist I, Boomsma D, Heid IM, et al. (2009) Loci influencing lipid levels and coronary heart disease risk in 16 European population cohorts. Nat Genet 41: 47–55. [PMC free article] [PubMed]
67. Middelberg RP, Ferreira MA, Henders AK, Heath AC, Madden PA, et al. (2011) Genetic variants in LPL, OASL and TOMM40/APOE-C1-C2-C4 genes are associated with multiple cardiovascular-related traits. BMC Med Genet 12: 123. [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science