Cited by 120Alzheimer's disease research and treatmentsDementia and Cognitive Impairment ResearchNeuroinflammation and Neurodegeneration Mechanisms

The DOI matches its Crossref record; 105 of the 105 references it cites were found in a registry or an index; 73 of the 85 sources its claims cite were checked for support.

The APOE-R136S mutation protects against APOE4-driven tau pathology, neurodegeneration, and neuroinflammation

Open access

Abstract

Apolipoprotein E4 ( APOE4 ) is the strongest genetic risk factor for late-onset Alzheimer’s disease (LOAD), leading to earlier age of clinical onset and exacerbating pathologies. There is a critical need to identify protective targets. Recently, a rare APOE variant, APOE3-R136S (Christchurch), was found to protect against early-onset AD in a PSEN1-E280A carrier. In this study, we sought to determine if the R136S mutation also protects against APOE4-driven effects in LOAD. We generated tauopathy mouse and human iPSC-derived neuron models carrying human APOE4 with the homozygous or heterozygous R136S mutation. We found that the homozygous R136S mutation rescued APOE4-driven Tau pathology, neurodegeneration and neuroinflammation. The heterozygous R136S mutation partially protected against APOE4-driven neurodegeneration and neuroinflammation but not Tau pathology. Single-nucleus RNA sequencing revealed that the APOE4-R136S mutation increased disease-protective and diminished disease-associated cell populations in a gene dose-dependent manner. Thus, the APOE-R136S mutation protects against APOE4-driven AD pathologies, providing a target for therapeutic development against AD.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

The APOE4-R136S mutation protects against APOE4-driven tau pathology, neurodegeneration, and neuroinflammation in PS19 tauopathy mice and human iPSC-derived neurons, with homozygous R136S generally required for full rescue. Heterozygous R136S partially protects against neurodegeneration and neuroinflammation but not tau pathology, and snRNA-seq shows gene-dose-dependent increases in disease-protective and decreases in disease-associated cell populations.

Key findings

  • Homozygous E4-S/S reduced hippocampal p-Tau: PS19-E4 mice had roughly two-fold higher AT8+ coverage than PS19-E3 mice, while PS19-E4-S/S fell to PS19-E3 levels; PS19-E4-R/S did not show this reduction. In human neurons, E4-S/S showed ~40% lower p-Tau by western blot and ~35% lower PHF1+ p-Tau by immunostaining versus E4, whereas E4-R/S showed no reduction.
  • Neurodegeneration was rescued by homozygous R136S: PS19-E4-S/S hippocampal atrophy was reduced to WT/PS19-E3 levels, with partial rescue in PS19-E4-R/S; DG GC layer thickness increased ~20% in PS19-E3 and PS19-E4-S/S over PS19-E4, while PS19-E4-R/S showed a non-significant 12% increase. R136S copy number correlated with hippocampal volume and posterior lateral ventricle volume.
  • Neuroinflammation was reduced in a gene-dose-dependent manner: PS19-E4-S/S and PS19-E3 had lower GFAP+ and S100β+ astrocytosis and lower Iba1+ and CD68+ microgliosis than PS19-E4; PS19-E4-R/S showed intermediate protection, with significantly higher Iba1+ and CD68+ cell numbers than PS19-E4-S/S.
  • snRNA-seq of 175,080 nuclei and 38 clusters showed E4-S/S increased disease-protective neuronal clusters 1, 6, 7, and 28 and homeostatic astrocyte subcluster 3 and microglia subclusters 2 and 11, while decreasing disease-associated oligodendrocytes (cluster 9), astrocytes (subclusters 5 and 7), and microglia (subcluster 8). E4-R/S also reduced DAA subcluster 7 (LOR −1.57 ± 0.57) and DAM subcluster 8 (LOR −2.58 ± 1.02), but less than E4-S/S (subcluster 7 LOR −2.75 ± 0.62; subcluster 8 LOR −3.36 ± 1.17).

Methods

  • CRISPR–Cas9 was used to introduce R136S into human APOE4 knock-in mice, which were crossed with PS19 Tau-P301S mice to generate PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S groups; mice were analyzed at ~10 months, with some 6-month comparisons, using both sexes.
  • Isogenic hiPSC lines from an APOE4/4 Alzheimer’s disease patient were edited to E4-S/S and E4-R/S, alongside E4, E3, and APOE-knockout lines, and differentiated into neurons for western blot, immunostaining, Tau-488 uptake flow cytometry with/without heparin, long-term heparin treatment, and E4-conditioned medium experiments.
  • Hippocampal snRNA-seq from 10-month PS19 mice was analyzed with Seurat clustering, differential expression, KEGG pathway enrichment, and generalized linear mixed-effects models for associations with genotype and histopathology.

Limitations

  • Stated in the paper: the models do not exhibit extensive Aβ pathology, so the study does not address whether APOE-R136S protects against Aβ-mediated neuroinflammatory and neurodegenerative effects.
  • Stated in the paper: the study focused on APOE4; future studies are needed to test whether R136S also protects in the context of APOE3.
  • Evident from the text: mouse tissue studies used one cohort; snRNA-seq had 3–4 mice per group and one PS19-E3 sample was excluded after quality control; hiPSC experiments used biological replicates over two or three independent differentiations; no randomization was used for group assignment, and snRNA-seq investigators were not blinded.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1038/s41593-023-01480-8

Cites

