Nature genetics. Apr/30/2011; 43(5): 436-441

Published online Apr/2/2011

PMID: 21460841

PMC: 3090745

doi: 10.1038/ng.801

Common variants in MS4A4/MS4A6E, CD2uAP, CD33, and EPHA1 are associated with late-onset Alzheimer’s disease

Badri Narayan Vardarajan

Jacqueline Buros

Paul J Gallins

Joseph D Buxbaum

Gail P Jarvik+146 authors

Abstract

The Alzheimer Disease Genetics Consortium (ADGC) performed a genome-wide association study (GWAS) of late-onset Alzheimer disease (LOAD) using a 3 stage design consisting of a discovery stage (Stage 1) and two replication stages (Stages 2 and 3). Both joint and meta-analysis analysis approaches were used. We obtained genome-wide significant results at MS4A4A [rs4938933; Stages 1+2, meta-analysis (P_M) = 1.7 × 10⁻⁹, joint analysis (P_J) = 1.7 × 10^−9; Stages 1–3, P_M = 8.2 × 10⁻¹²], CD2AP (rs9349407; Stages 1–3, P_M = 8.6 × 10⁻⁹), EPHA1 (rs11767557; Stages 1–3 P_M = 6.0 × 10⁻¹⁰), and CD33 (rs3865444; Stages 1–3, P_M = 1.6 × 10⁻⁹). We confirmed that CR1 (rs6701713; P_M = 4.6×10⁻¹⁰, P_J = 5.2×10⁻¹¹), CLU (rs1532278; P_M = 8.3 × 10⁻⁸, P_J = 1.9×10⁻⁸), BIN1 (rs7561528; P_M = 4.0×10⁻¹⁴; P_J = 5.2×10⁻¹⁴), and PICALM (rs561655; P_M = 7.0 × 10⁻¹¹, P_J = 1.0×10⁻¹⁰) but not EXOC3L2 are LOAD risk loci^1–3.

Alzheimer Disease (AD) is a neurodegenerative disorder affecting more than 13% of individuals aged 65 years and older and 30%–50% aged 80 years and older^4–5. Early work identified mutations in APP, PSEN1, and PSEN2 that cause early-onset autosomal dominant AD^6–9 and variants in APOE that affect LOAD susceptibility¹⁰. A recent GWAS identified CR1, CLU, PICALM, and BIN1 as LOAD susceptibility loci^1–3. However, because LOAD heritability estimates are high (h² ≈ 60–80%)¹¹, much of the genetic contribution remains unknown.

To identify genetic variants associated with risk for AD, the ADGC assembled a discovery dataset [Stage 1; 8,309 LOAD cases, 7,366 cognitively normal controls (CNEs)] using data from eight cohorts and a ninth newly assembled cohort from the 29 NIA-funded Alzheimer Disease Centers (ADCs) (Supplementary Tables 1 and 2, Supplementary Note) with data coordinated by the National Alzheimer Coordinating Center (NACC) and samples coordinated by the National Cell Repository for Alzheimer Disease (NCRAD). For the Stage 2 replication, we used four additional datasets and additional samples from the ADCs (3,531 LOAD cases, 3,565 CNEs). The Stage 3 replication used the results of association analyses provided by three other consortia (Hollingworth et al.¹²; 7,650 LOAD cases, 25,839 mixed-age controls). For Stages 1 and 2, we used both a meta-analysis (M) approach that integrates results from association analyses of individual datasets; and a joint analysis (J) approach where genotype data from each study are pooled. The latter method has improved power over meta-analysis in the absence of between-study heterogeneity¹³ and more direct correction for confounding sampling bias¹⁴. We were limited to meta-analysis for Stage 3.

Because cohorts were genotyped using different platforms, we used imputation to generate a common set of 2,324,889 SNPs. We applied uniform stringent quality control measures to all datasets to remove low-quality and redundant samples and problematic SNPs (Supplementary Tables 3, 4, and Online Methods). We performed association analysis assuming an additive model on the log odds ratio scale with adjustment for population substructure using logistic regression for case-control data and generalized estimating equations (GEE) with a logistic model for family data. Results from individual datasets were combined in the meta-analysis using the inverse variance method, applying a genomic control to each dataset. The joint analysis was performed using GEE and incorporated terms to adjust for population substructure and site-specific effects (Online Methods). For both approaches, we also examined an extended model of covariate adjustment that adjusted for age (age at onset or death in cases; age at exam or death in controls), sex, and number of APOE ε4 alleles (0, 1, or 2). Genomic inflation factors (λ) for both the discovery meta-analysis and the joint analysis and extended models were less than 1.05, indicating that there was not substantial inflation of the test statistics (Supplementary Table 3, Supplementary Figure 1). Association findings from meta-analysis and joint analysis were comparable.

