Cell-type Specific Expression Quantitative Trait Loci Associated with Alzheimer Disease in Blood and Brain Tissue
=================================================================================================================

* Devanshi Patel
* Xiaoling Zhang
* John J. Farrell
* Jaeyoon Chung
* Thor D. Stein
* Kathryn L. Lunetta
* Lindsay A. Farrer

## ABSTRACT

Because regulation of gene expression is heritable and context-dependent, we investigated AD-related gene expression patterns in cell-types in blood and brain. Cis-expression quantitative trait locus (eQTL) mapping was performed genome-wide in blood from 5,257 Framingham Heart Study (FHS) participants and in brain donated by 475 Religious Orders Study/Memory & Aging Project (ROSMAP) participants. The association of gene expression with genotypes for all cis SNPs within 1Mb of genes was evaluated using linear regression models for unrelated subjects and linear mixed models for related subjects. Cell type-specific eQTL (ct-eQTL) models included an interaction term for expression of “proxy” genes that discriminate particular cell type. Ct-eQTL analysis identified 11,649 and 2,533 additional significant gene-SNP eQTL pairs in brain and blood, respectively, that were not detected in generic eQTL analysis. Of note, 386 unique target eGenes of significant eQTLs shared between blood and brain were enriched in apoptosis and Wnt signaling pathways. Five of these shared genes are established AD loci. The potential importance and relevance to AD of significant results in myeloid cell-types is supported by the observation that a large portion of GWS ct-eQTLs map within 1Mb of established AD loci and 58% (23/40) of the most significant eGenes in these eQTLs have previously been implicated in AD. This study identified cell-type specific expression patterns for established and potentially novel AD genes, found additional evidence for the role of myeloid cells in AD risk, and discovered potential novel blood and brain AD biomarkers that highlight the importance of cell-type specific analysis.

## INTRODUCTION

Recent expression quantitative trait locus (eQTL) analysis studies suggest that changes in gene expression have a role in the pathogenesis of AD 1, 2. However, regulation of gene expression, as well as many biological functions, has been shown to be context-specific (e.g., tissue and cell-types, developmental time point, sex, disease status, and response to treatment or stimulus) 3-6. One study of 500 healthy subjects found over-representation of T cell-specific eQTLs in susceptibility alleles for autoimmune disease and AD risk alleles polarized for monocyte-specific eQTL effects 7. In addition, disease and trait-associated cis-eQTLs were more cell type specific than average cis-eQTLs 7. Another study classified 12% of more than 23000 eQTLs in blood as cell-type specific 5. Large eQTL studies across multiple human tissues have been conducted by the GTEx consortium, with a study on genetic effects on gene expression levels across 44 human tissues collected from the same donors characterizing patterns of tissue specificity recently published 8.

Microglia, monocytes and macrophages share a similar developmental lineage and are all considered to be myeloid cells 9. It is believed that a large proportion of AD genetic risk can be explained by genes expressed in myeloid cells and not other cell-types 10. Several established AD genes are highly expressed in microglia 9, 11 and a variant in the AD-associated locus *CELF1* has been associated with lower expression of *SPI1* in monocytes and macrophages 10. AD risk alleles have been shown to be enriched in myeloid specific epigenomic annotations and in active enhancers of monocytes, macrophages, and microglia 12, and to be polarized for cis-eQTL effects in monocytes 7. These findings suggest that a cell-type specific analysis in blood and brain tissue may identify novel and more precise AD associations that may help elucidate regulatory networks. In this study, we performed a genome-wide *cis* ct-eQTL analysis in blood and brain, respectively, then compared eQTLs and cell-type specific eQTLs (ct-eQTLs) between brain and blood with a focus on genes, loci, and cell-types previously implicated in AD risk by genetic approaches.

## MATERIALS, SUBJECTS AND METHODS

### Study cohorts

#### Framingham Heart Study (FHS)

The FHS is a multigenerational study of health and disease in a prospectively followed community-based and primarily non-Hispanic white sample. Procedures for assessing dementia and determining AD status in this cohort are described elsewhere 13. Clinical, demographic, and pedigree information, as well as 1000 Genomes Project Phase 1 imputed SNP genotypes and Affymetrix Human Exon 1.0 ST array gene expression data from whole blood, were obtained from dbGaP ([https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000007.v31.p12](https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000007.v31.p12)). Requisite information for this study was available for 5,257 participants. Characteristics of these subjects are provided in Table S1.

#### Religious Orders Study (ROS)/ Memory and Aging Project (MAP)

