CoChat
Sign Up
PreprintCited by 169Neuroinflammation and Neurodegeneration MechanismsAlzheimer's disease research and treatmentsSingle-cell and spatial transcriptomics

The DOI matches its Crossref record; 67 of the 67 references it cites were found in a registry or an index.

A multi-pronged human microglia response to Alzheimer’s disease Aβ pathology

Open access
Mancuso, Renzo, Fattorelli, Nicola, Martinez-Muriana, Anna, Davis, Emma, Wolfs, Leen, Van Den Daele, Johanna, Geric, Ivana, Preman, Pranav, Serneels, Lutgarde, Poovathingal, Suresh, Balusu, Sriram, Verfaille, Catherine, Fiers, Mark, De Strooper, Bart·2022·Nature Neuroscience
DOI 10.1101/2022.07.07.499139PMC PMC11089003
Open this paper in CoChat — free accountOpen PDF

Abstract

Microglial activation and neuroinflammation are initial steps in the pathogenesis of Alzheimer’s disease (AD). However, studies in mouse models and human postmortem samples have yielded divergent results regarding microglia cell states relevant to AD. Here, we investigate 127,000 single cell expression profiles of human microglia isolated freshly from a xenotransplantation model for early AD. While human microglia adopt a disease-associated (DAM) profile, they display a much more pronounced HLA-cell state related to antigen presentation in response to amyloid plaques. In parallel, a distinctive pro-inflammatory cytokine and chemokine CRM response is mounted against oligomeric amyloid-β. TREM2 and, to a lesser extent, APOE polymorphisms, modulate the response of microglia to amyloid-β plaques, in contrast with the response to oligomeric Aβ. Specific polygenic risk genes are enriched in each branch of these multi-pronged response of human microglia to amyloid pathology (ARM). ARM responses can be captured in post-mortem studies when reanalyzed in light of this novel, comprehensive data set. In conclusion, therapeutic strategies targeting microglia in AD need to carefully assess how they affect the different cell states, as the overall balance between distinct microglial profiles might determine a protective or damaging outcome.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1101/2022.07.07.499139

Cites