In Stage 1, the strongest signal was from the APOE region (e.g., rs4420638, P_M =1.1 × 10⁻²⁶⁶, P_J =1.3 × 10⁻²⁵³; Supplementary Table 5). Excluding the APOE region, SNPs at nine distinct loci yielded a P_M or P_J ≤ 10⁻⁶ (Table 1; all SNPs with P < 10⁻⁴ are in Supplementary Table 5). SNPs from these nine loci were carried forward to Stage 2. Five of these had not previously been associated with LOAD at a genome-wide significance level of P ≤ 5.0 × 10⁻⁸ (MS4A, EPHA1, CD33, ARID5B, and CD2AP). Because Hollingworth et al.¹² identified SNPs at ABCA7 as a novel LOAD locus, we included ABCA7 region SNPs in Stage 2 and provided the results to Hollingworth et al.¹². For all loci in Table 1, we did not detect evidence for effect heterogeneity (Supplementary Fig. 2). One novel locus (MS4A) was significant in the Stage 1+2 analysis. Four other loci approached but did not reach genome-wide significance in the Stage 1+2 analyses and were carried forward to Stage 3. For three of these (CD33, EPHA1, and CD2AP), Stage 3 analysis strengthened evidence for association. However, Stages 2 and 3 results did not support Stage 1 results for ARID5B 2 (Table 2).

Stage 1+2 analysis identified the MS4A gene cluster as a novel LOAD locus (P_M = 1.7 × 10⁻⁹, P_J = 1.7 × 10⁻⁹)(Table 1, Fig. 1A). The minor allele (MAF = 0.39) was protective with identical odds ratios (ORs) from both meta-analysis and joint analysis (OR_M and OR_J = 0.88, 95% CI: 0.85–0.92). In the Stage 1+2 analysis, other SNPs gave smaller P values when compared to discovery SNP rs4938933, with the most significant SNP being rs4939338 (P_M = 2.6 × 10⁻¹¹, P_J = 4.6 × 10⁻¹¹; OR_M and OR_J = 0.87, 95% CI: 0.84–0.91) (Supplementary Table 5). In the accompanying manuscript¹², genome-wide significant results were also obtained at the MS4A locus (rs670139, P_M = 5.0 × 10⁻¹²) using an independent sample. In a combined analysis of ADGC results and those from Hollingworth et al.¹², the evidence for this locus at rs4938933 increased to P_M = 8.2 × 10⁻¹² (Table 3: OR_M = 0.89, 95% CI: 0.87–0.92; Fig. 1A).

SNPs in the CD2AP locus also met our Stage 1 criteria for additional analysis (Fig. 1B). Stage 2 data modestly strengthened this association, but the results did not reach genome-wide significance. Stage 3 analysis yielded a genome-wide significance result for rs9349407 (P_M = 8.6 × 10⁻⁹), identifying CD2AP as a novel LOAD locus. The minor allele (MAF = 0.27) at this SNP increased risk for LOAD (OR_M = 1.11, 95% CI: 1.07–1.15) (Table 2, Fig. 1B).

Another locus studied further in Stages 2 and 3 centered on EPHA1. Previous work provided suggestive evidence that this is a LOAD risk locus, although the associations did not reach genome-wide significance (P = 1.7 × 10⁻⁶)². Here, results from Stages 1 and 2 for SNP rs11767557, located in the promoter region of EPHA1, reached genome-wide significance in the joint analysis. The addition of Stage 3 results increased evidence for association (P_M = 6.0 × 10⁻¹⁰, Table 2, Fig. 1C). The minor allele (MAF = 0.19) for this SNP is protective (OR_M = 0.90, 95% CI: 0.86–0.93). We observed no evidence for heterogeneity at this locus (Supplementary Fig. 2D, heterogeneity P = 0.58).

In Stages 1 and 2, strong evidence for association was also obtained for SNPs in CD33, a gene located approximately 6Mb from APOE, but the results did not reach genome-wide significance. The addition of Stage 3 data confirmed that CD33 is a LOAD risk locus (rs3865444; Stages 1–3, P_M = 1.6 × 10⁻⁹). The minor allele (MAF = 0.30) is protective (OR_M = 0.91, 95% CI: 0.88–0.93; Tables 1,2, Fig. 1D). A single SNP (rs3826656) in the 5’ region of CD33, was previously reported as an AD-related locus using a family-based approach as genome-wide significant (P = 6.6 × 10⁻⁶) ¹⁵. We were unable to replicate this finding (P_M = 0.73; P_J = 0.39, Stage 1 analysis for rs3826656). Though rs3826656 is only 1,348 bp from our top SNP (rs3865444), these 2 sites display only weak LD (r² = 0.13).

Hollingworth et al¹² report highly significant evidence for the association of an ABCA7 SNP rs3764650 with LOAD (P_M = 4.5 × 10⁻¹⁷) that included data from our study. In our Stage 1+2 analysis, we obtained suggestive evidence for association with ABCA7 SNP rs3752246 (P_M = 5.8 × 10⁻⁷, and P_J = 5.0 × 10⁻⁷), which is a missense variant (G1527A) that may alter the function of the ABCA7 protein (see Supplementary Table 6 for functional SNPs in LD with SNPs yielding P_M or P_J < 10⁻⁴).

Our Stage 1+2 analyses also confirmed the association of previously reported loci (BIN1, CR1, CLU, and PICALM) with LOAD (Table 1). For each locus, supporting evidence was P ≤ 5.0 × 10⁻⁸ in one or both types of analysis.