ROS enrolled older nuns and priests from across the US, without known dementia for longitudinal clinical analysis and brain donation and MAP enrolled older subjects without dementia from retirement homes who agreed to brain donation at the time of death 14, 15. RNA-sequencing brain gene expression and whole-genome sequencing (WGS) genotype data were obtained from the AMP-AD knowledge portal ([https://www.synapse.org/#!Synapse:syn3219045](https://www.synapse.org/#!Synapse:syn3219045)) 16.

#### Data processing

Generation and initial quality control (QC) procedures of the FHS GWAS and expression data are described elsewhere and include all genotype QC and pre-adjustment of gene expression levels for batch effects and other technical covariates 13. ROSMAP gene expression data were log-normalized and adjusted for known and hidden variables detected by surrogate variable analysis (SVA) 17 in order to determine which of these variables should be included as covariates in analysis models for eQTL discovery. Additional filtering steps of FHS and ROSMAP GWAS and gene expression data included eliminating subjects with missing data, restricting gene expression data to protein coding genes, and retaining common variants (MAF□≥ 0.05) with good imputation quality (R2□≥□0.3).

#### Cis eQTL mapping

Cis-eQTL mapping was performed using a genome-wide design (Fig. S1). The association of gene expression with SNP genotypes for all cis SNPs within1 Mb of protein-coding genes was evaluated using linear mixed models adjusting for family structure in FHS and linear regression models for unrelated individuals in ROSMAP. In FHS, lmekin function in the R kinship package (version 1.1.3) 18 was applied assuming an additive genetic model with covariates for age and sex, and family structure modeled as a random-effects term for kinship - a matrix of kinship coefficients calculated from pedigree structures. The linear model for analysis of FHS can data be expressed as follows: ![Formula][1]</img>  where *Y**i* is the expression value for gene *i, G**j* is the genotype dosage for cis SNP j, *Aij* and *S**ij* are the covariates for age and sex respectively, *U**ij* is the random effect for family structure, and *β**1*, *β**2*, and *β**3* are regression coefficients.

ROSMAP data were analyzed using the lm function in the base stats package in R 19. The regression model, which included covariates for age, sex, post-mortem interval (PMI), study (ROS or MAP), and a term for a surrogate variable (SV1) derived from analysis of high dimensional data, can be expressed as: ![Formula][2]</img>  where *Y**i* is the expression value for gene *i, G**j* is the genotype dosage for cis SNP j, *Aij, S**ij*, *PM**ij*, *S2**ij*, and *SV1**ij* are the covariates for age, sex, PMI, study and SV1 respectively, □ij is the residual error, and the *β*s are regression coefficients.

#### Cis ct-eQTL mapping

Models testing associations with cell type-specific eQTLs (ct-eQTLs) included an interaction term for expression levels of “proxy” genes that represent cell types. Proxy genes representing 10 cell types in whole blood 5 and five cell types in brain 20-22 were incorporated in cell type-specific models (Table S2). These proxy genes for cell types in blood were established previously using BLUEPRINT expression data to validate cell-type-specific expression in each cell-type module 5 and the proxy genes for brain cell types have been incorporated in several studies 20-22. Cell type-specific expression analyses in blood of FHS participants were conducted using the following model: ![Formula][3]</img>  where in each eQTLij pair, *Y**i* is the eQTL expression value for gene *i, G**j* is the genotype dosage for cis SNP j, P is the proxy gene, **P**G******j*** is the interaction term representing the effect of genotype in a particular cell type, *Aij* and *S**ij* are covariates for age and sex respectively, *U**ij* is the random effect for family structure, and *β*s are regression coefficients. Models with significant interaction terms indicate cell type specific eQTLs.

The following model was used to evaluate cell type-specific expression in brain in ROSMAP: ![Formula][4]</img>  where in each eQTLij pair, variables *Y**i*, *G**j*, P, *Aij, S**ij*, *P**ij*, □ij and *β*s are as described above, and *PM**ij*, *S2**ij*, and *SV1**ij* are covariates for PMI, study, and SV1 respectively.

A Bonferroni correction was applied to determine the significance threshold for each analysis (Table S3).

#### Selection of eQTLs in established AD loci and gene-set pathway enrichment analysis

eGenes (genes whose expression levels are associated with variation at a particular eSNP) were matched to 88 genes located near 80 distinct uncorrelated SNPs that have been associated with AD or AD-related traits by genetic association or experimental approaches (Table S4) and eSNPs (SNPs that significantly influence gene expression) under the 80 significant association peaks. Gene-set enrichment analysis was performed using the PANTHER (Protein ANalysis THrough Evolutionary Relationships) software tool 23 to determine if the unique genes in the significant eQTL/ct-eQTL pairs shared by both brain and blood datasets are associated with a specific biological process or molecular function. Significance of the pathways was determined by the Fisher’s Exact test with False discovery rate (FDR) multiple test correction.

#### Colocalization analyses

Assessment of causal variants shared by adjacent GWAS and eQTL signals was performed using a Bayesian colocalization approach implemented in the R package *coloc* 24. This analysis incorporated information about significantly associated variants for AD risk obtained from a recent large GWAS 25 and lead eQTL variants each defined as the eSNP showing the strongest association with gene expression. Following recommended guidelines, the variants were deemed to be colocalized by a high posterior probability that a single shared variant is responsible for both signals (PP4 > 0.8) 24, 26. A lower threshold for statistical significance with a false discovery rate (FDR) < 0.05 for eQTL significant results was applied to maximize detection of colocalized pairs. Regional plots were constructed with LocusZoom 27.

## RESULTS

A total of 847,429 eQTLs and 30,405 ct-eQTLs in blood, and 173,857 eQTLs and 51,098 ct-eQTLs in brain were significant after Bonferroni correction (Table S3 and Supplemental Resources). Among these results, 11,649 ct-eQTLs in brain 2,533 ct-eQTLs in blood involved SNP-eGene pairs that were not detected in eQTL analysis (Fig. 1A). Of note, 24,028 significant eQTLs were shared between blood and brain. The 386 distinct eGenes among these shared eQTLs (Table S5) are most enriched in the apoptosis signaling *(P=*0.023) and Wnt signaling (*P*=0.036) pathways (Table S6). Five of these eGenes (*HLA-DRB5, HLA-DRB1, ECHDC3, CR1*, and *WWOX*) were previously associated with AD 25, 28. Three eSNPs in eQTLs involving *HLA-DRB1*/*HLA-DRB5* (rs9271058) and *ARL17A*/*LRRC37A2* (rs2732703 and rs113986870, which are near *KANSL1* and *MAPT*) were previously associated with AD risk at the genome-wide significance level 25, 29.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/24/2020.11.23.20237008/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/11/24/2020.11.23.20237008/F1)

Figure 1. Significant gene-SNP eQTLs and ct-eQTLs in blood and brain tissue genome-wide.
**A)** Venn diagram shows the number of overlapping eQTLs and ct-eQTLs in blood and brain. Gold color indicates significant eQTLs that are cell-type specific. Orange color indicates significant eQTLs that are shared between blood and brain. **B)** Cell-type distributions of significant genome-wide ct-eQTL results in blood and brain.

eQTLs involving *CR1, ECHDC3* and *WWOX* were much more significant in brain than blood, whereas *HLA-DRB5* and *HLA-DRB1* were more significant in blood when comparing the effect sizes (Table 1). *ECHDC3* was a significant eGene in blood and brain eQTLs (specifically in neurons). *HLA-DRB5* and *HLA-DRB1* were the only eGenes ascribed to significant ct-eQTLs in both blood and brain noting that of the 10 distinct lead eSNPs, five are unique to each tissue. Although the eQTLs involving these genes with the largest effect were observed in blood across multiple cell types, the total number of significant eSNP-eGene combinations was far greater in brain (particularly in microglia and neurons). The only instance in which the lead eSNP is also associated with AD risk at the GWS level was observed in the blood eQTL pair of *HLA-DRB1* with eSNP rs9271058 (Table 1A). Among the AD-associated SNPs at the GWS level, rs9271058 is a significant eSNP for *HLA-DRB1* in both blood and brain cell types (the most significant association by p-value was observed in anti-bacterial cells and microglia) and rs9271192 is a significant ct-eQTL for the gene in multiple brain cell types (Table 1). Both of these SNPs are also eSNPs for *HLA-DB5* in the brain in neurons only.

