|Home | About | Journals | Submit | Contact Us | Français|
Revised the article for important intellectual content: EGH AVS BKC GHSB RAJ XS T. Aspelund T. Aung PNB EB CYC CMD GE VG TH AWH MI FJ BEKK LL XL GL TL P. McElduff BM PM BMP ER JIR RJS WT KT YYT AGU AV SX WTCCC2 JRV CCWK EST DS RK MFC TYW JA JJW. Approved final version to be published: EGH AVS BKC GHSB RAJ XS T. Aspelund T. Aung PNB EB CYC CMD GE VG TH AWH MI FJ BEKK LL XL GL TL P. McElduff BM P. Mitchell BMP ER JIR RJS WT KT YYT AGU AV SX WTCCC2 JRV CCWK EST DS RK MFC TYW JA JJW. Conceived and designed the experiments: JRV CCWK EST DS RK MFC TYW JA JJW. Performed the experiments: EGH AVS BKC GHSB RAJ XS. Analyzed the data: EGH AVS BKC GHSB RAJ XS. Contributed reagents/materials/analysis tools: T. Aspelund T. Aung PNB EB CYC CMD GE VG TH AWH MI FJ BEKK LL XL GL TL P. McElduff BM P. Mitchell BMP ER JIR RJS WT KT YYT AGU AV SX WTCCC2. Wrote the paper: EGH JA JJW.
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.
Age-related macular degeneration (AMD) is the most common cause of irreversible blindness in the elderly in many countries . Clinically, the disease is believed to develop via a series of progressive stages . 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 . 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 . 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 .
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 , . 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 . For example, heritability estimates from twin studies are 0.45 for early AMD , but 0.71 for late AMD .
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) , , , , , . 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 , , , , , .
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 .
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.
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) , the Atherosclerosis Risk in Communities (ARIC) study , the Cardiovascular Health Study (CHS) , the Blue Mountains Eye Study (BMES)  and three distinct cohorts from the Rotterdam Study (RS) : 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)  and the Singapore Malay Eye Study (SiMES) . 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.
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.
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 . 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 , 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.
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.
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  for AGES and ARIC cohorts, mach2dat  for BMES, R software for CHS and GRIMP software for the RS cohorts .
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 . Meta-analysis genomic control inflation factors (λGC) were calculated as previously described , as were standardised values for a sample of 1000 cases and 1000 controls (λ1000) .
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  formulated as
where , and , 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.
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.
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).
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.
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.
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.
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.
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) , , , . 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 , , , , , 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 . Notably, Shh signalling is required for the development of Muller glial stem cells , a retinal progenitor cell (RPC) with both proliferative and neurogenic ability . 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 (24). Muller glia isolated from the adult human retinal marginal region can also re-enter the cell cycle and demonstrate proliferative potential in vitro . The presence of a Muller cell-specific regulatory region upstream of the HTRA1 promoter , and the upregulation of glial fibrillary acidic protein (GFAP) immunoreactivity in the Müller cells of donor retina with drusen  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 , as zinc supplements have been shown reduce the risk of progression from early to late AMD .
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  and deficiencies of melanin in RPE cells have been associated with AMD , . Variants in TYR, including SNPs correlated with our lead SNP, have previously shown strong association with skin, hair and eye colour , , tanning ability , vitiligo  and melanoma risk ,  via GWAS. In addition to protecting from sunlight-induced damage , melanin is also an efficient antioxidant, reducing the oxidative stress resulting from lipid peroxidation and reactive oxygen species (ROS) generation . Specifically, elevated melanin content in RPE cells has been shown to reduce accretion of the auto-oxidant lipofuscin , a key waste product of photoreceptor outer segment phagocytosis implicated in extracellular drusen formation . This antioxidant function of melanin is consistent with evidence supporting supplemental lutein and zeaxanthin in modifying the course of AMD .
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 , , , , . Our most strongly associated SNPs have also previously shown unequivocal association with Alzheimer’s disease diagnosis, pathologic features and biomarkers , , , , , as well as lipid levels and cardiovascular traits ,  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) . 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.
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/
Supplementary Figures S1– S5.
Supplementary Tables S1– S3.
Supplementary Tables S4– S6.