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

Previous genome-wide association studies (GWAS), conducted by our group and others, have identified loci that harbor risk variants for neurodegenerative diseases, including Alzheimer's disease (AD). Human disease variants are enriched for polymorphisms that affect gene expression, including some that are known to associate with expression changes in the brain. Postulating that many variants confer risk to neurodegenerative disease via transcriptional regulatory mechanisms, we have analyzed gene expression levels in the brain tissue of subjects with AD and related diseases. Herein, we describe our collective datasets comprised of GWAS data from 2,099 subjects; microarray gene expression data from 773 brain samples, 186 of which also have RNAseq; and an independent cohort of 556 brain samples with RNAseq. We expect that these datasets, which are available to all qualified researchers, will enable investigators to explore and identify transcriptional mechanisms contributing to neurodegenerative diseases.

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...

Mutations in the progranulin gene (GRN) are an important cause of frontotemporal lobar degeneration (FTLD) with ubiquitin and TAR DNA-binding protein 43 (TDP43)-positive pathology. The clinical presentation associated with GRN mutations is heterogeneous and may include clinical probable Alzheimer's disease. All GRN mutations identified thus far cause disease through a uniform disease mechanism, i.e. the loss of functional GRN or haploinsufficiency. To determine if expression of GRN in plasma could predict GRN mutation status and could be used as a biological marker, we optimized a GRN ELISA and studied plasma samples of a consecutive clinical FTLD series of 219 patients, 70 control individuals, 72 early-onset probable Alzheimer's disease patients and nine symptomatic and 18 asymptomatic relatives of GRN mutation families. All FTLD patients with GRN loss-of-function mutations showed significantly reduced levels of GRN in plasma to about one third of the levels observed in non-GRN carriers and control individuals (P < 0.001). No overlap in distributions of GRN levels was observed between the eight GRN loss-of-function mutation carriers (range: 53–94 ng/ml) and 191 non-GRN mutation carriers (range: 115–386 ng/ml). Similar low levels of GRN were identified in asymptomatic GRN mutation carriers. Importantly, ELISA analyses also identified one probable Alzheimer's disease patient (1.4%) carrying a loss-of-function mutation in GRN. Biochemical analyses further showed that the GRN ELISA only detects full-length GRN, no intermediate granulin fragments. This study demonstrates that using a GRN ELISA in plasma, pathogenic GRN mutations can be accurately detected in symptomatic and asymptomatic carriers. The ∼75% reduction in full-length GRN, suggests an unbalanced GRN metabolism in loss-of-function mutation carriers whereby more GRN is processed into granulins. We propose that plasma GRN levels could be used as a reliable and inexpensive tool to identify all GRN mutation carriers in early-onset dementia populations and asymptomatic at-risk individuals.

Genetic variants that modify brain gene expression may also influence risk for human diseases. We measured expression levels of 24,526 transcripts in brain samples from the cerebellum and temporal cortex of autopsied subjects with Alzheimer's disease (AD, cerebellar n = 197, temporal cortex n = 202) and with other brain pathologies (non–AD, cerebellar n = 177, temporal cortex n = 197). We conducted an expression genome-wide association study (eGWAS) using 213,528 cisSNPs within ±100 kb of the tested transcripts. We identified 2,980 cerebellar cisSNP/transcript level associations (2,596 unique cisSNPs) significant in both ADs and non–ADs (q<0.05, p = 7.70×10−5–1.67×10−82). Of these, 2,089 were also significant in the temporal cortex (p = 1.85×10−5–1.70×10−141). The top cerebellar cisSNPs had 2.4-fold enrichment for human disease-associated variants (p<10−6). We identified novel cisSNP/transcript associations for human disease-associated variants, including progressive supranuclear palsy SLCO1A2/rs11568563, Parkinson's disease (PD) MMRN1/rs6532197, Paget's disease OPTN/rs1561570; and we confirmed others, including PD MAPT/rs242557, systemic lupus erythematosus and ulcerative colitis IRF5/rs4728142, and type 1 diabetes mellitus RPS26/rs1701704. In our eGWAS, there was 2.9–3.3 fold enrichment (p<10−6) of significant cisSNPs with suggestive AD–risk association (p<10−3) in the Alzheimer's Disease Genetics Consortium GWAS. These results demonstrate the significant contributions of genetic factors to human brain gene expression, which are reliably detected across different brain regions and pathologies. The significant enrichment of brain cisSNPs among disease-associated variants advocates gene expression changes as a mechanism for many central nervous system (CNS) and non–CNS diseases. Combined assessment of expression and disease GWAS may provide complementary information in discovery of human disease variants with functional implications. Our findings have implications for the design and interpretation of eGWAS in general and the use of brain expression quantitative trait loci in the study of human disease genetics.