Molecular Vision 2010; 16:2811-2821 <http://www.molvis.org/molvis/v16/a301>
Received 22 December 2009 | Accepted 13 December 2010 | Published 17 December 2010

Genome-wide association analyses of genetic, phenotypic, and environmental risks in the age-related eye disease study

Euijung Ryu,1 Brooke L. Fridley,1 Nirubol Tosakulwong,1 Kent R. Bailey,1 Albert O. Edwards2

1Department of Health Sciences Research, Mayo Clinic, Rochester, MN; 2Institute for Molecular Biology, University of Oregon, Eugene, OR

Correspondence to: Albert O. Edwards, Institute of Molecular Biology, 1229 University of Oregon, Eugene, OR, 97403-1229; Phone: (514) 346-2532; FAX: (541) 346-4854; email: edwardsa@uoregon.edu

Abstract

Purpose: To present genome-wide association analyses of genotypic and environmental risks on age-related macular degeneration (AMD) using 593 subjects from the age-related eye disease study (AREDS), after adjusting for population stratification and including questionable controls.

Methods: Single nucleotide polymorphism (SNP) associations with AMD for the non-Hispanic white population were investigated using a log-additive model after adjusting for population stratification. Replication of possible SNP-disease association was performed by genotyping an independent group of 444 AMD case and 300 control subjects. Logistic regression models were used to assess interaction effects between smoking and SNPs associated with AMD. Independent genetic risk effects among the disease-associated SNPs were also investigated using multiple logistic regression models.

Results: Population stratification was observed among the individuals having a self-reported race of non-Hispanic white. Risk allele frequencies at established AMD loci demonstrated that questionable control subjects were similar to control subjects in the AREDS, suggesting that they could be used as true controls in the analyses. Genetic loci (complement factor H [CFH], complement factor B [CFB], the age-related maculopathy susceptibility 2 locus containing the hypothetical gene [LOC387715]/the high-temperature requirement A-1 [HTRA1], and complement component 3 [C3]) that were already known to be associated with AMD were identified. An additional 26 novel SNPs potentially associated with AMD were identified, but none were definitely replicated in a second independent group of subjects. Smoking did not interact with known AMD loci, but was associated with late AMD. Statistically independent genetic signals were observed within the Pleckstrin homology domain-containing family A member 1 (PLEKHA1) region near LOC387715/HTRA1 and within a haplotype spanning exon 19 of the C3 gene.

Conclusions: Population stratification among Caucasian subjects from the multicentered AREDS was observed, suggesting that it should be adjusted for in future studies. The AREDS questionable control subjects can be used as control subjects in the AREDS genome-wide association study (GWAS). Smoking was an independent risk factor for advanced AMD in the AREDS subjects. There continues to be evidence that the 10q26 (age-related maculopathy susceptibility 2 gene [ARMS2]) locus spanning PLEKHA1-LOC387715-HTRA1 and the C3 gene may contain multiple independent genetic risks contributing to AMD.

Introduction

Age-related macular degeneration (AMD) is a progressive neurodegenerative disease that is the leading cause of vision loss in older individuals [1]. The more advanced forms of AMD are associated with atrophy of the retinal pigment epithelium (geographic atrophy) or abnormal angiogenesis (exudation) [2]. AMD is inherited as a complex disease with multiple genetic and environmental risks [3]. Environmental and lifestyle risk factors that may modify an individual’s risk of developing AMD include smoking, alcohol consumption, obesity, a lower level of physical activity, and diet [4]. The strongest identifiable risk factors for AMD continue to be age, family history, genetics, and smoking [5].

Currently, two biologic pathways contain genetic variants contributing to AMD risk. Genetic variations in multiple complement genes [6-12] and in a locus on chromosome 10q26 [13] are major risks for the development of AMD. Variants in other regions such as the human leukocyte antigen (HLA) locus, mitochondrial genome, interleukin 8 (IL8), apolipoprotein E (ApoE), and factor I (FI), may alter the risk of AMD, but additional studies are needed due to inconsistent or absent replication in independent study groups [14-17].

The genetic epidemiology of AMD is currently being investigated using genome-wide association studies (GWASs). The 593 subjects from the age-related eye disease study (AREDS) were genotyped using the Illumina Human-1 Bead-Chip panel containing 109,270 single nucleotide polymorphisms (SNPs). The data set was deposited into the National Center for Biotechnology Information database of Genotypes and Phenotypes (dbGaP). The AREDS subjects have been used for replication studies [11,18-20], meta-analyses [21], gentic epidemiology [22-24], and candidate gene selection by us and other investigators [12,21,25,26], but an analysis of the GWAS that includes strict quality control and population stratification has not yet been presented.

The AREDS was initiated in 1990 to learn more about the clinical course, prognosis, and risk factors of AMD and cataract. Over 4,700 people, ages 55 to 80 years, were enrolled in the AREDS cohort. These individuals were randomized for treatment with placebo, antioxidants, zinc, or antioxidants plus zinc, and were followed for a minimum of five years [27]. Results of the clinical trial suggested that people with a high risk of developing advanced AMD lowered their risk of AMD by 25% when they were treated with a high dose combination of vitamin C, vitamin E, beta-carotene, and zinc [27]. The AREDS GWAS is a valuable data set because of the extensive phenotypic and epidemiological variables collected during the study.