View this table:
[Table 1:](http://medrxiv.org/content/early/2020/11/24/2020.11.23.20237008/T1)

Table 1: 
eQTLs and ct-eQTLs in established AD loci appearing in both blood and brain

There were 657 gene-SNP eQTL pairs comprising 16 unique eGenes that were significant in blood and brain overall as well as in specific cell types in both blood and brain (Table S7). None of these eGenes were observed in significant pathways enriched for AD genes, however, they included AD-associated genes *HLA-DRB1* and *HLA-DRB5*.

Slightly more than half (42/80 = 52.5%) of the established AD associations (Table S3) are eGene targets for significant eQTLs in blood (Table S8). By comparison, only seven established AD loci were eGene targets for significant eQTLs in brain, among which *OARD1* was significant in endothelial cells only (Table S8). Many GWS SNPs for AD risk are eSNPs affecting expression of the nearest gene, which is usually recognized as the causative gene, but several GWS SNPs target other genes (Table S9). For example, AD-associated eSNPs rs113986870 and rs2732703 in the *MAPT/KANSL1* region target *ARL17A* in blood, but are paired in seven of eight eQTLs and ct-eQTLs with *LRRC37A2* in brain (Table S9). *HLA-DRB1* is the only AD gene with a significant ct-eQTL in blood, whereas many AD genes have significant blood eQTLs. In brain, only four AD loci (*CR1, HLA-DRB1/DRB5, IQCK* and *MAPT/KANSL1*) have significant brain eQTLs of which *HLA-DRB1/DRB5* and *MAPT/KANSL1* are the only brain ct-eQTLs, noting that all are significant in microglia, neurons and endothelial cells.

Next, we evaluated whether the most significant eSNPs and SNPs genome-wide significantly associated with AD status (i.e. AD-SNPs) co-localize and thus to identify a single shared variant responsible for both signals (posterior probability of shared signals (PP4) > 0.8). This analysis revealed eight eQTL/ct-eQTL signals that colocalized with seven AD GWAS signals and half of the co-localized signals involved a ct-eQTL (Table 2 and Fig. S2). Two different eSNPs for *CD2AP*, rs4711880 (eQTL *P*=1.4×10−104) and rs13201473 (ct-eQTL *P*= 1.47×10−9), flank *CD2AP* GWAS SNP rs10948363 which is also the second most significant eQTL (*P*=2.32×10−104) and the second most significant ct-eQTL in NK cells / CD8+ T-Cells (*P*=2.66×10−9). These three SNPs span a 9.0 kb region in intron 2 and are in complete linkage disequilibrium (LD, r2□=1.0), indicating that any one or more of them could affect the function of target gene CD2AP. Rs6557994 is the most significant eSNP for and located in *PTK2B* (blood ct-eQTL *P*=2.58×10−9) and is moderately correlated with the *PTK2B* GWAS SNP (rs28834970, r2=0.78, *P*=1.58×10−9). Thus, it is not surprising that rs6557994 is also significantly associated with AD risk (*P*=8.19×10−7). Rs6557994 is also correlated with a GWAS SNP in *CLU*, located approximately 150 kb from PTK2B, that is not significantly associated with expression of any gene. Because *PTK2B* and *CLU* are independent AD risk loci 28, it is possible that this eSNP has an effect on AD pathogenesis through independent pathways (Fig. S2). The most significant eSNP in *MADD* (rs35233100, *P*=2.88×10−10) was predicted to have functional consequences because it is a stop-gained mutation. This brain eQTL is colocalized (PP4=0.95) and weakly correlated with a GWAS SNP (*P*=1.91×10−5) in *CELF1* rs10838725 (r2=0.12).

View this table:
[Table 2:](http://medrxiv.org/content/early/2020/11/24/2020.11.23.20237008/T2)

Table 2: 
Colocalized AD GWAS/lead eQTL SNP pairs

Examination of the distribution of the significant ct-eQTL results genome-wide showed that nearly two-thirds of the ct-eQTLs in blood occurred in interferon response/anti-bacterial cells in blood, whereas brain ct-eQTLs are highly represented in endothelial cells, neurons and microglia (Fig. 1B, Table S10). Further examination of significant results within myeloid cell lineages (i.e., microglia and monocytes/ macrophages) which account for a large proportion of the genetic risk for late-onset AD 10 revealed that 3,234 or 10.6% of all significant ct-eQTLs in blood were in monocytes/macrophages. This subset includes 128 unique eGenes which are significantly enriched in the AD amyloid secretase pathway (FDR *P*=0.013, Table S11). A total of 974 or 30.1% of ct-eQTLs including 4 of the 20 most significant eGenes in monocytes/macrophages are located within 1 Mb of established AD loci (Table 4A). One of the eGenes in this top-ranked group (*HLA-DRB5*) is an established AD locus, and three others that are near established AD loci (*DLG2* near *PICALM* 30, *C4BPA* near *CR1* 31, and *MYO1E* near *ADAM10* 32) are reasonable AD gene candidates based on evidence using non-genetic approaches. Microglia accounted for 15,560 (30.5%) of significant ct-eQTLs in the brain (Table S10) which involved 304 unique eGenes. Approximately 52% of significant ct-eQTLs in microglia are located in AD regions including five of the 20 most significant ct-eQTLs in this group (Table 4B). One of these five eGenes is an established AD locus (*HLA-DRB1*) and two others (*ALCC* 33 and *WNT3* 34) have been linked to AD in previous studies.

View this table:
[Table 3:](http://medrxiv.org/content/early/2020/11/24/2020.11.23.20237008/T3)

Table 3: 
Top-ranked ct-eQTLs in myeloid cell types

View this table:
[Table 4:](http://medrxiv.org/content/early/2020/11/24/2020.11.23.20237008/T4)

Table 4: Overlap of ct-eQTLs in myeloid cell types in brain and blood

Considering significant eGene-eSNP pairs in myeloid cell types, 251 pairs including five distinct eGenes (*BTNL3, FAM118A, HLA-DOB, HLA-DRB1*, and *HLA-*DRB5) are shared between microglia and monocytes/macrophages (Table 4A and Fig. 2A). Three of these pairs involving eSNPs rs3763355, rs3763354, and rs1183595100 have the same target gene *HLA-DOB* and occur only in microglia and monocytes/macrophages (Table 4B).

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/24/2020.11.23.20237008/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/11/24/2020.11.23.20237008/F2)

Figure 2. Intersection of significant gene-SNP eQTL pairs between cell-types in blood and brain tissue.
**A)** Venn diagram showing overlap of ct-eQTL pairs in myeloid cell types (microglia and monocytes/macrophages). **B)** Number of significant eQTLs unique to and that overlap cell-types in blood and brain. Bar chart on the left side indicates the number of significant eQTLs involving each cell-type and bar chart above the matrix indicates the number of significant eQTLs that are unique to each cell type and set of cell-types. Pink colored bar indicates the number of eQTLs pairs that are unique to microglia and monocytes/macrophages.

Among the significant ct-eQTLs in brain, the cell types with the largest proportion that were also significant in monocytes/macrophages were microglia (1.6%) and neurons (1.3%) (Table 4C). Conversely, among the significant ct-eQTLs in blood, the cell types with the largest proportion that were also significant in microglia were NK/CD+ T-cells (12.9%) and monocytes/macrophages (7.8%). Among ct-eQTLs which are significant only for one cell-type each in blood and one in brain, monocytes/ macrophages shared three ct-eQTLs with microglia but with no other brain cell-types (Fig. 2B, Table 4C). By comparison, microglia shared 63 ct-eQTLs with interferons/ anti-bacterial cells, but with no other blood cell types. The much larger number of ct-eQTLs in microglia that were common with interferons/bacterial cells than monocytes/ macrophages may reflect the substantially greater proportion of significant eQTLs in blood involving interferons/ antibacterial cells (64%) than monocytes/macrophages (10.6%) (Table S10). The only other ct-eQTLs that were unique to a pair of cell types in brain and blood cell type involved neurons paired with neutrophils (n=3) and with interferons/anti-bacterial cells (n=65) (Fig. 2B).

## DISCUSSION

We identified several novel AD-related eQTLs that highlight the importance of cell-type dependent context. It is noteworthy that there were more significant ct-eQTLs in brain (n=51,098) than blood (n=30,405) even though the dataset containing expression data from blood (FHS) is several times larger than the brain expression dataset (ROSMAP). This could be due to greater cell type heterogeneity in brain, the enrichment of AD cases in the ROSMAP dataset who may show different patterns of gene expression compared to persons without AD, or highly variable gene expression across cell-types in the nervous system 35. Because expression studies in brain are often constrained by the small number specimens compared to studies in other tissues, post-mortem changes that may affect gene expression in brain 36, and the growing recognition that AD is a systemic disease 37-39, incorporating expression data from multiple tissues can enhance discovery of additional genetic influences on AD risk and pathogenesis.

Although most significant findings were tissue-specific, the 386 distinct eGenes among more than 24,000 significant gene-SNP eQTL pairs that were shared between blood and brain were enriched in the apoptosis signaling pathway that contributes to much of the underlying pathology associated with AD 40, 41. Five established AD genes (*CR1, ECHDC3, HLA-DRB1, HLA-DRB5*, and *WWOX* 25, 28) were shared eGenes in brain and blood and could be playing a key role in the systemic AD mechanisms. The complement receptor 1 (*CR1*) gene encodes a transmembrane glycoprotein functioning in the innate immune system by promoting phagocytosis of immune complexes, cellular debris, and Aβ 42. *CR1* is an eGene for several eSNPs, including AD GWAS peak SNP rs6656401 located within the gene, in brain and blood eQTLs and the effects on *CR1* expression are opposite in blood and brain. There are multiple possible explanations for the effect direction differences across tissues. The effect of eSNP rs6656401 on *CR1* expression may be developmental, noting that the average age of the FHS subjects (group with expression data in blood) is more than 30 years younger than the ROSMAP subjects (group with expression data in brain). The difference between brain and blood may also reflect post-mortem changes in brain that are not indicative of expression *in vivo*.

Alternatively, these effects may be related to AD because few FHS subjects were AD cases at the time of blood draw whereas 60% of subjects in the ROSMAP sample are AD cases. This idea is supported by the observation of a larger and positive effect of rs6656401 on *CR1* expression in AD (β=0.020) compared to control brains (β=-0.0086). Opposite effect directions of expression in brain and blood from AD patients compared to controls has been observed for several ribosomal genes 43. GWS variants located in the region spanning *ECHDC3* and *USP6NL* have previously been associated with AD 44. We found that *ECHDC3* is the target gene for eSNP rs866770710 located in its promoter region, and this eQTL was significant in brain and specifically in neurons. Altered *ECHDC3* expression in AD brains 45 supports the idea that this gene has a role in AD. Knockout of *WWOX* in mice leads to aggregation of amyloid-β (Aβ) and Tau, and subsequent cell death 46, 47.

The human leukocyte antigen (HLA) region is the key susceptibility locus in many immunological diseases and many associations have been reported between neurodegenerative diseases and HLA haplotypes 48. In addition, the most widely used marker to examine activated microglia in normal and diseased human brains is *HLA-DR* and microglia activation increases with the progression of AD 49, 50. *HLA-DRB5* and *HLA-DRB1* have been implicated in numerous GWAS studies as significantly associated with AD risk 25, 28 and appeared frequently among significant results in blood and brain in this study. Rs9271058, which is located approximately 17.8 kb upstream of *HLA-DRB1*, is significantly associated with AD risk (p=5.1 × 10−8 25) and when paired with *HLA-DRB1* is a significant eQTL and ct-eQTL in multiple cell types in blood and brain including myeloid lineage cells (i.e., monocytes/macrophages and microglia). This eSNP is also a significant eQTL in brain and specifically in neurons when paired with *HLA-DRB5*. Rs9271192, which is adjacent to rs9271058 and also significantly associated with AD risk (*P*=2.9 × 10−12) 28, is a significant eQTL and ct-QTL with multiple cell types in brain but not blood when paired with *HLA-DRB5 and HLA-DRB1*.

Significant associations for AD have been reported with variants spanning a large portion of the major histocompatibility (MHC) region in HLA class I, II and III loci 48, 51, 52. While the strongest statistical evidence for association in this region is with variants in *HLA-DRB1* 25, fine mapping in this region suggests that a class I haplotype (spanning the *HLA-A* and *HLA-B* loci) and a class II haplotype (including variants in *HLA-DRB1, HLA-DQA1* and *HLA-DQB1*) are more precise markers of AD risk. Given the complexity of the MHC region and extensive LD, further work is needed to confirm whether this is a true eQTL or a signal generated from a specific HLA allele or HLA haplotype. Although functional studies may be required to discern which HLA variants have AD-relevant consequences relevant to AD and develop methods to detect the differential expression of the HLA alleles, our findings support a role for the immune system in AD 37, 53 and the hypothesis that a large proportion of AD risk can be explained by genes expressed in myeloid cells 10.

The potential importance and relevance to AD of eQTLs and ct-eQTLs in myeloid cell-types is supported by the observation that a large portion of GWS ct-eQTLs we identified map within 1 Mb of established AD loci, and 58% (12/20 in monocytes/ macrophages and 11/20 in microglia) of the most significant eGenes have been previously implicated in AD. *DLG2* encodes a synaptic protein whose expression was previously reported as down-regulated in an AD proteome and transcriptome network 54 and inversely associated with AD Braak stage 30. Genome-wide significant associations of AD risk with *PTPRG* was observed in a family-based GWAS 55 and with *CLNK* in a recent large GWAS for which the evidence was derived almost entirely with a proxy AD phenotype in the UK Biobank 56. *NFXL1* is a novel putative substrate for *BACE1*, an important AD therapeutic target 57. *FCRL5* may interact with the *APOE*E2* allele and also modifies AD age of onset 58. *C4BPA* was shown to be a consistently down-regulated in MCI and AD patients, and the protein encoded by this gene accumulates in Aβ plaques in AD brains 31, 59. Lower levels of the *PAM* have been observed in the brains and CSF of AD patients compared to healthy controls 60 and *MYO1E* is expressed by anti-inflammatory disease associated microglia 32. As a calcium channel protein, *CACNB2* may affect AD risk by altering calcium levels which could cause mitochondrial damage and then induce apoptosis 61, 62.

Likewise, several eGenes of top-ranked ct-eQTLs in microglia that are not established AD loci may have a role in the disease. It was shown that copy number variants (CNVs) near *HNRNPCL1* overlapped the coding portion of the gene in AD cases but not controls 63. A region of epigenetic variation in *ALLC* was associated with AD neuropathology 33. *FAM21B*, a retromer gene in the endosome-to-Golgi retrieval pathway, was associated with AD in a candidate gene study 64. Vacuolar sorting proteins genes in this pathway including *SORL1* have been functionally linked to AD through trafficking of Aβ 65. One study demonstrated that *WNT3* expression in the hippocampus was increased by exercise and alleviated AD-associated memory loss by increasing neurogenesis 34. Expression of *RPL9* is downregulated in severe AD 66 and significantly differs by sex among persons with the *APOE* □4 allele 67. Significant evidence of association with a *TRIM49B* SNP was found in a genome-wide pleiotropy GWAS of AD and major depressive disorder (MDD) 68.

*HLA-DOB*, which is one of the five distinct eGenes (*BTNL3, FAM118A, HLA-DOB, HLA-DRB1*, and *HLA-DRB5)* for significant ct-eQTLs shared between microglia and monocytes/macrophages, and is the target gene for three eSNPs (rs3763355, rs3763354, and rs1183595100) that were evident only in these myeloid cell types.

These eSNPs have similar eQTL p-values in both cell types, but have slightly larger effect sizes in monocytes (Fig. 2). The effect of rs3763355 on expression is in opposite directions in monocytes and microglia which suggests *HLA-DOB* may be acting in different immune capacitates in AD in blood and brain. Though the functions of the genes *BTNL3* and *FAM118A* are unknown, a *BTNL8-BTNL3* deletion has been correlated with TNF and ERK1/AKT pathways, which have an important role in immune regulation inducing inflammation, apoptosis, and proliferation, suggesting the deletion could be correlated to inflammatory disease 69. This suggests that the majority of the shared myeloid cell types genes-the *HLA* genes and possibly *BTNL3*, are all immune-related. Ct-eQTLs involving microglia and monocytes/macrophages had a larger proportion of total intersection, an isolated set interaction and a statistically significant overlap (*P*<1.0E-314), demonstrating a stronger connection than other brain/blood cell types in this study and thus providing further evidence for importance of the immune system in AD.

The proportions of significant ct-eQTLs in NK cells/CD8+T cells, monocytes/ macrophages, and eosinophils are comparable to those observed in reference blood tissue 70, 71. Similarly, significant eQTL distributions in endothelial cells, neurons, and glia are consistent with reference brain tissue 72. The majority of significant blood eQTLs were type I interferon response cells which cross-regulate with pro-inflammatory cytokines that drive pathogenesis of autoimmune diseases including AD and certain heart diseases 73-75and the enrichment of interferon ct-eQTLs in this study could possibly be due to the high proportion of subjects these diseases in the FHS dataset. In contrast, the proportion of significant ct-eQTLs among glial cells is much lower in astrocytes and oligodendrocytes and three-fold higher in microglia than in reference brain tissue 72. Because many AD risk genes are expressed in myeloid cells including microglia 10, the large number of microglia ct-eQTLs is consistent with the high proportion of AD subjects in the ROSMAP dataset.

Several SNPs previously reported to be associated with AD at the GWS level were associated with eGenes that differ from genes ascribed to AD loci and thus may have a role in AD. Karch et al. observed that the expression of *PILRB*, which is involved in immune response and is the activator receptor to its inhibitory counterpart *PILRA*, an established AD gene 76, 77, was highest in microglia 11. *CNN2*, the eGene for eSNP rs4147929 located near the end of *ABCA7*, significantly alters extracellular Aβ levels in human induced pluripotent stem cell-derived neurons and astrocytes 78. Rs4147929 also targeted *HMHA1* which plays several roles in the immune system in an HLA-dependent manner 79. The eSNP/GWAS SNP rs3740688 located in *SPI1* also affects expression of *MYBPC3* and *MADD. MYBPC3* was recently identified as a target gene for eSNPs located in *CELF1* and *MS64A6A* in a study of eQTLs in blood for GWS AD loci 80. *MADD* is expressed in neurons 11, is involved in neuronal cell death in the hippocampus 81, and was shown to be a tau toxicity modulator 82. Although eSNP rs113986870 in *KANSL1* when paired with the nearby eGene *LRRC37A2* was a significant brain eQTL and ct-eQTL, *LRRC37A2* encodes a leucine rich repeat protein that is expressed primarily in testis and has no apparent connection to AD. However, rs113986870 also significantly influenced expression of another gene in this region, *ARL17A*, that was previously linked to progressive supranuclear palsy by analysis of gene expression and methylation 83.

Our study has several noteworthy limitations. The proxy genes individually or collectively may not accurately depict cell-type specific context. In addition, the comparisons of gene expression in blood and brain may yield false results because they are based on independent groups ascertained from a community-based longitudinal study of health (FHS – blood) and multiple sources for studies of AD (ROSMAP – brain) which were not matched for age, sex, ethnicity and other factors which may affect gene expression. Moreover, the FHS cohort contains many elderly but relatively few AD cases, whereas nearly 60% of the ROSMAP participants in this autopsy sample are AD cases. Although the dataset for eQTL analysis in blood was much larger than the dataset derived from brain, the effect sizes associated with many of the eQTLs common to both tissues were similar. Also, findings in brain may reflect post-mortem changes unrelated to disease or cell-type different expression 36. Another limitation of our findings is that some cell types are vastly under-represented compared to others.

Because myeloid cell types are represented in only a small proportion of the total cell populations in brain and blood (generally ∼10%), it is difficult to identify myeloid-specific signals 12. Despite this limitation, many of the most significant and noteworthy results of this study involved monocytes/macrophages and microglia. Finally, some of the findings may not be AD-related because expression was not compared between AD cases and controls.

## Conclusion

Our observation of cell-type specific expression patterns for established and potentially novel AD genes, finding of additional evidence for the role of myeloid cells in AD risk, and discovery of potential novel blood and brain AD biomarkers highlight the importance of cell-type specific analysis. Future studies that compare cell type specific differential gene expression among AD cases and controls using single cell RNA-sequencing or cell count data, as well as functional experiments, are needed to validate and extend our findings.

## Supporting information

Supplemental materials [[supplements/237008_file04.docx]](pending:yes)

Supplemental resources [[supplements/237008_file05.xlsx]](pending:yes)

## Data Availability

This research was conducted in part using data and resources from the FHS of the National Heart, Lung, and Blood Institute (NHLBI) of the US National Institutes of Health (NIH) and the Boston University School of Medicine. The analyses reflect intellectual input and resource development from the FHS investigators participating in the SNP Health Association Resource (SHARe) project and in the Systems Approach to Biomarker Research in Cardiovascular Disease (SABRe) project. The results published here are in whole or in part based on data obtained from the AD Knowledge Portal (https://adknowledgeportal.synapse.org). Study data were provided by the Rush Alzheimers Disease Center, Rush University Medical Center, Chicago. Data collection was supported through funding by NIA grants P30AG10161 (ROS), R01AG15819 (ROSMAP; genomics and RNAseq), R01AG17917 (MAP), R01AG30146, R01AG36042 (5hC methylation, ATACseq), RC2AG036547 (H3K9Ac), R01AG36836 (RNAseq), R01AG48015 (monocyte RNAseq) RF1AG57473 (single nucleus RNAseq), U01AG32984 (genomic and whole exome sequencing), U01AG46152 (ROSMAP AMP-AD, targeted proteomics), U01AG46161(TMT proteomics), U01AG61356 (whole genome sequencing, targeted proteomics, ROSMAP AMP-AD), the Illinois Department of Public Health (ROSMAP), and the Translational Genomics Research Institute (genomic). Additional phenotypic data can be requested at www.radc.rush.edu.

## CONFLICT OF INTEREST

The authors declare no conflicts of interest.

## ACKNOWLEDGEMENTS

This study was supported by NIH grants RF1-AG057519, 2R01-AG048927 U01-AG058654, P30-AG13846, 3U01-AG032984, U01-AG062602 and U19-AG068753. Framingham brain bank data was supported by grants 75N92019D00031 and HHSN2682015000011. Collection of study data provided by the Rush Alzheimer’s Disease Center, Rush University Medical Center, Chicago was supported through funding by NIA grants P30AG10161, R01AG15819, R01AG17917, R01AG30146, R01AG36836, U01AG32984, U01AG46152, U01AG61358, a grant from the Illinois Department of Public Health, and the Translational Genomics Research Institute.

*   Received November 23, 2020.
*   Revision received November 23, 2020.
*   Accepted November 24, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Rao, S. et al. An APOE-independent cis-eSNP on chromosome 19q13.32 influences tau levels and late-onset Alzheimer’s disease risk. Neurobiol Aging 66, 178.e1-178.e8 (2018).
    
    

2.  2.Zou, F. et al. Gene expression levels as endophenotypes in genome-wide association studies of Alzheimer disease. Neurology 74, 480–486 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0b013e3181d07654&link_type=DOI) 

3.  3.Jonkers, I. H. & Wijmenga, C. Context-specific effects of genetic variants associated with autoimmune disease. Hum Mol Genet 26, R185–R192 (2017).
    
    

4.  4.Dobbyn, A. et al. Landscape of conditional eQTL in dorsolateral prefrontal cortex and co-localization with schizophrenia GWAS. Am J Hum Genet 102, 1169–1184 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2018.04.011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29805045&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

5.  5. Daria V Zhernakova et al. Identification of context-dependent expression quantitative trait loci in whole blood. Nat Genet 49, 139–145 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3737&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27918533&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

6.  6.[https://www.hindawi.com/journals/ijad/2012/314185/abs/](https://www.hindawi.com/journals/ijad/2012/314185/abs/).
    
    

7.  7.Raj, T. et al. Polarization of the effects of autoimmune and neurodegenerative risk alleles in leukocytes. Science 344, 519–523 (2014).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNDQvNjE4My81MTkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMS8yNC8yMDIwLjExLjIzLjIwMjM3MDA4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

8.  8.Battle, A., Brown, C. D., Engelhardt, B. E. & Montgomery, S. B. Genetic effects on gene expression across human tissues. Nature 550, 204–213 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature24277&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29022597&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000412829500039&link_type=ISI) 

9.  9.Efthymiou, A. G. & Goate, A. M. Late onset Alzheimer’s disease genetics implicates microglial pathways in disease risk. Mol Neurodegen 12 (2017).
    
    

10. 10.Huang, K. et al. A common haplotype lowers PU.1 expression in myeloid cells and delays onset of Alzheimer’s disease. Nat Neurosci 20, 1052–1061 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn.4587&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28628103&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

11. 11.Karch, C. M., Ezerskiy, L. A., Bertelsen, S. & Goate, A. M. Alzheimer’s disease risk polymorphisms regulate gene expression in the ZCWPW1 and the CELF1 loci. PLoS One 11 (2016).
    
    

12. 12.Novikova, G. et al. Integration of Alzheimer’s disease genetics and myeloid genomics reveals novel disease risk mechanisms. bioRxiv, 694281 (2019).
    
    

13. 13.Zhang, X. et al. Identification of common genetic variants controlling transcript isoform variation in human whole blood. Nat Genet 47, 345–352 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3220&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25685889&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

14. 14.Bennett, D. A., Schneider, J. A., Arvanitakis, Z. & Wilson, R. S. Overview and findings from the religious orders study. Curr Alzheimer Res 9, 628–645 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2174/156720512801322573&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22471860&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308864500003&link_type=ISI) 

