Diego. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California. We thank Drs. D. Stephen Snyder and Marilyn Miller from NIA who are ex-officio ADGC members. EADI. This work has been developed and supported by the LABEX (laboratory of excellence program investment for the future) DISTALZ grant (Development of Innovative Strategies for a Transdisciplinary approach to ALZheimer's disease) including funding from MEL (Metropole européenne de Lille), ERDF (European Regional Development Fund) and Conseil Régional Rotterdam, Netherlands Organization for the Health Research and Development (ZonMw), the Research Institute for Diseases in the Elderly (RIDE), the Ministry of Education, Culture and Science, the Ministry for Health, Welfare and Sports, the European Commission (DG XII), and the Municipality of Rotterdam. The authors are grateful to the study participants, the staff from the Rotterdam Study and the participating general practitioners and pharmacists. The generation and management of GWAS genotype data for the Rotterdam Study (RS-I, RS-II, RS-III) was executed by the Human Genotyping Facility of the Genetic Laboratory of the

Introduction We identified rare coding variants associated with Alzheimer’s disease (AD) in a 3-stage case-control study of 85,133 subjects. In stage 1, 34,174 samples were genotyped using a whole-exome microarray. In stage 2, we tested associated variants (P<1×10-4) in 35,962 independent samples using de novo genotyping and imputed genotypes. In stage 3, an additional 14,997 samples were used to test the most significant stage 2 associations (P<5×10-8) using imputed genotypes. We observed 3 novel genome-wide significant (GWS) AD associated non-synonymous variants; a protective variant in PLCG2 (rs72824905/p.P522R, P=5.38×10-10, OR=0.68, MAFcases=0.0059, MAFcontrols=0.0093), a risk variant in ABI3 (rs616338/p.S209F, P=4.56×10-10, OR=1.43, MAFcases=0.011, MAFcontrols=0.008), and a novel GWS variant in TREM2 (rs143332484/p.R62H, P=1.55×10-14, OR=1.67, MAFcases=0.0143, MAFcontrols=0.0089), a known AD susceptibility gene. These protein-coding changes are in genes highly expressed in microglia and highlight an immune-related protein-protein interaction network enriched for previously identified AD risk genes. These genetic findings provide additional evidence that the microglia-mediated innate immune response contributes directly to AD development.

IntroductionLate-onset Alzheimer's disease (LOAD, onset age > 60 years) is the most prevalent dementia in the elderly 1 , and risk is partially driven by genetics 2 . Many of the loci responsible for this genetic risk were identified by genome-wide association studies (GWAS) [3][4][5][6][7][8] . To identify additional LOAD risk loci, the we performed the largest GWAS to date (89,769 individuals), analyzing both common and rare variants. We confirm 20 previous LOAD risk loci and identify four new genome-wide loci (IQCK, ACE, ADAM10, and ADAMTS1). Pathway analysis of these data implicates the immune system and lipid metabolism, and for the first time tau binding proteins and APP metabolism. These findings show that genetic variants affecting APP and Aβ processing are not only associated with early-onset autosomal dominant AD but also with LOAD. Analysis of AD risk genes and pathways show enrichment for rare variants (P = 1.32 x 10 -7 ) indicating that additional rare variants remain to be identified. Main TextOur previous work identified 19 genome-wide significant common variant signals in addition to APOE 9 , that influence risk for LOAD. These signals, combined with 'subthreshold' common variant associations, account for ~31% of the genetic variance of LOAD 2 , leaving the majority of genetic risk uncharacterized 10 . To search for additional signals, we conducted a GWAS metaanalysis of non-Hispanic Whites (NHW) using a larger sample (17 new, 46 total datasets) from our group, the International Genomics of Alzheimer's Project (IGAP) (composed of four AD consortia: ADGC, CHARGE, EADI, and GERAD). This sample increases our previous discovery sample (Stage 1) by 29% for cases and 13% for controls (N=21,982 cases; 41,944 controls) ( Supplementary Table 1 and 2, and Supplementary Note). To sample both common and rare variants (minor allele frequency MAF ≥ 0.01, and MAF < 0.01, respectively), we imputed the discovery datasets using a 1000 Genomes reference panel consisting of . CC-BY-NC-ND 4.0 International license peer-reviewed) is the author/funder. It is made available under a 11 36,648,992 single-nucleotide variants, 1,380,736 insertions/deletions, and 13,805 structural variants. After quality control, 9,456,058 common variants and 2,024,574 rare variants were selected for analysis (a 63% increase from our previous common variant analysis in 2013).Genotype dosages were analyzed within each dataset, and then combined with meta-analysis ( Supplementary Figures 1 and 2 and Supplementary Table 3). The Stage 1 discovery metaanalysis was first followed by Stage 2 using the I-select chip we previously developed in Lambert et al (including 11,632 variants, N=18,845) and finally stage 3A (N=6,998). The final sample was 33,692 clinical AD cases and 56,077 controls.Meta-analysis of Stages 1 and 2 produced 21 associations with P ≤ 5x10 -8 (Table 1 and Figure 1). Of these, 18 were previously reported as genome-wide significant and three of them are signals not initially described in Lambert et al: the rare R47H TREM2 coding va...