We also examined SNPs with statistically significant GWAS results reported by others (GAB2¹⁶, PCDH11X¹⁷, GOLM1¹⁸, and MTHFD1L¹⁹, Supplementary Table 7). Stage 1 data were used except for PCDH11X where Stage 1+2 data were used because Affymetrix platforms do not contain the appropriate SNP. Only SNPs in the APOE, CR1, PICALM, and BIN1 loci demonstrated P < 10⁻⁶. For MTHFD1L¹⁹, at rs11754661 (previously reported P = 4.7 × 10⁻⁸) we obtained modest independent association evidence (OR_M = 1.16, 95% CI: 1.04–1.29, P_M = 0.006; OR_J = 1.19, 95% CI: 1.08–1.32, P_J = 7.5 × 10⁻⁴). For the remaining sites, only nominal evidence (P < 0.05) or no evidence was obtained. For the GAB2 locus¹⁶ at rs10793294 (previously reported P = 1.60 × 10⁻⁷), we obtained nominal statistical significance results (P_M = 0.017; P_J = 0.029). The association for rs5984894 in the PCDH11X locus¹⁷ (previously reported P = 3.9 × 10⁻¹²), did not replicate (P_M = 0.89, P_J = 0.26). Likewise, findings at GOLM1¹⁸ for rs10868366 (previously reported P = 2.40 × 10⁻⁴) did not replicate (P_M = 0.71; P_J = 0.62). Another gene consistently implicated in LOAD is SORL1²⁰ where at rs3781835 (previously reported P = 0.006), we obtained modest evidence for association (OR_M = 0.72, 95% CI: 0.60–0.86, P_M = 2.9 × 10⁻⁴; OR_J = 0.78, 95% CI: 0.59–0.86; P_J = 3.8 × 10⁻⁴).

We examined the influence of the APOE ε4 allele on the loci in Table 1, stratified by and in interactions with APOE ε4 allele carrier status. After adjustment, all loci had similar effect sizes to the unadjusted analyses with some showing a modest reduction in statistical significance. We previously reported evidence for a PICALM-APOE²¹ interaction using a dataset that largely overlaps with the Stage 1 dataset used here. However, using the Stage 1+2 data, we do not replicate this finding or see evidence of SNP-APOE interactions with Table 1 loci (data not shown).

Previous work reported an association between LOAD and chromosome 19 SNP rs597668, located 7.2 kb proximal to EXOC3L2 and 296 kb distal of APOE². While we did observe a signal for this SNP (Stage 1, P_M = 1.5 × 10⁻⁹; P_J = 7.7 × 10⁻¹⁰) and other SNPs in the EXOC2L3-MARK4 region, evidence was completely extinguished for all SNPs after adjustment for APOE (Online Methods, Supplementary Table 8), suggesting that signal in this region is from APOE.

Our observation of genome-wide significant associations at MS4A4A, CD2AP, EPHA1, and CD33 extend our understanding of the genetic architecture of LOAD and confirm the emerging consensus that common genetic variation plays a significant role in the etiology of LOAD. With our findings and those by Hollingsworth et al.¹², there are now ten LOAD susceptibility loci (APOE, CR1, CLU, PICALM, BIN1, EPHA1, MS4A, CD33, CD2AP, and ABCA7). Examining the amount of genetic effect attributable to these candidate genes, the most strongly associated SNPs at each locus other than APOE demonstrated population attributable fractions (PAFs) between 2.72–5.97% (Supplemental Table 9), with a cumulative PAF for non-APOE loci estimated to be as much as 35%; however, these estimates may vary widely between studies²², and the actual effect sizes are likely to be much smaller than those estimated here because of the ‘winner’s curse’. Also the results do not account for interaction among loci, and are not derived from appropriate population-based samples.

A recent review of GWAS studies²³ noted that risk alleles with small effect sizes (0.80 < OR < 1.2) likely exist for complex diseases such as LOAD but remain undetected, even with thousands of samples, because of insufficient power²⁴. Our discovery dataset (Stage 1; 8,309 cases and 7,366 controls), was well-powered to detect associations exceeding the statistical significance threshold of P < 10⁻⁶ (Supplementary Table 9). If there are many loci of more modest effects, some, but not all, will likely be detected in any one study. This likely explains the genome-wide statistical significance for the ABCA7 locus in the accompanying manuscript¹², which reaches only modest statistical significance in our dataset (rs3752246; P_M = 1.0 × 10⁻⁵, P_J = 1.9 × 10⁻⁵). Finding additional LOAD loci will require larger studies with increased depth of genotyping to test for the effects of both common and rare variants.

Supplementary Material

NIHMS287699-supplement-1.pdf

NIHMS287699-supplement-2.xls

NIHMS287699-supplement-3.xls

NIHMS287699-supplement-4.xls

NIHMS287699-supplement-01.pdf

Figure 1