15. 15.[http://www.eurekaselect.com/99959/article](http://www.eurekaselect.com/99959/article).
    
    

16. 16.[https://www.synapse.org/](https://www.synapse.org/).
    
    

17. 17.Leek, J. T., Johnson, W. E., Parker, H. S., Jaffe, A. E. & Storey, J. D. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics 28, 882–883 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bts034&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22257669&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000301972900020&link_type=ISI) 

18. 18.Therneau, T. The Lmekin Function. (2018).
    
    

19. 19.[http://www.R-project.org/](http://www.R-project.org/).
    
    

20. 20.Penney, J., Ralvenius, W. T. & Tsai, L. Modeling Alzheimer’s disease with iPSC-derived brain cells. Mol Psychiatry 25, 148–167 (2020).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

21. 21.McKenzie, A. T. et al. Brain cell type specific gene expression and co-expression network architectures. Sci Rep 8, 8868 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-018-27293-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29892006&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

22. 22.Ren, Y. et al. TMEM106B haplotypes have distinct gene expression patterns in aged brain. Mol Neurodegen 13 (2018).
    
    

23. 23.Mi, H. et al. PANTHER version 11: expanded annotation data from gene ontology and reactome pathways, and data analysis tool enhancements. Nucleic Acids Res 45, D183–D189 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkw1138&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27899595&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000396575500028&link_type=ISI) 

24. 24.Giambartolomei, C. et al. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLOS Genet 10, e1004383 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1004383&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24830394&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

25. 25.Kunkle, B. W. et al. Genetic meta-analysis of diagnosed Alzheimer’s disease identifies new risk loci and implicates Aβ, tau, immunity and lipid processing. Nat Genet 51, 414–430 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0358-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30820047&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

26. 26.Wu, Y. et al. Colocalization of GWAS and eQTL signals at loci with multiple signals identifies additional candidate genes for body fat distribution. Hum Mol Genet 28, 4161–4172 (2019).
    
    

27. 27.Pruim, R. J. et al. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 26, 2336–2337 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq419&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20634204&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281714100054&link_type=ISI) 

