Genome-Wide Association Studies for Cerebrospinal Fluid Soluble TREM2 in Alzheimer’s Disease

Alzheimer’s disease (AD) is the most common form of dementia. Rare variants in triggering receptor expressed on myeloid cells 2 (TREM2) have been identified as risk factors for AD. Soluble TREM2 (sTREM2) in the cerebrospinal fluid (CSF) is a potential and novel biomarker of neuroinflammation implicated in the onset and progression of AD. To explore the roles of CSF sTREM2 on the pathogenesis of AD, we performed genome-wide association studies (GWAS) by using the data from Alzheimer’s Disease Neuroimaging Initiative (ADNI). We found CSF sTREM2 levels were elevated with the disease stages, but there was no significant difference between that of AD patients and normal participants. CSF sTREM2 was positively correlated with CSF total tau and phosphorylated-tau levels (ρ > 0.35, p < 1e-06; ρ > 0.32, p < 1e-05, respectively) for all disease states. We identified the most significant CSF sTREM2 related locus was rs7232 (FDR = 3.01e-08), a missense variant in MS4A6A gene of chromosome 11. Moreover, we also detected rs7232 was highly associated with MS4A6A gene expression (FDR = 1.37e-18). In addition, our pathway analysis for our significant GWAS results showed that biological processes for regulation of viruses and immune response were highly overrepresented or enriched. Our study suggests that CSF sTREM2 plays an informative role in AD progression. Moreover, CSF sTREM2 and AD is highly related to viral infections and immune response.

Introduction

The most common neurodegenerative disorder, Alzheimer’s disease (AD), is characterized by accumulation of amyloid-β peptide (Aβ) in senile plaques and hyperphosphorylated tau protein in neurofibrillary tangles (NFT) in brain, neuroinflammation, and progressive decline in cognition and memory loss finally (Price and Sisodia, 1998; Mukaetova-Ladinska et al., 2015). Although AD is a complex, multigenic and multifactorial disorder, most genetic, physiological and pathological studies suggest that the balance between production and clearance of Aβ contributes to its development (Hardy and Higgins, 1992). Microglia are the brain-resident effective phagocytes for the uptake and proteolytic clearance of both soluble Aβ oligomers and insoluble Aβ fibrils (Lee and Landreth, 2010). Many studies have detected that microglia surrounded Aβ plaques in the brain of AD patients (Perlmutter et al., 1990) and AD mouse models (Stalder et al., 1999) to reduce their sizes and subsequent toxicity (Condello et al., 2015). The protein triggering receptor expressed on myeloid cells 2 (TREM2) is an innate immune receptor expressed on microglia and on myeloid cells outside the brain (Schmid et al., 2002). It plays crucial roles in microglial phagocytosis of apoptotic neurons, damaged myelin, and Aβ plaques (Jay et al., 2017). In addition, TREM2 regulates microglial proliferation, survival (Zheng et al., 2017), cytokine release (Zhong et al., 2015), and the accumulation around Aβ plaques (Yuan et al., 2016). Several studies reported that the rare TREM2 mutation, arginine 47 to histidine (p.R47H, rs75932628) substitution in the extracellular immunoglobulin domain, significantly increases the risk for AD with odds ratios similar to those of carrying an apolipoprotein E (APOE) ε4 allele (Guerreiro et al., 2013; Jonsson et al., 2013), which is present in about half of late-onset AD (LOAD) patients and has been convincingly demonstrated to affect its risk (Bertram et al., 2010).

TREM2 protein undergoes proteolytic cleavage by ADAM proteases, releasing its ectodomain into the extracellular space as the form of soluble TREM2 (sTREM2) (Wunderlich et al., 2013), which is abundantly detected in human cerebrospinal fluid (CSF) and plasma. Although elevated CSF sTREM2 levels were initially reported in neuroinflammatory conditions such as multiple sclerosis (Piccio et al., 2008), exploring the relationship between CSF sTREM2 and other biomarkers of AD has recently become of great interest. Recent study reported that CSF sTREM2 exerts functional roles in microglia by promoting inflammatory responses and shielding them from apoptosis (Zhong et al., 2017). Furthermore, an AD mouse model study demonstrates a protective role of sTREM2 against amyloid pathology and related toxicity (Zhong et al., 2019). However, the roles of CSF sTREM2 on the pathogenesis of AD remain unclear. To gain insight into the mechanism of AD development according to CSF sTREM2, in this study, we analyzed the correlation of CSF sTREM2 levels with AD status and other important CSF and clinical biomarkers from Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort. Moreover, we performed genome-wide association studies (GWAS) to identify novel variants and genes associated with sTREM2 level and AD by using CSF sTREM2 and gene expression data as phenotypes. Additionally, we implemented pathway analysis for sTREM2 related genes to get a better understanding of AD pathology.

Materials and Methods

Data used in this study were obtained from the ADNI database¹. 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.

We applied for and were granted permission to obtain data from the ADNI cohort² for performing the analyses described in this paper.

Subjects

In this work, analyses were restricted to ADNI subjects with CSF sTREM2 data available. The study sample (N = 1001) included 224 healthy normal (NL), 72 significant memory concern (SMC), 234 early mild cognitive impairment (EMCI), 277 late mild cognitive impairment (LMCI), and 194 AD participants. Table 1 shows selected demographic and clinical characteristics of these subjects at baseline.

TABLE 1

Table 1. Selected demographic and clinical characteristics of 1001 ADNI participants at baseline.

CSF and Clinical Biomarkers

Cerebrospinal fluid sTREM2 measurements were done with the ELISA protocol previously established by the Haass’ group with minor changes. The assay is based on the MSD platform and it is comprehensively described in their previous publications (Suárez-Calvet et al., 2016a, b). The CSF levels of amyloid-β 1-42 peptide (Aβ₄₂), total tau (tau) and tau phosphorylated at the threonine 181 (p-tau) were determined using the fully automated Roche Elecsys immunoassay platform (Seibyl et al., 2017). 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). ADAS13 scores were automatically calculated on the electronic case report form based on item level data entered. Semi-automated hippocampal volumetry was carried out using a commercially available high dimensional brain mapping tool (Medtronic Surgical Navigation Technologies (SNT), Louisville, CO, United States) for MRI data of subjects, that has previously been validated and compared to manual tracing of the hippocampus (Hsu et al., 2002).

