Skip to main content

TREM2 is associated with increased risk for Alzheimer’s disease in African Americans

Abstract

Background

TREM2 encodes for triggering receptor expressed on myeloid cells 2 and has rare, coding variants that associate with risk for late-onset Alzheimer’s disease (LOAD) in Caucasians of European and North-American origin. This study evaluated the role of TREM2 in LOAD risk in African-American (AA) subjects. We performed exonic sequencing and validation in two independent cohorts of >800 subjects. We selected six coding variants (p.R47H, p.R62H, p.D87N, p.E151K, p.W191X, and p.L211P) for case–control analyses in a total of 906 LOAD cases vs. 2,487 controls.

Results

We identified significant LOAD risk association with p.L211P (p = 0.01, OR = 1.27, 95%CI = 1.05-1.54) and suggestive association with p.W191X (p = 0.08, OR = 1.35, 95%CI = 0.97-1.87). Conditional analysis suggests that p.L211P, which is in linkage disequilibrium with p.W191X, may be the stronger variant of the two, but does not rule out independent contribution of the latter. TREM2 p.L211P resides within the cytoplasmic domain and p.W191X is a stop-gain mutation within the shorter TREM-2V transcript. The coding variants within the extracellular domain of TREM2 previously shown to confer LOAD risk in Caucasians were extremely rare in our AA cohort and did not associate with LOAD risk.

Conclusions

Our findings suggest that TREM2 coding variants also confer LOAD risk in AA, but implicate variants within different regions of the gene than those identified for Caucasian subjects. These results underscore the importance of investigating different ethnic populations for disease risk variant discovery, which may uncover allelic heterogeneity with potentially diverse mechanisms of action.

Background

TREM2 is a transmembrane receptor expressed on microglia, macrophages, monocyte-derived dendritic cells, and osteoclasts [1,2]. In the human brain, TREM2 is primarily expressed on microglia and controls two signaling pathways: regulation of phagocytosis and suppression of inflammatory reactivity. Previous studies have shown that homozygous loss-of-function mutations in TREM2 are associated with the polycystic lipomembranous osteodysplasia with sclerosing leukoencephalopathy (PLOSL; also called Nasu-Hakola disease) characterized by progressive early-onset dementia and bone lesions leading to fractures [3-5]. Association between TREM2 and Alzheimer’s disease (AD) was initially reported by two independent groups, who identified a rare heterozygous missense mutation in TREM2, p.R47H (rs75932628-T), that increased AD risk by about 2–4.5 fold [6,7]. One of these studies also reported additional, rare coding variants in exon 2 of TREM2 that were collectively enriched in AD cases compared to controls [6].

Since then, association between p.R47H and AD risk was replicated in Spanish and French-Caucasian populations [8-10] and other North American-Caucasian series, [11,12] but not in Asian populations, where this variant is either very rare or absent, highlighting the need to study diverse populations [13,14]. Several studies suggest that TREM2-p.R47H could also be a risk factor for Parkinson’s disease, frontotemporal dementia (FTD), and amyotrophic lateral sclerosis [15-18]. These observations clearly highlight two common themes in neurodegenerative diseases: (1) homozygous and heterozygous mutations in TREM2 have pleiotropic effects on clinically distinct disorders and (2) neuroinflammation plays a key role in neurodegenerative disease pathogenesis.

To identify additional AD risk variants in TREM2, several deep re-sequencing studies have been performed to comprehensively screen the coding regions of this gene. Cuyvers et al. sequenced TREM2 in a Belgian population and found coding variants [19], some of which were not previously reported [6]. Although none of the observed variants were significantly associated with AD risk, Cuyvers et al. found an enrichment of TREM2 variants in both AD and FTD patients compared to controls [19]. Recently, our group performed deep re-sequencing of TREM2 coding regions in approximately 4,000 individuals of European descent [20]. We identified sixteen non-synonymous variants, six of which were not identified in previous AD studies [20]. Besides p.R47H, we found that p.R62H was significantly associated with AD risk.