Regional association plots from the three-stage meta-analysis with LOAD. P_M values for association are shown for: (A) MS4A gene cluster, (B) CD2AP, (C) EPHA1, and (D) CD33. For each locus, the genomic position (NCBI Build 37.1) is plotted on the X-axis against –log₁₀(P-value) on the Y-axis. For the SNP with the lowest P-value at each locus in Stage 1 analyses, three P-values for association are shown: P₁ meta-analysis of the ADGC Discovery (Stage 1) dataset (highlighted with a black diamond), P₁₊₂ meta-analysis of the Combined ADGC Discovery and Replication (Stages 1 + 2) datasets (highlighted with a blue diamond), and P₁₊₂₊₃ meta-analysis of the combined ADGC dataset and the external replication (Stages 1 + 2 + 3) datasets (highlighted with a red diamond). Computed estimates of linkage disequilibrium (r²) with the most significant SNP at each locus are shown as an orange diamond for r² ≥ 0.8, a yellow diamond for 0.5 ≤ r² < 0.8, a grey diamond for 0.2 ≤ r² < 0.5, and a white diamond for r² < 0.2. Genes in each region are indicated at the bottom of each panel. The length and the direction of the arrowhead represent the scaled size and the direction of the gene, respectively.

Table 1

Genome-wide Association Results for LOAD in the ADGC Stage 1 and Stage 2 datasets

Association signals represent SNPs with the strongest associations within each locus demonstrating P ≤ 10⁻⁶ in the Stage 1 dataset or in/near previously reported genes, excluding the APOE region (Supplementary Table 5).

SNP	CH:MB	Nearest Gene	MA	MAF	# SNPs	ADGC Discovery (Stage 1)				ADGC Replication (Stage 2)				Combined Analysis (Stages 1+2)
SNP	CH:MB	Nearest Gene	MA	MAF	# SNPs	OR_M (95% CI)	P_M	OR_J (95% CI)	P_J	OR_M (95% CI)	P_M	OR_J (95% CI)	P_J	OR_M (95% CI)	P_M	OR_J (95% CI)	P_J
rs6701713	1:207.8	CR1*	A	0.20	7	1.18 1.11–1.25	1.4×10⁻⁸	1.19 1.12–1.26	3.5×10⁻⁹	1.13 1.04–1.23	0.004	1.13 1.04–1.24	0.004	1.16 1.11–1.22	4.6×10⁻¹⁰	1.17 1.12–1.23	5.2×10⁻¹¹
rs7561528	2:127.9	BIN1*	A	0.35	10	1.18 1.13–1.24	2.9×10⁻¹¹	1.18 1.12–1.24	7.7×10⁻¹¹	1.15 1.07–1.24	1.4×10⁻⁴	1.15 1.07–1.24	1.0×10⁻⁴	1.17 1.13–1.22	4.2×10⁻¹⁴	1.17 1.12–1.22	5.2×10⁻¹⁴
rs9349407	6:47.5	CD2AP	C	0.27	1	1.14 1.08–1.21	1.2×10⁻⁶	1.14 1.08–1.20	5.3×10⁻⁶	1.07 0.98–1.17	0.118	1.08 0.99–1.18	0.074	1.12 1.07–1.18	1.0×10⁻⁶	1.12 1.07–1.17	2.1×10⁻⁶
rs11767557	7:143.1	EPHA1‡	C	0.19	1	0.85 0.80–0.90	7.3×10⁻⁸	0.84 0.79–0.89	3.1×10⁻⁸	0.94 0.86–1.03	0.169	0.93 0.85–1.02	0.133	0.87 0.83–0.92	2.4×10⁻⁷	0.87 0.83–0.91	4.9×10⁻⁸
rs1532278	8:27.5	CLU*	T	0.36	2	0.90 0.85–0.95	5.6×10⁻⁵	0.89 0.85–0.94	2.0×10⁻⁵	0.87 0.81–0.94	2.6×10⁻⁴	0.87 0.81–0.94	2.7×10⁻⁴	0.89 0.85–0.93	8.3×10⁻⁸	0.89 0.85–0.92	1.9×10⁻⁸
rs2588969	10:63.6	ARID5B	A	0.37	0	0.88 0.84–0.93	1.1×10⁻⁶	0.88 0.84–0.93	6.9×10⁻⁷	1.05 0.97–1.13	0.234	1.05 0.98–1.13	0.189	0.93 0.89–0.97	0.001	0.93 0.89–0.97	7.7×10⁻⁴
rs4938933	11:60.0	MS4A4A	C	0.39	22	0.88 0.84–0.92	5.2×10⁻⁸	0.87 0.83–0.92	4.5×10⁻⁸	0.90 0.84–0.97	0.005	0.90 0.84–0.97	0.004	0.88 0.85–0.92	1.7×10⁻⁹	0.88 0.85–0.92	1.7×10⁻⁹
rs561655	11:85.8	PICALM*	G	0.34	36	0.88 0.84–0.92	1.2×10⁻⁷	0.88 0.84–0.93	4.6×10⁻⁷	0.86 0.80–0.93	8.4×10⁻⁵	0.86 0.80–0.92	3.7×10⁻⁵	0.87 0.84–0.91	7.0×10⁻¹¹	0.87 0.84–0.91	1.0×10⁻¹⁰
rs3752246	19:1.1	ABCA7%	G	0.19	2	1.16 1.08–1.24	1.0×10⁻⁵	1.15 1.08–1.23	1.9×10⁻⁵	1.13 1.03–1.24	0.012	1.13 1.03–1.25	0.009	1.15 1.09–1.21	5.8×10⁻⁷	1.15 1.09–1.21	5.0×10⁻⁷
rs3865444	19:51.7	CD33#	A	0.30	1	0.88 0.84–0.93	8.2×10⁻⁷	0.88 0.84–0.93	1.9×10⁻⁶	0.91 0.85–0.99	0.021	0.92 0.85–0.99	0.029	0.89 0.86–0.93	1.1×10⁻⁷	0.89 0.86–0.93	2.0×10⁻⁷