Genotyping Data

The single-nucleotide polymorphisms (SNPs) 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; Shen et al., 2010). The SNPs shown in both arrays were used for the following analyses.

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) Hardy-Weinberg equilibrium (HWE) test of p-value >10^–6 (absolute value of z-score <4.753424). After QC analysis, 2,379,855 SNPs remained for the subsequent analyses.

Gene Expression Data

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

Statistical and Genetic Analyses

The CSF sTREM2 levels for different disease stages and their correlation with other CSF and clinical biomarkers were determined by using R software. To identify the association between SNPs and CSF sTREM2 levels, we performed (quantitative trait locus) QTL analysis using the R package Matrix expression quantitative trait loci (eQTL) (Shabalin, 2012). In this analysis, age, gender (1 for male and 2 for female) and diagnosis (1 for NL, 2 for SMC, 3 for EMCI, 4 for LMCI, and 5 for AD) at baseline were considered as covariates. Manhattan plots of QTL results were generated using the R package qqman (Turner, 2014). The regional plot for QTL results was obtained from LocusZoom⁴ (Pruim et al., 2010). The linkage disequilibrium (LD) plot was generated through R package gaston. The association between significant SNPs (FDR < 0.05) of our QTL results and gene expression data was detected by performing eQTL analysis. eQTL analysis was also conducted using the R package Matrix eQTL for the filtered SNPs, with the same covariates as our QTL analysis. Here, we considered both cis-eQTL analysis (local, distance < 1Mb) and trans-eQTL analysis (distant, distance ≥ 1Mb or the SNPs even locate on different chromosomes).

Pathway Analyses

For the significant genes (FDR < 0.05) of SNP-gene pairs form our eQTL results, gene ontology (GO) analyses were performed using the Protein Analysis Through Evolutionary Relationships (PANTHER) statistical over-representation test v14.1⁵, which used data from the Gene Ontology Consortium (GOC)⁶. PANTHER utilizes a binomial distribution test to calculate overrepresentation of candidate genes, relative to background, for different GO terms (Mi et al., 2019). Additionally, we also conducted gene set enrichment analysis (GSEA) for these genes through the software GSEA Desktop v3.0⁷ (Subramanian et al., 2005). Here, we applied R package limma (Ritchie et al., 2015) to make the pre-ranked gene list as the input for GSEA, according to the gene expression data for NL and AD subjects.

Results

Correlations of sTREM2 With AD Biomarkers

In ADNI cohort, we saw that CSF sTREM2 levels increased with the disease progression from SMC group to AD group (mean: 3801.212, 3891.837, 4176.13, and 4216.861 pg/mL for SMC, EMCI, LMCI, and AD subjects, respectively; Figure 1A). However, only the difference between EMCI state and AD state was significant (p < 0.05; Figure 1A). Though the mean CSF sTREM2 levels for AD patients was a little higher than that of normal subjects (4206.102 pg/mL), there was no significant difference.

FIGURE 1

Figure 1. Cerebrospinal fluid sTREM2 levels in the five clinical disease stages and the correlation between sTREM2 and other AD biomarkers. (A) Violin plots with boxplots for the comparison of CSF sTREM2 levels in disease states. Statistical significance was determined by Wilcoxon rank sum test. +: mean. ^∗p < 0.05. The correlation plots between log transformed CSF sTREM2 and (B) log transformed TREM2 gene expression in blood samples, (C) log transformed CSF Aβ₄₂, (D) log transformed CSF tau, (E) log transformed CSF p-tau, and (F) log transformed ADAS13 scores for each group. Black straight lines are the regression lines. Shaded areas around regression lines represent the pointwise 95% confidence intervals (CI). ρ: Spearman’s rank correlation coefficient (rho).

Levels of sTREM2 in the CSF did not correlate with TREM2 gene expression in the blood samples (Spearman’s ρ from −0.1007 to 0.1746, p > 0.05; Figure 1B). Next, we investigated the correlation of its CSF levels with other AD highly associated CSF and clinical biomarkers to explore its role in AD pathology. CSF sTREM2 levels showed highest correlation with CSF Aβ₄₂ in SMC group (ρ = 0.3307, p = 0.0046; Figure 1C), followed by AD group (ρ = 0.2639, p = 0.0002; Figure 1C). The overall correlation between CSF sTREM2 and CSF tau were high for all the five disease states (ρ from 0.3543 to 0.6546, p < 1e-06; Figure 1D). SMC group showed the highest correlation (ρ = 0.6546, p = 4.465e-10) while AD group had the lowest one (ρ = 0.3543, p = 4.292e-07; Figure 1D). Similarly, CSF sTREM2 also correlated with CSF p-tau for all the five categories (ρ from 0.3213 to 0.6398, p < 1e-05; Figure 1E). Same as the case for CSF tau, the correlation of SMC group ranked first while that of AD was the smallest. For ADAS13 scores, the correlation between them and CSF sTREM2 were low generally (ρ from −0.0137 to 0.1612; Figure 1F). Moreover, there was no strong correlation between CSF sTREM2 levels and hippocampus volumes (ρ from −0.1483 to 0.0082, p > 0.05; Supplementary Figure S1).

GWAS of sTREM2 Levels in CSF

