Abstract
TREM2 is an Alzheimer's disease (AD) risk gene expressed in microglia. To study the role of Trem2 in a mouse model of β-amyloidosis, we compared PS2APP transgenic mice versus PS2APP mice lacking Trem2 (PS2APP;Trem2ko) at ages ranging from 4 to 22 months. Microgliosis was impaired in PS2APP;Trem2ko mice, with Trem2-deficient microglia showing compromised expression of proliferation/Wnt-related genes and marked accumulation of ApoE. Plaque abundance was elevated in PS2APP;Trem2ko females at 6–7 months; but by 12 or 19–22 months of age, it was notably diminished in female and male PS2APP;Trem2ko mice, respectively. Across all ages, plaque morphology was more diffuse in PS2APP;Trem2ko brains, and the Aβ42:Aβ40 ratio was elevated. The amount of soluble, fibrillar Aβ oligomers also increased in PS2APP;Trem2ko hippocampi. Associated with these changes, axonal dystrophy was exacerbated from 6 to 7 months onward in PS2APP;Trem2ko mice, notwithstanding the reduced plaque load at later ages. PS2APP;Trem2ko mice also exhibited more dendritic spine loss around plaque and more neurofilament light chain in CSF. Thus, aggravated neuritic dystrophy is a more consistent outcome of Trem2 deficiency than amyloid plaque load, suggesting that the microglial packing of Aβ into dense plaque is an important neuroprotective activity.
SIGNIFICANCE STATEMENT Genetic studies indicate that TREM2 gene mutations confer increased Alzheimer's disease (AD) risk. We studied the effects of Trem2 deletion in the PS2APP mouse AD model, in which overproduction of Aβ peptide leads to amyloid plaque formation and associated neuritic dystrophy. Interestingly, neuritic dystrophies were intensified in the brains of Trem2-deficient mice, despite these mice displaying reduced plaque accumulation at later ages (12–22 months). Microglial clustering around plaques was impaired, plaques were more diffuse, and the Aβ42:Aβ40 ratio and amount of soluble, fibrillar Aβ oligomers were elevated in Trem2-deficient brains. These results suggest that the Trem2-dependent compaction of Aβ into dense plaques is a protective microglial activity, limiting the exposure of neurons to toxic Aβ species.
- Alzheimer's disease
- amyloid plaque
- microglia
- microgliosis
- neuritic dystrophy
- TREM2
Introduction
Since the discovery of TREM2 (triggering receptor expressed on myeloid cells 2) variants as genetic risk factors for Alzheimer's disease (AD) (Guerreiro et al., 2013; Jonsson et al., 2013), TREM2 biology has become a focal point in research efforts to better understand how the innate immune system impacts AD and other neurodegenerative diseases (Jay et al., 2017b; Ulrich et al., 2017; Yeh et al., 2017). However, whether Trem2 exerts protective or detrimental functions in mouse models of AD-related neuropathology has been rather unclear (Gratuze et al., 2018; Hansen et al., 2018; Ulland and Colonna, 2018).
In transgenic models of cerebral β-amyloidosis, plaque load has been reportedly increased, decreased, or unchanged in mice that lack Trem2, depending on the model, age, and brain region being analyzed (Jay et al., 2015, 2017a; Wang et al., 2015, 2016; Yuan et al., 2016; Parhizkar et al., 2019). Trem2 deletion in amyloidosis models has also been reported to either increase or decrease phosphorylation of the endogenous tau protein (Jay et al., 2015; Wang et al., 2016). Similarly, studies in neurodegeneration models driven by transgenic expression of the human tau protein have suggested disparate roles of Trem2. In the hTau model (Andorfer et al., 2003), Trem2 deletion increased the amounts of tau phosphorylation and aggregation detected (Bemiller et al., 2017). In contrast, in the PS19 model (Yoshiyama et al., 2007), Trem2 deletion had a protective effect, preventing tau-driven synaptic loss and atrophy in the hippocampus and entorhinal cortex, respectively (Leyns et al., 2017).
Transcriptional profiling studies have defined a disease/damage-associated microglial (DAM) activation state that is commonly observed in the brains of neurodegeneration models (Deczkowska et al., 2018; Friedman et al., 2018). The acquisition of the DAM state is Trem2-dependent in mouse models of AD, amyotrophic lateral sclerosis, and demyelinating disease (Poliani et al., 2015; Wang et al., 2015; Keren-Shaul et al., 2017; Krasemann et al., 2017). Some have argued that this state of microglial activation or alarm is fundamentally protective (Keren-Shaul et al., 2017), whereas others have argued that this state is damaging and that returning microglia to their normal “homeostatic” state would be beneficial (Krasemann et al., 2017). It is conceivable that microglial TREM2 activity may be either protective or detrimental, depending on the disease stage and types of pathology present.
To clarify the role of Trem2 in β-amyloid-driven AD models, we studied the effects of Trem2 deletion on microglial activation, plaque accumulation, and neuronal pathology in the PS2APP model across a wide range of ages and in both sexes. PS2APP mice develop amyloid plaque and attendant gliosis pathologies that increase with age, with female mice accumulating the pathology more rapidly than males (Ozmen et al., 2009). Here we report that the effects of Trem2 deficiency on plaque load varied with age and sex, but notably, plaque accumulation was reduced at older ages in both female and male Trem2 KO (PS2APP;Trem2ko) mice compared with age-matched Trem2 WT (PS2APP;Trem2wt) mice. We consistently observed that the Aβ42:Aβ40 ratio was elevated, plaque morphology was more diffuse, and neuritic dystrophy histopathology was more marked in PS2APP;Trem2ko mice, even at older ages when total plaque was reduced. Additional analyses of PS2APP;Trem2ko mice at 12 months of age revealed ApoE-laden microglia, increased levels of soluble fibrillar oligomeric Aβ, and elevated neurofilament-L in the CSF. By RNAseq, we observed that proliferation-related transcripts in PS2APP;Trem2ko microglia were reduced, particularly those encoding certain components and regulators of Wnt-related signaling. Collectively, our data indicate that Trem2-dependent microglial proliferation and activation attenuate the toxic effects of Aβ toward neurons (i.e., the DAM state is mainly protective) and that measurements of neuronal pathology are more informative than plaque load as readouts of microglial modulation in models of β-amyloidosis.
Materials and Methods
Animals.
All animal care and handling procedures were reviewed and approved by the Genentech Institutional Animal Care and Use Committee and were conducted in full compliance with regulatory statutes, Institutional Animal Care and Use Committee policies, and National Institutes of Health guidelines. Animals were housed in specific pathogen-free conditions with 12 h light/12 h dark/d and maintained on regular chow diets. The Trem2tm1(KOMP)Vlcg null allele (C57BL/6N background) was crossed into the PS2APP model (C57BL/6J background). To generate experimental cohorts, all breeding mice were homozygous carriers of the PS2APP transgene and heterozygous carriers of the Trem2-null allele to allow maximal use of littermates between Trem2 WT (Trem2wt) and Trem2 KO (Trem2ko) PS2APP groups. We designed our study to analyze female cohorts at early (4 months), intermediate (6–7 months), and late (12 months) stages of pathology. We also analyzed males at the age of 6–7 months, a common age for us to examine PS2APP histopathology, to check whether any observed effects of Trem2 deletion were sex-specific. The 4 and 6–7 month cohorts were processed together; the 12 and 19–22 month cohorts were separate batches. When we observed less plaque in 12 month PS2APP;Trem2ko females, we decided to also analyze their still aging male counterparts to learn whether that effect was reproducible. (By that time, the numbers of aged males were somewhat depleted, so we included available Trem2het mice to round out the analysis.) For dendritic spine analysis, animals also carrying the Thy1:GFP-M transgene (Jackson ImmunoResearch Laboratories, stock #007788) were used.
Flow cytometry and FACS.
Animals were anesthetized with ketamine/xylazine and transcardially perfused with 30 ml of ice-cold PBS before dissection of cortex + hippocampus. Care was taken to remove the choroid plexus and as much of the meninges as possible before dissociation and sorting. Tissues were dissociated and cell suspensions prepared as described previously (Srinivasan et al., 2016). All steps were performed on ice or at 4°C to prevent artifactual microglial activation.
For quantifying CD45 immunoreactivity and collecting brain-resident myeloid cells by FACS, cell suspensions from 7 PS2APP;Trem2wt (5 females, 2 males) and 7 PS2APP;Trem2ko (6 females, 1 male) mice at 14–15 months of age were stained with the following antibodies in Hibernate-A medium for 20 min at 4°C on a rotator: APC-conjugated anti-CD11b (BD Biosciences, 561690, 1:200), PE/Cy7-conjugated anti-CD45 (BD Biosciences, 552848, 1:500), FITC-conjugated anti-Ly6g (Tonbo Biosciences, 35-5931, 1:200), and PE-conjugated anti-Ccr2 (R&D Systems, FAB5538P, 1:200). Samples were briefly washed and stained with DAPI before FACS sorting. Myeloid cells were selected by gating live (DAPI-negative) cells for CD11b and CD45 immunoreactivity. To avoid the presence of peripheral myeloid cells in the flow analysis and FACS collections, Ccr2+ cells (peripheral monocytes/macrophages) and Ly6g+ cells (neutrophils) were excluded. The total population of brain-resident myeloid cells (defined as CD11b+CD45+Ccr2−Ly6g−), consisting almost entirely of microglia but including CD45high perivascular macrophages (∼2% of total), from each sample was collected in Hibernate-A. Collected cells were pelleted at 5000 rcf for 8 min, and RNA was extracted from cell pellets using QIAGEN RNeasy Micro kits. To measure CD45 immunoreactivity of brain myeloid populations during FACS, CD45low and CD45high gates were defined using a nontransgenic animal as a control to identify microglia in their normal state (CD45low). The same gates were transposed onto cells from PS2APP animals sorted on the same day in the same machine to ensure accuracy and consistency in determining percentages of CD45low and CD45high populations. Generally, animal pairs, including one PS2APP;Trem2wt and one PS2APP;Trem2ko animal, were processed together from perfusion through sorting. Data for nontransgenic animals in Figure 1F include the nontransgenic animals used to define the CD45 gates used during PS2APP cell sorting, as well as additional age-matched control animals from another cohort processed at another time (due to animal availability constraints).
For measuring microglial cells with β-amyloid content, animals at ∼12 months (n = 3 PS2APP;Trem2wt and 2 PS2APP;Trem2ko) were intraperitoneally injected with methoxy-X04 (10 mg/kg) 24 h before tissue collection. Animals were processed as described above, and dissociated cells were stained for CD11b. Cells were also incubated with Calcein-AM (eBioscience, 65-0853-39, 1:1000) just before flow cytometry to label live cells using the 488 nm excitation channel. Live CD11b+ cells were gated as X04+ or X04− using the DAPI excitation channel to determine the percentage of microglia with ingested amyloid content. Transgenic animals injected with PBS were also used as negative controls (data not shown).
RNA sequencing, differential expression, and gene set analysis.
RNA samples from sorted brain myeloid cells of 7 PS2APP;Trem2wt (5 females, 2 males) and 6 PS2APP;Trem2ko (5 females, 1 male) mice (14–15 months old) were selected for sequencing. The concentration of RNA samples was determined using DS-11 spectrophotometer (DeNovix), and the integrity of RNA was determined by 2100 Bioanalyzer (Agilent Technologies). Approximately 1–5 ng of total RNA was used as an input material for the library generation using SMART-seq v4 Ultra Low Input RNA kit (Clontech). Size of the libraries was confirmed using 4200 TapeStation and High Sensitivity D1K screen tape (Agilent Technologies), and their concentration was determined by qPCR-based method using KAPA Library Quantification Kit. The libraries were multiplexed and then sequenced on HiSeq4000 (Illumina) to generate 30 m of single-end 50 bp reads. Sorted cell RNA-Seq data were analyzed as described previously (Srinivasan et al., 2016). Briefly, Illumina adapters, low-quality sequences, and rRNA reads were first discarded. Remaining reads were aligned to the GRCm38 genome with GSNAP aligner (Wu et al., 2016), and reads overlapping each gene were quantified. Normalization was based on the nRPKM method, which is proportional to size factor normalization of DESeq (Love et al., 2014). Differential expression was performed using voom+limma (Law et al., 2014). Raw RNA-Seq data have been deposited to NCBI GEO under accession number GSE140744.
Heat maps in Figure 2A and Fig. 2-3 were generated using gene sets of interest to compare transcriptional responses in multiple datasets. z scores were calculated as follows. First, log2-scale expression matrices were calculated as max(log2(nRPKM), −4). Then, each gene was centered and scaled to give z scores: for a given gene/sample combination, the z score represents distance of nRPKM value in SDs from the mean log2-scale expression value for that gene across all samples within a dataset. Rows (genes) were organized hierarchically using the Euclidean distance function. Columns (sorted microglia samples) were organized by project and genotype.
Gene ontology (GO) query was submitted on the PANTHER Classification System version 14.1 (Mi et al., 2019) at the website (http://pantherdb.org/) using the following inputs: Gene list = genes differentially expressed between Trem2ko and Trem2wt PS2APP microglia (fold change ≥ 2, adjusted p value ≤ 0.05); Organism = Mus musculus; Analysis = statistical overrepresentation test; Annotation set = GO biological process complete; Reference list = Mus musculus whole genome genes; Test type = Fisher's exact; Correction = Calculate false discovery rate (FDR).
For gene set enrichment analyses in Figure 2D and Fig. 2-2), each sample was assigned a gene set score using log2(nRPKM) values for each gene in the set. A sample's gene set score reflected the average difference, for all genes in the set, between that sample's measured log2(nRPKM) value for a given gene and the average log2(nRPKM) value for the same gene across all samples. In cases when one or more samples had no transcripts detected for a given gene, an imputed log2(nRPKM) value was assigned equal to one log2 step below the lowest log2(nRPKM) value detected for that gene in that sample set. We compared gene set scores between genotype groups using two-tailed t tests assuming unequal variance between groups.
Sectioning, histological and immunological staining.
Single-sex cohorts of animals used for histological and biochemical analyses included 4 month females (n = 11 nontransgenic Trem2wt, 14 nontransgenic Trem2ko, 10 PS2APP;Trem2wt, and 12 PS2APP;Trem2ko), 6–7 month females (n = 11 PS2APP;Trem2wt and 13 PS2APP;Trem2ko), 6–7 month males (n = 16 PS2APP;Trem2wt and 14 PS2APP;Trem2ko), 12 month females (n = 15 PS2APP;Trem2wt and 15 PS2APP;Trem2ko), and 19–22 month males (n = 12 PS2APP;Trem2wt, 8 PS2APP;Trem2het, and 7 PS2APP;Trem2ko). Animals were deeply anesthetized with 2.5% tribromoethanol (0.5 ml/25 g body weight) and transcardially perfused with PBS. One brain hemisphere was drop-fixed in 4% PFA for 2 d at 4°C with agitation and then transferred to PBS for histopathological analyses. The other hemisphere was subdissected into cortical and hippocampal tissues that were frozen and stored at −80°C for biochemical assays. Immersion-fixed hemi-brains were cryoprotected, embedded up to 40 per block in a solid matrix, and coronally sectioned at 35 μm (MultiBrain processing by NeuroScience Associates) as previously described (Wang et al., 2011; Kallop et al., 2014). Sheets of sections were stored in cryoprotectant (30% glycerol, 30% ethylene glycol in PBS) at −20°C until use.
Immunohistochemical (IHC) stains for Iba1, CD68, and Gfap were performed at NeuroScience Associates as described previously (Wang et al., 2011), and CD68-stained sections were counterstained with Nissl (0.05% thionine/0.08 m acetate buffer, pH 4.5). Silver stains for amyloid plaque (Campbell-Switzer stain) (Switzer et al., 1993) and neuronal damage/degeneration (amino cupric silver) or “disintegrative degeneration” stain (de Olmos et al., 1994) were also performed at NeuroScience Associates. The bases for these silver stains are given by Switzer (2000) and described at the NeuroScience Associates website (https://www.neuroscienceassociates.com/technologies/staining/). IHC and silver stains spanned a broad rostral-caudal range, including 8–11 sections per animal. Stained slides were returned to Genentech for imaging and quantitation, and unused sections were also returned to Genentech for cryoprotected storage until used for additional stains.
For X-34 stains, sheets were mounted onto slides and completely dried. Slides were incubated with 10 μm X-34 in PBS containing 40% ethanol and 0.02 N NaOH for 10 min, followed by 3 quick washes in PBS, differentiation in 80% ethanol for 1 min, and additional 3 quick PBS washes. After applying ProLong Diamond Antifade Mountant (Thermo Fisher Scientific, P36961), slides were covered with no. 1 coverslips. Two sections per animal were stained, with all cohorts stained and analyzed simultaneously.
For costaining of plaque, microglia, and ApoE or dystrophic axons, sheets encompassing 2 or 3 sections per animal containing regions of the rostral and caudal hippocampus were washed in PBS and then PBS plus Triton X-100 (PBST, 0.1%) and then blocked in PBST (0.3%) with 5% BSA and 5% normal donkey serum, then incubated overnight with primary antibodies diluted in PBST (0.3%) plus 1% BSA at 4°C. Microglia were labeled with rabbit anti-Iba1 (Wako, 019-19741, 1:000) or goat anti-Iba1 (Abcam, ab5076, 1:1000), ApoE with a rabbit monoclonal antibody (Abcam, ab183597, 1:4000), and dystrophic neurites with rat anti-Lamp1 (Abcam, ab25245, 1:2000). Primary antibody incubation was followed by three 10 min washes in PBST, followed by incubation with secondary antibodies for 2 h at room temperature. Donkey anti-rabbit IgG-Alexa555, anti-rat IgG-Alexa647, and anti-goat IgG-Alexa647 (Thermo Fisher Scientific, 1:500) were used as secondary detection reagents. Following the stain, tissue sheets went through three 10 min washes in PBST (0.1%) and three quick washes in PBS. Sheets were mounted onto slides with 0.1% gelatin in PBS and allowed to dry and adhere to the slide at room temperature. To label plaque, slides were then incubated with 10 μm methoxy-X04 in 40% ethanol in PBS for 10 min, washed briefly in PBS, differentiated in 0.2% NaOH in 80% ethanol for 2 min, washed, and then allowed to dry. Slides were coverslipped with added ProLong Gold Antifade Mountant (Thermo Fisher Scientific, P36961). All cohorts were stained and analyzed simultaneously.
Imaging and quantitation of stained sections.
Brain tissue samples processed by NeuroScience Associates were imaged on the SCN400 whole-slide scanning system (Leica Microsystems) at 200× magnification. MATLAB (The MathWorks) running on a high-performance computing cluster was used for all whole-slide image analyses performed in a blinded manner. Quantification of CD68 or Iba1 staining and enlarged dark cluster areas was performed using morphometric-based methods as previously described (Le Pichon et al., 2013; Kallop et al., 2014). The large dark “cluster” of CD68 or Iba1+ cells coincided with the presence of amyloid plaques. Analysis of amino cupric staining was performed using color thresholds and morphological operations. Plaque area was analyzed from slides stained using the Campbell-Switzer method with plaques appearing with a black or amber hue. Multiple color classifiers spanning narrow ranges in RGB and HSV space were created for positive and negative features. Plaques were segmented using these classifiers and applying adaptive thresholding, Euclidean distance transform, morphological operations, and reconstruction. The percentage plaque load, amino cupric, Iba1, CD68, or Gfap positivity for the entire section were calculated by normalizing the positive pixel area to tissue section area and averaged from 8–11 sections/animal. All images, segmentation overlays, and data were reviewed by a pathologist.
Image acquisition of immunofluorescent slides costained for plaque, microglia, and either ApoE or dystrophic neurites was performed at 200× magnification using the Nanozoomer S60 or XR (Hamamatsu) digital whole-slide scanner. Ideal exposure for each channel was determined based on samples with the brightest intensity and set for the whole set of slides to run as a batch. Total tissue area was detected by thresholding on the Iba1 signal and merging and processing of the binary masks by morphological operations. Methoxy-X04, Lamp1, ApoE, and Iba1 staining was analyzed using a top-hat filter and local threshold followed by morphological opening and closing. For Lamp1 and methoxy-X04 staining, shape factor, roundness, and solidity features were used to eliminate elongated objects. In addition, a minimum size of 34 μm2 was applied to exclude small areas of staining. The detected plaques were used as markers in a marker-controlled watershed segmentation to create watershed lines of separation. The plaque mask was then dilated by 17 μm but constrained to be within watershed lines to prevent merging of plaques in close proximity during dilation. Total Lamp1-positive staining was normalized to the whole tissue area. Plaque-associated Lamp1 and Iba1 staining was constrained to be within the mask of dilated area around plaque and was normalized to the same area. Plaque-associated ApoE staining was constrained to be within the mask of plaque + dilated area and normalized to plaque area. Data were averaged from 2 or 3 sections per animal.
For X-34 stains, images were collected with a confocal laser scanning microscope LSM780 (Carl Zeiss) using Zen 2.3 SP1 software (Carl Zeiss). Eleven z-stack images at 1 μm intervals were collected with Plan-Apochromat 20×/0.8 M27, and maximum intensity projection images were created using Zen software. Images were collected and processed blind to genotypes. Image analysis was performed using MATLAB in a blinded fashion on the maximum intensity projection of the confocal z-stack. Control images that did not have X-34 positive staining were used to determine an initial threshold to exclude background. A threshold that is >99.99% of all pixel intensities in the control images was applied to all images to determine an initial segmentation mask. The binary masks were then smoothed out using morphological opening and closing. A minimum size of 9 μm2 was applied to exclude small areas of staining. A threshold corresponding to the 80th and 50th intensity percentile for the pixels within the segmentation mask of all positive images was applied to analyze compact and diffuse area, respectively. Post-threshold morphological operations and size exclusion were performed as described above. The plaque diffuseness index was calculated as follows: (Areadiffuse+compact − Areacompact)/Areadiffuse+compact. For each animal, data were averaged from two sections per animal, with 3 or 4 images per section consisting of two fields from cortex: one field from dorsal subiculum and/or one field from dentate gyrus molecular layer.
For analysis of methoxy-X04 or ApoE colocalization within Iba1+ microglia, images were collected from the costained slides described above for the 12 month female cohort of animals using confocal laser scanning microscope LSM780 (Carl Zeiss) with Zen 2.3 SP1 software (Carl Zeiss). Eleven z-stack images at 1 μm intervals were collected from the cortex with Plan-Apochromat 20×/0.8 M27. To determine colocalization of ApoE and Iba1 staining, or methoxy-X04 and Iba1 staining, we calculated the Mander's colocalization coefficients using the ImageJ plugin JACoP, as described previously (Dunn et al., 2011). The same thresholds were consistently used to identify the Iba1, ApoE, or methoxy-X04 channel across animal samples. Calculations were performed on the entire z-stack of images.
Aβ peptide measurements.
Frozen hippocampal tissues, described above, were homogenized in 10 volumes of TBS (50 mm Tris pH 7.5, 150 mm NaCl) including complete EDTA-free protease inhibitor mixture (Roche Diagnostics) with aprotinin (20 μg/ml) and leupeptin (10 μg/ml) in a QIAGEN TissueLyser II (3 min at 30 Hz). Samples were then centrifuged at 20,000 × g for 20 min at 4°C. Supernatants were collected as the “TBS fraction” and stored at −80°C until analyzed. The pellet was then homogenized in 10 volumes of 5 m guanidine HCl using the TissueLyser II and then placed on a rotisserie for 3 h at room temperature. Samples were diluted 1:10 in a casein buffer (0.25% casein/5 mm EDTA, pH 8.0, in PBS), including aprotinin (20 μg/ml) and leupeptin (10 μg/ml), vortexed and centrifuged at 20,000 × g for 20 min at 4°C. Supernatants were collected as “GuHCl fractions.” Aβ40 and Aβ42 concentrations in mouse hippocampal samples were measured using an ELISA. Briefly, rabbit polyclonal antibody specific for the C terminus of Aβ40 or Aβ42 (Millipore) was coated onto plates, and biotinylated monoclonal anti-Aβ1-16 (Covance, clone 6E10) was used for detection.
For dot blot analyses, ∼10 μg in 1 μl of lysate (TBS soluble fraction of homogenized mouse hippocampus) was blotted onto nitrocellulose membranes (#LC2001, Invitrogen) and incubated for at least 1 h at room temperature to ensure that the blots were dry. The membrane was blocked with Blocking Buffer (MB-070, Rockland Immunochemicals) with added 0.01% Tween 20, for 1 h at room temperature. The membrane was incubated with Amyloid Fibrils OC (Sigma Millipore, AB2286), Oligomer A11 (Thermo Fisher Scientific, AHB0052), 4G8 (BioLegend, 800703), or 6E10 (BioLegend, 803015) primary antibody diluted 1:1000 in Blocking Buffer for 1 h at room temperature. Total protein was normalized with mouse anti-β-actin (Cell Signaling Technology, 8H10D10, 1:10,000) or rabbit anti-GAPDH (Novus Biological, NB300-323, 1:10,000). After primary antibody incubation, membranes were washed 3 times (10 min each) with TBST (50 mm Tris, 0.5 m NaCl, 0.01% Tween 20). The membrane was incubated with secondary antibodies in Blocking Buffer at 1:15,000 dilution (IRDye 800CW donkey anti-rabbit IgG and IRDye 680LT donkey anti-mouse IgG, LI-COR Biosciences, 926-32213, 926-32212, 926-68022, and 926-68023) for 1 h at room temperature. Membrane was washed 3 times (10 min each) in TBST on rocker. Blots were scanned on Odyssey/LICOR scanner for signals followed by image analysis in Image Studio (version 5.2.5, LI-COR Biosciences).
CSF collection and neurofilament light chain (NfL) analysis.
A separate, mixed sex cohort of 12 month PS2APP;Trem2wt (n = 8) and PS2APP;Trem2ko mice (n = 8) was anesthetized, and CSF was collected from the cisterna magna and placed on ice; then blood was collected from terminal cardiac puncture, placed into EDTA collection tubes, and centrifuged at 20,000 × g for 2 min. (The two genotype groups in this analysis were not littermates since we had not collected plasma or CSF from our original cohorts, and we assembled this cohort just for NfL measurements due to a recommendation received during peer review.) Plasma was collected into tubes and stored at −80°C until transfer. CSF samples were diluted 1:10 in 0.1% BSA in TBS and then stored at −80°C until transfer. Plasma and CSF samples were sent to Quanterix for NfL measurements using the Simoa NF-Light Advantage Kit (product 103186). The Simoa assay is a 2 step digital immunoassay, which measures the quantity of NfL in samples using the Simoa HD-1 Analyzer and Single Molecule Array (Simoa) technology.
Two-photon imaging of plaque and dendritic spine measurements.
The somatosensory cortex from PS2APP mice carrying the Thy1:GFP-M transgene and different Trem2 genotypes was imaged ex vivo via 2-photon microscopy. Single-sex cohorts used for this purpose included 6 month PS2APP females (n = 7 for per Trem2 genotype, used for both plaque counts and dendritic spine measurements) and 8 month PS2APP males (n = 6 per Trem2 genotype, used only for plaque counts). At 24 h before brain collection, animals received intraperitoneal injections of methoxy-X04 (10 mg/kg) to label amyloid structures (Klunk et al., 2002). Animals were anesthetized using isoflurane and transcardially perfused with 10 ml PBS followed by 10 ml of 4% PFA + 10% sucrose in PBS, and the collected brains were fixed overnight in 4% PFA + 10% sucrose in PBS at 4°C. After fixation, brains were mounted in agarose and immersed in PBS. Imaging and analysis were performed under blinded conditions.
Apical dendrites and their spines in somatosensory cortex upper layers were imaged en bloc via a two-photon laser-scanning microscope (Ultima In Vivo Multiphoton Microscopy System, Prairie Technologies) using a Ti:sapphire laser (MaiTai DeepSee Spectra Physics, Newport) tuned to 840 nm and a 60× numerical aperture 1.0 immersion objective lens (Olympus) with pixel resolution of 0.1 μm/pixel across a 1024 × 1024 pixel FOV using 1.0 μm steps, with stack depth determined by the slant of the dendritic branch being imaged. For comparison of spine density relative to plaques in PS2APP animals, an FOV containing a dendrite and nearby plaque within 20 μm was considered “near plaque” and an FOV containing only a dendrite with no visible plaque was considered “away from plaque.” To meet the “away from plaque” criteria, we confirmed that no plaque was present in the FOV and at least 100 μm outside of the containing FOV. From each brain, at least five dendrites per condition (near plaque, away from plaque) were imaged. Dendritic spine density and size measurements were generated using custom, semiautomated image analysis routines in MATLAB (The MathWorks). Spine density was estimated as the total number of visible dendritic spines divided by the corresponding length of dendrite. Relative spine volumes were estimated for each detected spine based on the number of corresponding GFP+ pixels in x, y, z dimensions above a local threshold applied as part of an automated image segmentation algorithm. For en bloc plaque measurements, larger-volume stacks were collected using a 20× immersion objective lens across a 1024 × 1024 pixel FOV with 2 μm steps (∼200 μm depth). Plaque density was quantified by a threshold-based MATLAB routine designed to automatically identify methoxy-X04-labeled plaques.
Statistical analysis.
All values are expressed as mean ± SEM. Statistical analysis was performed using the JMP (version 14.2, SAS Institute) or Prism (version 8.3.0 for Mac, GraphPad) software packages. To compare differences between PS2APP;Trem2wt and PS2APP;Trem2ko groups, we performed unpaired t tests. For comparisons of three or more groups, we performed one-way ANOVA followed by Tukey's multiple-comparisons test. The 6–7 month cohort of male and female mice was purposely analyzed separately to determine the effects of Trem2 deficiency in each sex since female mice have accelerated amyloid pathology compared with males, and we did not want an analysis of interactions between sex and Trem2 genotype to be confounded by the age-dependent differences in pathology between males and females.
Results
Plaque-associated microgliosis is impaired in PS2APP;Trem2ko microglia
PS2APP mice express transgenes encoding familial AD mutations in human presenilin 2 (PS2 N141I) and APP (APP K670N/M671L). By 4 months of age, the first deposits of β-amyloid plaque are detected, with age-dependent plaque accumulation occurring faster in females than in males (Ozmen et al., 2009; Kallop et al., 2014). To determine the role of Trem2 in the progression of amyloid disease pathology, we crossed PS2APP mice with Trem2-deficient mice and examined plaque-related phenotypes in single-sex groups at various stages of pathology. We found a stark reduction in Iba1+ microglial clusters in PS2APP;Trem2ko mice at 4, 6–7, 12, and 19–22 months compared with PS2APP;Trem2wt (Fig. 1A,B). In addition, PS2APP;Trem2ko brains showed reduced staining for CD68 marking active microglial lysosomes by IHC (Fig. 1C,D), reduced percentage of CD45high (“activated”) microglia by flow cytometry (Fig. 1E,F), and reduced percentage of methoxy-X04+ (amyloid-containing) microglia by flow cytometry (Fig. 1G,H). Analysis of confocal z-stack images of cortical tissue from 12 month PS2APP;Trem2ko animals also found a significant reduction in the fraction of Iba1 and methoxy-X04 signals that colocalized with each other (Fig. 1-1). We also observed reductions in total Iba1 and Gfap staining at the later ages, indicative of reduced extents of microgliosis and astrogliosis, respectively, in PS2APP;Trem2ko mice (Fig. 1-1). These observations were consistent with reports of Trem2 deletion in other β-amyloidosis models (Jay et al., 2015; Wang et al., 2015, 2016; Yuan et al., 2016; Parhizkar et al., 2019) and suggested that Trem2 deficiency impairs the ability of microglia to engage plaques and phagocytose Aβ filbrils/aggregates.
Figure 1-1
Trem2-dependent induction of the proliferation and neurodegeneration-related gene expression modules
To further characterize the attenuated microglial response to β-amyloid pathology in PS2APP;Trem2ko mice, we FACS-isolated the resident myeloid cell population from the cortex+hippocampus of 14–15 month PS2APP;Trem2wt versus PS2APP;Trem2ko mice and compared their transcriptomic profiles by RNA sequencing (raw RNA-Seq data deposited in NCBI GEO under accession number GSE140744). Although not affording single-cell resolution, our approach provided certain overall advantages (genome-wide analysis, robust detection of low-copy transcripts, and avoidance of artifactual gene expression that occurs during warm-temperature dissociations) compared with other approaches for transcription profiling of Trem2ko microglia in β-amyloid models that used different cell isolation techniques and/or different RNA detection methods, such as microarray, Nanostring, or single-cell RNAseq (Wang et al., 2015; Keren-Shaul et al., 2017; Krasemann et al., 2017; Griciuc et al., 2019).
Applying cutoffs of ≥2-fold change and adjusted p value ≤ 0.05, we observed only 7 transcripts with increased abundance in PS2APP;Trem2ko versus PS2APP;Trem2wt microglia (excluding Treml1, an artifact of the KO cassette insertion) (Kang et al., 2018). In contrast, 144 transcripts (excluding Trem2) showed reduced abundance in PS2APP;Trem2ko compared with PS2APP;Trem2wt microglia using the same cutoffs (Fig. 2A; for extended data table of genome-wide expression values for each sample and summary statistics for differential gene expression, see Fig. 2-1. The majority of these transcripts showed upregulation in microglial expression profiles from the PS2APP model (Friedman et al., 2018) and other models of β-amyloid pathology (Orre et al., 2014; Wang et al., 2015) compared with nontransgenic mice (Fig. 2A). Therefore, their reduced expression in PS2APP;Trem2ko microglia is another manifestation of the impaired microglial response to β-amyloid pathology. The dependence of these transcripts on Trem2 for their induction in PS2APP microglia was approximately concordant with published data from sorted microglial populations from the 5xFAD model (Wang et al., 2015) (Fig. 2A).
Figure 2-1
Figure 2-2
Figure 2-3
Four of the eight most starkly reduced transcripts (in terms of fold change) in PS2APP;Trem2ko microglia are regulators of canonical Wnt signaling or proliferation (Dkk2, Wif1, Ctnna3, and Asb11) (Fig. 2B). Dkk2, Wif1, and Ctnna3 can all negatively regulate Wnt activity (Hsieh et al., 1999; Busby et al., 2004; Gage et al., 2008), whereas Asb11 is important for maintaining progenitor cell activity in multiple cell types (Diks et al., 2006; Tee et al., 2012). Although the lack of induction for negative Wnt regulators might suggest that Wnt-related signaling was enhanced in Trem2ko microglia, another possibility is that Wnt-related signaling was impaired since Dkk2 is also a context-dependent activator of the pathway (Wu et al., 2000; Mao and Niehrs, 2003; Devotta et al., 2018) and since induction of regulators, including Wif1, can occur downstream of active β-catenin as negative feedback (Diep et al., 2004; Boerboom et al., 2006). Supporting this interpretation, the GO knowledge base identified “positive regulation of canonical Wnt signaling pathway” as a biological process overrepresented (fold enrichment = 10.3, FDR = 0.026) among the 144 transcripts with ≥2-fold reduced abundance in PS2APP;Trem2ko microglia (Fig. 2-2A, with six positive factors in Wnt signaling showing reduced expression including Fzd9 (Karasawa et al., 2002), Sulf2 (Lai et al., 2010), Bambi (Lin et al., 2008), Ptk7 (Berger et al., 2017), and Aspm (Buchman et al., 2011) along with Dkk2 (Fig. 2C). We also analyzed recently published microglia RNA-Seq expression profiles from the 5xFAD model (Griciuc et al., 2019) and observed similar Trem2-dependent induction for the nine above-mentioned genes, with the exception of Ptk7 (Fig. 2-2B. The reduced expression of Wnt-related signaling components and regulators in PS2APP;Trem2ko microglia may be consistent with previous reports of coordinated signaling between Trem2 and β-catenin pathways within microglia (Zheng et al., 2017; Zulfiqar and Tanriöver, 2017).
Other GO biological processes implicated as being downregulated in PS2APP;Trem2ko microglia included positive regulation of bone resorption (fold enrichment = 24.6, FDR = 0.030), protein kinase B signaling (fold enrichment = 16.2, FDR = 0.047), negative regulation of tumor necrosis factor production (fold enrichment = 14.1, FDR = 0.026), positive regulation of smooth muscle cell migration (fold enrichment = 13.9, FDR = 0.025), transmembrane receptor protein tyrosine kinase signaling pathway (fold enrichment = 5.9, FDR = 0.0037), and actin cytoskeleton reorganization (fold enrichment = 3.9, FDR = 0.044) (Fig. 2-2A. The Trem2-dependent genes identified in these processes may underlie described roles for Trem2 in osteoclast function (Cella et al., 2003; Paloneva et al., 2003), AKT and mTOR signaling (Ulland et al., 2017), attenuation of proinflammatory macrophage activation (Turnbull et al., 2006), chemotaxis (Mazaheri et al., 2017), and DAP12 signaling (Bouchon et al., 2001).
We recently defined a number of gene expression modules that can be used to characterize the diverse ways that microglia respond to environmental and genetic perturbations (Friedman et al., 2018). In PS2APP compared with nontransgenic microglia, several of these gene sets were upregulated, including the Neurodegeneration-related, Interferon-related, Proliferation, and LPS-related modules, whereas the Microglia and Brain Myeloid modules that typify microglia in their “homeostatic” or “resting” state were modestly but significantly downregulated (Fig. 2D; Fig. 2-3. The entirety of these changes is approximately equivalent to the so-called DAM (Keren-Shaul et al., 2017) or MgND (Krasemann et al., 2017) microglial activation profiles. These findings underscore the utility of these gene modules in characterizing microglial activation states. For instance, even though only 5 of 82 genes in the Proliferation module were upregulated strongly enough in PS2APP microglia to reach genome-wide significance (adjusted p ≤ 0.05), the overall expression of the module was clearly enriched compared with microglia from nontransgenic mice (Fig. 2D; Fig. 2-3.
We next analyzed the degree to which these modular changes in microglial gene expression in PS2APP mice depended on Trem2. Again, despite only 2 of 82 genes (Ccna2 and Aspm) showing significant reduction in transcript abundance after correction for genome-wide analysis, the overall induction of the Proliferation module was compromised in PS2APP;Trem2ko microglia (Fig. 2D; Fig. 2-3). We also observed less induction of the Proliferation module in Trem2ko microglia expression profiles from the 5xFAD model (Griciuc et al., 2019) (Fig. 2-2C). The reduced expression of proliferation-related genes was consistent with the notion mentioned above that Wnt-related signaling was impaired and also corroborated reports from other β-amyloid models that proliferation markers, such as Ki67 or BrdU, were observed less frequently in Trem2-deficient microglia (Wang et al., 2016; Jay et al., 2017a).
As expected, the Neurodegeneration-related gene set, which overlaps with the so-called DAM (Keren-Shaul et al., 2017) and MGnD (Krasemann et al., 2017) genes but relates more specifically to neurodegenerative disease models, was notably impaired in PS2APP;Trem2ko microglia (Fig. 2D; Fig. 2-3). Of the 134 genes in this set, 80 showed upregulation in PS2APP versus normal microglia (Friedman et al., 2018), and approximately half of these showed impaired induction in PS2APP;Trem2ko microglia (adjusted p values ≤ 0.05) (Fig. 2E). Unlike the clear requirement of Trem2 for induction of many genes in the Neurodegeneration-related gene set, the downregulation of the Microglia and Brain Myeloid modules that normally occurs during any CNS challenge (Friedman et al., 2018) was not prevented by Trem2 deletion since expression of these modules was similar in PS2APP;Trem2wt and PS2APP;Trem2ko microglia (Fig. 2D; Fig. 2-3). No effect of Trem2 deletion on the induction of the LPS-related and Interferon-related modules was observed. Overall, our results are similar to previous analysis of the 5xFAD model, in which the microglial induction of many DAM genes showed substantial Trem2 dependence whereas the downregulation of so-called microglial “homeostatic” genes appeared largely Trem2-independent (Keren-Shaul et al., 2017).
Plaque load is reduced in aged PS2APP;Trem2ko mice
In the 5xFAD mouse model, plaque load was reportedly unchanged in 4 month Trem2ko mice but increased in 8 month Trem2ko hippocampus (Wang et al., 2015, 2016; Yuan et al., 2016). In the APPPS1 model, plaque load was reduced in Trem2ko brains at 2 months, reduced or unchanged at 4 months, and increased in the cortex at 8 months (Jay et al., 2015, 2017a).
The reductions in microglial clustering and amyloid engulfment in PS2APP;Trem2ko mice (Fig. 1) suggested that more plaque might accumulate over time in these brains, relative to PS2APP;Trem2wt mice. Indeed, using the Campbell-Switzer silver stain to label amyloid plaque (Campbell et al., 1987; Switzer, 2000), we observed that plaque burden was increased in 6–7 month PS2APP;Trem2ko females, and trending upward but not reaching significance in 6–7 month PS2APP;Trem2ko males (Fig. 3A,B). We observed similar results in distinct cohorts of 6 month females and 8 month males by in vivo labeling of plaque using methoxy-X04 injection, followed by fixation and two-photon imaging of intact somatosensory cortex (Fig. 3-1). We did not observe any effect of Trem2 deletion on the low levels of plaque deposition detected at the earliest stage examined (4 months).
Figure 3-1
We expected to see further exacerbation of amyloid plaque load in Trem2-deficient brains at more advanced ages, but, to our surprise, plaque load was reduced in 12 month PS2APP;Trem2ko females and in 19–22 month PS2APP;Trem2ko males compared with PS2APP;Trem2wt mice (Fig. 3A,B). An intermediate reduction in plaque load was observed in Trem2 heterozygous (PS2APP;Trem2het) mice at 19–22 months (the only age where heterozygous mice were analyzed) (Fig. 3A). At this age, the Campbell-Switzer stain also revealed a distinctive pattern of “amber core” amyloid staining in PS2APP brains (Fig. 3C), thought to represent a mature form of highly condensed plaque. These amber cores were reduced in PS2APP;Trem2het and nearly absent in PS2APP;Trem2ko brains, indicating that Trem2-dependent microglial activity is essential for the formation of these particular amyloid structures.
Overall, our results are reminiscent of a recent study of APPPS1 mice, which showed that Trem2 deletion produced increased seeding of amyloid plaques at early ages but slower rates of amyloid plaque accumulation at later ages (Parhizkar et al., 2019).
Reduced plaque consolidation, elevated neurotoxic Aβ species, and ApoE-laden microglia in Trem2-deficient brains
In contrast to the effects of Trem2 deletion on total plaque burden that varied with age or sex, we observed consistent changes in plaque compaction and composition across ages and sexes. We used the X-34 stain and confocal microscopy to visualize plaque morphology. Although we were blinded to Trem2 genotype, there was an obvious difference in plaque appearance between PS2APP;Trem2wt and PS2APP;Trem2ko brains, with X-34+ structures in PS2APP;Trem2ko brains looking more splayed and less compact (Fig. 4A), similar to descriptions of Trem2-dependent plaque alterations in other β-amyloid models (Wang et al., 2016; Yuan et al., 2016). Using an algorithm based on X-34 signal intensity to quantify the degree of plaque diffuseness, we observed that plaque morphologies were significantly more diffuse in PS2APP;Trem2ko brains in both sexes at all ages tested (Fig. 4B).
Because stains, such as Campbell-Switzer, methoxy-X04, and X-34, only label fibrillar amyloid structures, we also measured the total abundance of Aβ40 and Aβ42 peptides in soluble (TBS) and insoluble (GuHCl) hippocampal fractions by ELISA. The abundance of Aβ peptides rose markedly with age while the Aβ42:Aβ40 ratio declined, particularly in the GuHCl fraction (Fig. 5). Notably, the Aβ42:Aβ40 ratio was higher in PS2APP;Trem2ko than in PS2APP;Trem2wt brains, in both TBS and GuHCl fractions across ages (Fig. 5A,B). The elevated Aβ42:Aβ40 ratio in PS2APP;Trem2ko brains resulted more from reduced abundance of Aβ40 than from increased abundance of Aβ42 (Fig. 5-1), although Aβ42 abundance was elevated in the 6–7 month females, coinciding with the increased plaque deposition we observed in that group (Fig. 3A). Together with our observation that total amyloid plaque is reduced in PS2APP;Trem2ko brains at later ages (Fig. 3A), these results suggest that the elevated Aβ42:Aβ40 ratio in Trem2-deficient brains may increase plaque seeding at younger ages (since Aβ42 is more prone to aggregate and deposit than Aβ40) (Klein, 2002) while reducing the incorporation and compaction of Aβ into existing plaques at older ages (since Aβ40 permeates dense core structures more readily than Aβ42) (Condello et al., 2015). Thus, Trem2 may both restrict the initial seeding of plaques and promote sequestration and compaction of Aβ into existing plaques.
Figure 5-1
To determine whether the elevated Aβ42:Aβ40 ratio was accompanied by altered abundance of soluble, fibrillar Aβ oligomers, we performed nondenaturing dot blots of hippocampal TBS homogenate supernatants from 12 month females. Using the conformation-specific OC antibody (Tomic et al., 2009), we detected significantly higher levels of soluble fibrillar Aβ oligomers in the soluble fraction from PS2APP;Trem2ko mice (Fig. 5C,D), and we observed a positive correlation between the amount of OC+ fibrillar oligomers in this fraction and the Aβ42:Aβ40 ratio (Fig. 5E). In contrast, when we used the pan-reactive 6E10 Aβ antibody to detect total Aβ species in this fraction, the abundance trended slightly downward in PS2APP;Trem2ko mice and correlated negatively with Aβ42:Aβ40 ratio (Fig. 5F–H). Similar respective trends were also seen when staining with the prefibrillar Aβ oligomer antibody A11 and pan-reactive 4G8 Aβ antibody (data not shown). As a control, we spotted hippocampal TBS homogenate supernatant from a PS2APP;Bace1ko animal (Meilandt et al., 2019), in which soluble APP is still produced by α-cleavage while β-cleavage and thus Aβ production are prevented, and demonstrated that the OC antibody had minimal detection, whereas 6E10 still had substantial signal (compare # symbols in Fig. 5C,F and dotted lines in Fig. 5D,G from a PS2APP;Bace1ko mouse), consistent with the ability of 6E10 to detect both soluble APP and Aβ peptides. These results suggest that, in the absence of Trem2, the increased Aβ42:Aβ40 ratio enhances the potential shift of soluble Aβ to a fibrillar oligomeric form. Alternatively, the increased abundance of fibrillar Aβ oligomers in the TBS-soluble fraction could result from reduced incorporation into highly condensed plaques in brains with Trem2-deficient microglia.
Apoe is one of the most highly induced genes in mouse microglia in response to neurodegenerative stimuli (Deczkowska et al., 2018). A recent report in the APPPS1 model found that plaque-associated ApoE was reduced in Trem2-deficient mice (Parhizkar et al., 2019). To learn whether Trem2 deletion altered ApoE localization in the PS2APP model and whether such alteration correlated with observed changes in plaque abundance, morphology, or composition, we costained for ApoE, plaques (methoxy-X04), and microglia (Iba1). In contrast to the findings of Parhizkar et al. (2019), quantification of plaque-associated ApoE in whole-brain sections found a significant increase in female PS2APP;Trem2ko mice at 6–7 and 12 months of age, but no significant changes in male cohorts at any age tested (Fig. 6A). The ApoE immunostaining pattern was especially prominent in, but not limited to, the hippocampal subiculum where plaque tends to first deposit (Fig. 6B). While the incongruous results between studies may be explained by differences between mouse models, other variations in sampling and technical procedures could also account for the differences. Parhizkar et al. (2019) examined 4 month males; we only examined females at that age, and our findings of elevated plaque-associated ApoE in Trem2-deficient females were only observed at older ages. The studies used different ApoE monoclonal antibodies that likely interact with distinct epitopes or configurations of ApoE (Kim et al., 2012) and whose binding may be differentially impacted by variations in staining procedure. We immunostained for ApoE first, followed by plaque labeling, which involves treatments with ethanol and NaOH. Parhizkar et al. (2019) did the plaque labeling first, followed by the ApoE immunostains, so the nature of the ApoE at the time of immunostaining was somewhat different between studies.
Parhizkar et al. (2019) also reported that colocalization of ApoE with Iba1+ microglia was reduced in Trem2-deficient mice. To assess this finding in the PS2APP model, we performed confocal imaging on the ApoE/Iba1/methoxy-X04 costains from the cortex of 12 month females and again observed a contrasting result. Instead of ApoE labeling being diminished, we observed microglia in PS2APP;Trem2ko mice to be markedly laden with ApoE (Fig. 6C,D), suggesting that Trem2-deficient microglia exhibit deficits in lipid clearance. At the mRNA level, we did not detect a difference in Apoe expression between PS2APP;Trem2wt and PS2APP;Trem2ko microglia (Fig. 6E; see also Fig. 2E). Although two groups have reported reductions in microglial Apoe expression in the APPPS1 model when Trem2 is deleted (Krasemann et al., 2017; Parhizkar et al., 2019), this does not appear to be a typical feature of Trem2 deficiency in β-amyloidosis models since previous analyses in the 5xFAD model observed Trem2-independent Apoe induction in either bulk microglia microarray (Fig. 6F) (Wang et al., 2015) or single-cell microglia RNA-Seq profiles (Keren-Shaul et al., 2017). In a very recent dataset of bulk microglia RNA-Seq profiles from 5xFAD mice (Griciuc et al., 2019), we did see a twofold decrease in the extent of Apoe induction in Trem2-deficient microglia, but the gene was still highly induced relative to the expression level in microglia from nontransgenic mice (Fig. 6F). Together, our evidence indicates that microglial Apoe expression is induced by Aβ-driven neuropathology in a largely Trem2-independent manner, and that Trem2-deficient microglia accumulate disproportionately large amounts of ApoE compared with the smaller ApoE puncta observed in normal PS2APP microglia (Fig. 6C).
Axonal dystrophy, dendritic spine loss, and CSF NfL detection are exacerbated in PS2APP;Trem2ko mice
The elevated Aβ42:Aβ40 ratio and fibrillar Aβ oligomers in PS2APP;Trem2ko brains would seem to be detrimental for neuronal health since Aβ42 oligomers are commonly understood to be the more toxic form of Aβ (Klein, 2002; Haass and Selkoe, 2007). Alternatively, the reduced plaque load in aged Trem2ko brains suggested a possible benefit of Trem2 deficiency. Therefore, we turned to measures of neuronal dystrophy to better understand the potential consequences of loss of Trem2 function.
First, we looked at neuritic dystrophy around plaque (D'Amore et al., 2003) by fluorescent costaining using methoxy-X04 to label plaque, anti-Iba1 to label microglia, and anti-Lamp1 to label dystrophic axons (Gowrishankar et al., 2015) (Fig. 7A). The methoxy-X04/Iba1 costain showed that microglial association with plaque was severely compromised in PS2APP;Trem2ko brains at all examined ages (Fig. 7A,B), corroborating our earlier Iba1 immunohistochemical stains that measured microglial clustering. The Lamp1 immunolabeling, which stains dilated dystrophic axons, revealed two important findings. First, on a per plaque basis, axonal dystrophy was exacerbated from 7 months onward in PS2APP;Trem2ko mice (Fig. 7C), similar to findings in other β-amyloid models (Wang et al., 2016; Yuan et al., 2016). This is consistent with the idea that the diffuse plaque structures and elevated Aβ42:Aβ40 ratio in Trem2ko brains are more damaging to surrounding axons than the more compacted, Aβ40-enriched plaques in PS2APP;Trem2wt brains. Second, the total Lamp1+ area was also increased from 7 months onward in PS2APP;Trem2ko brain sections (Fig. 7D), indicating that total axonal damage was exacerbated at later ages despite the reduced plaque burden. These data strongly suggest that the Trem2-dependent clustering of microglia and their functions around plaque serve to mitigate the neurotoxic effects of Aβ.
Another method we used to visualize neuronal pathology was an amino-cupric-silver stain or “disintegrative degeneration” stain that labels damaged or degenerating neurons (de Olmos et al., 1994; Switzer, 2000). Overall, the staining pattern appeared very similar to the Lamp1 stain, with “bouquets” of argyrophilic structures presumably surrounding plaques throughout the cortex and hippocampus. Staining was also observed in relevant white matter tracts, such as the corpus callosum, perforant path, and fornix, suggesting that degenerating axonal processes are not restricted to dystrophic neurites around plaques. Degenerating neurites detected by this stain were more abundant in PS2APP;Trem2ko brains of the 12 month female and 19–22 month male cohorts (Fig. 8A,B). Thus, again, axonal damage was exacerbated in Trem2-deficient mice from older ages despite the fact that plaque accumulation was reduced. The intermediate effect of Trem2 heterozygosity on microglial clustering around plaque (Figs. 1A, 7B) and Campbell-Switzer plaque staining (Fig. 3A), but not on plaque diffuseness (Fig. 4B), Aβ42:Aβ40 ratio (Fig. 5A), or neuritic dystrophy (Figs. 7C,D, 8A), suggested that the form rather than the amount of plaque correlates with neuronal damage, and that sufficient microgliosis occurs in PS2APP;Trem2het mice to enable plaque compaction and neuroprotection.
NfL measured in the CSF or plasma has recently emerged as a potential biomarker of neurodegeneration in human patients and in mouse disease models (Bacioglu et al., 2016; Khalil et al., 2018). To determine whether Trem2 deficiency altered NfL levels, we collected plasma and CSF from a separate, mixed sex cohort of 12 month PS2APP;Trem2wt and PS2APP;Trem2ko mice. PS2APP;Trem2ko mice had significantly greater NfL levels in the CSF compared with PS2APP;Trem2wt mice (Fig. 8C), consistent with the measures of increased axonal dystrophy that were observed by histopathology in Trem2-deficient females at this age. Perhaps surprisingly, we did not observe a difference in plasma NfL levels between genotypes (Fig. 8C), suggesting that CSF NfL measurements better represent ongoing neuronal damage or degeneration in the CNS than plasma measurements.
Finally, we looked at whether another feature of neuronal pathology observed in AD tissues and β-amyloid mouse models, reduced synaptic density, particularly near plaque (Spires and Hyman, 2004; Tsai et al., 2004; Spires et al., 2005), was altered in Trem2-deficient mice. In β-amyloid models, the reduction in synapse number requires the presence of microglia since depleting the microglial cell population largely prevents loss of synaptic density (Olmos-Alonso et al., 2016; Spangenberg et al., 2016). Therefore, we asked whether Trem2 deletion and the resulting lack of activated, plaque-associated microglia would prevent the dendritic spine loss from occurring or would worsen it. To answer this, we crossed the Thy1:GFP-M line, which labels a sparse population of excitatory neurons, into the PS2APP model and analyzed a cohort of 6 month female mice with different Trem2 genotypes. Spine density loss in the proximity of plaques was not rescued by Trem2 deletion, but was actually further exacerbated in PS2APP;Trem2ko mice compared with PS2APP;Trem2wt or PS2APP;Trem2het groups (Fig. 9). Thus, Trem2 is not required for microglia-mediated dendritic spine loss around plaque; and indeed, Trem2 seems to hedge against synapse loss. Overall, the exacerbated axonal dystrophy and dendritic spine loss observed around plaque in PS2APP;Trem2ko mice imply that Trem2-dependent microglial activity is fundamentally neuroprotective in β-amyloid-driven models of AD-like pathology.
Discussion
In this study, we examined the role of Trem2 in microglial activation, plaque accumulation, and neuronal dystrophy in the PS2APP model of β-amyloidosis. We observed both age- and sex-dependent effects of Trem2 deletion on plaque abundance assessed using the Campbell-Switzer silver stain, with slightly more plaque in PS2APP;Trem2ko females at 6–7 months of age but markedly less plaque in both female and male PS2APP;Trem2ko mice at later ages. Together with a recent similar report (Parhizkar et al., 2019), these results suggest that Trem2-dependent microglial activity may both restrain the formation/seeding of plaques at an early stage of pathology, conceivably through uptake and degradation of soluble Aβ species, but also enhance the sequestration of Aβ into existing plaque structures, particularly at later stages of pathology.
At all ages examined, and in both sexes, microglial clustering around plaque and other measurements of microglial activation were sharply reduced in PS2APP;Trem2ko mice. Transcriptional induction of the neurodegeneration-related modules, the proliferation module, and certain genes related to Wnt regulation were impaired in PS2APP;Trem2ko microglia. Presumably as a consequence of the impaired microglial response around plaque, plaques in PS2APP;Trem2ko brains displayed a more diffuse morphology than in PS2APP;Trem2wt brains. From 6–7 months onward, axonal injury was magnified in PS2APP;Trem2ko mice, even at later ages when the abundance of argyrophilic amyloid plaques was diminished. Our finding that plaque-proximal dendritic spine loss was exacerbated in PS2APP;Trem2ko mice further underscored that the Trem2-dependent microglial activation around plaque is protective for neurons and is distinct from microglia-mediated, complement-mediated activities that contribute to synapse loss in β-amyloid models (Fonseca et al., 2004; Hong et al., 2016; Shi et al., 2017; Wu et al., 2019). It will be important for future studies of microglial modulation in mouse models to analyze plaque-associated neuritic dystrophy and not assume that decreases (or increases) in amyloid plaque burden are evidence of beneficial (or detrimental) effects.
While the elevated Aβ42:Aβ40 ratio we observed in the soluble fraction of PS2APP;Trem2ko brains could reflect direct deficits in the Trem2-mediated binding and clearance of oligomeric Aβ42 by microglia (Yeh et al., 2016; Lessard et al., 2018; Zhao et al., 2018), it could also reflect enhanced clearance of Aβ40 through other means, such as vascular efflux. The resulting elevation in soluble Aβ42:Aβ40 ratio may give rise to the increased plaque formation we observed in 6–7 month females, similar to the increased plaque seeding activity observed by others at a similar age in another model when Trem2 was deleted (Parhizkar et al., 2019). Given that Aβ40 is more easily incorporated into existing dense plaque structures than Aβ42 (Condello et al., 2015), our observation that the Aβ42:Aβ40 ratio declines as the mice age, particularly in the insoluble fraction, is consistent with the accumulation of compact plaques in size and number over time in the brains of PS2APP mice. The reduction of this plaque accumulation activity in PS2APP;Trem2ko brains may underlie the elevated levels of fibrillar oligomeric Aβ we detected in the soluble fraction, which may in turn be a source of increased neuronal injury in these animals.
The weakness of correlation between plaque abundance and cognitive status in humans, along with several unsuccessful clinical trials aimed at preventing cognitive decline by reducing brain Aβ levels, has led some to question the relevance of Aβ in the etiology and progression of AD (Krstic and Knuesel, 2013; Morris et al., 2014; Itzhaki et al., 2016). Moreover, many have proposed that chronic microglial activation is a key damaging agent that contributes to the neurotoxic environment in AD (Heneka et al., 2015; Park et al., 2018). Our results suggest the opposite, since preventing the microglial response to Aβ pathology via Trem2 deletion reversed neither axonal dystrophy nor dendritic spine loss around plaque and indeed made both of these pathologies worse. Trem2 deletion also increased the amount of NfL detected in CSF. Evidence in humans and mice supports the mechanistic model that β-amyloid pathology accelerates the accumulation of tau pathology or enhances its spreading (Pooler et al., 2015; He et al., 2018; Jack et al., 2018; Jacobs et al., 2018), and a new report indicates that this effect is further magnified in Trem2-deficient mice (Leyns et al., 2019). Together, these findings suggest that the form of microglial activation brought on by Aβ-related pathologies (the DAM state, or at least the Trem2-dependent component of it) protects neurons by limiting Aβ-induced neuronal injury.
Supporting the notion that microglial activation is primarily beneficial in the context of AD pathology, elevated PET signal for ligands of the “neuroinflammation marker” TSPO have predicted better cognitive measures and slower AD progression in mice with β-amyloid pathology and human patients, respectively (Hamelin et al., 2016; Focke et al., 2019). The fact that plaques in Trem2-deficient mice are more injurious to adjacent neurites but show weaker labeling with molecular probes (Thioflavin S, methoxy-X04 and X-34) related to those used in the clinic (11C-PiB, 18F-Florbetapir) helps explain why cognitive decline correlates better with tau pathology and synapse loss than with brain amyloid detection. The form of β-amyloid in the brain is more critical than the amount, and Trem2-mediated microgliosis facilitates the consolidation of β-amyloid into a highly compacted, less damaging form. Therefore, therapeutics that enhance this microglial activity may prevent AD or delay its progression while simultaneously (and perhaps counterintuitively) leading to increases in PET signals for amyloid content and microgliosis. Clinical biomarkers of neuronal degeneration and spreading AD pathology, such as neurofilament-L and tau, should be more informative indicators of whether a microglia-directed therapy is achieving efficacy.
Footnotes
- Received August 1, 2019.
- Revision received December 8, 2019.
- Accepted December 23, 2019.
We thank Melissa Gonzales Edick and Joanna Yung for imaging support; members of the Genentech animal care staff, FACS laboratory, and RNA sequencing laboratory for other research support; and Chris Bohlen for helpful critique of the manuscript.
All authors are or were employees of Genentech, Inc., with interests in developing novel therapeutics for neurodegenerative diseases at the time of performing work for this manuscript.
- Correspondence should be addressed to David V. Hansen at hansen.david{at}gene.com or William J. Meilandt at meilandt.william{at}gene.com
- Copyright © 2020 Meilandt et al.
This is an open-access article distributed under the terms of the Creative Commons Attribution License Creative Commons Attribution 4.0 International, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.