- Open Access
Genome-wide pleiotropy analysis of neuropathological traits related to Alzheimer’s disease
Alzheimer's Research & Therapyvolume 10, Article number: 22 (2018)
Simultaneous consideration of two neuropathological traits related to Alzheimer’s disease (AD) has not been attempted in a genome-wide association study.
We conducted genome-wide pleiotropy analyses using association summary statistics from the Beecham et al. study (PLoS Genet 10:e1004606, 2014) for AD-related neuropathological traits, including neuritic plaque (NP), neurofibrillary tangle (NFT), and cerebral amyloid angiopathy (CAA). Significant findings were further examined by expression quantitative trait locus and differentially expressed gene analyses in AD vs. control brains using gene expression data.
Genome-wide significant pleiotropic associations were observed for the joint model of NP and NFT (NP + NFT) with the single-nucleotide polymorphism (SNP) rs34487851 upstream of C2orf40 (alias ECRG4, P = 2.4 × 10−8) and for the joint model of NFT and CAA (NFT + CAA) with the HDAC9 SNP rs79524815 (P = 1.1 × 10−8). Gene-based testing revealed study-wide significant associations (P ≤ 2.0 × 10−6) for the NFT + CAA outcome with adjacent genes TRAPPC12, TRAPPC12-AS1, and ADI1. Risk alleles of proxy SNPs for rs79524815 were associated with significantly lower expression of HDAC9 in the brain (P = 3.0 × 10−3), and HDAC9 was significantly downregulated in subjects with AD compared with control subjects in the prefrontal (P = 7.9 × 10−3) and visual (P = 5.6 × 10−4) cortices.
Our findings suggest that pleiotropy analysis is a useful approach to identifying novel genetic associations with complex diseases and their endophenotypes. Functional studies are needed to determine whether ECRG4 or HDAC9 is plausible as a therapeutic target.
Alzheimer’s disease (AD) is the most common type of dementia in persons aged 65 years and older [1, 2]. Pathologically, it is characterized primarily by the appearance of both neuritic plaques (NPs) containing oligomers of β-amyloid and neurofibrillary tangles (NFTs), accompanied by a progressive loss of neurons in the brain [3, 4]. Also, cerebral amyloid angiopathy (CAA), which is caused by aggregates of β-amyloid in walls of blood vessels in the brain, is found in as many as 90% of autopsy-confirmed AD cases . Previously, Beecham et al. identified multiple significant gene associations in a genome-wide association study (GWAS) for several AD-related neuropathological traits, including NP, NFT, and CAA measured in brains from subjects with pathologically confirmed AD cases and from control subjects with no evidence of neurological disease . We hypothesized that additional novel associations could be identified in models allowing a genetic variant to influence more than one trait (i.e., pleiotropy). In this study, we performed genome-wide pleiotropy analyses of joint models of NP, NFT, and CAA using summary data from the previous study .
Study population, neuropathological trait selection, and data processing
We obtained summary statistics from univariate GWAS of NP, NFT, and CAA . These results were derived from meta-analyses of 12 studies including 3598 subjects (3135 AD cases, 463 controls) of European ancestry. Neuropathological data for the entire sample were reviewed and harmonized by one neuropathologist for consistency across studies . Although Beecham et al.  also evaluated Lewy body disease, hippocampal sclerosis, and vascular brain disease, we limited our present analyses to neuropathological outcomes most directly linked to AD and moderately correlated with each other (i.e., NP, NFT, and CAA). Uncorrelated traits are unlikely to show significant pleiotropic associations, and results from pleiotropy analysis will be similar to those from univariate models (i.e., single phenotype) if the traits are highly correlated. Details of subject recruitment, genotyping, genotype imputation, quality control procedures, population substructure analysis, and statistical methods for association analyses of individual traits were reported previously [6, 7]. Sample demography of the 3598 subjects with autopsied brains and genotypes (3135 cases and 463 controls) is described in Additional file 1: Table S1.
Univariate genome-wide association analyses
Results from the association tests by Beecham et al.  in each dataset for each neuropathologic trait with genotypes imputed using the 1000 Genomes Project reference panel (GRCh37 at December 2010) for a genome-wide set of single-nucleotide polymorphisms (SNPs) were obtained using ordinal logistic regression models including the first three principal components of ancestry as covariates to account for population substructure . NP and NFT measures were analyzed in well-established ordinal rankings (NPs: none, sparse, moderate, and frequent by Consortium to Establish a Registry for Alzheimer’s Disease [“CERAD”] scoring ; NFT: none, transentorhinal, limbic, and isocortical by Braak and Braak staging ), and CAA was analyzed as a binary trait (present or absent). Full details of these analyses are reported elsewhere . We used GWAS meta-analysis summary statistics (β and SE) of the three neuropathologic traits for 6.5 million imputed SNPs after omitting SNPs from studies if the minor allele frequency was ≤ 1%, imputation quality (R2) was ≤ 0.4, or dosage variance was ≤ 0.02.
Genome-wide pleiotropy analyses
We conducted a genome-wide pleiotropy analysis for each pair of the three neuropathological traits using the O’Brien method [10, 11], which is implemented in an R library (“CUMP”) . This method combines univariate test statistics (Z-scores from β and SE values) of all SNPs from separate GWASs for individual phenotypes to compute a test statistic that follows a multivariate normal distribution. The covariance matrix of the distribution was approximated by the sample covariance matrix of the test statistics of all SNPs. Under the null hypothesis, an SNP is not associated with any of the phenotypes. The alternative hypothesis is that an SNP is associated with at least one of the phenotypes. We defined a SNP as having a pleiotropic effect on two phenotypes when the P value for the O’Brien test statistic from the joint model of association of two phenotypes (Pjoint) with the SNP is at least one order of magnitude more significant than the P values (Punivariate) for both phenotypes and the univariate P values are at least nominally significant (Punivariate < 0.05). As a supplementary analysis, we also conducted a trivariate pleiotropy genome-wide analysis for the three neuropathological traits. The genome-wide significance (GWS) threshold for these analyses was set at P < 5.0 × 10−8.
We performed genome-wide gene-based tests for each joint model using results from individual SNP tests. SNPs within 30 kb of the transcription start and end sites were included in each gene-based test. These analyses were carried out using the versatile gene-based test (“VEGAS”) method , which computes an empirical P value through Monte Carlo simulations based on linkage disequilibrium patterns of the European ancestry population in the 1000 Genomes Project (GRCh37 released March 2012). The GWS level for the gene-based tests was set at 2.7 × 10−6, which was calculated as the nominal significance level 0.05 divided by the total number of genes tested (n = 18,500).
Expression quantitative trait locus analysis
The association of SNP genotypes with gene-level expression (i.e., expression quantitative trait loci [eQTLs]) was evaluated using version 6 of the GTEx Portal database (http://www.gtexportal.org/; ) and data from the Mayo Clinic brain expression GWAS (eGWAS) (https://www.synapse.org/#!Synapse:syn3157249 or http://alois.med.upenn.edu/niagads; ). The GTEx Portal provides eQTL association summary statistics (β and P values) across 43 different tissues from 175 subjects. The Mayo Clinic brain eGWAS data were generated from the cerebellum (CER; 197 AD and 177 non-AD control subjects) and temporal cortex (TCX; 202 AD and 197 non-AD control subjects) regions. Gene expression measures for 24,526 probes were generated with the Illumina Whole Genome DASL array (Illumina, San Diego, CA, USA). SNP genotype data for the Mayo Clinic eGWAS were obtained from the Mayo Clinic late-onset AD GWAS . AD cases were diagnosed as definite AD according to National Institute of Neurological and Communicative Disorders and Stroke/Alzheimer’s Disease and Related Disorders Association criteria, whereas non-AD controls had other neuropathologies. For each brain region, association of gene expression and imputed SNP genotype (GRCh36) was evaluated using linear regression, including covariates for AD status, apolipoprotein E (APOE) ε4 dosage (0, 1, or 2), age at death, sex, plate, RNA integrity number (RIN), and adjusted RIN (RIN − RINmean2). Analyses were also conducted for AD cases and controls separately.
Differential gene expression analysis
Differential gene expression (DGE) analysis was performed using publicly available brain whole-transcriptome RNA-sequencing (RNA-Seq) data  and microarray data (Gene Expression Omnibus accession number [GEO:GSE44772] ). The RNA-Seq data include DGE summary statistics for the CER and TCX derived from 86 patients with AD and 80 control subjects (https://www.synapse.org). Following a quality control step, 80 AD and 76 control brains were analyzed. All subjects underwent RNA-Seq using the Illumina HiSeq 2000 sequencing system (101 bp, paired-end sequencing) at the Mayo Clinic Genomic Core Facility. All AD and some of the control brains were from the Mayo Clinic Brain Bank, whereas other control brains were from the Banner Sun Health Research Institute (Sun City, AZ, USA). Following quality control, raw read counts normalized according to conditional quantile normalization (CQN) employing the Bioconductor package were used in the analyses. For DGE comparing AD with controls, multivariable linear regression analyses were conducted in R, using CQN normalized gene expression measures and including age at death, sex, RIN, brain tissue source, and flow cell as biological and technical covariates. To account for any CNS cell-population changes that occur as a consequence of disease pathology, we also included cell-specific gene levels as covariates, using the expression levels for the five central nervous system (CNS)-specific genes as follows: ENO2 for neurons [ENCODE:ENSG00000111674], GFAP for astrocytes [ENCODE:ENSG00000131095], CD68 for microglia [ENCODE:ENSG00000129226], OLIG2 for oligodendrocytes [ENCODE:ENSG00000205927], and CD34 for endothelial cells [ENCODE:ENSG00000174059].
The microarray gene expression data were generated from autopsied brains collected from dorsolateral prefrontal cortex (DLPFC), visual cortex (VCX), and CER regions of 129 AD patients and 101 control subjects. Samples were profiled on a custom-made Agilent 44K array (Agilent Technologies, Santa Clara, CA, USA) containing 40,638 human genes. Gene expression data were normalized using Rosetta Resolver gene expression analysis software as previously described . The association between expression of each gene (outcome) and AD status (predictor) was tested using linear regression adjusting for RIN, postmortem interval, batch, preservation method, pH in tissues, age, sex, and the five cell-type markers.
NP, NFT, and CAA were moderately correlated (NP-NFT, r = 0.68; NP-CAA, r = 0.56; NFT-CAA, r = 0.40; P < 2.2 × 10−16 for each pair of traits), indicating a potential for discovery of novel associations in pleiotropy analysis.
Bivariate GWAS results
There was no inflation in P values for the GWAS of the three neuropathological traits analyzed individually (genomic control parameter, λ = 1.00, 1.01, and 0.96 for NP, NFT, and CAA, respectively) or as joint outcomes (Additional file 1: Figure S1). Results of the pleiotropy GWAS are shown in Additional file 1: Figure S2. As reported previously, with the exception of APOE, only 15 of 25 previously known AD loci attained at least a nominal association with NP, NFT, or CAA . Three of the previously established AD loci—BIN1, HLA region, and PICALM—were moderately associated (P < 10−4) in the pleiotropy analysis for NP and NFT at a significance level of at least one order of magnitude smaller compared with the results from univariate analyses (Additional file 1: Table S2). Two novel GWS associations were detected in the pleiotropy analyses (Table 1, Fig. 1). rs34487851, an SNP located approximately 40 kb upstream of C2orf40, was associated with the joint model of NP and NFT (Pjoint = 2.0 × 10−8). An intronic SNP in HDAC9, rs79524815, was associated with the joint model of NFT and CAA (Pjoint = 1.1 × 10−8). The major allele A of rs34487851 and the minor allele G of rs79524815 are associated with increased NP and NFT and with increased NFT and the presence of CAA, respectively. Both of these findings were at least one order of magnitude more significant than for the univariate traits (Table 1) and were supported by evidence from multiple SNPs at those locations (Fig. 1, Additional file 1: Table S3).
Bivariate gene-based pleiotropy analysis results
Three contiguous novel genes on chromosome 2p25.3 (TRAPPC12, TRAPPC12-AS1, and ADI1) were associated with the joint model of NFT and CAA at a gene-wide significant level (P ≤ 2.0 × 10−6) (Table 2 and Additional file 1: Figure S3). Of note, one SNP in this region (rs35067331 in TRAPPC12) was associated with the NFT-CAA outcome at nearly the GWS level (Pjoint = 5.8 × 10−8) (Additional file 1: Table S4).
Trivariate GWAS and gene-based pleiotropy analysis results
We conducted trivariate GWAS and gene-based association analyses to identify genetic factors common to NP, NFT, and CAA. There was no evidence for genomic inflation (λ = 0.99) in the results from the trivariate model (Additional file 1: Figure S4). GWS association was observed only for APOE isoform SNP rs429358 (P = 2.1 × 10−47), whereas associations at C2orf40, HDAC9, and TRAPPC12 were attenuated (Additional file 1: Table S5).
We performed eQTL association analysis to examine whether the expression levels of the five GWS significant genes identified in the pleiotropy association tests differed between carriers and noncarriers of the risk alleles from those loci. Because information about the two GWS SNPs was not available in the GTEx Portal database or in the Mayo Clinic brain eGWAS, we analyzed proxy SNPs that are in high linkage disequilibrium (LD; D′ ≥ 80) with the GWS SNPs. According to GTEx, rs34487851 proxy SNP rs1232803 is a cis-acting eQTL, and the major allele A, which is associated with higher NP and NFT, is also significantly associated with decreased expression of C2orf40 in several tissues, including the esophagus (P = 3.5 × 10−5), transverse colon (P = 4.7 × 10−4), and tibial artery (P = 1.7 × 10−3), but not in any brain regions. In the Mayo Clinic brain eGWAS, proxy SNPs for rs34487851 were not cis-acting eQTLs for C2orf40. In GTEx, proxy SNPs for rs79524815 were not associated with the expression of HDAC9. However, in the brain eGWAS, the minor alleles of proxy SNPs for rs79524815, which are associated with higher NFT and CAA, were significantly associated with lower HDAC9 levels in the CER (probe ID: ILMN_1803563; best eQTL, rs4721721; P = 0.003) but not in the TCX (Additional file 1: Table S3). According to GTEx, rs35067331 is a cis-acting eQTL, and its major allele C, which is associated with higher NFT and CAA, is significantly associated with increased expression of TRAPPC12-AS1 in several brain regions (best P = 2.1 × 10−7 in cortex) and ADI1 (P = 0.03) in the caudate nucleus, but not with differential expression of TRAPPC12 in any brain regions (Additional file 1: Table S4). In the Mayo Clinic brain eGWAS data, rs35067331 and its proxy SNPs were not cis-acting eQTLs for ADI1 or TRAPPC12. Unfortunately, information about TRAPPC12-AS1 was unavailable in the brain eGWAS.
Differential gene expression analysis
We investigated whether the expression levels of C2orf40, HDAC9, and TRAPPC12/TRAPPC12-AS1/ADI1 differed in AD brains compared with non-AD control brains in the publicly available RNS-Seq and microarray datasets (Table 3 and Fig. 2). There were no significant differences in C2orf40 expression between subjects with AD and control subjects in the TCX or CER in the Mayo Clinic RNA-Seq DGE profiling. However, C2orf40 was significantly downregulated in subjects with AD compared with control subjects in the CER (P = 1.6 × 10−3), DLPFC (P = 0.04), and VCX (P = 2.7 × 10−3) in the microarray brain expression data. HDAC9 was significantly downregulated in subjects with AD compared with control subjects in several brain regions, including the TCX (P = 1.5 × 10−4) and CER (P = 0.04) in the RNA-Seq profiling data and in the DLPFC (P = 7.9 × 10−3) and VCX (P = 5.6 × 10−4) in the microarray expression data. ADI1 expression was downregulated in subjects with AD in the CER in the microarray data (P = 4.9 × 10−4). The RNA-Seq DGE profiling indicated that TRAPPC12-AS1 expression was significantly increased in subjects with AD in the TCX (P = 1.3 × 10−3). In contrast, expression of TRAPPC12 was significantly lower in subjects with AD than in control subjects in the TCX (P = 0.01) in the RNA-Seq and in the CER (P = 1.1 × 10−3) and VCX (P = 3.2 × 10−4) in the microarray data. Information about expression of TRAPPC12-AS1 was not available in the microarray data. The majority of these nominally significant findings survived correction for multiple testing (P < 2.27 × 10−3).
To contextualize our findings, we evaluated DGE among AD cases and controls for 26 previously established AD genes (Additional file 1: Table S6) [19,20,21]. With notable exceptions of EPHA1 (P = 2.8 × 10−7) and SLC24A4 (P = 7.0 × 10−5) in CER in the RNA-Seq data and ABCA7 in DLPFC (P = 2.5 × 10−7) and VCX (P = 4.9 × 10−7) in the microarray data, none of the results for the other 23 genes were significant after correcting for 127 tests (P < 3.94 × 10−4).
A previous GWAS of neuropathologic traits including NP, NFT, and CAA identified GWS associations with APOE only . Our pleiotropy analysis of all pairwise combinations of these traits identified GWS associations with APOE and three regions not previously reported with any neuropathologic traits or AD risk, including C2orf40 for the joint model comprising NP and NFT, as well as HDAC9 and TRAPPC12/TRAPPC12-AS1/ADI1 for the joint model comprising NFT and CAA. Our DGE study found that HDAC9 is significantly downregulated in several brain regions in subjects with AD compared with control subjects. Moreover, we observed that the G allele of HDAC9 SNP rs79524815 is associated with a higher level of the joint outcome of NFT and CAA, and proxy SNPs for rs79524815 (which are suggestively associated with the joint outcome of NFT and CAA) are associated with decreased HDAC9 expression in subjects with AD (Additional file 1: Table S3). The pleiotropy analysis also revealed that 4 (BIN1, HLA, PICALM, and APOE) of the 25 previously reported GWS AD risk loci [7, 19,20,21] were at least one order of magnitude more significantly associated with the joint model of NP + NFT than each of these traits analyzed separately, suggesting that these genes are involved in pathways leading to both plaques and tangles [22,23,24,25].
It is notable that pleiotropy analyses for the model including NP, NFT, and CAA did not yield any GWS findings. Moreover, the GWS associations identified in the bivariate models were attenuated in the trivariate model. These results suggest that the mechanisms or pathways underlying the bivariate associations probably do not encompass all three traits, and this conclusion may generalize genome-wide.
C2orf40, also known as esophageal cancer-related gene 4 (ECRG4), is a tumor suppressor gene  that encodes a peptide hormone that is involved in NFT formation in transgenic mice , senescence of precursor cells in the CNS during aging , and activation of microglia and peripheral mononuclear leukocytes . We observed that rs34487851 allele A is associated with higher NP and NFT and lower expression of ERCG4, albeit not in the brain. We also found that ERCG4 expression was significantly lower in AD cases than in controls in several brain regions. Abnormal downregulation of C2orf40 was previously reported in brain injury  as well as in several cancer cell types [31, 32].
HDAC9 encodes a member of class IIa histone genes that deacetylate histones, thereby remodeling chromatin structure and controlling gene expression [33, 34] that has previously been linked to epigenetic mechanisms  and memory loss  in AD and also has been proposed as a possible therapeutic target [37,38,39]. GWS association of ischemic stroke with an HDAC9 variant was identified by GWAS . Structural variants including deletions and copy number variants in HDAC9 have been identified in patients with schizophrenia and patients with autism [41, 42]. MEF2C, one of the well-established AD risk loci , stimulates HDAC9 expression, but HDAC9 suppresses MEF2C transcription, resulting in a negative feedback loop . In a previously reported coexpression network study in AD and control brains, HDAC9 and MEF2C were clustered together in the top fourth module ranked by relevance to AD pathology, and expression of HDAC9 and MEF2C was inversely correlated with Braak stage (HDAC9, r = −0.71; MEF2C, r = −0.65) and frontal atrophy (HDAC9, r = −0.57, MEF2C, r = −0.51) . These findings are consistent with our observation that HDAC9 expression is reduced in subjects with AD and in the subjects with HDAC9 SNP alleles associated with higher NFT and CAA. Decreased HDAC9 expression has also been linked to increased neuronal apoptosis [44, 45]. Collectively, findings from our and other studies indicate that MEF2C and HDAC9 may participate in a pathway leading to NFT formation and brain atrophy.
Gene-based analyses identified significant associations with three adjacent loci—TRAPPC12, TRAPPC12-AS1, and ADI1—in a gene-rich region near the end of the short arm of chromosome 2p. ADI1, encoding acireductone dioxygenase 1, is involved in methionine salvage and prostate cancer  and has no known relationship to AD. TRAPPC12 is a subunit of a trafficking protein particle complex that has a role in vesicle trafficking in endoplasmic reticulum (ER) to Golgi . TRAPPC12-AS1 is an antisense (noncoding) RNA that contains a 1168 transcript from TRAPPC12. We previously established that regulation of vesicular trafficking in the ER to Golgi by several VPS10 receptor domain receptor genes, including SORL1, and by other genes encoding members of the retromer complex is an important pathway leading to AD [48,49,50]. Of the genes in this region, only TRAPPC12-AS1 showed a pattern of expression in subjects with AD and control subjects that is consistent with the effect direction of the TRAPPC12 rs35067331 allele’s influence on NFT and CAA. It should be noted that TRAPPC12 expression was significantly lower in AD cases than in control subjects in the TCX (P = 0.01) in the RNA-Seq data and the CER (P = 1.1 × 10−3) and VCX (P = 3.2 × 10−4) in the microarray data, which could be due to negative feedback by the antisense TRAPPC12-AS1 transcript .
Our study has several potential caveats. The GWS associations identified in the pleiotropy analysis have moderate supportive evidence for association from other SNPs under the association peaks, probably because of low LD with the top SNPs. The two GWS SNPs near C2orf40 and HDAC9 were not associated with AD risk in one of the largest GWASs for AD (rs34487851, P = 0.07; rs79524815, P = 0.73) [19, 21]. However, approximately 87% of the autopsy samples used in this pleiotropy analysis (as well as in the Beecham et al. study ) were from patients with AD. This may indicate that our findings are more relevant with neuropathological progression after onset of AD clinical symptoms. Alternatively, because our study was focused on endophenotypes that might be more proximal than disease diagnosis to effects of the genetic variants [52, 53], our analyses might have more power to detect those novel associations. Finally, to our knowledge, additional large late-onset AD cohorts with neuropathological and genotype data are not currently available for replication of our association findings. Therefore, validation of the role of these loci in AD will likely require experimental evidence.
Our findings suggest that genome-wide pleiotropy analysis is a useful approach to identifying novel genetic associations with complex diseases and their endophenotypes. Functional studies are needed to determine whether C2orf40 or HDAC9 is a plausible therapeutic target.
Jorm AF, Jolley D. The incidence of dementia: a meta-analysis. Neurology. 1998;51(3):728–33.
Corrada MM, Brookmeyer R, Paganini-Hill A, Berlau D, Kawas CH. Dementia incidence continues to increase with age in the oldest old: the 90+ study. Ann Neurol. 2010;67(1):114–21.
Shoji M, Golde TE, Ghiso J, Cheung TT, Estus S, Shaffer LM, Cai XD, McKay DM, Tintner R, Frangione B, et al. Production of the Alzheimer amyloid β protein by normal proteolytic processing. Science. 1992;258(5079):126–9.
Trojanowski JQ, Schmidt ML, Shin RW, Bramblett GT, Rao D, Lee VM. Altered tau and neurofilament proteins in neuro-degenerative diseases: diagnostic implications for Alzheimer’s disease and Lewy body dementias. Brain Pathol. 1993;3(1):45–54.
Love S. Contribution of cerebral amyloid angiopathy to Alzheimer’s disease. J Neurol Neurosurg Psychiatry. 2004;75(1):1–4.
Beecham GW, Hamilton K, Naj AC, Martin ER, Huentelman M, Myers AJ, Corneveaux JJ, Hardy J, Vonsattel JP, Younkin SG, et al. Genome-wide association meta-analysis of neuropathologic features of Alzheimer’s disease and related dementias. PLoS Genet. 2014;10(9):e1004606.
Naj AC, Jun G, Beecham GW, Wang LS, Vardarajan BN, Buros J, Gallins PJ, Buxbaum JD, Jarvik GP, Crane PK, et al. Common variants at MS4A4/MS4A6E, CD2AP, CD33 and EPHA1 are associated with late-onset Alzheimer’s disease. Nat Genet. 2011;43(5):436–41.
Mirra SS, Hart MN, Terry RD. Making the diagnosis of Alzheimer’s disease: a primer for practicing pathologists. Arch Pathol Lab Med. 1993;117(2):132–44.
Braak H, Braak E. Neuropathological stageing of Alzheimer-related changes. Acta Neuropathol. 1991;82(4):239–59.
O’Brien PC. Procedures for comparing samples with multiple endpoints. Biometrics. 1984;40(4):1079–87.
Yang Q, Wu H, Guo CY, Fox CS. Analyze multivariate phenotypes in genetic association studies by combining univariate association tests. Genet Epidemiol. 2010;34(5):444–54.
Yang Q, Wang Y. Methods for analyzing multivariate phenotypes in genetic association studies. J Probab Stat. 2012;2012:652569.
Liu JZ, McRae AF, Nyholt DR, Medland SE, Wray NR, Brown KM, Investigators AMFS, Hayward NK, Montgomery GW, Visscher PM, et al. A versatile gene-based test for genome-wide association studies. Am J Hum Genet. 2010;87(1):139–45.
GTEx Consortium. The Genotype-Tissue Expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science. 2015;348(6235):648–60.
Zou F, Chai HS, Younkin CS, Allen M, Crook J, Pankratz VS, Carrasquillo MM, Rowley CN, Nair AA, Middha S, et al. Brain expression genome-wide association study (eGWAS) identifies human disease-associated variants. PLoS Genet. 2012;8(6):e1002707.
Carrasquillo MM, Zou F, Pankratz VS, Wilcox SL, Ma L, Walker LP, Younkin SG, Younkin CS, Younkin LH, Bisceglio GD, et al. Genetic variation in PCDH11X is associated with susceptibility to late-onset Alzheimer’s disease. Nat Genet. 2009;41(2):192–8.
Allen M, Carrasquillo MM, Funk C, Heavner BD, Zou F, Younkin CS, Burgess JD, Chai HS, Crook J, Eddy JA, et al. Human whole genome genotype and transcriptome data for Alzheimer’s and other neurodegenerative diseases. Sci Data. 2016;3:160089.
Zhang B, Gaiteri C, Bodea LG, Wang Z, McElwee J, Podtelezhnikov AA, Zhang C, Xie T, Tran L, Dobrin R, et al. Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer’s disease. Cell. 2013;153(3):707–20.
Jun G, Ibrahim-Verbaas CA, Vronskaya M, Lambert JC, Chung J, Naj AC, Kunkle BW, Wang LS, Bis JC, Bellenguez C, et al. A novel Alzheimer disease locus located near the gene encoding tau protein. Mol Psychiatry. 2016;21(1):108–17.
Jun GR, Chung J, Mez J, Barber R, Beecham GW, Bennett DA, Buxbaum JD, Byrd GS, Carrasquillo MM, Crane PK, et al. Transethnic genome-wide scan identifies novel Alzheimer’s disease loci. Alzheimers Dement. 2017;13(7):727–38.
Lambert JC, Ibrahim-Verbaas CA, Harold D, Naj AC, Sims R, Bellenguez C, DeStafano AL, Bis JC, Beecham GW, Grenier-Boley B, et al. Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat Genet. 2013;45(12):1452–8.
Liu CC, Kanekiyo T, Xu H, Bu G. Apolipoprotein E and Alzheimer disease: risk, mechanisms and therapy. Nat Rev Neurol. 2013;9(2):106–18.
Minett T, Classey J, Matthews FE, Fahrenhold M, Taga M, Brayne C, Ince PG, Nicoll JA, Boche D, MRC CFAS. Microglial immunophenotype in dementia with Alzheimer’s pathology. J Neuroinflammation. 2016;13(1):135.
Rosenthal SL, Kamboh MI. Late-onset Alzheimer’s disease genes and the potentially implicated pathways. Curr Genet Med Rep. 2014;2:85–101.
Tan MS, Yu JT, Tan L. Bridging integrator 1 (BIN1): form, function, and Alzheimer’s disease. Trends Mol Med. 2013;19(10):594–603.
Li LW, Li YY, Li XY, Zhang CP, Zhou Y, Lu SH. A novel tumor suppressor gene ECRG4 interacts directly with TMPRSS11A (ECRG1) to inhibit cancer cell growth in esophageal carcinoma. BMC Cancer. 2011;11:52.
Woo JM, Park SJ, Kang HI, Kim BG, Shim SB, Jee SW, Lee SH, Sin JS, Bae CJ, Jang MK, et al. Characterization of changes in global gene expression in the brain of neuron-specific enolase/human Tau23 transgenic mice in response to overexpression of Tau protein. Int J Mol Med. 2010;25(5):667–75.
Kujuro Y, Suzuki N, Kondo T. Esophageal cancer-related gene 4 is a secreted inducer of cell senescence expressed by aged CNS precursor cells. Proc Natl Acad Sci U S A. 2010;107(18):8259–64.
Podvin S, Miller MC, Rossi R, Chukwueke J, Donahue JE, Johanson CE, Baird A, Stopa EG. The orphan C2orf40 gene is a neuroimmune factor in Alzheimer’s disease. JSM Alzheimers Dis Relat Dement. 2016;3(1):1020.
Podvin S, Gonzalez AM, Miller MC, Dang X, Botfield H, Donahue JE, Kurabi A, Boissaud-Cooke M, Rossi R, Leadbeater WE, et al. Esophageal cancer related gene-4 is a choroid plexus-derived injury response gene: evidence for a biphasic response in early and late brain injury. PLoS One. 2011;6(9):e24609.
Li LW, Yu XY, Yang Y, Zhang CP, Guo LP, Lu SH. Expression of esophageal cancer related gene 4 (ECRG4), a novel tumor suppressor gene, in esophageal cancer and its inhibitory effect on the tumor growth in vitro and in vivo. Int J Cancer. 2009;125(7):1505–13.
Mori Y, Ishiguro H, Kuwabara Y, Kimura M, Mitsui A, Kurehara H, Mori R, Tomoda K, Ogawa R, Katada T, et al. Expression of ECRG4 is an independent prognostic factor for poor survival in patients with esophageal squamous cell carcinoma. Oncol Rep. 2007;18(4):981–5.
Strahl BD, Allis CD. The language of covalent histone modifications. Nature. 2000;403(6765):41–5.
Sugo N, Oshiro H, Takemura M, Kobayashi T, Kohno Y, Uesaka N, Song WJ, Yamamoto N. Nucleocytoplasmic translocation of HDAC9 regulates gene expression and dendritic growth in developing cortical neurons. Eur J Neurosci. 2010;31(9):1521–32.
Millan MJ. The epigenetic dimension of Alzheimer’s disease: causal, consequence, or curiosity? Dialogues Clin Neurosci. 2014;16(3):373–93.
Agis-Balboa RC, Pavelka Z, Kerimoglu C, Fischer A. Loss of HDAC5 impairs memory function: implications for Alzheimer’s disease. J Alzheimers Dis. 2013;33(1):35–44.
Fischer A. Targeting histone-modifications in Alzheimer’s disease: what is the evidence that this is a promising therapeutic avenue? Neuropharmacology. 2014;80:95–102.
Cuadrado-Tejedor M, Oyarzabal J, Lucas MP, Franco R, Garcia-Osta A. Epigenetic drugs in Alzheimer’s disease. Biomol Concepts. 2013;4(5):433–45.
Benito E, Urbanke H, Ramachandran B, Barth J, Halder R, Awasthi A, Jain G, Capece V, Burkhardt S, Navarro-Sala M, et al. HDAC inhibitor-dependent transcriptome and memory reinstatement in cognitive decline models. J Clin Invest. 2015;125(9):3572–84.
International Stroke Genetics Consortium (ISGC), Wellcome Trust Case Control Consortium 2 (WTCC2), Bellenguez C, Bevan S, Gschwendtner A, Spencer CC, Burgess AI, Pirinen M, Jackson CA, Traylor M, et al. Genome-wide association study identifies a variant in HDAC9 associated with large vessel ischemic stroke. Nat Genet. 2012;44(3):328–33.
Pinto D, Delaby E, Merico D, Barbosa M, Merikangas A, Klei L, Thiruvahindrapuram B, Xu X, Ziman R, Wang Z, et al. Convergence of genes and cellular pathways dysregulated in autism spectrum disorders. Am J Hum Genet. 2014;94(5):677–94.
Tam GW, van de Lagemaat LN, Redon R, Strathdee KE, Croning MD, Malloy MP, Muir WJ, Pickard BS, Deary IJ, Blackwood DH, et al. Confirmed rare copy number variants implicate novel genes in schizophrenia. Biochem Soc Trans. 2010;38(2):445–51.
Haberland M, Arnold MA, McAnally J, Phan D, Kim Y, Olson EN. Regulation of HDAC9 gene expression by MEF2 establishes a negative-feedback loop in the transcriptional circuitry of muscle differentiation. Mol Cell Biol. 2007;27(2):518–25.
Morrison BE, Majdzadeh N, Zhang X, Lyles A, Bassel-Duby R, Olson EN, D’Mello SR. Neuroprotection by histone deacetylase-related protein. Mol Cell Biol. 2006;26(9):3550–64.
Salian-Mehta S, Xu M, McKinsey TA, Tobet S, Wierman ME. Novel interaction of class IIb histone deacetylase 6 (HDAC6) with class IIa HDAC9 controls gonadotropin releasing hormone (GnRH) neuronal cell survival and movement. J Biol Chem. 2015;290(22):14045–56.
Oram SW, Ai J, Pagani GM, Hitchens MR, Stern JA, Eggener S, Pins M, Xiao W, Cai X, Haleem R, et al. Expression and function of the human androgen-responsive gene ADI1 in prostate cancer. Neoplasia. 2007;9(8):643–51.
Scrivens PJ, Noueihed B, Shahrzad N, Hul S, Brunet S, Sacher M. C4orf41 and TTC-15 are mammalian TRAPP components with a role at an early stage in ER-to-Golgi trafficking. Mol Biol Cell. 2011;22(12):2083–93.
Reitz C, Jun G, Naj A, Rajbhandary R, Vardarajan BN, Wang LS, Valladares O, Lin CF, Larson EB, Graff-Radford NR, et al. Variants in the ATP-binding cassette transporter (ABCA7), apolipoprotein E ε4, and the risk of late-onset Alzheimer disease in African Americans. JAMA. 2013;309(14):1483–92.
Rogaeva E, Meng Y, Lee JH, Gu Y, Kawarai T, Zou F, Katayama T, Baldwin CT, Cheng R, Hasegawa H, et al. The neuronal sortilin-related receptor SORL1 is genetically associated with Alzheimer disease. Nat Genet. 2007;39(2):168–77.
Vardarajan BN, Bruesegem SY, Harbour ME, Inzelberg R, Friedland R, St George-Hyslop P, Seaman MN, Farrer LA. Identification of Alzheimer disease-associated variants in genes that regulate retromer function. Neurobiol Aging. 2012;33(9):2231. e15–e30
Pelechano V, Steinmetz LM. Gene regulation by antisense transcription. Nat Rev Genet. 2013;14(12):880–93.
Cannon TD, Keller MC. Endophenotypes in the genetic analyses of mental disorders. Annu Rev Clin Psychol. 2006;2:267–90.
Flint J, Munafo MR. The endophenotype concept in psychiatric genetics. Psychol Med. 2007;37(2):163–80.
The authors express appreciation to Dr. Qiong Yang for help with using the R library (‘CUMP’ package).
Members of the Alzheimer Disease Genetics Consortium
Marilyn S. Albert1, Roger L. Albin2−–4, Liana G. Apostolova5, Steven E. Arnold6, Sanjay Asthana7−–9, Craig S. Atwood7,9, Clinton T. Baldwin10, Michael M. Barmada11, Lisa L. Barnes12,13, Thomas G. Beach14, James T. Becker15, Duane Beekly16, David A. Bennett12, Eileen H. Bigio17,18, Thomas D. Bird19,20, Deborah Blacker21,22, Bradley F. Boeve23, James D. Bowen24, Adam Boxer25, James R. Burke26, Joseph Buxbaum108, Nigel J. Cairns27, Laura B. Cantwell28, Chuanhai Cao29, Chris S. Carlson30, Cynthia M. Carlsson8, Regina M. Carney31, Minerva M. Carrasquillo32, Steven L. Carroll33, Helena C. Chui34, David G. Clark35, Paul Crane80, David H. Cribbs36, Elizabeth A. Crocco31, Carlos Cruchaga37, Charles DeCarli38, F. Yesim Demirci11, Philip De Jager109, Malcolm Dick39, Dennis Dickson32, Ranjan Duara40, Nilufer Ertekin-Taner32,41, Kelley M. Faber42, Kenneth B. Fallon33, Martin R. Farlow43, Steven Ferris44, Tatiana M. Foroud42, Matthew Frosch110, Douglas R. Galasko45, Marla Gearing46,47, Daniel H. Geschwind48, Bernardino Ghetti85, John Gilbert70, Jonathan D. Glass49, Alison M. Goate37, Neill R. Graff-Radford32,41, Robert C. Green50, John H. Growdon51, Hakon Hakonarson52, Ronald L. Hamilton53, Kara L. Hamilton-Nelson70, John Hardy111, Lindy E. Harrell35, Elizabeth Head54, Lawrence S. Honig55, Matthew J. Huentelman112, Christine M. Hulette56, Bradley Hyman51, Gail P. Jarvik57,58, Gregory A. Jicha59, Lee-Way Jin60, Gyungah Jun10,61,62, M. Ilyas Kamboh11,75, Anna Karydas25, John S.K. Kauwe63, Jeffrey A. Kaye64,65, Ronald Kim66, Neil W. Kowall67,68, Joel H. Kramer69, Patricia Kramer64, Walter Kukull16,113, Brian W. Kunkle70, Frank M. LaFerla71, James J. Lah49, Eric Larson80, James B. Leverenz72, Allan I. Levey49, Ge Li73, Andrew P. Lieberman74, Chiao-Feng Lin28, Oscar L. Lopez75, Kathryn L. Lunetta61, Constantine G. Lyketsos76, Wendy J. Mack77, Daniel C. Marson35, Eden R. Martin17, Frank Martiniuk78, Deborah Mash107, Eliezer Masliah45,79, Wayne C. McCormick80, Susan M. McCurry81, Andrew N. McDavid30, Ann C. McKee67,68, Marsel Mesulam18,82, Bruce L. Miller25, Carol A. Miller83, Joshua W. Miller60, John C. Morris27,84, Shubhabrata Mukherjee80, Jill R. Murrell42,85, Amanda J. Myers31, Adam C. Naj114, John M. Olichney38, Vernon S. Pankratz86, Joseph E. Parisi87, Amanda Partch28, Henry L. Paulson88, William Perry70, Elaine Peskind73, Ronald C. Petersen23, Aimee Pierce36, Wayne W. Poon39, Huntington Potter89, Joseph F. Quinn64, Ashok Raj29, Murray Raskind73, Eric Reiman115,116, Barry Reisberg44,90, Christiane Reitz55,91,92, John M. Ringman5, Erik D. Roberson35, Ekaterina Rogaeva93, Howard J. Rosen25, Roger N. Rosenberg94, Mark A. Sager8, Mary Sano95, Andrew J. Saykin42,96, Julie Schneider12,117, Lon S. Schneider34,97, William W. Seeley25, Amanda G. Smith29, Joshua A. Sonnen98, Salvatore Spina85, Peter St George-Hyslop93,99, Robert A. Stern67, Rudolph E. Tanzi51, Tricia A. Thornton-Wells100, John Trojanowski28, Juan C. Troncoso101, Debby W. Tsuang20,73, Otto Valladares28, Vivianna M. Van Deerlin28, Linda J. Van Eldik102, Badri N. Vardarajan55,91,92, Harry V. Vinters5,103, Jean-Paul Vonsattel118, Li-San Wang28, Sandra Weintraub18,104, Kathleen A. Welsh-Bohmer26,105, Jennifer Williamson55, Sarah Wishnek70, Randall L. Woltjer106, Clinton B. Wright107, Steven Younkin32, Chang-En Yu80, and Lei Yu12
1Department of Neurology, Johns Hopkins University, Baltimore, Maryland; 2Department of Neurology, University of Michigan, Ann Arbor, Michigan; 3Geriatric Research, Education and Clinical Center (GRECC), VA Ann Arbor Healthcare System (VAAAHS), Ann Arbor, Michigan; 4Michigan Alzheimer Disease Center, Ann Arbor, Michigan; 5Department of Neurology, University of California, Los Angeles, Los Angeles, California; 6Department of Psychiatry, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania; 7Geriatric Research, Education and Clinical Center (GRECC), University of Wisconsin, Madison, Wisconsin; 8Department of Medicine, University of Wisconsin, Madison, Wisconsin; 9Wisconsin Alzheimer’s Institute, Madison, Wisconsin; 10Department of Medicine (Biomedical Genetics), Boston University, Boston, Massachusetts; 11Department of Human Genetics, University of Pittsburgh, Pittsburgh, Pennsylvania; 12Department of Neurological Sciences, Rush University Medical Center, Chicago, Illinois; 13Department of Behavioral Sciences, Rush University Medical Center, Chicago, Illinois; 14Civin Laboratory for Neuropathology, Banner Sun Health Research Institute, Phoenix, Arizona; 15Departments of Psychiatry, Neurology, and Psychology, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania; 16National Alzheimer’s Coordinating Center, University of Washington, Seattle, Washington; 17Department of Pathology, Northwestern University Feinberg School of Medicine, Chicago, Illinois; 18Cognitive Neurology and Alzheimer’s Disease Center, Northwestern University Feinberg School of Medicine, Chicago, Illinois; 19Department of Neurology, University of Washington, Seattle, Washington; 20VA Puget Sound Health Care System/GRECC, Seattle, Washington; 21Department of Epidemiology, Harvard School of Public Health, Boston, Massachusetts; 22Department of Psychiatry, Massachusetts General Hospital/Harvard Medical School, Boston, Massachusetts; 23Department of Neurology, Mayo Clinic, Rochester, Minnesota; 24Swedish Medical Center, Seattle, Washington; 25Department of Neurology, University of California, San Francisco, San Francisco, California; 26Department of Medicine, Duke University, Durham, North Carolina; 27Department of Pathology and Immunology, Washington University, St. Louis, Missouri; 28Department of Pathology and Laboratory Medicine, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania; 29USF Health Byrd Alzheimer’s Institute, University of South Florida, Tampa, Florida; 30Fred Hutchinson Cancer Research Center, Seattle, Washington; 31Department of Psychiatry and Behavioral Sciences, Miller School of Medicine, University of Miami, Miami, Florida; 32Department of Neuroscience, Mayo Clinic, Jacksonville, Florida; 33Department of Pathology, University of Alabama at Birmingham, Birmingham, Alabama; 34Department of Neurology, University of Southern California, Los Angeles, Los Angeles, California; 35Department of Neurology, University of Alabama at Birmingham, Birmingham, Alabama; 36Department of Neurology, University of California, Irvine, Irvine, California; 37Department of Psychiatry and Hope Center Program on Protein Aggregation and Neurodegeneration, Washington University School of Medicine, St. Louis, Missouri; 38Department of Neurology, University of California, Davis, Sacramento, California; 39Institute for Memory Impairments and Neurological Disorders, University of California, Irvine, Irvine, California; 40Wien Center for Alzheimer’s Disease and Memory Disorders, Mount Sinai Medical Center, Miami Beach, Florida; 41Department of Neurology, Mayo Clinic, Jacksonville, Florida; 42Department of Medical and Molecular Genetics, Indiana University, Indianapolis, Indiana; 43Department of Neurology, Indiana University, Indianapolis, Indiana; 44Department of Psychiatry, New York University, New York, New York; 45Department of Neurosciences, University of California, San Diego, La Jolla, California; 46Department of Pathology and Laboratory Medicine, Emory University, Atlanta, Georgia; 47Emory Alzheimer’s Disease Center, Emory University, Atlanta, Georgia; 48Neurogenetics Program, University of California, Los Angeles, Los Angeles, California; 49Department of Neurology, Emory University, Atlanta, Georgia; 50Division of Genetics, Department of Medicine and Partners Center for Personalized Genetic Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, Massachusetts; 51Department of Neurology, Massachusetts General Hospital/Harvard Medical School, Boston, Massachusetts; 52Center for Applied Genomics, Children’s Hospital of Philadelphia, Philadelphia, Pennsylvania; 53Department of Pathology (Neuropathology), University of Pittsburgh, Pittsburgh, Pennsylvania; 54Sanders-Brown Center on Aging, Department of Molecular and Biomedical Pharmacology, University of Kentucky, Lexington, Kentucky; 55Taub Institute on Alzheimer’s Disease and the Aging Brain, Department of Neurology, Columbia University, New York, New York; 56Department of Pathology, Duke University, Durham, North Carolina; 57Department of Genome Sciences, University of Washington, Seattle, Washington; 58Department of Medicine (Medical Genetics), University of Washington, Seattle, Washington; 59Sanders-Brown Center on Aging, Department Neurology, University of Kentucky, Lexington, Kentucky; 60Department of Pathology and Laboratory Medicine, University of California, Davis, Sacramento, California; 61Department of Biostatistics, Boston University, Boston, Massachusetts; 62Department of Ophthalmology, Boston University, Boston, Massachusetts; 63Department of Biology, Brigham Young University, Provo, Utah; 64Department of Neurology, Oregon Health & Science University, Portland, Oregon; 65Department of Neurology, Portland Veterans Affairs Medical Center, Portland, Oregon; 66Department of Pathology and Laboratory Medicine, University of California, Irvine, Irvine, California; 67Department of Neurology, Boston University, Boston, Massachusetts; 68Department of Pathology, Boston University, Boston, Massachusetts; 69Department of Neuropsychology, University of California, San Francisco, San Francisco, California; 70The John P. Hussman Institute for Human Genomics, University of Miami, Miami, Florida; 71Department of Neurobiology and Behavior, University of California, Irvine, Irvine, California; 72Cleveland Clinic Lou Ruvo Center for Brain Health, Cleveland Clinic, Cleveland, Ohio; 73Department of Psychiatry and Behavioral Sciences, University of Washington School of Medicine, Seattle, Washington; 74Department of Pathology, University of Michigan, Ann Arbor, Michigan; 75University of Pittsburgh Alzheimer’s Disease Research Center, Pittsburgh, Pennsylvania; 76Department of Psychiatry, Johns Hopkins University, Baltimore, Maryland; 77Department of Preventive Medicine, University of Southern California, Los Angeles, California; 78Department of Medicine–Pulmonary, New York University School of Medicine, New York, New York; 79Department of Pathology, University of California, San Diego, La Jolla, California; 80Department of Medicine, University of Washington, Seattle, Washington; 81Northwest Research Group on Aging, University of Washington School of Nursing, Seattle, Washington; 82Department of Neurology, Northwestern University Feinberg School of Medicine, Chicago, Illinois; 83Department of Pathology, University of Southern California, Los Angeles, Los Angeles, California; 84Department of Neurology, Washington University, St. Louis, Missouri; 85Department of Pathology and Laboratory Medicine, Indiana University, Indianapolis, Indiana; 86Department of Biostatistics, Mayo Clinic, Rochester, Minnesota; 87Department of Laboratory Medicine and Pathology, Mayo Clinic, Rochester, Minnesota; 88Michigan Alzheimer’s Disease Center, Department of Neurology, University of Michigan, Ann Arbor, Michigan; 89Department of Neurology, University of Colorado School of Medicine, Aurora, Colorado; 90Alzheimer’s Disease Center, New York University, New York, New York; 91Gertrude H. Sergievsky Center, Columbia University, New York, New York; 92Department of Neurology, Columbia University, New York, New York; 93Tanz Centre for Research in Neurodegenerative Disease, University of Toronto, Toronto, Ontario, Canada; 94Department of Neurology, University of Texas Southwestern, Dallas, Texas; 95Department of Psychiatry, Mount Sinai School of Medicine, New York, New York; 96Department of Radiology and Imaging Sciences, Indiana University, Indianapolis, Indiana; 97Department of Psychiatry, University of Southern California, Los Angeles, California; 98Department of Pathology, University of Washington, Seattle, Washington; 99Cambridge Institute for Medical Research and Department of Clinical Neurosciences, University of Cambridge, Cambridge, United Kingdom; 100Center for Human Genetics and Research, Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, Tennessee; 101Department of Pathology, Johns Hopkins University, Baltimore, Maryland; 102Sanders-Brown Center on Aging, Department of Anatomy and Neurobiology, University of Kentucky, Lexington, Kentucky; 103Department of Pathology & Laboratory Medicine, University of California, Los Angeles, Los Angeles, California; 104Department of Psychiatry, Northwestern University Feinberg School of Medicine, Chicago, Illinois; 105Department of Psychiatry & Behavioral Sciences, Duke University, Durham, North Carolina; 106Department of Pathology, Oregon Health & Science University, Portland, Oregon; 107Evelyn F. McKnight Brain Institute, Department of Neurology, Miller School of Medicine, University of Miami, Miami, Florida; 108Department of Psychiatry, Mount Sinai Hospital, New York, New York; 109Department of Neurology and Psychiatry, Brigham and Women’s Hospital/Harvard Medical School, Boston, Massachusetts; 110C.S. Kubik Laboratory for Neuropathology, Massachusetts General Hospital, Charlestown, Massachusetts; 111Department of Molecular Neuroscience, University College London, London, United Kingdom; 112Neurogenomics Division, Translational Genomics Research Institute, Phoenix, Arizona; 113Department of Epidemiology, University of Washington, Seattle, Washington; 114Division of Epidemiology, University of Pennsylvania, Philadelphia, Pennsylvania; 115Arizona Alzheimer’s Consortium, Banner Alzheimer’s Institute, Phoenix, Arizona; 116Department of Psychiatry, University of Arizona, Phoenix, Arizona; 117Department of Pathology (Neuropathology), Rush University Medical Center, Chicago, Illinois; 118New York Brain Bank, Columbia University, New York, New York.
The Alzheimer’s Disease Genetics Consortium supported the collection of samples used in this study through National Institute on Aging (NIA) grants U01-AG032984 and RC2AG036528. Data for this study were prepared, archived, and distributed by the National Institute on Aging Alzheimer’s Disease Data Storage Site (NIAGADS) at the University of Pennsylvania and funded by NIA grant U24-AG041689-01. This work was also supported by NIA grants R01-AG048927 (to LAF), R01-AG032990 and RF-AG051504 (to NET), and U01-AG046139 (to NET, SGY, TEG, and NDP), as well as by National Institute of Neurological Disorders and Stroke (NINDS) grant R01-NS080820 (to NET). The Brain and Body Donation Program at Sun Health Research Institute (Sun City, AZ, USA) is supported by NINDS grant U24-NS072026 for the National Brain and Tissue Resource for Parkinson’s Disease and Related Disorders; NIA grant P30-AG19610 for the Arizona Alzheimer’s Disease Core Center; Arizona Department of Health Services contract 211002 for the Arizona Alzheimer’s Research Center; Arizona Biomedical Research Commission contracts 4001, 0011, 05-901, and 1001 to the Arizona Parkinson’s Disease Consortium; and the Michael J. Fox Foundation for Parkinson’s Research.
Availability of data and materials
Summarized genome-wide association study results are available from the National Institute on Aging Genetics of Alzheimer’s Disease Data Storage Site (https://www.niagads.org/).
Ethics approval and consent to participate
This study was approved by the Boston University Institutional Review Board.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Sample characteristics. Table S2. Association P values of Alzheimer disease loci previously established by GWAS in univariate and pleiotropy association tests of neuropathological features. Table S3. Association of cis-eQTL for HDAC9 in the Mayo Clinic brain expression genome-wide association study (eGWAS). rs79524815 was not available in the Mayo Clinic brain eGWAS, so proxy SNPs that are in LD (D′ > 0.90) with rs79524815 were used for the eQTL test with HDAC9 expression. Table S4. Association of expression of SNPs for TRAPPC12-AS1 and ADI1 with neuropathological traits and gene expression in the GTEx portal database. Table S5. Association results from the trivariate pleiotropy model of neuritic plaque (NP), neurofibrillary tangles (NFT), and cerebral amyloid angiopathy (CAA) for study-wide significant SNPs in the bivariate pleiotropy model. Table S6. Results of differential gene expression analysis by brain region among AD cases and controls for AD loci previously established by GWAS in RNA-Seq and microarray analysis. Figure S1. Quantile-quantile plots of observed (y-axis) vs. expected (x-axis) P values of all SNPs (black dots) and after excluding SNPs in APOE region (blue dots) for the pleiotropy analysis of (a) NP and NFT, (b) NP and CAA, and (c) NFT and CAA using the O’Brien method . Figure S2. Manhattan plots showing genome-wide pleiotropy analyses of (a) NP and NFT, (b) NP and CAA, and (c) NFT and CAA using the O’Brien method . Red dashed line represents the genome-wide significance threshold of P < 5.0 × 10−−8. Loci achieving genome-wide significance are highlighted in red, and known AD genes that attained at least a moderate significance level (P < 10−−4) are highlighted in gold. Figure S3. Regional association plots of genes, including TRAPPC12, TRAPPC12-AS1, and ADI1, on chromosome 2 from the joint model of NFT and CAA. Figure S4. Genome-wide trivariate pleiotropy analysis of NP, NFT, and CAA. (a) Quantile-quantile plot. (b) Manhattan plot. (DOCX 1028 kb)