CH:MB, chromosome:position (in mega base pairs, build 19); MA, minor allele; MAF, minor allele frequency; # SNPs, the number of SNPs for which P ≤ 1 × 10⁻⁶ in meta-analysis from the combined analysis in Stage 1+2; OR_M, odds ratio in meta-analysis; P_M, P-value in meta-analysis; OR_J, odds ratio in joint analysis; P_J, P-value in joint analysis.

Genes with previous case-control genome-wide statistically significant associations: CR1³, CLU ^1,3, PICALM¹, BIN1². Gene with a previous association not meeting genome-wide statistical significance: EPHA1². Family-based association study with reported genome-wide statistical significance: CD33¹⁵.

Genes with previously published case-control association signals at P ≤ 5.0 × 10⁻⁸ are denoted with *

the case-control locus that did not meet this level of statistical significance is denoted with ^‡

the locus previously reported in a family-based association study as genome-wide significant with ^#

locus identified in Hollingworth et al.¹² with genome-wide significant evidence for association with. ^%

Table 2

Meta-Analysis of Stage 1+2 with Stage 3 (CHARGE/GERAD/EADI1 Consortia ²) GWAS Results

Meta-analysis using an external replication case-control sample (Stage 3) for SNPs from novel loci at which associations did not exceed the genome-wide statistical significance threshold (P = 5.0 × 10⁻⁸) in the ADGC meta-analysis (Stage 1+2). Results for MS4A are also included to show association results from the ADGC and accompanying manuscript¹². The external replication dataset did not include results from TGEN, ADNI, and MAYO cohorts (Supplementary Tables 1 and 2).

Gene:SNP	Cases	Controls	Total	OR_M (95% CI)	P_M	OR_j (95% CI)	P_J
CD2AP: rs9349407
ADGC	11840	10931	22771	1.12 (1.07–1.18)	1.0 × 10⁻⁶	1.12 (1.07–1.17)	2.1 × 10⁻⁶
External	6922	18896	25818	1.09 (1.03–1.15)	0.002	-	-
ADGC + External	18762	29827	48589	1.11 (1.07–1.15)	8.6 × 10⁻⁹	-	-

EPHA1: rs11767557
ADGC	11840	10931	22771	0.87 (0.83–0.92)	2.4 × 10⁻⁷	0.87 (0.83–0.91)	4.9 × 10⁻⁸
External	6922	24666	31588	0.91 (0.87–0.96)	2.9 × 10⁻⁴	-	-
ADGC + External	18762	35597	54359	0.90 (0.86–0.93)	6.0 × 10⁻¹⁰	-	-

ARID5B: rs2588969
ADGC	11840	10931	22771	0.93 (0.89–0.97)	0.001	0.93 (0.89–0.97)	7.8 × 10⁻⁴
External	6922	18896	25818	1.06 (1.01–1.11)	0.018	-	-
ADGC + External	18762	29827	48589	0.99 (0.95–1.02)	0.362	-	-

MS4A4A: rs4938933
ADGC	11840	10931	22771	0.88 (0.85–0.92)	1.7 × 10⁻⁹	0.88 (0.85–0.92)	1.7 × 10⁻⁹
External	6922	18896	25818	0.92 (0.88–0.97)	5.4 × 10⁻⁴	-	-
ADGC + External	18762	29827	48589	0.89 (0.87–0.92)	8.2 × 10⁻¹²	-	-

CD33: rs3865444
ADGC	11840	10931	22771	0.89 (0.86–0.93)	1.1 × 10⁻⁷	0.89 (0.86–0.93)	2.0 × 10⁻⁷
External	6922	18896	25818	0.92 (0.88–0.97)	0.002	-	-
ADGC + External	18762	29827	48589	0.91 (0.88–0.93)	1.6 × 10⁻⁹	-	-

Footnotes

URLs.

The Alzheimer Disease Genetics Consortium (ADGC), http://alois.med.upenn.edu/adgc/about/overview.html; ADNI database, (www.loni.ucla.edu/ADNI); ADNI investigators, http://www.loni.ucla.edu/ADNI/Collaboration/ADNI_Manuscript_Citations.pdf; APOE Genotyping kit from TIB MOLBIOL, http://www.roche-as.es/logs/LightMix%C2%AE_40-0445-16_ApoE-112-158_V080904.pdf; PLINK, http://pngu.mgh.harvard.edu/~purcell/plink/; PREST, http://utstat.toronto.edu/sun/Software/Prest/; MACH, http://www.sph.umich.edu/csg/abecasis/mach/; EIGENSTRAT, http://genepath.med.harvard.edu/~reich/EIGENSTRAT.htm; The R Project for Statistical Computing, http://www.r-project.org/; Package GWAF in R, http://cran.r-project.org/web/packages/GWAF/index.html; Package gee in R, http://cran.r-project.org/web/packages/gee/index.html; UCSC Genome Browser, http://genome.ucsc.edu/; METAL, http://www.sph.umich.edu/csg/abecasis/Metal/; FUGUE, http://www.sph.umich.edu/csg/abecasis/fugue/.