We then performed QTL analysis to study the association between genotype (SNPs) and CSF sTREM2 levels. According to our analysis, the significant loci concentrated in chromosome 11 (Figure 2A). The most significant one was rs7232 (p = 1.32e-14, FDR = 3.01e-08; Table 2 and Figure 2), a missense variant in MS4A6A gene locus. The regional association plot of rs7232 is shown in Figure 2B. The second most significant SNP was rs1582763 (p = 3.85e-14, FDR = 4.41e-08; Table 2), which locates near the MS4A4E gene. The top 10 most significant SNPs were all located on chromosome 11. Moreover, they were either within or near to MS4A gene family (Table 2). Besides identified SNPs on chromosome 11, there were also SNPs showed significant association on other chromosomes, such as rs3799468 (p = 4.10e-12, FDR = 2.85e-07) located in MAP7 gene on chromosome 6, rs181768270 (p = 4.19e-12, FDR = 2.85e-07) proximal to gene RNF187 on chromosome 1, and rs116087048 (p = 4.19e-12, FDR = 2.85e-07; Supplementary Table S1) within the gene ASIC4 on chromosome 2. The list for all identified SNPs with p < 0.05 from our QTL analysis can be found in Supplementary Table S1.

FIGURE 2

Figure 2. Manhattan plot and regional plot of the results from QTL analysis for CSF sTREM2 levels. (A) Manhattan plot of –log₁₀ (p-value) from the results of QTL analysis. (B) Regional plot for the most significant SNP rs7232 identified by the QTL analysis. The r² measures the linkage disequilibrium of each SNP with the most significant SNP rs7232 according to hg19/1000 Genomes Nov 2014 AMR population. Points with gray color indicates that their r² values were not available in the reference genome.

TABLE 2

Table 2. Top 10 most significant SNPs identified by the QTL analysis for CSF sTREM2.

Next, for exploring the interactions among the identified significant SNPs on chromosome 11, we selected the SNPs with FDR < 0.05 and located within ±1Mb from rs7232. 40 SNPs including rs7232 remained for the analysis. The LD pattern for these 40 SNPs is illustrated in Figure 3. rs7232 showed strong LD with all the rest of top 10 identified most significant SNPs (Figure 3). Among them, the weakest LD occurred between rs7232 and rs1582763 (r² = 0.68), while the LD between rs7232 and rs12453 was the strongest (r² = 0.83; Figure 3). The strong LD between rs7232 and rs12453 may be due to that both of them locate in same gene MS4A6A. The LD of the SNP pairs rs1151065-rs558788, rs11230160-rs920573, and rs6591559-rs1530914 were extremely high (r² = 1.00 for all the three SNP pairs; Figure 3).

FIGURE 3

Figure 3. The linkage disequilibrium pattern for the identified SNPs close to rs7232 in our study. Here the LD measure is r².

For we wanted to study how the SNPs associated with CSF sTREM2 levels regulated gene expression, we chose the SNPs with FDR < 0.05 identified by our QTL analysis for the following eQTL analysis. There were 240 SNPs left for this analysis. In our cis-eQTL analysis (local eQTL analysis, the distance between SNPs and associated genes is less than 1Mb), the top 10 identified most significant SNPs were all located on chromosome 11, and they all associated with MS4A6A gene (Table 3). rs12453 ranked first (p = 9.02e-23, FDR = 8.14e-19; Table 3) in the results. Another SNP within MS4A6A gene, rs7232, showed the third most significant association (p = 4.57e-22, FDR = 1.37e-18) with it. Furthermore, rs1582763 was also significantly associated with MS4A6A gene expression (p = 3.58e-17, FDR = 1.9e-14; Supplementary Table S2). For the results of trans-eQTL analysis (distant eQTL analysis, the distance between SNPs and associated genes is bigger than 1Mb, or they even locate on different chromosomes), the association between rs113239743 and CPT1C gene was the most significant (p = 6.53e-44, FDR = 7.73e-37; Table 4). Moreover, rs113239743 also showed significant association with TESMIN gene (p = 3.91e-30, FDR = 9.27e-24), and CCND2 gene (p = 2.43e-25, FDR = 2.39e-19). The results for all identified SNP-gene pairs with p < 0.05 from our cis-eQTL analysis and those with FDR < 0.05 from trans-eQTL analysis are demonstrated in Supplementary Tables S2, S3, respectively.

TABLE 3

Table 3. Top 10 most significant SNP-gene pairs identified by our cis-eQTL analysis.

TABLE 4

Table 4. Top 10 most significant SNP-gene pairs identified by our trans-eQTL analysis.

Pathway Studies

For the SNP-gene pairs identified from both cis-eQTL and trans-eQTL analysis for CSF sTREM2, we chose the genes from the results with FDR < 0.05 for the following pathway analysis. We found 4295 qualified genes. Firstly, we performed statistical overrepresentation test for these genes. According to the value of fold enrichment, the highest one was negative regulation of viral process (fold enrichment = 2.27, p = 5.16e-03; Table 5). What’s more, top three categories were all related to virus regulation. In this list of top 10 categories, four of them were associated with viral process and one was in connection with immune response (neutrophil activation involved in immune response; Table 5). In addition, the most significant category was neutrophil activation (fold enrichment = 1.85, p = 3.76e-11; Table 5). The results for all identified GO categories with corrected p-value less than 0.05 from this analysis are illustrated in Supplementary Table S4.

TABLE 5

Table 5. Top 10 identified gene ontology categories from PANTHER overrepresentation test of our eQTL analysis results, according to fold enrichment.

In addition to overrepresentation test, we also implemented GSEA for these over 4000 significant genes identified by our eQTL analysis. For the gene sets (GO categories) upregulated (positively enriched) in AD, 10 of them were significant at FDR < 0.25. On the other hand, there were only 3 such categories for gene sets downregulated (negatively enriched) in AD (Figure 4). The top positively enriched GO category was immune response (NES = 2.4469, p < 0.001, FDR = 0.0859), while the top negatively enriched gene set was myeloid cell development (NES = −2.6369, p < 0.001, FDR = 0.0624; Figure 4). There was another GO category related to immune response (regulation of immune response, NES = 2.4966, p < 0.001, FDR = 0.1136) in the top 10 significant positively enriched gene sets, which ranked third. Moreover, two of the significant downregulated gene sets were for myeloid cells. The full lists for upregulated and downregulated results with p < 0.05 of GSEA are in Supplementary Tables S5, S6, respectively.

FIGURE 4