Despite increasing numbers of reports on TREM2 associations in subjects of European descent, to our knowledge, there are no comprehensive studies of TREM2 coding variants in African-Americans (AAs). Furthermore, to date, no deep re-sequencing studies have been conducted to comprehensively catalogue TREM2 rare variants in AAs. Given this knowledge gap and notable differences of AD risk variants between AAs vs. Caucasians [21], assessment of TREM2 in this ethnic group can bring novel insight into the role of this gene in AD. In this study, we hypothesized that coding TREM2 variants also affect AD risk in AAs and sought to catalogue such variants in AAs. To test this hypothesis, we used targeted next-generation sequencing of pooled-DNA and Sanger sequencing of individual DNA samples to search for TREM2 variants in two independent AA cohorts. Follow-up genotyping of six potentially functional variants was conducted in a large late-onset AD (LOAD) case–control series comprised of >3,300 AA subjects from five institutions. Finally, we conducted a joint analysis to evaluate the effect of TREM2 variants on AD risk.

Results

Pooled sequencing in 202 AD cases and 136 controls from WUSM identified eleven non-synonymous TREM2 variants in this AA cohort (Additional file 1: Table S7 and Figure 1A), one of which (p.A105V) has not been identified in previous AD studies [6,19,20]. Sanger sequencing of an independent AA cohort from Mayo, followed by direct genotyping in 179 AD cases and 334 controls identified seven non-synonymous TREM2 variants (Additional file 1: Table S8 and Figure 1A). In these sequenced samples, we identified one AD subject heterozygous for p.T66M. This variant was previously identified in its homozygous state in the proband of a Turkish pedigree with an FTD-like syndrome [16] and also observed in its heterozygous state in a different Caucasian AD subject [6].

Figure 1
figure1

Protein structures and conservation of TREM2 with marked variants. (A) The top panel shows the protein structure of TREM2 (based on ENST00000373113), a type-I transmembrane receptor that is encoded by a gene containing 5 exons. The isoform ENST00000373122 encodes a different protein coding sequence after exon 3 (gradient fill rectangle) compared to ENST00000373113. The soluble form of TREM2 (ENST00000338469) lacks exon 4, which encodes the transmembrane domain, and contains a coding region after exon 5 (texture fill rectangle). Figures shown below include the structure of three different TREM2 transcripts, the location of confirmed variants in the most common TREM2 transcript (ENST00000373113), and the location of confirmed variants only in the sTREM2 transcript (ENST00000338469). Most of the variants in the transmembrane form of TREM2 are located in the extracellular domain with three exceptions, located in the cytoplasmic tail. We identified two variants that are located near the C-terminus of the soluble form of TREM2. (B) The protein conservation analysis of confirmed TREM2 variants. Variants are shown with an arrow identifying the corresponding amino acid position. Protein sequences were downloaded from UniProt. The entries used for each species are as follows: Q9NZC2 (human), Q99NH8 (Mouse), D3ZZ89 (Rat), H2QSZ0 (Chimp), F7CW35 (Frog), Q2YHU4 (Chicken), and E2RP46 (Dog).

Based on the sequencing data from the WUSM and Mayo datasets, we selected six variants for follow-up genotyping in additional cohorts (Table 1): All variants were non-synonymous coding variants that were observed in both cohorts, except p.D87N, which was only observed in WUSM, but included as it was previously shown to be associated with AD risk [6]. Of the selected variants, p.R47H and p.R62H have been identified as risk factors for AD in previous studies [6,7,20]. Although p.T96K was observed in both sequenced cohorts, this was not chosen for follow-up genotyping due to its almost perfect linkage disequilibrium (LD) with p.L211P (D’ = 0.99, R2 = 0.96), which did not have the polyallelic variability problem of the p.T96 locus (p.T96K/M/R).

Table 1 Results of the joint-analysis of 6 follow-up SNPs in TREM2