28. 28.Lambert, J. C. et al. Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat Genet. 45, 1452–1458 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2802&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24162737&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

29. 29.Jun, G. et al. A novel Alzheimer disease locus located near the gene encoding tau protein. Mol. Psychiatry 21, 108–117 (2016).
    
    

30. 30.Hondius, D. C. et al. Profiling the human hippocampal proteome at all pathologic stages of Alzheimer’s disease. Alzheimers Dement 12, 654–668 (2016).
    
    

31. 31.Sjöberg, A. P., Trouw, L. A. & Blom, A. M. Complement activation and inhibition: a delicate balance. Trends Immunol. 30, 83–90 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.it.2008.11.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19144569&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264279300005&link_type=ISI) 

32. 32.Rangaraju, S. et al. Identification and therapeutic modulation of a pro-inflammatory subset of disease-associated-microglia in Alzheimer’s disease. Mol Neurodegen 13, 24 (2018).
    
    

33. 33.Smith, A. R. et al. Parallel profiling of DNA methylation and hydroxymethylation highlights neuropathology-associated epigenetic variation in Alzheimer’s disease. Clin Epigenetics 11, 52 (2019).
    
    

34. 34.Kim, D., Jung, S., Kim, K. & Kim, C. Treadmill exercise ameliorates Alzheimer disease-associated memory loss through the Wnt signaling pathway in the streptozotocin-induced diabetic rats. J Exerc Rehabil 12, 276–283 (2016).
    
    