Author Contributions

Sample collection, phenotyping, and data management: J.D.Buxbaum, G.P.J., P.K.C., E.B.L., T.D.B., B.F.B., N.R.G., P.L.D., D.E., J.A.Schneider, M.M.C., N.E., S.G.Y., C.C., J.S.K.K., P.N., P.K., J.H., M.J.H., A.J.M., M.M.B., F.Y.D., C.T.B., R.C.G., E.R., P.S.G.-H., S.E.A., R.B., T.B., E.H.B., J.D.Bowen, A.B., J.R.B., N.J.C., C.S.C., S.L.C., H.C.C., D.G.C., J.C., C.W.C., J.L.C., C.D., S.T.D., R.D.-A., M.D., D.W.D., W.G.E., K.M.F., K.B.F., M.R.F., S.F., M.P.F., D.R.G., M.Ganguli, M.Gearing, D.H.G., B.Ghetti, J.R.G., S.G., B.Giordani, J.G., J.H.G., R.L.H., L.E.H., E.H., L.S.H., C.M.H., B.T.H., G.A.J., L.-W.J., N.J., J.K., A.K., J.A.K., R.K., E.H.K., N.W.K., J.J.L., A.I.L., A.P.L., O.L.L., W.J.M., D.C.Marson, F.M., D.C.Mash, E.M., W.C.M., S.M.M., A.N.M., A.C.M., M.M., B.L.M., C.A.M., J.W.M., J.E.P., D.P.P., E.P., R.C.P., W.W.P., J.F.Q., M.R., B.R., J.M.R., E.D.R., R.N.R., M.S., L.S.S., W.S., M.L.S., M.A.S., C.D.S., J.A.Sonnen, S.S., R.A.S., R.E.T., J.Q.T., J.C.T., V.M.V., H.V.V., J.P.V., S.W., K.A.W., J.W., R.L.W., L.B.C., B.A.D., D.Beekly, M.I.K., A.J.S., E.M.R., D.A.B., A.M.G., W.A.K., T.M.F., J.L.H., R.M., M.A.P., L.A.F. Study management and coordination: L.B.C., D.Beekly, D.A.B., J.C.M., T.J.M., A.M.G., D.Blacker, D.W.T., H.H., W.A.K., T.M.F., J.L.H., R.M., M.A.P., L.A.F., G.D.S. Statistical methods and analysis: A.C.N., G.J., G.W.B., L.-S.W., B.N.V., J.B., P.J.G., R.M.C., R.A.R., M.A.S., K.L.L., E.R.M., J.L.H., M.A.P., L.A.F. Interpretation of results: A.C.N., G.J., G.W.B., L.-S.W., B.N.V., J.B., P.J.G., R.A.R., M.A.S., K.L.L., E.R.M., M.I.K., A.J.S., E.M.R., D.A.B., J.C.M., T.J.M., A.M.G., D.Blacker, D.W.T., H.H., W.A.K., T.M.F., J.L.H., R.M., M.A.P., L.A.F., G.D.S. Manuscript writing group: A.C.N., G.J., G.W.B., L.-S.W., B.N.V., J.B., P.J.G., J.L.H., R.M., M.A.P., L.A.F., G.D.S. Study design: D.A.B., J.C.M., T.J.M., A.M.G., D.Blacker, D.W.T., H.H., W.A.K., T.M.F., J.L.H., R.M., M.A.P., L.A.F., G.D.S.

Competing Financial Interests

T.D.B. received licensing fees from and is on the speaker's bureau of Athena Diagnostics, Inc. M.R.F. receives research funding from BristolMyersSquibb Company, Danone Research, Elan Pharmaceuticals, Inc., Eli Lilly and Company, Novartis Pharmaceuticals Corporation, OctaPharma AG, Pfizer Inc., and Sonexa Therapeutics, Inc; Receives honoraria as scientific consultant from Accera, Inc., Astellas Pharma US Inc., Baxter, Bayer Pharmaceuticals Corporation, BristolMyersSquibb, Eisai Medical Research, Inc., GE Healthcare, Medavante, Medivation, Inc., Merck & Co., Inc., Novartis Pharmaceuticals Corp., Pfizer, Inc., Prana Biotechnology Ltd., QR Pharma., Inc., The sanofi-aventis Group, and Toyama Chemical Co., Ltd.; and is speaker for Eisai Medical Research, Inc., Forest Laboratories, Pfizer Inc. and Novartis Pharmaceuticals Corporation. A.M.G. has research funding from AstraZeneca, Pfizer and Genentech, and has received remuneration for giving talks at Pfizer and Genentech. R.C.P. is on the Safety Monitory Committee of Pfizer, Inc. (Wyeth) and a consultant to the Safety Monitoring Committee at Janssen Alzheimer's Immunotherapy Program (Elan), to Elan Pharmaceuticals, and to GE Healthcare. R.E.T. is a consultant to Eisai, Japan in the area of Alzheimer's genetics and a shareholder in, and consultant to Pathway Genomics, Inc, San Diego, CA.

Acknowledgements

