Humanized APOE genotypes influence lifespan independently of tau aggregation in the P301S mouse model of tauopathy
Acta Neuropathologica Communications volume 11, Article number: 99 (2023)
Apolipoprotein (APOE) E4 isoform is a major risk factor of Alzheimer’s disease and contributes to metabolic and neuropathological abnormalities during brain aging. To provide insights into whether APOE4 genotype is related to tau-associated neurodegeneration, we have generated human P301S mutant tau transgenic mice (PS19) that carry humanized APOE alleles (APOE2, APOE3 or APOE4). In aging mice that succumbed to paralysis, PS19 mice homozygous for APOE3 had the longest lifespan when compared to APOE4 and APOE2 homozygous mice (APOE3 > APOE4 ~ APOE2). Heterozygous mice with one human APOE and one mouse Apoe allele did not show any variations in lifespan. At end-stage, PS19 mice homozygous for APOE3 and APOE4 showed equivalent levels of phosphorylated tau burden, inflammation levels and ventricular volumes. Compared to these cohorts, PS19 mice homozygous for APOE2 showed lower induction of phosphorylation on selective epitopes, though the effect sizes were small and variable. In spite of this, the APOE2 cohort showed shorter lifespan relative to APOE3 homozygous mice. None of the cohorts accumulated appreciable levels of phosphorylated tau compartmentalized in the insoluble cell fraction. RNAseq analysis showed that the induction of immune gene expression was comparable across all the APOE genotypes in PS19 mice. Notably, the APOE4 homozygous mice showed additional induction of transcripts corresponding to the Alzheimer’s disease-related plaque-induced gene signature. In human Alzheimer’s disease brain tissues, we found no direct correlation between higher burden of phosphorylated tau and APOE4 genotype. As expected, there was a strong correlation between phosphorylated tau burden with amyloid deposition in APOE4-positive Alzheimer’s disease cases. Overall, our results indicate that APOE3 genotype may confer some resilience to tauopathy, while APOE4 and APOE2 may act through multiple pathways to increase the pathogenicity in the context of tauopathy.
Both sporadic and familial forms of Alzheimer’s disease (AD) are associated with the presence of the E4 allele of Apolipoprotein E (APOE) [1, 2]. Indeed, relative to the most prevalent E3 allele representing the ‘reference allele’, the E4 allele significantly increases the risk of AD and lowers the age of dementia onset in a gene-dose dependent manner. On the other hand, the E2 allele is associated with increased longevity and protection from AD , though presence of E2 allele can often lead to hyperlipoproteinemia and metabolic pathologies [4, 5].
APOE is a lipoprotein that functions as a cholesterol carrier and its three major isoforms—E2, E3 and E4—differ in sequence at two positions . The pathogenic effects of APOE4 in aging and Alzheimer’s dementia is thought to be mediated by its role in various metabolic and immune pathways . APOE4 is associated with higher burden of amyloid β (Aβ) burden . However, there are some divergent datasets regarding the direct involvement of APOE alleles in tau pathogenesis. Previous studies in human tau transgenic PS19 mice showed that APOE4 triggers tau-associated neurodegeneration and unique phosphorylation patterns [7,8,9], while another report found that APOE2 exacerbates pathologic tau inclusions in a hyper-expression model of tau . In the K18-tau aggregate seeded PS19 mice, we showed that APOE3 homozygosity (E3H) exacerbates the prevalence of hyperphosphorylated tau (phospho-tau) in the early stages of seeded tauopathy (E3H > E4H ~ E2H) . Given these disparate observations, we wanted to examine whether presence of a particular human APOE isoform alters survival and tau pathology in aged PS19 mice. We found that PS19 mice carrying E3 allele (PS/E3H) survived the longest (E3H > E4H ~ E2H) and that the burden of tau pathology was equivalent between PS/E3H and PS/E4H (PS19 mice homozygous for APOE4) mice at end-stage. RNAseq analysis showed that PS/E4H mice showed AD-associated profiles, relative to PS/E3H mice and PS19 mice homozygous for APOE2 (PS/E2H). A survey using limited human AD cases demonstrated that phospho-tau levels were similar between APOE4 and non-APOE4 individuals. Interestingly, phospho-tau burden showed positive correlation with amyloid β (Aβ) deposits in APOE4 individuals, but not in non-APOE4 AD patients. Collectively, our data suggests that APOE3 may influence resilience to tau-associated phenotype in PS19 mice while presence of APOE4 was pathological in the same context.
Materials and methods
Mice and study design
Mouse husbandry and experimental procedures were performed in accordance with the protocols and policies approved by the Institutional Animal Care and Use Committee at the University of Florida. All mice were maintained under a 12-h light/dark cycle and had access to water and food ad libitum. The mice were maintained under specific pathogen-free conditions. PS19 mice were obtained from Jackson Labs and maintained on a B6/C3H background as heterozygotes for P301S tau transgene, developing age-progressive hindlimb paralysis at 9–12 months of age . APOE targeted replacement (TR) mice were obtained from Duke University and were maintained as homozygotes on C57BL6 background [13,14,15]. PS19 mice were mated with APOE homozygote mice to produce N1 cohorts that were heterozygous for APOE. N2 cohorts of PS19 mice that are heterozygous or homozygous for APOE were generated by back-crossing the N1 generation mice with corresponding APOE homozygous mice. APOE and tau levels in these N1 and N2 cohorts were described in a previous study . Cohorts were set up to be aged until the animals reached hindlimb paralysis. At endpoint, all euthanasia was performed with intracardiac perfusion of cold phosphate buffered saline (PBS) containing heparin and the brains were fixed in 10% normal buffered formalin (Fisher Scientific). The fixed brains were then sliced sagitally and processed for paraffin embedding.
Human brain samples—demographics and APOE genotyping protocol
All brains were obtained from the University of Florida Neuromedicine Human Brain and Tissue Bank. Information on individual patients are summarized in Additional file 6: Table S1.
For APOE genotyping, DNA was extracted from formalin fixed paraffin embedded (FFPE) sections or frozen brain tissue using DNA extraction kits (GeneRead DNA FFPE Kit, Qiagen; and QIAamp DNA Mini Kit, Qiagen). Two single nucleotide polymorphisms (SNPs) within the APOE gene, rs7412 and rs429358, were detected using the Human TaqMan™ SNP Genotyping Assay (Applied Biosystems, Cat. no. 4351376) and TaqMan™ Genotyping Master Mix (Applied Biosystems, Cat. no. 4371355). The relative fluorescent units (RFU) were measured by performing qRT-PCR. The frequencies of APOE alleles and corresponding genotypes were obtained by calculating the ratio and combination of RFUs for both SNPs.
Immunohistochemical analysis of mouse brain sections
FFPE slides with mouse brain sections were deparaffinized and probed with primary antibodies as described before . For antigen retrieval, slides receiving AT8 (1:1000; Invitrogen MN1020), AT100 (1:1000; Invitrogen MN1060), AT270 (1:1000; Invitrogen MN1050), AT180 (1:1000; Invitrogen MN1040) and glial fibrillary acidic protein (GFAP) (1:3000; Cell Signaling) antibodies were steamed in water at high pressure for 15 min, while Iba-1 (1:2000; Wako) antibody slides were steamed in citrate buffer pH 6.0 (Target Retrieval Solution, Dako). Slides were incubated in 3% hydrogen peroxide (Fisher Scientific) for 20 min to block endogenous peroxidase activity and then washed three times in PBS for 5 min each. Slides were then blocked in 2% fetal bovine serum (FBS) (Hyclone, GE) for 45 min before incubating in primary antibody diluted in block solution overnight at 4 °C. The following day, slides were washed and appropriate secondary antibody (ImmPRESS Polymer Reagent, Vector Labs) was applied for 30 min at room temperature. Following PBS washes, color was developed using 3,3’-diaminobenzidine (Vector DAB, Vector Labs) and slides were counterstained with haematoxylin (Vector Labs). Next, brain sections were dehydrated in a series of ethanol, cleared in xylene, mounted in Cytoseal-60 media (Fisher Scientific) and coverslipped.
Analysis of histochemical images from mice brain
Immunostained images were captured using a Scanscope XT image scanner (Aperio, Vista, CA). Percent immunoreactivity was computed using the Aperio Positive Pixel Count program (Aperio, Vista, CA). The data is shown as the average %immunoreactivity ± S.E.M. per group. Statistical comparisons were conducted using 1-way ANOVA (GraphPad Prism 7).
Estimation of ventricular volume
FFPE tissues were sectioned sagittally, starting 1 mm from the midline. Initial slides were stained with haematoxylin and eosin, their location identified using Paxinos and Franklin’s mouse brain atlas and further slides were prepared according to desired bregma locations. 5 µm sections were collected spaced at regular intervals of 25 µm followed by an interval of 250 µm. This pattern was repeated until reaching the target area, which was approximately 3.2 mm from the midline bregma. Slides were stained with haematoxylin and eosin at target area. Assessments of ventricle volume were made at two distances from bregma. One measurement of ventricle volume was made at 2.8–2.76 mm from bregma, and a second measurement was made at 3.2 mm from bregma. Using an Aperio Image Scope application to view the images, the entire perimeter of the lateral ventricles that were visible in the images was outlined and the area of the marked areas were recorded. In more medial sections where the ventricles were discontinuous, the total area of the ventricles was summed.
Analysis of mouse brains using immunoblotting
Frozen hemibrains were cryo-pulverized and used for tau extraction as detailed earlier . Briefly, brain fractions S1 and S3 containing tau with differential soluble properties were obtained by treating the brain homogenates with a sequential series of buffers with increasing detergent concentrations—Tris buffered saline, high salt buffer containing 10% sucrose and finally 1% sarkosyl solution. The S1 fraction refers to clarified homogenate prepared in Tris buffered saline while S3 fraction refers to the pellet obtained following 1% sarkosyl incubation. The S3 pellet was further solubilized in 4 M urea and 2% SDS and used in western blotting. All steps in the protocol included addition of protease and phosphatase inhibitors (Pierce Protease and Phosphatase Inhibitor Mini Tablets, Thermo Scientific). Immunoblotting methods were followed as described earlier . Membranes were incubated with primary antibodies [Iba-1 (1:250; Novus Biologicals), GFAP (1:1000; Cell Signaling), GAPDH (1:1000; Abcam), AT8 (1:1000; Invitrogen MN1020), AT100 (1:1000; Invitrogen MN1060), AT270 (1:1000; Invitrogen MN1050), AT180 (1:1000; Invitrogen MN1040), PHF1 (1:1000; Peter Davies), total tau (1:1000; Abcam ab254256] and signals detected using multiplex Li-Cor Odyssey Infrared Imaging system (Li-Cor Biosciences, Lincoln, NE). Relative band intensity was quantified using ImageJ software (NIH).
Analysis of histochemical images from human brain
Human brain tissue sections were received from the UF Neuromedicine Human Brain and Tissue Bank (UF HBTB) following protocols approved by the Institutional Review Board of the University of Florida. Staging of AD neuropathologic changes and associated co-pathologies was conducted per guidelines published by the National Institute on Aging-Alzheimer’s Association  by a board-certified neuropathologist (SP). 5 μm thick tissue sections of FFPE brain tissue specimens were deparaffinized and rehydrated. Heat-induced epitope retrieval was performed in a pressure cooker (Tintoretriever, Bio SB) for 15 min at high pressure in a 0.05% Tween-20 solution. Endogenous peroxidase was quenched by incubation of sections in 1.5% hydrogen peroxide/0.005% Triton-X-100 diluted in sterile PBS (Invitrogen) for 20 min, following multiple washes in tap water and subsequently, 0.1 M Tris, pH 7.6. Non-specific antibody binding was minimized with sections incubated in 2% FBS/0.1 M Tris, pH 7.6. Primary antibodies were diluted in 2% FBS/0.1 M Tris, pH 7.6 (anti-Aβ antibody Ab5 [Todd E. Golde], 1:1000; anti-tau (AT8), 1:5000, invitrogen MN1020; anti-α-synuclein (94-3A10) [Benoit I. Giasson], 1:10,000). Sections were incubated with primary antibody over night at 4 °C, washed one time in 0.1 M Tris, pH 7.6, followed by 2% FBS/0.1 M Tris, pH 7.6 for 5 min, incubated in goat anti-rabbit IgG horseradish peroxidase (HRP) conjugated secondary antibody (Millipore Sigma) for 1 h, additionally washed one time in 0.1 M Tris, pH 7.6, followed by 2% FBS/0.1 M Tris, pH 7.6 for 5 min, and incubated in VectaStain ABC Peroxidase HRP Kit (diluted in 2% FBS/0.1 M Tris, pH 7.6 at 1:1000) for 1 h. After a final wash in 0.1 M Tris, pH 7.6 for 5 min, immunocomplexes were visualized using 3,3′-diaminobenzidine (DAB). Tissue sections were counterstained with hematoxylin (Sigma Aldrich, St. Louis) for 2 min, dehydrated and cover slipped using Cytoseal 60 mounting medium (Thermo Scientific).
For analysis of stains, slides were scanned on an Aperio AT2 scanner (Leica biosystems) at 20 × magnification. Digital slides were analyzed using the QuPath platform (version 0.3.2) . Cortex and white matter were annotated for regional analysis. After exclusion of tissue and staining artifacts, the ‘Positive Pixel Detection’ tool (e.g., Downsample factor 2, Gaussian sigma 1 μm, Hematoxylin threshold (‘Negative’) 1.5 OD units, DAB threshold (‘Positive’) 0.2 OD units) was used to determine the percentage of area covered by pathological staining. ‘Positive Pixel Detection’ thresholds were adjusted for individual stains.
RNA isolation from frozen brains and bulk Rnaseq analysis
The number of samples were 3 mice each for parental E2H, E3H and E4H cohorts and 4 mice each for PS/E2H, PS/E3H and PS/E4H cohorts. All mice were male, except 1 female for E4H, 1 female for PS/E2H and 2 females for PS/E4H cohort. RNA was extracted from frozen, pulverized forebrains using TRIzol reagent (Invitrogen). Extracted RNA was cleaned using RNeasy mini extraction kit with an on-column DNase treatment (Qiagen). RNA quality was quantified with the Qubit RNA HS assay (Invitrogen) and transferred to LC Sciences (Houston, TX) for library construction and Rnaseq. Poly(A) RNA sequencing library was prepared following Illumina’s TruSeq-stranded-mRNA sample preparation protocol. RNA integrity was checked with Agilent Technologies 2100 Bioanalyzer. Poly(A) tail-containing mRNAs were purified using oligo-(dT) magnetic beads with two rounds of purification. After purification, poly(A) RNA was fragmented using divalent cation buffer in elevated temperature. Quality control analysis and quantification of the sequencing library were performed using Agilent Technologies 2100 Bioanalyzer High Sensitivity DNA Chip. Paired-ended sequencing was performed on Illumina’s NovaSeq 6000 sequencing system. Firstly, Cutadapt  and perl scripts in house were used to remove the reads that contained adaptor contamination, low quality bases and undetermined bases. Then sequence quality was verified using FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). We used HISAT2  to map reads to the genome of ftp://ftp.ensembl.org/pub/release-101/fasta/mus_musculus/dna/. The mapped reads of each sample were assembled using StringTie . Then, all transcriptomes were merged to reconstruct a comprehensive transcriptome using perl scripts and gffcompare. After the final transcriptome was generated, StringTie ] and ballgown (http://www.bioconductor.org/packages/release/bioc/html/ballgown.html) was used to estimate the expression levels of all transcripts. StringTie  was used to perform expression level for mRNAs by calculating FPKM (fragments per kilobase of exon per million mapped fragments). mRNAs differential expression analysis were performed by R package DESeq2  between two different groups (and by R package edgeR  between two samples). The mRNAs with the parameter of false discovery rate (FDR) below 0.05 and absolute fold change ≥ 2 were considered differentially expressed mRNAs. Gene Set Enrichment Analysis (GSEA) with GO and KEGG databases were done using gene sets with |NES|≥ 1 & NOM p-value < 0.05 & FDR q-value < 0.25, where NES = normalized enrichment score, NOM = NOM and FDR = false discovery rate.
Cell profile and microglial profile analysis
Gene lists identifying cell-types within the brain were described previously [24,25,26,27,28,29,30]. The geometric mean of the FPKM for genes identified for each cell or cellular sub-type was calculated from each animal. Group means were calculated, and between group significance values determined by 1-way ANOVA. Outliers were removed from the group if their value fell outside of 1.5 times the inter-quartile range (IQR).
Immunohistochemical data was analyzed using 1-way ANOVA unless otherwise indicated in figure legend. Outliers, if applicable, were removed using Rout’s test with Q = 1%. The neuropathological scoring was done in a blinded manner. All data was assembled using Adobe Photoshop Elements.
APOE isoform dependent survival in PS19 mice homozygous for APOE
As described earlier , mice carrying humanized APOE alleles [13,14,15] were bred with P301S tau transgenic Line PS19 mice  to generate tau transgenic mice that are heterozygous for APOE (PS/E2h, PS/E3h and PS/E4h) or homozygous (PS/E2H, PS/E3H and PS/E4H) for human APOE. PS19 mice show progressive paralysis as a result of accumulation of tau pathology in the spinal cord and brain . In our colony, the median age to humane endpoint of bilateral hindlimb paralysis for PS19 mice on the B6/C3H background, carrying mouse Apoe, is 308 days (n = 7 males, 9 females) (Fig. 1A). We examined whether the presence of specific human APOE or the gene dosage of APOE influenced lifespan in PS19 mice (Fig. 1A). The median number of days to humane endpoint of bilateral hindlimb paralysis for PS/E2H, PS/E3H and PS/E4H mice was 332, 381 and 318 days, respectively (Fig. 1A). The longer average survival for PS/E3H mice was statistically different from PS/E2H mice (p = 0.015), PS/E4H mice (p = 0.008) and PS19 mice (p = 0.003). There was no statistical difference in median age to paralysis between PS/E2H, PS/E4H and PS19 cohorts. Previous studies, using coronal brain sections of PS19 crossed with humanized APOE mice, showed that APOE4 exacerbates ventricular enlargement . In our study, we used sagittal sections from two selected brain areas, one a medial and another a lateral location with respect to midline to explore changes in ventricular volume. In our cohorts of PS19xAPOE mice, we observed high variability in ventricular volumes and did find any statistically significant differences in ventricular volume between the three genotypes (Additional file 1: Fig. S1A, B).
Phosphorylated tau burden in PS19 mice homozygous for APOE
We investigated the levels of tau and phospho-tau in brains of these paralyzed mouse cohorts (Fig. 1, Additional file 2: Fig. S2). Using an antibody against total tau, we observed that all of our mouse lines showed comparable tau levels in both the detergent soluble (S1) and detergent insoluble (S3) fractions (Fig. 1B). Using AT8 (pSer202/Thr205) antibody staining in brain sections, we observed lower levels of phospho-tau immunoreactivity in PS/E2H mice compared to PS/E3H (p < 0.05) and PS/E4H (p < 0.01) mice in the hippocampus (Fig. 1C). Cortical AT8 immunoreactivity was somewhat lower in PS/E2H mice compared to PS/E3H and PS/E4H mice, though this did not reach statistical significance due to high variability (Fig. 1C). Immunoblotting of soluble and detergent-insoluble fractions from hemibrains showed that AT8-immunopositive detergent-soluble tau was equivalent between the three APOE genotypes (Fig. 1D). Notably, we failed to detect insoluble tau showing AT8-site phosphorylation by immunoblotting in these mice (Additional file 2: Fig. S2G). These findings suggest that most of the AT8 immunoreactivity seen in the cortex of PS19 mice is from phosphorylated soluble tau.
Next, we investigated whether other tau phosphorylation epitopes showed any differential association with APOE genotypes. Immunohistochemistry analysis of pThr212/Ser214 (AT100), pThr181 (AT270) and pThr231 (AT180) burden showed equivalent levels between PS/E2H, PS/E3H and PS/E4H mice (Additional file 2: Fig. S2A-B). Immunoblotting of soluble and detergent-insoluble fractions with these antibodies showed variable patterns of immunoreactivity (Additional file 2: Fig. S2C–G). Consistent with the immunohistological staining, the level of AT100 immunostaining in soluble fractions was low (Additional file 2: Fig. S2C). In the limited samples we analyzed here, soluble AT100 levels were somewhat lower in PS/E2H mice as compared to other genotypes (Additional file 2: Fig. S2C; p < 0.05 relative to PS/E3H and PS/E4H). Robust levels of AT270 and AT180 immunoreactivity were detected in soluble fractions from all genotypes (Additional file 2: Fig. S2D, E). We also analyzed biochemical levels of PHF1 (Ser396/404), finding robust staining in soluble fractions as well as detectable amounts in detergent-insoluble fractions (Additional file 2: Fig. S2F). In the detergent-insoluble fractions, AT100-, AT180- and AT8-positive phospho-tau were undetectable (Additional file 2: Fig. S2C, E, G). Presence of total tau on these same blots confirmed that these phospho-tau epitopes were specifically not induced. Levels of silver-positive Gallyas staining showed high variability in forebrains of these mice and thus were not considered for analysis (data not shown). Collectively, these data suggest that most of the cortical phospho-tau present in the brains of these mice was present as detergent soluble tau. Since our data shows that PS/E3H and PS/E4H mice had equivalent phospho-tau pathology, but PS/E4H mice reached end-stage (hindlimb paralysis) earlier than PS/E3H mice, this suggests that the E3 allele imparts some resilience to tau-mediated pathogenesis.
Gliosis in PS19 mice homozygous for APOE
We examined hippocampal and cortical microgliosis and astrogliosis in PS19 mice carrying different APOE alleles, using Iba-1 and GFAP antibodies respectively. Hippocampal burden of Iba-1 immunoreactivity was somewhat lower in PS/E2H mice compared to PS/E3H and PS/E4H mice (PS/E2H vs PS/E3H: p < 0.05), with the latter two genotypes showing equivalent levels (Fig. 2A). There was no difference in cortical microglia staining in the three cohorts (Fig. 2A). Immunoblotting analyses of Iba-1 levels from the whole forebrain showed equivalent levels among the three APOE genotypes of PS19 mice (Additional file 3: Fig. S3A), indicating region-specific microgliosis in these mice.
We observed lower hippocampal burden of GFAP in PS/E2H mice compared to PS/E3H mice (p < 0.01), with no difference observed between PS/E3H and PS/E4H mice (Fig. 2B). There was no difference in cortical astroglia staining among the three cohorts of bigenic mice (Fig. 2B). GFAP immunoblotting showed similar levels in PS19 mice bearing mouse Apoe or human APOE2, APOE3 or APOE4 (Additional file 3: Fig. S3B). Notably, the PS/E3H bigenic mice showed higher astrogliosis compared to the corresponding parental E3H mice (Additional file 3: Fig. S3B; p < 0.05), while the other group comparisons (PS/E2H vs E2H and PS/E4H vs E4H) did not reach statistical significance. Our results indicate that presence of the tau transgene induces variable levels of astrogliosis in all parental APOE genotypes, and most consistently in E3H mice.
Survival, tau burden and gliosis in PS19 mice heterozygous for APOE
We investigated the lifespan of PS19 mice carrying one human APOE allele and one endogenous Apoe allele. In these heterozygous APOE mice, we did not observe any differential time to paralysis based on APOE genotype. The median ages of survival of the PS/E2h, PS/E3h and PS/E4h mice were 358, 325 and 365 days respectively (Additional file 4: Fig. S4A, B). Interestingly, the PS/E2h mice had a longer lifespan compared to its APOE homozygous counterpart, PS/E2H mice (p = 0.028) (Additional file 4: Fig. S4B), showing the detrimental influence of increased APOE2 gene dose on age to paralysis. Likewise, the PS/E4h mice had a longer lifespan compared to PS/E4H mice (p = 0.013) (Additional file 4: Fig. S4B). Although the median age to paralysis for PS/E3h mice was substantially shorter than that of PS/E3H mice (325 vs 381 days), the difference was not statistically significant. Within the APOE heterozygous PS19 cohorts, we did not observe any statistically significant changes in the cortical or hippocampal burden of phospho-tau or Iba-1 (Additional file 4: Fig. S4C, D). Cortical (but not hippocampal) levels of GFAP was lower in PS/E2h mice compared to PS/E3h mice (p < 0.05), with no difference noted between PS/E3h and PS/E4h mice (Additional file 4: Fig. S4E). Whether this difference is biologically meaningful is uncertain, given the high variability in reactivity within these cohorts.
RNAseq reveals APOE isoform-based differences in PS19 mice
To illuminate the potential mechanisms by which APOE could influence the neuropathology in PS19 mice, we conducted bulk RNAseq analysis of the forebrains of PS19 mice homozygous for human APOE and corresponding parental APOE homozygous cohorts (Additional file 6: Table S2). Comparison of the three parental APOE homozygous cohorts showed clear APOE genotype-based differences in gene expression (Fig. 3; Additional file 6: Table S2A-C). Relative to E3H mice, the chemokine genes, Ccl21b and Ccl21d, were downregulated and the cysteine protease, Capn11, was upregulated in E2H mice. Gene Ontology (GO) enrichment analysis showed that multiple GO terms associated with immune function (chemokine receptor binding, chemotaxis and cellular response to cytokines) were differentially altered in E2H mice. Analysis of KEGG pathways also showed consistent enrichment in immune pathways in E2H mice relative to E3H mice (Fig. 3B–D). Fewer differences were observed between E4H and E3H mice at this age. Notable genes that were upregulated in E4H mice (relative to E3H mice) were Serpina3n (cysteine protease inhibitor) and Eps8l1 (signaling adaptor). Enrichment of GO terms and KEGG pathways showed few significant events – most of which were related to metabolic pathways and cellular response to DNA damage (Fig. 3F–G). Comparison of E4H to E2H revealed several immune genes upregulated in E4H mice, such as Ccl121b and Ccl21d (Fig. 3H). The GO term and KEGG pathway profiles in E4H vs E2H comparison were similar to E2H vs E3H comparison, with immune pathways identified as enriched in both the analyses (Fig. 3I, J; C, D). Notably, while immune genes were downregulated in E2H mice (relative to E3H mice), these genes were upregulated in E4H mice (relative to E2H mice).
There were substantial differences in PS19 mice carrying the three APOE alleles (Additional file 5: Fig. S5A). Genes that were differentially upregulated in the comparisons of parental APOE lines were also found in the PS19 tau transgenics carrying the corresponding APOE allele (examples include Eps8l1 in E2H vs E3H, Eps8l1 in E4H vs E3H and Ccl21d in E4H vs E2H comparisons) (Fig. 3B, E, H vs Additional file 5: Fig. S5 B, C, D; Additional file 6: Table S2D-F). We normalized the within-group differences to gain insights into the differential gene expression (DEG) patterns that are determined by APOE genotype or presence of tau transgene (Additional file 5: Fig. S5E–H). The common DEG list showed few genes that are altered and few common pathways were found to be enriched when we investigated differences between E2H vs E3H mice or E4H vs E3H mice, after allowing for the effect of PS19 genotype to be normalized (Additional file 5: Fig. S5E, G). For example, relative to APOE3, NADP + transhydrogenase activity and myeloid cell chemotaxis were downregulated in APOE2 genotypes, whereas no upregulated pathways were identified in this comparison. On the other hand, relative to APOE3, dUMP metabolic process was down regulated and translational fidelity pathway was upregulated in APOE4 mice. We further tested for the effect of the PS19 genotype, after normalizing for the differential contribution of APOE. The presence of tau transgene led to upregulation of a larger number of genes in these mice, with 127 genes commonly upregulated in E2H vs E3H comparison and 252 genes upregulated in E4H vs E3H comparison (Additional file 5: Fig. S5F–H). Expression of tau transgene leads to a robust immune response and stress response that are upregulated in both APOE2 (relative to APOE3) and APOE4 (relative to APOE3) genotypes (Additional file 5: Fig. S5F–H). This suggests that tau expression imparts an overwhelming immune phenotype in these APOE homozygous mice. There were fewer genes that were commonly downregulated in these comparisons (Additional file 5: Fig. S5F–H), one being the keratinocyte activation pathway that was downregulated in both these DEGs (E2H vs E3H and E4H vs E3H).
Tau transgene induced DEG in APOE4 and APOE3 homozygous PS19 mice
We wanted to compare the individual DEG in PS19 mice within each APOE genotype to understand how tau transgene expression is related to eventual tau-mediated neurodegeneration phenotype. Comparing the PS/E4H to E4H mice, we noted that 216 genes were upregulated while 40 genes were down regulated (Fig. 4A; Additional file 6: Table S2G). Genes corresponding to glial activation, such as Gfap, C1qa, C1qb and Mpeg1, were highly upregulated (Fig. 4A). Notable genes that have been reported to be part of the APOE-associated pathogenic signature in AD, Trem2 (↑3.7x, p = 1.46E−22) and Clec7A (↑36.2x, p = 3.1E−25), were upregulated in PS/E4H mice compared to E4H mice. APOE, Trem1 and Treml1 remained unchanged, but Treml2 was upregulated (↑4.7x, p = 4.22E−05). We also examined the cell types most perturbed in these mice. There was a slight decrease in neuronal gene expression (p = 0.052) and robust increase in glial genes (astrocyte, p = 0.0081; microglia, p < 0.0001; oligodendrocyte, p = 0.001) (Fig. 4B). We further examined whether a specific AD-associated astrocyte profile could be correlated with the early death phenotype. We found that pan-reactive astrocyte profile was upregulated (p = 0.0053), consistent with Gfap upregulation, but neither A1 or A2 astrocyte profiles were specifically induced. Among other aging/AD-associated gene expression signatures identified from human brains, we could identify plaque induced gene (PIG, p = 0.0342) and the H2M (p < 0.0001) homeostatic microglia signatures in the PS/E4H mice (Fig. 4B). The upward trends in microglia neurodegenerative signatures (MGnD and DAM) did not survive multiple testing correction (Fig. 4B). We performed pathway analysis using enriched gene sets. Using the GO and KEGG databases, we identified immune processes as the main pathways emerging in the PS/E4H vs E4H comparison. Examples of enriched pathways in PS/E4H mice include Il-6 and Ifn-γ signaling, antigen processing and presentation, Toll-like receptor and chemokine signaling pathway (Fig. 4C, D).
Compared to age-matched E3H mice, the PS/E3H mice showed 155 transcripts upregulated and 22 transcripts downregulated (Fig. 5A, Additional file 6: Table S2H). Mostly immune system mediators, such as Itgax, Fcrls, Lag3 and Gfap were upregulated in the PS/E3H mice relative to age-matched E3H mice (Fig. 5A). Similar to changes observed in PS/E4H cohorts, Trem2 (↑2.8x, p = 6.32E−13) and Clec7A (↑19.5x, p = 7.4E−18) were also upregulated in PS/E3H mice compared to E3H mice, though the strength of association (based on fold change and p values) was lower in the in PS/E3H mice. APOE, Trem1 and Treml1 levels also remained similar between PS/E3H and E3H mice, while Treml2 was upregulated (↑4.1x, p = 8.7E−04). The cell types that were most perturbed in the PS/E3H mice were microglia (p = 0.0025) and oligodendrocytes (p = 0.0102) (Fig. 5B). Similar to PS/E4H vs E4H comparison, we observed global changes in astrocytic gene expression (p = 0.0295) in PS/E3H mice, though the A1/A2 markers were not changed (Fig. 5B). Among the aging/AD-associated gene expression profiles examined in the PS/E3H vs E3H comparison, only the H2M signature was upregulated (p < 0.0001), with no changes in ARM, DAM, MGnD and PIG (Fig. 5B).
As there was a distinct difference in lifespan of the PS/E4H and PS/E3H mice, we were curious whether DEG analysis between these two cohorts could provide insights into tau-mediated pathological effects in E4H vs E3H background. A total of 45 genes were significantly altered, with 14 genes upregulated in PS/E4H mice (Additional file 5: Fig. S5C; Additional file 6: Table S2E). Analysis of the enriched gene sets revealed gene ontogeny and pathways related to pyrimidine metabolism, prolactin signaling, histone demethylation and immunoglobulin pathways most enriched in PS/E4H relative to PS/E3H mice (Additional file 5: Fig. S5I–J). Most of these pathway changes involved only a few genes that were enriched in the comparison analysis, ranging from 1 to 3 genes and the strength of the association (p values and gene numbers enriched in each pathway) were also weak. While our expectation was that the pathogenic interaction of APOE-TREM2 would be exacerbated in PS/E4H mice, surprisingly, we did not note any differential regulation (measured by fold change and p-value of significance) of Trem family members in the PS/E4H mice compared to PS/E3H mice (Additional file 5: Fig. S5K).
Tau transgene induced DEG in APOE2 homozygous PS19 mice
Relative to age-matched E2H mice, the PS/E2H showed 82 upregulated and 11 down regulated genes (Fig. 6, Additional file 6: Table S2I). Genes that were most perturbed were microglia or astrocyte-specific, such as C1qa, Grn, Itgax and Gfap (Fig. 6A), reiterating the commonalities of tau overexpression paradigm observed previously in PS/E3H and PS/E4H mice (Fig. 4 and 5). PS/E2H mice showed upregulated transcripts of Trem2 (↑2.2x, p = 3.7E−07), Clec7a (↑16.6x, p = 2.4E−15) and Treml2 (↑2.9x, p = 0.018), whereas APOE, Trem1 and Treml1 were unchanged. We observed only modest upregulation in microglial cell profile (p = 0.019), with other cell types (neurons, astrocytes and oligodendrocytes) remaining unaffected (Fig. 6B). The pan-reactive astrocyte gene expression profile (p = 0.0038), but not the A1/A2, was also identified (Fig. 6B). Among the aging/AD-associated transcriptional signatures, we found robust upregulation of H2M (p < 0.001) and modestly increased PIG (p = 0.0533) signature (Fig. 6B). Other AD-associated profiles, such as ARM, DAM or MGnD were not affected (Fig. 6B). Pathway analysis revealed immune system processes and defense response being enriched overall, with specific involvement of antigen processing, complement cascades, and NOD-like receptor signaling cascades in PS/E2H mice (Fig. 6C, D).
APOE-based differences in neuropathology in human AD cases
In order to understand the relationship between tau phosphorylation and APOE genotype, we next turned our attention to human AD cases that were stratified for APOE4 or non-APOE4 genotypes (Fig. 7; Additional file 6: Table S1). All of the AD cases had a primary neuropathological diagnosis of AD . Some of the AD patients also exhibited co-pathology in form of α-synuclein positive Lewy bodies alongside canonical phospho-tau and amyloid β (Aβ) deposits (Additional file 6: Table S1). These patients were classified as having ‘mixed’ pathologies (open diamonds, Fig. 7). The age-matched healthy controls with no history of dementia were all non-E4 genotype (Additional file 6: Table S1). Immunostaining with pan Aβ and AT8 antibodies demonstrated that AD patients of both E4 and non-E4 genotypes showed the expected robust patterns of Aβ deposits and intracellular tau tangles relative to age-matched controls (****p < 0.0001 for Aβ; **p < 0.01 for phospho-tau) (Fig. 7A, B). We did not observe any obvious differences in the Aβ and tau burden between the E4 and non-E4 individuals within the neuropathologically-characterized AD group (Fig. 7A, B). Interestingly, we found that AD patients with E4 allele had slightly higher α-synuclein burden relative to non-E4 AD patients (p = 0.030) and controls (p = 0.0082) (Fig. 7C). We then wanted to investigate whether there was a correlation between tau and Aβ or tau and α-synuclein based on APOE4 genotype within AD patients. Phospho-tau burden was correlated to Aβ burden in the E4 AD group (p = 0.019) but not in the non-E4 AD group (Fig. 7D). On the other hand, tau and α-synuclein burden showed no correlation in either non-E4 or E4 AD cohorts (Fig. 7E).
Here, we report on life expectancy and end-stage tau pathology in tau transgenic mice with different doses of the three different human APOE alleles. In aging cohorts of PS19 mice that reached end-of-life hindlimb paralysis, APOE3 mice lived the longest, whereas both APOE4 and APOE2 homozygotes showed significantly shorter lifespan. While the APOE4 and APOE3 homozygotes had equivalent levels of phospho-tau burden (AT8, AT270, AT180 and PHF1) at end-stage, the APOE4 mice had significantly reduced lifespan compared to APOE3 homozygotes, suggesting that APOE3 conferred some resiliency to the degenerative phenotype inherent in PS19 mice. When we compared PS19 mice that are heterozygous or homozygous for human APOE4 or APOE2 (i.e., mouse carrying one human and one mouse allele), we observed that human APOE homozygous mice generally developed paralysis sooner than heterozygous counterparts. Thus, the presence of one copy of mouse Apoe modified human APOE4- and human APOE2-dependent effects on lifespan of PS19 mice. This finding suggests that either the effect of human APOE is dose-dependent, or that mouse Apoe could have modulatory effects on tau-associated neurodegeneration relative to human APOE, or that even subtle differences in mouse genetic background could have a regulatory effect on tau-APOE interactions. Overall, our data shows that presence of two alleles of human APOE3 substantially extends lifespan of PS19 mice relative to mice with heterozygous or homozygous for human APOE4, human APOE2 or mouse Apoe.
Previous studies have highlighted how APOE genotype affected tauopathy in tau transgenic animals , in viral-mediated tau overexpression cohorts  and in K18-tau seed injected tau transgenic animals . Collectively, these data on APOE genotype association with tau neuropathology seem to be disparate between the studies. The first report on PS19 transgenic animals showed that APOE4 exacerbated tau burden and neurodegeneration  while the second study provides evidence for a highly pathologic role of APOE2 on tau pathology . Additionally, in the K18-tau aggregate seeding model, induction of AT8-phospho-tau was exacerbated in APOE3 mice, but the levels of early tangle pathology did not seem to be APOE genotype-dependent . In this current study, we establish the relationship of motor neuron degeneration (leading to paralysis) and tau neuropathology in PS19 transgenic mice carrying different human APOE alleles. Our data generally supports the observations of Shi and colleagues who reported that APOE4 exacerbates tau-mediated neurodegeneration in PS19 mice . However, there are several differences between our current observations and this previous study. Namely, we observed equivalent phospho-tau burden in brains of PS/E4H and PS/E3H mice at end stage using several different phospho-tau antibodies, whereas the former study reported focused on AT8-tau profiles associated with APOE genotype around 9.5 months of age. A caveat to our study is that it is possible we may have missed a critical window of APOE genotype-associated differential tau pathology as we examined mice at physiological end-of-life when our cohorts are expected to reach a ‘ceiling effect’ in terms of phospho-tau burden concurrent with paralysis. Additionally, our colony of PS/E4H mice did not display consistent ventricular enlargement relative to other APOE genotypes. It is possible that our modest sample size precluded us from reaching statistically significant data on APOE genotype-related cortical atrophy due to the high degree of inter-animal variability in ventricle size within our colony. Other studies have noted inherent phenotypic variability PS19 mice . Additionally, variability in phenotypes between mouse colonies maintained in different facilities have been attributed to background strain effects  and microbiome variabilities unique to different vivarium conditions . Importantly, the accelerated motor neuron degeneration phenotype leading to paralysis in PS/E4 mice relative to PS/E3 mice mirrors the cortical degeneration phenotype reported by Shi et al. in their PS19 mice humanized for APOE4 .
We also observed an unexpected negative association of PS19 mice lifespan with APOE2 homozygosity relative to APOE3, in spite of modest reductions in selective soluble phospho-tau species (AT8 and AT100) in the PS/E2H mice. Notably, the immune profile determined from transcriptomic changes in PS/E2H mice were similar to PS/E3H and PS/E4H mice. It is possible that the early death phenotype in PS/E2H mice could be related to pathways independent of somatodendritic tau accumulation. In humans, APOE2 is associated with higher longevity, lower Aβ and reduced risk of AD . At least in mice, the biologically beneficial effects of APOE2 on brain health has been shown to be related to cholesterol handling in the CNS  as well as reduced amyloid precursor protein (APP) transcription  relative to the reference APOE3 allele. However, during aging, APOE2 homozygous mice, by themselves, show hyperlipoproteinemia and spontaneous atherosclerosis , consistent with phenome-wide association studies revealing harmful effects of APOE2 homozygosity on peripheral vascular homeostasis in humans . Together, our data are consistent with the idea that APOE3 confers some level of systemic resilience relative to APOE2.
AD patients are characterized by the presence of Aβ and tau pathologies, in addition to some cases showing abundant mixed pathologies such as α-synuclein and TDP-43 inclusions . This does not allow for easy interpretation of APOE genotype-tau association data or establishing the individual contribution of a specific APOE genotype on Aβ, α-synuclein and tau pathogenesis separately. Using a cohort of AD patient brains that have undergone extensive post-mortem neuropathology analysis, we show that tau and Aβ levels are individually upregulated in these patients, irrespective of APOE genotype, though interestingly, α-synuclein pathology is slightly more predominant in APOE4 + individuals. Further analysis in this cohort of human AD patients appear to indicate a strong association between phospho-tau and Aβ levels, but not between phospho-tau and α-synuclein levels, in APOE4 + individuals. This implies a pathologic synergy between tau and Aβ directed by the presence of APOE4 genotype. The direct association between APOE4 and increased Aβ deposition has been corroborated over the last three decades of work in humans and mice [1, 37,38,39], though the direct link between APOE4 and tau tangles remain controversial [7, 40,41,42]. On the other hand, there is also considerable evidence that development and progression of tau tangle pathology is directly related to the extent and distribution of Aβ in AD patients and rodent models [43,44,45,46,47]. Thus, broadly our data is in agreement with the theory that higher levels of tau correspond to APOE4 + individuals that also have higher levels of Aβ [40, 48]. However, our data does not agree with neuropathology data derived from several cohorts that showed a direct association between APOE4 and tau neuroimaging measures in specific brain areas . Further experiments are needed to experimentally establish if specific APOE genotypes are directly associated with tau and tau-associated neurodegeneration when other AD-typical pathologies such as Aβ are also present.
Our RNAseq data indicates that there are only modest overall differences between the three parental lines of APOE homozygous mice, especially between APOE4 and APOE3 mice where the major differences centered on fatty acid metabolism and cytosolic DNA sensing pathway (inflammasome and Type 1 IFN response). Additionally, APOE4 and APOE2 comparisons showed differential regulation of immune signaling pathways relative to APOE3. Notably, Serpina3n is upregulated in APOE4 mice relative to both APOE3 and APOE2 mice, which is a key mediator of reactive astrocytosis . Erdr1, a stress-induced cell survival factor, was downregulated in APOE4 mice relative to APOE3 and APOE2 mice while Eps8l1, a tauopathy-related Bin1-responsive gene, was downregulated in APOE3 mice relative to APOE4 and APOE2 mice [51, 52]. Interestingly, the presence of tau transgene led to the emergence of overwhelmingly robust immune transcriptome, irrespective of the APOE genotype. Indeed, when we analyzed the differential transcriptome in PS19 mice carrying APOE4 or APOE2 or APOE3 relative to their corresponding parental APOE lines, we found that immune response and stress response pathways were commonly upregulated in PS19 mice irrespective of the APOE genotype. We also found that tau overexpression and accompanying tau neuropathology increased transcription of Trem family members in these PS19 mice, consistent with previous evidence of a pathologic relationship between Trem2 and Apoe , though there was no clear APOE genotype-dependent outcome noticed in our study. Originally these Apoe/Trem2-associated DAM microglial phenotypes were characterized in aging mouse brains associated with Aβ- or demyelinated-associated neurodegeneration . Subsequent work showed that loss of Trem2 worsened tau pathology, microgliosis and neurodegeneration in APOE4 mice  but not in mice carrying mouse Apoe . Together, these previous studies suggested that there could be a pathologic relationship between APOE4, tau and Trem2-microglia function, though there is a lack of consensus on underlying mechanisms. In our Rnaseq analyses on paralyzed cohorts, we did not observe a clear differential upregulation of the Trem2-associated DAM microglial profile in any particular APOE genotypes. On the other hand, our transcriptome studies showed that the homeostatic H2M signature was upregulated in all cohorts of PS19 mice, independent of APOE genotype. It is possible that the microglia in our PS19xAPOE cohorts are persisting in the homeostatic phase, pausing prior to maturing into DAM or other degenerative phenotype in the absence of additional insults. This could possibly suggest that the PS/E4H mice could require additional triggers, such as Tgfβ or other AD-related factors, for gradual differentiation into progressively dystrophic Apoe/Trem2-associated DAM state . We also noted that the APOE4 mice uniquely showed modest upregulation of another AD-associated profile, the PIG, that was identified from APP/Aβ model and human AD brains and is related to neurodegenerative stimulus . Future studies could focus on relating the relative abundance of these transcriptional signatures to APOE genotype-dependent biological outcomes.
In conclusion, our studies demonstrate that relative to human APOE2 and APOE4, the APOE3 allele delays paralysis in PS19 mice. At end-stage, the burden of phospho-tau pathology was similar among the three APOE genotypes, with some select phospho-tau epitopes being under-represented in APOE2 mice. Robust immune changes in the transcriptome was noted in PS19 mice, irrespective of the APOE genotype. Our findings indicate that APOE genotype plays a regulatory role in tauopathy, though the exact pathogenic mechanism remains elusive.
Availability of data and materials
All data and code are available from the corresponding author upon reasonable request following publication. The RNAseq data is available from GEO database (GSE226093).
Amyloid precursor protein
Damage associated microglia
Differential expression of genes
APOE2 targeted replacement mouse
APOE3 targeted replacement mouse
APOE4 targeted replacement mouse
Fetal bovine serum
Formalin fixed, paraffin embedded
Fragments per kilobase of exon per million mapped fragments
Glial fibrillary acidic protein
Neurodegeneration associated microglia
Normalized enrichment score
- NOM p- value:
Nominal p value
Plaque induced genes
PS19 mice Homozygous for APOE2
PS19 mice Homozygous for APOE3
PS19 mice Homozygous for APOE4
PS19 mice heterozygous for APOE2
PS19 mice heterozygous for APOE3
PS19 mice heterozygous for APOE4
Progressive supranuclear palsy
Relative fluorescent units
Single nucleotide polymorphism
Yamazaki Y et al (2019) Apolipoprotein E and Alzheimer disease: pathobiology and targeting strategies. Nat Rev Neurol 15(9):501–518
Long JM, Holtzman DM (2019) Alzheimer disease: an update on pathobiology and treatment strategies. Cell 179(2):312–339
Reiman EM et al (2020) Exceptionally low likelihood of Alzheimer’s dementia in APOE2 homozygotes from a 5,000-person neuropathological study. Nat Commun 11(1):667
Raber J et al (2000) Apolipoprotein E and cognitive performance. Nature 404(6776):352–354
Lumsden AL et al (2020) Apolipoprotein E (APOE) genotype-associated disease risks: a phenome-wide, registry-based, case-control study utilising the UK Biobank. EBioMedicine 59:102954
Parhizkar S, Holtzman DM (2022) APOE mediated neuroinflammation and neurodegeneration in Alzheimer’s disease. Semin Immunol 59:101594
Shi Y et al (2017) ApoE4 markedly exacerbates tau-mediated neurodegeneration in a mouse model of tauopathy. Nature 549(7673):523–527
Shi Y et al (2019) Microglia drive APOE-dependent neurodegeneration in a tauopathy mouse model. J Exp Med 216(11):2546–2561
Wang C et al (2021) Selective removal of astrocytic APOE4 strongly protects against tau-mediated neurodegeneration and decreases synaptic phagocytosis by microglia. Neuron 109(10):1657–1674
Zhao N et al (2018) APOE epsilon2 is associated with increased tau pathology in primary tauopathy. Nat Commun 9(1):4388
Tristan Williams AJR, Ruiz AM, Vo Q, Tsering W, Xu G, McFarland K, Giasson BI, Sullivan P, Borchelt DR, Chakrabarty P (2022) Impact of APOE genotype on prion-type propagation of tauopathy. Acta Neuropathol Commun 10(1):57
Yoshiyama Y et al (2007) Synapse loss and microglial activation precede tangles in a P301S tauopathy mouse model. Neuron 53(3):337–351
Knouff C et al (1999) Apo E structure determines VLDL clearance and atherosclerosis risk in mice. J Clin Invest 103(11):1579–1586
Sullivan PM et al (1998) Type III hyperlipoproteinemia and spontaneous atherosclerosis in mice resulting from gene replacement of mouse Apoe with human Apoe*2. J Clin Invest 102(1):130–135
Sullivan PM et al (1997) Targeted replacement of the mouse apolipoprotein E gene with the common human APOE3 allele enhances diet-induced hypercholesterolemia and atherosclerosis. J Biol Chem 272(29):17972–17980
Koller EJ et al (2019) Combining P301L and S320F tau variants produces a novel accelerated model of tauopathy. Hum Mol Genet 28(19):3255–3269
Hyman BT et al (2012) National Institute on Aging-Alzheimer’s Association guidelines for the neuropathologic assessment of Alzheimer’s disease. Alzheimers Dement 8(1):1–13
Bankhead P et al (2017) QuPath: open source software for digital pathology image analysis. Sci Rep 7(1):16878
Martin M (2011) Cutadapt removes adapter sequences from high-throughput sequencing reads. Embnet J 17:1
Kim D, Langmead B, Salzberg SL (2015) HISAT: a fast spliced aligner with low memory requirements. Nat Methods 12(4):357–360
Pertea M et al (2015) StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol 33(3):290–295
Love MI, Huber W, Anders S (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15(12):550
Robinson MD, McCarthy DJ, Smyth GK (2010) edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26(1):139–140
Zhang Y et al (2014) An RNA-sequencing transcriptome and splicing database of glia, neurons, and vascular cells of the cerebral cortex. J Neurosci 34(36):11929–11947
Chen WT et al (2020) Spatial transcriptomics and in situ sequencing to study Alzheimer’s disease. Cell 182(4):976–991
Friedman BA et al (2018) Diverse brain myeloid expression profiles reveal distinct microglial activation states and aspects of Alzheimer’s disease not evident in mouse models. Cell Rep 22(3):832–847
Hammond TR et al (2019) Single-cell RNA sequencing of microglia throughout the mouse lifespan and in the injured brain reveals complex cell-state changes. Immunity 50(1):253–271
Keren-Shaul H et al (2017) A unique microglia type associated with restricting development of Alzheimer’s disease. Cell 169(7):1276–1290
Krasemann S et al (2017) The TREM2-APOE pathway drives the transcriptional phenotype of dysfunctional microglia in neurodegenerative diseases. Immunity 47(3):566–581
Sala Frigerio C et al (2019) The major risk factors for Alzheimer’s disease: age, sex, and genes modulate the microglia response to abeta plaques. Cell Rep 27(4):1293–1306
Woerman AL et al (2017) Kinetics of human mutant tau prion formation in the brains of 2 transgenic mouse lines. JAMA Neurol 1:1
Neuner SM et al (2019) Harnessing genetic complexity to enhance translatability of Alzheimer’s disease mouse models: a path toward precision medicine. Neuron 101(3):399–411
Burberry A et al (2020) C9orf72 suppresses systemic and neural inflammation induced by gut bacteria. Nature 582(7810):89–94
Shinohara M et al (2016) APOE2 eases cognitive decline during Aging: Clinical and preclinical evaluations. Ann Neurol 79(5):758–774
Huang YA et al (2017) ApoE2, ApoE3, and ApoE4 differentially stimulate APP transcription and abeta secretion. Cell 168(3):427–441
Jack CR Jr et al (2018) NIA-AA research framework: toward a biological definition of Alzheimer’s disease. Alzheimers Dement 14(4):535–562
Gottesman RF et al (2016) The ARIC-PET amyloid imaging study: brain amyloid differences by age, race, sex, and APOE. Neurology 87(5):473–480
Jansen WJ et al (2015) Prevalence of cerebral amyloid pathology in persons without dementia: a meta-analysis. JAMA 313(19):1924–1938
Sabbagh MN et al (2013) The influence of Apolipoprotein E genotype on regional pathology in Alzheimer’s disease. BMC Neurol 13:44
Hohman TJ et al (2018) Sex-specific association of apolipoprotein E with cerebrospinal fluid levels of tau. JAMA Neurol 75(8):989–998
Mattsson N et al (2018) Greater tau load and reduced cortical thickness in APOE epsilon4-negative Alzheimer’s disease: a cohort study. Alzheimers Res Ther 10(1):77
La Joie R et al (2021) Association of APOE4 and clinical variability in Alzheimer disease with the pattern of tau- and amyloid-PET. Neurology 96(5):e650–e661
Jack CR et al (2019) The bivariate distribution of amyloid-beta and tau: relationship with established neurocognitive clinical syndromes. Brain 142(10):3230–3242
Ramanan VK et al (2019) Association of apolipoprotein E varepsilon4, educational level, and sex with tau deposition and tau-mediated metabolic dysfunction in older adults. JAMA Netw Open 2(10):e1913909
Dincer A et al (2022) APOE epsilon4 genotype, amyloid-beta, and sex interact to predict tau in regions of high APOE mRNA expression. Sci Transl Med 14(671):l7646
Koller EJ et al (2022) Combinatorial model of amyloid beta and tau reveals synergy between amyloid deposits and tangle formation. Neuropathol Appl Neurobiol 48(2):e12779
Lee WJ et al (2022) Regional Abeta-tau interactions promote onset and acceleration of Alzheimer’s disease tau spreading. Neuron 110(12):1932–1943
Jack CR Jr et al (2017) Age-specific and sex-specific prevalence of cerebral beta-amyloidosis, tauopathy, and neurodegeneration in cognitively unimpaired individuals aged 50–95 years: a cross-sectional study. Lancet Neurol 16(6):435–444
Therriault J et al (2020) Association of apolipoprotein E epsilon4 with medial temporal tau independent of amyloid-beta. JAMA Neurol 77(4):470–479
Zamanian JL et al (2012) Genomic analysis of reactive astrogliosis. J Neurosci 32(18):6391–6410
Dormer P, Spitzer E, Moller W (2004) EDR is a stress-related survival factor from stroma and other tissues acting on early haematopoietic progenitors (E-Mix). Cytokine 27(2–3):47–57
Ponnusamy M et al (2023) Loss of forebrain BIN1 attenuates hippocampal pathology and neuroinflammation in a tauopathy model. Brain 146(4):1561–1579
Chen Y, Colonna M (2021) Microglia in Alzheimer’s disease at single-cell level: Are there common patterns in humans and mice? J Exp Med 218(9):1
Gratuze M et al (2023) TREM2-independent microgliosis promotes tau-mediated neurodegeneration in the presence of ApoE4. Neuron 111(2):202–219
Gratuze M et al (2020) Impact of TREM2R47H variant on tau pathology-induced gliosis and neurodegeneration. J Clin Invest 130(9):4954–4968
We thank the University of Florida Neuromedicine Human Brain and Tissue Bank for sharing human brain samples (funded in part by NIA P30AG066506). We thank Dr. Todd E. Golde and Dr. Benoit I. Giasson for sharing Ab5 and 94-3A10 antibodies respectively. We posthumously thank Dr. Peter Davies for sharing the PHF1 antibody.
The work was supported by NIA RF1AG057933 (DRB, PC), R01AG055798 (PC) and NIA T32 AG061892 (TW).
Ethics approval and consent to participate
All rodent studies were approved by the University of Florida IACUC. All human autopsy materials were handled with the approval of the Institutional Review Board (IRB) at University of Florida. All patients provided written informed consent.
Consent for publication
All authors have approved the manuscript and agree with its submission.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Analysis of aging cohorts of PS19xAPOE mice. Representative images of hematoxylin & eosin-stained brains and volumes of ventricles at two locations (bregma +2.6mm denoted as ‘lateral’, A, and bregma +3.2mm denoted as ‘medial’, B) in PS19 mice homozygous for APOE2, APOE3 or APOE4 shown. Male mice indicated in blue; female mice indicated in pink. N = 6 (PS/E2H), n = 6 (PS/E3H), n = 8 (PS/E4H).
. Analysis of phospho-tau epitopes in paralyzed PS19 mice homozygous for APOE2, APOE3 or APOE4. A-B. Representative images from hippocampus and corresponding quantitation of whole brain burden of phospho-tau in PS19 mice homozygous for APOE2 (triangle), APOE3 (circle) or APOE4 (diamond). 3 different phospho-tau epitopes were tested: AT100, AT270 and AT180. 1-way Anova. N = 10 (PS/E2H), n = 8 (PS/E3H), n = 9 (PS/E4H). C-G. Biochemical analysis of phospho-tau using various antibodies in detergent soluble (S1) and detergent-insoluble (S3 pellet) fractionated brains of paralyzed mice. Band intensities of phospho-tau in S1 fraction have been normalized to GAPDH, while band intensities of phospho-tau in S3 fraction (if detectable) have been normalized to total tau. Numerals on the left side of each blot denote molecular weight standards in kDa. As no phospho-tau was detected for AT100 (C), AT180 (E) and AT8 (G), S3-phospho-tau values have not been plotted. S1- and S3- phospho-tau levels have been presented for AT270 (D) and PHF1 (F) epitopes. Total tau levels in S3 fraction corresponding to the AT8-stained blot (G) is shown in Fig. 1D. The loading order has been denoted by consecutive numbering and the key is presented (H). 1-way Anova; *p<0.05. n = 3 mice/genotype.
. Inflammation markers in paralyzed PS19 mice homozygous for APOE2, APOE3 or APOE4. Immunoblots representing Iba-1 (A) and GFAP (B) levels from the whole forebrain and corresponding quantitation of protein bands are presented for PS19 mice carrying mouse ApoE, APOE mice (with no tau) or PS19 mice carrying human APOE alleles. Individual protein bands were normalized to actin and data from 3 different blots were averaged for the graph. N = 3 mice/group. 1-way Anova, **p<0.01, *p<0.05.
. Lifespan analysis and neuropathological characterization of PS19 mice heterozygous for APOE. PS/E2h, PS/E3h and PS/E4h mice were aged to paralysis as indicated (A). Median age to paralysis of these APOE heterozygous ‘h’ mice are compared to APOE-genotype matched homozygous ‘H’ mice - PS/E2H, PS/E3H and PS/E4H mice (shown in Fig. 1) (B). Log-rank (Mantel-Cox) test. Significant p values are in bold. n = 34 (PS/E2h), n = 31 (PS/E3h), n = 43 (PS/E4h). Paralyzed PS/E2h, PS/E3h, PS/E4h mice were analyzed for phospho-tau, microgliosis and astrocytosis levels using AT8 (C), Iba-1 (D) and GFAP (E) antibodies respectively. Quantification of AT8, Iba-1, and GFAP immunostaining is presented as % immunoreactivity in the cortex and hippocampus underneath corresponding antibody-stained panels. n = 12 (PS/E2h), n = 11 (PS/E3h), n = 26 (PS/E4h). 1-way ANOVA, *p<0.05. Scale bar: 100 µm.
. Comparative RNAseq analysis of PS/E2H, PS/E3H and PSE4H mice. Heat map (A) and individual volcano plots (B-D) of differentially expressed genes in PS/E2H, PS/E3H and PSE4H mice. Q value (qval) indicates false discovery rate. N = 4 mice/group. Asterisks in volcano plot refer to predicted genes. Commonly enriched gene sets (upregulated or downregulated) were identified from a comparison of PS/E2H vs PS/E3H and E2H vs E3H to enumerate genes that are driven by the variation in APOE allele. Total number of significantly altered genes are indicated by numerals within the Venn diagram. The comonly altered genes (in the overlapping Venn diagram) were used to map onto KEGG pathways for this comparison (E). Similar KEGG pathway mapping was done for PS/E3H vs E3H and PS/E2H vs E2H (F), PS/E4H vs PS/E3H and E4H vs E3H (G), PS/E4H vs E4H and PS/E3H vs E3H (H) using overlapping gene sets from the Venn diagram. Enriched gene sets from comparison of PS/E4H vs PS/E3H is represented as KEGG pathways (I) and GO pathways (J). The bubble plots (I, J) are colored by p-value and sized by number of genes in enriched gene sets. Tabulation of gene expression levels of Trem family genes (K) from PS/E4H vs PS/E3H analysis. BaseMean, means of normalized counts of all samples; FC, Fold change; lFcSE, log of Fold change Standard Error; stat, Wald statistics. N = 3-4 mice/group.
Detailed information of human brains used in the study. Individually characterized brains were obtained from the University of Florida Neuromedicine Human Brain and Tissue Bank. Experimental Group denotes the designation used in this study (Fig. 7). Dx1, primary diagnosis based on neuropathology; Dx2, Dx3, Dx4, secondary diagnoses based on neuropathology; Thal phase, burden of immunostained amyloid deposits in cortical and subcortical area; CERAD score, neuritic plaque frequency; AD, Alzheimer's disease; ARTAG, Aging-related tau astrogliopathy; CAA, cerebral amyloid angiopathy; CVD, cardiovascular disease; LATE-NC, limbic-predominant age-related TDP-43 encephalopathy neuropathological change; LBD, Lewy body dementia; PART, primary age-related tauopathy. Table S2. Transcriptome analysis of forebrains of PS19 mice homozygous for APOE2, APOE3 or APOE4 and corresponding APOE TR lines. RNAseq data has been deposited in GEO database. A selection of genes that are altered in each experimental group comparisons are presented. The corresponding tables show the DEG analysis from the following comparisons: Table S2A, E2H vs E3H; Table S2B, E4H vs E3H; Table S2C, E4H vs E2H; Table S2D, PS/E2H vs PS/E3H; Table S2E, PS/E4H vs PS/E3H; Table S2F, PS/E4H vs PS/E2H; Table S2G, PS/E4H vs E4H; Table S2H, PS/E3H vs E3H; Table S2I, PS/E2H vs E2H.
About this article
Cite this article
Williams, T., Bathe, T., Vo, Q. et al. Humanized APOE genotypes influence lifespan independently of tau aggregation in the P301S mouse model of tauopathy. acta neuropathol commun 11, 99 (2023). https://doi.org/10.1186/s40478-023-01581-2