Integrated biology approach reveals molecular and pathological interactions among Alzheimer’s Aβ42, Tau, TREM2, and TYROBP in Drosophila models
- Michiko Sekiya†1,
- Minghui Wang†2, 3,
- Naoki Fujisaki1, 4,
- Yasufumi Sakakibara1,
- Xiuming Quan1,
- Michelle E. Ehrlich2, 5, 6,
- Philip L. De Jager7, 8,
- David A. Bennett9,
- Eric E. Schadt2, 3,
- Sam Gandy5, 10, 11, 12,
- Kanae Ando13,
- Bin Zhang2, 3, 14Email author and
- Koichi M. Iijima1, 4Email author
© The Author(s). 2018
Received: 15 August 2017
Accepted: 12 March 2018
Published: 29 March 2018
Cerebral amyloidosis, neuroinflammation, and tauopathy are key features of Alzheimer’s disease (AD), but interactions among these features remain poorly understood. Our previous multiscale molecular network models of AD revealed TYROBP as a key driver of an immune- and microglia-specific network that was robustly associated with AD pathophysiology. Recent genetic studies of AD further identified pathogenic mutations in both TREM2 and TYROBP.
In this study, we systematically examined molecular and pathological interactions among Aβ, tau, TREM2, and TYROBP by integrating signatures from transgenic Drosophila models of AD and transcriptome-wide gene co-expression networks from two human AD cohorts.
Glial expression of TREM2/TYROBP exacerbated tau-mediated neurodegeneration and synergistically affected pathways underlying late-onset AD pathology, while neuronal Aβ42 and glial TREM2/TYROBP synergistically altered expression of the genes in synaptic function and immune modules in AD.
The comprehensive pathological and molecular data generated through this study strongly validate the causal role of TREM2/TYROBP in driving molecular networks in AD and AD-related phenotypes in flies.
Alzheimer’s disease (AD) is the leading cause of neurodegeneration and dementia. At the level of neuropathology, AD is characterized by aggregation and accumulation of two proteins, β-amyloid peptides (Aβ) and the microtubule-associated protein tau . It is accompanied by the activation of multiple neuroinflammatory pathways . Lines of evidence from laboratories and clinics worldwide support the concept that accumulation of Aβ peptides can be an initiating factor and can lie upstream of tau to drive synaptic dysfunction, neuron death and cognitive impairment [3–7].
A new model was developed to account for the fact that up to one-third of patients with clinically diagnosed AD have no evidence of amyloidosis on brain amyloid imaging . Alternatively, some older individuals with neuropathological AD were asymptomatic during their lifetime . These clinicopathological studies indicate that disease progression is a complex process resulting from the interplay of a number of genetic and environmental factors, some of which modulate accumulation of neuropathology while others modulate synaptic and neuronal resilience . System-level analyses of large datasets from patients have emerged as powerful tools for understanding complex diseases such as AD. Gene expression datasets, along with genomic and clinical information from multiple studies, continue to accumulate and data interpretation is becoming a difficult challenge in these “omics” approaches.
Gene regulatory network analysis is a powerful tool in identifying gene modules pathologically related to human complex diseases including AD [11, 12]. We employed an integrative multiscale network analysis approach to identify key molecular interactions of cellular pathways and causal regulators underlying pathological changes in AD. This approach identified TYROBP (tyrosine kinase binding protein, also known as DAP12), the intracellular adaptor of TREM2 (triggering receptor expressed on myeloid cells 2), as a key driver of immune- and microglia-specific networks that are associated with LOAD pathology [11, 13]. Genome-wide association studies (GWAS) revealed that TREM2, a TYROBP-binding protein, is a risk gene for late-onset sporadic AD [14–16]. More recently, TYROBP-coding sequence genetic variants were found to contribute to an increased risk of early-onset AD . Moreover, TREM2/TYROBP signaling is upregulated by plaque-associated myeloid cells in AD brains and in APP transgenic mice [18–20]. An ectodomain of TREM2 is cleaved and released into the extracellular space as a soluble form (sTREM2) and sTREM2 levels in CSF are elevated in the early symptomatic phase of AD [21–23]. Interestingly, this cleavage of TREM2 is reduced by pathogenic mutations for AD [24, 25]. These reports underscore the role of TREM2/TYROBP in AD pathogenesis.
TREM2 encodes a receptor expressed exclusively in the immune cells in the brain [26, 27]. Studies with TREM2-deficient AD model mice suggest that TREM2 may influence phagocytosis of Aβ-lipid complexes as well as microglial survival  and metabolic fitness . Microglia forms a barrier to restrict Aβ plaque growth and diffusion of soluble Aβ oligomers , thereby ameliorating tau pathology in AD mouse models [31–33]. TREM2 deficiency or the AD-associated R47H mutation in TREM2 significantly reduced accumulation of microglia around Aβ plaques [28, 34]. A recent study also shows that some effects of TREM2 on Aβ pathology may be disease-stage-dependent .
By contrast, activation of microglia can play not only beneficial but also detrimental roles in plaque-related neuropathology. Microglia in adult brains engulf spines and other synaptic processes; exposure to Aβ may inappropriately activate this process to mediate synapse loss . This “synaptophagy” involves complement and CR3, which, like TREM2, can provide an ectodomain protein that interacts with TYROBP . Intriguingly, TYROBP deficiency in APPswe/PS1dE9 mice reduces plaque-associated microglia with improved electrophysiological and learning behavior effects . In addition, TREM2 overexpression failed to improve neuropathology and cognitive impairment in aged APPswe/PS1dE9 mice . More recent studies demonstrate that TREM2 pathway promotes the transition from homeostatic to disease-associated microglia in brains of AD model mice [39, 40].
These reports are consistent with a model wherein TREM2/TYROBP signaling is activated as a protective response against Aβ pathology; however, sustained TREM2/TYROBP activation may ultimately aggravate inflammatory and synapse-related pathologies, thereby driving AD progression. Thus, elucidation of the molecular basis for various interrelationships involving Aβ, TREM2/TYROBP, and tau may fill some gaps in our understanding of AD pathogenesis.
In this study, we aimed to decipher molecular interactions among Aβ, TREM2/TYROBP, and tau by integrating gene expression signatures associated with TREM2/TYROBP from AD Drosophila models and transcriptome-wide gene co-expression networks from two human AD cohorts including Harvard Brain Tissue Resource Center (HBTRC) , and the Religious Orders Study and the Rush Memory and Aging Project (ROSMAP) [41, 42]. The impact of an AD-associated TREM2 R47H variant on these molecular interactions was also analyzed. Our data demonstrate that co-expression of neuronal Aβ42 with glial TREM2R47H/TYROBP led to synergistic downregulation of genes associated with synaptic function modules in fly brains. Moreover, glial expression of both TREM2WT/TYROBP and TREM2R47H/TYROBP exacerbated tau toxicity and synergistically affected the pathways implicated in AD-related neurodegeneration. Thus, gene regulatory networks highlighted by this unbiased, cross species analysis appeared to recapitulate some key features of AD progression and support a key driver role for TREM2/TYROBP in AD pathogenesis.
Flies were maintained in standard cornmeal media at 25 °C. Complementary DNA (cDNA) encoding the full length of TREM2 (NM_018965, RC221132) and TYROBP (NM_198125, RC203771) with Myc-DDK tag were obtained from OriGene Technologies, Inc. These constructs were subcloned into a pJFRC19-13XLexAop2 vector (Addgene #26224). TREM2R47H mutation was introduced by using site-directed mutagenesis kit (Takara Bio Inc.). Transgenic flies were generated by PhiC31 integrase-mediated transgenesis systems (Best Gene Inc.). Transgenic fly lines carrying UAS-Aβ42 and UAS-tau were previously described [43–45]. Repo-LexA (#67096), Elav-GAL4 (#458), GMR-GAL4 (#1104), UAS-para RNAi (#31626), and UAS-mcherry RNAi (#35785) were obtained from the Bloomington Stock Center. For RNA sequencing (RNA-seq), around seven-day-old male flies were used. All experiments were performed using age-matched male flies.
Fly heads for each genotype were homogenized in appropriate buffer and subjected to western blotting or co-immunoprecipitation. Details for western blotting and sequential extractions of Aβ42 were performed as previously described . Anti-FLAG (Sigma-Aldrich), anti-TREM2 (Cell signaling), anti-tau (Millipore), anti-non phospho tau (Merck Millipore), anti-pThr231 tau (Thermo Fisher Scientific), anti-pSer262 tau (Abcam), and anti-tubulin (Sigma-Aldrich) for western blotting were purchased.
Heads of male or female flies were fixed in 4% paraformaldehyde for 24 h at 4 °C and embedded in paraffin. Serial sections (6-μm thickness) through the entire heads were prepared, stained with hematoxylin and eosin (Sigma-Aldrich), and examined by bright-field microscopy. Images of the sections were captured with AxioCam 105 color (Carl Zeiss); the vacuole area was measured using Image J (NIH).
Approximately 25 flies were placed in an empty plastic vial. The vial was then gently tapped to knock all of the flies to the bottom. The numbers of flies in the top, middle, or bottom thirds of the vial were scored after 10 s. The percentages of flies that stayed at the bottom were subjected to statistical analyses. Experiments were repeated more than three times and a representative result was shown.
Courtship-conditioning assay was performed using the method by Ishimoto et al. [46, 47]. Unreceptive, mated-females were prepared as “trainers” one day before the conditioning. For training, a three- to five-day-old virgin male was placed with a trainer female in the courtship chamber (15 mm in diameter × 5 mm in depth) for 1 h. Trained males and non-trained naïve males were tested with freeze-killed virgin females as a courtship target in the courtship chamber 1 h after training. The courtship index (CI) was defined as the proportion of time spent in courtship behavior during 10 min observation period. We used more than 60 flies for each genotype. CIs for conditioned males and naïve controls were analyzed by Mann–Whitney U test. To compare the memory performances of each genotype, experimental data are presented as the performance index (PI), which was calculated using the following formula. PI = 100 × (CI average for naïve – CI conditioned)/CI average for naïve, after CIs were subjected to arcsine square root transformation to approximate normal distributions.
Reverse transcription polymerase chain reaction (RT-PCR) and quantitative reverse transcription polymerase chain reaction (qRT-PCR)
RNA extraction was described below (RNA-seq and analyses). Total RNA was reverse-transcribed using PrimeScript RT-PCR kit (TaKaRa Bio) and the resulting cDNA was used as a template for PCR (Veriti, Applied Biosystems). PCR products were analyzed by 1% agarose gel.
qRT-PCR was performed using PowerSYBR (Thermo Fisher Scientific) on a CFX96 real-time PCR detection system (Bio-Rad Laboratories). The average threshold cycle value (CT) was calculated from at least three replicates per sample. Expression of genes of interest was standardized relative to GAPDH1. Relative expression values were determined by the ΔΔCT method. Primers were designed using Primer-Blast (NIH) or FlyPrimerBank  as described in Additional file 1: Table S11.
RNA sequencing and analyses
More than 100 flies for each genotype were collected and frozen. Heads were mechanically isolated and total RNA was extracted using TRIzol Reagent (Invitrogen, Thermo Fisher Scientific) according to the manufacturer’s protocol with an additional centrifugation step (16,000 × g for 10 min) to remove cuticle membranes before the addition of chloroform. Total RNA was purified using phenol-chloroform reagents after treatment with DNAaseI.
Preparation of samples for RNA-seq analysis was performed using the TruSeq RNA Sample Preparation Kit v2 (Illumina). Briefly, ribosomal RNA was depleted from total RNA using the Ribo-Zero rRNA Removal Kit (Human/Mouse/Rat) (Illumina) to enrich for coding RNA and long non-coding RNA. The cDNA was synthesized using random hexamers, end-repaired, and ligated with appropriate adaptors for sequencing. The library then underwent size selection and purification using AMPure XP beads (Beckman Coulter). The appropriate Illumina recommended 6-bp bar-code bases are introduced at one end of the adaptors during PCR amplification step. The size and concentration of the RNA-seq libraries were measured by Bioanalyzer (Agilent) and Qubit fluorometry (Life Technologies, Thermo Fisher Scientific) before loading onto the sequencer. The Ribo-Zero libraries were sequenced on the Illumina HiSeq 2500 System with 100 nucleotide single-end reads, according to the standard manufacturer’s protocol (Illumina).
Single-ended RNA-seq data were generated with the Illumina HiSeq 2500 platform following the Illumina protocol. The raw sequencing reads were aligned to fly genome BDGP6 using star aligner (version 2.5.0b). Following read alignment, featureCounts  was used to quantify the gene expression at the gene level based on Ensembl gene model. Genes with at least 5 reads in at least one sample were considered expressed and hence retained for further analysis, otherwise removed. The gene level read counts data were normalized using trimmed mean of M-values normalization (TMM) method  to adjust for sequencing library size difference.
Differential gene expression between different genotypes was predicted by linear model analysis using Bioconductor package LIMMA . RNA integrity number (RIN) score was incorporated as a covariate in the linear model to control for sample quality. To adjust for multiple tests, false discovery rate (FDR) of the differential expression test was estimated using the Benjamini–Hochberg (BH) method . Genes with FDR < 0.05 and log2 fold change > 1 or < − 1 were considered significant.
Functional enrichment analysis (FEA)
For FEA of differential expression gene signatures, the gene ontology (GO) annotations were obtained from the flyBase database. Then the enrichment analysis was carried out using Fisher’s exact test (FET), assuming the genes in different sets were identically independently sampled from the genome-wide genes profiled. The BH approach was employed to constrain the FDR.
ROSMAP AD cohort network analysis
We utilized large-scale RNA-seq data of the ROSMAP AD cohort to build a gene co-expression network to capture the coordinated regulation of gene expression traits in brain samples. This dataset profiled gene expression of postmortem brain samples from two longitudinal studies of aging and AD [41, 42]. In both studies, participants enroll without dementia and agree to annual clinical evaluation and organ donation at death. As a result, most decedents are old, without dementia and few participants reach end stage dementia before death. They differ from the types of cases obtained in tertiary care clinics . In the ROSMAP dataset, there were 1059 samples including 362 AD cases and 697 controls. We downloaded preprocessed RNA-seq FPKM gene expression abundance data, SNP genotype data, and DNA methylation data from Synapse (https://doi.org/10.7303/syn3219045). We downloaded preprocessed RNA-seq FPKM gene expression abundance data from Synapse (https://doi.org/10.7303/syn3388564). Genes with at least 1 FPKM in at least 10% of the samples were selected and then the data were corrected for confounding factors including batch, PMI, sex, and RIN score. Co-expression network was constructed by using R package WINA , which implements a computationally optimized procedure for weighted gene co-expression network analysis (WGCNA) . In WINA analysis, we used power β = 6 with other parameters set by default. Thirty-five modules (i.e. clusters of gene showing highly correlated expression profiles across samples) were identified, which were annotated by the mostly enriched gene ontology/canonical pathway term. The modules were rank sorted in relation to AD pathology by multiple sorting features computed from the ROSMAP data, including module-trait correlations and enrichment for AD-related disease gene signatures including DEGs and trait-correlated genes (TCGs) regarding neuropathological/clinical traits such as Braak staging, cognitive score, CERAD neuropathological category, and NIA-Reagan score.
Bayesian causal network was constructed by integrating genome-wide gene expression, SNP genotype, DNA methylation, and known transcription factor (TF)-target relationships. Briefly, we first computed expression quantitative trait loci (eQTLs) and then employed a formal statistical causal inference test (CIT)  to infer the causal probability between gene pairs associated with the same eQTL. With a similar strategy, we also computed causal probability of gene pairs mediated or regulated by a DNA methylation site. The causal relationships inferred were combined with TF-target relationships, and together they were subsequently used as priors for building a causal network through a Monte Carlo Markov Chain (MCMC) simulation-based procedure .
Statistical analysis for biological assay
All results were expressed as mean ± SEM. Unpaired Student’s t-test (Prism7, GraphPad Software Inc.) was used to determine statistical significance as indicated in the figure legends. * indicates p < 0.05, ** indicates p < 0.01 and *** indicates p < 0.001 for biological assays throughout the manuscript.
RNA-seq raw data have been deposited in the Gene Expression Omnibus (GEO) database under accession number GSE99012.
The hypothesis underlying this work was that TREM2/TYROBP plays a causal role in driving molecular networks in AD . To test this hypothesis, we used Drosophila to identify molecular interactions between neuronal expression of Aβ42 or tau and glial expression of TREM2/TYROBP.
Establishment of transgenic flies co-expressing TREM2 (TREM2WT) or TREM2 with pathogenic R47H variant (TREM2R47H) with TYROBP in glial cells
We found that, while the expression of TYROBP proteins was readily detectable by western blotting, TREM2 protein levels were undetectable, raising the possibility that ectopically expressed human TREM2 proteins may be unstable in fly glial cells perhaps because a binding partner that is required to stabilize TREM2 protein was absent (Fig. 2b). Indeed, when TREM2 and TYROBP transgenes were combined and co-expressed in fly glial cells, TREM2 proteins became readily detectable (Fig. 2c). A prior report showed that R47H mutation reduced the stability of TREM2 proteins . We compared protein levels of TREM2WT and TREM2R47H in fly brains and found no significant difference between them (Fig. 2c). When TYROBP was immunoprecipitated from the lysate of bigenic fly brains co-expressing TYROBP and TREM2, TREM2 proteins (both TREM2WT and TREM2R47H) were also precipitated, indicating that TREM2 and TYROBP proteins interact and stabilize each other in fly glial cells (Fig. 2d).
In mammalian cells, TREM2 is cleaved by α-secretase, which results in production of N- and C-terminal fragments of TREM2 . The N-terminal fragments of TREM2 are secreted (sTREM2) and promote inflammatory responses , while C-terminal fragments of TREM2 are further processed by γ-secretase . Western blotting using an anti-TREM2 antibody detected the C-terminal fragment of both TREM2WT and TREM2R47H (Fig. 2e), suggesting that TREM2 is processed and that sTREM2 is produced in fly glial cells in a manner similar to that observed in mammalian cells.
FEA revealed that significant overlap between molecular pathways affected by neuronal expression of Aβ42 and those affected by glial expression of TREM2WT/TYROBP in fly brains
Since TREM2/TYROBP signaling is known to promote survival of microglial cells , we evaluated whether ectopic expression of TREM2 and TYROBP induced any structural changes in the fly brain and/or significantly altered the number of glial cells or neurons. No significant alteration in the size or gross morphology of brain structures was observed in TREM2WT/TYROBP bigenic flies (Additional file 1: Figure S1). In addition, immunostaining of fly brains against a glial marker protein, Repo, or a neuronal marker, Elav, revealed that the numbers of glial cells or neurons were not significantly different between control and TREM2WT/TYROBP bigenic fly brains (Additional file 1: Figure S1). These results suggest that gene expression changes induced by ectopic expression of TREM2WT/TYROBP are not due to either structural defects or altered number of neurons or glial cells in the fly brain.
To identify the biological pathways that are affected by glial expression of TREM2WT/TYROBP, we performed functional enrichment analysis (FEA) for the DEG signatures using GO annotation. The genes upregulated by the expression of TREM2WT/TYROBP were significantly enriched (multiple testing corrected FET p value < 0.05) for pathways designated as “oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen,” “intracellular membrane-bounded organelle,” “electron carrier activity,” “heme binding,” “iron ion binding,” “oxidation-reduction process,” and “cellular response to heat” (Fig. 3b and Additional file 3: Table S2). In contrast, the genes downregulated by the expression of TREM2WT/TYROBP were enriched (corrected FET p < 0.05) in the pathways including “integral component of plasma membrane,” “extracellular region,” “extracellular matrix,” “structural constituent of chitin-based cuticle,” “extracellular space,” “potassium ion transport,” and “myosin light chain kinase activity” (Fig. 3b and Additional file 3: Table S2).
RNA sequence analysis in our Aβ42 fly model identified that neuronal expression of Aβ42 upregulated 437 genes and downregulated 485 genes in heads as compared to control flies (Fig. 3a, Additional file 2: Table S1). The upregulated DEGs were enriched in pathways including “endomembrane system,” “endoplasmic reticulum,” and “oxidation-reduction process” (Fig. 3b and Additional file 3: Table S2). By contrast, the downregulated DEGs were enriched in “electron carrier activity,” “oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen,” “intracellular membrane-bounded organelle,” “heme binding,” “oxidation-reduction process,” “iron ion binding,” “extracellular space,” “transferase activity, transferring phosphorus-containing groups,” “extracellular region,” “metabolic process,” “carboxylic ester hydrolase activity,” “cuticle pigmentation,” and “melanin biosynthetic process” (Fig. 3b and Additional file 3: Table S2).
Interestingly, this analysis revealed that eight of the 15 pathways (13 pathways for DEGs downregulated by Aβ42, three pathways for DEGs upregulated by Aβ42, one pathway is overlapped) enriched for the Aβ42 DEGs were also enriched for the TREM2WT/TYROBP DEGs in the same or opposite direction (Fig. 3b and Additional file 3: Table S2). For example, the pathways “extracellular space” and “extracellular region” were enriched for the genes downregulated by Aβ42 and by TREM2WT/TYROBP, while “oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen,” “oxidation-reduction process,” “electron carrier activity,” “heme binding,” “iron ion binding,” and “intracellular membrane-bounded organelle” were enriched for the genes downregulated by Aβ42 and for the genes upregulated by TREM2WT/TYROBP.
These unbiased analyses revealed that gene expression changes induced by neuronal expression of Aβ42 and glial expression of TREM2WT/TYROBP merged onto the same molecular pathways. However, in fly genome, there is no clear ortholog of either TREM2 or TYROBP. One possibility could be that glial cells sense Aβ42 and/or associated neuronal damages and then induce gene expression changes through endogenous signaling pathways. Ectopically expressed human TREM2/TYROBP may sense these damage-associated signals and impact the overlapping molecular pathways.
TREM2R47H/TYROBP induces gene expression changes similar to those by TREM2WT/TYROBP in the fly brains
TREM2 variants were originally identified as causative mutations in patients with Nasu-Hakola disease . However, recent genetic analysis revealed that R47H variant of TREM2 is associated with a three- to fourfold increased risk for AD [14, 15]. To examine the impact of glial expression of TREM2R47H/TYROBP in the fly brains at the molecular level, we generated RNA-seq data from the head samples in flies with the co-expression of TREM2R47H and TYROBP (TREM2R47H/TYROBP) as described above . Expression of TREM2R47H/TYROBP resulted in 290 upregulated genes and 365 downregulated ones (Fig. 3a and Additional file 2: Table S1). No significant alteration in either the size, gross morphology of brain structures, the numbers of glial cells, or neurons was observed in TREM2R47H /TYROBP bigenic flies (Additional file 1: Figure S1), suggesting that gene expression changes induced by ectopic expression of TREM2R47H/TYROBP is not due to either structural changes or altered number of neurons or glial cells in the fly brain.
At the pathway level, FEA using GO annotation revealed that genes upregulated or downregulated by TREM2R47H/TYROBP were enriched in the same categories as those with TREM2WT/TYROBP (Fig. 3b and Additional file 3: Table S2), although some of the pathways that were enriched in the DEGs in TREM2WT/TYROBP, such as “iron ion binding,” “myosin light chain kinase activity,” and “cellular response to heat,” were not significantly enriched in the DEGs in TREM2R47H/TYROBP (Fig. 3b and Additional file 3: Table S2). The DEG signatures from TREM2WT/TYROBP and TREM2R47H/TYROBP flies shared about half of their members and that 98% of those overlapped genes changed in the same direction (corrected FET p = 2.0 × 10−246, 25.1-fold for downregulated genes; corrected FET p = 1.9 × 10−133, 27.4-fold for upregulated genes) (Fig. 3c).
To quantify differences in gene expression induced by TREM2WT/TYROBP and TREM2R47H/TYROBP, we directly compared mRNA expression levels between these two groups and identified 145 upregulated genes and 157 downregulated genes in TREM2R47H/TYROBP compared to TREM2WT/TYROBP (Fig. 3a). Interestingly, at an FDR of 5%, the upregulated DEGs are significantly enriched for “odorant binding,” “sensory perception of chemical stimulus,” “defense response,” and “response to pheromone,” suggesting that these functional pathways were activated by R47H mutation in TREM2 (Fig. 3b and Additional file 3: Table S2). Among these genes, Drosophila Toll-4 gene (the closest ortholog of human TLR7) detected in the “defense response” module is of particular interest, since TREM2 family proteins are known to modulate Toll-like receptor signaling in mammals [67, 68].
We also compared molecular pathways affected by neuronal expression of Aβ42 and those affected by glial expression of TREM2R47H/TYROBP (Fig. 3a, Additional file 2: Table S1). At the pathway levels, four of the above 15 pathways enriched in the Aβ42 DEGs were enriched in the TREM2R47H/TYROBP DEGs. At the gene level, a significant overlap between Aβ42 DEGs and TREM2 R47H/TYROBP DEGs was observed (Fig. 3c, corrected p value ≤ 1.0 × 10−32, ≥ 7-fold).
Taken all together, biological pathways affected by glial expression of TREM2WT/TYROBP or TREM2R47H/TYROBP in fly heads are similar but about 300 genes show significant difference in mRNA expression. Moreover, glial expression of TREM2R47H/TYROBP impacts several common molecular pathways affected by neuronal expression of Aβ42, though TREM2WT/TYROBP appears to impact many more other common pathways affected by neuronal expression of Aβ42.
Expression of TREM2/TYROBP in glial cells modifies molecular signatures induced by Aβ42 expression in neurons in fly brains
To investigate the effects of TREM2WT/TYROBP on phenotypes as well as gene expression signatures induced by Aβ42, we achieved neuronal expression of Aβ42 and glial expression of TREM2WT/TYROBP in fly brains by using two tissue-specific transgenes expression systems in Drosophila (Additional file 1: Figure S2A). Phenotypic characterization revealed that glial overexpression of TREM2WT/TYROBP did not affect either Aβ42 accumulation levels (Fig. 4a), courtship learning and memory (Fig. 4b), or Aβ42-mediated neurodegeneration (Fig. 4c); however, some exacerbation of Aβ42-mediated behavioral deficits was observed (Fig. 4d and Additional file 1: Figure S2B).
We next analyzed the effects of glial expression of TREM2WT/TYROBP on gene expression signatures induced by neuronal expression of Aβ42 in fly brains. RNA sequence analyses revealed that expression of Aβ42/TREM2WT/TYROBP resulted in upregulation of 533 genes and downregulation of 727 genes compared to control flies (Additional file 1: Figure S2C). Interestingly, comparison of FEA results between Aβ42/TREM2WT/TYROBP and Aβ42 flies revealed that seven of the 15 pathways enriched for the Aβ42 DEGs disappeared when TREM2WT/TYROBP was expressed in glia (Additional file 1: Figure S2D and Additional file 4: Table S3). These pathways include “electron carrier activity,” “oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen,” and “heme binding,” which were enriched in the DEGs downregulated by Aβ42 and in the DEGs upregulated by TREM2WT/TYROBP (Fig. 3b and Additional file 3: Table S2).
Using the same strategy, we also analyzed the effects of glial expression of TREM2R47H/TYROBP on phenotypes as well as gene expression signatures induced by neuronal expression of Aβ42 in fly brains (Additional file 1: Figure S2A). Similar to TREM2WT/TYROBP, glial overexpression of TREM2R47H/TYROBP did not affect Aβ42 accumulation levels (Fig. 4a), courtship learning and memory (Fig. 4b), or Aβ42-mediated neurodegeneration (Fig. 4c); although there is a trend toward subtle exacerbation of Aβ42-mediated behavioral deficits (Fig. 4d and Additional file 1: Figure S2B).
RNA sequence analyses revealed that expression of Aβ42/TREM2R47H/TYROBP resulted in upregulation of 661 genes and downregulation of 846 genes (Additional file 1: Figure S2C and Additional file 2: Table S1). Comparison of FEA results between Aβ42/TREM2R47H/TYROBP and Aβ42 flies revealed that nine of the 15 pathways enriched in the Aβ42 DEGs disappeared following glial expression of TREM2R47H/TYROBP (Additional file 1: Figure S2D and Additional file 4: Table S3). Among these nine pathways, six pathways were also disappeared following glial expression of TREM2WT/TYROBP, suggesting that the effects of TREM2R47H/TYROBP on Aβ42 were similar to that of TREM2WT/TYROBP by this analysis.
Taken together, glial expression of TREM2/TYROBP modifies molecular signatures induced by neuronal expression of Aβ42. Since glial expression of TREM2/TYROBP did not reduce Aβ42 levels (Fig. 4a), the observed changes in FEA results are not simply due to reduced response to Aβ42 in fly brains. In addition, since TREM2/TYROBP proteins are expressed in glial cells and Aβ42 peptides are expressed in neurons (Additional file 1: Figure S2A), this modulatory action likely reflects non-cell autonomous effects by TREM2/TYROBP.
Neuronal Aβ42 and glial TREM2R47H/TYROBP synergistically downregulated genes associated with synaptic and immune function modules of the co-expressed gene networks from human AD brains
To validate these findings, we performed qPCR analyses in five DEGs from the “synaptic transmission” module with known functions related to neuronal activity; Sh, SK, Shab, para, and Nmdar2 (fly orthologs for potassium voltage-gated channel subfamily A, potassium calcium-activated channel subfamily N, potassium voltage-gated channel subfamily B, sodium voltage-gated channel alpha subunits, and glutamate ionotropic receptor NMDA type subunits, respectively). We found that expression levels of Sh were slightly downregulated by Aβ42 expression alone, while those of Sh and Nmdar2 were slightly downregulated by expression of TREM2R47H/TYROBP alone (Fig. 5b). By contrast, expression levels of all five genes were significantly downregulated in Aβ42/TREM2R47H/TYROBP flies.
We further examined whether downregulation of para, a fly ortholog for sodium voltage-gated channel alpha subunits, modifies neuronal dysfunction in Aβ42 flies. Neuronal knockdown of para by RNAi significantly worsened Aβ42-induced locomotor deficits (Fig. 5c). Moreover, neuronal knockdown of para by itself caused modest decline in locomotor functions in flies (Fig. 5d).
Taken altogether, these network analysis results suggest that genes associated with synaptic transmission were synergistically downregulated by co-expression of Aβ42 and TREM2R47H/TYROBP, which may lead to neuronal dysfunction.
To further explore the association of the DEG signatures with early phases of AD, we intersect them with the co-expression network modules from an independent cohort in the ROSMAP study. Again, the modules were annotated by the GO/pathways that the modules were most enriched for. The result is summarized in Fig. 5a and Additional file 5: Table S4b. The genes downregulated by Aβ42/TREM2WT/TYROBP were enriched for an “inflammatory response” module (salmon; corrected FET p = 0.005, 2.6-fold) and an “organic acid metabolism” module (cyan; corrected FET p = 0.014, 2.5-fold). Note that the “inflammatory response” module salmon was ranked number 7 in relation to AD pathology after ranking ROSMAP modules using multiple sorting features, including module-trait correlations and enrichment for genes correlated with or differentially expressed regarding neuropathological and clinical traits Braak staging, global cognition, CERAD neuropathological category (and, by extension, NIA-Reagan score). Importantly, this inflammatory module is not enriched for the DEGs by either Aβ42 alone or TREM2WT/TYROBP alone, suggesting interactions between Aβ42 and TREM2WT/TYROBP at the level of gene expression.
When gene expression levels were compared between Aβ42/TREM2/TYROBP and Aβ42 alone, the downregulated genes in Aβ42/TREM2WT/TYROBP were enriched in the “inflammatory response” module (salmon; corrected FET p = 0.015, 3.5-fold), while downregulated genes in Aβ42/TREM2R47H/TYROBP were enriched in the same “inflammatory response” module (salmon; corrected FET p = 5 × 10−4, 3.4-fold) and a “locomotion” module (yellow; ranked number 9; corrected FET p = 0.034, 1.9-fold).
In summary, for both TREM2WT/TYROBP and TREM2R47H/TYROBP, interaction with Aβ42 affected the “inflammatory response” pathway in flies. This is an interesting observation since neuroinflammation is implicated as a significant contributor to AD pathogenesis and is also consistent with the proposed anti-inflammatory consequences of TREM2 signaling in human microglia.
To test if the enrichment of synaptic transmission and inflammatory response modules was biased by the conserveness of these two pathways between fly and human, we analyzed the enrichment of known GO categories (based on the MSigDB gene sets) in the human orthologs of fly genes. As shown in Additional file 6: Table S5, the mostly enriched gene sets are big pathways including cytoplasm, metabolic process, nucleus, organelle part, and macromolecular complex, which account for 14.6%, 12.2%, 10%, 8.8%, and 7.2% of the 10,938 high confidence human orthologous genes (DIOPT score > 1; http://www.flyrnai.org/cgi-bin/DRSC_orthologs.pl) accordingly, with FDR adjusted FET p value < 1.2E-58. In contrast, the immune system genes only account for < 1.1% of the orthologs and were not enriched (FET p value ≥ 0.54), while synaptic transmission accounted for 1.0% of the orthologs and was marginally enriched (FET p value = 0.004). Thus, it is unlikely that the significant correlation with the inflammatory and synaptic modules in fly signatures were caused by an artifact of overrepresentation of these pathways in the fly-human orthologous genes.
In summary, since Aβ42 accumulation, TREM2/TYROBP activation, altered inflammatory response, and synaptic dysfunctions are all implicated in early phases of AD pathogenesis, Aβ42/TREM2/TYROBP flies may recapitulate some molecular signatures relevant to early stages of AD.
Molecular pathways affected by neuronal expression of tau do not overlap with those affected by glial expression of TREM2/TYROBP in fly brains
In the pathogenesis of AD, abnormal accumulation and toxicity of tau is believed to play a critical role in neurodegeneration. Thus, identification of molecular signatures induced by simultaneous activation of TREM2/TYROBP axis and accumulation of tau may provide important information underlying neurodegenerative process in AD.
We first compared molecular pathways affected by neuronal expression of tau and those affected by glial expression of TREM2/TYROBP in fly heads. We performed RNA sequence analyses and characterized gene expression signatures using an established fly model of human tau toxicity  in which expression of human tau causes progressive degeneration of photoreceptor neurons in the retina . Expression of tau in photoreceptor neurons using GMR-GAL4 driver upregulated 384 genes and downregulated 418 genes in the heads compared to control flies (Additional file 1: Figure S3A and Additional file 2: Table S1). The upregulated DEGs in tau fly heads were associated with “endosome transport via multivesicular body sorting pathway,” “ESCRT III complex,” and “vacuolar transport” (Additional file 1: Figure S3B and Additional file 7: Table S6). In contrast, downregulated DEGs in tau fly heads were significantly enriched in “rhabdomere” and “striated muscle thin filament” (Additional file 1: Figure S3B and Additional file 7: Table S6).
We also analyzed gene expression changes caused by pan-glial expression of TREM2WT/TYROBP or TREM2R47H/TYROBP in the same genetic background carrying the GMR-GAL4 driver. Glial expression of TREM2WT/TYROBP upregulated 448 genes and downregulated 306 genes while TREM2R47H/TYROBP upregulated 475 genes and downregulated 426 genes (Additional file 1: Figure S3A and Additional file 2: Table S1). There were 29–52 genes common between tau DEGs and TREM2/TYROBP DEGs (Additional file 1: Figure S3C, corrected p value ≤ 10−10, ≥ 3.9-fold). However, we observed no pathway that was enriched in both DEG signatures (Additional file 1: Figure S3B and Additional file 7: Table S6).
Taken together, these results suggest that molecular signatures induced by expression of tau are dissimilar to those induced by TREM2/TYROBP in fly heads at the functional pathway level.
Glial expression of TREM2/TYROBP exacerbated tau-mediated neurodegeneration
Next, we examined the effects of glial expression of TREM2/TYROBP on gene expression signatures as well as neurodegenerative phenotypes induced by tau expression. In order to achieve expression of tau in photoreceptor neurons and expression of the TREM2/TYROBP complex in glial cells simultaneously, we utilized two tissue-specific transgenes expression systems in Drosophila (Additional file 1: Figure S4A).
Analysis of the gene regulatory network in AD brains revealed that tau and TREM2/TYROBP synergistically downregulated genes overrepresented in the modules related to immune systems associated with AD pathogenesis
We generated RNA-seq data from tau/TREM2WT/TYROBP and tau/TREM2R47H/TYROBP flies and identified gene expression signatures in comparison with control flies. Expression of tau/TREM2WT/TYROBP upregulated 377 genes and downregulated 476 genes, while expression of tau/TREM2R47H/TYROBP upregulated 596 genes and downregulated 601 genes (Additional file 1: Figure S4B and Additional file 2: Table S1).
Most of the pathways enriched in these DEGs (Additional file 1: Figure S4C and Additional file 8: Table S7) were the same as those detected in either TREM2/TYROBP alone or tau alone (Additional file 1: Figure S3B and Additional file 7: Table S6). However, we observed that “proteolysis” and “UDP-glycosyltransferase activity” were uniquely enriched in the DEGs downregulated by tau/TREM2WT/TYROBP and tau/TREM2R47H/TYROBP, respectively (Additional file 1: Figure S3B and Additional file 8: Table S7). The “proteolysis” pathway contains proteases including angiotensin-converting enzyme (ACE), which have been associated with AD , and UDP-glycosyltransferases, enzymes associated with oligodendrocyte myelination, disruption of which has been implicated in neurodegeneration in AD .
To further explore the relevance of the gene expression signatures in tau, TREM2WT/TYROBP, TREM2R47H/TYROBP, tau/TREM2WT/TYROBP, and tau/TREM2R47H/TYROBP flies to AD, we investigated their association with the 111 co-expressed gene modules derived from co-regulation analyses of brain gene expression in the HBTRC AD and controls . The enrichment analysis shows that no module was enriched for tau DEG signature. Moreover, neither the TREM2WT/TYROBP nor the TREM2R47H/TYROBP DEG signature showed significant enrichment in any of the HBTRC modules.
In the co-expression network from ROSMAP, the DEGs downregulated by tau/TREM2WT/TYROBP and tau/TREM2R47H/TYROBP significantly overlapped with three modules (Fig. 7a, Additional file 9: Table S8b): the “inflammatory response” module (salmon; ranked number 7) (corrected FET p = 9.0 × 10−4, 3.3-fold and corrected FET p = 3.7 × 10−4, 3.0-fold, respectively); “locomotion” (yellow; ranked number 9) (corrected FET p = 4.4 × 10−3, 2.1-fold and corrected FET p = 0.04, 1.7-fold, respectively); and “organic acid metabolism” (cyan; ranked number 24) (corrected FET p = 5.7 × 10−3, 3.1-fold and corrected FET p = 0.02, 2.5-fold, respectively). These three modules also significantly overlapped with the Aβ42-related DEG signatures, as described above (Fig. 5a). The “locomotion” and “organic acid metabolism” modules were also enriched for DEGs downregulated by expression of tau alone (corrected FET p = 0.01, 2.0-fold and corrected FET p = 0.02, 2.8-fold, respectively), while the “inflammatory response” module salmon was enriched in the DEGs downregulated by expression of TREM2WT/TYROBP alone (corrected FET p = 4.1 × 10−5, 4.9-fold). The “inflammatory response” module (or the salmon module, ranked number 7) was also enriched in the DEGs from downregulated by tau/TREM2WT/TYROBP in comparison with tau (corrected FET p = 0.01, 3.6-fold), or by tau/TREM2R47H/TYROBP in comparison with tau (corrected FET p = 0.006, 3.4-fold). In the ROSMAP data (Additional file 10: Table S9), this salmon module had three members downregulated in AD brains, including GADD45A, FABP5, and BAALC-AS1, the first two of which were also downregulated in the present tau/TREM2WT/TYROBP and tau/TREM2R47H/ TYROBP flies (FET p = 9.2 × 10−5, 14.7-fold), consistent with the existence of substantial network consistency when human data and fly data are compared.
Of particular interest, the DEGs upregulated by tau/TREM2WT/TYROBP were enriched for another “inflammatory response” module (lightcyan; ranked number 8; corrected FET p = 0.04, 2.9-fold) in the ROSMAP network. This lightcyan module contained five AD GWAS loci, including CD33, INPP5D, MS4A4A/MS4A6A, RIN3, and TREM2 (Additional file 10: Table S9). Moreover, TYROBP was a member of this ROSMAP lightcyan module. This module was not enriched with DEGs upregulated by either tau alone or TREM2WT/TYROBP alone, suggesting that genetic interactions between tau and TREM2WT/TYROBP may induce this gene expression signature. Moreover, significant enrichment was not observed with the DEG signatures in tau/TREM2R47H/TYROBP flies, suggesting that the TREM2R47H variant may have weaker impact on this module than does TREM2WT.
Taken together, these results revealed that different components of the immune response system were either activated or inhibited by the tau/TREM2/TYROBP pathway. Since both the salmon and lightcyan modules were highly ranked for their predicted relationship to AD pathology, the present results highlighted the importance of inflammatory response subnetworks as potential targets for disease intervention.
As shown in Additional file 1: Figure S5, the salmon and lightcyan modules in the ROSMAP network were adjacent to each other in the cluster dendrogram, indicating that the two inflammatory response modules were highly related in the human data, even though they were regulated differently in tau/TREM2/TYROBP flies. Therefore, our fly models provide valuable biological insights into the human data that were not otherwise evident. To investigate the causal regulatory relationships among the inflammatory response module genes, we combined the genes from the two inflammatory response modules and overlaid the combined gene set onto a Bayesian causal network constructed from the ROSMAP data by using an approach described in our previous study . Figure 7b shows the network structure of the 604 inflammatory response genes of which 270 genes have fly orthologs. Over one-third (96) of the 270 fly orthologs were differentially expressed in at least one of the fly transgenic models analyzed here, resulting in a 1.3-fold enrichment (p value = 2.5 × 10−4). TYROBP was highlighted as a key regulator for controlling a large number of downstream genes in this inflammatory response network: 17, 36, and 79 genes were in the immediate first, second, and third layer downstream of TYROBP, respectively. Therefore, the causal network analysis further validated the causal role of TYROBP and informed other novel key regulators, which modulate the inflammatory response pathways, such as, LAPTM5, MYO1F, CLIC1, and CSF1R, which were highlighted by a large node size in Fig. 7b.
There is an increasing appreciation that immunological mechanisms play important roles in AD pathogenesis, as evidenced by the identification of a number of genes expressed in immune cells of the central nervous system (CNS) carrying genetic variants associated with increased risk for late-onset AD, including CD33 , TREM2 [14, 15], and CR1 . Thus, dysregulation of immune response genes and/or pathways are believed to be key factors in the cause and/or progression of AD. The present transcriptomic analysis indicated an overlap between glial expression of TREM2/TYROBP and neuronal expression of Aβ42. FEA of the DEGs suggested a strong overlap of the common pathways regulated by TREM2WT/TYROBP and by Aβ42. More than half of the pathways detected in Aβ42 DEGs were also detected in TREM2WT/TYROBP DEGs with the same or opposite regulation direction (Fig. 3b). In addition, more than half of the pathways regulated by Aβ42 disappeared by glial expression of TREM2WT/TYROBP (Additional file 1: Figure S2D), suggesting that the changes in the TREM2WT/TYROBP signaling pathway might represent a defense reaction to Aβ42 toxicity. This is consistent with the proposed role of TREM2 as a component of the microglial reaction to Aβ-related pathology [18–20].
In order to identify potential modules associated with AD pathogenesis, we overlaid DEGs onto two independent human AD co-expression networks, one from the HBTRC AD cohort and the other from the ROSMAP AD cohort. Synaptic transmission modules from the HBTRC network and inflammatory response modules from the ROSMAP network were among the most interesting subnetworks enriched with different sets of DEGs (Fig. 5a). Different modules emerged from the two networks, possibly due to differences in the distribution of AD severity within each cohort. The ROSMAP cohort contains normal individuals, and patients with mild cognitive impairment (MCI) to mild to moderate stages of dementia with very few individuals with advanced dementia [41, 42], while the HBTRC AD cohort samples were concentrated in more advanced stages of the disease (CDR 3.0 and higher ).
When TREM2/TYROBP is expressed in glial cells, three out of the five ROSMAP modules enriched with the DEGs identified from fly models of tau toxicity were also enriched with the DEGs detected from fly models of Aβ42 toxicity (Figs. 5a and 7a). This suggests that changes in these pathways maybe part of the pathological interaction between Aβ and tau toxicity and therefore may have implications for elucidation of the pathogenesis of early phases of AD. In turn, identification of molecules that play roles in early phases of AD may point to novel sites of intervention where progression of AD may be slowed or arrested.
We also found that, while the “inflammatory response” module salmon was enriched in downregulated DEGs from both Aβ42/TREM2WT/TYROBP and tau/TREM2WT/TYROBP genotypes, the “inflammatory response” lightcyan module was enriched in upregulated DEGs from the tau/TREM2WT/TYROBP genotype (Fig. 7a). Activation of this subnetwork of the inflammatory response pathway may represent an event linked to late stages of AD characterized by tau toxicity and upregulation of TREM2/TYROBP signaling. Interestingly, the lightcyan module was not detected in the tau/TREM2R47H/TYROBP genotype, suggesting that the pathogenic R47H variant may have a negative impact on activating this inflammatory response. The finding has significant implication for selectively and differentially targeting subnetworks of inflammatory response for potential therapeutic intervention. The lightcyan module contained TYROBP, as well as several AD GWAS gene loci, including CD33, INPP5D, MS4A4A/MS4A6A, RIN3, and TREM2. In addition, this module was highly enriched for various AD signatures and hence was the top ranked module in relation to AD pathology among all ROSMAP modules. Taken together, these results highlighted the lightcyan module as an interesting target for potential disease intervention from both genetic and the molecular pathway perspectives.
Innate immune response is a conserved biological process that multicellular organisms use for their defense against pathogens and toxic stimuli. In AD brains, there is a sustained increase in innate immune activity. In fruit fly, immune response relies on combined action of both cellular processes, such as the phagocytosis of invading microbials, and humoral immune responses, such as the secretion of antimicrobial peptides (AMPs) into the hemolymph . NF-κB signaling pathways play paramount roles in modulating humoral immune response. We noted that DEGs in the present Aβ or Tau flies with or without TREM2/TYROBP showed a significant overlap with the fly Rel/NF-κB perturbation signatures induced by Rel mutation or Rel overexpression in Pal et al.  (Additional file 11: Table S10a). In addition, we found that Rel signatures were enriched for the immune response modules in both HBTRC and ROSMAP data (Additional file 11: Table S10b). For example, the Rel overexpression genes were enriched for the “yellow” (response to biotic stimulus) module in the HBTRC dataset (2.5-fold, BH adjusted FET p value 0.035) and the “salmon” (inflammatory response) module in the ROSMAP dataset (6.4-fold, BH adjusted FET p value 0.001). This suggests that fly is a promising model for studying the NF-κB-controlled immune signaling pathways that are implicated in the Aβ or tau pathologies of AD.
Besides the impact on the immune response modules, we systematically examined the impact of the TREM2 R47H variant on neuronal expression of Aβ42 or tau. Overall, the gene expression changes induced by TREM2 R47H were similar to those induced by TREM2 WT in terms of GO function and co-expression network enrichment under conditions with or without Aβ42 or tau. However, we noted a significant enrichment for the synaptic transmission modules with downregulated DEGs in Aβ42/TREM2R47H/TYROBP flies, but not with Aβ42/TREM2WT/TYROBP flies (Fig. 5a and Additional file 5: Table S4a), suggesting a potential role of R47H variant involved in dysregulation of neuronal activities. Direct comparison of mRNA levels between TREM2 R47H and TREM2 WT under various genotype configurations revealed several consistent GO categories, including upregulation of “odorant binding,” “extracellular region,” “defense response,” and “response to pheromone,” downregulation of “phosphatidate phosphatase activity.” Among these categories, “extracellular region” was consistently identified to differ except under the tau expression background. “Extracellular region” is of particular interest because the R47H variant is located in the extracellular region of TREM2 protein. It is postulated that the amino acid change by this mutation interferes the normal biological function of TREM2, such as the binding to its ligands, its receptor function and its processing by proteases, leading to impaired biological pathways implicated in the pathogenesis of AD . In addition, “defense response” contains Toll-4, a fly ortholog of mammalian TLR7, suggesting that TREM2 R47H variant may have distinct impact on Toll-like receptor singling. We anticipate that the gene signatures and pathways identified in this study will be a starting point for a complete identification of the exact molecular mechanisms underlying how risk for AD is specified by TREM2 R47H .
In summary, we constructed novel transgenic fly models of AD in order to study the genetic interactions between glial expression of TREM2/TYROBP and the neuronal expression of Aβ42 or tau, the two hallmark proteins for the characterization of AD neuropathology. Using these novel transgenic fly models of AD, we also investigated the impact of a TREM2 pathogenic R47H variant (rs75932628), for which the observed effect size has been estimated to be comparable to that of the APOE ε4 allele . To the best of our knowledge, we are the first to systematically analyze phenotypic and genome-wide gene expression changes associated with overexpression of the WT and R47H mutant type TREM2/TYROBP and their interaction with Aβ42- or tau-related pathobiology in vivo. A recent work reports that R47H mutation impairs TREM2-mediated microglial response to Aβ pathology , while our results demonstrate that TREM2 R47H is capable of promoting tau-mediated neurodegeneration. The comprehensive pathological and molecular data generated through this study strongly validate the causal role of TREM2/TYROBP in driving molecular networks in AD and AD-related phenotypes in flies and also provides insight into the role of R47H variant TREM2 in AD pathogenesis.
We thank Dr. T. Awasaki and Dr. H. Ishimoto for their technical advice and the Bloomington Stock Center for fly stocks.
This work was supported by the NIA/NIH grant U01AG046170 (to KMI, BZ, ES, SG, and ME), the Research Funding for Longevity Science from National Center for Geriatrics and Gerontology, Japan, grant number 28-26 (to KMI), Takeda Science Foundation (JP) (to KMI), the NIA/NIH grants RF1AG054014 (to BZ), RF1AG057440 (to BZ), R01AG057907 (to BZ), U01AG052411 (to BZ), U01AG46152 (to DB and PDJ), P30AG10161 (to DB and PDJ), RF1AG15819 (to DB and PDJ), and R01AG36836 (to DB and PDJ). U01AG046170 and U01AG46152 are components of the AMP-AD Target Discovery and Preclinical Validation Project.
Availability of data and materials
RNA-sequencing raw data have been deposited in the Gene Expression Omnibus (GEO) database under accession number GSE99012.
Conceptualization: MS, MW, ES, SG, KA, BZ, and KMI; investigation: MS, MW, NF, YS, XQ, BZ, and KMI; writing (original draft): MS, MW, SG, KA, BZ, and KMI; writing (review and editing): ME, PDJ, DB, and ES; supervision: BZ and KMI; funding acquisition: ME, PDJ, DB, ES, SG, BZ, and KMI. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hardy J, Selkoe DJ. The amyloid hypothesis of Alzheimer’s disease: progress and problems on the road to therapeutics. Science. 2002;297(5580):353–6.View ArticlePubMedGoogle Scholar
- Heneka MT, Carson MJ, El Khoury J, Landreth GE, Brosseron F, Feinstein DL, et al. Neuroinflammation in Alzheimer’s disease. Lancet Neurol. 2015;14(4):388–405.View ArticlePubMedGoogle Scholar
- Lewis J, Dickson DW, Lin WL, Chisholm L, Corral A, Jones G, et al. Enhanced neurofibrillary degeneration in transgenic mice expressing mutant tau and APP. Science. 2001;293(5534):1487–91.View ArticlePubMedGoogle Scholar
- Oddo S, Billings L, Kesslak JP, Cribbs DH, LaFerla FM. Abeta immunotherapy leads to clearance of early, but not late, hyperphosphorylated tau aggregates via the proteasome. Neuron. 2004;43(3):321–32.View ArticlePubMedGoogle Scholar
- Caccamo A, Oddo S, Sugarman MC, Akbari Y, LaFerla FM. Age- and region-dependent alterations in Abeta-degrading enzymes: implications for Abeta-induced disorders. Neurobiol Aging. 2005;26(5):645–54.View ArticlePubMedGoogle Scholar
- Chabrier MA, Blurton-Jones M, Agazaryan AA, Nerhus JL, Martinez-Coria H, LaFerla FM. Soluble abeta promotes wild-type tau pathology in vivo. J Neurosci. 2012;32(48):17345–50.View ArticlePubMed CentralPubMedGoogle Scholar
- Bennett DA, Schneider JA, Wilson RS, Bienias JL, Arnold SE. Neurofibrillary tangles mediate the association of amyloid load with clinical Alzheimer disease and level of cognitive function. Arch Neurol. 2004;61(3):378–84.View ArticlePubMedGoogle Scholar
- Jack CR Jr, Wiste HJ, Weigand SD, Knopman DS, Lowe V, Vemuri P, et al. Amyloid-first and neurodegeneration-first profiles characterize incident amyloid PET positivity. Neurology. 2013;81(20):1732–40.View ArticlePubMed CentralPubMedGoogle Scholar
- Iacono D, Resnick SM, O’Brien R, Zonderman AB, An Y, Pletnikova O, et al. Mild cognitive impairment and asymptomatic Alzheimer disease subjects: equivalent beta-amyloid and tau loads with divergent cognitive outcomes. J Neuropathol Exp Neurol. 2014;73(4):295–304.View ArticlePubMed CentralPubMedGoogle Scholar
- Arnold SE, Louneva N, Cao K, Wang LS, Han LY, Wolk DA, et al. Cellular, synaptic, and biochemical features of resilient cognition in Alzheimer’s disease. Neurobiol Aging. 2013;34(1):157–68.View ArticlePubMedGoogle Scholar
- Zhang B, Gaiteri C, Bodea LG, Wang Z, McElwee J, Podtelezhnikov AA, et al. Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer’s disease. Cell. 2013;153(3):707–20.View ArticlePubMed CentralPubMedGoogle Scholar
- Wang M, Roussos P, McKenzie A, Zhou X, Kajiwara Y, Brennand KJ, et al. Integrative network analysis of nineteen brain regions identifies molecular signatures and networks underlying selective regional vulnerability to Alzheimer’s disease. Genome Med. 2016;8(1):104.View ArticlePubMed CentralPubMedGoogle Scholar
- Forabosco P, Ramasamy A, Trabzuni D, Walker R, Smith C, Bras J, et al. Insights into TREM2 biology by network analysis of human brain gene expression data. Neurobiol Aging. 2013;34(12):2699–714.View ArticlePubMed CentralPubMedGoogle Scholar
- 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(2):117–27.Google Scholar
- 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(2):107–16.View ArticlePubMedGoogle Scholar
- Lill CM, Rengmark A, Pihlstrom L, Fogh I, Shatunov A, Sleiman PM, et al. The role of TREM2 R47H as a risk factor for Alzheimer’s disease, frontotemporal lobar degeneration, amyotrophic lateral sclerosis, and Parkinson’s disease. Alzheimers Dement. 2015;11(12):1407–16.View ArticlePubMed CentralPubMedGoogle Scholar
- Pottier C, Ravenscroft TA, Brown PH, Finch NA, Baker M, Parsons M, et al. TYROBP genetic variants in early-onset Alzheimer’s disease. Neurobiol Aging. 2016;48:2220. e229–2. e215.View ArticleGoogle Scholar
- Frank S, Burbach GJ, Bonin M, Walter M, Streit W, Bechmann I, et al. TREM2 is upregulated in amyloid plaque-associated microglia in aged APP23 transgenic mice. Glia. 2008;56(13):1438–47.View ArticlePubMedGoogle Scholar
- Melchior B, Garcia AE, Hsiung BK, Lo KM, Doose JM, Thrash JC, et al. Dual induction of TREM2 and tolerance-related transcript, Tmem176b, in amyloid transgenic mice: implications for vaccine-based therapies for Alzheimer’s disease. ASN Neuro. 2010;2(3):e00037.View ArticlePubMed CentralPubMedGoogle Scholar
- Jay TR, Miller CM, Cheng PJ, Graham LC, Bemiller S, Broihier ML, et al. TREM2 deficiency eliminates TREM2+ inflammatory macrophages and ameliorates pathology in Alzheimer’s disease mouse models. J Exp Med. 2015;212(3):287–95.View ArticlePubMed CentralPubMedGoogle Scholar
- Suarez-Calvet M, Kleinberger G, Araque Caballero MA, Brendel M, Rominger A, Alcolea D, et al. sTREM2 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. 2016;8(5):466–76.View ArticlePubMed CentralPubMedGoogle Scholar
- Heslegrave A, Heywood W, Paterson R, Magdalinou N, Svensson J, Johansson P, et al. Increased cerebrospinal fluid soluble TREM2 concentration in Alzheimer’s disease. Mol Neurodegener. 2016;11:3.View ArticlePubMed CentralPubMedGoogle Scholar
- Piccio L, Deming Y, Del-Aguila JL, Ghezzi L, Holtzman DM, Fagan AM, et al. Cerebrospinal fluid soluble TREM2 is higher in Alzheimer disease and associated with mutation status. Acta Neuropathol. 2016;131(6):925–33.View ArticlePubMed CentralPubMedGoogle Scholar
- Schlepckow K, Kleinberger G, Fukumori A, Feederle R, Lichtenthaler SF, Steiner H, et al. An Alzheimer-associated TREM2 variant occurs at the ADAM cleavage site and affects shedding and phagocytic function. EMBO Mol Med. 2017;9(10):1356–65.View ArticlePubMed CentralPubMedGoogle Scholar
- Thornton P, Sevalle J, Deery MJ, Fraser G, Zhou Y, Stahl S, et al. TREM2 shedding by cleavage at the H157-S158 bond is accelerated for the Alzheimer’s disease-associated H157Y variant. EMBO Mol Med. 2017;9(10):1366–78.View ArticlePubMed CentralPubMedGoogle Scholar
- Klesney-Tait J, Turnbull IR, Colonna M. The TREM receptor family and signal integration. Nat Immunol. 2006;7(12):1266–73.View ArticlePubMedGoogle Scholar
- Linnartz B, Neumann H. Microglial activatory (immunoreceptor tyrosine-based activation motif)- and inhibitory (immunoreceptor tyrosine-based inhibition motif)-signaling receptors for recognition of the neuronal glycocalyx. Glia. 2013;61(1):37–46.View ArticlePubMedGoogle Scholar
- Wang Y, Cella M, Mallinson K, Ulrich JD, Young KL, Robinette ML, et al. TREM2 lipid sensing sustains the microglial response in an Alzheimer’s disease model. Cell. 2015;160:1061–71.View ArticlePubMed CentralPubMedGoogle Scholar
- Ulland TK, Song WM, Huang SC, Ulrich JD, Sergushichev A, Beatty WL, et al. TREM2 maintains microglial metabolic fitness in Alzheimer’s disease. Cell. 2017;170(4):649–63. e613.View ArticlePubMedGoogle Scholar
- Ulrich JD, Ulland TK, Colonna M, Holtzman DM. Elucidating the role of TREM2 in Alzheimer’s disease. Neuron. 2017;94(2):237–48.View ArticlePubMedGoogle Scholar
- Jiang T, Tan L, Zhu XC, Zhang QQ, Cao L, Tan MS, et al. Upregulation of TREM2 ameliorates neuropathology and rescues spatial cognitive impairment in a transgenic mouse model of Alzheimer’s disease. Neuropsychopharmacology. 2014;39(13):2949–62.View ArticlePubMed CentralPubMedGoogle Scholar
- Yeh FL, Wang Y, Tom I, Gonzalez LC, Sheng M. TREM2 binds to apolipoproteins, including APOE and CLU/APOJ, and thereby facilitates uptake of amyloid-beta by microglia. Neuron. 2016;91(2):328–40.View ArticlePubMedGoogle Scholar
- Jay TR, Hirsch AM, Broihier ML, Miller CM, Neilson LE, Ransohoff RM, et al. Disease progression-dependent effects of TREM2 deficiency in a mouse model of Alzheimer’s disease. J Neurosci. 2017;37(3):637–47.View ArticlePubMed CentralPubMedGoogle Scholar
- Song WM, Joshita S, Zhou Y, Ulland TK, Gilfillan S, Colonna M. Humanized TREM2 mice reveal microglia-intrinsic and -extrinsic effects of R47H polymorphism. J Exp Med. 2018;215:745–60.View ArticlePubMedGoogle Scholar
- Hong S, Beja-Glasser VF, Nfonoyim BM, Frouin A, Li S, Ramakrishnan S, et al. Complement and microglia mediate early synapse loss in Alzheimer mouse models. Science. 2016;352(6286):712–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Squarzoni P, Oller G, Hoeffel G, Pont-Lezica L, Rostaing P, Low D, et al. Microglia modulate wiring of the embryonic forebrain. Cell Rep. 2014;8(5):1271–9.View ArticlePubMedGoogle Scholar
- Haure-Mirande JV, Audrain M, Fanutza T, Kim SH, Klein WL, Glabe C, et al. Deficiency of TYROBP, an adapter protein for TREM2 and CR3 receptors, is neuroprotective in a mouse model of early Alzheimer’s pathology. Acta Neuropathol. 2017;134:769–88.View ArticlePubMed CentralPubMedGoogle Scholar
- Jiang T, Wan Y, Zhang YD, Zhou JS, Gao Q, Zhu XC, et al. TREM2 Overexpression has no improvement on neuropathology and cognitive impairment in aging APPswe/PS1dE9 mice. Mol Neurobiol. 2017;54(2):855–65.View ArticlePubMedGoogle Scholar
- Keren-Shaul H, Spinrad A, Weiner A, Matcovitch-Natan O, Dvir-Szternfeld R, Ulland TK, et al. A unique microglia type associated with restricting development of Alzheimer’s disease. Cell. 2017;169(7):1276–90.View ArticlePubMedGoogle Scholar
- Krasemann S, Madore C, Cialic R, Baufeld C, Calcagno N, El Fatimy R, et al. The TREM2-APOE pathway drives the transcriptional phenotype of dysfunctional microglia in neurodegenerative diseases. Immunity. 2017;47(3):566–81. e569.View ArticlePubMedGoogle Scholar
- Bennett DA, Schneider JA, Arvanitakis Z, Wilson RS. Overview and findings from the religious orders study. Curr Alzheimer Res. 2012;9(6):628–45.View ArticlePubMed CentralPubMedGoogle Scholar
- Bennett DA, Schneider JA, Buchman AS, Barnes LL, Boyle PA, Wilson RS. Overview and findings from the rush Memory and Aging Project. Curr Alzheimer Res. 2012;9(6):646–63.View ArticlePubMed CentralPubMedGoogle Scholar
- Iijima K, Liu HP, Chiang AS, Hearn SA, Konsolaki M, Zhong Y. Dissecting the pathological effects of human Abeta40 and Abeta42 in Drosophila: a potential model for Alzheimer’s disease. Proc Natl Acad Sci U S A. 2004;101(17):6623–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Iijima K, Chiang HC, Hearn SA, Hakker I, Gatt A, Shenton C, et al. Abeta42 mutants with different aggregation profiles induce distinct pathologies in Drosophila. PLoS One. 2008;3(2):e1703.View ArticlePubMed CentralPubMedGoogle Scholar
- Sekiya M, Maruko-Otake A, Hearn S, Sakakibara Y, Fujisaki N, Suzuki E, et al. EDEM function in ERAD protects against chronic ER proteinopathy and age-related physiological decline in drosophila. Dev Cell. 2017;41(6):652–64. e655.View ArticlePubMedGoogle Scholar
- Ishimoto H, Sakai T, Kitamoto T. Ecdysone signaling regulates the formation of long-term courtship memory in adult Drosophila melanogaster. Proc Natl Acad Sci U S A. 2009;106(15):6381–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Ishimoto H, Wang Z, Rao Y, Wu CF, Kitamoto T. A novel role for ecdysone in Drosophila conditioned behavior: linking GPCR-mediated non-canonical steroid action to cAMP signaling in the adult brain. PLoS Genet. 2013;9(10):e1003843.View ArticlePubMed CentralPubMedGoogle Scholar
- Hu Y, Sopko R, Foos M, Kelley C, Flockhart I, Ammeux N, et al. FlyPrimerBank: an online database for Drosophila melanogaster gene expression analysis and knockdown evaluation of RNAi reagents. G3 (Bethesda). 2013;3(9):1607–16.View ArticleGoogle Scholar
- Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.View ArticlePubMedGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.View ArticlePubMedGoogle Scholar
- Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.View ArticlePubMed CentralPubMedGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995;B 57:289–300.Google Scholar
- Schneider JA, Aggarwal NT, Barnes L, Boyle P, Bennett DA. The neuropathology of older persons with and without dementia from community versus clinic cohorts. J Alzheimers Dis. 2009;18(3):691–701.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005;4(1):Article 17.View ArticleGoogle Scholar
- Millstein J, Zhang B, Zhu J, Schadt EE. Disentangling molecular relationships with a causal inference test. BMC Genet. 2009;10(1):23.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhu J, Wiener MC, Zhang C, Fridman A, Minch E, Lum PY, et al. Increasing the power to detect causal associations by combining genotypic and expression data in segregating populations. PLoS Comput Biol. 2007;3(4):e69.View ArticlePubMed CentralPubMedGoogle Scholar
- Lai SL, Lee T. Genetic mosaic with dual binary transcriptional systems in Drosophila. Nat Neurosci. 2006;9(5):703–9.View ArticlePubMedGoogle Scholar
- 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(243):243ra286.View ArticleGoogle Scholar
- 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(46):33027–36.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhong L, Chen XF, Wang T, Wang Z, Liao C, Wang Z, et al. Soluble TREM2 induces inflammatory responses and enhances microglial survival. J Exp Med. 2017;214(3):597–607.PubMed CentralPubMedGoogle Scholar
- Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3.View ArticlePubMedGoogle Scholar
- Gandy S. The role of cerebral amyloid beta accumulation in common forms of Alzheimer disease. J Clin Invest. 2005;115(5):1121–9.PubMed CentralPubMedGoogle Scholar
- Crowther DC, Kinghorn KJ, Miranda E, Page R, Curry JA, Duthie FA, et al. Intraneuronal Abeta, non-amyloid aggregates and neurodegeneration in a Drosophila model of Alzheimer’s disease. Neuroscience. 2005;132(1):123–35.View ArticlePubMedGoogle Scholar
- Schilling S, Zeitschel U, Hoffmann T, Heiser U, Francke M, Kehlen A, et al. Glutaminyl cyclase inhibition attenuates pyroglutamate Abeta and Alzheimer’s disease-like pathology. Nat Med. 2008;14(10):1106–11.View ArticlePubMedGoogle Scholar
- Casas-Tinto S, Zhang Y, Sanchez-Garcia J, Gomez-Velazquez M, Rincon-Limas DE, Fernandez-Funez P. The ER stress factor XBP1s prevents amyloid-beta neurotoxicity. Hum Mol Genet. 2011;20(11):2144–60.View ArticlePubMed CentralPubMedGoogle Scholar
- 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(3):656–62.View ArticlePubMed CentralPubMedGoogle Scholar
- Hamerman JA, Jarjoura JR, Humphrey MB, Nakamura MC, Seaman WE, Lanier LL. Cutting edge: inhibition of TLR and FcR responses in macrophages by triggering receptor expressed on myeloid cells (TREM)-2 and DAP12. J Immunol. 2006;177(4):2051–5.View ArticlePubMedGoogle Scholar
- Turnbull IR, Gilfillan S, Cella M, Aoshi T, Miller M, Piccio L, et al. Cutting edge: TREM-2 attenuates macrophage activation. J Immunol. 2006;177(6):3520–4.View ArticlePubMedGoogle Scholar
- Hu Y, Flockhart I, Vinayagam A, Bergwitz C, Berger B, Perrimon N, et al. An integrative approach to ortholog prediction for disease-focused and other functional studies. BMC Bioinformatics. 2011;12:357.View ArticlePubMed CentralPubMedGoogle Scholar
- Wittmann CW, Wszolek MF, Shulman JM, Salvaterra PM, Lewis J, Hutton M, et al. Tauopathy in Drosophila: neurodegeneration without neurofibrillary tangles. Science. 2001;293(5530):711–4.View ArticlePubMedGoogle Scholar
- Iijima-Ando K, Sekiya M, Maruko-Otake A, Ohtake Y, Suzuki E, Lu B, et al. Loss of axonal mitochondria promotes tau-mediated neurodegeneration and Alzheimer’s disease-related tau phosphorylation via PAR-1. PLoS Genet. 2012;8(8):e1002918.View ArticlePubMed CentralPubMedGoogle Scholar
- Chotard C, Salecker I. Glial cell development and function in the Drosophila visual system. Neuron Glia Biol. 2007;3(1):17–25.View ArticlePubMed CentralPubMedGoogle Scholar
- Leyns CEG, Ulrich JD, Finn MB, Stewart FR, Koscal LJ, Remolina Serrano J, et al. TREM2 deficiency attenuates neuroinflammation and protects against neurodegeneration in a mouse model of tauopathy. Proc Natl Acad Sci U S A. 2017;114(43):11524–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Chou PS, Wu MN, Chou MC, Chien I, Yang YH. Angiotensin-converting enzyme insertion/deletion polymorphism and the longitudinal progression of Alzheimer’s disease. Geriatr Gerontol Int. 2017;17:1544–50.View ArticlePubMedGoogle Scholar
- Griciuc A, Serrano-Pozo A, Parrado AR, Lesinski AN, Asselin CN, Mullin K, et al. Alzheimer’s disease risk gene CD33 inhibits microglial uptake of amyloid beta. Neuron. 2013;78(4):631–43.View ArticlePubMed CentralPubMedGoogle Scholar
- Crehan H, Holton P, Wray S, Pocock J, Guerreiro R, Hardy J. Complement receptor 1 (CR1) and Alzheimer’s disease. Immunobiology. 2012;217(2):244–50.View ArticlePubMedGoogle Scholar
- Lemaitre B, Hoffmann J. The host defense of Drosophila melanogaster. Annu Rev Immunol. 2007;25(1):697–743.View ArticlePubMedGoogle Scholar
- Pal S, Wu J, Wu LP. Microarray analyses reveal distinct roles for Rel proteins in the Drosophila immune response. Dev Comp Immunol. 2008;32(1):50–60.View ArticlePubMedGoogle Scholar
- Park J-S, Ji IJ, Kim D-H, An HJ, Yoon S-Y. The Alzheimer’s disease-associated R47H variant of TREM2 has an altered glycosylation pattern and protein stability. Front Neurosci. 2016;10:618.View ArticlePubMedGoogle Scholar