Genome-Wide Association and Mechanistic Studies Indicate That Immune Response Contributes to Alzheimer’s Disease Development

Alzheimer’s disease (AD) is the most common cause of dementia. Although genome-wide association study (GWAS) have reported hundreds of single-nucleotide polymorphisms (SNPs) and genes linked to AD, the mechanisms about how these SNPs modulate the development of AD remain largely unknown. In this study, we performed GWAS for three traits in cerebrospinal fluid (CSF) and one clinical trait in the Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort. Our analysis identified five most significant AD related SNPs (FDR < 0.05) within or proximal to APOE, APOC1, and TOMM40. One of the SNPs was co-inherited with APOE allele 4, which is the most important genetic risk factor for AD. Three of the five SNPs were located in promoter or enhancer regions, and transcription factor (TF) binding affinity calculations showed dramatic changes (| Log2FC| > 2) of three TFs (PLAG1, RREB1, and ZBTB33) for two motifs containing SNPs rs2075650 and rs157580. In addition, our GWAS showed that both rs2075650 and rs157580 were significantly associated with the poliovirus receptor-related 2 (PVRL2) gene (FDR < 0.25), which is involved in spreading of herpes simplex virus (HSV). The altered regulation of PVRL2 may increase the susceptibility AD patients to HSV and other virus infections of the brain. Our work suggests that AD is a type of immune disorder driven by viral or microbial infections of the brain during aging.

Introduction

Alzheimer’s disease (AD) is the most common neurodegenerative disorder. It is characterized by progressive memory loss and cognitive decline, cerebral accumulation of amyloid-β peptide (Aβ) in senile plaques and hyper-phosphorylated tau in neurofibrillary tangles (NFT) (Price and Sisodia, 1998; Mukaetova-Ladinska et al., 2015). Since AD is a complex and multifactorial disease, large datasets with multiple data types have been critical to identify its genetic risk factors (Harold et al., 2009). For several decades, only the allele 4 of Apolipoprotein E (APOE), which is present in about half of late-onset AD (LOAD) patients, has been convincingly demonstrated to affect risk for LOAD (Bertram et al., 2010).

Studies have shown that the levels of amyloid-β 1-42 peptide (Aβ₄₂), total tau (T-tau), and tau phosphorylated at threonine 181 (p-tau) in cerebrospinal fluid (CSF) samples can be used as AD diagnostic biomarkers (Shaw et al., 2009; Hampel et al., 2010). AD patients show lower levels of CSF Aβ₄₂ (Noguchi et al., 2005), which is negatively correlated with Aβ plaque counts in brain samples (Ikonomovic et al., 2008). The CSF levels of T-tau and p-tau are increased in AD patients (Price and Morris, 1999). Increased CSF T-tau levels are also found in stroke (Hesse et al., 2001) and traumatic brain injury (Ost et al., 2006). However, elevated CSF p-tau levels appear to be specific to AD (Buerger et al., 2006). The 13-item version of the Alzheimer’s Disease Assessment Scale-Cognitive subscale (ADAS13) was developed to measure memory and cognition for patients with mild to moderate AD (Podhorna et al., 2016). It’s one of the most frequently applied tests in experimental studies and clinical trials for new drugs and other interventions. A normal ADAS13 score for a person who does not have AD is 5 (Graham et al., 2004), while 31.2 is the average score for those who have been diagnosed with AD or mild cognitive impairment (The Alzheimer’s Disease Neuroimaging Initiative et al., 2015).

Most genetic association studies analyze at a single marker level or focus on detecting risk factors for AD, but ignore the mechanisms and functions associated with what they find. In our study, we conducted quantitative trait locus (QTL) analysis of Aβ₄₂, T-tau/Aβ₄₂ ratio, p-tau/Aβ₄₂ ratio and ADAS13 as quantitative traits on SNPs to identify significantly AD-associated SNPs. We then performed expression quantitative trait loci (eQTL) analysis for these AD-associated SNPs to locate the related genes. After that, we computed the linkage disequilibrium (LD) pattern, allele distributions and transcription factor motifs binding affinity for SNPs in regulatory regions (promoter regions or enhancer regions) to study how they modulate the target genes (Figure 1). Overall, our analysis may contribute to understand the mechanism and etiology of AD.

FIGURE 1

FIGURE 1. The flow chart of the analysis conducted in this study.

Materials and Methods

Data used in this study were obtained from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (adni.loni.usc.edu). ADNI was launched in 2003 by the National Institute on Aging (NIA), the National Institute of Biomedical Imaging and Bioengineering (NIBIB), the Food and Drug Administration (FDA), and by private pharmaceutical companies and non-profit organizations, as a public-private partnership. The principal investigator of ADNI is Michael W. Weiner, MD. The primary goal of ADNI has been to test whether serial magnetic resonance imaging (MRI), positron emission tomography (PET), biological markers, and clinical and neuropsychological assessment can be combined together to measure the progression of AD.