35. 35.Xu, X., Wells, A. B., O’Brien, D. R., Nehorai, A. & Dougherty, J. D. Cell type-specific expression analysis to identify putative cellular mechanisms for neurogenetic disorders. J. Neurosci. 34, 1420–1431 (2014).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjk6IjM0LzQvMTQyMCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzExLzI0LzIwMjAuMTEuMjMuMjAyMzcwMDguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

36. 36.Blair, J. A. et al. Individual case analysis of postmortem interval time on brain tissue preservation. PLoS One 11 (2016).
    
    

37. 37.Le Page, A. et al. Role of the peripheral innate immune system in the development of Alzheimer’s disease. Exp. Gerontol. 107, 59–66 (2018).
    
    

38. 38.Goldeck, D., Witkowski, J. M., Fülop, T. & Pawelec, G. Peripheral immune signatures in Alzheimer disease. Curr Alzheimer Res 13, 739–749 (2016).
    
    

39. 39.Harris, S. A. & Harris, E. A. Herpes simplex virus type 1 and other pathogens are key causative factors in sporadic Alzheimer’s disease. J Alzheimers Dis. 48, 319–353 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3233/JAD-142853&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26401998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

40. 40.Obulesu, M. & Lakshmi, M. J. Apoptosis in Alzheimer’s disease: an understanding of the physiology, pathology and therapeutic avenues. Neurochem. Res. 39, 2301–2312 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11064-014-1454-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25322820&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