After completion of direct-genotyping of these variants in an additional 525 LOAD cases and 2,017 controls from Indiana, WHICAP, and Emory, we performed association analyses for TREM2 variants in the combined cohort of 906 AA LOAD cases and 2,487 elderly AA controls. The primary results of the combined analysis of all cohorts are shown in Table 1, and analyses of individual cohorts are depicted in Additional file 1: Tables S7, S8 and Additional file 2: Tables S9, S10 and S11. Surprisingly, we did not find any significant association for the previously reported risk variants p.R47H (p = 0.61, odds ratio [OR] = 1.83 [0.31-10.98]) and p.R62H (p = 1, OR = 0.75 [0.20-2.69]) (Table 1), possibly due to their rarity in AAs. Likewise, the rare variant p.D87N was not significant in this study (p = 1, OR = 1.39 [0.13-15.37]). In contrast, p.L211P in exon 4 of TREM2 is significantly associated with increased risk of AD (p = 0.01, OR = 1.27 [1.05-1.54]; Table 1). Additionally, p.W191X, a stop-gain variant that is only present in the shortest transcript of TREM2, shows a trend toward association with AD risk (p = 0.08, OR = 1.35 [0.97-1.87]; Table 1). We did not observe any significant heterogeneity in the disease association of p.W191X and p.L211P between cohorts (PBreslow-Day = 0.95 and 0.25 for p.W191X and p.L211P respectively; Additional file 3: Figures S1 and Additional file 4: Figure S2). TREM2 p.L211P and p.W191X are in tight LD but only weakly correlated (D’ = 0.98, R2 = 0.25). When we adjusted for p.L211P, the association of p.W191X became insignificant (p = 0.63, OR = 1.10 [0.75-1.61]; Additional file 5: Table S12). Nonetheless, after adjusting for p.W191X, p.L211P remained significantly associated with AD risk (p = 0.05, OR = 1.25 [1.00-1.57]; Additional file 5: Table S12). This suggests that of these two polymorphisms p.L211P is more likely to be the underlying risk variant or the stronger proxy for another functional variant. Finally, TREM2 p.E151K did not achieve significance in the analysis of the combined series (p = 0.22, OR = 2.14 [0.74-6.17], Table 1).

Discussion

This is the first investigation of TREM2 for coding variant discovery, validation and AD risk association in a sizable cohort of AAs (n = 3,393). In this study, we identified 12 TREM2 coding region variants, six of which were genotyped in a combined cohort of 906 AA LOAD cases and 2,487 elderly AA controls. Amongst these six coding variants, two of them, p.L211P and p.W191X, showed nominally significant or suggestive association with risk of LOAD in the combined AA cohort. TREM2 p.L211P was previously observed in two of the three LOAD risk variant discovery studies [6,19,20] performed in Caucasian subjects of European or North American origin. Based on these studies as well as the NHLBI GO Exome Sequencing Project (ESP) Exome Variant Server (EVS) [19], TREM2 p.L211P is a rare variant in the Caucasian population with an allele frequency of ~0.14% in comparison to African-Americans that have a MAF ~ 13%. Likewise, p.W191X, which was observed only one time in the published TREM2 sequencing studies of Caucasians [20] and has a MAF of 0.06% in the European-American (EA) cohort from EVS, has an allele frequency of ~3% in the AA controls in this study and ~4% in the AA LOAD subjects. These two variants are in tight LD with a D’ of 0.98, but R2 of 0.25 due to different allele frequencies. Logistic regression analysis conditioned on either variant revealed a stronger influence of p.L211P, suggesting that this variant is more likely to be functional than p.W191X. That said, p.W191X still has an estimated OR indicative of increased LOAD risk after conditioning on p.L211P, but was no longer significant.

We note that p.L211P is in tight LD with p.T96K, which was not genotyped in all cohorts because this locus is polyallelic (p.T96K/M/R). P.T96K resides within exon 2 where all previously reported TREM2 disease risk variants are localized. Given these and the fact that p.T96K has a predicted damaging functional outcome; we recognize that the AD risk association observed by p.L211P may be driven by this variant. Collectively, our findings are consistent with a model where p.L211P, or p.T96K, is a stronger functional LOAD risk variant, with p.W191X also harboring a smaller but independent effect on LOAD risk. It will also be important to consider an alternative model, where p.L211P is in LD with several functional variants, including p.W191X and/or p.T96K, thus leading to enhanced risk estimates for p.L211P. In this alternative scenario, the other functional risk variants would be non-coding variants that would not be identified in our sequencing paradigm, rarer coding variants not tested in our study or a combination of both. Finally, it should be noted that p.W191X is in LD with the strongest TREM2 region variant observed in a LOAD risk GWAS of AAs (rs7748513: OR ± SE, 1.16 ± 0.05, P = 0.001. LD between p.W191X and s7748513 estimated in the 274 Mayo Clinic samples that were included in the AA GWAS was D’ = 1, R2 = 0.06) [22]. Given that rs7748513 is deep within intron 2 of TREM2, the observed association in the AA GWAS is likely a consequence of tagging functional coding variants.