The first phase of ADNI (ADNI-1) launched in 2004 included 400 subjects diagnosed with late mild cognitive impairment (LMCI), 200 subjects with AD, and 200 elderly cognitively normal (CN) subjects. ADNI was extended in 2009 during the ADNI-GO phase, which assessed the existing ADNI-1 cohort along with 200 new participants with early mild cognitive impairment (EMCI). In 2011, assessing participants from ADNI-1/ADNI-GO phases, ADNI-2 began with new subject groups: 150 CN, 100 EMCI, 150 LMCI and 150 AD.

We were granted permission to obtain data from the ADNI cohort (www.adni-info.org) for performing the analysis described in this paper.

Subjects

We analyzed the CSF Aβ₄₂, T-tau and p-tau levels, ADAS13 scores, quality-controlled SNP data and gene expression data for 812 ADNI subjects, including 281 CN, 235 EMCI, 249 LMCI and 47 AD cases.

Genotyping Data

The SNP data of ADNI-1, ADNI-GO, and ADNI-2 cohorts were collected from either the Illumina 2.5-M array or the Illumina OmniQuad array (Saykin et al., 2010). The SNPs shown in both arrays were used for the following analysis.

Quality control (QC) analysis was conducted by using R package snpStats (Clayton, 2012) in R software (R Core Team, 2013). In the QC, we excluded any SNPs that did not meet any of the following criteria: (1) SNPs on chromosome 1-22; (2) call rate per SNP > 95%; (3) minor allele frequency (MAF) > 5%; (4) Hardy-Weinberg equilibrium (HWE) test of p-value > 10^-6 (absolute value of z-score < 4.753424). After QC analysis, 575353 SNPs remained for the subsequent analysis.

CSF Biomarkers and ADAS13

The CSF levels of Aβ₄₂, T-tau and p-tau were determined using the fully automated Roche Elecsys immunoassay platform (Seibyl et al., 2017). The ADAS13 was conducted by a certified ADAS13 rater. ADAS13 scores were automatically calculated on the electronic case report form based on item level data entered.

Gene Expression

Gene expression profiles of peripheral blood samples from ADNI participants were performed at Bristol-Myers Squibb (BMS) laboratories. The Affymetrix Human Genome U219 Array (www.affymetrix.com) was used for expression profiling, which contains 530467 probes for 49293 transcripts. Raw expression values obtained directly from CEL files were pre-processed using the RMA (Robust Multi-chip Average) normalization method (Vawter et al., 2004).

We also used the dataset GSE28146 (Blalock et al., 2011) in Gene Expression Omnibus [GEO (Barrett et al., 2012)]. The gene expression data was collected from laser-captured hippocampus tissue gray matter from formalin fixed, paraffin embedded specimens by Affymetrix Human Genome U133 Plus 2.0 Array. The dataset includes eight controls and seven cases of severe AD.

Regulatory Information for SNP Location

The regulatory region information (promoter, enhancer and TF binding regions) of the identified AD-associated SNPs was obtained from Genome Browser (Kent et al., 2002) and HaploReg (Ward and Kellis, 2012). We used TF binding motifs data in Motif browser (Kheradpour and Kellis, 2014) and R package MotifDb (Shannon, 2014) with motifs in JASPAR (Khan et al., 2018), SwissRegulon (Pachkov et al., 2012) and other databases, for the computation of binding affinity.

Analysis

To determine the SNPs associated with Aβ₄₂ level, T-tau/Aβ₄₂ ratio, p-tau/Aβ₄₂ ratio in CSF and ADAS13 score, we performed QTL analysis using the R package MatrixEQTL (Shabalin, 2012). In this analysis, age, gender (1 for male and 2 for female) and diagnosis (1 for CN, 2 for EMCI, 3 for LMCI and 4 for AD) at baseline were considered as covariates. Manhattan plots of QTL results were generated using the R package qqman (Turner, 2014). Following analysis was restricted to the associated SNPs on the chromosome 19 (for APOE on chromosome 19) with the false discovery rate (FDR) < 0.05. eQTL analysis was also conducted using the R package MatrixEQTL for the filtered SNPs, with the same covariates as QTL. The cis-eQTL results (local, distance < 1 Mb) with FDR < 0.25 were considered as significant. The LD pattern was plotted using the R package LDheatmap (Shin et al., 2006).

TF motifs binding affinity were calculated using the R package PWMEnrich for the SNPs in the regulatory regions. We used the regions located +/- 20 base pair (bp) of the SNPs as the DNA strings to perform the computation of binding affinity. To determine the changes in binding affinity, we calculated the log2 fold change (Log2FC) of binding affinity for major alleles against the binding affinity for minor alleles. TF motifs with binding affinity of both sequences with major allele and minor allele less than 1 (unlikely to bind on the sequences) were excluded for the following analysis. Remaining TF motifs and related SNPs with absolute value of Log2FC greater than 2 (|Log2FC| > 2) were used for the functional and mechanism analysis.

Gene differential expression between CN group and AD group was conducted by the R package limma (Ji et al., 2014; Ritchie et al., 2015). We checked the PVRL2 gene expression for both ADNI and dataset GSE28146 from GEO.

Results