We present a comprehensive analysis of the AREDS GWAS to search for novel genetic loci for AMD, their interactions with smoking, and multiple independent genetic risks within AMD loci. In doing so, we have applied genotype quality-control filters and assessed genetic association using additive models accounting for population stratification. We have also compared the questionable control group (see Methods section) to the control group to determine if they have allele frequencies at AMD risk loci that are more similar to true controls or to cases. We attempted replication of all SNP effects that exceeded a significance level (p value) of less than 0.0001 in an independent case-control study of AMD, using logistic regression under an additive genetic model.

Methods

AREDS subjects and GWAS

We downloaded genotype data from the AREDS database (NIH-dbGaP), controlled through dbGaP accession number phs000001.v1.p1. Of the 593 subjects included in the AREDS genome-wide scan, 395 cases and 198 controls were successfully genotyped. The distribution of AMD subtypes, control subtypes, smoking exposure, and race are shown in Table 1. The cases were graded using standard clinical definitions of AMD [27]. The control category included subjects with AMD category 1 in both eyes at all visits. AMD category 1 was defined as any drusen present being less than 63 μm and the total area of drusen being less than the area of a 125 μm circle. The questionable control categories ranged from 1 to 4. The questionable control category 1 is for subjects who have AMD category 2 at only one visit except the last, and the questionable control category 4 is for those who have AMD category 2 (small to intermediate drusen) in one eye at the last vist, but AMD category 1 for both eyes at all other visits.

Replication

SNPs from the AREDS genome-wide scan were analyzed using logistic regression under an additive genetic model. Those with p values of less than 0.0001 from loci not known to contribute to AMD were considered for replication using an independent sample of AMD cases and controls. The replication subjects consisted of 744 individuals including 444 AMD cases and 300 controls without AMD. The use of these subjects was approved by the institutional review board of the Mayo Clinic, written informed consent was obtained from all subjects, and the study was performed in accordance with the Helsinki declaration. Diagnosis was determined by review of fundus photographs as previously described [6]. Briefly, all subjects diagnosed with AMD had large drusen or more advanced findings, and controls had five or fewer hard drusen without pigment changes or more advanced findings. SNPs for replication were genotyped using either Illumina or TaqMan assays (Applied Biosystem, Foster City, CA) as previously described [11,21,26]. A log-additive model using logistic regression was used for single SNP association analysis.

Quality-control filters

Quality-control filters were applied to the SNPs from the AREDS GWAS before the assessment of association with AMD. SNPs with a call rate less than 95%, HWE p values based on controls <1×10−7, and SNPs with a minor allele frequency less than 1% were removed from further analyses. Among 593 subjects, three subjects were further excluded, as they were considered possible first-degree relatives to the remaining subjects, based on pairwise identities determined by state calculations. Ninety-five percent of the SNPs passed the quality-control filters, leaving 103,895 SNPs for analysis.

Population stratification

Assessment of population stratification was performed using a principal components method [28] with the SNPs that passed our quality-control thresholds across the genome found on the 100K Illumina SNP panel. Principal component analysis was first completed using all cases and controls for all races (n=590), followed by an analysis including only subjects whose self-reported race was white, non-Hispanic (n=576). All further analyses were corrected for the first two principal components.

Genome-wide analysis

Single SNP analyses were performed using R, PLINK, and SAS version 9.1 (SAS Institute Inc., Cary, NC). Log-additive models using logistic regression were performed to examine the association between single SNPs and AMD status, as this protected against finding positive associations due only to differences in the heterozygote genotype frequencies between the AMD cases and controls [26]. To visualize the results, quantile-quantile (Q-Q) plots were drawn to compare the distribution of observed p values to the expected distribution under the null hypothesis of no association. Plots of the level of significance of association with AMD status versus chromosomal position (Manhattan plots) were presented to show the significance of the SNPs across genomes.

Correction for known SNPs

We wished to determine if any SNPs would be associated with AMD after adjusting for known loci. The process of choosing SNPs to exclude known loci was as follows: First, we calculated single SNP p values for all SNPs in each known locus (complement factor H [CFH], complement factor I [CFI], complement component 2/complement factor B [C2/CFB], age-related maculopathy susceptibility 2 [ARMS2], and complement component 3 [C3] loci) after adjusting for the first two principal components, and we found the most significant SNP for each locus. Second, we calculated the extent of linkage disequilibrium (using r2) between the most significant SNP in each locus and the rest of the SNPs in that locus to determine if there were any significant SNPs not in high linkage disequilibrium (r2 less than 0.2). Third, we performed logistic regression to see if such SNPs had any additional effect after conditioning on the most significant SNP for each locus. Seven SNPs were thus obtained, statistically accounting for the genetic risk contributed to AMD by the known loci as follows: the CFH locus, rs1853883; CFI locus, rs6848178; C2/CFB locus, rs541862; ARMS2 locus, rs2014307 and rs4565845; and C3 locus, rs2230199 and rs433594. The two SNPs from the ARMS2 locus (and the C3 locus) were selected to represent that locus, as they were in low linkage disequilibrium (r2<0.2), they statistically significantly contributed to AMD, and they each made an independent contribution to AMD conditioning on the other SNP. To understand the independent genetic effects within the ARMS2 locus (and C3 locus), haplotype analyses were performed for using the score test implemented in R package haplo.stats [29]. We then repeated the GWAS analysis, performing single SNP analyses with adjustments for these seven SNPs and the first two principal components.