105 of 105 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Fig. 1 Homozygous R136S mutation rescues APOE4-promoted Tau pathology in tauopathy mice. a, Schematic of CRISPR–Cas-9-mediated gene editing strategy to generate human APOE4-R136S knock-in mice. b, Representative images of p-Tau immunostaining in hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice with the AT8 monoclonal antibody. c, Quantification of the percent AT8 coverage area in the hippocampus of these mice (PS19-E4, n = 29; PS19-E3, n = 22; PS19-E4-S/S, n = 20; PS19-E4-R/S, n = 22) and WT mice (n = 11). d, Representative western blot images with p-Tau-specific AT8 or PHF1 antibody. TUJ1 was used as a loading control. e,f, Quantification of AT8+ (e) and PHF1+ (f) p-Tau levels in hippocampal lysates of PS19-E4 (n = 8), PS19-E3 (n = 7), PS19-E4-S/S (n = 7) and PS19-E4-R/S (n = 8) mice. p-Tau levels were normalized to TUJ1 first and then to those of PS19-E4 mice. g, Representative images of four AT8 staining patterns in the hippocampus. h, Distribution of four p-Tau staining patterns in the hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice (PS19-E4, n = 29; PS19-E3, n = 22; PS19-E4-S/S, n = 20; PS19-E4-R/S, n = 22). Scale bars in b and g, 500 µm. Throughout, data are expressed as mean ± s.e.m. Differences between groups were determined by Welch’s ANOVA followed by Dunnett’s T3 multiple comparison test (c) or ordinary one-way ANOVA followed by Dunnett’s multiple comparison test (e,f). Comparisons of P ≤ 0.05 are labeled on the graph. Source data
Fig. 1 Homozygous R136S mutation rescues APOE4-promoted Tau pathology in tauopathy mice. a, Schematic of CRISPR–Cas-9-mediated gene editing strategy to generate human APOE4-R136S knock-in mice. b, Representative images of p-Tau immunostaining in hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice with the AT8 monoclonal antibody. c, Quantification of the percent AT8 coverage area in the hippocampus of these mice (PS19-E4, n = 29; PS19-E3, n = 22; PS19-E4-S/S, n = 20; PS19-E4-R/S, n = 22) and WT mice (n = 11). d, Representative western blot images with p-Tau-specific AT8 or PHF1 antibody. TUJ1 was used as a loading control. e,f, Quantification of AT8+ (e) and PHF1+ (f) p-Tau levels in hippocampal lysates of PS19-E4 (n = 8), PS19-E3 (n = 7), PS19-E4-S/S (n = 7) and PS19-E4-R/S (n = 8) mice. p-Tau levels were normalized to TUJ1 first and then to those of PS19-E4 mice. g, Representative images of four AT8 staining patterns in the hippocampus. h, Distribution of four p-Tau staining patterns in the hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice (PS19-E4, n = 29; PS19-E3, n = 22; PS19-E4-S/S, n = 20; PS19-E4-R/S, n = 22). Scale bars in b and g, 500 µm. Throughout, data are expressed as mean ± s.e.m. Differences between groups were determined by Welch’s ANOVA followed by Dunnett’s T3 multiple comparison test (c) or ordinary one-way ANOVA followed by Dunnett’s multiple comparison test (e,f). Comparisons of P ≤ 0.05 are labeled on the graph. Source dataCC BY · PubMed Central
Fig. 2 Homozygous R136S mutation protects against APOE4-induced p-Tau accumulation in human neurons. a–d, Representative western blot images (a) and quantification of APOE (b), PHF1+ p-Tau (c) and AT8+ p-Tau (d) levels in lysates of E4, E3, E4-S/S or E4-R/S neurons. In b, APOE levels were normalized to those of E4. TUJ1 was used as loading control (E4, n = 32; E3, n = 31; E4-S/S, n = 64 (n = 24 from E4-S/S-A; n = 8 from E4-S/S-B; n = 32 from E4-S/S-C); E4-R/S, n = 59 (n = 31 from E4-R/S-A; n = 28 from E4-R/S-B)). In c, PHF1+ p-Tau levels were normalized to those of E4. TUJ1 was used as loading control (E4, n = 32; E3, n = 31; E4-S/S, n = 64 (n = 24 from E4-S/S-A; n = 8 from E4-S/S-B; n = 32 from E4-S/S-C); E4-R/S, n = 59 (n = 31 from E4-R/S-A; n = 28 from E4-R/S-B)). In d, AT8+ p-Tau levels were normalized to those of E4. TUJ1 was used as loading control (E4, n = 28; E3, n = 27; E4-S/S, n = 55 (n = 16 from E4-S/S-A; n = 11 from E4-S/S-B; n = 28 from E4-S/S-C); E4-R/S, n = 59 (n = 31 from E4-R/S-A; n = 28 from E4-R/S-B)). e, Representative images showing immunostaining of p-Tau (PHF1), total Tau and MAP2 in E4, E3, E4-S/S or E4-R/S human neurons, with some PHF1+ puncta (white arrowheads in insets). f,g, Quantification of fraction of the PHF1+ area over MAP2+ area (f) and fraction of PHF1+ puncta area over MAP2+ area (g) (E4, n = 25 fields of view; E3, n = 25 fields of view; E4-S/S, n = 38 fields of view; E4-R/S, n = 33 fields of view). The ratio of PHF1+ area (f) and PHF1+ puncta area (g) over MAP2 area was normalized to that of E4. Western blot data (b–d) were made up of at least three independent rounds of differentiation, and all data were combined. Scale bars (e), 20 µm. In b–d, n = biological replicates. Throughout, data are expressed as mean ± s.e.m. Differences between groups were determined by Welch’s ANOVA followed by Dunnett’s T3 multiple comparison test (b–d,g) or ordinary one-way ANOVA followed by Tukey’s multiple comparison test (f). Comparisons of P ≤ 0.05 are labeled on the graph. Source data
Fig. 2 Homozygous R136S mutation protects against APOE4-induced p-Tau accumulation in human neurons. a–d, Representative western blot images (a) and quantification of APOE (b), PHF1+ p-Tau (c) and AT8+ p-Tau (d) levels in lysates of E4, E3, E4-S/S or E4-R/S neurons. In b, APOE levels were normalized to those of E4. TUJ1 was used as loading control (E4, n = 32; E3, n = 31; E4-S/S, n = 64 (n = 24 from E4-S/S-A; n = 8 from E4-S/S-B; n = 32 from E4-S/S-C); E4-R/S, n = 59 (n = 31 from E4-R/S-A; n = 28 from E4-R/S-B)). In c, PHF1+ p-Tau levels were normalized to those of E4. TUJ1 was used as loading control (E4, n = 32; E3, n = 31; E4-S/S, n = 64 (n = 24 from E4-S/S-A; n = 8 from E4-S/S-B; n = 32 from E4-S/S-C); E4-R/S, n = 59 (n = 31 from E4-R/S-A; n = 28 from E4-R/S-B)). In d, AT8+ p-Tau levels were normalized to those of E4. TUJ1 was used as loading control (E4, n = 28; E3, n = 27; E4-S/S, n = 55 (n = 16 from E4-S/S-A; n = 11 from E4-S/S-B; n = 28 from E4-S/S-C); E4-R/S, n = 59 (n = 31 from E4-R/S-A; n = 28 from E4-R/S-B)). e, Representative images showing immunostaining of p-Tau (PHF1), total Tau and MAP2 in E4, E3, E4-S/S or E4-R/S human neurons, with some PHF1+ puncta (white arrowheads in insets). f,g, Quantification of fraction of the PHF1+ area over MAP2+ area (f) and fraction of PHF1+ puncta area over MAP2+ area (g) (E4, n = 25 fields of view; E3, n = 25 fields of view; E4-S/S, n = 38 fields of view; E4-R/S, n = 33 fields of view). The ratio of PHF1+ area (f) and PHF1+ puncta area (g) over MAP2 area was normalized to that of E4. Western blot data (b–d) were made up of at least three independent rounds of differentiation, and all data were combined. Scale bars (e), 20 µm. In b–d, n = biological replicates. Throughout, data are expressed as mean ± s.e.m. Differences between groups were determined by Welch’s ANOVA followed by Dunnett’s T3 multiple comparison test (b–d,g) or ordinary one-way ANOVA followed by Tukey’s multiple comparison test (f). Comparisons of P ≤ 0.05 are labeled on the graph. Source dataCC BY · PubMed Central
Fig. 3 Homozygous R136S mutation protects against APOE4-induced p-Tau accumulation by reducing Tau uptake via the HSPG pathway. a, Diagram of Tau-488 uptake assay. Neurons treated with either Tau-488 alone (left) or Tau-488 together with 100 µg ml−1 heparin (right) before flow cytometry analysis. b, Measurement of individual neuronal Tau-488 uptake (25 nM, 1-h incubation) based on MFI per cell in human neurons. c, Measurement of Tau-488 uptake (25 nM, 1-h incubation) based on percent Tau-488+ human neurons. In b,c, n = independent experiments and normalized to E4 MFI (b) or uptake (%) (c). E4, n = 13; E4+heparin, n = 4; E3, n = 4; E3+heparin, n = 4; E4-S/S, n = 4; E4-S/S+heparin, n = 4; E4-R/S, n = 4; E4-R/S+heparin, n = 4; EKO, n = 3; EKO+heparin, n = 4. Analysis was performed on a live cell population of estimated 5,000 cells for each sample. d, Experimental design for long-term heparin treatment of human neurons. e,f, Representative western blot images (e) and quantification of PHF1+ p-Tau levels (f) in lysates of E4, E3, E4-S/S or E4-R/S neurons under long-term heparin treatment. In f, PHF1+ p-Tau levels were normalized to those of E4. TUJ1 was used as loading control (E4, n = 16; E4+heparin, n = 16; E3, n = 16; E3+heparin, n = 15; E4-S/S, n = 32; E4-S/S+heparin, n = 32; E4-R/S, n = 16; E4-R/S+heparin, n = 16). g, Experimental design for E4 neuron-conditioned medium treatment of human neurons with different APOE genotypes. h,i, Representative western blot images (h) and quantification of PHF1+ p-Tau levels (i) in lysates of E4, E3, E4-S/S or E4-R/S neurons after E4 neuron-conditioned medium treatment. In i, PHF1+ p-Tau levels were normalized to those of E4. TUJ1 was used as loading control (E4, n = 5; E3, n = 5; E4-S/S, n = 10; E4-R/S, n = 10). In f,i, n =biological replicates. Throughout, data are expressed as mean ± s.e.m. Differences between groups were determined by two-way ANOVA followed by Tukey’s multiple comparison test (b,c,f) or ordinary one-way ANOVA followed by Tukey’s multiple comparison test (i). Some comparisons of P ≤ 0.05 are labeled on the graph. Hep, heparin; fluor, fluorescence. Source data
Fig. 3 Homozygous R136S mutation protects against APOE4-induced p-Tau accumulation by reducing Tau uptake via the HSPG pathway. a, Diagram of Tau-488 uptake assay. Neurons treated with either Tau-488 alone (left) or Tau-488 together with 100 µg ml−1 heparin (right) before flow cytometry analysis. b, Measurement of individual neuronal Tau-488 uptake (25 nM, 1-h incubation) based on MFI per cell in human neurons. c, Measurement of Tau-488 uptake (25 nM, 1-h incubation) based on percent Tau-488+ human neurons. In b,c, n = independent experiments and normalized to E4 MFI (b) or uptake (%) (c). E4, n = 13; E4+heparin, n = 4; E3, n = 4; E3+heparin, n = 4; E4-S/S, n = 4; E4-S/S+heparin, n = 4; E4-R/S, n = 4; E4-R/S+heparin, n = 4; EKO, n = 3; EKO+heparin, n = 4. Analysis was performed on a live cell population of estimated 5,000 cells for each sample. d, Experimental design for long-term heparin treatment of human neurons. e,f, Representative western blot images (e) and quantification of PHF1+ p-Tau levels (f) in lysates of E4, E3, E4-S/S or E4-R/S neurons under long-term heparin treatment. In f, PHF1+ p-Tau levels were normalized to those of E4. TUJ1 was used as loading control (E4, n = 16; E4+heparin, n = 16; E3, n = 16; E3+heparin, n = 15; E4-S/S, n = 32; E4-S/S+heparin, n = 32; E4-R/S, n = 16; E4-R/S+heparin, n = 16). g, Experimental design for E4 neuron-conditioned medium treatment of human neurons with different APOE genotypes. h,i, Representative western blot images (h) and quantification of PHF1+ p-Tau levels (i) in lysates of E4, E3, E4-S/S or E4-R/S neurons after E4 neuron-conditioned medium treatment. In i, PHF1+ p-Tau levels were normalized to those of E4. TUJ1 was used as loading control (E4, n = 5; E3, n = 5; E4-S/S, n = 10; E4-R/S, n = 10). In f,i, n =biological replicates. Throughout, data are expressed as mean ± s.e.m. Differences between groups were determined by two-way ANOVA followed by Tukey’s multiple comparison test (b,c,f) or ordinary one-way ANOVA followed by Tukey’s multiple comparison test (i). Some comparisons of P ≤ 0.05 are labeled on the graph. Hep, heparin; fluor, fluorescence. Source dataCC BY · PubMed Central
Fig. 4 The R136S mutation ameliorates APOE4-driven neurodegeneration in tauopathy mice. a, Representative images of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mouse brain sections stained with Sudan black to enhance hippocampal visualization (scale bar, 1 mm). b,c, Quantification of hippocampal volume (WT, n = 11; PS19-E4, n = 31; PS19-E3, n = 23; PS19-E4-S/S, n = 20; PS19-E4-R/S, n = 23; n = mice) (b) and posterior lateral ventricle volume (PS19-E4, n = 30; PS19-E3, n = 23; PS19-E4-S/S, n = 20; PS19-E4-R/S, n = 23; n = mice) (c). d,e, Correlation of the APOE4-R136S gene copy number with the average of hippocampal volume (PS19-E4, n = 31; PS19-E4-S/S, n = 20; PS19-E4-R/S, n = 23; n = mice) (d) or the average of posterior lateral ventricle volume (PS19-E4, n = 30; PS19-E4-S/S, n = 20; PS19-E4-R/S, n = 23; n = mice) (e) in PS19-E4 mice with 0, 1 or 2 copies of the APOE4-R136S gene mutation. f–i, Correlations between percent AT8 coverage area and hippocampal volume in 10-month-old PS19-E4 (n = 29) (f), PS19-E3 (n = 21) (g), PS19-E4-S/S (n = 19) (h) and PS19-E4-R/S (n = 22) (i) mice. j, Representative images of the DG stained for NeuN in 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice (scale bar, 100 μm). k, Quantification of DG GC layer thickness in 10-month-old PS19-E4 (n = 30), PS19-E3 (n = 23), PS19-E4-S/S (n = 20) and PS19-E4-R/S (n = 23) mice. Throughout, data are expressed as mean ± s.e.m. except for correlation plots. Differences between groups were determined by ordinary one-way ANOVA followed by Tukey’s multiple comparison test (b,k) or Welch’s ANOVA followed by Dunnett’s T3 multiple comparison test (c). Comparisons of P ≤ 0.05 are labeled on the graph. Pearson’s correlation analysis (two-sided). LV, lateral ventricle. Source data
Fig. 4 The R136S mutation ameliorates APOE4-driven neurodegeneration in tauopathy mice. a, Representative images of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mouse brain sections stained with Sudan black to enhance hippocampal visualization (scale bar, 1 mm). b,c, Quantification of hippocampal volume (WT, n = 11; PS19-E4, n = 31; PS19-E3, n = 23; PS19-E4-S/S, n = 20; PS19-E4-R/S, n = 23; n = mice) (b) and posterior lateral ventricle volume (PS19-E4, n = 30; PS19-E3, n = 23; PS19-E4-S/S, n = 20; PS19-E4-R/S, n = 23; n = mice) (c). d,e, Correlation of the APOE4-R136S gene copy number with the average of hippocampal volume (PS19-E4, n = 31; PS19-E4-S/S, n = 20; PS19-E4-R/S, n = 23; n = mice) (d) or the average of posterior lateral ventricle volume (PS19-E4, n = 30; PS19-E4-S/S, n = 20; PS19-E4-R/S, n = 23; n = mice) (e) in PS19-E4 mice with 0, 1 or 2 copies of the APOE4-R136S gene mutation. f–i, Correlations between percent AT8 coverage area and hippocampal volume in 10-month-old PS19-E4 (n = 29) (f), PS19-E3 (n = 21) (g), PS19-E4-S/S (n = 19) (h) and PS19-E4-R/S (n = 22) (i) mice. j, Representative images of the DG stained for NeuN in 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice (scale bar, 100 μm). k, Quantification of DG GC layer thickness in 10-month-old PS19-E4 (n = 30), PS19-E3 (n = 23), PS19-E4-S/S (n = 20) and PS19-E4-R/S (n = 23) mice. Throughout, data are expressed as mean ± s.e.m. except for correlation plots. Differences between groups were determined by ordinary one-way ANOVA followed by Tukey’s multiple comparison test (b,k) or Welch’s ANOVA followed by Dunnett’s T3 multiple comparison test (c). Comparisons of P ≤ 0.05 are labeled on the graph. Pearson’s correlation analysis (two-sided). LV, lateral ventricle. Source dataCC BY · PubMed Central
Fig. 5 The R136S mutation reduces APOE4-driven gliosis in tauopathy mice. a,b, Representative images of GFAP immunostaining of astrocytes in the hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice (a) and quantification of the number of GFAP+ cells per mm2 (b) in the hippocampus of these mice and WT mice. c,d, Representative images of S100β immunostaining of astrocytes in the hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice (c) and quantification of the number of S100β+ cells per mm2 (d) in the hippocampus of these mice and WT mice. e–h, Correlations between GFAP+ cells per mm2 and hippocampal volume in PS19-E4 (n = 31) (e), PS19-E3 (n = 23) (f), PS19-E4-S/S (n = 20) (g) and PS19-E4-R/S (n = 23) (h) mice. i–l, Correlations between S100β+ cells per mm2 and hippocampal volume in PS19-E4 (n = 31) (i), PS19-E3 (n = 23) (j), PS19-E4-S/S (n = 20) (k) and PS19-E4-R/S (n = 23) (l) mice. m,n, Representative images of Iba1 immunostaining of microglia in the hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice (m) and quantification of the number of Iba1+ cells per mm2 (n) in the hippocampus of these mice and WT mice. o,p, Representative images of CD68 immunostaining of microglia in the hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice (o) and quantification of the number of CD68+ cells per mm2 (p) in the hippocampus of these mice and WT mice. q–t, Correlations between Iba1+ cells per mm2 and hippocampal volume in PS19-E4 (n = 31) (q), PS19-E3 (n = 23) (r), PS19-E4-S/S (n = 20) (s) and PS19-E4-R/S (n = 23) (t) mice. u–x, Correlations between CD68+ cells per mm2 and hippocampal volume in PS19-E4 (n = 31) (u), PS19-E3 (n = 23) (v), PS19-E4-S/S (n = 20) (w) and PS19-E4-R/S (n = 23) (x) mice. For all quantifications in b,d,n,p, WT, n = 11; PS19-E4, n = 31; PS19-E3, n = 24; PS19-E4-S/S, n = 21; PS19-E4-R/S, n = 23. Scale bars, 500 µm in a,c,m,o. Throughout, data are expressed as mean ± s.e.m. In b,d,n,p, differences between groups were determined by Welch’s ANOVA followed by Dunnett’s T3 multiple comparison test. Comparisons of P ≤ 0.05 are labeled on the graph. Pearson’s correlation analysis (two-sided). Source data
Fig. 5 The R136S mutation reduces APOE4-driven gliosis in tauopathy mice. a,b, Representative images of GFAP immunostaining of astrocytes in the hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice (a) and quantification of the number of GFAP+ cells per mm2 (b) in the hippocampus of these mice and WT mice. c,d, Representative images of S100β immunostaining of astrocytes in the hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice (c) and quantification of the number of S100β+ cells per mm2 (d) in the hippocampus of these mice and WT mice. e–h, Correlations between GFAP+ cells per mm2 and hippocampal volume in PS19-E4 (n = 31) (e), PS19-E3 (n = 23) (f), PS19-E4-S/S (n = 20) (g) and PS19-E4-R/S (n = 23) (h) mice. i–l, Correlations between S100β+ cells per mm2 and hippocampal volume in PS19-E4 (n = 31) (i), PS19-E3 (n = 23) (j), PS19-E4-S/S (n = 20) (k) and PS19-E4-R/S (n = 23) (l) mice. m,n, Representative images of Iba1 immunostaining of microglia in the hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice (m) and quantification of the number of Iba1+ cells per mm2 (n) in the hippocampus of these mice and WT mice. o,p, Representative images of CD68 immunostaining of microglia in the hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S and PS19-E4-R/S mice (o) and quantification of the number of CD68+ cells per mm2 (p) in the hippocampus of these mice and WT mice. q–t, Correlations between Iba1+ cells per mm2 and hippocampal volume in PS19-E4 (n = 31) (q), PS19-E3 (n = 23) (r), PS19-E4-S/S (n = 20) (s) and PS19-E4-R/S (n = 23) (t) mice. u–x, Correlations between CD68+ cells per mm2 and hippocampal volume in PS19-E4 (n = 31) (u), PS19-E3 (n = 23) (v), PS19-E4-S/S (n = 20) (w) and PS19-E4-R/S (n = 23) (x) mice. For all quantifications in b,d,n,p, WT, n = 11; PS19-E4, n = 31; PS19-E3, n = 24; PS19-E4-S/S, n = 21; PS19-E4-R/S, n = 23. Scale bars, 500 µm in a,c,m,o. Throughout, data are expressed as mean ± s.e.m. In b,d,n,p, differences between groups were determined by Welch’s ANOVA followed by Dunnett’s T3 multiple comparison test. Comparisons of P ≤ 0.05 are labeled on the graph. Pearson’s correlation analysis (two-sided). Source dataCC BY · PubMed Central
Fig. 6 snRNA-seq reveals protective effects of the R136S mutation on APOE4-driven neuronal and oligodendrocytic deficits in mice. a, UMAP projection of 38 distinct cell clusters in hippocampi of 10-month-old PS19-E4 (n = 4), PS19-E3 (n = 3), PS19-E4-S/S (n = 4) and PS19-E4-R/S (n = 4) mice. b, Feature plot showing relative levels of normalized human APOE gene expression across all 38 hippocampal cell clusters by APOE genotype (PS19-E4, n = 4; PS19-E3, n = 3; PS19-E4-S/S, n = 4; PS19-E4-R/S, n = 4; n = mice). c, UMAP projection highlighting hippocampal cell clusters 1, 6, 7, 9 and 28 for each genotype group. d, Box plot of the proportion of cells from each sample in clusters 1, 6, 7, 9 and 28 in PS19-E4 (n = 4), PS19-E3 (n = 3), PS19-E4-R/S (n = 4), and PS19-E4-S/S (n = 4) mice. The lower, middle and upper hinges of the box plots correspond to the 25th, 50th and 75th percentiles, respectively. The upper whisker of the box plot extends from the upper hinge to the largest value no further than 1.5 Å~ IQR from the upper hinge. IQR, interquartile range, or distance between the 25th and 75th percentiles. The lower whisker extends from the lower hinge to the smallest value at most 1.5 Å~ IQR from the lower hinge. The LORs are the mean ± s.e.m. estimates of LOR for these clusters, which represents the change in the log odds of cells per sample from PS19-E3, PS19-E4-R/S or PS19-E4-S/S mice belonging to the respective clusters compared to the log odds of cells per sample from PS19-E4 mice. e,f, KEGG pathway enrichment dot plot of top 20 pathways significantly enriched for DE genes of neuronal cluster 6 in PS19-E4-S/S (e) or PS19-E4-R/S (f) versus PS19-E4 mice. P values are based on a two-sided hypergeometric test and are adjusted for multiple testing using the Benjamini–Hochberg method. Gene ratio represents the proportion of genes in the respective gene set that are deemed to be DE using the two-sided Wilcoxon rank-sum test as implemented in the FindMarkers function in Seurat. g, KEGG pathway enrichment dot plot of top 20 pathways significantly enriched for DE genes of oligodendrocyte cluster 9 versus oligodendrocyte cluster 2. h, Heat map plot of LOR per unit change in each pathological measurement for clusters 1, 6, 7, 9 and 28. The LOR represents the mean estimate of the change in the log odds of cells per sample from a given animal model, corresponding to a unit change in a given histopathological parameter. Associations with pathologies are colored (negative associations, blue; positive associations, red). P values in d are from fits to a GLMM_AM, and P values in h are from fits to a GLMM_histopathology; the associated tests are two-sided. All error bars represent s.e.m. Ex neuron, excitatory neuron; In neuron, inhibitory neuron.
Fig. 6 snRNA-seq reveals protective effects of the R136S mutation on APOE4-driven neuronal and oligodendrocytic deficits in mice. a, UMAP projection of 38 distinct cell clusters in hippocampi of 10-month-old PS19-E4 (n = 4), PS19-E3 (n = 3), PS19-E4-S/S (n = 4) and PS19-E4-R/S (n = 4) mice. b, Feature plot showing relative levels of normalized human APOE gene expression across all 38 hippocampal cell clusters by APOE genotype (PS19-E4, n = 4; PS19-E3, n = 3; PS19-E4-S/S, n = 4; PS19-E4-R/S, n = 4; n = mice). c, UMAP projection highlighting hippocampal cell clusters 1, 6, 7, 9 and 28 for each genotype group. d, Box plot of the proportion of cells from each sample in clusters 1, 6, 7, 9 and 28 in PS19-E4 (n = 4), PS19-E3 (n = 3), PS19-E4-R/S (n = 4), and PS19-E4-S/S (n = 4) mice. The lower, middle and upper hinges of the box plots correspond to the 25th, 50th and 75th percentiles, respectively. The upper whisker of the box plot extends from the upper hinge to the largest value no further than 1.5 Å~ IQR from the upper hinge. IQR, interquartile range, or distance between the 25th and 75th percentiles. The lower whisker extends from the lower hinge to the smallest value at most 1.5 Å~ IQR from the lower hinge. The LORs are the mean ± s.e.m. estimates of LOR for these clusters, which represents the change in the log odds of cells per sample from PS19-E3, PS19-E4-R/S or PS19-E4-S/S mice belonging to the respective clusters compared to the log odds of cells per sample from PS19-E4 mice. e,f, KEGG pathway enrichment dot plot of top 20 pathways significantly enriched for DE genes of neuronal cluster 6 in PS19-E4-S/S (e) or PS19-E4-R/S (f) versus PS19-E4 mice. P values are based on a two-sided hypergeometric test and are adjusted for multiple testing using the Benjamini–Hochberg method. Gene ratio represents the proportion of genes in the respective gene set that are deemed to be DE using the two-sided Wilcoxon rank-sum test as implemented in the FindMarkers function in Seurat. g, KEGG pathway enrichment dot plot of top 20 pathways significantly enriched for DE genes of oligodendrocyte cluster 9 versus oligodendrocyte cluster 2. h, Heat map plot of LOR per unit change in each pathological measurement for clusters 1, 6, 7, 9 and 28. The LOR represents the mean estimate of the change in the log odds of cells per sample from a given animal model, corresponding to a unit change in a given histopathological parameter. Associations with pathologies are colored (negative associations, blue; positive associations, red). P values in d are from fits to a GLMM_AM, and P values in h are from fits to a GLMM_histopathology; the associated tests are two-sided. All error bars represent s.e.m. Ex neuron, excitatory neuron; In neuron, inhibitory neuron.CC BY · PubMed Central
Fig. 7 The APOE4-R136S mutation increases disease-protective and decreases disease-associated astrocyte subpopulations. a, UMAP projection of 12 astrocyte subclusters after subclustering hippocampal cell clusters 13 and 36 (Fig. 6a) from 10-month-old mice with different APOE genotypes. b, UMAP projection highlighting astrocyte subclusters 3, 5 and 7 for each genotype group (PS19-E4, n = 4; PS19-E3, n = 3; PS19-E4-S/S, n = 4; PS19-E4-R/S, n = 4; n = mice). c, Box plot of the proportion of cells from each sample in astrocyte subclusters 3, 5 and 7 in PS19-E4 (n = 4), PS19-E3 (n = 3), PS19-E4-R/S (n = 4), and PS19-E4-S/S (n = 4) mice. The lower, middle and upper hinges of the box plots correspond to the 25th, 50th and 75th percentiles, respectively (see Fig. 6d for details). The LORs are the mean ± s.e.m. estimates of LOR for these clusters, which represents the change in the log odds of cells per sample from PS19-E3, PS19-E4-R/S or PS19-E4-S/S mice belonging to the respective clusters compared to the log odds of cells per sample from PS19-E4 mice. LOR versus PS19-E4 for subcluster 3: PS19-E3, 0.67 ± 0.30; PS19-E4-S/S, 0.71 ± 0.27; subcluster 5: PS19-E4-S/S, −0.74 ± 0.34; subcluster 7: PS19-E4-R/S, −1.57 ± 0.57; PS19-E4-S/S, −2.75 ± 0.62. d, Dot plot of normalized average expression of selected homeostatic and DAA marker genes for astrocyte subclusters 3, 5 and 7. e, Volcano plot for top 30 DE genes of astrocyte subcluster 5 versus other astrocyte subclusters. f, Volcano plot for top 30 DE genes of astrocyte subcluster 5 in PS19-E4-S/S versus PS19-E4 mice. g, Volcano plot for top 30 DE genes of astrocyte subcluster 7 versus other astrocyte subclusters. h, Volcano plot for top 30 DE genes of astrocyte subcluster 7 in PS19-E4-S/S versus PS19-E4 mice. i, Representative images of Nkain2+GFAP+ astrocytes in the hippocampus of 10-month-old PS19-E4 (n = 9) and PS19-E4-S/S (n = 10) mice. j, Quantification of the number of Nkain2+GFAP+ cells (per mm2) within the molecular layer of hippocampus. k, Representative images of Id3+GFAP+ astrocytes in the hippocampus of 10-month-old PS19-E4 (n = 10) and PS19-E4-S/S (n = 10) mice. l, Quantification of the number of Id3+GFAP+ cells (per mm2) within the molecular layer of hippocampus. m, Heat map plot of LOR per unit change in each pathological measurement for astrocyte subclusters 3, 5 and 7. P values in c are from fits to a GLMM_AM, and P values in m are from fits to a GLMM_histopathology; the associated tests are two-sided. In e–h, horizonal dashed line indicates P = 0.05, and vertical dashed lines indicate log2 fold change = 0.4. The unadjusted P values and log2 fold change values used were generated from the gene set enrichment analysis using the two-sided Wilcoxon rank-sum test as implemented in the FindMarkers function of the Seurat package. Gene names highlighted in red text indicate that they are selected marker genes for DAAs. Scale bars in i and k, 50 µm. All error bars represent s.e.m. Differences between groups in j and l were determined by unpaired, two-sided Welch’s t-test. AS, astrocyte; NS, not significant; FC, fold change. Source data
Fig. 7 The APOE4-R136S mutation increases disease-protective and decreases disease-associated astrocyte subpopulations. a, UMAP projection of 12 astrocyte subclusters after subclustering hippocampal cell clusters 13 and 36 (Fig. 6a) from 10-month-old mice with different APOE genotypes. b, UMAP projection highlighting astrocyte subclusters 3, 5 and 7 for each genotype group (PS19-E4, n = 4; PS19-E3, n = 3; PS19-E4-S/S, n = 4; PS19-E4-R/S, n = 4; n = mice). c, Box plot of the proportion of cells from each sample in astrocyte subclusters 3, 5 and 7 in PS19-E4 (n = 4), PS19-E3 (n = 3), PS19-E4-R/S (n = 4), and PS19-E4-S/S (n = 4) mice. The lower, middle and upper hinges of the box plots correspond to the 25th, 50th and 75th percentiles, respectively (see Fig. 6d for details). The LORs are the mean ± s.e.m. estimates of LOR for these clusters, which represents the change in the log odds of cells per sample from PS19-E3, PS19-E4-R/S or PS19-E4-S/S mice belonging to the respective clusters compared to the log odds of cells per sample from PS19-E4 mice. LOR versus PS19-E4 for subcluster 3: PS19-E3, 0.67 ± 0.30; PS19-E4-S/S, 0.71 ± 0.27; subcluster 5: PS19-E4-S/S, −0.74 ± 0.34; subcluster 7: PS19-E4-R/S, −1.57 ± 0.57; PS19-E4-S/S, −2.75 ± 0.62. d, Dot plot of normalized average expression of selected homeostatic and DAA marker genes for astrocyte subclusters 3, 5 and 7. e, Volcano plot for top 30 DE genes of astrocyte subcluster 5 versus other astrocyte subclusters. f, Volcano plot for top 30 DE genes of astrocyte subcluster 5 in PS19-E4-S/S versus PS19-E4 mice. g, Volcano plot for top 30 DE genes of astrocyte subcluster 7 versus other astrocyte subclusters. h, Volcano plot for top 30 DE genes of astrocyte subcluster 7 in PS19-E4-S/S versus PS19-E4 mice. i, Representative images of Nkain2+GFAP+ astrocytes in the hippocampus of 10-month-old PS19-E4 (n = 9) and PS19-E4-S/S (n = 10) mice. j, Quantification of the number of Nkain2+GFAP+ cells (per mm2) within the molecular layer of hippocampus. k, Representative images of Id3+GFAP+ astrocytes in the hippocampus of 10-month-old PS19-E4 (n = 10) and PS19-E4-S/S (n = 10) mice. l, Quantification of the number of Id3+GFAP+ cells (per mm2) within the molecular layer of hippocampus. m, Heat map plot of LOR per unit change in each pathological measurement for astrocyte subclusters 3, 5 and 7. P values in c are from fits to a GLMM_AM, and P values in m are from fits to a GLMM_histopathology; the associated tests are two-sided. In e–h, horizonal dashed line indicates P = 0.05, and vertical dashed lines indicate log2 fold change = 0.4. The unadjusted P values and log2 fold change values used were generated from the gene set enrichment analysis using the two-sided Wilcoxon rank-sum test as implemented in the FindMarkers function of the Seurat package. Gene names highlighted in red text indicate that they are selected marker genes for DAAs. Scale bars in i and k, 50 µm. All error bars represent s.e.m. Differences between groups in j and l were determined by unpaired, two-sided Welch’s t-test. AS, astrocyte; NS, not significant; FC, fold change. Source dataCC BY · PubMed Central
Fig. 8 The APOE4-R136S mutation increases disease-protective and decreases disease-associated microglial subpopulations. a, UMAP projection of 15 microglia subclusters after subclustering hippocampal cell clusters 17 and 19 (Fig. 6a) from 10-month-old mice with different APOE genotypes. b, UMAP projection highlighting microglia subclusters 2, 8 and 11 for each mouse genotype group (PS19-E4, n = 4; PS19-E3, n = 3; PS19-E4-S/S, n = 4; PS19-E4-R/S, n = 4; n = mice). c, Box plot of the proportion of cells from each sample in microglia subclusters 2, 8, and 11 in PS19-E4 (n = 4), PS19-E3 (n = 3), PS19-E4-R/S (n = 4), and PS19-E4-S/S (n = 4) mice. The lower, middle and upper hinges of the box plots correspond to the 25th, 50th and 75th percentiles, respectively (see Fig. 6d for details). The LORs are the mean ± s.e.m. estimates of LOR for these clusters, which represents the change in the log odds of cells per sample from PS19-E3, PS19-E4-R/S or PS19-E4-S/S mice belonging to the respective clusters compared to the log odds of cells per sample from PS19-E4 mice. LOR versus PS19-E4 for subcluster 2: PS19-E4-S/S, 3.02 ± 1.28; subcluster 8: PS19-E4-R/S, −2.58 ± 1.02; PS19-E4-S/S, −3.36 ± 1.17; subcluster 11: 2.69 ± 1.33. d, Dot plot of normalized average expression of selected homeostatic and DAM marker genes for microglia subclusters 2, 8 and 11. e, Volcano plot for top 30 DE genes of microglia subcluster 8 versus other microglia subclusters. f, Volcano plot for top 30 DE genes of microglia subcluster 8 in PS19-E4-S/S versus PS19-E4 mice. g, Representative images of Gpnmb+Iba1+ microglia in the hippocampus of 10-month-old PS19-E4 (n = 10) and PS19-E4-S/S (n = 10) mice. Scale bars, 50 µm. h, Quantification of the number of Gpnmb+Iba1+ cells (per mm2) within the DG of hippocampus. Difference between groups in h was determined by unpaired, two-sided Welch’s t-test. i, Heat map plot of LOR per unit change in each pathological measurement for microglia subclusters 2, 8 and 11. P values in c are from fits to a GLMM_AM, and P values in i are from fits to a GLMM_histopathology; the associated tests are two-sided. In e,f, horizonal dashed line indicates P = 0.05, and vertical dashed lines indicate log2 fold change = 0.4. The unadjusted P values and log2 fold change values used were generated from the gene set enrichment analysis using the two-sided Wilcoxon rank-sum test as implemented in the FindMarkers function of the Seurat package. Gene names highlighted in red text indicate that they are selected marker genes for DAMs. All error bars represent s.e.m. AS, astrocyte; FC, fold change; MG, microglia; NS, not significant. Source data
Fig. 8 The APOE4-R136S mutation increases disease-protective and decreases disease-associated microglial subpopulations. a, UMAP projection of 15 microglia subclusters after subclustering hippocampal cell clusters 17 and 19 (Fig. 6a) from 10-month-old mice with different APOE genotypes. b, UMAP projection highlighting microglia subclusters 2, 8 and 11 for each mouse genotype group (PS19-E4, n = 4; PS19-E3, n = 3; PS19-E4-S/S, n = 4; PS19-E4-R/S, n = 4; n = mice). c, Box plot of the proportion of cells from each sample in microglia subclusters 2, 8, and 11 in PS19-E4 (n = 4), PS19-E3 (n = 3), PS19-E4-R/S (n = 4), and PS19-E4-S/S (n = 4) mice. The lower, middle and upper hinges of the box plots correspond to the 25th, 50th and 75th percentiles, respectively (see Fig. 6d for details). The LORs are the mean ± s.e.m. estimates of LOR for these clusters, which represents the change in the log odds of cells per sample from PS19-E3, PS19-E4-R/S or PS19-E4-S/S mice belonging to the respective clusters compared to the log odds of cells per sample from PS19-E4 mice. LOR versus PS19-E4 for subcluster 2: PS19-E4-S/S, 3.02 ± 1.28; subcluster 8: PS19-E4-R/S, −2.58 ± 1.02; PS19-E4-S/S, −3.36 ± 1.17; subcluster 11: 2.69 ± 1.33. d, Dot plot of normalized average expression of selected homeostatic and DAM marker genes for microglia subclusters 2, 8 and 11. e, Volcano plot for top 30 DE genes of microglia subcluster 8 versus other microglia subclusters. f, Volcano plot for top 30 DE genes of microglia subcluster 8 in PS19-E4-S/S versus PS19-E4 mice. g, Representative images of Gpnmb+Iba1+ microglia in the hippocampus of 10-month-old PS19-E4 (n = 10) and PS19-E4-S/S (n = 10) mice. Scale bars, 50 µm. h, Quantification of the number of Gpnmb+Iba1+ cells (per mm2) within the DG of hippocampus. Difference between groups in h was determined by unpaired, two-sided Welch’s t-test. i, Heat map plot of LOR per unit change in each pathological measurement for microglia subclusters 2, 8 and 11. P values in c are from fits to a GLMM_AM, and P values in i are from fits to a GLMM_histopathology; the associated tests are two-sided. In e,f, horizonal dashed line indicates P = 0.05, and vertical dashed lines indicate log2 fold change = 0.4. The unadjusted P values and log2 fold change values used were generated from the gene set enrichment analysis using the two-sided Wilcoxon rank-sum test as implemented in the FindMarkers function of the Seurat package. Gene names highlighted in red text indicate that they are selected marker genes for DAMs. All error bars represent s.e.m. AS, astrocyte; FC, fold change; MG, microglia; NS, not significant. Source dataCC BY · PubMed Central
Extended Data Fig. 1 Generation of human APOE4-R136S knock-in mice and APOE4-R136S hiPSC lines by CRISPR/Cas-9-mediated gene editing. a, Schematic of generating human APOE4-R136S knock-in (E4-S/S-KI) mice using CRISPR/Cas-9-mediated gene editing. b, Schematic of gene editing strategy to generate the R136S mutation in human E4-KI mice (with the loxP sites, unused in current gene-editing strategy) and to generate APOE4-R136S (E4-S/S) hiPSC lines (without the loxP sites). c, DNA sequences of WT human APOE4 loci encoding for R136, designed sgRNA, and single-stranded oligodeoxynucleotides donor repair template encoding for S136 and silent mutation at PAM site for generating E4-S/S-KI mice. d, Sanger DNA sequencing of WT APOE4 and APOE4-S/S at and near the site encoding for residue 136 in E4-S/S-KI mice. e, Summary of on-target R136S editing in APOE4 and potential off-target mutation screening for knock-in mice. f, Representative immunofluorescent images of APOE (green) and GFAP (red) in CA1 hippocampal subfield in E4-KI and E4-S/S-KI mice at 12 months of age (scale bar, 50 μm). g, DNA sequences of WT APOE4 loci encoding for R136, designed sgRNA, and single-stranded oligodeoxynucleotides donor repair template encoding for S136 and silent mutation at PAM site for generating E4-S/S hiPSC lines. h, Sanger DNA sequencing of WT E4, E4-S/S, and E4-R/S at and near the site encoding for residue 136 in hiPSC lines. i, Summary of on-target R136S editing in APOE4 and potential off-target mutation screening in hiPSC lines. Experiments depicted in representative images in f were performed on n = 3 mice per genotype using 2 brain sections per mouse, with reproducible data. WT, wildtype; sgRNA, single guide RNA; ssODN, single-stranded oligodeoxynucleotide.
Extended Data Fig. 1 Generation of human APOE4-R136S knock-in mice and APOE4-R136S hiPSC lines by CRISPR/Cas-9-mediated gene editing. a, Schematic of generating human APOE4-R136S knock-in (E4-S/S-KI) mice using CRISPR/Cas-9-mediated gene editing. b, Schematic of gene editing strategy to generate the R136S mutation in human E4-KI mice (with the loxP sites, unused in current gene-editing strategy) and to generate APOE4-R136S (E4-S/S) hiPSC lines (without the loxP sites). c, DNA sequences of WT human APOE4 loci encoding for R136, designed sgRNA, and single-stranded oligodeoxynucleotides donor repair template encoding for S136 and silent mutation at PAM site for generating E4-S/S-KI mice. d, Sanger DNA sequencing of WT APOE4 and APOE4-S/S at and near the site encoding for residue 136 in E4-S/S-KI mice. e, Summary of on-target R136S editing in APOE4 and potential off-target mutation screening for knock-in mice. f, Representative immunofluorescent images of APOE (green) and GFAP (red) in CA1 hippocampal subfield in E4-KI and E4-S/S-KI mice at 12 months of age (scale bar, 50 μm). g, DNA sequences of WT APOE4 loci encoding for R136, designed sgRNA, and single-stranded oligodeoxynucleotides donor repair template encoding for S136 and silent mutation at PAM site for generating E4-S/S hiPSC lines. h, Sanger DNA sequencing of WT E4, E4-S/S, and E4-R/S at and near the site encoding for residue 136 in hiPSC lines. i, Summary of on-target R136S editing in APOE4 and potential off-target mutation screening in hiPSC lines. Experiments depicted in representative images in f were performed on n = 3 mice per genotype using 2 brain sections per mouse, with reproducible data. WT, wildtype; sgRNA, single guide RNA; ssODN, single-stranded oligodeoxynucleotide.CC BY · PubMed Central
Extended Data Fig. 2 Histopathological analyses of 10-month-old WT and 6-month-old PS19-E3, PS19-E4, PS19-E4-S/S, and PS19-E4-R/S mice. a, Representative images of 10-month-old WT mouse (n = 11) brain sections stained with AT8 monoclonal antibody to visualize p-Tau (scale bar, 500 µm), Sudan black to enhance hippocampal visualization (scale bar, 1 mm), GFAP and S100β to measure astrocytosis (scale bar, 500 µm), and Iba1 and CD68 to measure microgliosis (scale bar, 500 µm). b,c, Representative images of 6-month-old PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mouse brain sections stained with AT8 antibody for p-Tau (scale bar, 500 µm) (b) or Sudan black (scale bar, 1 mm) (c). d, Representative images of GFAP and S100β immunostaining for astrocytes and reactive astrocytes, respectively, as well as Iba1 and CD68 immunostaining for microglia and reactive microglia, respective, in the hippocampus of 6-month-old PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice. Scale bar, 500 µm. e–j, Quantification of % AT8 coverage area (e), hippocampal volume (f), % GFAP coverage area (g), % S100β coverage area (h), % Iba1 coverage area (i), and % CD68 coverage area (j). In a–j, PS19-E4, n = 12; PS19-E3, n = 11; PS19-E4-S/S, n = 12; PS19-E4-R/S, n = 9; n=mice. Experiments depicted in representative images in a–d were performed using 2 brain sections per mouse, with reproducible data. Throughout, data are expressed as mean ± s.e.m. Differences between groups were determined by ordinary one-way ANOVA followed with Dunnett’s multiple comparison test (e,f) or Welch’s ANOVA followed with Dunnett T3 multiple comparison test (g-j); comparisons of p ≤ 0.05 were labeled on graph. Source data
Extended Data Fig. 2 Histopathological analyses of 10-month-old WT and 6-month-old PS19-E3, PS19-E4, PS19-E4-S/S, and PS19-E4-R/S mice. a, Representative images of 10-month-old WT mouse (n = 11) brain sections stained with AT8 monoclonal antibody to visualize p-Tau (scale bar, 500 µm), Sudan black to enhance hippocampal visualization (scale bar, 1 mm), GFAP and S100β to measure astrocytosis (scale bar, 500 µm), and Iba1 and CD68 to measure microgliosis (scale bar, 500 µm). b,c, Representative images of 6-month-old PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mouse brain sections stained with AT8 antibody for p-Tau (scale bar, 500 µm) (b) or Sudan black (scale bar, 1 mm) (c). d, Representative images of GFAP and S100β immunostaining for astrocytes and reactive astrocytes, respectively, as well as Iba1 and CD68 immunostaining for microglia and reactive microglia, respective, in the hippocampus of 6-month-old PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice. Scale bar, 500 µm. e–j, Quantification of % AT8 coverage area (e), hippocampal volume (f), % GFAP coverage area (g), % S100β coverage area (h), % Iba1 coverage area (i), and % CD68 coverage area (j). In a–j, PS19-E4, n = 12; PS19-E3, n = 11; PS19-E4-S/S, n = 12; PS19-E4-R/S, n = 9; n=mice. Experiments depicted in representative images in a–d were performed using 2 brain sections per mouse, with reproducible data. Throughout, data are expressed as mean ± s.e.m. Differences between groups were determined by ordinary one-way ANOVA followed with Dunnett’s multiple comparison test (e,f) or Welch’s ANOVA followed with Dunnett T3 multiple comparison test (g-j); comparisons of p ≤ 0.05 were labeled on graph. Source dataCC BY · PubMed Central
Extended Data Fig. 3 Characterization of E4-S/S and E4-R/S hiPSC lines and neuronal differentiation. a, Karyogram of E4-S/S-A, E4-S/S-B, E4-S/S-C, E4-R/S-A, and E4-R/S-B hiPSC lines. For all cell lines, metaphase was examined for n = 20 cells per line. b, Representative immunofluorescent images of pluripotent stem cell markers Nanog, OCT3/4, SOX2, hNuclei, and TRA-1-60 from E4-S/S-A, E4-S/S-B, E4-S/S-C, E4-R/S-A, and E4-R/S-B hiPSC lines (scale bar, 100 μm). c, Representative immunofluorescent images of mature neuronal marker MAP2 (green) and DAPI (blue) in 3-week-old neurons differentiated from E4-S/S-A, E4-S/S-B, E4-S/S-C, E4-R/S-A, and E4-R/S-B hiPSC lines (scale bar, 20 μm). Representative images were selected from n = 3 fields of view from each cell line.
Extended Data Fig. 3 Characterization of E4-S/S and E4-R/S hiPSC lines and neuronal differentiation. a, Karyogram of E4-S/S-A, E4-S/S-B, E4-S/S-C, E4-R/S-A, and E4-R/S-B hiPSC lines. For all cell lines, metaphase was examined for n = 20 cells per line. b, Representative immunofluorescent images of pluripotent stem cell markers Nanog, OCT3/4, SOX2, hNuclei, and TRA-1-60 from E4-S/S-A, E4-S/S-B, E4-S/S-C, E4-R/S-A, and E4-R/S-B hiPSC lines (scale bar, 100 μm). c, Representative immunofluorescent images of mature neuronal marker MAP2 (green) and DAPI (blue) in 3-week-old neurons differentiated from E4-S/S-A, E4-S/S-B, E4-S/S-C, E4-R/S-A, and E4-R/S-B hiPSC lines (scale bar, 20 μm). Representative images were selected from n = 3 fields of view from each cell line.CC BY · PubMed Central
Extended Data Fig. 4 Measurement of neuronal uptake of Tau-488 by flow cytometry. a, Internal controls of Tau-488 uptake assay as measured by median fluorescent intensity at 488 nm of live cell population (left) or % 488-positive live cell population (right) via flow cytometry. All values normalized to E4 cells at 37 °C, collected over separate experiments. Samples treated with Tau-488 are colored in green and samples treated with unlabeled tau are colored in white. E4 at 37 °C with Tau-488, n = 10; E4 at 37 °C with unlabeled Tau, n = 4; E4 at 4 °C with Tau-488, n = 18; E4 at 37 °C no Tau CTL, n = 7; n=independent experiment with unique biological samples. b, Gating strategy for Tau uptake assay. First, cells were gated on forward scatter/side scatter (FSC/SSC). Cells were then gated on forward scatter height (FSC-H) versus area (FSC-A) to discriminate doublets. Dead cells were removed from the analysis using nuclear stain with DAPI, and positive cells were determined by gating on a control (no Tau-488 added) population. c, Scatter plots of flow cytometry analysis of live cell population at 488 nm (FITC-A versus FSC-A) for E4, E3, E4-S/S, or E4-R/S neuronal cultures with or without the treatment with 100 µg/ml heparin. Roughly 1×105 to 5×105 events were recorded in each experiment. The live cell population analyzed was roughly 5000 cells for each sample. In a, data are expressed as mean ± s.e.m. Source data
Extended Data Fig. 4 Measurement of neuronal uptake of Tau-488 by flow cytometry. a, Internal controls of Tau-488 uptake assay as measured by median fluorescent intensity at 488 nm of live cell population (left) or % 488-positive live cell population (right) via flow cytometry. All values normalized to E4 cells at 37 °C, collected over separate experiments. Samples treated with Tau-488 are colored in green and samples treated with unlabeled tau are colored in white. E4 at 37 °C with Tau-488, n = 10; E4 at 37 °C with unlabeled Tau, n = 4; E4 at 4 °C with Tau-488, n = 18; E4 at 37 °C no Tau CTL, n = 7; n=independent experiment with unique biological samples. b, Gating strategy for Tau uptake assay. First, cells were gated on forward scatter/side scatter (FSC/SSC). Cells were then gated on forward scatter height (FSC-H) versus area (FSC-A) to discriminate doublets. Dead cells were removed from the analysis using nuclear stain with DAPI, and positive cells were determined by gating on a control (no Tau-488 added) population. c, Scatter plots of flow cytometry analysis of live cell population at 488 nm (FITC-A versus FSC-A) for E4, E3, E4-S/S, or E4-R/S neuronal cultures with or without the treatment with 100 µg/ml heparin. Roughly 1×105 to 5×105 events were recorded in each experiment. The live cell population analyzed was roughly 5000 cells for each sample. In a, data are expressed as mean ± s.e.m. Source dataCC BY · PubMed Central
Extended Data Fig. 5 Quantifications and correlations between % coverage area of glial markers and hippocampal volume in tauopathy mice. a,b, Quantification of % GFAP (a) or % S100β (b) coverage area in hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice. c–f, Correlation between % GFAP coverage area and hippocampal volume in PS19-E4 (n = 31) (c), PS19-E3 (n = 23) (d), PS19-E4-S/S (n = 20) (e), and PS19-E4-R/S (n = 23) (f) mice at 10 months of age. g–j, Correlation between % S100β coverage area and hippocampal volume in PS19-E4 (n = 31) (g), PS19-E3 (n = 23) (h), PS19-E4-S/S (n = 20) (i), and PS19-E4-R/S (n = 23) (j) mice at 10 months of age. k,l, Quantification of % Iba1 (k) or % CD68 (l) coverage area in hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice. m–p, Correlation between % Iba1 coverage area and hippocampal volume in PS19-E4 (n = 31) (m), PS19-E3 (n = 23) (n), PS19-E4-S/S (n = 20) (o), and PS19-E4-R/S (n = 23) (p) mice at 10 months of age. q–t, Correlation between % CD68 coverage area and hippocampal volume in PS19-E4 (n = 31) (q), PS19-E3 (n = 23) (r), PS19-E4-S/S (n = 20) (s), and PS19-E4-R/S (n = 23) (t) mice at 10 months of age. For all quantifications in a,b,k,l, WT, n = 11; PS19-E4, n = 31; PS19-E3, n = 24; PS19-E4-S/S, n = 21; PS19-E4-R/S, n = 23. In a,b,k,l, data are expressed as mean ± s.e.m. and differences between groups were determined by Welch’s ANOVA followed with Dunnett T3 multiple comparison test; comparisons of p ≤ 0.05 were labeled on graph. Pearson’s correlation analysis (two-sided). Source data
Extended Data Fig. 5 Quantifications and correlations between % coverage area of glial markers and hippocampal volume in tauopathy mice. a,b, Quantification of % GFAP (a) or % S100β (b) coverage area in hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice. c–f, Correlation between % GFAP coverage area and hippocampal volume in PS19-E4 (n = 31) (c), PS19-E3 (n = 23) (d), PS19-E4-S/S (n = 20) (e), and PS19-E4-R/S (n = 23) (f) mice at 10 months of age. g–j, Correlation between % S100β coverage area and hippocampal volume in PS19-E4 (n = 31) (g), PS19-E3 (n = 23) (h), PS19-E4-S/S (n = 20) (i), and PS19-E4-R/S (n = 23) (j) mice at 10 months of age. k,l, Quantification of % Iba1 (k) or % CD68 (l) coverage area in hippocampus of 10-month-old PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice. m–p, Correlation between % Iba1 coverage area and hippocampal volume in PS19-E4 (n = 31) (m), PS19-E3 (n = 23) (n), PS19-E4-S/S (n = 20) (o), and PS19-E4-R/S (n = 23) (p) mice at 10 months of age. q–t, Correlation between % CD68 coverage area and hippocampal volume in PS19-E4 (n = 31) (q), PS19-E3 (n = 23) (r), PS19-E4-S/S (n = 20) (s), and PS19-E4-R/S (n = 23) (t) mice at 10 months of age. For all quantifications in a,b,k,l, WT, n = 11; PS19-E4, n = 31; PS19-E3, n = 24; PS19-E4-S/S, n = 21; PS19-E4-R/S, n = 23. In a,b,k,l, data are expressed as mean ± s.e.m. and differences between groups were determined by Welch’s ANOVA followed with Dunnett T3 multiple comparison test; comparisons of p ≤ 0.05 were labeled on graph. Pearson’s correlation analysis (two-sided). Source dataCC BY · PubMed Central
Extended Data Fig. 6 Quality control measures in snRNA-seq analysis of tauopathy mice with different APOE genotypes. a, Dot-plot depicting normalized average expression of selected cell identity marker genes for all 38 distinct hippocampal cell clusters from mice at 10 months of age. b, Cluster identity of 38 distinct cell types. c, The number of cells per cluster. d, The number of genes per cell in each cluster (± s.e.m). e, The number of nUMI (unique molecular identified) per cell in each cluster (± s.e.m). f, The % mitochondrial genes per cell in each cluster (± s.e.m). Number of cells (n) for each cell cluster can be found in c. For details of the analyses, see Methods.
Extended Data Fig. 6 Quality control measures in snRNA-seq analysis of tauopathy mice with different APOE genotypes. a, Dot-plot depicting normalized average expression of selected cell identity marker genes for all 38 distinct hippocampal cell clusters from mice at 10 months of age. b, Cluster identity of 38 distinct cell types. c, The number of cells per cluster. d, The number of genes per cell in each cluster (± s.e.m). e, The number of nUMI (unique molecular identified) per cell in each cluster (± s.e.m). f, The % mitochondrial genes per cell in each cluster (± s.e.m). Number of cells (n) for each cell cluster can be found in c. For details of the analyses, see Methods.CC BY · PubMed Central
Extended Data Fig. 7 snRNA-seq analysis of disease-protective neuronal clusters and disease-associated oligodendrocyte cluster. a,b, KEGG pathway enrichment dot-plot of top 20 pathways significantly enriched for DE genes of neuronal cluster 1 in PS19-E4-S/S (a) or PS19-E4-R/S (b) versus PS19-E4 mice at 10 months of age. c,d, KEGG pathway enrichment dot-plot of top 20 pathways significantly enriched for DE genes of neuronal cluster 7 in PS19-E4-S/S (c) or PS19-E4-R/S (d) versus PS19-E4 mice at 10 months of age. e,f, KEGG pathway enrichment dot-plot of top 20 pathways significantly enriched for DE genes of neuronal cluster 28 in PS19-E4-S/S (e) or PS19-E4-R/S (f) versus PS19-E4 mice at 10 months of age. g, Dot-plot normalized average expression of selected homeostatic and disease-associated oligodendrocyte (DAO) marker genes for oligodendrocyte clusters 2, 3, and 9. h, Volcano plot for top DE genes of oligodendrocyte cluster 9 versus cluster 2. i, Volcano plot for top DE genes of oligodendrocyte cluster 9 in PS19-E4-S/S versus PS19-E4 mice at 10 months of age. j, Representative images of Kirrel3+ Olig2+ oligodendrocytes in the hippocampus of 10-month-old PS19-E4 and PS19-E4-S/S mice at 10 months of age. Scale bars are 25 µm. k, Quantification of the number of Kirrel3+ Olig2+ oligodendrocytes (per mm2) within the dentate gyrus of hippocampus in 10-month-old PS19-E4 (n = 10) and PS19-E4-S/S (n = 10) mice at 10 months of age. In a–f, P-values are based on a two-sided hypergeometric test and are adjusted for multiple testing using the Benjamini-Hochberg method. Gene ratio represents the proportion of genes in the respective gene set that are deemed to be differentially expressed using the two-sided Wilcoxon Rank-Sum test as implemented in the FindMarkers function in Seurat. In h,i, horizonal dashed line indicates p-value = 0.05 and vertical dashed lines indicate log2 fold change = 0.4. All error bars represent s.e.m. Difference between groups in k was determined by unpaired, two-sided Welch’s t-test. DE, differentially expressed; NS, nonsignificant; FC, fold change. For details of the analyses, see Methods. Source data
Extended Data Fig. 7 snRNA-seq analysis of disease-protective neuronal clusters and disease-associated oligodendrocyte cluster. a,b, KEGG pathway enrichment dot-plot of top 20 pathways significantly enriched for DE genes of neuronal cluster 1 in PS19-E4-S/S (a) or PS19-E4-R/S (b) versus PS19-E4 mice at 10 months of age. c,d, KEGG pathway enrichment dot-plot of top 20 pathways significantly enriched for DE genes of neuronal cluster 7 in PS19-E4-S/S (c) or PS19-E4-R/S (d) versus PS19-E4 mice at 10 months of age. e,f, KEGG pathway enrichment dot-plot of top 20 pathways significantly enriched for DE genes of neuronal cluster 28 in PS19-E4-S/S (e) or PS19-E4-R/S (f) versus PS19-E4 mice at 10 months of age. g, Dot-plot normalized average expression of selected homeostatic and disease-associated oligodendrocyte (DAO) marker genes for oligodendrocyte clusters 2, 3, and 9. h, Volcano plot for top DE genes of oligodendrocyte cluster 9 versus cluster 2. i, Volcano plot for top DE genes of oligodendrocyte cluster 9 in PS19-E4-S/S versus PS19-E4 mice at 10 months of age. j, Representative images of Kirrel3+ Olig2+ oligodendrocytes in the hippocampus of 10-month-old PS19-E4 and PS19-E4-S/S mice at 10 months of age. Scale bars are 25 µm. k, Quantification of the number of Kirrel3+ Olig2+ oligodendrocytes (per mm2) within the dentate gyrus of hippocampus in 10-month-old PS19-E4 (n = 10) and PS19-E4-S/S (n = 10) mice at 10 months of age. In a–f, P-values are based on a two-sided hypergeometric test and are adjusted for multiple testing using the Benjamini-Hochberg method. Gene ratio represents the proportion of genes in the respective gene set that are deemed to be differentially expressed using the two-sided Wilcoxon Rank-Sum test as implemented in the FindMarkers function in Seurat. In h,i, horizonal dashed line indicates p-value = 0.05 and vertical dashed lines indicate log2 fold change = 0.4. All error bars represent s.e.m. Difference between groups in k was determined by unpaired, two-sided Welch’s t-test. DE, differentially expressed; NS, nonsignificant; FC, fold change. For details of the analyses, see Methods. Source dataCC BY · PubMed Central
Extended Data Fig. 8 snRNA-seq analysis of astrocyte and microglia subclusters. a, Feature plot showing relative levels of normalized human APOE gene expression across all 12 astrocyte subclusters by APOE genotype from mice at 10 months of age. b–d, KEGG pathway enrichment dot-plot of top 20 pathways significantly enriched for DE genes of astrocyte subclusters 3 (b), 5 (c), and 7 (d) versus other astrocyte subclusters. e, Feature plot showing relative levels of normalized human APOE gene expression across all 15 microglia subclusters by APOE genotype from mice at 10 months of age. f–h, KEGG pathway enrichment dot-plot of top 20 pathways significantly enriched for DE genes of microglia subclusters 2 (f), 8 (g), and 11 (h) versus other microglia subclusters. In b–d and f–h, P-values are based on a two-sided hypergeometric test and are adjusted for multiple testing using the Benjamini-Hochberg method. Gene Ratio represents the proportion of genes in the respective gene set that are deemed to be differentially expressed using the two-sided Wilcoxon Rank-Sum test as implemented in the FindMarkers function in Seurat.
Extended Data Fig. 8 snRNA-seq analysis of astrocyte and microglia subclusters. a, Feature plot showing relative levels of normalized human APOE gene expression across all 12 astrocyte subclusters by APOE genotype from mice at 10 months of age. b–d, KEGG pathway enrichment dot-plot of top 20 pathways significantly enriched for DE genes of astrocyte subclusters 3 (b), 5 (c), and 7 (d) versus other astrocyte subclusters. e, Feature plot showing relative levels of normalized human APOE gene expression across all 15 microglia subclusters by APOE genotype from mice at 10 months of age. f–h, KEGG pathway enrichment dot-plot of top 20 pathways significantly enriched for DE genes of microglia subclusters 2 (f), 8 (g), and 11 (h) versus other microglia subclusters. In b–d and f–h, P-values are based on a two-sided hypergeometric test and are adjusted for multiple testing using the Benjamini-Hochberg method. Gene Ratio represents the proportion of genes in the respective gene set that are deemed to be differentially expressed using the two-sided Wilcoxon Rank-Sum test as implemented in the FindMarkers function in Seurat.CC BY · PubMed Central
Extended Data Fig. 9 The R136S mutation reduces APOE4-induced neutral lipid accumulation in astrocytes and microglia in tauopathy mice. a, Representative images of neutral lipid visualized with BODIPY staining in GFAP+ astrocytes in the dentate gyrus of hippocampus in PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice at 10 months of age (example overlap shown in red arrows). b, Quantification of the fraction of BODIPY+ GFAP+ cells per total GFAP+ cells in the dentate gyrus of hippocampus in PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice at 10 months of age. c, Representative images of neutral lipid visualized with BODIPY staining in Iba1+ microglia in the dentate gyrus of hippocampus in PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice at 10 months of age (example overlap shown in white arrows). d, Quantification of the fraction of BODIPY+ Iba1+ cells per total Iba1+ cells in the dentate gyrus of hippocampus in PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice at 10 months of age. In a and c, scale bars are 50 µm. In b and d, PS19-E4, n = 10; PS19-E3, n = 10; PS19-E4-S/S, n = 10; PS19-E4-R/S, n = 10; n=mice. In b,d, data are expressed as mean ± s.e.m. and differences between groups were determined by Welch’s ANOVA followed with Dunnett T3 multiple comparison test; comparisons of p ≤ 0.05 were labeled on graph. DE, differentially expressed; AS, astrocyte; MG, microglia. For details of the analyses, see Methods. Source data
Extended Data Fig. 9 The R136S mutation reduces APOE4-induced neutral lipid accumulation in astrocytes and microglia in tauopathy mice. a, Representative images of neutral lipid visualized with BODIPY staining in GFAP+ astrocytes in the dentate gyrus of hippocampus in PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice at 10 months of age (example overlap shown in red arrows). b, Quantification of the fraction of BODIPY+ GFAP+ cells per total GFAP+ cells in the dentate gyrus of hippocampus in PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice at 10 months of age. c, Representative images of neutral lipid visualized with BODIPY staining in Iba1+ microglia in the dentate gyrus of hippocampus in PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice at 10 months of age (example overlap shown in white arrows). d, Quantification of the fraction of BODIPY+ Iba1+ cells per total Iba1+ cells in the dentate gyrus of hippocampus in PS19-E4, PS19-E3, PS19-E4-S/S, and PS19-E4-R/S mice at 10 months of age. In a and c, scale bars are 50 µm. In b and d, PS19-E4, n = 10; PS19-E3, n = 10; PS19-E4-S/S, n = 10; PS19-E4-R/S, n = 10; n=mice. In b,d, data are expressed as mean ± s.e.m. and differences between groups were determined by Welch’s ANOVA followed with Dunnett T3 multiple comparison test; comparisons of p ≤ 0.05 were labeled on graph. DE, differentially expressed; AS, astrocyte; MG, microglia. For details of the analyses, see Methods. Source dataCC BY · PubMed Central
Extended Data Fig. 10 PCA clustering of selected hippocampal cell clusters as well as astrocyte and microglia subclusters against all measured pathologies. Principal component analysis plot for hippocampal clusters 1, 6, 7, 9, 28, astrocyte subclusters 3, 5, 7, and microglia subclusters 2, 8, 11 against all measured pathologies (hippocampal volume and coverage areas of AT8, Iba1, CD68, GFAP, and S100β). PC1 and PC2 showed the overall relationship between clusters based on similarity of the estimated log odds ratio per unit change in six pathologies per cluster/subcluster. AS, astrocyte; MG, microglia. For details of the analyses, see Methods.
Extended Data Fig. 10 PCA clustering of selected hippocampal cell clusters as well as astrocyte and microglia subclusters against all measured pathologies. Principal component analysis plot for hippocampal clusters 1, 6, 7, 9, 28, astrocyte subclusters 3, 5, 7, and microglia subclusters 2, 8, 11 against all measured pathologies (hippocampal volume and coverage areas of AT8, Iba1, CD68, GFAP, and S100β). PC1 and PC2 showed the overall relationship between clusters based on similarity of the estimated log odds ratio per unit change in six pathologies per cluster/subcluster. AS, astrocyte; MG, microglia. For details of the analyses, see Methods.CC BY · PubMed Central

