
Trem2 is essential for microglia-mediated synaptic refinement, but whether Trem2 contributes to shaping neuronal development remains unclear. Tagliatti, Desiato, et al. demonstrate that microglial Trem2 regulates the metabolic fitness of hippocampal neurons in a region-specific manner. Lack of Trem2-mediated communication between microglia and neurons during early development links metabolic derangements with altered synapses and circuit maturation.
Erica Tagliatti,1,2,8 Genni Desiato,1,8 Sara Mancinelli,3 Matteo Bizzotto,1,3 Maria C. Gagliani,4 Elisa Faggiani,1
Rebeca Herna´ ndez-Soto,1 Andrea Cugurra,1 Paola Poliseno,1 Matteo Miotto,1 Rafael J. Arguello, € 5 Fabia Filipello,1,6
Katia Cortese,4 Raffaella Morini,1 Simona Lodato,1,3 and Michela Matteoli1,7,9, *
1 IRCCS Humanitas Research Hospital, via Manzoni 56, Rozzano, 20089 Milan, Italy
2 Department of Clinical and Experimental Epilepsy, UCL Queen Square Institute of Neurology, University College London, London, UK
3 Humanitas University, Department of Biomedical Sciences, Via Levi Montalicini 4, Pieve Emanuele 20072 Milan, Italy
4 Cellular Electron Microscopy Laboratory, Department of Experimental Medicine (DIMES), Human Anatomy, Universita` di Genova, Via Antonio de Toni 14, 16132 Genova, Italy
5 Aix Marseille Univ, CNRS, INSERM, CIML, Centre d’Immunologie de Marseille-Luminy, Marseille, France
6 Department of Psychiatry, Washington University School of Medicine, St. Louis, MO, USA
7 Institute of Neuroscience - National Research Council, 20139 Milan, Italy
8 These authors contributed equally
9 Lead contact
* Correspondence: michela.matteoli@hunimed.eu
https://doi.org/10.1016/j.immuni.2023.12.002
Triggering receptor expressed on myeloid cells 2 (Trem2) is a myeloid cell-specific gene expressed in brain microglia, with variants that are associated with neurodegenerative diseases, including Alzheimer’s disease. Trem2 is essential for microglia-mediated synaptic refinement, but whether Trem2 contributes to shaping neuronal development remains unclear. Here, we demonstrate that Trem2 plays a key role in controlling the bioenergetic profile of pyramidal neurons during development. In the absence of Trem2, developing neurons in the hippocampal cornus ammonis (CA)1 but not in CA3 subfield displayed compromised energetic metabolism, accompanied by reduced mitochondrial mass and abnormal organelle ultrastructure. This was paralleled by the transcriptional rearrangement of hippocampal pyramidal neurons at birth, with a pervasive alteration of metabolic, oxidative phosphorylation, and mitochondrial gene signatures, accompanied by a delay in the maturation of CA1 neurons. Our results unveil a role of Trem2 in controlling neuronal development by regulating the metabolic fitness of neurons in a region-specific manner.
Over the last few years, the outdated concept that the brain is an immunologically privileged organ has been replaced by the observation that a continuous crosstalk occurs between the nervous system and the immune system, especially during development and aging. Microglia, the main brain residential immune cells, represent the starring actor in these processes. Besides representing a line of defense against pathogenic insults, microglia are also emerging as centrally involved in physiological functions essential for correct CNS development and plasticity, regulating neuronal apoptosis, neurogenesis, myelin formation, and removal of supernumerary synapses during development. 1–9
Triggering receptor expressed on myeloid cells 2 (Trem2) is an immunoglobulin superfamily transmembrane receptor, expressed in the brain exclusively by microglia,10–15 which controls the functional microglia profile by regulating their energetic metabolism.16–20 Upon engagement with membrane-bound or soluble ligands, including lipids, lipoproteins, DNA, and bacterial products, Trem2 activates an intracellular signal transduction signaling via the adaptor proteins DAP10 and DAP12, and it enhances microglial phagocytosis of apoptotic neurons, cellular debris, bacterial products, and protein aggregates, including the neurotoxic -amyloid peptides. Further, Trem2 promotes bmyeloid cell survival and proliferation21–25 a nd modulates inflammatory signaling, controlling the switch from homeostatic to disease-associated microglia (DAM) state.26 Cleavage of membrane-bound Trem2 by the -secretases ADAM10 and aADAM17,27,28 or translation of an alternative spliced Trem2 transcript,29 leads to the release into the extracellular environment of soluble Trem2 (sTrem2), which maintains its biological activity.30–32 Genome-wide association studies demonstrated that Trem2 missense homozygous and heterozygous variants associate to neurodegenerative diseases, in particular Alzheimer disease (AD).33–37 Specifically, Trem2 variants deprive microglia of the dedicated sensory mechanisms required to detect damage within the CNS, also preventing them from acquiring the transcriptional and functional signatures needed to fight the pathological conditions.38,39 Accordingly, in AD models or in demyelination paradigms, lack of Trem2 results in increased neuritic dystrophy and axonal damage around amyloid plaques.18,40–42 The possibility that microglia, through Trem2, may shape the neuronal molecular profile and that the protein loss may force neurons to acquire altered transcriptional and molecular profiles has never been addressed.

Figure 1. Energetic metabolism is compromised in early post-natal neurons from Trem2–/– hippocampi
(A) Cartoon showing the experimental procedure.
(B) Average OCR kinetics showing the response of hippocampal slice punches from Trem2+/+ and Trem2/ newborn P1 pups to different mitochondrial inhibitors. Mean ± SEM.
© From left to right: basal, maximal, ATP-coupled respiration, and non-mitochondrial oxygen consumption. Box-and-whiskers plots of n = 6 vs. 12 punches from
N = 3 Trem2+/+ and Trem2/ P1 newborn pups. One sample t test, **p < 0.01 and *p < 0.05.
(D) Description of the experimental design.
(E) Flow cytometry gating strategy on Trem2+/+ and Trem2/ purified CD11b+CD45int microglia and CD11bCD45NeuN+ neurons.
(F) Puromycin incorporation following the treatment with metabolic inhibitors on microglia and neurons.
(G) Puromycin mean fluorescence intensity (MFI) in Trem2+/+ microglia and neurons.
(H and I) Puromycin MFI in microglia (H) and neurons (I) from Trem2+/+ and Trem2/ P1 hippocampi. Box-and-whiskers plots of N = 16 Trem2+/+, N = 12 Trem2/
P1 mice, of N = 5 independent experiments. Two-way ANOVA followed by Tukey’s test, **p < 0.01.
(J) Flow cytometry gating strategy on Trem2+/+ and Trem2/ purified CD11b+CD45int microglia.
(K and L) FACS Mitotracker and CellROS green dye MFI from Trem2+/+ and Trem2/ P1 microglia. Box-and-whiskers plots of N = 12 Trem2+/+, N = 10 Trem2/
from P1 mice, N = 2 independent experiments. Unpaired Student’s t test, **p < 0.01.
(M) Cartoon showing the experimental procedure.
(N) Left: average OCR kinetics showing the response of primary microglia from Trem2+/+ or Trem2/ newborn pups to different mitochondrial inhibitors.
(O) From left to right: basal, maximal, ATP-coupled respiration, and non-mitochondrial oxygen consumption. Mean ± SEM of at least n = 15 replicates from N = 4
independent Trem2+/+ and Trem2/ cultures. One sample t test.See also Figure S1.
Trem2 plays key roles also during neurodevelopment, when it controls the microglia-mediated process of supernumerary synapse elimination, influencing neuronal wiring and brain connectivity.43–45 At these stages, the protein is highly expressed and only modestly increases as a function of age.46,47 We focused on early developmental stages, when critical events of neuronal maturation occur, to investigate whether Trem2 contributes to shaping the neuronal profiles. We demonstrate that lack of microglial Trem2 functionally impacts neuronal metabolic signatures. We also demonstrate that this imbalanced pattern is associated with delayed developmental trajectory of Cornus Ammonis (CA)1 hippocampal pyramidal neurons.
Since Trem2 regulates circulating macrophage and microglial bioenergetics,16 we investigated the impact of Trem2 deficiency on the overall mitochondrial metabolism in P1 acute hippocampal sections from Trem2+/+ and mice by measuring oxygen consumption rates (OCRs) through the Seahorse extracellular flux analyzer platform (Figure 1A). OCR measurements showed that post-natal hippocampi display overall reduced mitochondrial respiration kinetics, compared with (Figure 1B), with significant reduction in basal, maximal, and ATP-dependent respiration (Figure 1C), suggesting an impaired ability of hippocampi to meet the overall energetic demand. Hippocampal slices from P1 mice also showed a lower non-mitochondrial-related respiration (Figure 1C), indicating a lower capacity of oxidizing by other forms of enzymatic activities.48
We assessed whether the observed reduction was associated to microglial alterations in the absence of Trem2. Differently from P18 where the density of microglia in CA1 region is lower, compared with Trem2+/+, 43 no significant differences occurred in the number of hippocampal Iba1-positive cells in either CA1 or CA3 regions of P mice (Figures S1A and S1C). Conversely, morphological analysis of P1 microglia revealed a higher number of branches and junctions per cell (Figures S1B, S1D, and S1E) and a higher expression of the homeostatic receptor P2yr12 (Figure S1F). These changes were statistically significant only in the CA1 region (Figures 1D– 1F, left).
To directly discriminate the contribution of different cell types to the altered metabolic phenotype, we exploited a flow cytometry-based method to functionally profile energy metabolism at single-cell resolution, namely SCENITH (single-cell energetic metabolism by profiling translation inHibition)49 (Figure 1D). This assay relies on the concept that most of the energy deriving from glucose, amino acids, and/or lipid catabolism in living cells is consumed by protein synthesis machinery.50 Thus, protein synthesis represents a proxy measure of global metabolic activity. We used the incorporation of puromycin combined with an anti-puro monoclonal antibody49 as a reliable readout for measuring protein synthesis51–53 and profiling the metabolic phenotype at single-cell resolution.
Both microglia and NeuN+ neurons purified from P1 hippocampi were processed (Figure 1E). We observed that P1 microglia equally rely on glycolysis and mitochondria for ATP-supplying protein synthesis, since both deoxy-glucose (DG) and oligomycin (O) treatments, which block glycolysis and mitochondrial ATP synthase, respectively, decreased puromycin incorporation (Figures 1F and 1G). Conversely, P1 neurons largely rely on mitochondria as an ATP source, as they did not display any changes in puromycin incorporation after DG (Figure 1G).
By comparing and microglia, we observed a mild reduction in their basal metabolism (Figure 1H). The metabolic defect in microglia was corroborated by mitotracker staining in microglia fluorescence-activated cell sorting (FACS)-isolated from P1 hippocampi (Figure 1J), which demonstrated a reduction of mitochondrial mass (Figure 1K), without changes in reactive oxygen species (ROS) production (Figure 1L). Furthermore, primary cultured microglia isolated from P1 brains and analyzed by Seahorse assay (Figure 1M) showed a negligible impairment in metabolic fitness relative to (Figure 1N), without gross changes in the main OCR parameters (Figure 1O).
Unexpectedly, we observed that post-natal hippocampal neurons, whose density was unaltered in either CA1 or CA3 regions (Figures S1G and S1H), displayed a significant impairment in mitochondrial metabolism (Figure 1I), suggesting that the lack of Trem2 impacts neuronal metabolism at early developmental stages.
To investigate how the lack of Trem2 in developing microglia impacts neuronal behavior, we transcriptionally profiled hippocampi from and littermates at P1, through single-cell RNA sequencing (scRNA-seq)54 (Figures 2A and 2B). Proportionally distributed between the 2 genotypes, 26,326 cells were profiled (Figure 2C), and unsupervised clustering analysis identified 15 clusters subsequently assigned to 13 distinct cell types, according to their canonical signatures (Figures 2C and S2B).
As expected, we identified neuronal progenitors (apical progenitor [AP] and intermediate progenitor [IPC]) and glial precursors (astrocyte and oligodendrocyte precursors). Distinct classes of postmitotic excitatory neurons of the CA (CA1- and CA3-Pyr) and the dentate gyrus (granules and mossy cells) and inhibitory GABAergic interneurons derived from medial and caudal ganglionic eminences (MGE and CGE) could be distinguished. An abundant cluster of immature CA-Pyr neurons was also identified cells), highlighting the dynamics of neuronal differentiation present at this stage.55 We also found microglial cells characterized by the expression of Trem2, Tmem119, Iba1, Cxc3r1, Csf1r, and Tyrobp (Figures S2A and S2B). As previously reported, Trem2 expression was found from early post-natal stages46,47 to adulthood (P90, Figure S2C).
By performing differential cell abundance testing for cell types,56 we did not observe significant changes in the composition of P1 hippocampi isolated from and littermates (Figures 2C, S2D, and S2E). However, differential gene expression analysis (Figure 2D) showed a significant difference in the transcriptional footprints of the excitatory pyramidal and granule cells (DG-Granule: 569 deregulated genes, 279 upregulated and 290 downregulated; CA1-Pyr: 487 deregulated genes), suggesting a prominent alteration of neuronal molecular signatures of the immature subtypes (DG-Granule and ${ \mathsf { C A - l m m } } $ and maturing CA neurons. In contrast, inhibitory neurons and non-neuronal cell types (including microglia), albeit less abundant, displayed limited transcriptional changes upon the missing Trem2 signaling (Figure 2D).
In agreement with our previously identified region-specific microglia defect in adult Trem2 mice,43 we observed that CAImm, CA1-Pyr, and CA3-Pyr showed a distinct pattern of transcriptional deregulation with limited gene overlap (Figure 2E). These results suggest that subtype-specific responses are triggered in the absence of Trem2 signaling. Indeed, the expression of most of the genes deregulated in CA neuronal classes isolated from Trem2 pups was selectively altered in each class (Figure 2E), and only a small fraction was shared among the other classes (Table S1). Pathways enrichment analysis (differentially expressed genes [DEGs]) in all CA neuronal clusters revealed pervasive deregulation of key biological processes in all neuronal types (Figure 2F; Table S2) related to translational processes (eukaryotic translation initiation factor 2 [EIF2] signaling and mammalian target of rapamycin [mTOR] pathways), supporting our previous results from SCENITH analysis (Figure 1I). Further underpinning the region-specific response, genes belonging to mTOR and mitochondrial dysfunction pathways displayed a selective enrichment only in CA1-Pyr and CA-Imm neurons but not in CA3-Pyr (Figures S2G and S2H). Also, metabolic alterations were identified as a prominent feature of their neuronal footprints: major dysregulations are on the mitochondrial machinery, without transcriptional involvement of the glycolysis (Figure 2I).
The oxidative phosphorylation (OXPHOS) pathway was severely deregulated in neurons. Indeed, we observed important changes in the behavior of genes associated with the respiratory chain, encoding all complexes from I to V, except complex II (Figure 2J). We found genes belonging to complex I (e.g., Ndufa5 and Ndufb2), to complex III (e.g., Uqcrq), and to complex IV (e.g., Cox6a1 and ), which were selectively downregulated in CA1 but not CA3 pyramidal neurons (Figure 2H). We also observed CA1 (both CA1-Pyr and CA-Imm) subtype-specific deregulation of genes belonging to complex V (e.g., Atp5e), coding the catalytic unit of the ATP synthase complex responsible for ATP production in the mitochondria (Figure 2H).
Overall, this evidence indicates that neurons, when microglial Trem2 is lacking, undergo a deep transcriptional rearrangement that selectively impacts the neuronal bioenergetics. This is accompanied by a predicted inhibition of ATP synthesis in neuronal cells committed to specific regional fates (i.e., CA1 neurons).
Since mitochondrial metabolism is a key driver of cell fate transitions in the brain,57–59 we tested whether the observed downregulation of genes involved in the OXPHOS pathways in the absence of microglial Trem2 could be involved in possible alteration of differentiation/developmental dynamics of hippocampal neuronal populations. To understand how the lack of Trem2 in microglia impacts neuronal commitment, we examined the developmental route of the neuronal lineages by inferring trajectory topology and pseudotime (Figure 3A). We first captured all the relevant cells in a standard model of hippocampus development and identified the portion of reduced dimensional space corresponding to the excitatory lineages, which included the cell types that showed the most pronounced transcriptional alteration upon lack of Trem2 in microglial cells (Figure 2D). We excluded inhibitory interneurons, as their trajectories could not be captured because the precursors are not localized within the hippocampus60 and mossy fiber cells, whose generation occurs earlier in hippocampal neurogenesis ( E10–E12.5).61

Figure 2. Single-cell transcriptomic analysis of Trem2+/+ and Trem2–/– P1 hippocampal tissue reveals metabolic reprogramming of devel oping neurons
(A) Schematic of hippocampal dissection and dissociation from P1 Trem2+/+ and Trem2/ mice.
(B) Uniform manifold approximation and projection (UMAP) plot colored by cell-type annotation.
© Beeswarm plot of the cell abundancy changes across genotypes.
(D) Bar plot showing the number of differentially expressed genes (DEGs) for each cell type (model-based analysis of single-cell transcriptomics [MAST]-re, p adjusted ).
(E) Venn diagram showing the DEGs overlap between the selected cell types.
(F) Dot plot showing the pathway enrichment analysis for individual cell types.
(G–I) Average expression heatmap of selected cell types showing deregulated genes related to selected pathways.
(J) Schematic showing predicted inhibition complexes and downregulated genes belonging to OXPHOS pathway.
See also Figure S2.

Figure 3. Mitochondrial defects in Trem2–/– newborn mice specifically disrupt the maturation of CA1 neurons
(A) UMAP plot showing the lineages trajectory. Density plot by cell types colored by pseudotime score value.
(B–D) Density and empirical cumulative distribution function plots show the distribution of the cells along the pseudotime between genotypes by individual cell type per lineage (KS test).
(E) Violin plots of selected DEGs associated with mitochondrial function in cells belonging to pseudotime interva split by genotype. MAST-re statistic test, *p adj , adj , and ***p adj .
(F) Representative confocal images of P1 CA1 (top) and CA3 (bottom) subfields of Trem2+/+ and Trem2/ hippocampi labeled with antibodies against doublecortin (Dcx) (green) and Tom20 (red). Scale bars, .
m(G and H) Tom20 mean fluorescence intensity in CA1 and CA3 pyramidal layers, respectively. Data are presented as mean SEM of Trem2+/+ and Trem2/ P1–2 pups. Data are normalized over Trem2+/+ counterpart. Unpaired Student’s t test, .
See also Figure S3.

Figure 4. Mitochondrial maintenance in neurons requires Trem2 signaling (A) Confocal images of P18 CA1 (left) and CA3 (right) subfield of Trem2+/+ and Trem2 $ / -$ hippocampi labeled with antibodies (Abs) against Map2 (magenta) and Tom20 (red). Scale bars, .
Upon trajectory inference, we detected three different excitatory lineages represented by multiple cell types (Figure 3A). From the common starting point (‘‘root’’) of the APs, the trajectory proceeds through IPCs to split into three distinct sublineages representing the developmental fate acquisition of the three different neuronal subtypes, which are defined by alternative ending branches (Figure 3A). No changes in the lineage topography were observed between and control littermates; however, the pseudotime analysis revealed differences in the cell distribution along the trajectory branches between conditions (Figures 3B–3D). The density distribution analysis (Figures S3A and S3B) revealed a significant difference in cell frequency along the pseudotime, specifically for CA1 lineage (Kolmogorov-Smirnov [KS] test: and not for the other two cell lineages (KS: (Figures S3A and S3B).
We also analyzed each lineage by parsing individual cell types for the two genotypes, focusing our attention on CA neurons. The analysis revealed that the common progenitor cells of the CA1 and CA3 lineages (AP and IPC) display similar distribution along the trajectory (Figures 3B and S3C). The analysis of the density distribution of the CA-ImmPyr exhibited a significant increase in both lineages (CA3 lineage, KS: , Figure 3C, top panel; CA1 lineage, KS: , Figure 3D, top panel) in the intermediate pseudotime interval (pseudotime score 17–22), which also included maturing CA1 neurons from , compared with . Of note, the accumulation of cells in the intermediate pseudotime interval coupled with a significant reduction of the terminally differentiated states (pseudotime score 22–25) only for CA1-Pyr neuronal cell types, (KS: , Figure 3C, bottom panel) but not for CA3-Pyr (KS: , Figure 3D, bottom panel). These data provide support for impairment in the dynamics of excitatory neuronal development in the absence of Trem2 and highlight a delayed maturation, specifically affecting the excitatory CA1 neurons.
Next, we investigated changes in gene expression among cells exhibiting different densities across the specific pseudotime value interval (17–22, gray) and identified, among the significantly deregulated genes between conditions, key mitochondrial genes: structural regulators of the outer mitochondrial membrane, crucial for the assembly of mitochondrial complexes (e.g., Tomm7), showed significant downregulation in cells, together with critical functional regulators of mitochondrial respiration machinery, such as Atp5e, Cox6a1, Cox7c, Cox7b, Cox7a2, Cox8a, and (Figure 3E; Table S1) (p adjusted ).
Finally, to assess whether the observed transcriptomic dysregulation of mitochondrial pathways in hippocampal neurons was reflected by an altered mitochondrial organization, we immunostained P1 sections from and brains against the translocase of the outer mitochondrial membrane complex subunit 20 (Tom20, Figure 3F), a reliable marker of mitochondrial mass.62 Confocal analysis showed a reduction in Tom20 mean intensity, significant only in CA1 region of hippocampi, compared with (Figures 3G and 3H).
Overall, these data indicate that lack of microglial Trem2 induces a shift in the developmental trajectory of pyramidal neurons, with a specific delay effect in the CA1 lineage, which is associated with alterations of gene programs controlling mitochondrial assembly.
Mitochondria are dynamic organelles that continuously remodel their morphology by fission and fusion to integrate bioenergetics and signaling in cells.63 The reduction in Tom20 mean intensity was detected in CA1 but not CA3 region of hippocampi, even at later developmental stages (P18, Figures 4A– 4C). The quantitative evaluation of the mitochondrial structure in CA1 and CA3 regions of or hippocampi by transmission electron microscopy (TEM) (Figure 4D) revealed a significant reduction in average mitochondria length in the soma of CA1 but not CA3 pyramidal neurons (Figures 4E and 4F).
In line with mitochondria ultrastructure reflecting their functionality, 64 OCR measurements in tissue punches obtained from CA1 or CA3 areas from and mice (Figure 4G) indicated that CA1 (Figures 4H and 4I) but not CA3 (Figures 4J and 4K) subfields display a selective impairment in all OCR parameters.
These findings indicate that the genetic lack of Trem2 in microglia induces a CA1-targeted selective impairment of OXPHOS chain accompanied by mitochondrial structural alterations, which is maintained at later developmental stages.
We next investigated the basis for the more prominent effect of Trem2 deficiency in CA1 pyramidal neurons. As previously shown,43 microglia at P18 were more enriched in the CA1 than in the CA3 region (Figures 5A and 5B, left panel), with no differences in cell volume revealed by Iba1 staining (Figures 5A and 5B, middle panel). Quantitative confocal analysis revealed that the expression of Trem2 relative to Iba1 is significantly higher in CA1 relative to CA3 (Figures 5A and 5B, right panel). Consistently, real-time quantitative PCR performed on 1-mm punches obtained from different regions of P18 hippocampal slices revealed a higher amount of Trem2 mRNA in the CA1 region (Figure 5C). Also, while microglial phagocytosis of supernumerary synapses is impaired in the CA1 region of hippocampi,43 we did not detect, in the CA3 region, any relevant difference in synapse engulfment (Figures 5D and 5E), frequency of miniature excitatory postsynaptic currents (mEPSCs, Figures 5F and 5G), and density of PSD95-labeled excitatory synapses (Figures 5H and 5I). Thus, differently from CA1, CA3 is not heavily impacted by the genetic lack of Trem2, suggesting that different mechanisms, independent from Trem2, may take place in this hippocampal region. Analysis of the somatosensory cortex (SSCx) revealed a significantly lower density of microglia, compared with hippocampal CA1 (Figure S5B), and a slightly lower expression of Trem2 quantified for single cell (Figure S5C, left). Quantitation of Trem2 volume in the entire fields confirmed that the protein expression was significantly higher in the CA1 region, with no difference between CA3 and SSCx (Figure S5C, right). Furthermore, similarly to the CA3 region, no differences in Tom20 expression were detected in the SSCx between and mice (Figure S5E). Thus, the CA1 region of the hippocampus seems to be more prominently dependent on Trem2-bearing microglia for its proper function.

Figure 5. Trem2 is highly expressed in the CA1 region, and its partial reduction is sufficient to deregulate neuronal metabolic fitness
(A) Confocal images and 3D reconstruction of microglial cells from P18 CA1 (up) and CA3 (down) subfields of hippocampi labeled with Abs against Iba1 (red) and Trem2 (green). Scale bars, .
(B) Microglial number per hippocampi (left), Iba1 volume (middle), and normalized Trem2 volume on Iba1 labeling (right) per cell in CA1 and CA3 subfields. Box-and-whiskers plots of N = 6 Trem2+/+ and Trem2−/− P18 mice. Unpaired Student’s t test, p < 0.05 and p < 0.01.
© Left: cartoon of the experimental procedure. Right: qPCR analysis of Trem2 mRNA in CA1 and CA3 punches from n = 8 P18 Trem2+/+ mice. Unpaired Student’s t test, p < 0.01.
(D) Representative confocal image and 3D reconstruction of microglial cells from P18 CA3 subfields of Trem2+/+ and Trem2−/− hippocampi labeled with Abs against Iba1 (green), CD68 (white), and Iba1 (red). Scale bars, 50 μm.
(E) Left: PSD95 (left), normalized CD68 volume on Iba1 staining (middle), and PSD95 volume engulfed in microglia (right). Mean ± SEM of N = 4 Trem2+/+ and N = 5 Trem2−/− P18 mice. Unpaired Student’s t test.
(F) Cartoon showing the experimental procedure.
(G) Left: representative traces of mEPSC recordings in the CA3 region of hippocampal slices from P18 Trem2+/+ and Trem2−/− mice. Scale bars, 20 pA, 200 ms. Down: quantification of mEPSC frequency (left) and amplitude (right) of CA3 pyramidal neurons. Box-and-whisker plots of n = 12 (Trem2+/+) and 13 (Trem2−/−) neurons from n = 3 Trem2+/+ and Trem2−/− mice. Unpaired t test.
(H) Confocal images of CA1 (left) and CA3 (right) subfields of P18 Trem2+/+ and Trem2−/− hippocampi labeled with Ab against PSD95. Scale bars, 10 μm.
(I) PSD95 volume fraction in CA1 (left) and CA3 (right) region. Mean ± SEM of n = 4 Trem2+/+ and Trem2−/− P18 mice. Unpaired t test.
(J) Confocal images of P1 CA1 (left) and CA3 (right) subfield of Trem2+/+, Trem2+/-, and Trem2−/− hippocampi labeled with Abs against Dcx (green) and Tom20 (red). Scale bars, 75 μm.
(K) Tom20 mean intensity in CA1 (left) and CA3 (right) pyramidal layers, respectively, from Trem2+/+, Trem2+/-, and Trem2−/− P1 mice. Mean ± SEM of N = 3 Trem2+/+ or Trem2+/- and N = 4 Trem2−/− P1 mice. One-way ANOVA test followed by Bonferroni, p < 0.05.
(L) Cartoon of the experimental procedure.
(M) Average OCR kinetics showing the response of hippocampal slice punches from Trem2+/+ and Trem2−/− newborn P1 pups to different mitochondrial inhibitors. Mean ± SEM.
(N) From left to right: basal, maximal, ATP-coupled respiration, and non-mitochondrial oxygen consumption. Box-and-whiskers plots of n = 8 vs. n = 10 punches from n = 3 Trem2+/+ and N = 4 Trem2−/− P1 newborn pups. One sample t test, p < 0.01 and p < 0.05.
See also Figure S5.
To assess whether the CA1 region strictly requires microglia carrying high Trem2 expression, we exploited (hemizygous) mice to test whether even partial reductions of Trem2 expression could be sufficient for inducing neuronal metabolic derangements. We found that both P1 (Figures 5J and 5K) and P18 (Figure S5F) hippocampi display reduced Tom20 mean intensity, specifically in the CA1 regions. Accordingly, OCR measurements revealed that P1 Trem2+/ hippocampi, similarly to (Figures 1A–1C), display overall reduced mitochondrial respiration kinetics (Figures 5M and 5N), suggesting their impaired ability to meet the overall energetic demand.
These data indicate that neurons in the CA1 hippocampal region require high microglial Trem2 for proper differentiation and functioning and that even a partial reduction of Trem2 is sufficient to derange the neuronal metabolic fitness.
We next investigated whether the metabolic alterations in P18 hippocampal neurons impact synaptic function and neurotransmission. P18 mice lacking Trem2 display a significantly higher number of excitatory synapses in the CA1 region, as contributed by the lack of supernumerary synapse elimination by microglia during the period of synaptic pruning.43 Here, we analyzed the ultrastructural features of these synapses. By TEM, we observed that synapses in the CA1 pyramidal layer contain mitochondria displaying a significant reduction in length (Figures 6A and S6A). This was accompanied by reduced synaptic area and increased synaptic vesicle density (Figures 6B and S6B). No differences were, instead, observed in active zone length and docked synaptic vesicle density (Figures 6B and S6B). Differently from CA1 synapses, CA3 synapses did not show alterations in mitochondria length, synaptic area, or synaptic vesicle density (Figures 6C, 6D, S6C, and S6D), also not in active zone length and docked vesicle density (Figures 6D and S6D).
To probe neuronal network functionality featuring both CA1 and CA3 regions in P18 mice lacking microglial Trem2, we employed multi-electrode arrays (MEAs) electrophysiology, which gives a high-throughput readout of field potential electrical activity of neuronal populations at once (Figures 6E and S6E). We measured mean firing rate (MFR) as a proxy of spontaneous neuronal activity in whole hippocampal slices (Figure 6E) and found that MFR in CA1 (Figure 6F) but not in CA3 (Figure 6G) was significantly higher in P18 slices. Such an effect was mainly due to an average increase in neuronal spontaneous activity and not in the number of firing cells65 (Figure S6F). Also, P18 Trem2 hippocampi showed a higher synchronicity and connectivity, compared with Trem2+/+ ones (Figures 6H and 6I). The increase in spontaneous neuronal activity is in line with the higher density of excitatory synaptic markers and miniatures events that we already demonstrated in mice.43

Figure 6. Trem2–/– hippocampal neurons display synaptic and functional alterations in a region-dependent manner
(A) Electron micrograph of synapses in the CA1 pyramidal layer of P18 Trem2+/+ and mice. Scale bars, . Box-and-whiskers plot of synaptic mitochondria length of (Trem2+/+) or 18 (Trem2 $ / -$ ) synapses from Trem2+/+ and Trem2 P18 mice. Mann-Whitney U test, .
(B) Synaptic vesicles (SV) density, synaptic area, docked SV density, and active zone length in CA1 pyramidal layer synapses of P18 Trem2+/+ and Trem2−/− mice. Box-and-whiskers plots of n = 56 (Trem2+/+), or 59 (Trem2−/−) synapses from active zone length in Trem2+/+ and Trem2−/− mice. Mann-Whitney U Test, p < 0.01 and **p < 0.001.
© Left: electron micrographs of synapses in the CA3 pyramidal layer of P18 Trem2+/+ and Trem2−/− mice. Scale bars, 200 nm. Right: box-and-whiskers plots of synaptic mitochondrial length, docked SV density, or active zone length in CA3 pyramidal layer synapses of P18 Trem2+/+ and Trem2−/− mice. Mann-Whitney U test. (D) SV density, synaptic area, docked SV density, (active zone) synapses from N = 3 Trem2+/+ and Trem2−/− P18 mice. Mean ± SEM. Box-and-whiskers plots of n = 58 (Trem2+/+) or 65 (Trem2−/−) synapses from N = 3 Trem2+/+ and Trem2−/− mice. Mann-Whitney U test.
(E) Up: image of a coronal hippocampal slice from a P18 Trem2+/+ mouse on the CMOS multi-electrode array. Scale bars, 500 μm. Down: raw data, filtered signal (350 and 0 Hz), and spike detection originating from an electrode.
(F, G) Cumulative frequency distributions and quantification of average spontaneous firing rates from neurons in CA1 (F) and CA3 (G) subfields of Trem2+/+ and Trem2−/− P18 mice. Box-and-whiskers plots of n > 600 neurons from n = 6 hippocampal slices from 3 mice. Dots represents individual firing units. Unpaired Student’s t test, **p < 0.0001.
(H) Cross-correlation matrix of neuron pairs during CA1 and CA3 spontaneous activity in Trem2+/+ and Trem2−/− P18 male.
(I) Average correlation coefficient from all neuron pairs in CA1 and CA3 regions in Trem2+/+ and Trem2−/− mice. Box-and-whiskers plots of n = 6 slices from N = 3 P18 Trem2+/+ and Trem2−/− mice. One-sided unpaired Student’s t test, p < 0.05.
(J, K) Confocal images of CA1 (J) and CA3 (K) subfields of P90 Trem2+/+ and Trem2−/− hippocampi labeled with Ab against PSD95. Scale bars, 10 μm.
(L and M) PSD95 volume fraction in CA1 (L) and CA3 (M) regions. Box-and-whiskers plots of N = 6 Trem2+/+ and N = 6 Trem2−/− P90 mice. Unpaired Student’s t test.
See also Figure S6.
We have previously shown that at more mature developmental stages (P90), the hippocampal amount of PSD95 protein in hippocampus undergoes a reduction, compared with P18, in the absence of changes in spine density.43 By specifically analyzing different hippocampal regions, we now demonstrate that the synaptic impoverishment at P90 occurs specifically in the CA1 but not the CA3 region of mice (Figures 6J–6M). These data suggest that the metabolic faultiness in the CA1 region may cause defective synapse strengthening in the adult (P90).
We next wondered whether the metabolic derangement detected in neurons from brains required the continuous proximity of microglial cells or relied on cell-autonomous mechanisms. To address this question, primary neurons were harvested from P1 hippocampi of and mice, maintained in culture for 4 days (4 days in vitro [DIV]) and then assayed for OCR measurements (Figure 7A). Results showed that the mitochondrial fitness of neurons is significantly altered (Figure 7B). Indeed, 4 DIV Trem2 neurons showed an impairment to meet their energetic demand (Figure 7C). Non-mitochondrial respiration was also defective in neurons (Figure 7C). By virtue of the transcriptional signature emerging from scRNA-seq, involving most of genes linked to OXPHOS and mitochondrial dysfunction pathways, these data indicate that the functional defects in metabolic fitness have been transcriptionally imprinted during earlier developmental times and that this signature intimately predetermines the mitochondrial functional outcome. To parallel the neuronal metabolic derangement observed with extracellular flux measurements in 4 DIV cultures, a significant reduction of Tom20 (Figure 7D) was detected in both the soma and neurites of neurons, compared with (Figures 7E and 7F). Furthermore, a reduction in mitochondrial ATP production at the neuronal soma was detected by live-imaging experiments (Figure 7G). Also, Sholl analysis (Figure 7H) provided evidence of an overall reduction in neuronal ramification (Figure 7I) and in primary neurite length (Figure 7J), supporting the hypothesis of a developmentally delayed neuronal maturation.
Given that the composition of and channels significantly changes along with the maturation of neuronal cells,66 an electrophysiological characterization of 4 DIV neurons was performed. We recorded current-voltage relationships to estimate inward and outward currents (Figure 7K). The curve demonstrated a significant reduction in but not in current density in neurons, compared with neurons (Figures 7L and 7M). Moreover, neurons displayed a more depolarized resting membrane potential than Trem2+/+ neurons (Figure 7N), which is generally considered a feature of neuronal immaturity.67,68 Finally, neurons showed slightly higher electrical resistivity and lower capacitance values than Trem2+/+ (Figure 7O). These data indicate that even in primary culture, neurons from Trem2/ hippocampi display defective metabolic properties and a delay in maturation, paralleling what was observed in the in vivo setting (Figures 1A–1C, 3A–3D, and 3F–3H).
We next asked whether the excitatory neurotransmission derangement that was observed in vivo at P18 (Figures 6E–6I) was also detectable in mature cultures of neurons from hippocampi. We have previously shown that the electrophysiological properties of 14 DIV neurons from or mice are not different.43 Also, neuronal resting potential and passive properties (resistivity; capacitance) were similar among the two populations, in both CA1 and CA3 subfields of P18 slices (Figures S6G and S6H) and 14 DIV cultured neurons (Figures S7D and S7E). However, and similarly to the ex vivo setting, mitochondria remained defective in neuronal somata and synaptic contacts of 14 DIV neurons from Trem2 mice (Figures S5A–S5C). Furthermore, live-imaging recordings at synaptic boutons of 14 DIV Trem2/ synapses, using the glutamate sensing fluorescent reporter (iGluSNFR), revealed a significant increase in glutamate release, compared with control cultures (Figures 7P and 7Q). To evaluate the dependency of these synapses on mitochondrial ATP, we recorded, from the very same synapses, the fluctuation in mitochondrial ATP in response to high-frequency stimulation and the synaptic ATP fluctuations using ATP-red dye and found that synapses fail to produce ATP in response to high-frequency stimulation (Figure 7R), in line with reduced mitochondrial functionality.

Figure 7. Trem2–/– neurons display mitochondrial defects and faulty differentiation in culture
(A) Cartoon showing the experimental design.
(B) Average OCR kinetics showing the response of 4 DIV Trem2+/+ or primary hippocampal neurons to different mitochondrial inhibitors.
© From left to right: basal, maximal, ATP-coupled respiration, and non-mitochondrial oxygen consumption rates. Mean ± SEM of n = 25 replicates from n = 3 independent Trem2+/+ and Trem2−/− cultures. One sample t test, **p < 0.0001.
(D) Confocal images of 4 DIV Trem2+/+ or Trem2−/− primary hippocampal neurons labeled with Abs against Tom20 (fire look-up table [LUT]) and MAP2 (red). Scale bars, 8 and 2 μm.
(E and F) Tom20 mean intensity at somata (E) and primary neurites (F), respectively. Mean ± SEM of n = 40 (Trem2+/+) and n = 40 (Trem2−/−) neurons from N = 6 independent Trem2+/+ and Trem2−/− neuronal cultures. Unpaired t test, **p < 0.0001.
(G) Representative ATP-red fluorescence variation upon pharmacological inhibition of mitochondrial complex 1 and 3 (rotenone/Aa) in 4 DIV Trem2+/+ neuronal somata. Scale bars, 2 μm. Down left: average ATP-red traces at soma (ΔF/F₀). Data are presented as mean ± SEM. Down right: plots showing ATP reduction after 15 min from the treatment (ΔF/F₀). Box-and-whiskers plots of n = 13 (Trem2+/+) and n = 14 (Trem2−/−) neurons from N = 4 Trem2+/+ and Trem2−/− neuronal cultures. Unpaired Student’s t test, p < 0.05.
(H) Confocal images of 4 DIV Trem2+/+ or Trem2−/− primary hippocampal neurons labeled with Ab against MAP2 and processed for Sholl analysis. Scale bars, 20 μm.
(I) Sholl analysis of 4 DIV Trem2+/+ or Trem2−/− primary hippocampal neurons. Mean ± SEM of n = 54 (Trem2+/+) and n = 41 (Trem2−/−) neurons from N = 5 independent Trem2+/+ or N = 4 Trem2−/− primary neuronal cultures. Two-way ANOVA followed by Sidak’s multi-comparison test, ***p < 0.0001.
(J) Primary neurite length in 4 DIV Trem2+/+ or Trem2−/− primary hippocampal neurons. Box-and-whiskers plots of n = 54 (Trem2+/+) and n = 41 (Trem2−/−) neurons from N = 5 independent N = 5 Trem2+/+ or N = 4 Trem2−/− primary neuronal cultures. Unpaired t test, p < 0.05.
(K and L) Cartoon of the experimental design.
(M) Mean inward Na+ and outward K+ currents of 4 DIV Trem2+/+ or Trem2−/− primary hippocampal neurons. n = 5–10 (Trem2+/+) and n = 6–9 (Trem2−/−) neurons from N = 3 Trem2+/+ or Trem2−/− primary neuronal cultures. Two-way ANOVA, **p < 0.001.
(N) Resting potential of 4 DIV Trem2+/+ or Trem2−/− primary hippocampal neurons. Box-and-whiskers plots of n = 14 (Trem2+/+) and n = 16 (Trem2−/−) neurons from N = 4 Trem2+/+ or Trem2−/− primary neuronal cultures. Mann-Whitney U test, p < 0.05.
(O) Resistivity (left) and capacitance (right) values of 4 DIV Trem2+/+ or Trem2−/− primary hippocampal neurons. Box-and-whiskers plots of n = 9 (Trem2+/+) and n = 10 (Trem2−/−) neurons from N = 3 Trem2+/+ and N = 2 Trem2−/− neuronal cultures. Unpaired t test.
§ Up: description of the iGluSNFR live-imaging experimental design. Down: representative iGluSNFR fluorescence images in a transfected Trem2+/+ neuron. Zoomed image of an axonal branch with synaptic boutons after a 10-Hz stimulation of 50 APs. Scale bars, 20 and 5 μm. Left: average traces of iGluSNFR response to AP Trains at 10 Hz in Trem2+/+ and Trem2−/− synapses. Mean ± SEM of n = 22 (Trem2+/+) and n = 21 (Trem2−/−) coverslips from N = 5
To start evaluating whether the defective metabolic phenotype of neurons developing in the absence of microglial Trem2 could result from the early neuronal exposure to an altered microglial secretome,69–72 pre-conditioned medium (CM) harvested from either or primary microglia was applied to primary neurons isolated from or hippocampi and maintained for 4 DIV (Figure S7F). Seahorse experiments (Figure S7G) revealed defective metabolic parameters in as compared with neurons. However, released soluble factors did not boost any of the metabolic parameters (Figures S7H–S7K), with a very modest, not significant, effect of the supernatant, regardless of whether it was derived from or microglia.
Taken together, these results indicate that the lack of Trem2- mediated communication between microglia and neurons during early developmental windows is sufficient to derange the forthcoming neuronal metabolism and mitochondria organization, which is maintained in neurons even after their isolation from the brain context. These defects are associated with a delay in the maturation of neurons, followed, both ex vivo and in vitro, by dysregulated glutamatergic neurotransmission. Furthermore, the secretome from microglia is not able to revert the defective metabolic profile of neurons, suggesting that more structured bidirectional microglia-neuron communication may be required.
Microglia are highly dynamic cells that patrol brain parenchyma, sensing the neuronal activity and regulating neuronal functions through the release of a wide spectrum of immunoregulatory factors. Through secretory and cell-to-cell contact-mediated processes, microglia dynamically regulate neuronal development, growth, formation of synapses, and physiological functions of brain cells.1,73 We have previously demonstrated that lack of the immune receptor Trem2 results in microglia’s inability to properly perform supernumerary synapse elimination during brain development.43 Also, we have shown that adult (P90) Trem2-deficient mice display altered brain connectivity and behavioral defects.43,45 We now demonstrate that Trem2 deficiency impairs neuronal transcriptomic and energetic profile, suggesting that early metabolic derangements in neurons may contribute to connectivity and behavioral defects occurring in adult mice.
Lack of Trem2 strongly impacts the expression of EIF2 signaling in neurons, which is heavily downregulated in CAImmPyr and upregulated in CA1-Pyr and CA3-Pyr. By controlling the eIF2.GTP.Met-tRNAiMet ternary complexes, EIF2 activity governs the rate of global translation, and its phosphorylation is recognized as implicated in nervous system development and memory consolidation.74,75 Indeed, genetic reduction of Elf2b activity induces delayed post-natal brain development, abnormal glial cell abundance, increased abundance of demyelinated axons, and impaired cerebral inflammatory response upon insults.76,77
In this study, we provide several lines of evidence indicating that mitochondrial and energetic assets of CA1 neurons are profoundly altered in the absence of microglial Trem2. Microglia react to the environment by flexibly reprogramming their metabolic pathways, in a process called immunometabolism.78,79 Indeed, while homeostatic microglia mostly rely upon OXPHOS for ATP production, inflammatory microglia reprogram their metabolism toward aerobic glycolysis.80–82 Further, Trem2 plays a key role in modulating cellular biosynthetic metabolism16 and is centrally involved in inducing transcriptomic and functional programs of DAMs, a phagocytic state detected in the brain during neurodegeneration.26,33
We found that in the early post-natal period microglia equally rely on glycolysis and OXPHOS for ATP supply. We also showed that, in line with the central role of Trem2 in the control of microglia metabolism,16 the lack of the gene causes, at birth, a reduction in mitochondrial mass without, however, major changes in the OCR. Conversely, hippocampal neurons from mice displayed impaired mitochondrial metabolism, with selective deterioration in the basal, maximal, and ATP-coupled respiration, accompanied by diminished mitochondrial mass and reduced mitochondria length. Thus, the lack of microglial Trem2 predominantly affects neuronal metabolic properties.
Differently from microglia, neurons are allowed only a limited degree of metabolic flexibility, which is required to ensure fate stability and sustained cellular functionality. During development, a continuous metabolic glycolytic rewiring of stem and precursor cells allows cell fate regulation, proliferation, and differentiation processes. Conversely, mature neurons, to prioritize network maintenance, mostly rely on OXPHOS and mitochondria, with a very minor proportion of their energy deriving from glycolytic processes83 deputed to the control of processes such as synaptic vesicle trafficking.84 It has been proposed that glycolysis or certain glycolytic metabolites are incompatible with neuronal identity and survival,83 causing neural cell dedifferentiation, loss of cell fate stability, defects in synaptic transmission, and ultimately, cell death. In line with the metabolic switch from glycolysis to OXPHOS representing an essential step in neuronal maturation, defects in the OXPHOS genes contribute toward neurodevelopmental disorders.85,86 The neuronal reduction in OXPHOS observed in mice likely represents an impasse to proper neuronal and synapse maturation and possibly even to brain circuit development.
We showed that the delayed CA1 pyramidal neuron maturation in mice is followed at P18 by the enhancement of excitatory neurotransmission. A similar phenomenon has been reported in different neurodevelopmental disorders, such as Rett syndrome, Down syndrome, and epileptic encephalopathies, where brain energy dysfunctions are associated with neuronal developmental delay and followed by an increase in excitatory transmission,87,88 correlating with enhanced network synchronicity.65,89 Furthermore, the ATP produced by glycolysis, but not mitochondria, is harnessed for glutamate accumulation into synaptic vesicles.90 These data provide a logical frame to our observation that neurons, which are defective in mitochondrial respiration, release a higher amount of glutamate upon stimulation. At later stages of brain maturation (P90), mice display an impoverishment of hippocampal synaptic markers,43 which specifically occurs in the CA1 region. This event associates with altered functional connectivity and significant defects in circuit maturation.43,45 Thus, mice seem unable to properly build and strengthen the correct synaptic setup during brain maturation, as a possible result of the neuronal metabolic impairment as well as the defective synapse elimination during the period of synaptic pruning.
Ample evidence has recently implicated mitochondria in early brain development. Mitochondria function is associated with neurogenesis, differentiation, and maturation of neurons. During mid-corticogenesis, mitochondria assembly, fusion/fission, and mobility dynamics are key processes that warrant progenitors’ self-renewal capability and maintenance.57–59 Alterations of the timing in the fusion/fission paradigm affect neurogenesis,91 and defects in mitochondria mobility and energy metabolism selectively impinge on neuronal differentiation and migration of cortical interneurons92 Moreover, mitochondrial biogenesis is critical for spine formation and maintenance in hippocampal neurons.93 To our knowledge, only cell-autonomous mechanisms have been so far identified in underpinning neuronal metabolic programs. In this study, we show a role for microglia in controlling neuronal maturation and metabolic fitness and identify Trem2 signaling as an early potential modulator of such microglia-neuron crosstalk.
Beyond DG neurons, we detected a selective vulnerability of specific classes of neurons and identify CA1 neurons as one of the neuronal subtypes most susceptible to the lack of Trem2. Consistently, mitochondrial defects and metabolic derangements occur specifically—and are even exacerbated at later (P18) developmental stages—in the CA1 and not the CA3 region. Consistently, assessment of the neuronal developmental trajectories in newborn mouse hippocampus by pseudotime analysis revealed a clear shift in the metabolic signature of CA1 neurons, associated with a delay in neuronal maturation. We have previously shown that also synaptic pruning depends on Trem2 predominantly in the CA1 as compared with the CA3 region.43 Thus, several Trem2 functions are achieved in a region-dependent manner. Although the basis for the selective susceptibility of the CA1 region to Trem2-depletion is not fully defined, it is relevant that CA1 microglia express higher Trem2. Furthermore, even partial reductions of Trem2 expression, as in Trem2 hemizygous mice, are sufficient to derange CA1 pyramidal neuron metabolism, indicating that the CA1 requires high microglial Trem2 to properly control pyramidal neuron differentiation. These data have important implications, opening the possibility that conditions characterized by reduced Trem2, including Trem2 hemizygous missense variants in patients,34,38,94 may suffer from neuronal metabolic dysfunctions.
How Trem2 modulates microglia-neuron communication remains to be defined. One possibility is that the lack of Trem2 and consequent modification of the microglial state could result in altered secretome, which in turn would differentially affect neuronal trajectories. Alternatively, or concomitantly, microgliato-neuron contacts occurring during early development may play a role. Specifically, Trem2 localized at the microglial membrane could bind to immature neuronal partners, influencing their metabolic trajectories. It has been found that throughout embryonic and adult neurogenesis, microglial processes can form specialized contacts with the cell bodies of developing neurons, characterized by the clustering of functional microglial P2yr12 receptors and enriched with mitochondria.95 In the absence of microglial P2yr12, erratic cortical cytoarchitecture occurs and is maintained in adulthood.95,96 The selective effect on excitatory cells, which are born locally in the hippocampus, where all the developing stages are potentially in contact with instructive microglia, is in line with the possibility of a cell-to-cell contactdependent mechanism. Our pseudotime analysis suggests that this putative interaction, occurring at early stages of maturation of the neuronal excitatory lineage, may exert a potential ‘‘priming’’ effect, possibly differently in CA1 vs. CA3 regions, where Trem2 is expressed at a different extent.
At present, we do not know whether microglia control the neuronal metabolic maturation only during development or are also active in different temporal windows. Neuronal changes in glucose metabolism occur during aging, when senescent neurons exhibit mitochondrial dysfunctions and increase their capability to oxidize glucose via glycolysis.97 A reduction in Trem2 expression also occurs during senescence,98 suggesting that these two mechanisms might contribute to the neuronal metabolic switch from OXPHOS to glycolysis during aging. Indeed, mouse brains analyzed by fluoro-2-deoxy-d-glucose micro positron emission tomography (FDG- ) display mreduced glucose metabolism,99 while a knockin mouse model for the disease-associated Trem2 T66M mutation presents a significant reduction in both microglial activity and overall brain glucose metabolism.38 Although the total FDG uptake measured by PET is the composite of neuronal, astroglial, and microglial cells’ metabolic profiles, it is possible that a crosstalk, like that occurring between microglia and neurons during early development, might be operative also during aging, and whenever compromised, it contributes to neurodegeneration. In this respect, it is particularly relevant that one of the functional consequences of the lack of Trem2 is the increase in glutamate release and spontaneous neuronal activity. It is known for many years that the hippocampus is importantly affected in and that an increase in hyperactive neurons in the CA1 region occurs before the formation of plaques .101,102 Also, thinning of the CA1 apical neuropil occurs in subjects with mild AD, compared with controls,103 while at later stages the CA1 region suffers more severe synaptic alterations.104 Importantly, we show that adult (P90) Trem2/ mice are characterized by a significant synaptic impoverishment, which occurs selectively in the CA1 region and accompanies the connectomic and behavioral deficits already described.43,45 Thus, the possibility that a faulty control of neuronal metabolism by defective microglial Trem2 results in deranged neuronal activity also during aging and neurodegeneration is a scenario worth investigating.
Although revealing a new role of TREM2 in early neurodevelopment, our study still lacks in identifying the molecular mechanisms by which Trem2-bearing microglia control the metabolic neuronal fate. More detailed analyses should be performed, addressing the possible roles of sTrem2 or microglia-derived extracellular vesicles105 and identifying the molecular code of the microglia-to-neuron communication, including whether CA1 and CA3 regions differently express putative Trem2 interactors. Furthermore, experiments of microglial Trem2 expression downregulation in early post-natal windows would univocally prove that the synaptic and circuits defects observed in the adult mouse result from a microglia-to-neuron crosstalk dysregulation occurring during the early development. Moreover, since the lack of Trem2 impairs both synapse refinement43 and neuronal metabolism (this study), it would be important in the future to define next at which extent the two processes relatively contribute to the synaptic setup during brain maturation.
Detailed methods are provided in the online version of this paper and include the following:
d KEY RESOURCES TABLE
RESOURCE AVAILABILITY
Lead contact Materials availability Data and Code availability
Mice
Primary hippocampal neuronal culture and transfection
Primary microglia culture
Live-imaging experiments
RNAseq analysis
Electrophysiology
Immunofluorescence analysis
Electron microscopy
Supplemental information can be found online at https://doi.org/10.1016/j.
immuni.2023.12.002.
We thank Prof. K. Volynski (Department of Experimental and Clinical Epilepsy) for kindly providing iGluSnFR and GFP plasmids. We thank C. Ferrari, D. Pozzi, E. Fraviga, C.A. Elia, and C. Saulle for help in setting up experimental methodologies. The present work was supported by ERC AdG MATILDA 101055323, EraNET Neuron JTC2021 InflASD and Telethon GGP20030 to M.M., by ERC StG IMPACT 101043003, Cariplo Giovani 2019-1785 to S.L., by a Grant D.R. 3404 on Heavy Equipment, HITACHI 120 kV TEM microscope HT7800 to K.C. E.T. and A.C. are supported by a High Profile PostDoctoral (HiPPO) Fellowship by Fondazione Humanitas per la Ricerca. Cartoons and GA have been created with Biorender.com.
E.T. performed primary neuronal and microglia cultures; performed and analyzed in vitro live-imaging experiments, ex vivo electron microscopy, and in vitro and ex vivo immunofluorescence experiments; and assisted in sample collection. G.D. performed primary neuronal cultures, performed and analyzed in vitro and ex vivo Seahorse experiments, and assisted in sample collection. M.B. performed primary microglia cultures, immunofluorescence experiments, conditioned media preparation, and sample collection. E.F. performed qPCR, colony maintenance, genotyping, and sample collection. A.C., R.J.A., and G.D. designed and A.C. performed SCENITH experiments. F.F. and A.C. set up and performed FACS experiments. S.M., M.M., and S.L.
performed, analyzed, and interpreted the single-cell sequencing data. M.C.G. and K.C. performed and analyzed ex vivo electron microscopy. R.M. performed in vitro electrophysiological recordings, immunofluorescence experiments, and assisted in sample preparation. R.H.-S. performed MEA recordings, and P.P. performed ex vivo immunofluorescence experiments. M.M., E.T., G.D., and R.M. conceived the study, designed the experiments, and analyzed the data. M.M., E.T., G.D., and S.L. wrote the original draft of the manuscript, and all authors revised the manuscript.
The authors declare no competing interests.
Received: September 15, 2022
Revised: July 17, 2023
Accepted: December 5, 2023
Published: December 29, 2023
KEY RESOURCES TABLE
| REAGENT or RESOURCE | SOURCE IDENTIFIER | |
| Antibodies | ||
| Monoclonal mouse anti-Tom20 | Santa-Cruz | Cat#sc-17764 |
| Polyclonal rabbit anti-P2YR12 | Anaspec | Cat#AS-55043A |
| Polyclonal guinea pig anti-IBA1 | Synaptic Systems | Cat#234 308 |
| Polyclonal rabbit anti-IBA1 | WAKO Chemicals | Cat#019-19741 |
| Polyclonal guinea pig anti-DCX | Merck-Milipore | Cat#AB2253 |
| Polyclonal guinea pig anti-MAP2 | Synaptic Systems | Cat#188 004 |
| Polyclonal rabbit anti-vGlut1 | Synaptic Systems | Cat#135 302 |
| Monoclonal rat anti-CD68 | Biolegend | Cat#137001 |
| Polyclonal rabbit anti-Neurod2 | Abcam | Cat#ab104430 |
| Polyclonal rat anti-Ctip2 | Abcam | Cat#ab18465 |
| Monoclonal rat anti-CD45 Pe | Biolegend | Cat#103106 |
| Monoclonal rat anti-CD45 PerCP | Biolegend | Cat#103130 |
| Monoclonal rat anti-CD11b Pe-Cy7 | Biolegend | Cat#101216 |
| Monoclonal rat anti-CD16.32 | Biolegend | Cat#101301 |
| Monoclonal mouse anti-NeuN 488 | Merck-Milipore | Cat#MAB377X |
| Anti-Puromycin 647 | Scenith -Auguello Lab | RRID:AB_2827926 |
| Polyclonal rabbit anti-PSD95 | Thermo-Fisher | Cat#51-6900 |
| Polyclonal rabbit anti-Trem2 | R&D Systems | Cat#AF1729 |
| Goat anti-Mouse IgG (H+L) Highly Cross-Adsorbed Secondary Antibody, Alexa Fluor 568 | Thermo-Fisher | Cat#A11031 |
| Goat anti-Rabbit IgG (H+L)Highly Cross-Adsorbed Secondary Antibody, Alexa Fluor 488 | Thermo-Fisher | Cat#A11034 |
| Goat nti-Rabit IgG (H+L) Highly Cross-AdsorbedSecondary Antibod Alexa Fluor 633 | Thermo-Fisher | Cat#A21071 |
| Goat anti-Rat IgG (H+L) Highly Cross-Adsorbed Secondary Antibody, Alexa Fluor 647 | Thermo-Fisher | Cat#A21247 |
| Goatanti-Guinea PigIgG(H+L)Highly Cross-Adsorbed SecondaryAntibody Alexa Fluor 488 | Thermo-Fisher | Cat#A11073 |
| Goat anti-Guinea Pig IgG (H+L) Highly Cross-Adsorbed SecondaryAntibody Alexa Fluor 568 | Thermo-Fisher | Cat.#A11075 |
| Goat anti-Guinea Pig IgG (H+L) Highly Cross-Adsorbed SecondaryAntibody | Thermo-Fisher | Cat#A21105 |
| Alexa Fluor 633 Donkey anti-Sheep IgG (H+L) Highly Cross-Adsorbed Secondary Antibody, | Thermo-Fisher | Cat#A-11015 |
| Alexa Fluor 488 Donkey anti-Rabbit IgG (H+L) Highly Cross-Adsorbed Secondary Antibody, | Thermo Fisher | Cat#A-10042 |
| Alexa Fluor 568 Fluorsave | Millipore | Cat. #345789 |
Continued
| Chemicals,peptides,and recombinant proteins | ||
| CNQX disodium salt | Tocris | Cat#1045/1 |
| D-AP5 | Tocris | Cat#0106/1 |
| FCCP | Agilent Technologies | Cat#103015-100 |
| Oligomycin | Agilent Technologies | Cat#103015-100 |
| Rotenone/Antimycin A | Agilent Technologies | Cat#103015-100 |
| Puromycin | Scenith/Auguello Lab | N/A |
| 2-Deoxy Glucose | Scenith/Auguell Lab | N/A |
| Oligomycin | Scenith/Auguello Lab | N/A |
| Neuromag reagent | OZ-Biosciences | Cat#KC30800 |
| Epoxy resin Poly/Bed@ 812 | Polysciences | Cat#08792-1 |
| 1% Uranyl Acetate | Electron Microscopy Sciences (EMS) | Cat#22400-1 |
| 4% Osmium Tetroxide | Electron Microscopy Sciences (EMS) | Cat#19140 |
| Propylene Oxide | TAAB | Cat#P021 |
| Hoechst-33342 | Thermo-Fisher | Cat#62249 |
| BioTracker ATP-Red Live Cel Dye | Merck-Milipore | Cat#SCT045 |
| MitoTracker Deep Red FM | Thermo-Fisher | Cat#M22426 |
| CellRosGreen reagent | Thermo Fisher | Cat#C1044 |
| Anti-MouseIgegativeCtro Compensation Particles Set | BD Biosciences | Cat#552843; RRID: AB_10051478 |
| Anti-Rat Ig, K/Negative Control | BD Biosciences | Cat#552844; |
| Compensation Particles Set Critical commercial assays | RRID: AB_10055784 | |
| Seahorse XF Cel Mito Stress Test it | Agilent Technologies | Cat#103015-100 |
| Fixation/Permeabilization Solution Kit | BD Bioscience | Cat#554714 |
| Foxp3/Transcription Factor Staining Buffer Set | Thermo-Fisher | Cat#00-5523-00 |
| Zombie NIRTm Fixable Viability Kit | Biolegend | Cat#423105 |
| Chromium Next GEM Single Cell3 Kt v3.1 | 10X Genomics | |
| Chromium Next GEM Chip G Single Cell Kit | 10X Genomics | Cat# PN-1000269 |
| Dual Index Kit TT Set A | 10X Genomics | Cat# PN-1000127 |
| Direct-ZolTm MiniPrep Isolation Kit | Cat# PN-1000215 | |
| High-Capacity cDNA RT kit | Zymo Research | Cat# R2050 |
| TaqMan Fast Universal PCR Master | Applied Biosystems | Cat# 4368814 |
| Mix(2x), no AmpErase UNG | Applied Biosystems | Cat# 4352042 |
| Deposited data | ||
| scRNA-seq data | This Paper | Geo number: GSE249036 |
| Experimental models: Organisms/strains | ||
| Mouse C56BL/6J | Charles River | Strain Code 632 |
| Mouse: Trem2-/- | Turnbull et al.15 | N/A |
| Recombinant DNA | ||
| pFU_GFP | Kirill Volynski, UCL | N/A |
| pAAV.hSynap.SF-iGluSnFR.A184V | Kill Volynski, UCL | Cat#106175; RRID: Addgene_106175 |
| Software and algorithms | ||
| GraphPad Prism 9 | GraphPad | https://www.graphpad.com; RRID:SCR_002798 |
| Fiji | NIH | https:/fij.sc/; RID:SCR_002285 |
| Adobe Ilustrator CC | Adobe | http://www.adobe.com/products/illustrator.html RRID:SCR_010279 |
| Radius 2.0 software | EMSIS | https://www.emsis.eu/products/radius |
Continued
| REAGENT or RESOURCE | SOURCE | IDENTIFIER |
| Seahorse Wave Version | Agilent Technologies | http://www.agilent.com/en-us/products/ cell-analysis-(seahorse)/software- download-for-wave-desktop; RRID:SCR_014526 |
| MetaFluor Fluorescence Ratio Imaging Software | Molecular Devices | RRID:SCR_014294 |
| Matlab R2018 or R2021b | Mathworks | https://t.mathworks.com/; RRID: SCR_001622 |
| FlowJo 10 | FlowJo | https://www.flowjo.com/; RRID: SCR_008520 |
| Bcl-convert 3.8.4 | Bcl-convert | https://emea.support.ilumina.com/ sequencing/sequencing_software/ bcl-convert.html |
| 10x Genomics Cell Ranger 6.1.1 | Cellranger | https://www.10xgenomics.com/; RRID:SCR_023221 |
| R 4.2.2 | R | https://ww.r-project.org/; RRID: SCR_002394 |
| Seurat 4.2.0 | Seurat | https://satijalab.org/seurat/articles/ get_started.html; RRID:SCR_016341 |
| Bioconductor Milor 1.2.0 | MiloR | https://www.bioconductor.org/packages/ release/bioc/html/miloR.html |
| Slingshot 2.6.0 | Slingshot | https://www.bioconductor.org/packages/ release/bioc/html/slingshot.html106 |
| pClamp-10 | Axon Instruments, Foster City, CA | https://www.moleculardevices.com; RRID:SCR_011323 |
| BrainWave 4 | 3Brain AG, Switzerland | https://www.3brain.com/products/ software/brainwave4 |
| Biorender | Biorender | http://biorender.com: RRID:SCR_018361 |
| Viia7 software system | Applied Biosystems | http:https://www.thermofisher.com/it/ en/home/life-science/pcr/real-time-pcr/ real-time-pcr-instruments/viia-7-real- |
| Imaris v. 9.7.2 or 7.2.3 | Oxford Instruments | time-pcr-system.html RRID:SCR_023358 http://www.bitplane.com/imaris/imaris RRID:SCR_007370 |
| QuPath v0.3.2 | QuPath | https://doi.org/10.1038/ s41598-017-17204-5 |
Further information and requests for resources and reagents should be directed to and will be fulfilled by lead contact Michela Matteoli (michela.matteoli@hunimed.eu).
This study did not generate new unique reagents.
d Single-cell RNA-seq data have been deposited at GEO and are publicly available as of the date of publication. Accession numbers are listed in the key resources table.
d This paper does not report original code.
d Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.
All experiments were performed in accordance with the guidelines established by the European Community Council (Directive 2010/ 63/EU of September 22nd,2010) and the Italian D.Lg. 26/2014. The study was approved by the Institutional Animal Care and Use Committee (IACUC) of the Humanitas Research Hospital and by the Italian Ministry of Health. C57BL/6J Trem2-/- mice, generated as previously described15 were provided by Bioxell-Cosmo Pharmaceutical (Milan, Italy).107 C57BL/6J mice were acquired by Charles River laboratories. Mice were housed in a Specific Pathogen Free (SPF) facility under constant temperature and humidity conditions with a 12 hrs- light/dark cycle and were provided with food and water ad libitum. P0-1 and P18-20 male and female animals were used for each experiment, unless otherwise stated.
Primary hippocampal neurons were isolated from either or P0-1 mice15 and cultured in complete Neurobasal A medium (cat.10888, Thermo-Fisher) supplemented with B27 (cat.17504044,Thermo Fisher Scientific). Briefly, hippocampi were dissected and dissociated by enzymatic digestion in trypsin for at and then triturated using a standard micropipette. Neurons were plated on poly-L-lysine ; cat. P2636, Sigma-Aldrich)–coated or diameter glass coverslips at a density of 100,000 to 120,000 cells per coverslip (for imaging and electrophysiological experiments) or on poly-L-lysine treated cartridges at a density of 12000 cells/well on XF-96 cell culture plates for Seahorse Extracellular Flux Analysis. Experiments were performed at 4-5 or 14-16 days in vitro (DIV). For both live imaging and immunofluorescence staining, hippocampal neurons were transfected at 5 DIV with either pFU_GFP or pAAV.hSynap.SF-iGluSnFR.A184V plasmid using Neuromag reagent (cat.KC30800; OZ Biosciences). This allowed expression of the iGluSnFR probe only in a small to ) subpopulation of neurons, which was essential for imaging of vesicular release in individual synaptic boutons. The pFU_GFP and the pAAV.hSynap.SF-iGluSnFR.A184V plasmids were kindly provided by K. Volynski, University College London, London, UK.
Primary microglia cultures were obtained from mixed glial cultures as previously described43,44). Briefly, meninges were removed from P2 newborn brains, then cortices and hippocampi were gently dissected in cold Hanks Balanced Salts Solution (HBSS, cat.14180046, Thermo Scientific). After digestion with trypsin at and mechanical dissociation, cells were plated in T-75 flasks and maintained in EMEM (cat.M5650, Sigma-Aldrich) containing of heat-inactivated Fetal Bovine Serum (FBS, cat.ECS0186L, Euroclone), glucose, 2mM sodium pyruvate (cat.ECM9542D, Euroclone) and penicillin/streptomycin (cat.15140-122, Thermo Fisher). Microglia were isolated around DIV 15 by shaking flasks at 245 rpm speed for 60 minutes.
For conditioned medium preparation, microglial cells were plated on poly-L Ornithine (cat. P8638, Sigma Aldrich) pre-coated 6-multiwell in Neurobasal-A/B27 medium at a density of cells/well. After 24 hours in vitro, Conditioned Medium (CM) was collected and centrifuged for 20 minutes, then used for subsequent experiments.
4DIV primary hippocampal neurons obtained from and newborn pups were loaded with the mitochondrial ATP biotracker ATP-red live cell dye (cat.SCT045, Merk-Millipore) for at in Neurobasal-A Medium and then imaged for measuring mitochondrial ATP production. Somatic neurons were identified by either sparse GFP transfection or by bright-light images of neuronal soma. Neurons were maintained in the extracellular buffer 1. The extracellular buffer (EB1) contained , 2.5 mM KCl, 2 mM MgCl2, 2 mM CaCl2, glucose, and Hepes (pH 7.4).108,109 After 2 mins recording of basal ATP-red fluorescence, mitochondrial complex 1 blocker Rotenone plus Antimycin A , cat.103015-100, Agilent Technologies) was mapplied to the cells and ATP decay measured for an additional 10 minutes. Time-lapse recording of ATP dynamics was performed with an acquisition rate of . A 18-mm recording chamber (cat.QR-40LP, Warner Instruments, Hamden, CT) was placed on the stage of an IX-71 inverted microscope (Olympus, Hamburg, Germany) equipped with an EMCCD (electron-multiplying CCD) camera (QuantEM-SC 512x512, Photometrics). Illumination was obtained using a light-emitting diode LED (Cairn research, Optoled Lite), coupled to the acquisition software (MetaFluor, Molecular Devices) to minimize photobleaching. A 40X (1,3 numerical aperture; N.A.) oil-immersion objective was used to identify multiple neuronal somas in the field of view pixel). ATP decays were measured as F/F0. Frotenone was then measured to determine mitochondrial ATP production. All values were normalized Dto Trem2+/+ neurons.
14-16 DIV hippocampal neurons obtained from and newborn pups and transfected with iGluSNFR probe were incubated for 45 mins at with the mitochondrial ATP biotracker ATP-red live cell in Neurobasal-A medium. Neurons were first stimulated with 5 stimuli at and recorded in the green channel to measure glutamate release at individual boutons. After 5 minutes of recovery, synaptic ATP-red signal was recorded in the red channel on the very same boutons. To evoke mitochondrial ATP production in synapses, neurons were challenged with 200 stimuli at and ATP dynamics recorded for 10 mins. Dual color imaging experiments were conducted in an open laminar flow field-stimulation chamber (cat.QR/RC-47FSLP, Warner Instruments, Hamden, CT) at 23 to . The extracellular buffer (EB1) contained 125 mM NaCl, 2.5 mM KCl, 2 mM MgCl2, 2 mM CaCl2, 30 mM glucose, and 25 mM Hepes (pH 7.4)108,109. To block recurrent activity, EB1 was supplemented with 6-Cyano-7-nitroquinoxaline-2,3- dione (CNQX) ; cat.479347-85-8, Tocris) and DL-AP5 ; cat.79055-68-8, Tocris). Electrical-evoked responses were induced by delivering 1 msec pulses with an amplitude using a train generation unit (Digitimer Ltd, DG2A) connected to a stimulus isolation unit (SIU-102; Warner Instruments, Hamden, CT). Recording chambers were placed on the stage of an IX-71 inverted microscope (Olympus, Hamburg, Germany) equipped with an EMCCD (electron-multiplying CCD) camera (QuantEM-SC 512x512, Photometrics). Illumination was obtained using a light-emitting diode LED (Cairn research, Optoled Lite), coupled to the acquisition software (MetaFluor, Molecular Devices) to minimize photobleaching. A 40X (1.3 numerical aperture; N.A.) oil-immersion objective was used to identify an axonal branch of an iGluSNFR positive neuron. Evoked iGluSNfr or ATP-red responses were recorded in an pixels) region of interest (ROI) typically containing 20 to 100 individual boutons. Time-lapse recording of glutamate dynamics was performed with an acquisition rate of for seconds, while recording of mitochondrial ATP was performed with an acquisition rate of for 12 minutes. Images were analyzed using ImageJ (NIH) and MATLAB (MathWorks) custom-made software scripts.109 Glutamate synaptic responses were identified by subtracting the resting iGluSNFR fluorescence from the peak fluorescence after a burst of . This was followed by the generation of ROIs that were used to measure the average fluorescence response in all identified boutons in each experiment. After subtracting the background, the data were normalized to the resting iGluSNFR or ATP-red signal ( F/F0).
The bioenergetic properties of primary hippocampal neurons and cultured microglia obtained from Trem2+/+ and newborn pups were determined using the XF-96 Seahorse extracellular flux analyzer (Seahorse Bioscience, CA). Neurons and microglia were grown at , with , up to 4 DIV (neurons) or for 24 hours (microglia). At indicated time points, neurons and microglia were washed and assayed in the XF Assay medium, XF DMEM pH 7.4 supplemented with glucose, 1 mM pyruvate sodium and L-Glutamine (cat.103575,103578,103577 and 103579, Agilent Technologies) for Mito Stress Test respiration measurements according to the manufacturer’s protocol. Measurements were obtained in real-time with no drug treatment (basal conditions) and with the sequential injection of Oligomycin, FCCP, Rotenone plus Antimycin A (Rot/AA) (cat.103015- m m m100, Agilent Technologies). Oxygen Consumption Rate (OCR) data were obtained as picomoles of consumption in the media and normalized against the control condition. OCR measurements under basal conditions in the absence of drugs represent the basal OCR. OCR measurements after FCCP injection represent the maximal OCR. OCR measurements after Oligomycin injection represent the ATP-coupled respiration. OCR measurements after Rot/AA injection represent the non-mitochondrial residual respiration. Three cycles of measurements were acquired for each condition. At least 15 replicates for each condition, for each experiment, have been assayed. OCR values, under basal condition and after the drugs’ additions were then plotted and run for statistical analysis using GraphPad Prism.
P1 hippocampi were carefully dissected in cold HBSS under a stereotactic microscope and mechanically digested in ice with cold buffers at as previously.110
Cells were first incubated with FC blocking CD16.32 (cat.101301, clone 93, 1:100, Biolegend) for 15’, then incubated in FACS buffer FBS, 1 mM EDTA in phosphate buffer saline) with the following fluorescently conjugated antibodies: CD45-Pe or PerCP (cat.103106 and 103130, clone 30-F11; 1:500, Biolegend), CD11b-Pe-Cy7 (cat.101216, clone M1-70; 1:200, Biolegend) for in ice. After surface staining, samples were washed twice with FACS buffer and centrifuged at for at . Viability dye staining was then performed for 15’ by resuspending the samples in phosphate buffer saline (PBS, cat.TL1006, Microgem) containing Zombie Nir (cat.423105, 1:800, Biolegend). For intracellular staining FOXP3 buffer kit (cat.00-5523-00, Thermo-Fisher) was used. Permeabilization and fixation, as well as subsequent washing steps incubations and centrifugations, were performed according to the protocol. Intracellular staining was performed for 1h on ice and the following antibodies were used: anti-NeuN 488 (cat.MAB377X, clone A60; 1:300 Merck-Millipore), anti-Puromycin AF647 (cat.AB_2827926, clone R4743L-E8; 1:400, Auguello Lab). For metabolic staining, the cells were incubated with CellROXgreen (cat.C1044, Thermo-Fisher) or with Mitotracker Deep Red FM (cat.M22426, ThermoFisher) for 35’ at . Metabolic staining was not ‘‘coupled’’ with intracellular staining. FACS analysis was performed on a FACS Fortessa machine or FACS Canto II (BD Biosciences), and data were analyzed with FlowJo Software (TreeStar).
SCENITH was performed according to Arguello et al., € 49 with minor modifications (www.scenith.com). After dissection, hippocampi were digested following a papain-based protocol (cat.LK003153, Worthington Biomedical) for minutes at under continuous agitation. After neutralizing the reaction with Neurobasal-A medium supplemented with FBS, digested hippocampi were gently homogenized with a P1000 pipette, filtered with a strainer, and centrifuged at for at RT. The obtained cell suspension mwas then resuspended and plated at a density of 150.000 cell/well in a 96 multi-well plate. Cells were allowed for an additional 45’ for recovery at and , after which they were incubated with a mix of different metabolic inhibitors for at and . The following metabolic inhibitors were used: 2-Deoxy Glucose, Oligomycin, 2-Deoxy Glucose combined with Olimgomycin (same concentrations as stated above, Scenith). Next, m Puromycin was applied to the samples, and left in incubation for at and m. Cells were then washed twice in FACS buffer and centrifuged at for and processed for FACS as described in the previous paragraph.
The bioenergetic properties of ex vivo acute brain slices obtained from and male offspring aged P1 or P18, were determined using the XF-96 Seahorse extracellular flux analyzer (Seahorse Bioscience, CA) with adapted protocol.111,112 Mice were deeply anesthetized by inhalation and then decapitated. Brains were rapidly removed and immediately immersed in ice-cold dissection aCSF solution previously saturated for at least 15 minutes with and . aCSF solution contained (in mM): 129 NaCl, , 10 glucose, , , , 2 . Coronal hippocampus slices were cut using a VT1000S vibratome (Leica Biosystems Nussloch GmbH) and maintained at m in a high-sucrose protective solution at for a minimum of 60 minutes before being used for the forthcoming processing.
For tissue subprocessing, 1mm-diameter -thick punches covering hippocampal CA area on P1, and both CA1 and CA3 subfield areas for P18 brain slices were manually obtained with a 1-mm Rapid-Core sampling puncher (Ted Pella, Inc). Punches were rapidly immersed in oxygenated Assay medium (in mM, , 5 KCl, , , , 25 HEPES, 0,01 Pyruvate and 6 Glucose, pH 7.4) and plated into a XF-96 cell culture plate, single punch/well, previously coated with poly-D-Lysine . Punches were washed and assayed in the XF Assay medium for Mito Stress Test respiration measurements. Measurements were obtained in real-time with no drug treatment (basal conditions) and with the sequential injection of Oligomycin, FCCP, Rotenone plus Antimycin A (Rot/AA) (cat.103015-100, Agilent Technologies). Oxygen Consumption Rate (OCR) data were obtained as picomoles of consumption in the media and normalized against the control condition within each experiment. Given the slower kinetics in sliced tissues, 8 cycles of measurement were used for most conditions. OCR measurements under basal conditions in the absence of drugs represent the basal OCR. OCR measurements after FCCP injection represent the maximal OCR. OCR measurements after Oligomycin injection represent the ATP-coupled respiration. OCR measurements after Rot/AA injection represent the non-mitochondrial residual consumption. OCR values, under basal condition and after the drugs’ additions, were normalized within each experiment over the counterpart, then plotted and run for statistical analysis using GraphPad Prism. Samples that did not display basal respiration values above the advisable ranges (>10 pmol/min) or did not respond to Seahorse drugs were discarded from the analysis.
P1 hippocampi from and pups from 3 independent litters embryos per sample per litter) were rapidly dissected in ice cold HBSS solution and pooled together. A balanced pool of males and females in both control and Trem2-/- mice were employed, to minimize the sex effect.
Dissected tissue was dissociated using Papain Dissociation System kit (Worthington, cat#LK003150) following the manufactured instructions as previously optimized. Dissociated hippocampi were suspended in PBS supplemented with BSA (Sigma Aldrich). Around 8,000 cells per sample were loaded into one channel of the Single Cell Chip using the Single Cell single cell reagent kit (10X Genomics) for Gel Bead Emulsion generation into the Chromium system. Following capture and lysis, cDNA was synthesized and amplified for 14 cycles following the manufacturer’s protocol (10X Genomics). 50 ng of the amplified cDNA were then used for each sample to construct Illumina sequencing libraries. Sequencing was performed on the NextSeq2000 Illumina sequencing platform following Genomics instruction for reads generation.
Raw sequencing data were converted to fastq files using Illumina bcl-convert tool, integrated into the CellRanger (10X Genomics) suite (version 6.1.1.113 The CellRanger analysis pipeline was used to generate a digital gene expression matrix starting from raw data. Pre-build mouse genome (version mm10-1.2.0) was used as genome reference. CellRanger count module was used to map reads with default settings setting and sequence length set to r1-length and– -length . At least 35,000 reads per cell were produced for each sample. The raw digital gene expression matrix (UMI counts per gene per cell) was imported in R version 4.4.2 using Seurat R package Briefly, UMI counts per gene per cell for each biological replicate was imported in Seurat and sample quality control was assessed by filtering out cells meeting any of the following criteria: less than 1,000 unique genes expressed, less than 1,000 UMIs, or more than of reads mapping to mitochondria. For each time point, data were normalized for cell feature expression by multiplying its total expression by a scale factor of 10,000, and -transformation of the result. This was followed by assignment of cell cycle scores to individual cells based on the expression of G2/M and S phase markers. We next scaled the expression values and identified the 2,000 most variable genes with FindVariableFeatures Seurat function (selection.method ’’vst’’, nfeatures ). We then performed principal component analysis (PCA) linear dimensionality reduction on the scaled data and clustered the cells with a graph-based clustering approach (RunPCA). We perform SCT-Normalization (negative binomial regression) without regressing out any variable (SCTransform(seurat, vars.to.regress NULL, return.only.var.genes ) for each individual sample. Due to differences in cell clustering due to the sex, we integrated the data using Seurat built-in SelectIntegrationFeatures(nFeatures $= \ 3 { , } 0 0 0 $ ), followed by PrepSCTIntegration, FindIntegrationAnchors(normalization.method “SCT”) and IntegrateData functions. To choose the number of principal components we used a graph-based method (elbow plot). We then clustered the cells using the Louvain algorithm (within Seurat) to iteratively group cells together. Optimal resolution was chosen by analyzing the cluster tree plot (clustree). To annotate clusters, we determined differentially expressed genes using FindAllMarkers function embedded in Seurat (Wilcoxon Rank Sum test with Bonferroni correction for multiple testing; adjusted ). We tested only genes that were detected in a minimum of of the cells within the cluster and that showed, on average, at least a 0.25-fold difference (log-scale) between the cells in the cluster and all remaining ones. By reviewing the resulting markers, as well as the expression of canonical marker genes and additional markers from the literature,115–117 we assigned a cell-type identity. We then filtered out vascular, mesenchymal and ependymal cells and re-clustered the remaining cells as previously described.
To test the differential cells abundances between genotypes, we used MiloR 1.2.0 function on the KNN-graph generated from the latent space embeddings with default parameters, where ‘genotype’ was used as a covariate.
To determine differentially expressed genes, we used FindMarkers from Seurat function implemented with default two-part hurdle model (MAST) with added random effect variables (MAST-re test; adjusted testing only genes found in at least of cells within the cluster and that showed, on average, at least a 0.10-fold difference (log-scale) across genotypes and setting the batch as random effect variable. The list of differentially expressed genes of CA-ImmPyr, CA1-Pyr and CA3-Pyr cell types was analyzed using the Ingenuity Pathway Analysis (IPA - Qiagen Ingenuity Systems) canonical pathways enrichment.
Pseudotime and trajectory tree were inferred using the Slingshot package in a two-step process. First, we identified the global lineage structure with a cluster-based minimum spanning tree (MST). Then, we fit simultaneous principal curves to describe each lineage onto the UMAP space using the slingshot function. Trajectory inference and visualization were performed on integrated pre-processed data object UMAP space and with previously assigned cell type labels. Apical Progenitor cells (AP) were chosen as the root of the trajectory. All the inferred pseudotime variables (one per lineage) are added to the metadata, individually. We next compared pseudotime distribution of cells in the three lineages by randomly subsetting 100 cells per each lineage and performing a two-sided Kolmogorov-Smirnov test of the empirical cumulative distribution functions. Finally, we repeated the same procedure to compare pseudotime distribution of cells for each genotype for individual cell types for each lineage.
Either frozen hippocampi from P1, P20, P90 or CA1/CA3 1-mm punches from Trem2+/+ P18 mice were gently thawed on ice and then subjected to mechanical disruption with homogenizing beads onto automatic Tissue Lyser, in 5 0 0 \mu \up TRI Reagent ^ \mathrm { \textregistered } (Zymo Research). mRNA was isolated with RNA Direct-Zol MiniPrep Isolation Kit (Zymo Research), according to the manufacturer guidelines. RNA was eventually eluted in 2 5 \mu \up DNAse/RNAse-free water and quantified with NANOdrop 2000c spectrophotometer (Thermo Fisher Scienmtific) for RNA concentration and optical density ratios.
500 ng RNA for each condition underwent reverse transcription into cDNA with High-Capacity cDNA RT kit (Applied Biosystems). Quantitative Real-time polymerase chain reaction (qRT-PCR) was performed with TaqMan detection kit (TaqMan Fast Universal PCR Master , no AmpErase UNG, ThermoFisher) onto qRT-PCR Viia7 software system (Applied Biosystems) in a final volume of 1 0 ~ \mu \mu \up . Each gene was subjected to at least duplicate measurements and data analyses were performed with the comparative method. mRNA measurements for each target gene were normalized to the housekeeping gene Gapdh. The following TaqMan assays were used (Applied Biosystems): mouse Trem2 FAM-MGB Mm04209424_g1; mouse GAPD(GAPDH) VIC-MGB Endogenous Control.
Whole cell voltage-clamp recordings were performed on 4DIV primary hippocampal neurons from and P1 mice. During recordings, neurons were mantained in extracellular buffer 2 in an recording chamber (cat.QR-40LP, Warner Instruments, Hamden, CT). The extracellular buffer (EB2) contained (in mM): 125 NaCl, 5 KCl, 1.2 , , , 6 glucose, and 25 HEPES-NaOH, pH 7.4. Recording pipettes (resistances of 3-5 M ) were filled with a standard intracellular solution containing (in mM): U-gluconate, 1 EGTA, 10 HEPES, 2 MgCl2,4 MgATP, and 0.3 Tris-GTP, . For voltage-dependent current activation experiments, neurons were held at and incremental depolarizing steps ( for , starting from -120 and heading to , were applied. current density was obtained by dividing peak inward current by cell capacitance (nA/pF). Resting potential was calculated at in current clamp configuration, whereas input resistance was calculated in voltage clamp configuration by using relationship slope of the steady state current measured at different depolarizing voltage steps (from to ). Recordings were performed in voltage clamp mode using a Multiclamp 700B amplifier and pClamp-10 software (Axon Instruments, Foster City, CA). Series resistance ranged from 10 to and was monitored all over the whole recordings. Signals were amplified, sampled at , filtered to 2 or U, and analyzed using the pClamp 10 data acquisition and analysis program. Cells in culture with leak currents pA were excluded from the analysis.
P18 Trem2+/+ and male mice were deeply anesthetized by inhalation and then decapitated. Brains were dissected out and immediately immersed in ice-cold dissection aCSF solution previously saturated for at least 15 minutes with O2 and CO2. aCSF solution contained (in mM): , , 10 glucose, , , 3 KCl, 21 . Coronal hippocampus slices were cut using a VT1000S vibratome (Leica Biosystems Nussloch GmbH) and maintained at in a high-sucrose protective solution at for a minimum of 60 minutes before being used for recordings. During the recordings, slices were placed on the multielectrode array and perfused with aCSF at 2 ml/min speed. After 15 minutes of slice spiking activity stabilization, spontaneous extracellular signals were collected for additional 5-10 minutes.
All extracellular multisite recordings were performed using a CMOS biosensor and acquisition system (3Brain AG, Switzerland). High-Density Multielectrode CMOS Arena chips integrating 4096 recording electrodes in a array configuration were employed. The acquisition was controlled by software BrainWave 4 (3Brain AG, Switzerland). Raw data were digitized at 18kHz.
For each experiment an image of the slice in the recording unit was acquired. Only the electrodes located in the CA1 and CA3 regions were analyzed. Custom-written MATLAB scripts were used for all MEAs analysis algorithms. At first, noise signals were identified and excluded from the analysis. Recordings were then high-pass filtered ) using a second-order Butterworth filter with BrainWave Software (3Brain Gmbh, Switzerland). A spike detection algorithm (Precise Timing Spike Detection, PTSD) was implemented to extract spike times from the 18 kHz signal. A threshold of 7 standard deviations (SDs) was set to filter out basal noise.
Spikes that occurred within the refractory period (set to 1 ms) were removed. Timestamps of every detected neuron were used to create raster plots. The average firing rate per neuron was calculated by dividing the total number of spikes by the recording session duration. To assess the interaction between the hippocampal subregions, a suitable cross-correlation measure was calculated using MATLAB xcorr function. To determine the degree of synchronization between different neuronal firing patterns, the maximum correlation at zero lag was calculated and a correlation matrix was built for each experiment and the mean correlation coefficient estimated for each experiment. The number of neurons analyzed, instead, corresponded to the number of active electrodes sampled for each hippocampal subregion.
P1 pups were sacrificed by decapitation and the brains quickly removed from the skull, washed in fresh PBS1X and fixed in paraformaldehyde (PFA) overnight. P18 were deeply anesthetized with xylazine/ketamine mixture and transcardially perfused with PFA and additionally post-fixed overnight. Coronal slices thick) were cut with a VT1000S vibratome (Leica Microsystems). Slices were washed in cold PBS1X, and permeabilized with Digitonin in PBS1X for 15 mins or TritonX-100 in PBS1X for 1 hour and nonspecific binding sites were blocked by incubation for with either normal goat serum (NGS), Bovine Serum Albumin (BSA) or donkey serum (DS) and TritonX-100 in PBS1X before an overnight incubation at 4C with primary antibodies.
Slices were then washed 3 times with PBS and incubated with secondary antibodies conjugated to Alexa Fluor (488, 555, or 633). The following primary antibodies were used: guinea pig anti-IBA1 (cat.234 308, 1:1000, Synaptic Systems), rabbit anti-IBA1 (cat.019- 19741, 1:1000, WAKO Chemicals), rabbit anti-P2yr12 (cat.AS-55043A, 1:400, Anaspec), mouse anti-Tom20 (cat.sc-17764, 1:200, Santa Cruz), guinea pig anti-DCX (cat.AB2253, 1:500, Merck-Millipore), guinea pig anti-MAP2 (cat.188 004, 1:1000, Synaptic System), rat anti-CD68 (cat.137001, 1:1000, Biolegend), rabbit anti-PSD95 (cat.51-6900, 1:200, Thermo-Fisher), sheep anti-Trem2 (cat.AF1729, 1:200, R&D Systems), rabbit anti-Neurod2 (cat.ab104430, 1:1000, Abcam), rat anti-Ctip2 (cat.ab18465, 1:100, Abcam). All slices were counterstained with Hoechst-33342 (cat.62249, 1:1000, Thermo-Fisher) and mounted with Fluorsave (cat.345759, Millipore). Images were acquired using a Leica SP8I or a SP8II laser scanning confocal microscope equipped with an HC PL APO CS2 objective and ACS APO or oil immersion objective. For Tom20 analysis, the entire pyramidal layer (P1) or a region of CA1 and CA3 hippocampal subfield (P18) were acquired using or oil immersion objective. Tom20 fluorescence intensity was estimated by generating a binary mask created on Dcx or Map2 staining of CA1 and CA3 pyramidal layer and normalized by area fraction. For staining related to Figure 5, a subset of control animals from Figures 3 and 4 was used. Neuronal density was estimated by counting NeuroD2-positive cells in the CA1 and CA3 pyramidal layer areas using Positive Cell Detection function implemented in QuPath software.119 Similarly, microglia density was estimated by counting Iba1 positive cells and normalized per CA1 and CA3 area. Hippocampal microglia morphology was assessed using Skeleton plugin in Fiji software (NIH, Bethesda), which enables a three-dimensional (3D) morphological analysis of skeletonized microglia present in tissue. P2y12r intensity at single hippocampal microglia cells was measured by generating a binary mask created on 3D Iba1 labeling. Microglia phagocytosis and Trem2 immunoreactivity were assessed by performing three-dimensional surface rendering and volumetric quantification using Imaris Software (Oxford Instruments, v. 9.7.2,).120 PSD95 intensity was measured as volume fraction of the whole image. For quantification of PSD95 engulfment, only the signal co-localized with CD68 within microglia was considered.
4DIV or 14DIV neurons were fixed with PFA in PBS1X and sucrose solution for and permeabilized with Saponin in PBS1X for 10 min. Nonspecific binding sites were blocked by incubation for 30 min with NGS and Saponin in PBS1X before a 2-h incubation at room temperature with primary antibodies. Neurons were washed 3 times with PBS1X and incubated with secondary antibodies conjugated to Alexa Fluor (488, 568, or 633) for at room temperature. Coverslips were washed 3 times in PBS1X, counterstained with Hoechst to reveal nuclei and mounted onto slides with Fluorsave. The following primary antibodies were used: mouse anti-Tom20 (cat.sc-17764; 1:200, Santa Cruz); guinea pig anti-Map2 (cat.188004; 1:1000, Synaptic Systems) and rabbit anti-VGlut1 (cat.135302; 1:1000, Synaptic Systems). Somatic and neuritic or synaptic fluorescence intensity was quantified using ImageJ (NIH). Somatic and neuritic Tom20 fluorescence intensity, corrected for background signal, was estimated by manual tracking of cell bodies and primary neurites while Tom20 synaptic intensity was quantified by generating a binary mask created from the vGlut1 staining over-imposing GFP positive neurons.
Sholl analysis of 4DIV sparsely distributed neurons was performed using the Sholl analysis Fiji plugin. Briefly, circles of increasing radius, originating from the center of the transfected cell body, were drawn, and the numbers of intersections between dendritic branches and the edges of the circles were used as an estimate of neuronal arborization. Images were acquired using a Leica SP8I confocal microscope equipped with an ACS APO or 63x oil immersion objective.
and male mice at deeply anesthetized with xylazine/ketamine mixture and transcardially perfused with fresh PBS1X. Brains were removed and placed in ice-cold solution containing the following (in mM): 129 NaCl, 1.25 NaH , 10 glucose,
, , 3 KCl, 21 NaH (pH 7.4), equilibrated with and . Coronal slices thick) were cut with a VT1000S vibratome (Leica Microsystems) from medial Prefrontal Cortex (PFC). Slices were washed in cold PBS1X and then fixed for 3 hours with glutaraldehyde in 0.1M Cacodylate buffer at room temperature. Slices were post-fixed in osmium tetroxide (cat.19140, Electron Microscopy Science, Hatfield, PA, USA) for 2 hours and in aqueous uranyl acetate (cat.22400-1, Electron Microscopy Science) for 1 hour at room temperature.
Slices were then dehydrated through a graded ethanol series with propylene oxide as a transition fluid (cat.P021, TAAB Laboratories Equipment, Aldermaston, UK) and embedded in epoxy resin (Poly-Bed; cat.08792-1, Polysciences, Warrington, PA, USA) overnight at e ultimately for 2 days at . Subsequent examination of the regions of interest (hippocampal CA1 and CA3 regions) was performed on semithin sections stained with Toluidine Blue. Ultrathin sections were then cut and counterstained with uranyl acetate. Electron micrographs were acquired using a HT7800 transmission electron microscope equipped with Megaview III digital camera and Radius 2.0 software (EMSIS, Muenster, Germany). Tile scan images were acquired to quantify the major axis length of mitochondria in pyramidal neurons somata. Single electron micrographs were used to acquire synapses in the pyramidal layer of both CA1 and CA3 regions. Morphometric analysis was done using Fiji software (NIH, Bethesda). Postsynaptic electron density was used to identify excitatory synapses. Structures with sagittal diameter between 20 and were classified as SVs. SVs touching AZ were classified as docked SVs.
The distribution of data in each set of experiments was first tested for normality using the Shapiro–Wilk test. The similarity of variances between each group of data was tested using the F test. Normally distributed data are presented as mean SEM; each plot also contains the individual data points. Unpaired Student’s t tests for group means, One-Way ANOVA or Two-Way ANOVA followed by post-Hoc multiple comparisons test was used as indicated. To compare datasets that failed the normality test, Mann–Whitney U test and ANOVA on Ranks followed by Dunn’s multiple comparison test were used as indicated. No statistical methods were used to predetermine sample sizes, but our sample sizes are like those reported in previous publications in the field. All statistical tests were performed using GraphPad9 (Prism software).