APOE ε4, the most significant genetic risk factor for Alzheimer disease (AD), may mask effects of other loci. We re-analyzed genome-wide association study (GWAS) data from the International Genomics of Alzheimer’s Project (IGAP) Consortium in APOE ε4+ (10,352 cases and 9,207 controls) and APOE ε4− (7,184 cases and 26,968 controls) subgroups as well as in the total sample testing for interaction between a SNP and APOE ε4 status. Suggestive associations (P<1x10−4) in stage 1 were evaluated in an independent sample (stage 2) containing 4,203 subjects (APOE ε4+: 1,250 cases and 536 controls; APOE ε4-: 718 cases and 1,699 controls). Among APOE ε4− subjects, novel genome-wide significant (GWS) association was observed with 17 SNPs (all between KANSL1 and LRRC37A on chromosome 17 near MAPT) in a meta-analysis of the stage 1 and stage 2 datasets (best SNP, rs2732703, P=5·8x10−9). Conditional analysis revealed that rs2732703 accounted for association signals in the entire 100 kilobase region that includes MAPT. Except for previously identified AD loci showing stronger association in APOE ε4+ subjects (CR1 and CLU) or APOE ε4− subjects (MS4A6A/MS4A4A/ MS4A6E), no other SNPs were significantly associated with AD in a specific APOE genotype subgroup. In addition, the finding in the stage 1 sample that AD risk is significantly influenced by the interaction of APOE with rs1595014 in TMEM106B (P=1·6x10−7) is noteworthy because TMEM106B variants have previously been associated with risk of frontotemporal dementia. Expression quantitative trait locus analysis revealed that rs113986870, one of the GWS SNPs near rs2732703, is significantly associated with four KANSL1 probes that target transcription of the first translated exon and an untranslated exon in hippocampus (P≤1.3x10−8), frontal cortex (P≤1.3x10−9), and temporal cortex (P≤1.2x10−11). Rs113986870 is also strongly associated with a MAPT probe that targets transcription of alternatively spliced exon 3 in frontal cortex (P=9.2x10−6) and temporal cortex (P=2.6x10−6). Our APOE-stratified GWAS is the first to show GWS association for AD with SNPs in the chromosome 17q21.31 region. Replication of this finding in independent samples is needed to verify that SNPs in this region have significantly stronger effects on AD risk in persons lacking APOE ε4 compared to persons carrying this allele, and if this is found to hold, further examination of this region and studies aimed at deciphering the mechanism(s) are warranted.

We report the development and validation of experimental methods, study designs, and analysis software for pooling-based genomewide association (GWA) studies that use high-throughput single-nucleotide-polymorphism (SNP) genotyping microarrays. We first describe a theoretical framework for establishing the effectiveness of pooling genomic DNA as a low-cost alternative to individually genotyping thousands of samples on high-density SNP microarrays. Next, we describe software called "GenePool," which directly analyzes SNP microarray probe intensity data and ranks SNPs by increased likelihood of being genetically associated with a trait or disorder. Finally, we apply these methods to experimental case-control data and demonstrate successful identification of published genetic susceptibility loci for a rare monogenic disease (sudden infant death with dysgenesis of the testes syndrome), a rare complex disease (progressive supranuclear palsy), and a common complex disease (Alzheimer disease) across multiple SNP genotyping platforms. On the basis of these theoretical calculations and their experimental validation, our results suggest that pooling-based GWA studies are a logical first step for determining whether major genetic associations exist in diseases with high heritability.