41. 41.Fu, H. et al. A tau homeostasis signature is linked with the cellular and regional vulnerability of excitatory neurons to tau pathology. Nat. Neurosci. 22, 47–56 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41593-018-0298-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30559469&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

42. 42.Crehan, H. et al. Complement receptor 1 (CR1) and Alzheimer’s disease. Immunobiology 217, 244–250 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.imbio.2011.07.017&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21840620&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

43. 43.L, R. et al. Differential expression of ribosomal genes in brain and blood of Alzheimer’s disease patients. Curr Alzheimer Res 12, 984–989 (2015).
    
    

44. 44.Jun, G. R. et al. Transethnic genome-wide scan identifies novel Alzheimer’s disease loci. Alzheimers Dement 13, 727–738 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jalz.2016.12.012&link_type=DOI) 

45. 45.Desikan, R. S. et al. Polygenic overlap between C-reactive protein, plasma lipids, and Alzheimer disease. Circulation 131, 2061–2069 (2015).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjExOiIxMzEvMjMvMjA2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzExLzI0LzIwMjAuMTEuMjMuMjAyMzcwMDguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

46. 46.Chang, J. & Chang, N. WWOX dysfunction induces sequential aggregation of TRAPPC6AΔ, TIAF1, tau and amyloid β, and causes apoptosis. Cell Death Discov 1, 15003 (2015).
    
    

47. 47.Liu, C. et al. WWOX phosphorylation, signaling, and role in neurodegeneration. Front Neurosci 12, 563 (2018).
    
    

48. 48.Steele, N. Z. R. et al. Fine-mapping of the human leukocyte antigen locus as a risk factor for Alzheimer disease: A case–control study. PLOS Medicine 14, e1002272 (2017).
    
    

49. 49.Xiang, Z., Haroutunian, V., Ho, L., Purohit, D. & Pasinetti, G. M. Microglia activation in the brain as inflammatory biomarker of Alzheimer’s disease neuropathology and clinical dementia. Dis. Markers 22, 95–102 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1155/2006/276239&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16410654&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249581800008&link_type=ISI) 

50. 50.Walker, D. G. & Lue, L. Immune phenotypes of microglia in human neurodegenerative disease: challenges to detecting microglial polarization in human brains. Alzheimers Res Ther 7, 56 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13195-015-0139-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26286145&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

51. 51.Neill, D. et al. Risk for Alzheimer’s disease in older late-onset cases is associated with HLA-DRB1*03. Neurosci. Lett. 275, 137–140 (1999).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10568518&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

52. 52.Lehmann, D. J. et al. HLA class I, I. & III genes in confirmed late-onset Alzheimer’s disease. Neurobiol. Aging 22, 71–77 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0197-4580(00)00180-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11164278&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000166728900010&link_type=ISI) 

53. 53.Lambert, J. et al. Implication of the immune system in Alzheimer’s disease: evidence from genome-wide pathway analysis. J. Alzheimers Dis. 20, 1107–1118 (2010).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20413860&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000279539500015&link_type=ISI) 

54. 54.Hallock, P. & Thomas, M. A. Integrating the Alzheimer’s disease proteome and transcriptome: a comprehensive network model of a complex disease. OMICS 16, 37–49 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/omi.2011.0054&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22321014&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

55. 55.Herold, C. et al. Family-based association analyses of imputed genotypes reveal genome-wide significant association of Alzheimer’s disease with OSBPL6, PTPRG, and PDCL3. Mol. Psychiatry 21, 1608–1612 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/mp.2015.218&link_type=DOI) 

