Microglia regroup into morphologically distinct cell aggregates in AD hippocampi.
To elucidate the alterations of microglia within the hippocampus of AD patients, we performed IHC analysis of the microglial morphological marker, ionized calcium-binding adaptor molecule 1 (Iba1),in a cohort of CTLs and AD post-mortem FFPE samples. We identified two distinct types of large Iba1-positive cellular accumulations. The most common were the typical Aβ plaque-associated microglia (PaM), characterized by a high number of microglia clustered around an Aβ core plaque in a rosette-like conformation. Less frequently encountered, but consistently present in our sample collection, were clusters of microglia with a polarised shape reminiscent of a 'soma' and 'tail' configuration, observed predominantly in the pyramidal layer (PL) of the CA1 and CA2 regions of the hippocampus (Fig. 1A-B). Across serial sections of hippocampi from 16 AD cases and 5 age-matched CTLs, we observed a ratio of approximately 1 of this type of accumulation to 11 PaM in the CA of AD patients, whereas it was absent in CTL samples, as well as in the prefrontal, temporal and entorhinal cortices of both CTL and AD patients.
To further characterise the three-dimensional (3D) morphological features and differences from PaM, we immunostained PFA-fixed 80 µm hippocampal sections from AD patients for Iba1 and 4G8 (Aβ). Using confocal and STED microscopy, we confirmed the enrichment of this peculiar type of Iba1- accumulation in the CA1 and CA2 PL of AD hippocampi (Fig. 1C) and showed that they were spatially separated from 4G8-positive Aβ plaques (Fig. 1D-G, Supplementary Fig. 1A-B). High and super-resolution microscopy revealed new features of their complexity with Iba1-cells showing tight boundaries and an encapsulating cellular organisation that was clearly distinct from PaM [9, 16] with a size and consistent orientation resembling that of pyramidal neurons (Supplementary Fig. 1, Supplementary Movie 1; Supplementary Movie 2). Co-staining with DRAQ7™ often showed nuclei that appeared to be engulfed by Iba1-positive cells (Supplementary Movie 1, Supplementary Fig. 2). These microglia collectively formed a tight shell around as yet-unidentified structures, resembling the glial embrace and coffin formation described in early twentieth-century literature [2, 24, 38], that we hereafter referred to as “coffin-like microglia” or CoM (Fig. 1A, C, E, G; Supplementary Fig. 1B, D).
Our high-resolution 3D analysis allowed us to distinguish between two types of CoM. We noticed some CoM formations were predominantly composed of bipolar or rod-shaped microglia, resulting in a thin and elongated structure in 3D, often traversing the PL, while only partially enveloping surrounding cells. Unlike the first CoM described, the other type did not seem to enclose any other structure. We termed them CoM-Rod (Fig. 1H, Supplementary Fig. 3). In addition, we often observed smaller groups of Iba1-positive cells 'clumped' together, located between other brain cells (Supplementary Fig. 3). Our previous findings indicated that morphological alterations of individual microglia in DLB resembled to those observed in AD, but with less severity [19]. Consistent with these observations, we detected PaM, CoM, CoM-Rod, and clumped microglia to a lesser extent in the hippocampi of DLB patients (Supplementary Fig. 1C-D).
We quantified CA1/CA2 PaM and CoM number of cells and their total volume from 3D confocal stacks (40x). We did not observe CoM in age-matched control samples (n = 8), however, we found few PaM and CoM-Rod. Using Imaris software, we quantified volumes and number of nuclei contained in 61 3D z-stack images of PaM, CoM and CoM-Rod across age-matched control (n = 1), AD (n = 8) and DLB patients (n = 4) samples (Fig. 1H). PaM presented a mean volume of 19,744µm3 with an average of 20 nuclei, CoM had a mean volume of 15,776µm3 and an average of 17 nuclei, and CoM-Rod a mean volume of 8,870µm3 and an average of 9 nuclei (Fig. 1H). Using the Kruskal-Wallis test by ranks, we showed that CoM-Rod are significantly smaller in volume (P < 0.05, P < 0.01 respectively) and in nuclei count (P < 0.01) than CoM and PaM. Even though PaM and CoM are very different in morphology, their volume and nuclei count did not present significant differences.
We then investigated the association of a typical marker of microglia upregulated in disease, cluster of differentiation 68 (CD68) [7, 27, 68]. We first quantified CD68 expression in hippocampal subfields of CTL and AD FFPE samples (n = 7 CTL, 8 AD). CD68 was expressed in microglia-like cells in CTL and AD in all brain regions examined, slightly increased in all subfields but more drastically in the CA1 AD (P < 0.05). We found CD68 associated with PaM and CoM-like structures in the CA hippocampus, showing a shared feature of PaM and CoM (Fig. 1I).
PaM and CoM exhibit distinct pathological and astrocytic microenvironments.
To elucidate the relationship of PaM and CoM with surrounding cells, tau and α-synuclein pathologies, we investigated the characteristics of their pathological microenvironment. Given the distinct association of PaM and CoM with Aβ plaques, we further investigated their interactions with pTau NFTs and pSyn aggregates in AD [19, 92]. Using high-resolution confocal imaging of immunostained sections for Iba1 and AT8 (pTau), we found that CoM can fully encompass NFT-bearing neurons, enveloping the entire 'flame' shape of the tangle and spreading to surrounding neuronal cells (Fig. 2A-B, Supplementary Movie 3). Using thiazine red (TR), which binds β-pleated sheet structures, we also observed CoM containing TR-positive structures resembling NFTs (Supplementary Fig. 4B, Supplementary Movie 4). Similarly, we found some cells harboring intracellular granular pSyn inclusions enveloped by CoM (Fig. 2C-D, Supplementary Movie 5). Notably, some CoM lacked AT8 tangles (Supplementary Fig. 4A) or pSyn aggregates.
We next examined the distribution of the P-Ser396 pTau (PS396), a preferred site of GSK-3 β phosphorylation activated by phospho-AKT1 (S473) [78], in combination with the S202 and S205 phosphorylated forms of tau labelled by AT8, 4G8 and Iba1 staining in IHC multiplex chromogenic experiments (n = 4 AD). We found relative heterogeneity of PS396 and AT8 tangles and neurites across samples, with some samples showing a distinct ratio of PS396 to AT8 loading (Fig. 2E-F). We observed single stained neurons, AT8 + or PS396+, and double stained AT8 + PS396+. The plaque microenvironment defined by 4G8 staining and Iba1 PaM showed a severe tau pathology, with both in AT8 and PS396 labelled swollen structures and often double-positive (Fig. 2E-F). Groups of microglia reminiscent of CoM were also found surrounding neurons with double positive tangles (Fig. 2F). To understand the compartmentalisation of neurodegenerative processes around Aβ plaques, we used multiplex chromogenic labelling for PS396 Tau, 4G8 and glial fibrillary acidic protein (GFAP), to label PaM astrocytes that were shown to form nets [9]. We found that astrocytes frequently encapsulate PS396 + structures and Aβ plaques (Fig. 2G).
We then investigated in more detail the relationship of PaM and CoM with surrounding astrocytes, whose morphological changes may indicate a specific interplay with microglia. Using chromogenic multiplex IHC, for Iba1, GFAP and 4G8 (n = 5 CTL, 8 AD), we quantified the number of PaM and PaM-associated astrocytes (PaM astrocytes), those polarised around PaM and forming a net, associated with Aβ plaques (n = 391) in the CA1 and measured the corresponding plaque area (Fig. 3A-C). As expected, more PaM than PaM astrocytes were attached to Aβ plaques per area (Fig. 3B). We found a positive correlation between CA1 plaque area, the number of PaM (r = 0.398, ***) and the number of astrocytes forming a net (r = 0.37, ***) (Fig. 3C) with a similar trend as in the cortex [9]. In total, 26.1% of Aβ plaques were encapsulated by both PaM and PaM astrocytes, 22% of Aβ plaques smaller than 2500 µm2, 53.3% of Aβ plaques larger than 2500 µm2. With this approach, there was no evidence of CoM encapsulation by astrocytes. Using 3D confocal microscopy, we compared the distribution and morphology of the astrocytes surrounding PaM or CoM. Around core Aβ plaques (4G8), hippocampal PaM (Iba1) were often surrounded by an outer sphere of hypertrophic and polarised astrocytes (GFAP, PaM astrocytes), as previously observed in AD cortical brain samples (Fig. 3D), but also in APP mouse models [9]. However, we found that astrocytes surrounding CoM ( CoM astrocytes) mainly exhibited a dysmorphic phenotype, with thin and irregular branches, lacking polarisation towards the CoM and showing lower reactivity, here assessed by GFAP intensity (Fig. 3E, Supplementary Fig. 4B). Some GFAP processes were found invading CoM boundaries suggesting a possible interaction at subcellular level different from that around Aβ plaques (Supplementary Movie 4).
Digital spatial profiling of 147 proteins highlights core differences between PaM and CoM and their microenvironment.
To delve deeper into the molecular signatures of PaM and CoM in relation to their surrounding microenvironment in AD, we used the Nanostring GeoMx Digital Spatial Profiler (DSP) platform. We quantified 147 proteins in spatially defined CoM (n = 5) and PaM (n = 7) Iba1-positive areas and their close surrounding in FFPE hippocampal sections from three AD patients (cases #49, #50 and #52; Table 1). Utilising fluorescent antibodies for microglia (Iba1), astrocytes (GFAP), Aβ-plaques (MOAB-2) and nucleic acid stain (SYTO™ 13), we accurately localised and identified PaM with astrocyte net and CoM. By manually delineating regions of interest for PaM and CoM in the hippocampal PL from the CA2 and CA1, we selected two AOIs: the Iba1-positive area and its immediate surrounding (Fig. 4A). After negative probes normalisation, only a limited number of DEPs reached significance FDR > 0.1, |Log2 fold change (FC)|>0.5, padj<0.05) between PaM and CoM and between their respective surroundings (Fig. 4B, Supplementary Excel 1). Her2 (receptor tyrosine-protein kinase erbB2), complement component 4B (C4B), and protein kinase B (Phospho-AKT1, S473) were enriched in PaM compared to CoM while stimulator of interferon genes (STING) exhibited higher levels in CoM. Granzyme A (GZMA), progesterone receptor (PR), and p53 were enriched in vicinity of PaM vs. CoM. Despite not reaching statistical significance, several proteins indicated potential differences in cellular signatures or microenvironments. Purinergic receptor P2RY12, Iba1, and microglia/immune cells associated proteins such as arginase 1 (ARG1) or lymphocyte activating 3 (LAG3) were elevated in CoM compared to PaM. Amyloid peptides, amyloid precursor protein (APP), GFAP, CD163, CD44 and apolipoprotein E (ApoE) were found at higher levels in both PaM and PaM surroundings in coherence with our data and literature [58, 62]. We further validated some of these findings with immunohistochemistry (IHC), in FFPE and/or in PFA AD hippocampal sections. We confirmed the enrichment of STING in CoM, as well as C4d, a cofactor of C4B, and ErbB4, a co-receptor of ErbB2 [43] in PaM and PaM surroundings (Fig. 4C-E). Specifically, STING expression was elevated in AD hippocampal samples (8 AD, 7 CTL, FFPE), often associated with blood vessels and microglia-like cells in the white matter of the stratum oriens, CA3, and particularly the PL of CA1 (Fig. 4C). Often microglia aggregated around pyramidal neurons of the CA1 were STING immunoreactive. C4d staining was enriched around Aβ plaques in the hippocampus (CA3, CA1 and DG), subiculum, entorhinal and temporal cortex in a set of our samples (3 AD, 3 CTL, FFPE), validating previous observations [85] (Fig. 4D). Our 3D confocal analysis revealed the presence of ErbB4-positive PaM encapsulated by GFAP-positive astrocytes in PFA-fixed hippocampal sections (Fig. 4E). These observations highlight distinct functional states and pathological microenvironments between PaM and CoM.
Digital spatial transcriptomic profiling reveals distinct and complex molecular and metabolic states of PaM vs. CoM and surrounding astrocytes.
To unravel the molecular profiles and functional states of PaM and CoM and their respective immediate surrounding cells, we used the GeoMx DSP transcriptomic platform to quantify gene expression (18,677 RNA targets) in Iba1-positive PaM and CoM, and their respective surrounding GFAP-positive astrocytes (Fig. 5A). In addition, we also profiled some Rod, defined by their elongated morphology, and sometimes clustered in CoM-Rod (Supplementary Fig. 5A). We assessed a total of 38 AOIs across two consecutive sections from three AD patient samples (cases #49, #50 and #52): 10 PaM-Iba1, 10 PaM-GFAP, 7 CoM-Iba1, 7 CoM-GFAP, and 4 Rod-Iba1 AOIs. Under UV illumination, the indexing oligonucleotide tags were collected from the probe mix and counted individually from each Iba1-positive and, their respective, GFAP-positive AOI. Briefly, all AOIs passed sequencing quality control and targets with expression in less than 5% of the segments above the LOQ were removed. The target count matrix was scaled to the geometric mean of the nuclei and normalised by the upper quartile (Q3) normalisation. A final dataset of 15,518 genes was analysed.
Comparative analysis of the 500 most highly expressed genes in each microglial and astrocyte subgroup revealed overlaps and unique gene sets. We found an overlap of 197 genes (20.6% of total genes) between all the three groups of microglia, and 290 genes (40.8%) between the two groups of astrocytes (Fig. 5B). Approximately 20% of genes were exclusive to one specific group for microglia and 30% for astrocytes. Among the highly expressed genes in all microglial subgroups were allograft inflammatory factor 1 (AIF-1), CD74, or colony-stimulating factor-1 receptor (CSF-1R).
In contrast, astrocytic subgroups exhibited expression of astrocyte-specific markers such as GFAP, aldolase C fructose-bisphosphate (ALDOC), and ezrin (EZR). We then analysed the DEGs between the microglia or astrocyte subgroups using unpaired t-test and multiple comparison adjustment (Benjamini-Hochberg with FDR of 0.1). We identified 1,238 DEGs in our microglia dataset and 895 DEGs in our astrocyte dataset. Unsupervised clustering resulted in a clear separation of all AOIs into morphologically predefined groups for the microglia and astrocyte datasets, as shown in the heatmaps (Fig. 5C). DEGs were displayed in volcano plots for all comparisons (Fig. 5D, Supplementary Fig. 5B-C for Rod) and the top 30 DEGs identified in heatmaps (most enriched and lowest p-value) (Fig. 5E, G, Supplementary Fig. 5D-E). When comparing PaM to CoM transcripts, we found RAB27B, RHOQ, PCED1B, TMED8, HCAR1, CLVS1, GOLPH3, PTPRT, DCHS1, TNFAIP1, FOXO1, OXCT2, and HMCES among the most enriched and statistically significant genes for PaM; while in CoM they were C9orf131, TBXT, CASP2, AICDA, E2F6, ATP5F1D, MAP2K5, GIPC3, PPP1R9A, FBXL15, ACADS, ZMAT2, GSDME DUS3L, GDF1, and BLOC1S3. Rod-shaped microglia also showed numerous DEGs compared to both CoM and PAM. TAAR1, BCKDHA, CTNNAL1, EDIL3, RYR1, CALML4, RAX2, SLC20A2 and SERPINB6 were among the top 20 most DEGs in both comparisons (Supplementary Fig. 5B-C, Supplementary Excel 1).
We then performed an over-representation analysis (ORA) for PaM and CoM DEGs. The Gene Ontology (GO) pathways over-represented in PaM revealed phosphorylation and lipid degradation activities, antioxidant and glutathione transferase and tyrosine kinase activities. Consistent with this, the ORA based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) showed associations with sphingolipid signaling pathways, insulin resistance, ErbB and Wnt signaling pathways or necroptosis (Supplementary Excel 2) (Fig. 5H). The top five GO pathways for CoM were associated with ubiquitin-like protein transferase, purinergic nucleotide receptor, oxidoreductase cyclin-dependent kinase activities (Fig. 5H). ORA KEGG and Wikipathways (WP) suggested associations with transforming growth factor (TGF)-β and NF-κB signaling, among others (Supplementary Excel 2).
Next, we analysed DEGs in PaM astrocytes vs. CoM astrocytes. PaM astrocytes were associated with distinct signatures, and we found GOLGA8F, PCSK9, HOXD13, PIWIL1, PIK3C3, FAM71E2, FLG2, GBA2, SLC10A5, ACE, FBXO10, ARHGAP29, CES1, TMEM273, and GALNT7 among the most enriched genes (Fig. 5F-G, Supplementary Excel 1). We found SERPINA5, TUBB2A, FAM171A1, PROB1, GP1BA, TMEM37, FNDC4, JUNB, SCN5A, TSC22D4, GFPT1, BASP1, and UBE2W among the most enriched and statistically significant genes for CoM astrocytes. We performed ORA on PaM astrocytes and CoM astrocytes DEGs. ORA GO of DEGs PaM astrocytes showed enrichment for fibroblast growth factor receptor signalling, response to TGF-β, glucose transmembrane transport, positive regulation of apoptotic process and amyloid beta clearance (Supplementary Excel 2). ORA KEGG and WP indicated activities associated with PI3K-Akt signalling pathway, cholesterol metabolism or cell adhesion molecules, fatty acid and lipoprotein transport or complement system. ORA GO highlighted the association of DEGs CoM astrocytes with metabolite and energy generation (OXPHOS), cytokine mediated pathways but also with adaptive immune response, JAK-STAT signalling and cellular response to misfolded proteins (Fig. 5H, Supplementary Excel 2). Our DSP analysis revealed specific cellular and molecular identities of PaM and CoM, and of their respective surrounding astrocytes.
Markers of the molecular signatures of PaM, CoM and surrounding astrocytes in the hippocampus of AD patients.
To validate the molecular signatures identified in PaM and PaM astrocytes by our DSP analysis, we performed IHC on selected proteins that embodied an enriched pathway in FFPE or PFA samples. First, we examined the complement component 1q (C1q) (n = 5 AD, 5 CTL) and complement component 3 (C3) (n = 3 AD, 2 CTL) association with PaM and the PaM environment or PaM astrocytes. In AD, C1q was strongly distributed around plaques, and on some hippocampal neuronal cell bodies and neurites in the CA3-CA1 subfields, but also in the DG (Fig. 6A) and entorhinal cortex (Supplementary Fig. 6A), and occasionally diffuse in areas reminiscent of astrocyte domains. C1q-positive hippocampal neurons were also found in CTL (Supplementary Fig. 6A). Similarly, C3 staining exhibited lower overall expression, primarily labeling the core and corona of Aβ plaques in AD (Fig. 6B, Supplementary Fig. 6B). Using high-resolution 3D confocal imaging, we observed C3 expressed within GFAP-positive PaM astrocytes territories, in cell bodies but also along processes (Fig. 6C). We then analysed the distribution of Proprotein convertase subtilisin/kexin type 9 (PCSK9) (n = 3 AD, 2 CTL) and Integrin subunit alpha 6 (ITGA6, also known as CD49f [5]) (n = 5 AD, 5 CTL), associated with the PaM astrocytes and PI3K-Akt pathways in our previous analysis. In FFPE sections, we observed a PCSK9 expression in neurons both in CTL and AD ( Supplementary Fig. 6), however in AD samples, PCSK9 staining was also associated with plaque microenvironments and with cells resembling PaM astrocytes in hippocampus and entorhinal cortex (Fig. 6D-E, Supplementary Fig. 6). ITGA6 was predominantly distributed along blood vessels in CTL, while in AD, expression around blood vessels was weaker and occasionally observed near plaques, apparently associated with PaM astrocytes. We confirmed a low expression of ITGA6 in the cell bodies and processes of PaM astrocytes with confocal microscopy (Fig. 6F). We then investigated the association of the PaM microenvironment with necroptosis identified through pathway analysis. We found, in bright-field and 3D confocal images, a clear enrichment of phosphorylated mixed lineage kinase domain-like protein (pMLKL), part of the necrosome complex, in PaM vicinity but no staining in CoM (Fig. 6G). PMLKL was frequently concentrated in vacuole-like accumulations within the cytoplasm of numerous pyramidal neurons in AD samples and, to a lesser extent, in CTL as previously described [39].
DSP analysis suggested that CoM was involved in different biological processes than PaM. To further elucidate these processes, we performed co-immunostaining against Iba1, GFAP and NF-κB -p65 in PFA sections revealing strong expression and translocation of NF-κB -p65 into CoM nuclei by confocal microscopy (Fig. 6H). The staining was not exclusive to CoM, but also detected in the surrounding microenvironment, without being predominantly expressed in astrocytes. In contrast, PaM and Pam astrocytes displayed heterogeneous NF-κB -p65 staining patterns (Fig. 6I). ORA analysis revealed strong ubiquitin ligase activity in CoM. Confocal microscopy analysis of E3 ubiquitin-protein ligase SMURF2 distribution in AD PFA-fixed samples revealed expression in CoM, although not significantly enriched compared to the surrounding microenvironment (Fig. 6J). Acyl-CoA dehydrogenase (ACADS) is associated with fatty acid oxidation and high metabolic activity. In CTL samples, ACADS was expressed in astrocytic cells in the temporal cortex, CA4 hippocampus and CA3-CA1 PL. In AD, astrocytes retained ACADS expression, but so did some microglia surrounding pyramidal neurons. 3D confocal analysis confirmed the co-distribution of ACADS in the CoM (Fig. 6K-L).
Macrophage and T cell association with PaM and CoM.
To determine whether infiltrating immune cells contribute to the PaM or CoM molecular markers and composition, we next examined the distribution of T cells and CNS-associated macrophages in CTL and AD hippocampal samples. Indeed CoM exhibited a cellular organisation reminiscent of nodules observed in various brain infections, often containing T cells [77, 90]. First, we quantified the number of CD8 + T cells, which have previously been detected near Aβ plaques in AD brain samples, in CTL and AD hippocampi (4CTL, 6 AD). CD8 + T cells were mainly found in blood vessels (BV) and less frequently in hippocampal parenchyma (Fig. 7A). CD8 + T cells were slightly increased in AD without reaching significance (Fig. 7B). To further characterise the properties and distribution of CD8 + T cells, we used multiplex chromogenic IHC to stain CD3 in combination with CD8, together with 4G8 for Aβ plaques and Iba1 for microglia in some CTL and AD hippocampal samples (CTL n = 3; AD n = 6). With this approach, we detected only single positive cells, mainly CD8 + and sometimes CD3+ (Fig. 7C-D). In CTL hippocampus, CD8 + cells were strongly localised in the BV, with occasional presence in the parenchyma, close to microglia. In AD, we again observed a slight increase in parenchymal CD8 + cells, sometimes infiltrating in hotspots, but with very little to no association with PaM, Aβ plaques or CoM (Fig. 7D). Similar observations were made in the entorhinal and temporal cortices of AD samples (Supplementary Fig. 7A).
We next examined the presence of macrophages as DSP revealed shared transcripts between PaM and peripheral immune cells. As CD163 was one of the proteins enriched in PaM and the PaM surroundings, we investigated its distribution using two multiplex chromogenic IHC strategies. When we coupled Iba1 staining with CD163, we observed, in CTL samples (n = 3), that CD163 + perivascular cells were regularly distributed along the BV, while Iba1 + cells populated the parenchyma. Few Iba1 + CD163 + cells were detected lining the capillaries (Fig. 7E). In contrast, AD samples (n = 4) exhibited a notable increase in CD163 + cells within the parenchyma, often regrouped among PaM-like accumulations with several CD163 + cells being double positive for Iba1(Fig. 7E). CD163 + Iba1 + were abundant in the white matter of the stratum oriens. To further understand their relationship with Aβ plaques, we combined CD163, Iba1 and 4G8 staining on CTL and AD samples (5 CTL, 8 AD). We found a recurrent co-distribution of CD163+, Iba1 + and Iba1 + CD163 + around Aβ plaques in the hippocampus (Fig. 7F), entorhinal and temporal cortices (Supplementary Fig. 7B). Using digital pathology, we have quantified the association of Iba1+, CD163+, and CD163 + Iba1 + with Aβ plaques and their respective areas in the CA1 of AD cases. The number of Iba1 + PaM associated with plaque area was significantly higher than CD163 + and CD163 + Iba1 + which were the less frequent. We found a positive correlation between CA1 plaque size and the number of PaM (r = 0.273, ****) but low positive correlation for CD163+ (r = 0.128, **) and CD163 + Iba1+ (r = 0.111, *).