Our results of QTL analysis for Aβ₄₂, T-tau/Aβ₄₂ ratio, p-tau/Aβ₄₂ ratio and ADAS13 (Supplementary Tables S1–S4) were illustrated in the Manhattan plots (Figure 2), showing that almost all SNPs significantly associated with AD (p-value < 10^-7) are located at the chromosome 19 which contains APOE gene. Table 1 shows the SNPs with FDR < 0.05 from the QTL analysis of Aβ₄₂, T-tau/Aβ₄₂ ratio, p-tau/Aβ₄₂ ratio and ADAS13, respectively. A total of five SNPs were identified according to the threshold FDR < 0.05, and they were also reported by previous studies (Li et al., 2017). Detected SNPs rs2075650, rs157580 and rs157582 are located within gene TOMM40. SNP rs769449 locates within gene APOE and rs4420638 is proximal to the downstream of gene APOC1. SNP rs4420638 is the most significant AD-associated SNP in our QTL analysis for Aβ₄₂ (FDR = 5.47E-15), T-tau/Aβ₄₂ ratio (FDR = 4.28E-23), and p-tau/Aβ₄₂ ratio (FDR = 1.78E-21; Figure 2 and Table 1). It is also the second most significant AD-associated SNP in our QTL analysis for ADAS13 (FDR = 8.73E-06). These imply that rs4420638 is the most relevant SNP to AD in ADNI according to our QTL analysis. The significant SNPs identified based on the T-tau/Aβ₄₂ ratio and p-tau/Aβ₄₂ ratio are identical due to the number and the order of SNPs (Figure 2 and Table 1). SNPs rs4420638, rs769449, and rs2075650 were identified to be significant for all of the four traits and ranked in the top three (Figure 2 and Table 1).

FIGURE 2

FIGURE 2. The Manhattan plots of the observed –log₁₀ (p-value) for the results of QTL analysis for (A) T-tau/Aβ₄₂ ratio, (B) p-tau/Aβ₄₂ ratio, (C) Aβ₄₂ concentration, and (D) ADAS13 scores. SNPs with FDR < 0.05 were labeled on each plot. Redline stood for p-value = 10^-7.

TABLE 1

TABLE 1. Results of QTL analysis.

We then performed eQTL analysis for the five identified SNPs (Supplementary Table S5). Table 2 shows the cis-eQTL results with FDR < 0.25. In order to include more SNPs and genes for the following mechanistic analysis, we set FDR < 0.25 as the threshold instead of the more classic cutoff FDR < 0.05 for our eQTL results. For we want to find the mechanism of AD, here we just focused on the cis-eQTL results and did not consider the trans-eQTL results. SNPs rs769449, rs2075650 and rs157580 were associated with PVRL2 gene (FDR = 0.19035, 0.0367, and 0.08866, respectively), while rs4420638 was associated with SYMPK gene (FDR = 0.22, Table 2). For three SNPs showed significant association with PVRL2 gene, we studied its function in the following. In addition to PVRL2 gene, rs2075650 was also associated with HIF3A gene (FDR = 0.19035, Table 2). There were no genes significantly associated with rs157582 according to our cis-eQTL results.

TABLE 2

TABLE 2. Results of eQTL analysis for the five SNPs closely related to AD.

We calculated allele distributions and minor allele frequencies (MAF) of the five identified SNPs for the four diagnosis groups (Table 3) to study their different features among different groups. Except for rs157580, MAFs were increased along with the progression of AD. For rs4420638, rs769449, rs2075650, and rs157582, the percentage of subjects carrying two minor alleles was greater in the AD group than that in the other three diagnosis groups. These genetic features of AD group highlight the contribution of the genotype to the development of AD.

TABLE 3

TABLE 3. Alleles distribution of the five AD-associated SNPs in four diagnostic groups.

Among these five identified SNPs, rs2075650 shows strong LD with rs769449 (D′ = 0.91, r² = 0.69; Figure 3) and rs157582 (D′ = 1.00, r² = 0.61; Figure 3). The strong LD between rs2075650 and rs157582 may be due to that both of them locate within gene TOMM40 closely. Moreover, rs4420638 shows strong LD with rs769449 (D′ = 0.97, r² = 0.56; Figure 3). Besides the LD pattern for these five SNPs, we checked that SNP rs4420638 was in strong LD with rs429358 (D′ = 0.94, r² = 0.65; by HaploReg), which contributes to define the allele types of APOE gene. We explored the relationship between rs4420638 alleles and the APOE allele 4 copy numbers (Figure 4). Our results showed that minor allele G of rs4420638 was highly correlated with APOE allele 4 copy numbers. Most subjects carrying one or two G alleles have one or two copies of APOE allele 4 among all the four diagnosis groups, respectively. This situation was much more obvious for AD group. All AD patients with two G alleles had two copies of APOE allele 4. All but one (95.65%) AD subjects carrying one allele G had just one APOE allele 4.

FIGURE 3

FIGURE 3. The LD pattern for the five significantly associated SNPs.

FIGURE 4

FIGURE 4. The relationship between APOE allele 4 copy numbers and rs4420638 alleles in four diagnostic groups.