TREM2 p.L211P resides within exon 4 of ENCODE transcript ID ENST00000373113 (RefSeq NM_018965), which encodes the full length TREM2 protein. Exon 4 is lacking from the shorter TREM2 transcripts (ENST00000373122 which is a predicted protein and ENCODE ID ENST00000338469, a.k.a. RefSeq NM_001271821 or TREM-2V; See Figure 1). ENCODE ID ENST00000338469 encodes a shorter form of TREM2 (sTREM2) that includes the first 143 of the 156 amino acids within the predicted secreted extracellular domain, in addition to exon 5 that partially encodes the cytoplasmic domain and an additional 54 amino acid C-terminal domain not shared with the full length protein. TREM2 p.W191X resides within this unique domain of TREM-2V and leads to a stop-gain mutation resulting in a predicted truncated TREM-2V. Theoretically, both of these coding changes could have consequences on the function of TREM2 via a number of mechanisms.

TREM2 is a type-I transmembrane protein, the extracellular domain of which was shown to undergo proteolytic processing via an ADAM or MMP protein family member [23] and in a more recent report, by ADAM10 [24], followed by intramembranous cleavage of its C-terminal fragment (CTF) by γ-secretase, which also cleaves amyloid precursor protein (APP) (reviewed [25]). Disruption of γ-secretase cleavage of TREM2 leads to accumulation of TREM2-CTFs, which in turn trap the TREM2 signalling adaptor protein DAP12, ultimately leading to reduced signaling [23]. Hence, a conformational change in TREM2, which alters cleavage by γ-secretase, may have downstream functional consequences. It will be important to test whether TREM2 p.L211P which resides within the predicted cytoplasmic domain of the full length TREM2 alters its γ-secretase cleavage.

In a recent in-vitro study of TREM2 coding variants that reside within the extracellular domain, the FTD and FTD-like syndrome risk variants were shown to reduce maturation of TREM2, resulting in reduced cell surface expression, lower sTREM2 release extracellularly and lower phagocytic activity [24]. Thus, whether p.L211P or p.T96K may also alter TREM2 maturation with similar downstream consequences is another potential mechanism that requires evaluation. It is plausible that sTREM2 per se can have functions in signaling [24] or the immune system [26]. Therefore, the effects of p.L211P, p.T96K and p.W191X on sTREM2 levels merit investigation. Finally, p.W191X may reduce sTREM2 protein levels via nonsense-mediated mRNA decay. The missense mutations p.T96K and p.L211P could theoretically result in reduced protein levels via inefficient translation or degradation of full length TREM2, as previously demonstrated for missense mutations in other genes [27].

One caveat of this study is that since we sequenced a relatively small number of individuals and only screened for exon 2 in the Mayo control samples, it is possible that some potential functional variants (e.g. p.P3T, p.R5W, p.F170L, p.A192T, p.D216Y, and p.T223I based on NHLBI Exome Variant Server [EVS] data base) could have been missed in this study. However, since the aforementioned variants are extremely rare (range from 0.02% to 0.11% according to EVS; see Additional file 6: Table S13) in African Americans (AAs), our study lacks the statistical power to reach any meaningful conclusions for their association with AD risk. A larger-scale sequencing study of the TREM2 region in AAs is needed to further determine the underlying variants driving the association with AD.

In our study of AA subjects we did not find a significant LOAD risk association with p.R47H, the most widely replicated TREM2 LOAD risk variant in Caucasians, or with the additional risk variants p.R62H or p.D87N [6,7]. However, p.R47H and p.D87N are more frequent in AA cases compared to controls, as would be expected. The lack of association of these variants may be due to their very low MAF in AAs. All three of these variants are rarer in AA vs. EA subjects (EVS AA and EA MAFs are p.R47H: 0.02% vs. 0.26%, p.R62H 0.27% vs. 1%, and p.D87N 0% vs. 0.14%, respectively). Given this, our cohort was likely underpowered to detect any effects from these variants. Another possibility is that since these TREM2 variants are risk factors rather than deterministic mutations, AA subjects may be lacking additional genetic and/or environmental factors, which exist in EAs and which may potentiate the effects of these variants on LOAD risk. Besides limited power, another weakness of our study is the lack of autopsy-confirmation, which is a universal problem with most LOAD case–control studies, but especially those conducted in non-Caucasian populations.

Conclusions