Figure 4. Bar plot for enriched gene ontology categories with FDR < 0.25 from GSEA of our eQTL analysis results. The numbers on the bars are the normalized enrichment scores for the corresponding gene ontology categories. NES: normalized enrichment score.

Discussion

In 2014, a study revealed a significant reduction of CSF sTREM2 in AD patients (n = 56) compared to control individuals (n = 88, p = 0.001; Kleinberger et al., 2014). However, 2 years later, another study reported that the CSF sTREM2 levels were significant higher in AD cases compared to controls [median and range: 1028 (244–2570) and 832 (163–2196) pg/ml, respectively; p = 0.015] in the Knight Alzheimer’s Disease Research Center (ADRC) cohort with 73 AD subjects and 107 controls (Piccio et al., 2016). Here, we detected that the difference of CSF sTREM2 levels between AD and normal participants was not significant in ADNI cohort (Figure 1A). These inconsistent results from different cohort studies may be caused by the timing for measuring the CSF sTREM2 levels in AD patients. Some studies for different cohorts report that the CSF sTREM2 increased in early symptomatic stages of AD, but it decreased in late stage of the disease (Heslegrave et al., 2016; Suárez-Calvet et al., 2016b; Liu D. et al., 2018). The mechanisms behind such findings still remain elusive. A possible reason would be that a part of sTREM2 released from microglia is retained within the barrier formed by microglia around plaques in the late stage of AD (Condello et al., 2015; Wang et al., 2016). Longitudinal data is required for further study about the relation between CSF sTREM2 levels and disease states. Additionally, another explanation for the inconsistency would be the different proportions of subjects with TREM2 variants in their studies. Some TREM2 variants affect CSF sTREM2 levels dramatically (Piccio et al., 2016; Wolfe et al., 2018).

CSF total tau and p-tau are passively released into CSF by dying neurons. CSF total tau levels reflect the intensity of neuronal damage and degeneration, with higher amounts of tau releasing to CSF for more intense neurodegeneration. Similar to total tau, some studies report correlation between high CSF p-tau and higher rate of cognitive decline (Blennow et al., 2010). We showed that CSF sTREM2 levels were positively correlated with CSF total tau and phosphorylated tau levels for all disease status in ADNI cohort (Figures 1D,E), which were also reported by some previous studies (Heslegrave et al., 2016; Piccio et al., 2016). The mechanism for CSF sTREM2 in pathogenesis of AD is still unclear based on these preliminary findings. Further studies are needed to reveal its role in neurodegeneration. As for CSF Aβ₄₂, some studies showed that its levels were not correlated with CSF sTREM2 level generally (Heslegrave et al., 2016; Piccio et al., 2016). However, our analyses by disease stages detected that there were positive correlations between them for all status except for EMCI group in ADNI. Furthermore, the correlation was strongest for SMC group (Figure 1C). These imply CSF sTREM2 may play a crucial role at the very early symptomatic stage of AD.

Our QTL analysis showed that rs7232 was the most significant SNP associated with CSF sTREM2 level. A recent study for Chinese Alzheimer’s Biomarker and Lifestyle (CABLE) cohort (Hou et al., 2019) also reported this association was significant (p = 0.00106). Moreover, one study showed that rs7232 correlated with atrophy rate of left middle temporal and minor T allele carriers had less loss in the volume of left middle temporal than A allele homozygotes subjects (Ma et al., 2016). In addition, rs7232 was identified as a protective locus for AD (Lambert et al., 2013) in combined dataset of International Genomics of Alzheimer’s Project (IGAP) stages 1 and 2 (OR: 0.90, 95%CI: 0.87–0.92, p = 2.621 × 10^–14). Our eQTL analysis detected that rs7232 was significantly associated with its located gene MS4A6A, which was supported by other studies (Proitsi et al., 2014). Membrane Spanning 4-Domains A6A (MS4A6A) has been identified as one of the significantly associated loci with AD (Hollingworth et al., 2011; Deng et al., 2012). MS4A6A expression levels were found to be associated with elevated Braak tangle and Braak plaque scores (Karch et al., 2012). Furthermore, one study reported that MS4A6A expression was significantly correlated to AD-related neurofibrillary pathology and tau phosphorylation (Martiskainen et al., 2015). Additionally, the expression of MS4A6A is highly correlated with the expression of TREM2 in the brain (Hou et al., 2019). For both MS4A6A and TREM2 are mainly expressed on microglia cells in the brain (Darmanis et al., 2015), we speculate that MS4A6A may regulate TREM2 expression and the levels of CSF sTREM2. The second most significant SNP associated with CSF sTREM2 was rs1582763 from our GWAS (Table 2), which was also highly associated with MS4A6A gene expression according to our eQTL analysis. These associations were justified by studies from other cohorts (Deming et al., 2018). Similar to rs7232, rs1582763 was identified to be associated with reduced AD risk (Jun et al., 2016). Our GWAS demonstrated that the top significant SNPs associated with CSF sTREM2 levels were located within or near the MS4A gene cluster and the genes of top significant SNP-gene pairs were also from MS4A gene family (Tables 2, 3). The MS4A gene cluster encodes cell membrane proteins. In addition to possibly being involved in the regulation of calcium signaling, MS4A gene cluster has also been reported to be involved in immune-system function (Ma et al., 2015). More recently, a study showed that MS4A4A co-localized with TREM2 in human macrophages and both proteins were upregulated in response to IL-4-mediated stimulation (Deming et al., 2019). They also found that antibody-mediated targeting of MS4A4A was sufficient to reduce sTREM2 levels in human macrophages. All these provide strong evidence of a biological relation between TREM2 and proteins in the MS4A gene cluster. Additional studies are required to understand the mechanisms for how MS4A gene cluster modulates TREM2, and affects sTREM2 levels.