67 of 67 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1 Human microglia display a complex, heterogeneous response to Aβ pathology. a, Experimental design used in this study. b, UMAP plot visualizing 138,577 single xenografted human microglial cells sorted from mouse brain (CD11b+hCD45+) ~6 months after transplantation. Cells are colored according to clusters identified: HM, RM, DAM, IRM, HLA, CRM-1 and CRM-2, and tCRM. The assignment of different clusters to distinct cell types or states is based on previous experimental data from our and other laboratories3,4,7,14 (see c,d and Extended Data Figs. 3 and 4). c, UMAP plots as in b, colored by the combined level of expression of groups of genes that characterize distinct microglial transcriptional states. d, Top ten most differentially expressed genes in each cluster (normalized expression scaled by gene is shown). FC, fold change; GO, gene ontology; MP, microglial precursor; SC, stem cell.
Fig. 1 Human microglia display a complex, heterogeneous response to Aβ pathology. a, Experimental design used in this study. b, UMAP plot visualizing 138,577 single xenografted human microglial cells sorted from mouse brain (CD11b+hCD45+) ~6 months after transplantation. Cells are colored according to clusters identified: HM, RM, DAM, IRM, HLA, CRM-1 and CRM-2, and tCRM. The assignment of different clusters to distinct cell types or states is based on previous experimental data from our and other laboratories3,4,7,14 (see c,d and Extended Data Figs. 3 and 4). c, UMAP plots as in b, colored by the combined level of expression of groups of genes that characterize distinct microglial transcriptional states. d, Top ten most differentially expressed genes in each cluster (normalized expression scaled by gene is shown). FC, fold change; GO, gene ontology; MP, microglial precursor; SC, stem cell.CC BY · PubMed Central
Fig. 2 Multipronged response of human microglia upon exposure to Aβ. a, Density plots displaying the average distribution of human microglia transplanted in AppNL-G-F (n = 23) and AppWT (n = 10) mice (density is normalized to sample size). b, Distribution and proportion of cells across all identified clusters (box plots are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum values. Each dot represents a single mouse. Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P < 0.05). c, Representative images and quantification of human microglia engrafted in the brain of AppNL-G-F mice at 6 months of age, labeled with human-specific antibodies for P2RY12 (HM microglia), CD9 (DAM), HLA-DR/DQ (HLA microglia), as well as X-34 for Aβ plaques. The shift in immunofluorescent signal in the proximity of Aβ plaques was performed using a modified Sholl analysis where the fluorescent intensity of microglial markers P2RY12, CD9 and HLA was measured through concentric rings (annuli) of increasing diameter surrounding the X-34 plaque center. Intensities of each channel were scaled for comparison using z-score normalization. Intensity over distance (µm) was plotted using Loess nonparametric regression in R with estimated standard error for each predicted value. For comparison of intensities near and distant from the plaque center, the means of the inner and outer three annuli were independently calculated (near, 0–10 mm, distant, 70–80 mm from the plaque; n = 3 per group; bar plots represent mean ± s.e.m.; unpaired t-test, one-tailed, alpha = 0.05, significance was set as P < 0.05). Scale bars, 50 μm. MFI, median fluorescent intensity.
Fig. 2 Multipronged response of human microglia upon exposure to Aβ. a, Density plots displaying the average distribution of human microglia transplanted in AppNL-G-F (n = 23) and AppWT (n = 10) mice (density is normalized to sample size). b, Distribution and proportion of cells across all identified clusters (box plots are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum values. Each dot represents a single mouse. Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P < 0.05). c, Representative images and quantification of human microglia engrafted in the brain of AppNL-G-F mice at 6 months of age, labeled with human-specific antibodies for P2RY12 (HM microglia), CD9 (DAM), HLA-DR/DQ (HLA microglia), as well as X-34 for Aβ plaques. The shift in immunofluorescent signal in the proximity of Aβ plaques was performed using a modified Sholl analysis where the fluorescent intensity of microglial markers P2RY12, CD9 and HLA was measured through concentric rings (annuli) of increasing diameter surrounding the X-34 plaque center. Intensities of each channel were scaled for comparison using z-score normalization. Intensity over distance (µm) was plotted using Loess nonparametric regression in R with estimated standard error for each predicted value. For comparison of intensities near and distant from the plaque center, the means of the inner and outer three annuli were independently calculated (near, 0–10 mm, distant, 70–80 mm from the plaque; n = 3 per group; bar plots represent mean ± s.e.m.; unpaired t-test, one-tailed, alpha = 0.05, significance was set as P < 0.05). Scale bars, 50 μm. MFI, median fluorescent intensity.CC BY · PubMed Central
Fig. 3 Human microglial transcriptional trajectory in response to Aβ pathology. a, Phenotypic trajectory followed by the human microglia after exposure to Aβ in vivo, obtained by an unbiased pseudotime ordering with Monocle 3. b, Distribution of cells from different host mouse genetic backgrounds (y axis) across the two main transcriptional trajectories, DAM and HLA (top panel) and CRM (bottom panel), colored by clusters shown in a. Note the shift in transcriptional states in AppNL-G-F versus AppWT mice (box plots are limited by lower and upper quartiles and midline indicates median; whiskers extend from the box to the smallest or largest value no further than 1.5 × interquartile range. Each dot represents a single cell, n = 10 mice (10,663 cells) for AppWT, n = 23 mice (32,436 cells) for AppNL-G-F. Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P < 0.05). c,d, Expression profile of selected genes across the HM–DAM–HLA pseudotime axis from the entire dataset (c), and divided into the different experimental groups: AppWT, AppNL-G-F at 3 months and AppNL-G-F at 6 months (d). e, Top ten most differentially expressed genes in each cluster divided by experimental group: AppWT, AppNL-G-F at 3 months and AppNL-F-G at 6 months (normalized expression scaled by gene is shown). f, Representative images of human microglia engrafted in the brain of AppNL-G-F mice at 6 months of age and labeled with human-specific CD9 (DAM) and HLA-DR/DQ (HLA microglia) antibodies. Scale bars, 25 μm (upper panel) and 50 μm (lower panel). APP–NLGF_3m, AppNL-G-F 3 months; APP–NLGF_6m, AppNL-G-F 6 months; APP–WT_6m, AppWT 6 months.
Fig. 3 Human microglial transcriptional trajectory in response to Aβ pathology. a, Phenotypic trajectory followed by the human microglia after exposure to Aβ in vivo, obtained by an unbiased pseudotime ordering with Monocle 3. b, Distribution of cells from different host mouse genetic backgrounds (y axis) across the two main transcriptional trajectories, DAM and HLA (top panel) and CRM (bottom panel), colored by clusters shown in a. Note the shift in transcriptional states in AppNL-G-F versus AppWT mice (box plots are limited by lower and upper quartiles and midline indicates median; whiskers extend from the box to the smallest or largest value no further than 1.5 × interquartile range. Each dot represents a single cell, n = 10 mice (10,663 cells) for AppWT, n = 23 mice (32,436 cells) for AppNL-G-F. Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P < 0.05). c,d, Expression profile of selected genes across the HM–DAM–HLA pseudotime axis from the entire dataset (c), and divided into the different experimental groups: AppWT, AppNL-G-F at 3 months and AppNL-G-F at 6 months (d). e, Top ten most differentially expressed genes in each cluster divided by experimental group: AppWT, AppNL-G-F at 3 months and AppNL-F-G at 6 months (normalized expression scaled by gene is shown). f, Representative images of human microglia engrafted in the brain of AppNL-G-F mice at 6 months of age and labeled with human-specific CD9 (DAM) and HLA-DR/DQ (HLA microglia) antibodies. Scale bars, 25 μm (upper panel) and 50 μm (lower panel). APP–NLGF_3m, AppNL-G-F 3 months; APP–NLGF_6m, AppNL-G-F 6 months; APP–WT_6m, AppWT 6 months.CC BY · PubMed Central
Fig. 4 TREM2 and APOE differentially modulate the transition to DAM and HLA states. a,d, Density plots displaying the average distribution of human H9-WT (n = 2) and H9-TREM2−/− (n = 3) (a), and human APOE3/3 (n = 6) and APOE−/− (n = 7) (d) microglia transplanted in AppNL-G-F mice. Density is normalized by sample size. b,e, Distribution and proportion of cells across all identified clusters for H9-WT (n = 2) and H9-TREM2−/− (n = 3) (b), and APOE3/3 (n = 6) and APOE−/− (n = 7) (e) transplanted microglia. Dots represent single mice. c,f, Phenotypic trajectory followed by H9-WT (n = 2 mice and 3,282 cells) and H9-TREM2−/− (n = 3 mice and 1,301 cells) (c), and APOE3/3 (n = 6 mice and 7,894 cells) and APOE−/− (n = 7 mice and 14,012 cells) (f) human microglia obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y axis) at different stages of the pseudotime trajectory (x axis), colored as shown in Fig. 1a. Dots represent single cells. g, Correlation analysis of the logFC in microglia transplanted in AppNL-G-F ApoE−/− versus AppNL-G-F mice (y axis from Fig. 4) and H9-TREM2−/− versus H9-WT from AppNL-G-F mice (x axis) (Pearson’s correlation, R = 0.17, differentially expressed genes adjusted using Bonferroni correction and colored as in Fig. 1b; ‘activated’ indicates that differential expression was performed comparing cell states reactive to pathology, excluding homeostatic or transitioning clusters). h, Correlation analysis of the logFC in TREM2−/− (x axis) and APOE−/− (y axis) microglia transplanted in AppNL-G-F mice (Pearson’s correlation, R = −0.11, differentially expressed genes adjusted using Bonferroni correction and colored as in Fig. 1b; ‘activated’ indicates that differential expression was performed comparing reactive cell states, excluding homeostatic or transitioning clusters). Box plots in b,c,e,f are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum value (b,e) or extend from the box to the smallest or largest value no further than 1.5× interquartile range (c,f). Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P < 0.05.
Fig. 4 TREM2 and APOE differentially modulate the transition to DAM and HLA states. a,d, Density plots displaying the average distribution of human H9-WT (n = 2) and H9-TREM2−/− (n = 3) (a), and human APOE3/3 (n = 6) and APOE−/− (n = 7) (d) microglia transplanted in AppNL-G-F mice. Density is normalized by sample size. b,e, Distribution and proportion of cells across all identified clusters for H9-WT (n = 2) and H9-TREM2−/− (n = 3) (b), and APOE3/3 (n = 6) and APOE−/− (n = 7) (e) transplanted microglia. Dots represent single mice. c,f, Phenotypic trajectory followed by H9-WT (n = 2 mice and 3,282 cells) and H9-TREM2−/− (n = 3 mice and 1,301 cells) (c), and APOE3/3 (n = 6 mice and 7,894 cells) and APOE−/− (n = 7 mice and 14,012 cells) (f) human microglia obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y axis) at different stages of the pseudotime trajectory (x axis), colored as shown in Fig. 1a. Dots represent single cells. g, Correlation analysis of the logFC in microglia transplanted in AppNL-G-F ApoE−/− versus AppNL-G-F mice (y axis from Fig. 4) and H9-TREM2−/− versus H9-WT from AppNL-G-F mice (x axis) (Pearson’s correlation, R = 0.17, differentially expressed genes adjusted using Bonferroni correction and colored as in Fig. 1b; ‘activated’ indicates that differential expression was performed comparing cell states reactive to pathology, excluding homeostatic or transitioning clusters). h, Correlation analysis of the logFC in TREM2−/− (x axis) and APOE−/− (y axis) microglia transplanted in AppNL-G-F mice (Pearson’s correlation, R = −0.11, differentially expressed genes adjusted using Bonferroni correction and colored as in Fig. 1b; ‘activated’ indicates that differential expression was performed comparing reactive cell states, excluding homeostatic or transitioning clusters). Box plots in b,c,e,f are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum value (b,e) or extend from the box to the smallest or largest value no further than 1.5× interquartile range (c,f). Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P < 0.05.CC BY · PubMed Central
Fig. 5 Human microglia display a differential response to Aβ plaques and Aβo. a,d, Density plots displaying the average distribution of human microglia transplanted in AppNL-G-F (n = 11) and AppNL-G-F ApoE−/− (n = 6) mice (a), and human microglia transplanted in AppWT mice and challenged with scrambled peptide (Scr, n = 7) or Aβo (n = 13) (d). Density is normalized by sample size. b,e, Distribution and proportion of cells across all identified clusters for microglia transplanted in AppNL-G-F (n = 11) and AppNL-G-F ApoE−/− (n = 6) mice (b), and Scr/Aβo injected mice (e). Dots represent single mice. c,f, Phenotypic trajectory followed by human microglia transplanted in AppNL-G-F (n = 11 and 14,649 cells) and AppNL-G-F ApoE−/− (n = 6 and 7,738 cells) mice (c), and challenged with scrambled peptide (Scr, n = 7 and 10,967 cells) or Aβo (n = 13 and 18,691 cells) (f), obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells from different mouse hosts (y axis) at different stages of the pseudotime trajectory (x axis), colored as shown in Fig. 1a. Dots represent single cells. g, Correlation analysis of the logFC in microglia either challenged with Aβο versus Src (y axis) or transplanted in AppNL-G-F ApoE−/− versus AppNL-G-F mice (x axis) (Pearson’s correlation, R = 0.05, differentially expressed genes were adjusted using Bonferroni correction and colored according to clusters in Fig. 1a). h, Correlation analysis of the logFC in microglia either challenged with Aβ soluble aggregates (Aβο) versus scrambled peptide (Src, y axis) or transplanted in AppNL-G-F mice (x axis) (Pearson’s correlation, R = 0.06, differentially expressed genes adjusted using Bonferroni correction). Box plots in b,c,e,f are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum values (b,e) or extend from the box to the smallest or largest value no further than 1.5× interquartile range (c,f). Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P < 0.05. Scr, scrambled peptide.
Fig. 5 Human microglia display a differential response to Aβ plaques and Aβo. a,d, Density plots displaying the average distribution of human microglia transplanted in AppNL-G-F (n = 11) and AppNL-G-F ApoE−/− (n = 6) mice (a), and human microglia transplanted in AppWT mice and challenged with scrambled peptide (Scr, n = 7) or Aβo (n = 13) (d). Density is normalized by sample size. b,e, Distribution and proportion of cells across all identified clusters for microglia transplanted in AppNL-G-F (n = 11) and AppNL-G-F ApoE−/− (n = 6) mice (b), and Scr/Aβo injected mice (e). Dots represent single mice. c,f, Phenotypic trajectory followed by human microglia transplanted in AppNL-G-F (n = 11 and 14,649 cells) and AppNL-G-F ApoE−/− (n = 6 and 7,738 cells) mice (c), and challenged with scrambled peptide (Scr, n = 7 and 10,967 cells) or Aβo (n = 13 and 18,691 cells) (f), obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells from different mouse hosts (y axis) at different stages of the pseudotime trajectory (x axis), colored as shown in Fig. 1a. Dots represent single cells. g, Correlation analysis of the logFC in microglia either challenged with Aβο versus Src (y axis) or transplanted in AppNL-G-F ApoE−/− versus AppNL-G-F mice (x axis) (Pearson’s correlation, R = 0.05, differentially expressed genes were adjusted using Bonferroni correction and colored according to clusters in Fig. 1a). h, Correlation analysis of the logFC in microglia either challenged with Aβ soluble aggregates (Aβο) versus scrambled peptide (Src, y axis) or transplanted in AppNL-G-F mice (x axis) (Pearson’s correlation, R = 0.06, differentially expressed genes adjusted using Bonferroni correction). Box plots in b,c,e,f are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum values (b,e) or extend from the box to the smallest or largest value no further than 1.5× interquartile range (c,f). Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P < 0.05. Scr, scrambled peptide.CC BY · PubMed Central
Fig. 6 Single-microglial nuclei from human postmortem brain. a–e, Human snRNA-seq datasets from Gerrits et al.20 (n = 176,136) (a), Sayed et al.35 (n = 28,767) (b), Zhou et al.18 (n = 3,978) (c), Olah et al.36 (n = 16,242) (d) and Sun et al.21 (n = 194,000) (e) re-analyzed to reproduce their original embeddings and annotated with our xenotransplanted microglial profiles as in Fig. 1b. f, Pairwise Pearson correlation between logFC of all differentially expressed genes (logFC cut-off set at 0.25, P < 0.05) of each microglial subtype and logFC of all differentially expressed genes (P < 0.05) of clusters from each human snRNA-seq study, with significance set at a P-adjusted value < 0.05 (*P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001). Only positive correlations are depicted here (r > 0). Additional correlations are shown in Extended Data Figs. 6 and 7. scRNA-seq, single cell RNA-sequencing.
Fig. 6 Single-microglial nuclei from human postmortem brain. a–e, Human snRNA-seq datasets from Gerrits et al.20 (n = 176,136) (a), Sayed et al.35 (n = 28,767) (b), Zhou et al.18 (n = 3,978) (c), Olah et al.36 (n = 16,242) (d) and Sun et al.21 (n = 194,000) (e) re-analyzed to reproduce their original embeddings and annotated with our xenotransplanted microglial profiles as in Fig. 1b. f, Pairwise Pearson correlation between logFC of all differentially expressed genes (logFC cut-off set at 0.25, P < 0.05) of each microglial subtype and logFC of all differentially expressed genes (P < 0.05) of clusters from each human snRNA-seq study, with significance set at a P-adjusted value < 0.05 (*P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001). Only positive correlations are depicted here (r > 0). Additional correlations are shown in Extended Data Figs. 6 and 7. scRNA-seq, single cell RNA-sequencing.CC BY · PubMed Central
Fig. 7 AD genetic risk modifies the response of human microglia to Aβ pathology. a,e, Density plots displaying the average distribution of human H9-WT (n = 2) and H9-TREM2R47H (n = 3) (a), and APOE2/2 (n = 6), APOE3/3 (n = 6) and APO4/4 (n = 6) human microglia (e) transplanted in AppNL-G-F mice. b,f, Distribution and percentage of cells across all identified clusters for H9-WT (n = 2) and H9-TREM2R47H (n = 3) (b), and APOE2/2 (n = 6), APOE3/3 (n = 6) and APO4/4 (n = 6) (f). Dots represent single mice. c,g, Phenotypic trajectory followed by H9-WT (n = 2) and H9-TREM2R47H (n = 3) (c), and APOE2/2 (n = 6), APOE3/3 (n = 6) and APO4/4 (n = 6) (g) human microglia obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y axis) over the binned pseudotime trajectory (x axis), colored by genotypes shown in a. Dots represent single mice (*P < 0.05). d,h, Correlation of the logFC in TREM2R47H versus H9-WT (y axis) and H9-TREM2−/− versus H9-WT (x axis) (d), and TREM2R47H versus H9-WT (y axis) and APOE4/4 versus APOE3/3 (x axis) (h) microglia transplanted in AppNL-G-F mice (Pearson’s correlation, R = 0.17, differentially expressed genes were adjusted using Bonferroni correction and colored as in Fig. 1b). i,j, Comparison of transcriptomic profiles of APOE4/4 (versus APOE3/3) versus APOE−/− (versus APOE3/3) (i), and APOE4/4 (versus APOE2/2) versus APOE−/− (versus APOE2/2) (j) microglia transplanted in AppNL-G-F mice (Pearson’s correlation, differentially expressed genes were adjusted using Bonferroni correction and colored according to clusters in Fig. 1b). Box plots in b,c,f,g are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum values (b,f) or extend from the box to the smallest or largest value no further than 1.5 × interquartile range (c,g). Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P < 0.05 (b,c); or one-way ANOVA with Tukey’s multiple comparisons as post hoc test, alpha = 0.05, significance was set as P-adjusted value < 0.05 (f,g). In d,h,i, ‘activated’ indicates that differential expression was performed comparing reactive cell states, excluding homeostatic or transitioning clusters.
Fig. 7 AD genetic risk modifies the response of human microglia to Aβ pathology. a,e, Density plots displaying the average distribution of human H9-WT (n = 2) and H9-TREM2R47H (n = 3) (a), and APOE2/2 (n = 6), APOE3/3 (n = 6) and APO4/4 (n = 6) human microglia (e) transplanted in AppNL-G-F mice. b,f, Distribution and percentage of cells across all identified clusters for H9-WT (n = 2) and H9-TREM2R47H (n = 3) (b), and APOE2/2 (n = 6), APOE3/3 (n = 6) and APO4/4 (n = 6) (f). Dots represent single mice. c,g, Phenotypic trajectory followed by H9-WT (n = 2) and H9-TREM2R47H (n = 3) (c), and APOE2/2 (n = 6), APOE3/3 (n = 6) and APO4/4 (n = 6) (g) human microglia obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y axis) over the binned pseudotime trajectory (x axis), colored by genotypes shown in a. Dots represent single mice (*P < 0.05). d,h, Correlation of the logFC in TREM2R47H versus H9-WT (y axis) and H9-TREM2−/− versus H9-WT (x axis) (d), and TREM2R47H versus H9-WT (y axis) and APOE4/4 versus APOE3/3 (x axis) (h) microglia transplanted in AppNL-G-F mice (Pearson’s correlation, R = 0.17, differentially expressed genes were adjusted using Bonferroni correction and colored as in Fig. 1b). i,j, Comparison of transcriptomic profiles of APOE4/4 (versus APOE3/3) versus APOE−/− (versus APOE3/3) (i), and APOE4/4 (versus APOE2/2) versus APOE−/− (versus APOE2/2) (j) microglia transplanted in AppNL-G-F mice (Pearson’s correlation, differentially expressed genes were adjusted using Bonferroni correction and colored according to clusters in Fig. 1b). Box plots in b,c,f,g are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum values (b,f) or extend from the box to the smallest or largest value no further than 1.5 × interquartile range (c,g). Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P < 0.05 (b,c); or one-way ANOVA with Tukey’s multiple comparisons as post hoc test, alpha = 0.05, significance was set as P-adjusted value < 0.05 (f,g). In d,h,i, ‘activated’ indicates that differential expression was performed comparing reactive cell states, excluding homeostatic or transitioning clusters.CC BY · PubMed Central
Fig. 8 AD genetic risk genes are differentially expressed in human microglial cell states and modulated by Aβ pathology. a, Analysis of genome-wide association studies (GWAS) genes enrichment in xenotransplanted microglia. The black bars represent the number of cells (in thousands) with detectable expression (≥1 read per cell) for each candidate gene from Bellenguez et al.42. The heatmap summarizes the deregulated expression (logFC, color scale) of these genes across cell states (each cluster compared with all others), as well as after exposure to Aβ plaque pathology, upon injection of soluble Aβ aggregates or altering the genetic background of the mice or the transplanted cells. The genes are ranked in rows based on hierarchical clustering. We identify three sets of genes that display a common profile across cell states (based on their enrichment in the specific microglial phenotypic transcriptional states HM, DAM and HLA, and CRM), Aβ pathology and genetic risk, and we group these profiles as: microglial homeostasis, plaque-induced genes and soluble aggregates-induced genes. The remaining genes did not show a clear enrichment in cell states or other conditions. All differential expressions were significant after adjusting P values using Bonferroni correction (FDR < 0.05). Only genes that are significantly changing in at least one of the tested conditions are reported, see Supplementary Table 4 and Extended Data Fig. 10 for further details. b, Illustration of the complex microglial AD genetics by cell-state profiles and driven by different Aβ pathologies as found in a and Extended Data Fig. 10.
Fig. 8 AD genetic risk genes are differentially expressed in human microglial cell states and modulated by Aβ pathology. a, Analysis of genome-wide association studies (GWAS) genes enrichment in xenotransplanted microglia. The black bars represent the number of cells (in thousands) with detectable expression (≥1 read per cell) for each candidate gene from Bellenguez et al.42. The heatmap summarizes the deregulated expression (logFC, color scale) of these genes across cell states (each cluster compared with all others), as well as after exposure to Aβ plaque pathology, upon injection of soluble Aβ aggregates or altering the genetic background of the mice or the transplanted cells. The genes are ranked in rows based on hierarchical clustering. We identify three sets of genes that display a common profile across cell states (based on their enrichment in the specific microglial phenotypic transcriptional states HM, DAM and HLA, and CRM), Aβ pathology and genetic risk, and we group these profiles as: microglial homeostasis, plaque-induced genes and soluble aggregates-induced genes. The remaining genes did not show a clear enrichment in cell states or other conditions. All differential expressions were significant after adjusting P values using Bonferroni correction (FDR < 0.05). Only genes that are significantly changing in at least one of the tested conditions are reported, see Supplementary Table 4 and Extended Data Fig. 10 for further details. b, Illustration of the complex microglial AD genetics by cell-state profiles and driven by different Aβ pathologies as found in a and Extended Data Fig. 10.CC BY · PubMed Central
Extended Data Fig. 1 Comparison of human microglia transplanted in AppWT and Rag2−/−hCSFKI mice, and primary microglia isolated from human surgical resections published before14. a, UMAP plot of the 39,538 cells passing quality control, coloured by annotated cell states after removal of CAMs, other myeloid, low quality and proliferating clusters. b, UMAP plots as in a, coloured by sample. c, UMAP plots as in a, split between the three experimental groups. d, Box plot showing proportion of cells across all identified clusters (box plots are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum value. Each dot represents a single replicate, n = 13 for AppWT, n = 1 for Rag2−/−/hCSF1KI, n = 7 for primary MG. One-way ANOVA with Tukey’s multiple comparisons as post-hoc test, alpha = 0.05, significance was set as P-adjusted value < 0.05). e, Heatmap displaying the top20 most upregulated genes in each cluster as in a. f, Heatmap depicting the expression levels of homeostatic microglial markers across all experimental groups (normalized expression scaled by gene is shown).
Extended Data Fig. 1 Comparison of human microglia transplanted in AppWT and Rag2−/−hCSFKI mice, and primary microglia isolated from human surgical resections published before14. a, UMAP plot of the 39,538 cells passing quality control, coloured by annotated cell states after removal of CAMs, other myeloid, low quality and proliferating clusters. b, UMAP plots as in a, coloured by sample. c, UMAP plots as in a, split between the three experimental groups. d, Box plot showing proportion of cells across all identified clusters (box plots are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum value. Each dot represents a single replicate, n = 13 for AppWT, n = 1 for Rag2−/−/hCSF1KI, n = 7 for primary MG. One-way ANOVA with Tukey’s multiple comparisons as post-hoc test, alpha = 0.05, significance was set as P-adjusted value < 0.05). e, Heatmap displaying the top20 most upregulated genes in each cluster as in a. f, Heatmap depicting the expression levels of homeostatic microglial markers across all experimental groups (normalized expression scaled by gene is shown).CC BY · PubMed Central
Extended Data Fig. 2 Extended analysis of the cell state transition in AppNL-G-F and AppWT mice. a, Volcano plot showing a paired comparison between AppNL-G-F and AppWT mice (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset). The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (ns = not significant). b, Volcano plot showing a paired comparison between AppNL-G-F and AppWT mice (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset), coloured by clusters as in Fig. 1 (HM = red squares, CRMs (CRM-1 & CRM-2) = ochre circles, DAM = olive triangles, HLA = cyan diamonds, other markers = grey circles). The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (ns = not significant). c, Phenotypic trajectory followed by the human microglia obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y-axis) present in each cell state as in Fig. 1a (rows) at different stages of the binned pseudotime trajectory (x-axis) in AppNL-G-F and AppWT mice (box plots are limited by lower and upper quartiles and midline indicates median; whiskers extend from the box to the smallest or largest value no further than 1.5*inter-quartile range. Each dot represents a single mouse, n = 10 for AppWT, n = 23 AppNL-G-F. Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P-value < 0.05; *p < 0.05, **p < 0.01, ****p < 0.0001). d, Representative images of human microglia engrafted in the brain of AppNL-G-F mice at 6 months of age and labelled with human specific antibodies FTH1 (RM). Human P2RY12 was stained in green, and amyloid-β plaques were stained with X34 in white. Scale bars are 25 μm (top row) and 50 μm (bottom row).
Extended Data Fig. 2 Extended analysis of the cell state transition in AppNL-G-F and AppWT mice. a, Volcano plot showing a paired comparison between AppNL-G-F and AppWT mice (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset). The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (ns = not significant). b, Volcano plot showing a paired comparison between AppNL-G-F and AppWT mice (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset), coloured by clusters as in Fig. 1 (HM = red squares, CRMs (CRM-1 & CRM-2) = ochre circles, DAM = olive triangles, HLA = cyan diamonds, other markers = grey circles). The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (ns = not significant). c, Phenotypic trajectory followed by the human microglia obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y-axis) present in each cell state as in Fig. 1a (rows) at different stages of the binned pseudotime trajectory (x-axis) in AppNL-G-F and AppWT mice (box plots are limited by lower and upper quartiles and midline indicates median; whiskers extend from the box to the smallest or largest value no further than 1.5*inter-quartile range. Each dot represents a single mouse, n = 10 for AppWT, n = 23 AppNL-G-F. Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P-value < 0.05; *p < 0.05, **p < 0.01, ****p < 0.0001). d, Representative images of human microglia engrafted in the brain of AppNL-G-F mice at 6 months of age and labelled with human specific antibodies FTH1 (RM). Human P2RY12 was stained in green, and amyloid-β plaques were stained with X34 in white. Scale bars are 25 μm (top row) and 50 μm (bottom row).CC BY · PubMed Central
Extended Data Fig. 3 Extended comparison of the response of human and mouse microglia to amyloid pathology. a, UMAP plot of the 3,803 mouse microglia isolated from AppNL-G-F and Rag2−/− Il2rγ−/− hCSF1KI AppNL-G-F passing quality control, coloured by annotated cell states and after removal of CAMs, other myeloid, low quality and proliferating clusters. b, Distribution and proportion of cells across all identified clusters (box plots are limited by lower and upper quartiles and mid line indicates median; whiskers show minimum to maximum value. Each dot represents a single mouse, n = 10 for xMG AppWT 6 m, n = 23 for xMG AppNL-G-F 6 m, n = 2 for AppNL-G-F). One-way ANOVA with Tukey’s multiple comparisons as post-hoc test, alpha = 0.05, significance was set as P-adjusted value < 0.05). For cell states that only had 2 groups (tCRM and CRM-1), we applied an unpaired t-test with Welch’s correction, alpha = 0.05, two-tailed, significance was set as P-value < 0.05. Note that whereas mouse cells polarize to DAM, HLA and CRM, proportions do not increase. c, Correlation analysis of the logFC in DAM2 vs. DAM1 clusters defined by Keren-Shaul et al.4 (y-axis) and HLA vs. DAM clusters defined in this study (x-axis). Pearson’s correlation, R = −0.44, differentially expressed genes were adjusted using Bonferroni correction and colored according to clusters in Fig. 1a. d, Overlap in the number of genes between our DAM signature and that from mouse systems (KarenShaul DAM and SalaFrigerio ARM), other human microglia xenotransplantation studies (Hasselman DAM), and postmortem single nuclei samples (Gerrits AD1, and Olah DAM). The arrow depicts the direction used to calculate the percentage overlap.
Extended Data Fig. 3 Extended comparison of the response of human and mouse microglia to amyloid pathology. a, UMAP plot of the 3,803 mouse microglia isolated from AppNL-G-F and Rag2−/− Il2rγ−/− hCSF1KI AppNL-G-F passing quality control, coloured by annotated cell states and after removal of CAMs, other myeloid, low quality and proliferating clusters. b, Distribution and proportion of cells across all identified clusters (box plots are limited by lower and upper quartiles and mid line indicates median; whiskers show minimum to maximum value. Each dot represents a single mouse, n = 10 for xMG AppWT 6 m, n = 23 for xMG AppNL-G-F 6 m, n = 2 for AppNL-G-F). One-way ANOVA with Tukey’s multiple comparisons as post-hoc test, alpha = 0.05, significance was set as P-adjusted value < 0.05). For cell states that only had 2 groups (tCRM and CRM-1), we applied an unpaired t-test with Welch’s correction, alpha = 0.05, two-tailed, significance was set as P-value < 0.05. Note that whereas mouse cells polarize to DAM, HLA and CRM, proportions do not increase. c, Correlation analysis of the logFC in DAM2 vs. DAM1 clusters defined by Keren-Shaul et al.4 (y-axis) and HLA vs. DAM clusters defined in this study (x-axis). Pearson’s correlation, R = −0.44, differentially expressed genes were adjusted using Bonferroni correction and colored according to clusters in Fig. 1a. d, Overlap in the number of genes between our DAM signature and that from mouse systems (KarenShaul DAM and SalaFrigerio ARM), other human microglia xenotransplantation studies (Hasselman DAM), and postmortem single nuclei samples (Gerrits AD1, and Olah DAM). The arrow depicts the direction used to calculate the percentage overlap.CC BY · PubMed Central
Extended Data Fig. 4 Extended pseudotime analysis, including transplanted microglia isolated from 3 months old AppNL-G-F mice. a, UMAP plot of the 38,360 mouse microglia isolated from 3 and 6 months AppNL-G-F and AppWT mice passing quality control, coloured by annotated cell states and after removal of CAMs, other myeloid, low quality and proliferating clusters. b, Phenotypic trajectory followed by the human microglia after exposure to amyloid-β in vivo, obtained by an unbiased pseudotime ordering with Monocle 3. c, Density plots displaying the average distribution of human microglia transplanted in 3 months old AppNL-G-F (n = 11) 6 months old AppNL-G-F (n = 9) and AppWT (n = 8) mice (density is normalized to sample size). d, Box plot showing proportion of cells across all identified clusters (eEach dot represents a single mouse, n = 8 for AppWT, n = 11 for AppNL-G-F 3 m, n = 9 for AppNL-G-F 6 m). e, Distribution of cells from different host mouse genetic backgrounds and ages (y-axis) across the DAM and HLA transcriptional trajectories colored by clusters shown in a (Each dot represents a single cell, n = 8 mice (8,974 cells) for AppWT, n = 11 mice (10,822 cells) for AppNL-G-F 3 m, n = 9 mice (12,197 cells) for AppNL-G-F 6 m). f, Phenotypic trajectory followed by the human microglia obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y-axis) present in each cell state as in a (rows) at different stages of the binned pseudotime trajectory (x-axis) in 3 and 6 months old AppNL-G-F and AppWT mice (Each dot represents a single mouse, n = 8 for AppWT, n = 11 for AppNL-G-F 3 m, n = 9 for AppNL-G-F 6 m; *p < 0.05, ****p < 0.001). Box plots in d, e and f are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum value (d) or extend from the box to the smallest or largest value no further than 1.5*inter-quartile range (e and f). One-way ANOVA with Tukey’s multiple comparisons as post-hoc test, alpha = 0.05, significance was set as P-adjusted value < 0.05.
Extended Data Fig. 4 Extended pseudotime analysis, including transplanted microglia isolated from 3 months old AppNL-G-F mice. a, UMAP plot of the 38,360 mouse microglia isolated from 3 and 6 months AppNL-G-F and AppWT mice passing quality control, coloured by annotated cell states and after removal of CAMs, other myeloid, low quality and proliferating clusters. b, Phenotypic trajectory followed by the human microglia after exposure to amyloid-β in vivo, obtained by an unbiased pseudotime ordering with Monocle 3. c, Density plots displaying the average distribution of human microglia transplanted in 3 months old AppNL-G-F (n = 11) 6 months old AppNL-G-F (n = 9) and AppWT (n = 8) mice (density is normalized to sample size). d, Box plot showing proportion of cells across all identified clusters (eEach dot represents a single mouse, n = 8 for AppWT, n = 11 for AppNL-G-F 3 m, n = 9 for AppNL-G-F 6 m). e, Distribution of cells from different host mouse genetic backgrounds and ages (y-axis) across the DAM and HLA transcriptional trajectories colored by clusters shown in a (Each dot represents a single cell, n = 8 mice (8,974 cells) for AppWT, n = 11 mice (10,822 cells) for AppNL-G-F 3 m, n = 9 mice (12,197 cells) for AppNL-G-F 6 m). f, Phenotypic trajectory followed by the human microglia obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y-axis) present in each cell state as in a (rows) at different stages of the binned pseudotime trajectory (x-axis) in 3 and 6 months old AppNL-G-F and AppWT mice (Each dot represents a single mouse, n = 8 for AppWT, n = 11 for AppNL-G-F 3 m, n = 9 for AppNL-G-F 6 m; *p < 0.05, ****p < 0.001). Box plots in d, e and f are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum value (d) or extend from the box to the smallest or largest value no further than 1.5*inter-quartile range (e and f). One-way ANOVA with Tukey’s multiple comparisons as post-hoc test, alpha = 0.05, significance was set as P-adjusted value < 0.05.CC BY · PubMed Central
Extended Data Fig. 5 Extended pseudotime analysis, including transplanted microglia isolated from 3 months old AppNL-G-F mice (continuation). a, b, Heatmaps displaying the cells ordered by the (a) HM-DAM-HLA and (b) HM-CRM pseudotime axis, coloured by single-cell levels of expression of the top markers per cluster. c, d, Expression profile of selected genes across the (c) HM-DAM-HLA and (d) CRM pseudotime axis averaging all experimental groups in the dataset. e, Comparison of the expression profile of selected CRM markers across the trajectory. f, Expression profile of selected genes across the CRM pseudotime axis as in d divided into the different experimental groups: AppWT, AppNL-G-F 3 M and AppNL-F-G 6 M.
Extended Data Fig. 5 Extended pseudotime analysis, including transplanted microglia isolated from 3 months old AppNL-G-F mice (continuation). a, b, Heatmaps displaying the cells ordered by the (a) HM-DAM-HLA and (b) HM-CRM pseudotime axis, coloured by single-cell levels of expression of the top markers per cluster. c, d, Expression profile of selected genes across the (c) HM-DAM-HLA and (d) CRM pseudotime axis averaging all experimental groups in the dataset. e, Comparison of the expression profile of selected CRM markers across the trajectory. f, Expression profile of selected genes across the CRM pseudotime axis as in d divided into the different experimental groups: AppWT, AppNL-G-F 3 M and AppNL-F-G 6 M.CC BY · PubMed Central
Extended Data Fig. 6 Extended analysis on the differential impact of TREM2 and APOE deficiency in the transplanted microglia. a, Phenotypic trajectory followed by TREM2 deficient and control microglia, obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y-axis) present in each cell state as in Fig. 1b (rows) at different stages of the binned pseudotime trajectory (x-axis), colored by host genotypes (each dot represents a single mouse, n = 2 mice H9-WT, n = 3 mice H9-TREM2−/−). b, c, Volcano plots showing a paired comparison between TREM2−/− and WT microglia transplanted in AppNL-G-F, coloured by (b) cells genotype and (c) cluster as in Fig. 1b. The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset, ns = not significant). d, Phenotypic trajectory followed by APOE deficient microglia compared to isogenic APOE3/3 controls, obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y-axis) present in each cell state as in Fig. 1b (rows) at different stages of the binned pseudotime trajectory (x-axis), colored by host genotypes (each dot represents a single mouse, n = 6 APOE3/3 and n = 7 for APOE−/−). e, f, Volcano plots showing a paired comparison between APOE−/− and APOE3/3 microglia transplanted in AppNL-G-F, coloured by (e) cells genotype and (f) cluster as in Fig. 1b. The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset, ns = not significant). Box plots in a and d, are limited by lower and upper quartiles and midline indicates median; whiskers extend from the box to the smallest or largest value no further than 1.5*inter-quartile range. Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P-value < 0.05, *p < 0.05, **p < 0.01.
Extended Data Fig. 6 Extended analysis on the differential impact of TREM2 and APOE deficiency in the transplanted microglia. a, Phenotypic trajectory followed by TREM2 deficient and control microglia, obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y-axis) present in each cell state as in Fig. 1b (rows) at different stages of the binned pseudotime trajectory (x-axis), colored by host genotypes (each dot represents a single mouse, n = 2 mice H9-WT, n = 3 mice H9-TREM2−/−). b, c, Volcano plots showing a paired comparison between TREM2−/− and WT microglia transplanted in AppNL-G-F, coloured by (b) cells genotype and (c) cluster as in Fig. 1b. The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset, ns = not significant). d, Phenotypic trajectory followed by APOE deficient microglia compared to isogenic APOE3/3 controls, obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y-axis) present in each cell state as in Fig. 1b (rows) at different stages of the binned pseudotime trajectory (x-axis), colored by host genotypes (each dot represents a single mouse, n = 6 APOE3/3 and n = 7 for APOE−/−). e, f, Volcano plots showing a paired comparison between APOE−/− and APOE3/3 microglia transplanted in AppNL-G-F, coloured by (e) cells genotype and (f) cluster as in Fig. 1b. The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset, ns = not significant). Box plots in a and d, are limited by lower and upper quartiles and midline indicates median; whiskers extend from the box to the smallest or largest value no further than 1.5*inter-quartile range. Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P-value < 0.05, *p < 0.05, **p < 0.01.CC BY · PubMed Central
Extended Data Fig. 7 Extended analysis on the differential effects of amyloid-β plaques and soluble aggregates on microglia cell states. a, Phenotypic trajectory followed by the human microglia transplanted in AppNL-G-F ApoE−/− and AppNL-G-F, obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y-axis) present in each cell state as in Fig. 1b (rows) at different stages of the binned pseudotime trajectory (x-axis), colored by host genotypes (dots represent single mice, n = 11 AppNL-G-F, n = 6 AppNL-G-F ApoE−/−). b, c, Volcano plots showing a paired comparison between microglia transplanted in AppNL-G-F ApoE−/− and AppNL-G-F mice, coloured by (b) host genotype and (c) cluster as in Fig. 1a. The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset, ns = not significant). d, Phenotypic trajectory followed by the human microglia after injection of amyloid-b oligomers (Aβο) or scrambled peptide (Scr), obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y-axis) present in each cell state as in Fig. 1a (rows) at different stages of the binned pseudotime trajectory (x-axis), colored by cells genotypes (dots represent single mice, n = 7 for Scr, n = 13 Aβο). e, f, Volcano plots showing a paired comparison between microglia transplanted in AppWT mice and treated with Aβο or Scr, coloured by (e) treatment and (f) cluster as in Fig. 1a. The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset, ns = not significant). Box plots in a and d, are limited by lower and upper quartiles and midline indicates median; whiskers extend from the box to the smallest or largest value no further than 1.5*inter-quartile range. Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P-value < 0.05, *p < 0.05, **p < 0.01, ***p < 0.001,****p < 0.0001.
Extended Data Fig. 7 Extended analysis on the differential effects of amyloid-β plaques and soluble aggregates on microglia cell states. a, Phenotypic trajectory followed by the human microglia transplanted in AppNL-G-F ApoE−/− and AppNL-G-F, obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y-axis) present in each cell state as in Fig. 1b (rows) at different stages of the binned pseudotime trajectory (x-axis), colored by host genotypes (dots represent single mice, n = 11 AppNL-G-F, n = 6 AppNL-G-F ApoE−/−). b, c, Volcano plots showing a paired comparison between microglia transplanted in AppNL-G-F ApoE−/− and AppNL-G-F mice, coloured by (b) host genotype and (c) cluster as in Fig. 1a. The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset, ns = not significant). d, Phenotypic trajectory followed by the human microglia after injection of amyloid-b oligomers (Aβο) or scrambled peptide (Scr), obtained by an unbiased pseudotime ordering with Monocle 3. Proportion of cells (y-axis) present in each cell state as in Fig. 1a (rows) at different stages of the binned pseudotime trajectory (x-axis), colored by cells genotypes (dots represent single mice, n = 7 for Scr, n = 13 Aβο). e, f, Volcano plots showing a paired comparison between microglia transplanted in AppWT mice and treated with Aβο or Scr, coloured by (e) treatment and (f) cluster as in Fig. 1a. The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset, ns = not significant). Box plots in a and d, are limited by lower and upper quartiles and midline indicates median; whiskers extend from the box to the smallest or largest value no further than 1.5*inter-quartile range. Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P-value < 0.05, *p < 0.05, **p < 0.01, ***p < 0.001,****p < 0.0001.CC BY · PubMed Central
Extended Data Fig. 8 Extended analysis of the single nuclei from human post-mortem brains (continuation). a, UMAP plots as in Fig. 6, coloured by the combined level of expression of groups of genes that characterise distinct microglial transcriptional states from xenotransplanted microglia. b, Full set of pairwise Pearson correlations between logFC of all DE genes (logFC cut off set at 0.25, P-value < 0.05) of each microglia sub-type and logFC of all DE genes (p < 0.05) of clusters from each human snRNA-seq study, with significance set at a P-adjusted value < 0.05 (*p < =0.05, **p < =0.01, ***p < =0.001). c, Representative scatter plots showing the correlation analysis of the logFC between xenotransplanted and single nuclei post-mortem microglia. The labels are coloured by cell state as in Fig. 1.
Extended Data Fig. 8 Extended analysis of the single nuclei from human post-mortem brains (continuation). a, UMAP plots as in Fig. 6, coloured by the combined level of expression of groups of genes that characterise distinct microglial transcriptional states from xenotransplanted microglia. b, Full set of pairwise Pearson correlations between logFC of all DE genes (logFC cut off set at 0.25, P-value < 0.05) of each microglia sub-type and logFC of all DE genes (p < 0.05) of clusters from each human snRNA-seq study, with significance set at a P-adjusted value < 0.05 (*p < =0.05, **p < =0.01, ***p < =0.001). c, Representative scatter plots showing the correlation analysis of the logFC between xenotransplanted and single nuclei post-mortem microglia. The labels are coloured by cell state as in Fig. 1.CC BY · PubMed Central
Extended Data Fig. 9 Extended analysis of the impact of clinical mutations on microglial cell states. a, b, Volcano plots showing a paired comparison between TREM2R47H and TREM2-WT microglia transplanted in AppNL-G-F mice and coloured by (a) cells genotype or (b) cluster as in Fig. 1a. The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset, ns = not significant). c, d, Distribution of TREM2R47H and TREM2-WT cells across (c) DAM-HLA and (d) CRM transcriptional trajectories. Note the shift in transcriptional states in both the DAM-HLA and CRM axes (each dot represents a single cell, n = 2 (3,282 cells) H9-WT, n = 3 (5,845 cells) H9-TREM2R47H). e, f, Volcano plots showing a paired comparison between APOE4/4 and APOE3/3 microglia transplanted in AppNL-G-F mice and coloured by (e) cells genotype or (f) cluster as in Fig. 1a. The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset, ns = not significant). g, h, Distribution of APOE4 and APOE3 cells across (g) DAM and HLA and (h) CRM transcriptional trajectories (each dot represents a single mouse, n = 6 (7,894 cells) for APOE3/3, n = 6 (10,151 cells) for APO4/4). i, Distribution and proportion of cells across all identified clusters as in Fig. 5e, coloured by the microglial APOE genotype (each dot represents a single mouse, n = 7 for Scr, n = 13 Aβο). Box plots in c, d, g, h and i, are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum value (i) or extend from the box to the smallest or largest value no further than 1.5*inter-quartile range (c, d, g and h). Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P-value < 0.05.
Extended Data Fig. 9 Extended analysis of the impact of clinical mutations on microglial cell states. a, b, Volcano plots showing a paired comparison between TREM2R47H and TREM2-WT microglia transplanted in AppNL-G-F mice and coloured by (a) cells genotype or (b) cluster as in Fig. 1a. The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset, ns = not significant). c, d, Distribution of TREM2R47H and TREM2-WT cells across (c) DAM-HLA and (d) CRM transcriptional trajectories. Note the shift in transcriptional states in both the DAM-HLA and CRM axes (each dot represents a single cell, n = 2 (3,282 cells) H9-WT, n = 3 (5,845 cells) H9-TREM2R47H). e, f, Volcano plots showing a paired comparison between APOE4/4 and APOE3/3 microglia transplanted in AppNL-G-F mice and coloured by (e) cells genotype or (f) cluster as in Fig. 1a. The number of significant genes per condition is reported in brackets in the legend (Log2 Fold-Change threshold > 0.1). Adjusted p-value threshold < 0.05 (Wilcoxon rank-sum test, P-values adjusted with Bonferroni correction based on the total number of genes in the dataset, ns = not significant). g, h, Distribution of APOE4 and APOE3 cells across (g) DAM and HLA and (h) CRM transcriptional trajectories (each dot represents a single mouse, n = 6 (7,894 cells) for APOE3/3, n = 6 (10,151 cells) for APO4/4). i, Distribution and proportion of cells across all identified clusters as in Fig. 5e, coloured by the microglial APOE genotype (each dot represents a single mouse, n = 7 for Scr, n = 13 Aβο). Box plots in c, d, g, h and i, are limited by lower and upper quartiles and midline indicates median; whiskers show minimum to maximum value (i) or extend from the box to the smallest or largest value no further than 1.5*inter-quartile range (c, d, g and h). Unpaired t-test with Welch’s correction, two-tailed, alpha = 0.05, significance was set as P-value < 0.05.CC BY · PubMed Central
Extended Data Fig. 10 Extended analysis of GWAS genes enrichment in xenotransplanted microglia. a, Summary heatmap with the expression levels of the Top10 marker genes for each cluster (as in Fig. 1c), compared across all experimental groups used in this study. b, Expanded gene list with the canonical AD genes according to the latest GWAS study of Fig. 8a, with the addition of other genes of interest, that is the genes selected in our previous publication14 (see Supplementary Table 4 for the complete list). Raw counts (log scale) are reported on the x axis, while standard deviation of counts is shown on the y-axis. Dots are coloured by SCT-normalised expression (see Methods) and size reflects percentage of cells in our dataset where the gene is detected (see all values in Supplementary Table 4). All genes from (a) that are significantly changing in at least one of the tested conditions in this manuscript are reported in an extended heat map build as in Fig. 8a. The black bars represent the number of cells (in thousands) with detectable expression (>=1 read per cell) for each candidate gene. The heat map summarizes the deregulated expression (LogFC, colour scale) of these genes across cell states (each cluster compared to all others), as well as after exposure to amyloid-β plaque pathology, upon injection of soluble aggregates of Aβ, or altering the genetic background of the mice or the transplanted cells. The genes are ranked in rows based on hierarchical clustering. We identify 3 sets of genes that display a common profile across cell states (based on their enrichment in the specific microglial phenotypic transcriptional states HM, DAM and HLA, and CRM), amyloid-β pathology and genetic risk, and we group these profiles as: microglia Homeostasis, plaque-induced genes, and soluble aggregates-induced genes. The remaining genes did not show a clear enrichment in cell states or other conditions. All differential expressions were significant after adjusting P-values using Bonferroni correction (FDR < 0.05).
Extended Data Fig. 10 Extended analysis of GWAS genes enrichment in xenotransplanted microglia. a, Summary heatmap with the expression levels of the Top10 marker genes for each cluster (as in Fig. 1c), compared across all experimental groups used in this study. b, Expanded gene list with the canonical AD genes according to the latest GWAS study of Fig. 8a, with the addition of other genes of interest, that is the genes selected in our previous publication14 (see Supplementary Table 4 for the complete list). Raw counts (log scale) are reported on the x axis, while standard deviation of counts is shown on the y-axis. Dots are coloured by SCT-normalised expression (see Methods) and size reflects percentage of cells in our dataset where the gene is detected (see all values in Supplementary Table 4). All genes from (a) that are significantly changing in at least one of the tested conditions in this manuscript are reported in an extended heat map build as in Fig. 8a. The black bars represent the number of cells (in thousands) with detectable expression (>=1 read per cell) for each candidate gene. The heat map summarizes the deregulated expression (LogFC, colour scale) of these genes across cell states (each cluster compared to all others), as well as after exposure to amyloid-β plaque pathology, upon injection of soluble aggregates of Aβ, or altering the genetic background of the mice or the transplanted cells. The genes are ranked in rows based on hierarchical clustering. We identify 3 sets of genes that display a common profile across cell states (based on their enrichment in the specific microglial phenotypic transcriptional states HM, DAM and HLA, and CRM), amyloid-β pathology and genetic risk, and we group these profiles as: microglia Homeostasis, plaque-induced genes, and soluble aggregates-induced genes. The remaining genes did not show a clear enrichment in cell states or other conditions. All differential expressions were significant after adjusting P-values using Bonferroni correction (FDR < 0.05).CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 67 references