Verification

Claims the paper makes

131 Own result · 61 Own interpretation · 34 Own method · 18 Partly supported · 15 Supported · 3 Mixed support · 2 Not supported · 2 No full text of the sources · 1 Couldn’t confirm · 1 Not checked yet · 1 No basis given · accuracy not yet measured

Show all 171 claims

  • It may state two separate claims as one. 52
  • Probably not a substantial claim. 53
98 of 98
  • The paper relies on this without support.

    The PS19 tauopathy mouse model crossed with human APOE4-KI mice recapitulates APOE4-driven AD tau pathology sufficiently to test protection by R136S.

    “we cross-bred isogenic E4-KI and E4-R136S-KI mice with a widely used tauopathy mouse model expressing Tau-P301S (PS19 line)53.”
  • The paper relies on this without support.

    The engineered R136S mutation in mouse and hiPSC APOE4 has the same functional consequence of defective HSPG binding as the human APOE3-R136S.

    “we engineered the R136S mutation into the APOE4 allele background both in vivo in human APOE4 knock-in (E4-KI) mice52 and in vitro in human induced pluripotent stem cells (hiPSCs) derived from a patient with AD with homozygous APOE4 (ref. 6).”
  • The paper relies on this without support.

    E4 neuron-conditioned medium increases p-Tau in E3 and E4-S/S neurons through exogenous APOE4 rather than other secreted factors.

    “Thus, the presence of exogenous APOE4, even with much higher endogenous APOE4-S/S levels (Fig. 2a,b), was sufficient to increase p-Tau accumulation in E4-S/S neurons to levels similar to those seen in E4 neurons.”
  • The paper relies on this without support.

    The top predicted off-target sites are sufficient to rule out off-target effects of CRISPR editing on phenotypes.

    “found no changes within the top predicted potential off-target sites”
  • The paper relies on this without support.

    Correlations between gliosis markers and hippocampal volume reflect a causal contribution of gliosis to hippocampal degeneration.

    “suggesting the contribution of astrocytosis to hippocampal degeneration”
  • The paper relies on this without support.

    The APOE4-expressing tauopathy mouse model and hiPSC-derived neurons from a patient with AD are valid models for APOE4-driven AD pathologies relevant to human disease.

    “In the present study, we demonstrate the protective power of the R136S mutation against APOE4-driven effects in a tauopathy mouse model and in hiPSC-derived neurons from a patient with AD, neither of which exhibit extensive Aβ pathology.”
  • The authors state this as an assumption.

    Data distribution was assumed to be normal, but this was not formally tested.

    “Data distribution was assumed to be normal, but this was not formally tested.”