The National Institutes of Health, National Institute on Aging (NIH-NIA) supported this work through the following grants: ADGC, U01 AG032984, RC2 AG036528; NACC, U01 AG016976; NCRAD, U24 AG021886; NIA LOAD, U24 AG026395, U24 AG026390; Boston University, P30 AG013846, R01 HG02213, K24 AG027841, U01 AG10483, R01 CA129769, R01 MH080295, R01 AG009029, R01 AG017173, R01 AG025259; Columbia University, P50 AG008702, R37 AG015473; Duke University, P30 AG028377; Emory University, AG025688; Indiana University, P30 AG10133; Johns Hopkins University, P50 AG005146, R01 AG020688; Massachusetts General Hospital, P50 AG005134; Mayo Clinic, P50 AG016574; Mount Sinai School of Medicine, P50 AG005138, P01 AG002219; New York University, P30 AG08051, MO1RR00096, and UL1 RR029893; Northwestern University, P30 AG013854; Oregon Health & Science University, P30 AG008017, R01 AG026916; Rush University, P30 AG010161, R01 AG019085, R01 AG15819, R01 AG17917, R01 AG30146; University of Alabama at Birmingham, P50 AG016582, UL1RR02777; University of Arizona/TGEN, P30 AG019610, R01 AG031581, R01 NS059873; University of California, Davis, P30 AG010129; University of California, Irvine, P50 AG016573, P50, P50 AG016575, P50 AG016576, P50 AG016577; University of California, Los Angeles, P50 AG016570; University of California, San Diego, P50 AG005131; University of California, San Francisco, P50 AG023501, P01 AG019724; University of Kentucky, P30 AG028383; University of Michigan, P50 AG008671; University of Pennsylvania, P30 AG010124; University of Pittsburgh, P50 AG005133, AG030653; University of Southern California, P50 AG005142; University of Texas Southwestern, P30 AG012300; University of Miami, R01 AG027944, AG010491, AG027944, AG021547, AG019757; University of Washington, P50 AG005136, UO1 AG06781, UO1 HG004610; Vanderbilt University, R01 AG019085; and Washington University, P50 AG005681, P01 AG03991. ADNI Funding for ADNI is through the Northern California Institute for Research and Education by grants from Abbott, AstraZeneca AB, Bayer Schering Pharma AG, Bristol-Myers Squibb, Eisai Global Clinical Development, Elan Corporation, Genentech, GE Healthcare, GlaxoSmithKline, Innogenetics, Johnson and Johnson, Eli Lilly and Co., Medpace, Inc., Merck and Co., Inc., Novartis AG, Pfizer Inc, F. Hoffman-La Roche, Schering-Plough, Synarc, Inc., Alzheimer's Association, Alzheimer's Drug Discovery Foundation, the Dana Foundation, and by the National Institute of Biomedical Imaging and Bioengineering and NIA grants U01 AG024904, RC2 AG036535, K01 AG030514. We thank Creighton Phelps, Marcelle Morrison-Bogorad, and Marilyn Miller from NIA who are ex-officio ADGC members. Support was also from the Alzheimer’s Association (LAF, IIRG-08-89720; MP-V, IIRG-05-14147) and the Veterans Affairs Administration. P.S.G.-H. is supported by Wellcome Trust, Howard Hughes Medical Institute, and the Canadian Institute of Health Research.

References