Smoking interaction

For assessment of the interaction effect of SNPs significantly associated with AMD and smoking status, a logistic regression model under a log-additive genetic model was used, in which each significant SNP and smoking status were included in the main effects. An interaction effect for SNP smoking was also included, after adjusting for the first two principal components. Smoking status was coded as 1 if the subject was a former smoker or current smoker at the beginning of the study, or was coded as 0 if the subject had never smoked. The distribution of smokers is shown in Table 1.

Results

Questionable controls had greater similarity to controls than to cases

The risk allele frequencies in the 27 questionable control subjects were compared to the other control and AMD subtypes using a two-sample test for equality of proportions (Figure 1). The risk allele frequencies in the questionable control group were very similar to those from the control group, but not from the cases. P values comparing questionable controls to control groups varied from 0.48 to 1, while those comparing them to cases varied from 1.04×10−5 to 0.06. Thus, the questionable controls could be combined with the controls for subsequent analyses.

Population stratification

To explore possible population stratification, a principal components approach was used to assess population stratification in all subjects and in those who self-reported their race as white, non-Hispanic (Table 1). Figure 2 shows a plot of the first two principal components, which discriminated blacks from non-Hispanic whites, with the exception of one self-defined black subject. The two individuals with unspecified race clustered close to the non-Hispanic whites. The figure shows that there was potential population stratification among the white, non-Hispanic subjects. For further analyses, we used subjects who self-reported as white, and corrected for population stratification among the white subjects.

Genome-wide analysis

The association between each SNP and AMD risk was investigated under the log-additive genetic model. Figure 3 shows Q-Q plots of the association p values under four different scenarios. The Q-Q plot using all subjects (3A) clearly shows population stratification with an unacceptably high genomic inflation factor of 1.23. Using only white subjects and after adjusting for the first two principal components (Figure 3C), the genomic inflation factor was reduced to an acceptable level of 1.014. Figure 3D shows that there is evidence for additional SNPs associated with AMD after correction for SNPs already known to be associated with AMD. All seven SNPs were incorporated into a logistic regression model to correct for known association. Figure 4 shows the distribution of p values across the genome; individual SNPs with p values less than 0.00001 are labeled, and additional information is provided in Table 2.

Association of individual SNPs with AMD

SNPs associated with AMD at a p value of less than 0.0001 before or after adjustment for SNPs already known to be associated with AMD are shown in Table 2. Note that we did not include any SNPs from known loci. There were 26 SNPs meeting these criteria before adjustment and nine SNPs after adjustment; six SNPs met both criteria. Replication in an independent group of subjects was attempted for all SNPs. Only one SNP reached nominal significance (CTRB locus) in the replication subject group, but the p value (0.02) would not be considered significant after Bonferroni correction for multiple testing. The value of incorporating population stratification and using the log-additive model is illustrated by the SNP (rs2230199) in C3 in Table 2, which reached the statistical significance threshold (p<0.0001) only after correction of population stratification and using a log-additive genetic model (Fisher’s exact test; p=0.00048).

Interaction between smoking and SNPs associated with AMD