56. 56.Jansen, I. E. et al. Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer’s disease risk. Nat Genet 51, 404–413 (2019).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

57. 57.Johnson, J. L., Chambers, E. & Jayasundera, K. Application of a bioinformatics-based approach to identify novel putative in vivo BACE1 substrates. Biomed Eng Comput Biol 5, 1–15 (2013).
    
    

58. 58.Vélez, J. I. et al. APOE*E2 allele delays age of onset in PSEN1 E280A Alzheimer’s disease. Molecular Psychiatry 21, 916–924 (2016).
    
    

59. 59.Song, F. et al. Plasma protein profiling of Mild Cognitive Impairment and Alzheimer’s disease using iTRAQ quantitative proteomics. Proteome Sci 12, 5 (2014).
    
    

60. 60.Wand, G. S. et al. Alzheimer’s disease: low levels of peptide alpha-amidation activity in brain and CSF. Neurology 37, 1057–1061 (1987).
    
    

61. 61.Liang, X. et al. Genomic convergence to identify candidate genes for Alzheimer disease on chromosome 10. Hum Mutat 30, 463–471 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/humu.20953&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19241460&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264011700024&link_type=ISI) 

62. 62.Pang, C. et al. Identification and analysis of Alzheimer’s candidate genes by an amplitude deviation algorithm. J Alzheimers Dis Parkinsonism 9, 460 (2019).
    
    

63. 63.Swaminathan, S. et al. Genomic Copy number analysis in Alzheimer’s disease and mild cognitive impairment: an ADNI study. Int J Alzheimers Dis 2011, 729478 (2011).
    
    

64. 64.Vardarajan, B. N. et al. Candidate gene study in the endosome-to-Golgi retrieval pathway reveals association of retromer genes with Alzheimer’s disease. Alzheimer Dement 6, S194 (2010).
    
    

65. 65.Rogaeva, E. et al. The neuronal sortilin-related receptor SORL1 is genetically associated with Alzheimer’s Disease. Nat Genet 39, 168–177 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng1943&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17220890&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000244063900014&link_type=ISI) 

66. 66.Kong, W. et al. Independent component analysis of Alzheimer’s DNA microarray gene expression data. Mol Neurodegener 4, 5 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1750-1326-4-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19173745&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

67. 67.Hsu, M., Dedhia, M., Crusio, W. E. & Delprato, A. Sex differences in gene expression patterns associated with the APOE4 allele. F1000Res 8, 387 (2019).
    
    

68. 68.Lutz, M. W., Sprague, D., Barrera, J. & Chiba-Falek, O. Shared genetic etiology underlying Alzheimer’s disease and major depressive disorder. Transl Psychiatry 10 (2020).
    
    

69. 69.Guo, Y. & Wang, A.Y. Novel immune check-point regulators in tolerance maintenance. Front Immunol. 6, 421 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fimmu.2015.00421&link_type=DOI) 

70. 70.[https://www.miltenyibiotec.com/US-en/resources/macs-handbook/human-cells-and-organs/human-cell-sources/blood-human.html#gref](https://www.miltenyibiotec.com/US-en/resources/macs-handbook/human-cells-and-organs/human-cell-sources/blood-human.html#gref).
    
    

71. 71.[https://www.stemcell.com/media/files/wallchart/WA10006-Frequencies\_Cell\_Types\_Human\_Peripheral\_Blood.pdf](https://www.stemcell.com/media/files/wallchart/WA10006-Frequencies\_Cell\_Types_Human_Peripheral_Blood.pdf).
    
    

72. 72.von Bartheld, C. S., Bahney, J. & Herculano-Houzel, S. The search for true numbers of neurons and glial cells in the human brain: A review of 150 years of cell counting. J Comp Neurol 524, 3865–3895 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cne.24040&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27187682&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

73. 73.Dinkins, M. B. et al. The 5XFAD mouse model of Alzheimer’s disease exhibits an age-dependent increase in anti-ceramide IgG and exogenous administration of ceramide further increases anti-ceramide titers and amyloid plaque burden. J Alzheimers Dis 46, 55–61 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3233/JAD-150088&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25720409&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

74. 74.Abou-Raya, A. & Abou-Raya, S. Inflammation: A pivotal link between autoimmune diseases and atherosclerosis. Autoimmunity Reviews 5, 331–337 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.autrev.2005.12.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16782558&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000238925700006&link_type=ISI) 

75. 75.Ivashkiv, L. B. & Donlin, L. T. Regulation of type I interferon responses. Nat Rev Immunol 14, 36–49 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nri3581&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24362405&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

76. 76.Patel, T. et al. Whole-exome sequencing of the BDR cohort: evidence to support the role of the PILRA gene in Alzheimer’s disease. Neuropathol. Appl. Neurobiol. (2017).
    
    

77. 77.Logue, M. W. et al. Search for age-related macular degeneration risk variants in Alzheimer disease genes and pathways. Neurobiol. Aging 35, 1510.e7-18 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neurobiolaging.2013.12.007&link_type=DOI) 

78. 78.Se, S. et al. Candidate-based screening via gene modulation in human neurons and astrocytes implicates FERMT2 in Aβ and TAU proteostasis. Hum Mol Genet 28, 718–735 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddy376&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30371777&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

79. 79.de Kreuk, B. et al. The Human Minor Histocompatibility Antigen1 Is a RhoGAP. PLoS One 8 (2013).
    
    

80. 80.Katsumata, Y., Nelson, P. T., Estus, S. & Fardo, D. W. Translating Alzheimer’s disease–associated polymorphisms into functional candidates: a survey of IGAP genes and SNPs. Neurobiol Aging 74, 135–146 (2019).
    
    

81. 81.Del Villar, K. & Miller, C. A. Down-regulation of DENN/MADD, a TNF receptor binding protein, correlates with neuronal cell death in Alzheimer’s disease brain and hippocampal neurons. Proc. Natl. Acad. Sci. U. S. A. 101, 4210–4215 (2004).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTAxLzEyLzQyMTAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMS8yNC8yMDIwLjExLjIzLjIwMjM3MDA4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

82. 82.Dourlen, P. et al. Functional screening of Alzheimer risk loci identifies PTK2B as an in vivo modulator and early marker of Tau pathology. Molecular Psychiatry 22, 874–883 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/mp.2016.59&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27113998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F24%2F2020.11.23.20237008.atom) 

83. 83.Allen, M. et al. Gene expression, methylation and neuropathology correlations at progressive supranuclear palsy risk loci. Acta Neuropathol 132, 197–211 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00401-016-1576-7&link_type=DOI)

 [1]: /embed/graphic-1.gif
 [2]: /embed/graphic-2.gif
 [3]: /embed/graphic-3.gif
 [4]: /embed/graphic-4.gif