1. HaroldDGenome-wide association study identifies variants at CLU and PICALM associated with Alzheimer's diseaseNat Genet2009411088-U61[Google Scholar]
2. SeshadriSGenome-wide analysis of genetic loci associated with Alzheimer diseaseJAMA201030318321840[PubMed][Google Scholar]
3. LambertJCGenome-wide association study identifies variants at CLU and CR1 associated with Alzheimer's diseaseNat Genet2009411094-U68[Google Scholar]
4. HebertLEScherrPABieniasJLBennettDAEvansDAAlzheimer disease in the US population - Prevalence estimates using the 2000 censusArch. Neurol20036011191122[PubMed][Google Scholar]
5. Alzheimer's Association, Alzheimer's Disease Facts and Figures2009Washington, D.C.
6. GoateASegregation of a missense mutation in the amyloid precursor protein gene with familial Alzheimer's diseaseNature1991349704706[PubMed][Google Scholar]
7. SherringtonRCloning of a gene bearing missense mutations in early-onset familial Alzheimer's diseaseNature1995375754760[PubMed][Google Scholar]
8. RogaevEIFamilial Alzheimer's disease in kindreds with missense mutations in a gene on chromosome 1 related to the Alzheimer's disease type 3 geneNature1995376775778[PubMed][Google Scholar]
9. Levy-LahadECandidate gene for the chromosome 1 familial Alzheimer's disease locusScience1995269973977[PubMed][Google Scholar]
10. CorderEHGene Dose of Apolipoprotein-E Type-4 Allele and the Risk of Alzheimer's Disease in Late Onset FamiliesScience1993261921923[PubMed][Google Scholar]
11. GatzMHeritability for Alzheimer's disease: The study of dementia in Swedish twinsJournals of Gerontology Series A - Biological Sciences and Medical Sciences199752M117M125[Google Scholar]
12. HollingworthPNature Genetics2011(in press)
13. SkolADScottLJAbecasisGRBoehnkeMJoint analysis is more efficient than replication-based analysis for two-stage genome-wide association studiesNat Genet200638209213[PubMed][Google Scholar]
14. IoannidisJPRosenbergPSGoedertJJO’BrienTRCommentary: meta-analysis of individual participants' data in genetic epidemiologyAm J Epidemiol2002156204210[PubMed][Google Scholar]
15. BertramLGenome-wide association analysis reveals putative Alzheimer's disease susceptibility loci in addition to APOEAm J Hum Genet200883623632[PubMed][Google Scholar]
16. ReimanEMGAB2 alleles modify Alzheimer's risk in APOE epsilon 4 carriersNeuron200754713720[PubMed][Google Scholar]
17. CarrasquilloMMGenetic variation in PCDH11X is associated with susceptibility to late-onset Alzheimer's diseaseNat Genet200941192198[PubMed][Google Scholar]
18. LiHCandidate single-nucleotide polymorphisms from a genomewide association study of Alzheimer diseaseArch. Neurol2008654553[PubMed][Google Scholar]
19. NajACBeechamGWMartinERGallinsPJPowellEHKonidariIWhiteheadPLCaiGHaroutunianVScottWKVanceJMSliferMAGwirtsmanHEGilbertJRHainesJLBuxbaumJDPericak-VanceMADementia Revealed: Novel Chromosome 6 Locus for Late-Onset Alzheimer Disease Provides Genetic Evidence for Folate-Pathway AbnormalitiesPlos Genetics20106[Google Scholar]
20. RogaevaEThe neuronal sortilin-related receptor SORL1 is genetically associated with Alzheimer diseaseNat. Genet200739168177[PubMed][Google Scholar]
21. JunGMeta-analysis Confirms CR1, CLU, and PICALM as Alzheimer Disease Risk Loci and Reveals Interactions With APOE GenotypesArch Neurol2010[Google Scholar]
22. RockhillBNewmanBWeinbergCUse and misuse of population attributable fractionsAm J Public Health1998881519[PubMed][Google Scholar]
23. KuCSLoyEYPawitanYChiaKSThe pursuit of genome-wide association studies: where are we now?J Hum Genet201055195206[PubMed][Google Scholar]
24. FlorezJCClinical review: the genetics of type 2 diabetes: a realistic appraisal in 2008J Clin Endocrinol Metab20089346334642[PubMed][Google Scholar]
25. PurcellSPLINK: A tool set for whole-genome association and population-based linkage analysisAm. J. Hum. Genet200781559575[PubMed][Google Scholar]
26. McPeekMSSunLStatistical tests for detection of misspecified relationships by use of genome-screen dataAm J Hum Genet20006610761094[PubMed][Google Scholar]
27. AbecasisGRChernySSCooksonWOCardonLRGRR: graphical representation of relationship errorsBioinformatics200117742743[PubMed][Google Scholar]
28. LiYAbecasisGRRapid haplotype reconstruction and missing genotype inferenceAmerican Journal of Human Genetics2006S792290[Google Scholar]
29. WittwerCTThe LightCycler: a microvolume multisample fluorimeter with rapid temperature controlBiotechniques199722176181[PubMed][Google Scholar]
30. AhmadianASingle-nucleotide polymorphism analysis by pyrosequencingAnal Biochem2000280103110[PubMed][Google Scholar]
31. HixsonJEVernierDTRestriction isotyping of human of human apolipoprotien E by gene amplification and cleavage with HhaIJ. Liped Res199031545548[Google Scholar]
32. LaiERileyJPurvisIRosesAA 4-Mb high-density single nucleotide polymorphism-based map around human APOEGenomics1998543138[PubMed][Google Scholar]
33. PritchardJKStephensMRosenbergNADonnellyPAssociation mapping in structured populationsAm. J. Hum. Genet200067170181[PubMed][Google Scholar]
34. PritchardJKStephensMDonnellyPJInference of population structure using multilocus genotype dataGenetics2000155945959[PubMed][Google Scholar]
35. PriceALPrincipal components analysis corrects for stratification in genome-wide association studiesNat. Genet200638904909[PubMed][Google Scholar]
36. CareyVJGEE: Generalized Estimation Equation Solver. 4.13–15edn2010[Google Scholar]
37. ChenMHYangQGWAF: an R package for genome-wide association analyses with family dataBioinformatics201026580581[PubMed][Google Scholar]
38. R Development Core Team. R: A language and environment for statistical computingR Foundation for Statistical Computing, V., Austria2009ISBN 3-900051-07-0, URL http://www.R-project.org
39. WillerCJLiYAbecasisGRMETAL: fast and efficient meta-analysis of genomewide association scansBioinformatics20102621902191[PubMed][Google Scholar]
40. HigginsJPThompsonSGQuantifying heterogeneity in a meta-analysisStat Med20022115391558[PubMed][Google Scholar]
41. HigginsJPThompsonSGDeeksJJAltmanDGMeasuring inconsistency in meta-analysesBMJ2003327557560[PubMed][Google Scholar]
42. AbecasisGRWiggintonJEHandling marker-marker linkage disequilibrium: Pedigree analysis with clustered markersAm. J. Hum. Genet200577754767[PubMed][Google Scholar]