Further analysis shows that rs769449, rs2075650, and rs157580 are located in promoter and enhancer regions for multiple types of brain tissues and these SNPs are also located in TF binding regions (annotations from Genome Browser and HaploReg). We calculated the corresponding TF motifs binding affinity for these three SNPs (Supplementary Table S6). SNP rs2075650 and rs157580 significantly altered the binding affinity (|Log2FC| > 2) of four TF motifs (Figure 5). SNP rs2075650 altered the binding affinity of TFs PLAG1 (Log2FC = -7.17919) and RREB1 (Log2FC = -2.517812 for motif MA0073.1 and Log2FC = -3.701917 for motif RREB1.SwissRegulon) while rs157580 only altered the binding affinity of TF ZBTB33 (Log2FC = -7.642652). All of the four TF motifs had higher binding affinity scores to the sequences containing the minor allele of the related SNPs, indicating that the regulatory function of these TFs may be enhanced for the subjects carrying minor allele of rs2075650 and rs157580. Through altering the regulatory function of these TFs, these SNPs can affect their target genes.

FIGURE 5

FIGURE 5. The transcription factor binding affinity for AD related SNPs in regulatory regions. Only binding affinity of TF motifs that are largely altered (|Log2FC| > 2) by the SNP genotypes are listed. PWM, position weight matrix. Log2FC, Log2 fold change for the binding affinity of sequence with major allele against the binding affinity of sequence with minor allele.

Discussion

APOE has three different alleles (allele 2, 3, and 4) that are defined by two SNPs: rs429358 (T/C, T as major allele) and rs7412 (C/T, C as major allele). APOE allele 2 [(T;T) for (rs429358 and rs7412)] is relatively rare and considered to be protective against AD. The most common allele is the APOE allele 3 (T;C) which is believed to have no effect on AD (Farrer, 1997). APOE allele 4 (C;C) increases the risk for developing AD (Liu et al., 2013). Due to the diploid nature of the human genome, a normal person can have zero, one, or two copies of APOE allele 4. Carrying one copy of APOE allele 4 can increase the risk of AD by 2 to 3 times, while carrying two such copies have a 25-fold increased risk for developing AD compared to people with two APOE allele 3 (Michaelson, 2014). SNP rs4420638 is the most significant AD-associated SNP identified in our study (Figure 2 and Table 1). Its minor allele G occurs almost simultaneously with APOE allele 4, especially in AD patients (Figure 4). More importantly, increased MAF for this SNP was observed in the AD patients (Table 3). These results suggest that rs4420638 (allele G) is co-inherited with APOE allele 4, which is supported by the previous studies (Li et al., 2008).

Pleomorphic adenoma gene 1 (PLAG1), encoding a zinc finger protein with 2 putative nuclear localization signals, has been shown to be consistently rearranged in pleomorphic adenomas of the salivary glands. As a transcription factor, its activation leads to the up-regulation of target genes (Van Dyck et al., 2007). Ras-Responsive Element-Binding Protein 1 is a zinc finger transcription factor encoded by RREB1 gene. It binds to RAS-responsive elements (RREs) of gene promoters and potentiates the transcriptional activity of Neurogenic differentiation 1 gene (NEUROD1) and other genes (Melani et al., 2008). SNP rs2075650 altered the binding affinity for some motifs of PLAG1 and RREB1, where minor allele G increased the binding affinity of these motifs for both TFs (Figure 5). Moreover, AD patients had the highest MAF for this SNP among all the four diagnosis groups (Table 3), implying increased function and regulation of these two TFs and increased activation of their target genes in the subjects of AD group. Zinc finger and BTB domain-containing protein 33 (ZBTB33) encodes a transcriptional regulator with bimodal DNA-binding specificity. It recruits the N-CoR repressor complex to promote histone deacetylation and the formation of repressive chromatin structures in target gene promoters. In addition, it may contribute to the repression of target genes (Daniel and Reynolds, 1999). By our analysis, ZBTB33 is more likely to bind to the string with minor allele G of rs157580 (Figure 5), compared with the string containing major allele A. The MAF of rs157580 in AD group was less than that in the CN group (Table 3), suggesting that the regulatory function of ZBTB33 in AD patients is weaker than that in the healthy people. Our eQTL analysis indicated that both rs2075650 and rs157580 were significantly associated with the PVRL2 gene (Table 2). PVRL2 appears to be the common target gene of the transcription factor PLAG1, RREB1, and ZBTB33. The altered binding affinity of these TFs in AD patients may potentiate the transcriptional activity of PVRL2. More importantly, PVRL2 was found to be upregulated in AD patients compared to CN in ADNI and another dataset (Table 4).

TABLE 4

TABLE 4. PVRL2 gene differential expression for AD group against CN group in ADNI and GSE28146.