As the most significant variant associated with AD risk in TREM2 coding regions, p.R47H (rs75932628) was also genotyped in ADNI cohort. However, there are no subjects carrying its minor allele T for the available SNP data. Besides p.R47H, the minor allele (T) of TREM2 rare variant p.R62H (rs143332484) was also reported to be significantly associated with increased AD risk (Cuyvers et al., 2014). But it was not genotyped in ADNI genotype data. These may be due to TREM2 coding variants present very low frequency among people [minor allele frequency (MAF) < 0.5%]. Furthermore, some other TREM2 variants have been identified to be associated with AD status (Jin et al., 2014; Guerreiro et al., 2013), including p.T96K (rs2234253), p.H157Y (rs2234255), p.L211P (rs2234256), and p.W191X (rs2234258). Our GWAS showed that p.T96K and p.L211P were associated with CSF sTREM2 levels (p = 0.005946, 0.009005, respectively; Supplementary Table S1), while p.H157Y and p.W191X were not as a result of their extreme low MAF. These imply that TREM2 variants may modify CSF sTREM2 levels.

For the pathway analyses, we found that some biological processes related to virus regulation were significantly enriched from our statistical overrepresentation test (Table 5). Interestingly, a previous study for a different AD cohort reported similar results (Piccio et al., 2016), which identified that the top two significant GO categories were receptor-mediated endocytosis of virus by host cell and endocytosis involved in viral entry into host cell (both with fold enrichment > 5 and p = 4.79e-06). On the other hand, our GSEA showed that some categories for immune response were significantly positively enriched (Figure 4). What’s more, our previous study (Liu C. et al., 2018) illustrated that how poliovirus receptor-related 2 (PVRL2), a gene located near the APOE locus and mediating the entry of herpes simplex virus (HSV), contributed to the progression of AD. Additionally, some studies reported that the accumulation of Aβ plaque deposits might be a consequence of the over-production of Aβ peptide during viral infection of the brain and Aβ peptide acted as a defense molecule of the innate immunity (Porcellini et al., 2010; Kumar et al., 2016). The viral hypothesis for AD development has been proposed almost immediately after the first case of AD was reported by Alois Alzheimer in 1907 (Lawrence, 2017). Researchers started to focus on the herpes simplex virus type 1 (HSV1) at the beginning of 1980s (Ball, 1982; Gannicliffe et al., 1986). So far, more than 100 studies have connected AD with some forms of pathogen (Itzhaki et al., 2016). People may argue that, if viruses do play a causal role in AD, then why there are so many individuals who are infected with the viruses but without developing AD? This may be due to the feature of many types of virus and bacteria that it is possible to be infected but without being affected. In recent studies, an important issue about the association of viruses with the substantial accumulation of Aβ in subjects who are cognitively normal is still not addressed (Eimer et al., 2018; Readhead et al., 2018). The researchers suggest that it may be the case that Aβ is protective in such cases of participants (Balin and Hudson, 2018; Eimer et al., 2018). If viruses or other microbes contribute to the development of AD, people may reason that could antiviral drug inhibit AD progression. Some research illustrate that antivirals have such effect. In one study, researchers found that three different antiviral agents reduced Aβ and p-tau accumulation for Vero cell cultures infected with HSV1 (Wozniak et al., 2011). A recent nationwide cohort study from Taiwan provides the first population-level evidence. This study enrolled 33,000 subjects that about one-quarter of whom were newly diagnosed with HSV infection and followed them for 15 years. The researchers found that subjects with HSV infection showed 2.5 times more likely to develop dementia than those without infection. More importantly, compared with HSV infected participants without receiving therapy, subjects with HSV who were treated with antivirals had a 10-fold reduction in the risk of dementia development (Tzeng et al., 2018). These studies show a potential causal role for viruses in AD. However, the exact nature of the link between viruses and AD is still ambiguous. More clinical trials of antiviral drugs are needed to evaluate their impacts on AD.

Conclusion

We started from CSF sTREM2 levels to explore the pathogenesis of AD according to the data in ADNI cohort. Firstly, we checked the CSF sTREM2 levels by disease status and found that there was no significant difference between AD and NL groups, though they were elevated with disease stages. Next, we studied the correlation between CSF sTREM2 levels and other AD highly associated CSF and clinical biomarkers, such as CSF tau, p-tau, Aβ₄₂ levels, ADAS13 scores and hippocampus volumes, demonstrating that CSF sTREM2 levels were significantly positively correlated with CSF total tau and phosphorylated-tau levels for all disease status. Furthermore, we performed QTL analysis by setting CSF sTREM2 levels as the phenotype and identified that the SNPs located within or near the MS4A gene cluster were significantly associated with them, with rs7232 and rs1582763 as the top two significant SNPs. After that, our eQTL analysis for these significant SNPs showed that they were also associated with the expression of the genes from MS4A gene family. Additionally, our pathway analyses for the important genes from the results of GWAS illustrated that the biological processes for virus regulation and immune response were highly associated with AD. According to our study, we speculated that the genetic architecture of AD patients might increase their susceptibility to viral infections of the brain. The more frequent immune activation and response against viral infections may result in progressive neurodegeneration and lead to AD finally.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: http://adni.loni.usc.edu.

Author Contributions

CL analyzed the data and wrote the manuscript. JY helped with the interpretation of the results, writing of the manuscript, the study design, and statistical analysis. Both authors oversaw the overall research plan, the study design, and statistical analysis, and read and approved the final manuscript.

Conflict of Interest

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.