Smoking was not associated with early AMD (OR=0.58, 95% CI=0.18–1.80. The p value equaled 0.34, but was associated with geographic atrophy (OR=1.62, 95% CI=1.07–2.44, p=0.02), exudative AMD (OR=1.51, 95% CI=1.00–2.26, p=0.05), and advanced AMD (OR=1.56, 95% CI=1.10–2.22, p=0.01) compared to control groups among the AREDS GWAS subjects. The replication sample showed similar results: early AMD (OR=0.87, 95% CI=0.60–1.25, p=0.45), geographic atrophy (OR=1.68, 95% CI=0.97–2.92, p=0.06), exudative AMD (OR=1.75, 95% CI=1.18–2.58, p=0.005), and advanced AMD (OR=1.73, 95% CI=1.22–2.46, p=0.002). All 59 SNPs associated with AMD at p values less than 0.0001 or already known to be associated were investigated for possible interactions with smoking status. Five SNPs (rs11208590, rs572515, rs7529589, rs12038333, and rs203674) reached nominal significance levels (with p values varying from 0.02 to 0.01), but none of the interactions would be considered significant after Bonferroni correction for multiple testing.

Independent genetic effects within the ARMS2 and C3 loci

We observed statistically independent effects of both rs4565845 and rs2014307 across the ARMS2 locus at chromosome 10q26 (see Correction for known SNPs in Methods section). SNP rs4565845 is located within the pleckstrin homology domain-containing family A member 1 (PLEKHA1) gene, while rs2014307 is located between LOC387715 and the high-temperature requirement A-1 (HTRA1) promoter. We inspected the haplotypes formed by the six SNPs genotyped in the AREDS GWAS (rs4146894, rs4752692, rs4565845, rs7097701, rs2292623, and rs2014307) across 62 kb of the 10q26 locus (PLEKHA1-LOC387715-HTRA1 promoter). The SNP most highly associated with AMD was rs2014307, which captures information from both risk and protective haplotypes across the ARMS2 locus. Formation of haplotypes with and without rs4752692 did not appear to lead to any refinement of the association of haplotypes with AMD. Specifically, the addition of rs4752692 split the haplotype formed by the other five SNPs (GC-TAC) that were highly associated with AMD (p=7.4×10−12) into two haplotypes (GCCTAC, 3.1×10−7 and GCATAC, 1.0×10−4) that were both associated with AMD.

We also observed independently, statistically significant effects of two SNPs (rs433594 and rs2230199) across the C3 locus. The non-synonymous SNP rs2230199 forms a high-risk haplotype with another non-synonymous SNP (rs1047286), and the effects of these two SNPs cannot be clearly distinguished [12]. SNP rs433594 is only modestly associated with AMD in the AREDS GWAS (p=0.01). However, four SNPs in this region near exon 19 (rs237554, rs423490, rs433594, and rs428453) form a haplotype (GGTC; with a frequency of 0.1 in cases and 0.18 in controls) that is associated with AMD (p=3×10−4) at a level similar to the non-synonymous SNP rs2230199 (3.4×10−5). The GGTC haplotype spanning exon 9 does not appear to be in linkage disequilibrium with the haplotype tagged by rs2230199 that spans exons 3–9.

Discussion

A genome-wide scan performed on 593 subjects from the AREDS was released in the database of Genotypes and Phenotypes (dbGaP) without an associated publication describing the genetic features of the data set. The AREDS genome-wide scan data set has been used extensively for genetic studies on AMD, as detailed in the introduction. Because of the ongoing analyses using the AREDS subjects and the AREDS genome-wide scan data set, we felt that a detailed analysis of the genome-wide scan data set might provide valuable guidance on its use in AMD research.

The AREDS data set has been studied previously in two reports [20,26]. These studies did not analyze the data set using the actual genotypes from the AREDS data set. Instead, both studies obtained p values from the AREDS dbGaP website to select the SNPs for replication studies. Herein, we have performed a primary analysis of the AREDS data set using standard approaches, including log-additive models and correction for population stratification.

Prior to initiating the analysis of the genome-wide scan data set, we sought to determine the heterogeneity within the various categories of control subjects within the AREDS. The phenotypic range of drusen that has been allowed for control subjects in AMD studies has been controversial. Many groups use an extreme phenotype approach in which only a few hard drusen (without more advanced features such as pigment changes) are allowed in control subjects, based on fundus photography [6,30,31]. Other groups have used clinical examination without fundus photography [32]. Some investigators have argued that even the extreme phenotypes used by many groups are not sufficiently stringent, and that different criteria for controls could explain the failure of other investigators to replicate genetic association in the TLR3 locus [33,34]. However, the high prevalence of a few hard drusen in aged individuals using both clinical examination and fundus photography suggests that more stringent criteria may not be possible. For example, 96% of all subjects (age 43 and older) had at least one drusen at the baseline in the Beaver Dam Eye Study [35].

A gold standard for defining a control could be based on clinical criteria (e.g., rate of progression to advanced disease) or some combination of biomarkers contributing to disease risk. Although genetic variation accounting for approximately 80% of the risk for AMD has now been identified, this information was not available at the time that the AREDS AMD severity scales were developed [36,37]. For the AREDS genome-wide scan data set, photographic grading was performed at a reading center using standardized approaches, and longitudinal follow-up detailing rates of progression allowed us to directly ask if “questionable control” subjects were more similar to cases or to controls. Figure 1 shows that the questionable control subjects were indistinguishable from the controls, based on risk allele frequencies at known AMD loci. Further, the rate of progression to advanced AMD for questionable control subjects was extremely low [27]. Thus, by both genetic and clinical criteria, questionable control subjects appear to be similar to controls. Our observation would be most valid in older subjects (≥60 years), and probably would not apply to younger individuals with a few hard drusen that are known to be at high risk for incident AMD [38]. Note that this observation will be useful to a variety of researchers studying the extensive clinical and epidemiologic data sets from the AREDS clinical trial.

We evaluated the AREDS GWAS for problems that could lead to faulty conclusions in genetic association studies. Genetic heterogeneity can inflate the significance of the association between SNPs and disease (type I errors) [28,39-41]. It is thus standard practice to perform genetic association analyses on homogeneous racial populations defined by self-reported race. Racial classifications defined using genetic data may be even more accurate than self-reported race. For example, we observed individuals that did not fit genetically with the self-reported race (Figure 2). Further, we observed residual substructure within the white population that artificially inflated the overall association with AMD in the AREDS GWAS (Figure 2 and Figure 3). Principal components methods statistically corrected for the substructure, reducing the genomic inflation to an acceptable level (Figure 3). Our analyses suggest that investigators using genetic data from the AREDS trial should account for population substructure in their studies.

Another problem in GWAS studies is the observation that differences in heterozygote frequencies between cases and controls can lead to type I errors, or false positive associations [26]. Thus, we performed here an analysis using a log-additive genetic model that is less sensitive to genotype distributions that differ primarily by heterozygote frequencies, and that is more sensitive to additive changes that are often used for analyses of complex traits. However, we acknowledge the limitation of the log-additive genetic model, which can be less powerful if the true model is not additive. Another problem can arise from relying on statistical tests of allele frequency differences between subject groups. While genotype tests are independent of the Hardy–Weinberg equilibrium (HWE), allele tests require the calculation of allele frequencies assuming HWE [42]. Deviation from HWE can arise from inaccurate genotyping or copy-number variation. Thus, we favor statistical models of inheritance based on genotype frequencies.

We observed statistically independent effects of two SNPs (rs4565845 and rs2014307) across the ARMS2 locus at chromosome 10q26 and of two SNPs (rs433594 and rs2230199) across the C3 locus. The independent effect of rs4565845 that we observed in the PLEKHA1 gene adjacent to the ARMS2 locus does not seem consistent with a well defined genomic region (e.g., haplotype). Others also have reported independent association with AMD in the PLEKHA1 region, but it has not yet been possible to convincingly demonstrate an independent biologic effect, given the modest association with AMD, compared to the strong association in the haplotype block spanning the hypothetical gene LOC387715 and the promoter of HTRA1 [43]. It is hoped that these questions will be resolved in the future, as the biology of the ARMS2 locus becomes better understood.

On the other hand, the independent effect of rs433594 in the C3 gene appears to arise due to the strong association between AMD and a haplotype spanning exon 19 of C3. The genetic risk in the exon 19 region appears to be independent of the established genetic risks of non-synonymous SNPs in exons 3 and 9 of C3. Thus, it is possible that the genomic region spanning exon 19 contains a variant on a specific haplotype (see Independent genetic effects within the ARMS2 and C3 loci in Results section) that increases the risk of AMD. Re-sequencing studies are underway that may help address this question.

We identified 29 SNPs with p values greater than 0.0001 either before or after correction for known loci associated with AMD (Table 2). We attempted replication by genotyping each of these SNPs or a proxy (a nearby SNP in high linkage disequilibrium) in a larger, independent group of controls without AMD and in cases with AMD of similar ethnic background. However, none of the SNPs were definitely replicated in the independent group of subjects. Larger studies are ongoing and will be helpful in determining if the modest effects we observed upon replication (Table 2) and in our prior work with the AREDS GWAS may make a minor contribution to AMD [26]. Future work is needed to validate these results and to determine the functional relevance of these SNPs, if any.

Smoking is well established as a risk factor for both incident early AMD and advanced AMD [44-47]. In the AREDS GWAS, smoking was associated with both forms of advanced AMD only, presumably due to the stronger effect of smoking on the development of geographic atrophy and exudation. There have been inconsistent reports that the variants in the ARMS2 locus may mediate [48-50] or interact [51,52] with the effect of smoking, but we did not observe an interaction between smoking and genetic risks in the AREDS subjects. Smoking appears to be an independent risk factor contributing to the development of AMD.

In summary, we reported the value of correcting for population stratification when working with the AREDS GWAS and have provided evidence supporting the use of questionable control subjects as controls in the AREDS GWAS study. Further, we presented further evidence that the contribution of smoking to the development of AMD is independent of genetic risks. Finally, we observed that the AREDS data set has evidence suggestive of a previously unknown genetic risk near exon 19 of the C3 gene. Additional studies will be needed to determine if this observation can be replicated. There continues to be evidence for independent risks in the ARMS2 locus.

Acknowledgments

We thank the AREDS participants and the AREDS Research Group for their valuable contribution to this research. This work was supported by the National Eye Institute (R01EY014467), Bethesda, MD; the Foundation Fighting Blindness, Owing Mills, MD, and the American Health Assistance Foundation, Clarksburg, MD.

References

  1. Friedman DS, O'Colmain BJ, Muñoz B, Tomany SC, McCarty C, de Jong PT, Nemesure B, Mitchell P, Kempen J, Eye Diseases Prevalence Research Group.. Prevalence of age-related macular degeneration in the United States. Arch Ophthalmol. 2004; 122:564-72. [PMID: 15078675]
  2. Klein R, Davis MD, Magli YL, Segal P, Klein BE, Hubbard L. The Wisconsin age-related maculopathy grading system. Ophthalmology. 1991; 98:1128-34. [PMID: 1843453]
  3. Edwards AO, Malek G. Molecular genetics of AMD and current animal models. Angiogenesis. 2007; 10:119-32. [PMID: 17372852]
  4. Risk factors for neovascular age-related macular degeneration. The Eye Disease Case-Control Study Group. Arch Ophthalmol. 1992; 110:1701-8. [PMID: 1281403]
  5. Klein R, Peto T, Bird A, Vannewkirk MR. The epidemiology of age-related macular degeneration. Am J Ophthalmol. 2004; 137:486-95. [PMID: 15013873]
  6. Edwards AO, Ritter R, , 3rd Abel KJ, Manning A, Panhuysen C, Farrer LA. Complement factor H polymorphism and age-related macular degeneration. Science. 2005; 308:421-4. [PMID: 15761121]
  7. Haines JL, Hauser MA, Schmidt S, Scott WK, Olson LM, Gallins P, Spencer KL, Kwan SY, Noureddine M, Gilbert JR, Schnetz-Boutaud N, Agarwal A, Postel EA, Pericak-Vance MA. Complement factor H variant increases the risk of age-related macular degeneration. Science. 2005; 308:419-21. [PMID: 15761120]
  8. Klein RJ, Zeiss C, Chew EY, Tsai JY, Sackler RS, Haynes C, Henning AK, SanGiovanni JP, Mane SM, Mayne ST, Bracken MB, Ferris FL, Ott J, Barnstable C, Hoh J. Complement factor H polymorphism in age-related macular degeneration. Science. 2005; 308:385-9. [PMID: 15761122]
  9. Yates JR, Sepp T, Matharu BK, Khan JC, Thurlby DA, Shahid H, Clayton DG, Hayward C, Morgan J, Wright AF, Armbrecht AM, Dhillon B, Deary IJ, Redmond E, Bird AC, Moore AT, Genetic Factors in AMD Study Group.. Complement C3 variant and the risk of age-related macular degeneration. N Engl J Med. 2007; 357:553-61. [PMID: 17634448]
  10. Maller JB, Fagerness JA, Reynolds RC, Neale BM, Daly MJ, Seddon JM. Variation in complement factor 3 is associated with risk of age-related macular degeneration. Nat Genet. 2007; 39:1200-1. [PMID: 17767156]
  11. Edwards AO. Genetics of age-related macular degeneration. Adv Exp Med Biol. 2008; 613:211-9. [PMID: 18188947]
  12. Park KH, Fridley BL, Ryu E, Tosakulwong N, Edwards AO. Complement component 3 (C3) haplotypes and risk of advanced age-related macular degeneration. Invest Ophthalmol Vis Sci. 2009; 50:3386-93. [PMID: 19234341]
  13. Rivera A, Fisher SA, Fritsche LG, Keilhauer CN, Lichtner P, Meitinger T, Weber BH. Hypothetical LOC387715 is a second major susceptibility gene for age-related macular degeneration, contributing independently of complement factor H to disease risk. Hum Mol Genet. 2005; 14:3227-36. [PMID: 16174643]
  14. Fagerness JA, Maller JB, Neale BM, Reynolds RC, Daly MJ, Seddon JM. Variation near complement factor I is associated with risk of advanced AMD. Eur J Hum Genet. 2009; 17:100-4. [PMID: 18685559]
  15. Ennis S, Jomary C, Mullins R, Cree A, Chen X, Macleod A, Jones S, Collins A, Stone E, Lotery A. Association between the SERPING1 gene and age-related macular degeneration: a two-stage case-control study. Lancet. 2008; 372:1828-34. [PMID: 18842294]
  16. Goverdhan SV, Howell MW, Mullins RF, Osmond C, Hodgkins PR, Self J, Avery K, Lotery AJ. Association of HLA class I and class II polymorphisms with age-related macular degeneration. Invest Ophthalmol Vis Sci. 2005; 46:1726-34. [PMID: 15851575]
  17. Jones MM, Manwaring N, Wang JJ, Rochtchina E, Mitchell P, Sue CM. Mitochondrial DNA haplogroups and age-related maculopathy. Arch Ophthalmol. 2007; 125:1235-40. [PMID: 17846364]
  18. Ross RJ, Bojanowski CM, Wang JJ, Chew EY, Rochtchina E, Ferris FL, , 3rd Mitchell P, Chan CC, Tuo J. The LOC387715 polymorphism and age-related macular degeneration: replication in three case-control samples. Invest Ophthalmol Vis Sci. 2007; 48:1128-32. [PMID: 17325155]
  19. Conley YP, Jakobsdottir J, Mah T, Weeks DE, Klein R, Kuller L, Ferrell RE, Gorin MB. CFH, ELOVL4, PLEKHA1 and LOC387715 genes and susceptibility to age-related maculopathy: AREDS and CHS cohorts and meta-analyses. Hum Mol Genet. 2006; 15:3206-18. [PMID: 17000705]
  20. Zhang H, Morrison MA, Dewan A, Adams S, Andreoli M, Huynh N, Regan M, Brown A, Miller JW, Kim IK, Hoh J, Deangelis MM. The NEI/NCBI dbGAP database: genotypes and haplotypes that may predispose to risk of neovascular age-related macular degeneration. BMC Med Genet. 2008; 9:51 [PMID: 18541031]
  21. Edwards AO, Chen D, Fridley BL, James KM, Wu Y, Abecasis G, Swaroop A, Othman M, Branham K, Iyengar SK, Sivakumaran TA, Klein R, Klein BE, Tosakulwong N. Toll-like receptor polymorphisms and age-related macular degeneration. Invest Ophthalmol Vis Sci. 2008; 49:1652-9. [PMID: 18385087]
  22. Seddon JM, Francis PJ, George S, Schultz DW, Rosner B, Klein ML. Association of CFH Y402H and LOC387715 A69S with progression of age-related macular degeneration. JAMA. 2007; 297:1793-800. [PMID: 17456821]
  23. Francis PJ, George S, Schultz DW, Rosner B, Hamon S, Ott J, Weleber RG, Klein ML, Seddon JM. The LOC387715 gene, smoking, body mass index, environmental associations with advanced age-related macular degeneration. Hum Hered. 2007; 63:212-8. [PMID: 17347568]
  24. Klein ML, Francis PJ, Rosner B, Reynolds R, Hamon SC, Schultz DW, Ott J, Seddon JM. CFH and LOC387715/ARMS2 genotypes and treatment with antioxidants and zinc for age-related macular degeneration. Ophthalmology. 2008; 115:1019-25. [PMID: 18423869]
  25. Park KH, Ryu E, Tosakulwong N, Wu Y, Edwards AO. Common variation in the SERPING1 gene is not associated with age-related macular degeneration in two independent groups of subjects. Mol Vis. 2009; 15:200-7. [PMID: 19169411]
  26. Edwards AO, Fridley BL, James KM, Sharma AS, Cunningham JM, Tosakulwong N. Evaluation of clustering and genotype distribution for replication in genome wide association studies: the age-related eye disease study. PLoS ONE. 2008; 3:e3813 [PMID: 19043567]
  27. A randomized, placebo-controlled, clinical trial of high-dose supplementation with vitamins C and E, beta carotene, and zinc for age-related macular degeneration and vision loss: AREDS report no. 8. Arch Ophthalmol. 2001; 119:1417-36. [PMID: 11594942]
  28. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006; 38:904-9. [PMID: 16862161]
  29. Schaid DJ, Rowland CM, Tines DE, Jacobson RM, Poland GA. Score tests for association between traits and haplotypes when linkage phase is ambiguous. Am J Hum Genet. 2002; 70:425-34. [PMID: 11791212]
  30. Majewski J, Schultz DW, Weleber RG, Schain MB, Edwards AO, Matise TC, Acott TS, Ott J, Klein ML. Age-related macular degeneration–a genome scan in extended families. Am J Hum Genet. 2003; 73:540-50. [PMID: 12900797]
  31. Hageman GS, Hancox LS, Taiber AJ, Gehrs KM, Anderson DH, Johnson LV, Radeke MJ, Kavanagh D, Richards A, Atkinson J, Meri S, Bergeron J, Zernant J, Merriam J, Gold B, Allikmets R, Dean M, AMD Clinical Study Group.. Extended haplotypes in the complement factor H (CFH) and CFH-related (CFHR) family of genes protect against age-related macular degeneration: characterization, ethnic distribution and evolutionary implications. Ann Med. 2006; 38:592-604. [PMID: 17438673]
  32. Abecasis GR, Yashar BM, Zhao Y, Ghiasvand NM, Zareparsi S, Branham KE, Reddick AC, Trager EH, Yoshida S, Bahling J, Filippova E, Elner S, Johnson MW, Vine AK, Sieving PA, Jacobson SG, Richards JE, Swaroop A. Age-related macular degeneration: a high-resolution genome scan for susceptibility loci in a population enriched for late-stage disease. Am J Hum Genet. 2004; 74:482-94. [PMID: 14968411]
  33. Edwards AO, Swaroop A, Seddon JM. Geographic atrophy in age-related macular degeneration and TLR3. N Engl J Med. 2009; 360:2254-5. [PMID: 19469037]
  34. Ambati J, Katsanis N, Zhang K. Geograpic atrophy in age-related macular degeneration and TLR. N Engl J Med. 2009; 360:2255-6.
  35. Klein R, Klein BE, Linton KL. Prevalence of age-related maculopathy. The Beaver Dam Eye Study. Ophthalmology. 1992; 99:933-43. [PMID: 1630784]
  36. Hecker LA, Edwards AO, Ryu E, Tosakulwong N, Baratz KH, Brown WL, Charbel Issa P, Scholl HP, Pollok-Kopp B, Schmid-Kubista KE, Bailey KR, Oppermann M. Genetic control of the alternative pathway of complement in humans and age-related macular degeneration. Hum Mol Genet. 2010; 19:209-15. [PMID: 19825847]
  37. Chen W, Stambolian D, Edwards AO, Branham KE, Othman M, Jakobsdottir J, Tosakulwong N, Pericak-Vance MA, Campochiaro PA, Klein ML, Tan PL, Conley YP, Kanda A, Kopplin L, Li Y, Augustaitis KJ, Karoukis AJ, Scott WK, Agarwal A, Kovach JL, Schwartz SG, Postel EA, Brooks M, Baratz KH, Brown WL. Complications of Age-Related Macular Degeneration Prevention Trial Research Group, Brucker AJ, Orlin A, Brown G, Ho A, Regillo C, Donoso L, Tian L, Kaderli B, Hadley D, Hagstrom SA, Peachey NS, Klein R, Klein BE, Gotoh N, Yamashiro K, Ferris Iii F, Fagerness JA, Reynolds R, Farrer LA, Kim IK, Miller JW, Cortón M, Carracedo A, Sanchez-Salorio M, Pugh EW, Doheny KF, Brion M, Deangelis MM, Weeks DE, Zack DJ, Chew EY, Heckenlively JR, Yoshimura N, Iyengar SK, Francis PJ, Katsanis N, Seddon JM, Haines JL, Gorin MB, Abecasis GR, Swaroop A. Genetic variants near TIMP3 and high-density lipoprotein-associated loci influence susceptibility to age-related macular degeneration. Proc Natl Acad Sci USA. 2010; 107:7401-6. [PMID: 20385819]
  38. Klein R, Klein BE, Tomany SC, Meuer SM, Huang GH. Ten-year incidence and progression of age-related maculopathy: The Beaver Dam eye study. Ophthalmology. 2002; 109:1767-79. [PMID: 12359593]
  39. Devlin B, Roeder K. Genomic control for association studies. Biometrics. 1999; 55:997-1004. [PMID: 11315092]
  40. Teo YY. Common statistical issues in genome-wide association studies: a review on power, data quality control, genotype calling and population structure. Curr Opin Lipidol. 2008; 19:133-43. [PMID: 18388693]
  41. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000; 155:945-59. [PMID: 10835412]
  42. Schaid DJ, Jacobsen SJ. Biased tests of association: comparisons of allele frequencies when departing from Hardy-Weinberg proportions. Am J Epidemiol. 1999; 149:706-11. [PMID: 10206619]
  43. Jakobsdottir J, Conley YP, Weeks DE, Mah TS, Ferrell RE, Gorin MB. Susceptibility genes for age-related maculopathy on chromosome 10q26. Am J Hum Genet. 2005; 77:389-407. [PMID: 16080115]
  44. Paetkau ME, Boyd TA, Grace M, Bach-Mills J, Winship B. Senile disciform macular degeneration and smoking. Can J Ophthalmol. 1978; 13:67-71. [PMID: 647465]
  45. Hyman LG, Lilienfeld AM. Ferris FLd, Fine SL. Senile macular degeneration: a case-control study. Am J Epidemiol. 1983; 118:213-27. [PMID: 6881127]
  46. Seddon JM, Willett WC, Speizer FE, Hankinson SE. A prospective study of cigarette smoking and age-related macular degeneration in women JAMA. 1996; 276:1141-6.
    see comments
    [PMID: 8827966]
  47. Tomany SC, Wang JJ, Van Leeuwen R, Klein R, Mitchell P, Vingerling JR, Klein BE, Smith W, De Jong PT. Risk factors for incident age-related macular degeneration: pooled findings from 3 continents. Ophthalmology. 2004; 111:1280-7. [PMID: 15234127]
  48. Schmidt S, Hauser MA, Scott WK, Postel EA, Agarwal A, Gallins P, Wong F, Chen YS, Spencer K, Schnetz-Boutaud N, Haines JL, Pericak-Vance MA. Cigarette smoking strongly modifies the association of LOC387715 and age-related macular degeneration. Am J Hum Genet. 2006; 78:852-64. [PMID: 16642439]
  49. Mori K, Horie-Inoue K, Kohda M, Kawasaki I, Gehlbach PL, Awata T, Yoneya S, Okazaki Y, Inoue S. Association of the HTRA1 gene variant with age-related macular degeneration in the Japanese population. J Hum Genet. 2007; 52:636-41. [PMID: 17568988]
  50. Seddon JM, George S, Rosner B, Klein ML. CFH Gene Variant, Y402H, and Smoking, Body Mass Index, Environmental Associations with Advanced Age-Related Macular Degeneration. Hum Hered. 2006; 61:157-65. [PMID: 16816528]
  51. Baird PN, Robman LD, Richardson AJ, Dimitrov PN, Tikellis G, McCarty CA, Guymer RH. Gene-environment interaction in progression of AMD: the CFH gene, smoking and exposure to chronic infection. Hum Mol Genet. 2008; 17:1299-305. [PMID: 18203751]
  52. Wang JJ, Ross RJ, Tuo J, Burlutsky G, Tan AG, Chan CC, Favaloro EJ, Williams A, Mitchell P. LOC387715 polymorphism, inflammatory markers, smoking, and age-related macular degeneration. A population-based case-control study. Ophthalmology. 2008; 115:693-9. [PMID: 17675241]