Poliovirus receptor-related 2 (PVRL2) gene, also known as nectin cell adhesion molecule 2 (NECTIN2) or herpesvirus entry mediator B (HVEB), is a human plasma membrane glycoprotein (Lopez et al., 1998). PVRL2 gene is located in close vicinity of the APOE locus on the chromosome 19. It is expressed in multiple types of cell and tissues, including neurons. It belongs to the plasma membrane components of adherens junctions (Takai, 2003). More importantly, it mediates the entry of herpes simplex virus (HSV) (Warner et al., 1998). Dysregulation for PVRL2 may have impact on the susceptibility of individuals to HSV infection of brain by affecting virus entry to cells and intercellular virus spreading (Spear, 2004). There is evidence showing that HSV infection in AD brains was observed (Itzhaki and Wozniak, 2008). A recent study reported that Aβ peptide might be involved in resistance to microbial infection in mouse and worm models of AD and acted as a defense molecule of the innate immunity (Kumar et al., 2016), which is compatible with the viral association with AD etiology and pathology. The accumulation of Aβ plaque deposits may be a consequence of the over-production of Aβ peptide during viral infection of the brain (Porcellini et al., 2010). Additionally, some previous studies showed that HSV could contribute to the development of Aβ plaques directly (Piacentini et al., 2011). What’s more, reactivated HSV1 can directly induce inflammatory damage, which may lead to an increase on the formation of Aβ and tau pathology (Wozniak et al., 2007). Last but the least, a new study of late-onset AD-associated virome found the evidence that links the activity of specific viral species (includes HSV1) with molecular, genetic, clinical and neuropathological aspects of AD (Readhead et al., 2018).

In summary, we identified five SNPs that were highly associated with AD by QTL analysis for multiple traits. Among them, the allele G of rs4420638 is co-inherited with the APOE allele 4. In addition, our eQTL analysis indicated that both rs2075650 and rs157580 were significantly associated with the PVRL2 gene. The transcription factors PLAG1, RREB1 and ZBTB33 had higher binding affinities to the motifs containing minor allele of the rs2075650 and rs157580, indicating an enhanced regulatory function of these TFs in the subjects with such minor alleles. The altered binding affinity of these TFs to these SNP regions increased the expression of PVRL2 in AD patients as compared to CN. Increased expression of PVRL2 may also increase the patients’ susceptibility to HSV and other viral infections of the brain. Overtime, the more frequent immune activation against infections may result in progressive neurodegeneration (Figure 6). Our findings suggest that AD is a type of immune disorder driven by virus or microbial infections of the brain. For microglia are the brain-resident immune cells and TREM2 gene plays a crucial role in microglia function (Keren-Shaul et al., 2017; Kober and Brett, 2017), we will study how they affect the progression of AD for further research.

FIGURE 6

FIGURE 6. The suggested AD mechanism according to our analysis.

Alzheimer’s Disease Neuroimaging Initiative

Data used in preparation of this article were obtained from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (adni.loni.usc.edu). As such, the investigators within the ADNI contributed to the design and implementation of ADNI and/or provided data but did not participate in analysis or writing of this report. A complete listing of ADNI investigators can be found at: http://adni.loni.usc.edu/wpcontent/uploads/how_to_apply/ADNI_Acknowledgement_List.pdf.

Author Contributions

CL analyzed the data and wrote the manuscript. JC, CS, and WZ helped with the interpretation of the results and writing of the manuscript. YX, ZJ, and HT helped with the study design and statistical analysis. XZ oversaw the overall research plan, the study design and statistical analysis. All authors read and approved the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

Data collection and sharing for this project was funded by the Alzheimer’s Disease Neuroimaging Initiative (ADNI) (National Institutes of Health grant U01 AG024904) and DOD ADNI (Department of Defense award number W81XWH-12-2-0012). ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, and through generous contributions from the following: AbbVie, Alzheimer’s Association; Alzheimer’s Drug Discovery Foundation; Araclon Biotech; BioClinica, Inc.; Biogen; Bristol-Myers Squibb Company; CereSpir, Inc.; Cogstate; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; EuroImmun; F. Hoffmann-La Roche Ltd., and its affiliated company Genentech, Inc.; Fujirebio; GE Healthcare; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC; Johnson & Johnson Pharmaceutical Research & Development LLC; Medpace, Inc.; Merck & Co, Inc.; Meso Scale Diagnostics, LLC.; NeuroRx Research; NeuroTrack Technologies; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Takeda Pharmaceutical Company; and Transition Therapeutics. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (www.fnih.org). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer’s Therapeutic Research Institute at the University of Southern California. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00410/full#supplementary-material

References

Barrett, T., Wilhite, S. E., Ledoux, P., Evangelista, C., Kim, I. F., Tomashevsky, M., et al. (2012). NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 41, D991–D995. doi: 10.1093/nar/gks1193

PubMed Abstract | CrossRef Full Text | Google Scholar

Bertram, L., Lill, C. M., and Tanzi, R. E. (2010). The genetics of Alzheimer Disease: back to the future. Neuron 68, 270–281. doi: 10.1016/j.neuron.2010.10.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Blalock, E. M., Buechel, H. M., Popovic, J., Geddes, J. W., and Landfield, P. W. (2011). Microarray analyses of laser-captured hippocampus reveal distinct gray and white matter signatures associated with incipient Alzheimer’s disease. J. Chem. Neuroanat. 42, 118–126. doi: 10.1016/j.jchemneu.2011.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Buerger, K., Ewers, M., Pirttila, T., Zinkowski, R., Alafuzoff, I., Teipel, S. J., et al. (2006). CSF phosphorylated tau protein correlates with neocortical neurofibrillary pathology in Alzheimer’s disease. Brain 129, 3035–3041. doi: 10.1093/brain/awl269