In summary, our investigation of TREM2 coding variants in AA subjects identified variants which confer LOAD risk, although the most strongly associated variants were not previously demonstrated as risk factors in Caucasian studies because they are much rarer in Caucasian samples. Our study has several strengths including variant discovery through deep-sequencing of two independent cohorts via different approaches, validation by genotyping of a relatively large combined LOAD case–control series from five centers, case–control association analyses and a careful analytic paradigm to discern independent effects of the top two variants. Our findings highlight the importance of investigating different ethnic populations, which may lead to discoveries of novel risk variants due to allelic heterogeneity. The strongest variants identified in our study appear to implicate different regions (cytoplasmic domain) and transcripts (TREM-2V) of TREM2 in LOAD risk than prior reports (exon 2/extracellular domain) in Caucasians. Although the functional consequences of these variants remain to be established, such findings pave the way to inquire alternative disease risk mechanisms and may therefore further our understanding of the role of TREM2 in LOAD and other neurodegenerative diseases.

Methods

In order to maximize the statistical power in detecting significant rare-variant association, we analyzed DNA samples from AA individuals evaluated at Washington University School of Medicine in St. Louis (WUSM; includes samples from Knight-Alzheimer’s Disease Research Center [Knight-ADRC] and NIA-LOAD), Mayo Clinic Jacksonville in Florida (Mayo), Indiana University School of Medicine (Indiana), Washington Heights-Inwood Community Aging Project (WHICAP), and the Emory University (Emory). All of these individuals were self-reported AAs. Characteristics of the individual cohorts are listed in Table 2. The overall study design is depicted in Figure 2. Further details are provided in the Additional file 7: Supplementary methods and Additional file 8: Tables S1-S6.

Table 2 Demographic information of studied African Americans
Figure 2
figure2

Study design for assessing TREM2 association with AD risk in African Americans.

Abbreviations

TREM2:

Triggering receptor expressed on myeloid cells 2

PLOSL:

Polycystic lipomembranous osteodysplasia with sclerosing leukoencephalopathy

AD:

Alzheimer’s disease

FTD:

Frontotemporal dementia

AA:

African-Americans

LOAD:

Late-onset Alzheimer’s disease

SNP:

Single nucleotide polymorphism

WUSM:

Washington University School of Medicine in St. Louis

WHICAP:

Washington Heights-Inwood Community Aging Project

OR:

Odds ratio

LD:

Linkage disequilibrium

EA:

European American

MAF:

Minor allele frequency

ESP:

Exome Sequencing Project

EVS:

Exome Variant Server

GWAS:

Genome-wide association study

ADAM:

A disintegrin and metalloproteinase

MMP:

Matrix metalloproteinase

CTF:

C-terminal fragment

APP:

Amyloic precursor protein