7 of 7

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.

Each contradiction between two passages listed here was confirmed by an independent check.

1 of 1

Citations

Cites

Show all 103 references

Cited by

Sorted by citations

Items 101–115 of 115

Metadata

Authors
Nelson, M. R., Liu, P., Agrawal, A., Yip, O., Blumenfeld, J., Traglia, M., Kim, M. J., Koutsodendris, N., Rao, A., Grone, B., Hao, Y., Yoon, S. Y., Xu, Q., De Leon, S., Choenyi, T., Thomas, R., Lopera, F., Quiroz, Y. T., Arboleda-Velasquez, J. F., Reiman, E. M., Mahley, R. W., & Huang, Y.
Affiliations
  • Gladstone Institutes (US)
  • Universidad de Antioquia (CO)
  • Harvard University (US)
  • University of Arizona (US)
Corresponding author
Huang, Yadong
Venue
Nature Neuroscience
Published
2023-11-13
Volume
26
Issue
12
Pages
2104-2121
ISSN
  • 1097-6256 (ISSN-L)
  • 1546-1726
OpenAlex
W4388625926
Cited by
120
References
103
Citations by year
2026: 34 · 2025: 47 · 2024: 37 · 2023: 2
FWCI
11.71
Type
journal article
Language
en
Open access
Locations
  • Nature Neuroscience (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • eScholarship (California Digital Library) (submittedVersion)
  • eScholarship (California Digital Library) (submittedVersion, cc-by)
  • Repositorio Digital - Academia Colombiana de Ciencias Exactas, Fisicas y Naturales (submittedVersion, other-oa)
  • UA Campus Repository (The University of Arizona) (submittedVersion, cc-by)
  • Greater South Information System
Funders
U.S. Department of Health & Human Services | NIH | National Institute on Aging (R01AG071697, RF1AG076647, P01AG073082, F31AG074690, F31AG074672, P3AG072980, RF1AG059751)
Topics
  • Life Sciences › Neuroscience › Cellular and Molecular Neuroscience › Alzheimer's disease research and treatments
  • Health Sciences › Medicine › Psychiatry and Mental health › Dementia and Cognitive Impairment Research
  • Life Sciences › Neuroscience › Neurology › Neuroinflammation and Neurodegeneration Mechanisms
Keywords
APOE ε4, Alzheimer's disease, tau pathology, neurodegeneration, neuroinflammation, iPSC-derived neurons, single-nucleus RNA sequencing, gene dosage, neuroprotection, therapeutic target
MeSH terms
Neuroinflammatory Diseases, Alzheimer Disease, Animals, Apolipoproteins E, Humans, Mutation, Tauopathies, Mice, Apolipoprotein E3, Apolipoprotein E4

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.

A language model drew the claims from the paper on 2026-10-31, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-31 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Cited sources are checked with a language model’s help (gpt-6-sol).

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