PubMed Abstract | CrossRef Full Text | Google Scholar

Clayton, D. (2012). snpStats: SnpMatrix and XSnpMatrix classes and methods. R Package.

Google Scholar

Daniel, J. M., and Reynolds, A. B. (1999). The Catenin p120 ctn Interacts with Kaiso, a novel BTB/POZ domain zinc finger transcription factor. Mol. Cell. Biol. 19, 3614–3623. doi: 10.1128/MCB.19.5.3614

PubMed Abstract | CrossRef Full Text | Google Scholar

Farrer, L. A. (1997). Effects of age, sex, and ethnicity on the association between apolipoprotein e genotype and Alzheimer disease: a meta-analysis. JAMA 278, 1349–1356. doi: 10.1001/jama.1997.03550160069041

PubMed Abstract | CrossRef Full Text | Google Scholar

Graham, D. P., Cully, J. A., Snow, A. L., Massman, P., and Doody, R. (2004). The Alzheimer’s Disease Assessment Scale-Cognitive subscale: normative data for older adult controls. Alzheimer Dis. Assoc. Disord. 18, 236–240.

Google Scholar

Hampel, H., Blennow, K., Shaw, L. M., Hoessler, Y. C., Zetterberg, H., and Trojanowski, J. Q. (2010). Total and phosphorylated tau protein as biological markers of Alzheimer’s disease. Exp. Gerontol. 45, 30–40. doi: 10.1016/j.exger.2009.10.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Harold, D., Abraham, R., Hollingworth, P., Sims, R., Gerrish, A., Hamshere, M. L., et al. (2009). Genome-wide association study identifies variants at CLU and PICALM associated with Alzheimer’s disease. Nat. Genet. 41, 1088–1093. doi: 10.1038/ng.440

PubMed Abstract | CrossRef Full Text | Google Scholar

Hesse, C., Rosengren, L., Andreasen, N., Davidsson, P., Vanderstichele, H., Vanmechelen, E., et al. (2001). Transient increase in total tau but not phospho-tau in human cerebrospinal fluid after acute stroke. Neurosci. Lett. 297, 187–190. doi: 10.1016/S0304-3940(00)01697-9

CrossRef Full Text | Google Scholar

Ikonomovic, M. D., Klunk, W. E., Abrahamson, E. E., Mathis, C. A., Price, J. C., Tsopelas, N. D., et al. (2008). Post-mortem correlates of in vivo PiB-PET amyloid imaging in a typical case of Alzheimer’s disease. Brain 131, 1630–1645. doi: 10.1093/brain/awn016

PubMed Abstract | CrossRef Full Text | Google Scholar

Itzhaki, R. F., and Wozniak, M. A. (2008). Herpes simplex virus type 1 in Alzheimer’s Disease: the enemy within. J. Alzheimers Dis. 13, 393–405. doi: 10.3233/JAD-2008-13405

CrossRef Full Text | Google Scholar

Ji, Z., Su, J., Liu, C., Wang, H., Huang, D., and Zhou, X. (2014). Integrating genomics and proteomics data to predict drug effects using binary linear programming. PLoS One 9:e102798. doi: 10.1371/journal.pone.0102798

PubMed Abstract | CrossRef Full Text | Google Scholar

Kent, W. J., Sugnet, C. W., Furey, T. S., Roskin, K. M., Pringle, T. H., Zahler, A. M., et al. (2002). The human genome browser at UCSC. Genome Res. 12, 996–1006. doi: 10.1101/gr.229102

PubMed Abstract | CrossRef Full Text | Google Scholar

Keren-Shaul, H., Spinrad, A., Weiner, A., Matcovitch-Natan, O., Dvir-Szternfeld, R., Ulland, T. K., et al. (2017). A unique microglia type associated with restricting development of Alzheimer’s Disease. Cell 169, 1276.e17–1290.e17. doi: 10.1016/j.cell.2017.05.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, A., Fornes, O., Stigliani, A., Gheorghe, M., Castro-Mondragon, J. A., van der Lee, R., et al. (2018). JASPAR 2018: update of the open-access database of transcription factor binding profiles and its web framework. Nucleic Acids Res. 46, D260–D266. doi: 10.1093/nar/gkx1126

PubMed Abstract | CrossRef Full Text | Google Scholar

Kheradpour, P., and Kellis, M. (2014). Systematic discovery and characterization of regulatory motifs in ENCODE TF binding experiments. Nucleic Acids Res. 42, 2976–2987. doi: 10.1093/nar/gkt1249

PubMed Abstract | CrossRef Full Text | Google Scholar

