Ulland 2017
Ulland 2017
Ulland 2017
Correspondence
[email protected]
In Brief
The Alzheimer’s disease risk factor
TREM2 regulates microglial function
through modulation of cellular
biosynthetic metabolism.
Highlights
d TREM2-deficient microglia undergo increased autophagy in
an AD mouse model
*Correspondence: [email protected]
http://dx.doi.org/10.1016/j.cell.2017.07.023
Cell 170, 649–663, August 10, 2017 ª 2017 Elsevier Inc. 649
Figure 1. Defect in TREM2 Enhances Autophagy In Vivo in the 5XFAD Mouse Model and in AD Patients
(A) TEM images of microglia (CD45+, CD11b+, F4/80+ cells) sorted from 8-month-old WT, Trem2–/–, 5XFAD, and Trem2–/– 5XFAD mice.
(B) Average number of multivesicular and multilamellar structures/cell (30 cells analyzed/genotype).
(C) Confocal images of plaque-bearing regions of the cortex (1.1 mm Bregma to 0.8 mm Bregma) of 8-month-old WT, Trem2–/–, 5XFAD, and Trem2–/– 5XFAD mice
show Iba-1+ microglia (red), methoxy X04+ plaques (blue), and LC3 (green). Z stacks composed of 30 images taken at 1.2-mm intervals were analyzed. Results
are reported as an average of 2 regions of interest (ROI) analyzed.
(D) Quantification of the % of microglia that are positive for LC3 puncta. 150–400 microglia/HPF were analyzed depending on the genotype of the animal.
(E) Confocal images of sections from post-mortem brains of R47H+ AD patients and case-matched controls (CV, common variant of TREM2) show Iba-1+
microglia (red), methoxy X04+ plaques (blue), and LC3 (green). 3 ROIs/donor were analyzed and between 400–700 microglia/individual were analyzed.
(F) Percentages of LC3+ microglia in post mortem specimens of AD patients with different genotypes.
Error bar represents mean ± SEM. ***p < 0.005, ****p < 0.001 by one-way ANOVA with Holm-Sidak’s multiple comparisons test. 15 cells from 2 separate mice
were visualized for TEM (A and B). Confocal images are representative of 3 female mice per group (C) or 7 R47H, 4 R62H, and 8 case-matched AD patients for
post-mortem specimens (E). Arrowheads indicate multilamellar and multivesicular structures (A) or LC3+ vesicles (C and E).
See also Figure S1 and Table S1.
decreased phosphorylation of 4EBP1, an mTORC1 effector, as with reduced mTOR activation and increased autophagy (Fig-
well as Akt (serine 473) and NDRG1, both mTORC2 effectors, ure 2A). Impaired mTOR signaling was associated with phos-
in the TREM2-deficient microglia (Figure 2A). Ulk1, a key inducer phorylation of AMP-activated protein kinase (AMPK), a sensor
of autophagy, which is inhibited by mTOR signaling through of low energy states (Figure 2A).
phosphorylation at serine 757 (Kim et al., 2011), was less phos- Providing further evidence for defective energetic and
phorylated in microglia from Trem2–/– 5XFAD mice, consistent anabolic metabolism associated with TREM2 deficiency,
Macrophages Lacking TREM2 Inadequately Signal TREM2 Deficiency Curtails Anabolic and Energetic
through mTOR and Undergo More Autophagy Metabolism in BMDMs
We then asked whether TREM2-deficiency could derail mTOR To directly demonstrate the impact of TREM2 deficiency on
signaling in bone marrow-derived macrophages (BMDMs) energetic and anabolic pathways in BMDMs, we performed
in vitro. To mimic the metabolic stress that occurs during dis- mass spectrometry to quantify cellular metabolites and RNA-
ease, we used growth-factor deprivation. BMDMs were cultured seq to quantify mRNA levels of metabolic enzymes. Analysis of
overnight in concentrations of CSF1-containing L cell-condi- metabolite data alone, or in combination with RNA-seq data by
tioned medium (LCCM) ranging from optimal to limiting (10% a systems-based algorithm (Sergushichev et al., 2016), revealed
to 0.5%). Trem2–/– BMDMs contained more autophagic vesicles widespread differences between WT and Trem2–/– BMDMs in
(Figures 3A and 3B) and had a higher LC3II/LC3I ratio than WT various metabolic pathways. Compared to WT BMDMs in
cells when CSF1 was limiting (Figure 3C). Addition of the lyso- optimal CSF1, Trem2–/– BMDMs cultured under the same condi-
somal inhibitor bafilomycin greatly increased LC3II in Trem2–/– tions exhibited (1) a marked decrease of key intermediates in the
BMDMs, confirming that the increase in autophagosomes was synthesis of nucleotides (e.g., phosphoribosyl pyrophosphate),
due to increased autophagic flux rather than reduced autopha- N-glycosylated proteins (e.g., UDP-glucose), and phospholipids
gosome degradation (Figures 3D and 3E). (e.g., CDP-ethanolamine); (2) a decrease in glycolytic metabo-
As observed in sorted microglia, autophagy in BMDMs was lites (e.g., glucose 6-phosphate and fructose bisphosphate)
linked to impaired mTOR signaling. Trem2–/– BMDMs had less and tricarboxylic acid (TCA) cycle intermediates (citrate and suc-
phosphorylated 4EBP1, Akt (serine 473), and NDRG1 (Figure 3F) cinate); and (3) an increase in catabolic products of amino acids
and more phosphorylated AMPK (Figure 3G) in both optimal (e.g., indolacetate) and phospholipids (e.g., glycerol 3-phos-
CSF1 and limiting CSF1 than did WT BMDMs. In limiting CSF1, phate) (Figures 4A and S4A). Moreover, a selective increase in
Trem2–/– BMDMs had decreased inhibitory phosphorylation of malate and fumarate suggested an enhanced malate-aspartate
Ulk1 at serine 757, while activating phosphorylation of Ulk1 at shuttle to sustain defective NADH oxidation and NAD regenera-
serine 317 increased (Figure 3H). Thus, lack of TREM2 sup- tion (Figure 4A). Unbiased network analysis combining metabolic
pressed mTOR activation and elicited compensatory AMPK and RNA-seq data highlighted defects in metabolites and
and Ulk1 activation and autophagy in BMDMs in response to enzymes involved in glycolysis, TCA cycle, and pentose phos-
metabolic stress, very similar to our observations of microglia phate pathway in Trem2–/– BMDMs (Figure 4B).
in 5XFAD mice. CSF1 reduction further deteriorated energy and anabolic
As TREM2 signaling adapters DAP12 and DAP10 have been metabolism in Trem2–/– BMDMs. Under these conditions, again
shown to activate PI3-K, which in turn can activate mTOR, we in comparison to WT BMDMs, Trem2–/– BMDMs underwent a
asked whether enhanced mTOR signaling in WT BMDM marked increase in symmetrical dimethyl arginine, indicative
compared to Trem2–/– BMDM was dependent on PI3-K. Inhibi- of protein catabolism, as well as an increase in ADP-ribose,
tion of PI3-K with Wortmannin or LY294002 caused a major indicative of NAD degradation (Figures 4C, S4B, and S4C).
reduction in phosphorylation of mTOR and its downstream Furthermore, stores of high-energy phosphates, such as
untreated Trem2–/– 5XFAD mice had a lower density than those did not moderate plaque accumulation or the engulfment of pla-
in 5XFAD mice as measured by methoxy-X04 staining intensity, que particulates by microglia, at least at this time point (Figures
the density of plaques in cyclocreatine-treated Trem2–/– 5XFAD S5G–S5I). As APP is known to accumulate as distinct rounded
mice resembled that of plaques in 5XFAD mice (Figure 7E), particles in dystrophic neurites, we used APP deposition around
although plaque shape complexity was not significantly altered plaques to assess neurite dystrophy (Masliah et al., 1996; Wang
(Figure S5F). Despite reducing plaque density, cyclocreatine et al., 2016; Yuan et al., 2016). Cyclocreatine treatment of
(F and G) Mitochondrial mass of WT and Trem2–/– BMDM assessed by Mito Tracker Green incorporation (F) and by the ratio of mitochondrial-to nuclear DNA (G).
Error bar represents mean ± SEM. *p < 0.05, **p < 0.01, or ****p < 0.001 by one-way ANOVA with Holm-Sidak’s multiple comparisons test (C) or Student’s t test
(F and G). Data are representative of at least 3 independent experiments.
See also Figure S4.
(C) Immunoblots performed on lysates of microglia sorted from the indicated genotype and treatment group of mice. Immunoblots for phosphorylated Akt (serine
473), NDRG1, total LC3, Akt, and actin.
(D) Quantification of the LC3II/LC3I ratio observed in immunoblots from 3 mice of each of the indicated genotypes and treatment groups.
(E) Average intensity of the plaques observed in the cortexes of mice from the indicated genotypes and treatment groups.
(F) Representative images depicting plaques (X04 in blue), nuclei (To-Pro3 in white), and N terminus APP (green) from the indicated mice and treatment groups.
Confocal images were taken as in Figure 1C.
(G) Quantification of the number of dystrophic neurites/plaque in the indicated mice and treatment group.
Error bar represents mean ± SEM. N.S. indicates not significant. *p < 0.05, and ****p < 0.001 by one-way ANOVA with Holm-Sidak’s multiple comparisons test
(A, C, D, and F) results pooled from 2 independent experiments representing a total of 5–8 male and female mice per treatment group.
See also Figure S5.
SUPPLEMENTAL INFORMATION Guerreiro, R., and Hardy, J. (2014). Genetics of Alzheimer’s disease. Neuro-
therapeutics 11, 732–737.
Supplemental Information includes five figures and two tables and can be Hara, T., Nakamura, K., Matsui, M., Yamamoto, A., Nakahara, Y., Suzuki-Mi-
found with this article online at http://dx.doi.org/10.1016/j.cell.2017.07.023. gishima, R., Yokoyama, M., Mishima, K., Saito, I., Okano, H., and Mizushima,
Further information and requests for resources and reagents should be directed to and will be fulfilled by the Lead Contact, Dr. Marco
Colonna ([email protected]).
Mice
Mice were of mixed sexes. Mice within experiments were age and sex matched. For studies using 5XFAD and Trem2–/– 5XFAD
animals all animals were 8 months of age at the time of use. For bone marrow and primary microglia mice were used from 6 weeks
of age until 12 weeks of age. Mice used in this study include WT C57BL/6J, 5XFAD, Trem2–/–, and Trem2–/– 5XFAD animals. All
animals were backcrossed until at least > 98% C57BL/6J confirmed by genotype wide microsatellite typing. Mice were housed under
specific pathogen free conditions. Mice from different genotypes were cohoused. Mice did not undergo any procedures prior to their
stated use. For cyclocreatine treatment mixed litters of sex matched mice were randomly assigned to experimental groups. All
studies performed on mice were done in accordance with the Institutional Animal Care and Use Committee at Washington University
in St. Louis approved all protocols used in this study.
METHOD DETAILS
Mice
The generation of Trem2–/– and Trem2–/– 5XFAD mice has been described previously (Oakley et al., 2006; Turnbull et al., 2006; Wang
et al., 2015). All mice were on a C57BL/6 background. Age and sex matched mice were used for all experiments; experimental
cohorts of mice were cohoused from birth to control for the microbiota. For in vivo cyclocreatine treatment 10-week old mice
were put on cyclocreatine-containing water, treatment was continued until mice reached 8 months of age (Santa Cruz
SC-217964 S). Desired intake of cyclocreatine was approximately 0.28 mg/g of body weight/day, which is approximately the
same as the standard creatine dose used in humans of 285 mg/kg of body weight/day (Kurosawa et al., 2012). Cyclocreatine was
administered in drinking water at a final concentration of 2.33 mg/ml. The Institutional Animal Care and Use Committee at Washington
University in St. Louis approved all protocols used in this study.
Microglia sorting
Microglia were isolated from the indicated animals as previously described (Wang et al., 2015). CD45+, CD11b+, F4/80+ (Biolegend
Cat. Number 103134, eBioscience Cat. Numbers 11-0112 and 17-4801) cells in the brain were fluorescence-activated cell-sorted
(FACS) directly into RLT-plus lysis buffer for microarray or 2% FBS in PBS for TEM or immunoblot lysates. For microarray RNA extrac-
tion was performed using a RNeasy micro kit (QIAGEN). Microarray hybridization (Affymetrix MoGene 1.0 ST array) and data process-
ing were performed at the Washington University Genome Center. For normalization, raw data was processed by Robust Multi-Array
(RMA) method and genes were pre-filtered for expression value R 120 expression units, a cut-off above which genes have a 95%
chance of expression demonstrated in Immgen dataset, which uses the same array platform (Wang et al., 2015). QIAGEN IPA anal-
ysis was preformed by comparing fold change and p values for all genes 5XFAD and Trem2-deficient 5XFAD microglia. Heatmaps
and hierarchical clustering were generated from preselected gene-lists using Morpheus. Microarray data has been deposited at
GEO: GSE65067.
Immunoblotting
BMDM or microglia were lysed in RIPA buffer (50 mM Tris, 150mM NaCl, 1% SDS, and 1% Triton X-100) containing PMSF, leupeptin,
activated sodium orthovanidate, apoprotinin, and phosphatase inhibitor cocktail 3 (Sigma Aldrich Cat. Number P0044). Lysates were
flash frozen on dry ice and stored at 80 C until use. Lysates were thawed and 4x LDS running buffer and 10% b-mercaptoethanol
were added. Lysates were heated to 95 C for 10 min and run on either a 15% polyacrylamide gel with a 4% stacking gel, a 12% bis-
tris gel (Nupage), or a 4%–12% bis-tris gel (Nupage). Proteins were transferred to nitrocellulose and blocked for 1 hr at RT in 5% milk
in Tris buffered saline with 0.05% Tween 20 (TBST). Membranes were incubated in primary antibody overnight at 4 C (Table S2).
Membranes were subsequently washed and incubated in Lienco anti-rabbit HRP for 1 hr at RT, washed, and developed using either
SuperSignal West Pico Chemiluminescent Substrate or a combination of SuperSignal West Pico Chemiluminescent Substrate and
SuperSignal West Femto Chemiluminescent Substrate.
RNA-seq analysis
Cells were cultured as described in the metabolite profiling by EIS-MS/MS section above. mRNA was extracted from cell lysates
using oligo-dT beads (Invitrogen). For cDNA synthesis, we used custom oligo-dT primer with a barcode and adaptor-linker sequence
(CCTACACGACGCTCTTCCGATCT-XXXXXXXX-T15). After first-strand synthesis, samples were pooled together based on Actb
qPCR values and RNA-DNA hybrids were degraded with consecutive acid-alkali treatment. Subsequently, a second sequencing
linker (AGATCGGAAGAGCACACGTCTG) was ligated with T4 ligase (NEB) followed by SPRI-beads (Agencourt AMPure XP,
BeckmanCoulter) clean-up. The mixture was enriched by PCR for 12 cycles and SPRI-beads (Agencourt AMPure XP,
BeckmanCoulter) purified to yield final strand-specific RNA-seq libraries. Libraries were sequenced using a HiSeq 2500 (Illumina)
using 50 bp x 25 bp pair-end sequencing. Second read (read-mate) was used for sample demultiplexing. Reads were aligned to
the GRCm38.p2 assembly of mouse genome using STAR aligner. Aligned reads were quantified using quant3p script (https://
github.com/ctlab/quant3p) to account for specifics of 30 sequencing. RefSeq genome annotation was used and DESeq2 was
used for differential gene expression analysis. RNaseq data has been deposited at GEO:GSE98563.
qRT-PCR
Total RNA was isolated with TRIzol Reagent (Invitrogen) and single-strand cDNA was synthesized with a High Capacity cDNA
Reverse Transcription Kit (Applied Biosystems). Genomic DNA was extracted using the QIAamp DNA micro kit (QIAGEN) to deter-
mine mtDNA/nDNA ratios. Real-time PCR was performed using SYBR Green real-time PCR master mix (Thermo-Fisher) and
LightCycler 96 detection system (Roche). mtDNA primers were to cytochrome c oxidase subunit 1 and nDNA primers were to
NADH:ubiquinone oxidoreductase core subunit V1.
Metabolism assays
For real-time analysis of extracellular acidification rates (ECAR) macrophages were analyzed using an XF96 Extracellular Flux
Analyzer (Agilent Technologies). Cells were incubated overnight in complete RPMI in the indicated concentration of LCCM with or
without cyclocreatine (10 mM). Measurements were taken under basal conditions and following the sequential addition of 1 mM
oligomycin and 1.5 mM fluoro-carbonyl cyanide phenylhydrazone (FCCP) (purchased from Sigma-Aldrich).
Data in figures are presented as mean ± SEM. Unless otherwise stated statistical analysis was performed using Prism (GraphPad).
Quantification of confocal microscopy, immunoblots, and electron microscopy images were performed using Imaris, ImageJ,
MATLAB, and FIJI. Differential metabolite expression was analyzed using limma. Pathway analysis of microarray data was performed
using IPA software. RNaseq analysis was performed by using Second read (read-mate) for sample demultiplexing. Reads were
aligned using STAR aligner and quantified using quant3p scirpt. RefSeq genome annotation was used and DESeq2 was used for
differential gene expression analysis. Combined RNaseq and metabolite network analysis was performed utilizing Shiny GAM. Sta-
tistical analysis to compare the mean values for multiple groups was performed using Prism by one-way ANOVA with Holm-Sidak’s
multiple comparisons test. Comparison of two groups was performed in Prism using a two-tailed unpaired t test (Mann Whitney).
Values were accepted as significant if p % 0.05. Intragroup variation compared between groups was similar in all experiments.