Research works printed but not linked: 2

  1. Sala Frigerio, C. et al. The Major Risk Factors for Alzheimer’s Disease: Age, Sex, and Genes Modulate the Microglia Response to Aβ Plaques. Cell Rep. 27, (2019). · DOI · Paper page
  2. 10.1038/s41593-019-0525-x · DOI · Paper page

Cited by

Sorted by citations

See all 169

Related

Metadata

Authors
Mancuso, R., Fattorelli, N., Martinez-Muriana, A., Davis, E., Wolfs, L., Van Den Daele, J., Geric, I., Preman, P., Serneels, L., Poovathingal, S., Balusu, S., Verfaille, C., Fiers, M., & De Strooper, B.
Affiliations
  • University of Antwerp (BE)
  • Vlaams Instituut voor Biotechnologie (BE)
  • UK Dementia Research Institute (GB)
Corresponding author
Mancuso, Renzo
ORCID
Venue
Nature Neuroscience
Publisher
openRxiv
Published
2022-07-07
Volume
27
Issue
5
Pages
886-900
ISSN
1097-6256 (ISSN-L)
DOI
10.1101/2022.07.07.499139
PMC
PMC11089003
OpenAlex
W4393225316
Cited by
169
References
67
Citations by year
2026: 76 · 2025: 57 · 2024: 30 · 2023: 3 · 2022: 3
FWCI
25.22
Type
journal article
Language
en
Open access
Licence
cc-by
Locations
  • Nature Neuroscience (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • UCL Discovery (University College London) (submittedVersion, cc-by)
  • Lirias (submittedVersion, cc-by)
Funders
  • Alzheimer's Association (AARF-20-684397)
  • BrightFocus Foundation (A2021034S)
  • UK Research and Innovation (MR/Y014847/1)
  • Alzheimer's Society
  • Queen Elisabeth Medical Foundation
  • European Commission (HORIZON2020, ERC-834682 CELLPHASE_AD, 834682, 2022-2025)
  • Fonds Wetenschappelijk Onderzoek (G0C9219N, G0K9422N, 1122423N, 12ZX621N, 1139520N, G056022N, S001221N)
  • KU Leuven
  • Vlaams Instituut voor Biotechnologie
  • Universiteit Antwerpen (FFB210307)
  • Vlaamse regering
  • UK Dementia Research Institute
  • Fondation pour la Recherche sur Alzheimer
  • Fondation Médicale Reine Elisabeth
  • Horizon 2020 Framework Programme
  • Medical Research Council (MR/Y014847/1)
  • H2020 Marie Skłodowska-Curie Actions (HORIZON2020)
  • HORIZON EUROPE European Research Council (101041867)
Topics
  • Life Sciences › Neuroscience › Neurology › Neuroinflammation and Neurodegeneration Mechanisms
  • Life Sciences › Neuroscience › Cellular and Molecular Neuroscience › Alzheimer's disease research and treatments
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
Keywords
microglia, Alzheimer's disease, amyloid-beta, xenotransplantation, single-cell transcriptomics, disease-associated microglia, HLA, antigen presentation, amyloid plaques, TREM2, apolipoprotein E, inflammatory cytokines, chemokines, amyloid-beta oligomers
MeSH terms
Alzheimer Disease, Animals, Apolipoproteins E, Brain, Humans, Membrane Glycoproteins, Mice, Transgenic, Receptors, Immunologic, Amyloid beta-Peptides, Microglia, Mice, Plaque, Amyloid, Transcriptome, Heterografts

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review