Kober, D. L., and Brett, T. J. (2017). TREM2-ligand interactions in health and disease. J. Mol. Biol. 429, 1607–1629. doi: 10.1016/j.jmb.2017.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, D. K. V., Choi, S. H., Washicosky, K. J., Eimer, W. A., Tucker, S., Ghofrani, J., et al. (2016). Amyloid-β peptide protects against microbial infection in mouse and worm models of Alzheimer’s disease. Sci. Transl. Med. 8:340ra72. doi: 10.1126/scitranslmed.aaf1059

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Wetten, S., Li, L., St. Jean, P. L., Upmanyu, R., Surh, L., et al. (2008). Candidate single-nucleotide polymorphisms from a genomewide association study of Alzheimer Disease. Arch. Neurol. 65, 45–53. doi: 10.1001/archneurol.2007.3

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Zhang, Q., Chen, F., Meng, X., Liu, W., Chen, D., et al. (2017). Genome-wide association and interaction studies of CSF T-tau/Aβ 42 ratio in ADNI cohort. Neurobiol. Aging 57, 247.e1–247.e8. doi: 10.1016/j.neurobiolaging.2017.05.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C.-C., Kanekiyo, T., Xu, H., and Bu, G. (2013). Apolipoprotein E and Alzheimer disease: risk, mechanisms and therapy. Nat. Rev. Neurol. 9, 106–118. doi: 10.1038/nrneurol.2012.263

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopez, M., Aoubala, M., Jordier, F., Isnardon, D., Gomez, S., and Dubreuil, P. (1998). The human poliovirus receptor related 2 protein is a new hematopoietic/endothelial homophilic adhesion molecule. Blood 92, 4602–4611.

PubMed Abstract | Google Scholar

Melani, M., Simpson, K. J., Brugge, J. S., and Montell, D. (2008). Regulation of cell adhesion and collective cell migration by hindsight and its human homolog RREB1. Curr. Biol. 18, 532–537. doi: 10.1016/j.cub.2008.03.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Michaelson, D. M. (2014). APOE ε4: the most prevalent yet understudied risk factor for Alzheimer’s disease. Alzheimers Dement. 10, 861–868. doi: 10.1016/j.jalz.2014.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Mukaetova-Ladinska, E. B., Abdel-All, Z., Mugica, E. S., Li, M., Craggs, L. J. L., Oakley, A. E., et al. (2015). Tau Proteins in the temporal and frontal cortices in patients with vascular dementia. J. Neuropathol. Exp. Neurol. 74, 148–157. doi: 10.1097/NEN.0000000000000157

PubMed Abstract | CrossRef Full Text | Google Scholar

Noguchi, M., Yoshita, M., Matsumoto, Y., Ono, K., Iwasa, K., and Yamada, M. (2005). Decreased β-amyloid peptide42 in cerebrospinal fluid of patients with progressive supranuclear palsy and corticobasal degeneration. J. Neurol. Sci. 237, 61–65. doi: 10.1016/j.jns.2005.05.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Ost, M., Nylen, K., Csajbok, L., Ohrfelt, A. O., Tullberg, M., Wikkelso, C., et al. (2006). Initial CSF total tau correlates with 1-year outcome in patients with traumatic brain injury. Neurology 67, 1600–1604. doi: 10.1212/01.wnl.0000242732.06714.0f

PubMed Abstract | CrossRef Full Text | Google Scholar

Pachkov, M., Balwierz, P. J., Arnold, P., Ozonov, E., and van Nimwegen, E. (2012). SwissRegulon, a database of genome-wide annotations of regulatory sites: recent updates. Nucleic Acids Res. 41, D214–D220. doi: 10.1093/nar/gks1145

PubMed Abstract | CrossRef Full Text | Google Scholar

Piacentini, R., Civitelli, L., Ripoli, C., Marcocci, M. E., De Chiara, G., Garaci, E., et al. (2011). HSV-1 promotes Ca2 + -mediated APP phosphorylation and Aβ accumulation in rat cortical neurons. Neurobiol. Aging 32, 2323.e13–2323.e26. doi: 10.1016/j.neurobiolaging.2010.06.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Podhorna, J., Krahnke, T., Shear, M., Harrison, J. E., and for the Alzheimer’s Disease Neuroimaging Initiative (2016). Alzheimer’s Disease assessment scale–cognitive subscale variants in mild cognitive impairment and mild Alzheimer’s disease: change over time and the effect of enrichment strategies. Alzheimers Res. Ther. 8:8. doi: 10.1186/s13195-016-0170-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Porcellini, E., Carbone, I., Ianni, M., and Licastro, F. (2010). Alzheimer’s disease gene signature says: beware of brain viral infections. Immun. Ageing 7:16. doi: 10.1186/1742-4933-7-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, D. L., and Sisodia, S. S. (1998). Mutant genes in familial Alzheimer’s disease and transgenic models. Annu. Rev. Neurosci. 21, 479–505.

Google Scholar

Price, J. L., and Morris, J. C. (1999). Tangles and plaques in nondemented aging and “preclinical” Alzheimer’s disease. Ann. Neurol. 45, 358–368.

Google Scholar

R Core Team (2013). R: A Language and Environment for Statistical Computing. Vienna: R Core Team.