Funding

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 was 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 the ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (http://www.fnih.org). The grantee organization is the Northern California Institute for Research and Education, and the study was 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/fnagi.2019.00297/full#supplementary-material

Footnotes

1. ^ http://adni.loni.usc.edu
2. ^ http://www.adni-info.org
3. ^ http://www.affymetrix.com
4. ^ http://locuszoom.org
5. ^ http://pantherdb.org
6. ^ http://www.geneontology.org
7. ^ http://software.broadinstitute.org/gsea/index.jsp

References

Balin, B. J., and Hudson, A. P. (2018). Herpes viruses and Alzheimer’s disease: new evidence in the debate. Lancet Neurol. 17, 839–841. doi: 10.1016/s1474-4422(18)30316-8

CrossRef Full Text | Google Scholar

Ball, M. J. (1982). Limbic predilection in Alzheimer Dementia: is reactivated herpesvirus involved? Can. J. Neurol. Sci. J. Can. Sci. Neurol. 9, 303–306. doi: 10.1017/S0317167100044115

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

Blennow, K., Hampel, H., Weiner, M., and Zetterberg, H. (2010). Cerebrospinal fluid and plasma biomarkers in Alzheimer disease. Nat. Rev. Neurol. 6, 131–144. doi: 10.1038/nrneurol.2010.4

PubMed Abstract | CrossRef Full Text | Google Scholar

Clayton, D. (2012). Snpstats: SnpMatrix and XSnpMatrix Classes and Methods. R Package.

Google Scholar

Condello, C., Yuan, P., Schain, A., and Grutzendler, J. (2015). Microglia constitute a barrier that prevents neurotoxic protofibrillar Aβ42 hotspots around plaques. Nat. Commun. 6:6176. doi: 10.1038/ncomms7176

PubMed Abstract | CrossRef Full Text | Google Scholar

Cuyvers, E., Bettens, K., Philtjens, S., Van Langenhove, T., Gijselinck, I., van der Zee, J., et al. (2014). Investigating the role of rare heterozygous TREM2 variants in Alzheimer’s disease and frontotemporal dementia. Neurobiol. Aging 35:26.e11-726.e19. doi: 10.1016/j.neurobiolaging.2013.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Darmanis, S., Sloan, S. A., Zhang, Y., Enge, M., Caneda, C., Shuer, L. M., et al. (2015). A survey of human brain transcriptome diversity at the single cell level. Proc. Natl. Acad. Sci. U.S.A. 112, 7285–7290. doi: 10.1073/pnas.1507125112

PubMed Abstract | CrossRef Full Text | Google Scholar

Deming, Y., Filipello, F., Cignarella, F., Cantoni, C., Hsu, S., Mikesell, R., et al. (2019). The MS4A gene cluster is a key modulator of soluble TREM2 and Alzheimer’s disease risk. Sci. Transl. Med. 11:eaau2291. doi: 10.1126/scitranslmed.aau2291

PubMed Abstract | CrossRef Full Text | Google Scholar

Deming, Y., Filipello, F., Cignarella, F., Hsu, S., Mikesell, R., Li, Z., et al. (2018). The MS4A Gene Cluster is a Key Regulator of Soluble TREM2 and Alzheimer Disease risk biorxiv. [preprint]. doi: 10.1101/352179.

CrossRef Full Text | Google Scholar

Deng, Y.-L., Liu, L.-H., Wang, Y., Tang, H.-D., Ren, R.-J., Xu, W., et al. (2012). The prevalence of CD33 and MS4A6A variant in Chinese Han population with Alzheimer’s disease. Hum. Genet. 131, 1245–1249. doi: 10.1007/s00439-012-1154-1156

PubMed Abstract | CrossRef Full Text | Google Scholar

Eimer, W. A., Vijaya Kumar, D. K., Navalpur Shanmugam, N. K., Rodriguez, A. S., Mitchell, T., Washicosky, K. J., et al. (2018). Alzheimer’s disease-associated β-Amyloid is rapidly seeded by herpesviridae to protect against brain infection. Neuron 99, 56.e3–63.e3. doi: 10.1016/j.neuron.2018.06.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Gannicliffe, A., Sutton, R. N. P., and Itzhaki, R. F. (1986). Viruses, brain and immunosuppression. Psychol. Med. 16, 247–249. doi: 10.1017/S0033291700009053

PubMed Abstract | CrossRef Full Text | Google Scholar

Guerreiro, R., Wojtas, A., Bras, J., Carrasquillo, M., Rogaeva, E., Majounie, E., et al. (2013). TREM2 variants in Alzheimer’s Disease. N. Engl. J. Med. 368, 117–127. doi: 10.1056/NEJMoa1211851

PubMed Abstract | CrossRef Full Text | Google Scholar

Hardy, J., and Higgins, G. (1992). Alzheimer’s disease: the amyloid cascade hypothesis. Science 256, 184–185. doi: 10.1126/science.1566067

PubMed Abstract | CrossRef Full Text | Google Scholar

Heslegrave, A., Heywood, W., Paterson, R., Magdalinou, N., Svensson, J., Johansson, P., et al. (2016). Increased cerebrospinal fluid soluble TREM2 concentration in Alzheimer’s disease. Mol. Neurodegener. 11:3. doi: 10.1186/s13024-016-0071-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hollingworth, P., Harold, D., Sims, R., Gerrish, A., Lambert, J. C., Carrasquillo, M. M., et al. (2011). Common variants at ABCA7, MS4A6A/MS4A4E, EPHA1, CD33 and CD2AP are associated with Alzheimer’s disease. Nat. Genet. 43, 429–435. doi: 10.1038/ng.803

PubMed Abstract | CrossRef Full Text | Google Scholar

Hou, X.-H., Bi, Y.-L., Tan, M.-S., Xu, W., Li, J.-Q., Shen, X.-N., et al. (2019). Genome-wide association study identifies Alzheimer’s risk variant in MS4A6A influencing cerebrospinal fluid sTREM2 levels. Neurobiol. Aging S0197-4580, 30143–30145. doi: 10.1016/j.neurobiolaging.2019.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsu, Y.-Y., Schuff, N., Du, A.-T., Mark, K., Zhu, X., Hardin, D., et al. (2002). Comparison of automated and manual MRI volumetry of hippocampus in normal aging and dementia. J. Magn. Reson. Imaging 16, 305–310. doi: 10.1002/jmri.10163

PubMed Abstract | CrossRef Full Text | Google Scholar

Itzhaki, R. F., Lathe, R., Balin, B. J., Ball, M. J., Bearer, E. L., Braak, H., et al. (2016). Microbes and Alzheimer’s Disease. J. Alzheimers Dis. 51, 979–984. doi: 10.3233/JAD-160152

PubMed Abstract | CrossRef Full Text | Google Scholar

Jay, T. R., von Saucken, V. E., and Landreth, G. E. (2017). TREM2 in Neurodegenerative Diseases. Mol. Neurodegener. 12:56. doi: 10.1186/s13024-017-0197-195

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, S. C., Benitez, B. A., Karch, C. M., Cooper, B., Skorupa, T., Carrell, D., et al. (2014). Coding variants in TREM2 increase risk for Alzheimer’s disease. Hum. Mol. Genet. 23, 5838–5846. doi: 10.1093/hmg/ddu277

PubMed Abstract | CrossRef Full Text | Google Scholar

Jonsson, T., Stefansson, H., Steinberg, S., Jonsdottir, I., Jonsson, P. V., Snaedal, J., et al. (2013). Variant of TREM2 associated with the risk of Alzheimer’s Disease. N. Engl. J. Med. 368, 107–116. doi: 10.1056/NEJMoa1211103

PubMed Abstract | CrossRef Full Text | Google Scholar

Jun, G., Ibrahim-Verbaas, C. A., Vronskaya, M., Lambert, J.-C., and Chung, J. (2016). A novel Alzheimer disease locus located near the gene encoding tau protein. Mol. Psychiatry 21, 108–117. doi: 10.1038/mp.2015.23

PubMed Abstract | CrossRef Full Text | Google Scholar

Karch, C. M., Jeng, A. T., Nowotny, P., Cady, J., Cruchaga, C., and Goate, A. M. (2012). Expression of Novel Alzheimer’s disease risk genes in control and Alzheimer’s Disease brains. PLoS One 7:e50976. doi: 10.1371/journal.pone.0050976

PubMed Abstract | CrossRef Full Text | Google Scholar

Kleinberger, G., Yamanishi, Y., Suarez-Calvet, M., Czirr, E., Lohmann, E., Cuyvers, E., et al. (2014). TREM2 mutations implicated in neurodegeneration impair cell surface transport and phagocytosis. Sci. Transl. Med. 6:243ra86. doi: 10.1126/scitranslmed.3009093

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

Lambert, J.-C., Ibrahim-Verbaas, C. A., Harold, D., Naj, A. C., Sims, R., Bellenguez, C., et al. (2013). Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat. Genet. 45, 1452–1458. doi: 10.1038/ng.2802

PubMed Abstract | CrossRef Full Text | Google Scholar

Lawrence, B. (2017). Are the infectious roots of alzheimers buried deep in the past? J. Mol. Pathol. Epidemol. 3:2.

Google Scholar

Lee, C. Y. D., and Landreth, G. E. (2010). The role of microglia in amyloid clearance from the AD brain. J. Neural. Transm. 117, 949–960. doi: 10.1007/s00702-010-0433-434

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Chyr, J., Zhao, W., Xu, Y., Ji, Z., Tan, H., et al. (2018). Genome-wide association and mechanistic studies indicate that immune response contributes to Alzheimer’s Disease development. Front. Genet. 9:410. doi: 10.3389/fgene.2018.00410

CrossRef Full Text | Google Scholar

Liu, D., Cao, B., Zhao, Y., Huang, H., McIntyre, R. S., Rosenblat, J. D., et al. (2018). Soluble TREM2 changes during the clinical course of Alzheimer’s disease: a meta-analysis. Neurosci. Lett. 686, 10–16. doi: 10.1016/j.neulet.2018.08.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, J., Yu, J.-T., and Tan, L. (2015). MS4A Cluster in Alzheimer’s Disease. Mol. Neurobiol. 51, 1240–1248. doi: 10.1007/s12035-014-8800-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, J., Zhang, W., Tan, L., Wang, H.-F., Wan, Y., Sun, F.-R., et al. (2016). MS4A6A genotypes are associated with the atrophy rates of Alzheimer&#x2019;s disease related brain structures. Oncotarget 7, :58779–58788. doi: 10.18632/oncotarget.9563

PubMed Abstract | CrossRef Full Text | Google Scholar

Martiskainen, H., Viswanathan, J., Nykänen, N.-P., Kurki, M., Helisalmi, S., Natunen, T., et al. (2015). Transcriptomics and mechanistic elucidation of Alzheimer’s disease risk genes in the brain and in vitro models. Neurobiol. Aging 36:1221.e15-1221.e28. doi: 10.1016/j.neurobiolaging.2014.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Mi, H., Muruganujan, A., Ebert, D., Huang, X., and Thomas, P. D. (2019). PANTHER version 14: more genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nucleic Acids Res. 47, D419–D426. doi: 10.1093/nar/gky1038

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

Perlmutter, L. S., Barron, E., and Chui, H. C. (1990). Morphologic association between microglia and senile plaque amyloid in Alzheimer’s disease. Neurosci. Lett. 119, 32–36. doi: 10.1016/0304-3940(90)90748-X

CrossRef Full Text | Google Scholar

Piccio, L., Buonsanti, C., Cella, M., Tassi, I., Schmidt, R. E., Fenoglio, C., et al. (2008). Identification of soluble TREM-2 in the cerebrospinal fluid and its association with multiple sclerosis and CNS inflammation. Brain 131, 3081–3091. doi: 10.1093/brain/awn217

PubMed Abstract | CrossRef Full Text | Google Scholar

Piccio, L., Deming, Y., Del-Águila, J. L., Ghezzi, L., Holtzman, D. M., Fagan, A. M., et al. (2016). Cerebrospinal fluid soluble TREM2 is higher in Alzheimer disease and associated with mutation status. Acta Neuropathol 131, 925–933. doi: 10.1007/s00401-016-1533-1535

CrossRef Full Text | Google Scholar

Podhorna, J., Krahnke, T., Shear, M., and Harrison, J. E. 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-175

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. doi: 10.1146/annurev.neuro.21.1.479

CrossRef Full Text | Google Scholar

Proitsi, P., Lee, S. H., Lunnon, K., Keohane, A., Powell, J., Troakes, C., et al. (2014). Alzheimer’s disease susceptibility variants in the MS4A6A gene are associated with altered levels of MS4A6A expression in blood. Neurobiol. Aging 35, 279–290. doi: 10.1016/j.neurobiolaging.2013.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Pruim, R. J., Welch, R. P., Sanna, S., Teslovich, T. M., Chines, P. S., Gliedt, T. P., et al. (2010). LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 26, 2336–2337. doi: 10.1093/bioinformatics/btq419

PubMed Abstract | CrossRef Full Text | 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–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

Schmid, C. D., Sautkulis, L. N., Danielson, P. E., Cooper, J., Hasel, K. W., Hilbush, B. S., et al. (2002). Heterogeneous expression of the triggering receptor expressed on myeloid cells-2 on adult murine microglia. J. Neurochem. 83, 1309–1320. doi: 10.1046/j.1471-4159.2002.01243.x

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

Shen, L., Kim, S., Risacher, S. L., Nho, K., Swaminathan, S., West, J. D., et al. (2010). Whole genome association study of brain-wide imaging phenotypes for identifying quantitative trait loci in MCI and AD: a study of the ADNI cohort. NeuroImage 53, 1051–1063. doi: 10.1016/j.neuroimage.2010.01.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Stalder, M., Phinney, A., Probst, A., Sommer, B., Staufenbiel, M., and Jucker, M. (1999). Association of microglia with amyloid plaques in brains of APP23 transgenic mice. Am. J. Pathol. 154, 1673–1684. doi: 10.1016/S0002-9440(10)65423-65425

PubMed Abstract | CrossRef Full Text | Google Scholar

Suárez-Calvet, M., Araque Caballero, M. A., Kleinberger, G., Bateman, R. J., Fagan, A. M., Morris, J. C., et al. (2016a). Early changes in CSF sTREM2 in dominantly inherited Alzheimers disease occur after amyloid deposition and neuronal injury. Sci. Transl. Med. 8, 369ra178. doi: 10.1126/scitranslmed.aag1767

PubMed Abstract | CrossRef Full Text | Google Scholar

Suárez-Calvet, M., Kleinberger, G., Araque Caballero, M. Á, Brendel, M., Rominger, A., Alcolea, D., et al. (2016b). Cerebrospinal fluid levels are a potential biomarker for microglia activity in early−stage Alzheimer’s disease and associate with neuronal injury markers. EMBO Mol. Med. 8, 466–476. doi: 10.15252/emmm.201506123

PubMed Abstract | CrossRef Full Text | Google Scholar

Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U.S.A. 102, 15545–15550. doi: 10.1073/pnas.0506580102

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

Tzeng, N.-S., Chung, C.-H., Lin, F.-H., Chiang, C.-P., Yeh, C.-B., Huang, S.-Y., et al. (2018). Anti-herpetic medications and reduced risk of dementia in patients with herpes simplex virus infections—a nationwide, population-based cohort study in taiwan. Neurotherapeutics 15, 417–429. doi: 10.1007/s13311-018-0611-x

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

Wang, Y., Ulland, T. K., Ulrich, J. D., Song, W., Tzaferis, J. A., Hole, J. T., et al. (2016). TREM2-mediated early microglial response limits diffusion and toxicity of amyloid plaques. J. Exp. Med. 213, 667–675. doi: 10.1084/jem.20151948

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolfe, C., Fitz, N., Nam, K., Lefterov, I., and Koldamova, R. (2018). The role of APOE and TREM2 in Alzheimer’s Disease—current understanding and perspectives. Int. J. Mol. Sci. 20:E81. doi: 10.3390/ijms20010081

PubMed Abstract | CrossRef Full Text | Google Scholar

Wozniak, M. A., Frost, A. L., Preston, C. M., and Itzhaki, R. F. (2011). Antivirals reduce the formation of key Alzheimer’s Disease molecules in cell cultures acutely infected with herpes simplex virus type 1. PLoS One 6:e25152. doi: 10.1371/journal.pone.0025152

PubMed Abstract | CrossRef Full Text | Google Scholar

Wunderlich, P., Glebov, K., Kemmerling, N., Tien, N. T., Neumann, H., and Walter, J. (2013). Sequential proteolytic processing of the triggering receptor expressed on myeloid cells-2 (TREM2) PROTEIN BY ECTODOMAIN SHEDDING AND γ-Secretase-dependent intramembranous cleavage. J. Biol. Chem. 288, 33027–33036. doi: 10.1074/jbc.M113.517540

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, P., Condello, C., Keene, C. D., Wang, Y., Bird, T. D., Paul, S. M., et al. (2016). TREM2 haplodeficiency in mice and humans impairs the microglia barrier function leading to decreased amyloid compaction and severe axonal dystrophy. Neuron 90, 724–739. doi: 10.1016/j.neuron.2016.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, H., Jia, L., Liu, C.-C., Rong, Z., Zhong, L., Yang, L., et al. (2017). TREM2 Promotes Microglial Survival By Activating Wnt/β-Catenin pathway. J. Neurosci. 37, 1772–1784. doi: 10.1523/JNEUROSCI.2459-16.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhong, L., Chen, X.-F., Wang, T., Wang, Z., Liao, C., Wang, Z., et al. (2017). Soluble TREM2 induces inflammatory responses and enhances microglial survival. J. Exp. Med 214, 597–607. doi: 10.1084/jem.20160844

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhong, L., Chen, X.-F., Zhang, Z.-L., Wang, Z., Shi, X.-Z., Xu, K., et al. (2015). DAP12 STABILIZEs the C-terminal fragment of the triggering receptor expressed on myeloid cells-2 (TREM2) and protects against LPS-induced pro-inflammatory response. J. Biol. Chem. 290, 15866–15877. doi: 10.1074/jbc.M115.645986

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhong, L., Xu, Y., Zhuo, R., Wang, T., Wang, K., Huang, R., et al. (2019). Soluble TREM2 ameliorates pathological phenotypes by modulating microglial functions in an Alzheimer’s disease model. Nat. Commun. 10:1365. doi: 10.1038/s41467-019-09118-9119

CrossRef Full Text | Google Scholar