References

  1. 1.

    Bouchon A, Dietrich J, Colonna M. Cutting edge: inflammatory responses can be triggered by TREM-1, a novel receptor expressed on neutrophils and monocytes. J Immunol. 2000;164:4991–5.

    Article  PubMed  CAS  Google Scholar 

  2. 2.

    Bouchon A, Hernandez-Munain C, Cella M, Colonna M. A DAP12-mediated pathway regulates expression of CC chemokine receptor 7 and maturation of human dendritic cells. J Elem Exp Med. 2001;194:1111–22.

    Article  CAS  Google Scholar 

  3. 3.

    Paloneva J, Manninen T, Christman G, Hovanes K, Mandelin J, Adolfsson R, et al. Mutations in two genes encoding different subunits of a receptor signaling complex result in an identical disease phenotype. Am J Hum Genet. 2002;71:656–62.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  4. 4.

    Paloneva J, Mandelin J, Kiialainen A, Bohling T, Prudlo J, Hakola P, et al. DAP12/TREM2 deficiency results in impaired osteoclast differentiation and osteoporotic features. J Elem Exp Med. 2003;198:669–75.

    Article  CAS  Google Scholar 

  5. 5.

    Klunemann HH, Ridha BH, Magy L, Wherrett JR, Hemelsoet DM, Keen RW, et al. The genetic causes of basal ganglia calcification, dementia, and bone cysts: DAP12 and TREM2. Neurology. 2005;64:1502–7.

    Article  PubMed  CAS  Google Scholar 

  6. 6.

    Guerreiro R, Wojtas A, Bras J, Carrasquillo M, Rogaeva E, Majounie E, et al. TREM2 variants in Alzheimer’s disease. N Engl J Med. 2013;368:117–27.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  7. 7.

    Jonsson T, Stefansson H, Steinberg S, Jonsdottir I, Jonsson PV, Snaedal J, et al. Variant of TREM2 associated with the risk of Alzheimer’s disease. N Engl J Med. 2013;368:107–16.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  8. 8.

    Benitez BA, Cooper B, Pastor P, Jin SC, Lorenzo E, Cervantes S, et al. TREM2 is associated with the risk of Alzheimer’s disease in Spanish population. Neurobiol Aging. 2013;34:1711.e15–7.

    Article  CAS  Google Scholar 

  9. 9.

    Pottier C, Wallon D, Rousseau S, Rovelet-Lecrux A, Richard AC, Rollin-Sillaire A, et al. TREM2 R47H variant as a risk factor for early-onset Alzheimer’s disease. J Alzheimers Dis. 2013;35:45–9.

    PubMed  CAS  Google Scholar 

  10. 10.

    Ruiz A, Dols-Icardo O, Bullido MJ, Pastor P, Rodriguez-Rodriguez E, Lopez De Munain A, et al. Assessing the role of the TREM2 p.R47H variant as a risk factor for Alzheimer’s disease and frontotemporal dementia. Neurobiol Aging. 2014;35:444.e1–4.

    Article  CAS  Google Scholar 

  11. 11.

    Gonzalez Murcia JD, Schmutz C, Munger C, Perkes A, Gustin A, Peterson M, et al. Assessment of TREM2 rs75932628 association with Alzheimer’s disease in a population-based sample: the Cache County Study. Neurobiol Aging. 2013;34:2889.e11–3.

    Article  CAS  Google Scholar 

  12. 12.

    Hooli BV, Parrado AR, Mullin K, Yip WK, Liu T, Roehr JT, et al. The rare TREM2 R47H variant exerts only a modest effect on Alzheimer disease risk. Neurology. 2014;83:1353–8.

    Article  PubMed  CAS  Google Scholar 

  13. 13.

    Yu JT, Jiang T, Wang YL, Wang HF, Zhang W, Hu N, et al. Triggering receptor expressed on myeloid cells 2 variant is rare in late-onset Alzheimer’s disease in Han Chinese individuals. Neurobiol Aging. 2014;35:937.e931–3.

    Google Scholar 

  14. 14.

    Miyashita A, Wen Y, Kitamura N, Matsubara E, Kawarabayashi T, Shoji M, et al. Lack of genetic association between TREM2 and late-onset Alzheimer’s disease in a Japanese population. Journal of Alzheimer’s disease: JAD. 2014;41:1031–8.

    PubMed  CAS  Google Scholar 

  15. 15.

    Benitez BA, Cruchaga C, United States-Spain Parkinson’s Disease Research G. TREM2 and neurodegenerative disease. N Engl J Med. 2013;369:1567–8.

    PubMed Central  PubMed  Google Scholar 

  16. 16.

    Guerreiro RJ, Lohmann E, Bras JM, Gibbs JR, Rohrer JD, Gurunlian N, et al. Using exome sequencing to reveal mutations in TREM2 presenting as a frontotemporal dementia-like syndrome without bone involvement. JAMA neurology. 2013;70:78–84.

    Article  PubMed Central  PubMed  Google Scholar 

  17. 17.

    Rayaprolu S, Mullen B, Baker M, Lynch T, Finger E, Seeley WW, et al. TREM2 in neurodegeneration: evidence for association of the p.R47H variant with frontotemporal dementia and Parkinson’s disease. Mol Neurodegener. 2013;8:19.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  18. 18.

    Cady J, Koval ED, Benitez BA, Zaidman C, Jockel-Balsarotti J, Allred P, et al. TREM2 Variant p.R47H as a Risk Factor for Sporadic Amyotrophic Lateral Sclerosis. JAMA neurology. 2014;71:449–53.

    Article  PubMed Central  PubMed  Google Scholar 

  19. 19.

    Cuyvers E, Bettens K, Philtjens S, Van Langenhove T, Gijselinck I, van der Zee J, et al. Investigating the role of rare heterozygous TREM2 variants in Alzheimer’s disease and frontotemporal dementia. Neurobiol Aging. 2014;35:726.e711–29.

    Article  CAS  Google Scholar 

  20. 20.

    Jin SC, Benitez BA, Karch CM, Cooper B, Skorupa T, Carrell D, et al. Coding variants in TREM2 increase risk for Alzheimer’s disease. Hum Mol Genet. 2014;23:5838–46.

    Article  PubMed  Google Scholar 

  21. 21.

    Reitz C, Jun G, Naj A, Rajbhandary R, Vardarajan BN, Wang LS, et al. Variants in the ATP-binding cassette transporter (ABCA7), apolipoprotein E 4, and the risk of late-onset Alzheimer disease in African Americans. JAMA. 2013;309:1483–92.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  22. 22.

    Reitz C, Mayeux R, Alzheimer’s Disease Genetics C. TREM2 and neurodegenerative disease. N Engl J Med. 2013;369:1564–5.

    Article  PubMed Central  PubMed  Google Scholar 

  23. 23.

    Wunderlich P, Glebov K, Kemmerling N, Tien NT, Neumann H, Walter J. Sequential proteolytic processing of the triggering receptor expressed on myeloid cells-2 (TREM2) protein by ectodomain shedding and gamma-secretase-dependent intramembranous cleavage. J Biol Chem. 2013;288:33027–36.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  24. 24.

    Kleinberger G, Yamanishi Y, Suarez-Calvet M, Czirr E, Lohmann E, Cuyvers E, et al. TREM2 mutations implicated in neurodegeneration impair cell surface transport and phagocytosis. Sci Transl Med. 2014;6:243ra286.

    Article  CAS  Google Scholar 

  25. 25.

    Hardy J, Selkoe DJ. The amyloid hypothesis of Alzheimer’s disease: progress and problems on the road to therapeutics. Science. 2002;297:353–6.

    Article  PubMed  CAS  Google Scholar 

  26. 26.

    Piccio L, Buonsanti C, Cella M, Tassi I, Schmidt RE, Fenoglio C, et al. Identification of soluble TREM-2 in the cerebrospinal fluid and its association with multiple sclerosis and CNS inflammation. Brain. 2008;131:3081–91.

    Article  PubMed Central  PubMed  Google Scholar 

  27. 27.

    Shankaran SS, Capell A, Hruscha AT, Fellerer K, Neumann M, Schmid B, et al. Missense mutations in the progranulin gene linked to frontotemporal lobar degeneration with ubiquitin-immunoreactive inclusions reduce progranulin production and secretion. J Biol Chem. 2008;283:1744–53.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgments

We are grateful to all the individuals who participated in the study and the scientists, clinicians, technicians, and administrative staff who helped in the implementation of this study.

We thank the Genome Technology Access Center in the Department of Genetics at Washington University School of Medicine for help with genomic analysis. The Center is partially supported by NCI Cancer Center Support Grant #P30 CA91842 to the Siteman Cancer Center and by ICTS/CTSA Grant # UL1TR000448 from the National Center for Research Resources (NCRR), a component of the National Institutes of Health (NIH), and NIH Roadmap for Medical Research.

Support for this research was provided by grants from the National Institutes of Health (P30 NS069329, R01 AG044546, and R01 AG035083); National Institute on Aging (R01 AG032990 and U01 AG046139 to NET; R01 AG09956 and P30 AG10133 to Indiana University); National Institute of Neurological Disorders and Stroke (R01 NS080820 to NET), Mayo Alzheimer’s Disease Research Center: (P50 AG0016574 to NRG-R and NET); Emory Alzheimer’s Disease Research Center (P50 AG025688); Department of Veterans Affairs (IBX001820A). The recruitment and clinical characterization of research participants at Washington University were supported by NIH P50 AG05681, P01 AG03991, and P01 AG026276. NIALOAD samples were collected under a cooperative agreement grant (U24 AG026395 and R01 AG041797) awarded by the National Institute on Aging. NIALOAD samples were obtained from the National Cell Repository for Alzheimer’s Disease (NCRAD), which receives government support under a cooperative agreement grant (U24 AG21886) awarded by the National Institute on Aging (NIA). African Americans in the WHICAP sample were recruited under a research award R01AG037212.

This project was also generously supported by the Robert and Clarice Smith and Abigail Van Buren Alzheimer’s Disease Research Program (NRG-R). MMC and NET are supported partly by GHR Foundation grants. MMC is also supported by an Alzheimer’s Association MNIRGD award. Additional funding was provided by the Alzheimer Association (NIRG-11-200110). CC was a recipient of a New Investigator Award in Alzheimer’s disease from the American Federation for Aging Research.

The authors would like to thank the NHLBI GO Exome Sequencing Project and its ongoing studies which produced and provided exome variant calls for comparison: the Lung GO Sequencing Project (HL-102923), the WHI Sequencing Project (HL-102924), the Broad GO Sequencing Project (HL-102925), the Seattle GO Sequencing Project (HL-102926) and the Heart GO Sequencing Project (HL-103010).

Author information

Affiliations

Authors