Google Scholar

Readhead, B., Haure-Mirande, J.-V., Funk, C. C., Richards, M. A., Shannon, P., Haroutunian, V., et al. (2018). Multiscale analysis of independent Alzheimer’s cohorts finds disruption of molecular, genetic, and clinical networks by human herpesvirus. Neuron 99, 64.e7–82.e7. doi: 10.1016/j.neuron.2018.05.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43:e47. doi: 10.1093/nar/gkv007

PubMed Abstract | CrossRef Full Text | Google Scholar

Saykin, A. J., Shen, L., Foroud, T. M., Potkin, S. G., Swaminathan, S., Kim, S., et al. (2010). Alzheimer’s Disease Neuroimaging Initiative biomarkers as quantitative phenotypes: genetics core aims, progress, and plans. Alzheimers Dement. 6, 265–273. doi: 10.1016/j.jalz.2010.03.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Seibyl, J., Shaw, L. M., Blennow, K., Widmann, M., Corradini, V., Wahl, S., et al. (2017). AMYLOID-PET concordance of Elecsys^® CSF biomarker immunoassays for Alzheimer’s Disease. Alzheimers Dement. 13, 199–200. doi: 10.1016/j.jalz.2017.07.062

PubMed Abstract | CrossRef Full Text | Google Scholar

Shabalin, A. A. (2012). Matrix eQTL: ultra fast eQTL analysis via large matrix operations. Bioinformatics 28, 1353–1358. doi: 10.1093/bioinformatics/bts163

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon, P. (2014). MotifDb: An annotated collection of protein-DNA binding sequence motifs. R Package Version 1.

Google Scholar

Shaw, L. M., Vanderstichele, H., Knapik-Czajka, M., Clark, C. M., Aisen, P. S., Petersen, R. C., et al. (2009). Cerebrospinal fluid biomarker signature in Alzheimer’s disease neuroimaging initiative subjects. Ann. Neurol. 65, 403–413. doi: 10.1002/ana.21610

PubMed Abstract | CrossRef Full Text | Google Scholar

Shin, J.-H., Blay, S., Graham, J., and McNeney, B. (2006). LDheatmap: an R function for graphical display of pairwise linkage disequilibria between single nucleotide polymorphisms. J. Stat. Softw. 16, 1–9. doi: 10.18637/jss.v016.c03

CrossRef Full Text | Google Scholar

Spear, P. G. (2004). Herpes simplex virus: receptors and ligands for cell entry. Cell. Microbiol. 6, 401–410. doi: 10.1111/j.1462-5822.2004.00389.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Takai, Y. (2003). Nectin and afadin: novel organizers of intercellular junctions. J. Cell Sci. 116, 17–27. doi: 10.1242/jcs.00167

PubMed Abstract | CrossRef Full Text | Google Scholar

The Alzheimer’s Disease Neuroimaging Initiative, Verma, N., Beretvas, S. N., Pascual, B., Masdeu, J. C., and Markey, M. K. (2015). New scoring methodology improves the sensitivity of the Alzheimer’s Disease Assessment Scale-Cognitive subscale (ADAS-Cog) in clinical trials. Alzheimers Res. Ther. 7, 49–52. doi: 10.1186/s13195-015-0151-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Turner, S. D. (2014). qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. bioRxiv [Preprint]. doi: 10.1101/005165

CrossRef Full Text | Google Scholar

Van Dyck, F., Declercq, J., Braem, C., and Van de Ven, W. (2007). PLAG1, the prototype of the PLAG gene family: versatility in tumour development (Review). Int. J. Oncol. 30, 765–774. doi: 10.3892/ijo.30.4.765

PubMed Abstract | CrossRef Full Text | Google Scholar

Vawter, M. P., Evans, S., Choudary, P., Tomita, H., Meador-Woodruff, J., Molnar, M., et al. (2004). Gender-Specific gene expression in post-mortem human brain: localization to sex chromosomes. Neuropsychopharmacology 29, 373–384. doi: 10.1038/sj.npp.1300337

PubMed Abstract | CrossRef Full Text | Google Scholar

Ward, L. D., and Kellis, M. (2012). HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 40, D930–D934. doi: 10.1093/nar/gkr917

PubMed Abstract | CrossRef Full Text | Google Scholar

Warner, M. S., Geraghty, R. J., Martinez, W. M., Montgomery, R. I., Whitbeck, J. C., Xu, R., et al. (1998). A cell surface protein with herpesvirus entry activity (HveB) confers susceptibility to infection by mutants of herpes simplex Virus Type 1, herpes simplex virus type 2, and pseudorabies virus. Virology 246, 179–189. doi: 10.1006/viro.1998.9218

PubMed Abstract | CrossRef Full Text | Google Scholar

Wozniak, M. A., Itzhaki, R. F., Shipley, S. J., and Dobson, C. B. (2007). Herpes simplex virus infection causes cellular β-amyloid accumulation and secretase upregulation. Neurosci. Lett. 429, 95–100. doi: 10.1016/j.neulet.2007.09.077

PubMed Abstract | CrossRef Full Text | Google Scholar