Corresponding authors

Correspondence to Carlos Cruchaga or Nilüfer Ertekin-Taner.

Additional information

Competing interests

Dr. Graff-Radford has multicenter treatment study grants from Lilly, TauRx and consulted for Cytox. Dr. Goate has received compensation as a member of the scientific advisory board of Cognition Therapeutics. She has also consulted for Amgen and received compensation. The authors declare that they have no competing interests.

Authors’ contributions

SCJ, MMC, CC, AMG and NET were involved with the conception, design, analysis and interpretation of data. BAB, TK, DC, DP, SL, SK, MK, CR and TSW performed the experiments. RM, TSW, JJL, AIL, JM, HH, TF, NRG, AMG and NET recruited and maintained the study cohorts. CC, AMG and NET provided overall supervision of the study and acquired funding. All authors contributed to the drafting and critical revision of the manuscript and have given final approval of the version to be published.

Sheng Chih Jin and Minerva M Carrasquillo contributed equally to this work.

Additional files

Additional file 1: Tables S7-S8.

TREM2 variants discovered and validated in African Americans. Tables show the variants discovered and validated in the (Table S7). Washington University in St. Louis and (Table S8). Mayo Clinic Florida cohorts. NA represents not applicable. The Fisher’s exact test was used to calculate p-values and odds-ratio (OR) for each variant using the default commands in PLINK. Multivariate logistic regression was performed to evaluate the association of p.T96K, p.W191X and p.L211P with AD risk adjusting for age, gender, and APOE genotype.

Additional file 2: Tables S9-S11.

TREM2 variants genotyped in African Americans. Tables show the variants genotyped at (Table S9). Indiana University (Table S10). WHICAP, and (Table S11). Emory University. NA represents not applicable. The Fisher’s exact test was used to calculate p-values and odds-ratio (OR) for each variant using the default commands in PLINK. Multivariate logistic regression was performed to evaluate the association of p.W191X and p.L211P with AD risk adjusting for age, gender, and APOE genotype.

Additional file 3: Figure S1.

Forest plot for p.W191X odds ratios across cohorts. Forest plot of multivariate logistic regression results generated using the R package ‘rmeta’.

Additional file 4: Figure S2.

Forest plot for p.L211P odds ratios across cohorts. Forest plot of multivariate logistic regression results generated using the R package ‘rmeta’.

Additional file 5: Table S12.

Conditional analyses for p.W191X and p.L211P. Conditional association testing was performed to evaluate the association of p.W191X and p.L211P with AD risk adjusting for age, gender, APOE genotype, cohort, and the conditioning SNP.

Additional file 6: Table S13.

TREM2 variants that are present in the African Americans according to NHLBI EVS. The amino acid change, exon location, dbSNP reference number, GVS function category, minor allele frequency (MAF) in 2,203 African Americans, and PolyPhen2 functional prediction of TREM2 variants are listed in each row. The amino acid position is annotated based on the longest protein isoform. &Rs2234258 (p.W191X) is only present in the coding region of the shortest isoform of TREM2 (ENST00000338469). Variants observed in the in 202 AD cases and 136 controls sequenced at Washington University, or in the 179 AD cases sequenced at Mayo Clinic, are shown in bold font. EVS does not list 5 of the missense variants (p.T66M, p.D87N, p.T96M, p.A105V, and p.E202D) that were observed in our study sample.

Additional file 7:

Supplementary Methods. Extended methods detailing the design, cohorts, and statistical methods applied.

Additional file 8: Tables S1-S6.

Demographic characteristics of African Americans. Tables detail the demographic characteristics of the various cohorts. Sample size (N), mean, standard deviation and range for age in years, percentage of female subjects and percentage of subjects that carry at least one APOE-ε4 allele for the AD cases and cognitively normal elderly controls from (Table S1). Knight-ADRC, (Table S2). NIA-LOAD, (Table S3). Mayo Clinic, (Table S4). Indiana University, (Table S5). WHICAP, (Table S6). Emory University. Ages are age-at-onset for ADs and age of last evaluation for controls.

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Jin, S.C., Carrasquillo, M.M., Benitez, B.A. et al. TREM2 is associated with increased risk for Alzheimer’s disease in African Americans. Mol Neurodegeneration 10, 19 (2015). https://doi.org/10.1186/s13024-015-0016-9

Download citation

Keywords

  • LOAD
  • African-American
  • TREM2
  • Coding variants
  • Case–control