APOE4 exacerbates synapse loss and neurodegeneration in Alzheimer’s disease patient iPSC-derived cerebral organoids

Zhao, Jing; Fu, Yuan; Yamazaki, Yu; Ren, Yingxue; Davis, Mary D.; Liu, Chia-Chen; Lu, Wenyan; Wang, Xue; Chen, Kai; Cherukuri, Yesesri; Jia, Lin; Martens, Yuka A.; Job, Lucy; Shue, Francis; Nguyen, Thanh Thanh; Younkin, Steven G.; Graff-Radford, Neill R.; Wszolek, Zbigniew K.; Brafman, David A.; Asmann, Yan W.; Ertekin-Taner, Nilüfer; Kanekiyo, Takahisa; Bu, Guojun

doi:10.1038/s41467-020-19264-0

Download PDF

Article
Open access
Published: 02 November 2020

APOE4 exacerbates synapse loss and neurodegeneration in Alzheimer’s disease patient iPSC-derived cerebral organoids

Jing Zhao^1,2,
Yuan Fu¹,
Yu Yamazaki ORCID: orcid.org/0000-0002-0806-4960¹,
Yingxue Ren³,
Mary D. Davis^1,2,
Chia-Chen Liu ORCID: orcid.org/0000-0001-7203-3926¹,
Wenyan Lu^1,2,
Xue Wang³,
Kai Chen¹,
Yesesri Cherukuri³,
Lin Jia¹,
Yuka A. Martens ORCID: orcid.org/0000-0002-5194-4269¹,
Lucy Job^1,2,
Francis Shue ORCID: orcid.org/0000-0003-1659-1952¹,
Thanh Thanh Nguyen ORCID: orcid.org/0000-0002-2281-6540¹,
Steven G. Younkin¹,
Neill R. Graff-Radford⁴,
Zbigniew K. Wszolek⁴,
David A. Brafman⁵,
Yan W. Asmann³,
Nilüfer Ertekin-Taner^1,4,
Takahisa Kanekiyo^1,2 &
…
Guojun Bu ORCID: orcid.org/0000-0003-3491-1016^1,2

Nature Communications volume 11, Article number: 5540 (2020) Cite this article

30k Accesses
152 Citations
66 Altmetric
Metrics details

Subjects

An Author Correction to this article was published on 05 May 2021

This article has been updated

Abstract

APOE4 is the strongest genetic risk factor associated with late-onset Alzheimer’s disease (AD). To address the underlying mechanism, we develop cerebral organoid models using induced pluripotent stem cells (iPSCs) with APOE ε3/ε3 or ε4/ε4 genotype from individuals with either normal cognition or AD dementia. Cerebral organoids from AD patients carrying APOE ε4/ε4 show greater apoptosis and decreased synaptic integrity. While AD patient-derived cerebral organoids have increased levels of Aβ and phosphorylated tau compared to healthy subject-derived cerebral organoids, APOE4 exacerbates tau pathology in both healthy subject-derived and AD patient-derived organoids. Transcriptomics analysis by RNA-sequencing reveals that cerebral organoids from AD patients are associated with an enhancement of stress granules and disrupted RNA metabolism. Importantly, isogenic conversion of APOE4 to APOE3 attenuates the APOE4-related phenotypes in cerebral organoids from AD patients. Together, our study using human iPSC-organoids recapitulates APOE4-related phenotypes and suggests APOE4-related degenerative pathways contributing to AD pathogenesis.

APOE2 mitigates disease-related phenotypes in an isogenic hiPSC-based model of Alzheimer’s disease

Article Open access 09 April 2021

Nicholas Brookhouser, Sreedevi Raman, … David A. Brafman

Neuronal APOE4 removal protects against tau-mediated gliosis, neurodegeneration and myelin deficits

Article Open access 20 February 2023

Nicole Koutsodendris, Jessica Blumenfeld, … Yadong Huang

Cell type-specific roles of APOE4 in Alzheimer disease

Article 08 January 2024

Jessica Blumenfeld, Oscar Yip, … Yadong Huang

Introduction

Alzheimer’s disease (AD) characterized by progressive neurodegeneration is the most common form of dementia¹. While formations of senile plaques and neurofibrillary tangles (NFTs) are two major neuropathological hallmarks in AD, these pathologies composed of amyloid-β peptides (Aβ) and phosphorylated tau often precede the onset of symptomatic dementia by decades^2,3.

The ε4 allele of APOE gene (APOE4) is the strongest genetic risk factor for AD among the three polymorphic alleles (APOE2, APOE3, and APOE4)^4,5. Apolipoprotein E (apoE) is produced primarily by astrocytes in the central nervous system as a carrier of cholesterol and other lipids to support membrane homeostasis, synaptic integrity, and injury repair^6,7. Increasing evidences from animal models and postmortem human brains have shown that APOE4 is associated with multiple aspects of AD pathogenesis⁸. In particular, mouse models carrying human APOE4 have accelerated Aβ seeding and suppressed Aβ clearance, as well as disturbed synaptic plasticity and blood–brain barrier integrity^7,9; however, whether these findings are translatable to humans is not clear. In addition, studies with human postmortem brains represent only the characteristics of end-stage disease^10,11,12. Thus, there is an urgent need to establish human-relevant models to define APOE4-related pathogenic pathways in AD.

The induced pluripotent stem cells (iPSCs) derived from human somatic cells with AD-linked mutations or risk alleles are promising in vitro models by recapitulating the earliest molecular and pathological changes in age-related disorders^{13,14,15,16,17}. Elevated levels of Aβ and phosphorylated tau, as well as increased cellular stress markers have been reported in iPSC-derived neurons from AD patients^14,15,16,18. For example, Aβ oligomers accumulate in iPSC-derived neurons and astrocytes in cells from patients with a familial amyloid precursor protein (APP) p.E693Δ mutation and sporadic AD, leading to endoplasmic reticulum (ER) and oxidative stress¹⁵. In addition to altered APP processing, an increase in levels of total and phosphorylated tau is observed in iPSC-derived neurons from AD patients carrying the APP p.V717I mutation or APOE4 risk allele^18,19. However, the relevance of the observations from two-dimensional (2-D) cell cultures to AD is questionable as AD pathology is intricate and involves diverse cell types, tissue structures, and cellular pathways. Thus, in this study, we utilize a large number of human iPSC lines from healthy subjects and AD patients carrying APOE ε3/ε3 or ε4/ε4 genotype, and investigate AD-related phenotypes using the iPSC-derived three-dimensional (3-D) cerebral organoid model system, which is highly reminiscent of human brain structure with diverse cell types^20,21,22. Here, we show that APOE4 aggravates neurodegeneration in iPSC-derived cerebral organoids from AD patients; however, its effects on tauopathy are significant in both healthy subject-derived and AD patient-derived cerebral organoids. Interestingly, the levels of Aβ are increased in AD organoids independent of APOE genotype. Importantly, isogenic conversion of APOE4 to APOE3 attenuates AD-related phenotypes in iPSC-derived cerebral organoids. These findings reveal APOE4-related pathways in 3-D models that are directly relevant to humans.

Results

Characterization of human iPSC-derived cerebral organoids

To study the effects of APOE genotype on AD-related pathways in a physiologically relevant environment, we generated 3-D cerebral organoid models using human iPSC lines from cognitively unimpaired individuals carrying APOE ε3/ε3 (Con-E3; N = 5), cognitively unimpaired individuals carrying APOE ε4/ε4 (Con-E4; N = 5), AD patients carrying APOE ε3/ε3 (AD-E3; N = 5) and AD patients carrying APOE ε4/ε4 (AD-E4; N = 5), collected from multiple sources (Supplementary Table 1). Among them, 5 iPSCs lines were generated de novo from fibroblasts by transfecting with three episomal vectors encoding five transcription factors (OCT3/4, SOX2, L-MYC, KLF4, and LIN28) and p53 shRNA²³. Expression of pluripotency stem cell-specific markers was confirmed by immunostainings for SSEA4, Nanog, and TRA-1-60 (Supplementary Fig. 1A). The pluripotency of the iPSC lines was also validated by their ability to differentiate into endodermal, mesodermal, and ectodermal origin cells upon immunostaining of Brachyury (Mesoderm marker), Sox17 (endoderm marker), and Nestin/Sox2 (ectoderm marker) (Supplementary Fig. 1B). The iPSC lines maintained a normal karyotype after reprogramming (Supplementary Fig. 1C). All other iPSC lines utilized in published studies^19,24,25,26 or from California Institute for Regenerative Medicine (CIRM) (FUJIFILM Cellular Dynamics, Inc.) have been fully validated and characterized previously.

To generate cerebral organoids from iPSCs with efficiency and reproducibility, we followed an optimized protocol developed by Stemcell Technologies. Cerebral organoid formation was initiated through an intermediate embryonic body (EB) formation step followed by expansion of neuroepithelia in a matrigel scaffold. On day 12, the iPSC-derived organoids were transferred to an orbital shaker in neuronal differentiation medium and maintained under rotary conditions for maturation (Fig. 1a). On week 4, cerebral organoids showed a predominantly dorsal forebrain region specification, containing fluid-filled ventricle-like structures aligned with Sox2-positive neural progenitors in a ventricular/subventricular-like zone (VZ/SVZ) and beta-tubulin III (Tuj1)-positive neuroblasts in an outer layer (Fig. 1b). A deep cortical layer marker Ctip2-positive neurons were detected as early as week 4 (Fig. 1c), whereas a superficial cortical layer marker Satb2-positive neurons emerged in the later stage at week 12 (Fig. 1d). These observations revealed the sequential emergence of different neuronal layers along with the differentiation, which is consistent with previous publication²⁰. Since apoE is mainly produced by astrocytes in the brain²⁷, we assessed the presence of astrocytes at different time points by immunostaining for glial fibrillary acidic protein (GFAP). We found that small clusters of GFAP-positive astrocytes started to emerge in some VZ area at week 4. GFAP-positive astrocytes showed an immature morphology with short processes, which were separated from surrounding Tuj1-positive neuronal cells (Fig. 1e). At week 12 of differentiation, GFAP-positive astrocytes increased in number and migrated within the neuronal layers, displaying typical mature astrocyte morphology with long processes (Fig. 1f). Furthermore, to determine the influence of the technical difference to the differentiation of cerebral organoids across different rounds, we compared the levels of GFAP and Tuj1 in the cerebral organoids from two rounds of differentiation. We found no significant differences in the expression of GFAP and Tuj1 at week 12 between the first and second round of experiments (Supplementary Fig. 9). These results indicate the successful development of iPSCs to cerebral organoids with cortical like organization composed of abundant mature neurons and astrocytes.

**Fig. 1: Generation and characterization of cerebral organoids from human iPSCs.**

APOE4 and AD status exacerbate neurodegeneration

With the establishment of the cerebral organoid culture system, we generated four groups of iPSC-derived cerebral organoids (Con-E3, Con-E4, AD-E3, and AD-E4). Cerebral organoid size was monitored and there were no evident differences observed among different groups at 2 and 12 weeks (Supplementary Fig. 2). To evaluate the effect of APOE4 and disease status on neuronal apoptosis/degeneration in the cerebral organoids, cleaved caspase-3 (CASP3) was analyzed by immunostaining at week 12 (Fig. 2a). To avoid the possible influences of necrosis observed in the core portion of organoids, the immunoreactivity was measured only in the surface neuronal layers. We found the increased cleaved CASP3 immunoreactivity in the AD-E4 organoids compared to other groups with an interactive effect between APOE4 and AD status (Fig. 2b). Consistently, western blotting also revealed higher cleaved CASP3/CASP3 ratio in AD-E4 organoids (Fig. 2c, d). These results indicate that APOE4 and AD status synergistically exacerbate apoptosis in late stage of organoid development. Presynaptic synaptophysin and postsynaptic PSD95 were decreased in AD organoids groups compared to healthy subject-derived cerebral organoids, whereas no significant APOE4 effect was observed (Fig. 2c, e, and f). In contrast, on week 4, reverse transcription-quantitative PCR (RT-qPCR) revealed that cerebral organoids from AD patients exhibited higher mRNA levels of mature neuronal markers including MAP2, CTIP2, and SATB2, but not GFAP, an astrocyte marker (Supplementary Fig. 3A–D). Synaptophysin and PSD95 were also upregulated in AD organoid groups by week 4 (Supplementary Fig. 3E–G). While synaptophysin levels were increased in AD-E4 organoids (Supplementary Fig. 3F), APOE4 did not influence PSD95 levels (Supplementary Fig. 3G). Together, these results suggest that organoid maturation and synaptic formations are accelerated in cerebral organoids from AD patients at the early stage.

**Fig. 2: *APOE4* enhances apoptosis and synaptic loss in cerebral organoids from AD patients.**

Increased Aβ amounts in cerebral organoids from AD-iPSCs

To investigate the impacts of APOE4 and disease status on Aβ accumulation and deposition, iPSC-derived cerebral organoids were utilized for the analyses at different time points (weeks 4, 8 and 12). Organoids were sequentially lysed in RIPA buffer and formic acid (FA), and subjected to the measurements for Aβ40 and Aβ42. Using ELISA, we detected Aβ at an earlier time point than that reported by Lin et al.²⁸, in which they assessed Aβ accumulation by western blotting. Higher levels of Aβ40 at week 8 (Supplementary Fig. 4A) and Aβ42 at weeks 4 and 8 (Supplementary Fig. 4B) were observed in the RIPA-soluble fraction from AD organoids as compared to those from healthy subject-derived cerebral organoids. The increases of Aβ40 (Fig. 3a), Aβ42 (Fig. 3b), and Aβ42/Aβ40 ratio (Fig. 3c) in AD organoid groups were magnified at week 12 independent of APOE4. Aβ40 and Aβ42 in detergent-insoluble FA fraction were undetectable at any time points. In addition, neither APOE4 nor disease status affected APP derivatives including sAPPα, sAPPβ, and CTF-β at week 12 (Fig. 3d–f), whereas APOE4 but not AD status was associated with lower levels of full-length APP (Fig. 3g, h). These results suggest that Aβ clearance mechanism rather than APP processing is altered in cerebral organoids from AD patients regardless of APOE4 status, resulting in enhanced Aβ accumulation and increased Aβ42/Aβ40 ratio.

**Fig. 3: Increased Aβ accumulation in iPSC-derived cerebral organoids from AD patients.**

Exacerbated tauopathy in cerebral organoids from APOE4-iPSCs

To examine the tau pathology in the healthy subject-derived and AD patient-derived cerebral organoids carrying different APOE genotypes, we stained cerebral organoids with anti-phosphorylated tau (p-tau) AT8 at week 12 and found more p-tau accumulation in organoids with APOE4 (Fig. 4a). Consistent with the results, western blotting with AT8 revealed that APOE4 and AD status were correlated with higher p-tau/tau ratio in the organoids (Fig. 4b–d). When RIPA lysates (Fig. 4e) and FA fractions (Fig. 4f) were subjected to ELISA, APOE4, and AD status were independently associated with p-tau upregulation in the organoids at week 12. While RIPA-soluble p-tau increased in all groups in time, APOE4 or AD patient-derived organoids showed higher p-tau levels as early as week 4 (Supplementary Fig. 4C). Together, these results indicate that APOE4 and AD status aggravate tau pathology in iPSC-derived organoids without interactive effects.

**Fig. 4: *APOE4* and AD status enhance p-tau levels in iPSC-derived cerebral organoids.**

When apoE levels were measured in the organoids by ELISA, higher soluble apoE levels in RIPA fraction was associated with APOE4 at week 12 (Fig. 5a); however, AD status rather than APOE4 was associated with increased insoluble apoE in the FA fraction at week 12 (Fig. 5b). Both APOE4 and AD status were associated with the increased soluble apoE levels at week 4 (Supplementary Fig. 4D). No significant differences in GFAP/Tuj1 ratio were observed among the four groups of cerebral organoids by immunostaining (Supplementary Fig. 5A–B), RT-qPCR (Supplementary Fig. 5C), and western blotting (Supplementary Fig. 5D), suggesting that APOE4 or AD status leads to the higher apoE level without increasing astrocyte population in the cerebral organoids. Among all groups, significant positive correlations were observed between p-tau and apoE levels in RIPA (Fig. 5c) and FA fractions (Fig. 5d). Although there were no evident correlations between Aβ and apoE levels (Fig. 5e, f), Aβ levels (Fig. 5g, h) were also positively correlated with p-tau levels in RIPA fraction. These results suggest that Aβ and apoE might associate with p-tau through different molecular pathways.

**Fig. 5: Positive correlation between apoE and p-tau levels in the iPSC-derived cerebral organoids.**

Altered transcriptomes in cerebral organoids from AD-iPSCs

To further address the impact of APOE4 and/or AD status on transcriptional profiles in iPSC-derived organoids, we performed RNA-sequencing (RNA-seq) at week 12. We identified 1331 and 717 differentially expressed genes (DEGs) between Con-E3 and Con-E4, and between AD-E3 and AD-E4, respectively, with an overlap of 302 genes (Supplementary Data 1, Supplementary Fig. 5A). Among the overlapped genes, 265 genes changed in the same direction. When comparing between healthy subject-derived and AD patient-derived cerebral organoids, we identified 1188 DEGs between Con-E3 and AD-E3, and 501 DEGs between Con-E4 and AD-E4, with an overlap of 317 genes 304 of which were changed in the same direction (Supplementary Data 1, Supplementary Fig. 6B). Furthermore, pathway analysis revealed “Amyloid proteins” as the top-ranked network affected by the disease status and APOE4 interaction (Supplementary Fig. 6C). The top-ranked pathways enriched by DEGs were “Actin filaments” in Con-E3 vs. Con-E4 comparison (Supplementary Fig. 6D), “Regulation of epithelial-to-mesenchymal transition” in AD-E3 vs. AD-E4 comparison (Fig. S6E), “Axonal guidance” in Con-E3 vs. AD-E3 comparison (Supplementary Fig. 6F), and “Endoplasmic reticulum stress pathway” in Con-E4 vs. AD-E4 comparison (Supplementary Fig. 6G). “Innate inflammatory response” also showed up in all of the top 10-ranked pathways, suggesting a potential involvement of inflammation in APOE4 effects and AD pathogenesis (Supplementary Fig. 6D–G). The expression of selected most significant changed genes (DND1, NR2F6, VPS9D1, CTH, and CDHR3) was further validated by RT-qPCR, the results of which were consistent with those from RNA-seq (Supplementary Fig. 6H–L).

To assess potential contributions of different brain cell types for these DEGs, we analyzed the RNA-seq data using CIBERSORT and CellCODE program. We found that genes detected in the cerebral organoids were mainly derived from neurons (median: 0.9231) and astrocytes (median: 0.06698), whereas oligodendrocytes compose negligible proportion of the cell population (Supplementary Fig. 7A). No significant differences in neuron and astrocyte proportion were observed among different groups (Supplementary Fig. 7B–C). As more DEGs were designated to neurons than those assigned to astrocytes in all comparisons (Supplementary Fig. 7D–G), those results suggest that neurons predominantly contributed to the phenotypes of cerebral organoids from different groups compared to astrocytes.

Weighted Gene Co-expression Network Analysis (WGCNA) identified two gene modules that were significantly associated with disease status: module magenta was positively correlated with AD, while module yellow was positively correlated with healthy subject-derived cerebral organoids (negatively correlated with AD) (Fig. 6a). Genes in the yellow module were enriched for RNA metabolism dysregulation (Fig. 6b), which included ERCC4, DGCR8, POLR3A, CLP1, HSPA4, PNO1, VPS18, RAD17, LCMT2, and RPUSD2 as top-ranked hub genes (Fig. 6c). We further validated expressions of selected genes by RT-qPCR and western blotting, and confirmed decreased mRNA and protein levels in AD-E4 organoids (Fig. 6d, f, j–l), Genes in the magenta module were enriched for those related to DNA and nucleosome metabolism pathways (Supplementary Fig. 8A). CHIC1, H2BC7, H2BC8, ASB3, H2BC18, H2BC21, H1-3, ARRDC3, and TBL1XR1-AS1 were identified as top-ranked hub genes (Supplementary Fig. 8B). We further confirmed the increased mRNA levels of CHIC1, ASB3, and ARRDC3 by RT-qPCR in AD organoids, in particular those with APOE4 (Supplementary Fig. 8C–E). Since some of the hub genes in the yellow module (PNO1, DGCR8, VPS18, and HSPA4) are closely related to lysosomal stress granules formation under cellular stress conditions^29,30, we immunostained the cerebral organoids for a stress granule marker Ras GTPase-activating protein-binding protein 1 (G3BP) in the cerebral organoids at week 12, and found an increase of G3BP-positive punctates in AD organoids (Fig. 6e). We further quantified G3BP, endosome marker protein early endosome antigen 1 (EEA1), and lysosome marker protein lysosomal-associated membrane protein 1 (LAMP1) by western blotting, and found the interactive effect between APOE4 and AD status; G3BP significantly increased in organoids from AD patients carrying APOE4 (Fig. 6f, g). No significant changes were observed in EEA1 (Fig. 6h) and LAMP1 levels (Fig. 6i). These results imply that the disruption of RNA metabolism in cerebral organoids from AD patients accelerates stress granule formation especially in the presence of APOE4.

**Fig. 6: Transcriptomics profiling of the iPSC-derived cerebral organoids by WGCNA.**

Conversion of APOE4 to APOE3 attenuates AD pathologies

Using a parental iPSC line with APOE ε4/ε4 from a sporadic AD patient (Par-E4/4) and the isogenic iPSC line with APOE ε3/ε3 (Iso-E3/3)¹⁹, we tested whether converting APOE4 to APOE3 could ameliorate APOE4-related phenotypes in iPSC-derived organoids. While no obvious size difference was observed between cerebral organoids from the parental line and isogenic line at week 12 (Fig. 7a), immunostaining showed the reduction of cleaved CASP3 immunoreactivity in the edge of cerebral organoids from Iso-E3/3 compared to Par-E4/4 organoids (Fig. 7b, Supplementary Fig. 10C). The levels of cleaved CASP3 in the center of the cerebral organoids at week 12 were measured to evaluate the possible influence of necrotic core. No significant differences in the center of cerebral organoids from the parental line and isogenic line were observed (Supplementary Fig. 10B). Correlation analysis showed no significant correlation between the levels of cleaved CASP3 at the center with those at the edge, suggesting the apoptosis at the edge of cerebral organoids is not driven by the necrotic core at the center (Supplementary Fig. 10D). Western blotting also showed a lower cleaved CASP3/CASP3 ratio in Iso-E3/3 organoids than Par-E4/4 organoids (Fig. 7c). In addition, Iso-E3/3 organoids showed the decreased Aβ40 (Fig. 7d) and Aβ42 (Fig. 7e) levels in the RIPA lysate, although the Aβ42/40 ratio (Fig. 7f) remained unaffected by APOE4 conversion. Whereas RIPA-soluble apoE levels were not affected (Fig. 7g), insoluble apoE in FA fraction was decreased in the Iso-E3/3 organoids (Fig. 7h). Immunostaining (Fig. 7i) and western blotting (Fig. 7j) showed that APOE4 conversion to APOE3 reduced p-tau levels and p-tau/tau ratios in both RIPA (Fig. 7k) and FA fractions (Fig. 7l) in Iso-E3/3 organoids compared with Par-E4/4 organoids, which were confirmed by ELISA. Converting APOE4 to APOE3 reduced the levels of G3BP in cerebral organoids when analyzed by western blotting (Fig. 7m). Taken together, these results indicate that APOE4-related phenotypes observed in AD patient-derived organoids can largely be reversed through genome editing to APOE3.

**Fig. 7: Isogenic conversion of *APOE4* to *APOE3* attenuates AD-related phenotypes in cerebral organoids.**

Discussion

While AD is neuropathologically diagnosed through postmortem assessments in patients with dementia³¹, recent research framework has established antemortem AD classification using biomarkers for Aβ deposition (A), pathologic tau (T), and neurodegeneration (N)³². Thus, to explore the complex AD pathogenesis and develop therapeutic interventions for the disease, the establishment of human models system recapitulating ATN phenotypes has become essential. Toward this, we have comprehensively investigated AD-related pathogenic pathways using iPSC-derived 3-D cerebral organoids from sporadic AD patients with or without APOE4. cerebral organoids not only simulate intrinsic spatial patterning, but also display acquisition of cell identity in a timed manner that closely mimics the temporal patterning with sequential neuronal layer formation, accompanied with matured astrocytes^21,22,33. Furthermore, the phenotypes related to aberrant extracellular protein aggregation can be recapitulated only in 3-D culture systems^28,34,35. Using iPSC-derived cerebral organoid system, we identified multiple APOE4- and/or AD disease status-dependent pathogenic pathways, revealing complex etiology associated with APOE4, the strongest genetic risk factor for late-onset AD.

Consistent with results from 2-D cultures of iPSC-derived neurons^13,18,28,36, we found the elevated Aβ40 and Aβ42 levels in RIPA-soluble fraction in AD patient-derived organoids regardless of APOE4 status. Since APP processing was not altered, AD brain organoids may have compromised machineries for enzymatic Aβ degradation and/or cellular Aβ clearance. Of note, a significant reduction of soluble Aβ levels was observed in the cerebral organoids from a sporadic AD patient upon isogenic conversion of APOE4 to APOE3. Thus, it is possible that APOE4 and other gene variants synergistically facilitate Aβ accumulation in AD organoids. Notably, insoluble Aβ and amyloid plaque pathology was undetectable by ELISA and immunostaining in our organoid models from sporadic AD patients at week 12. While Aβ plaque deposition is detected 2–4 months after the differentiation in iPSC-organoid models from familial AD and Down syndrome patients^35,37, those from sporadic AD patients likely have Aβ deposition 6 months after differentiation²⁸. Thus, longer differentiation duration may be required to assess amyloid deposition in iPSC-derived organoids from sporadic AD cases.

Another pathological hallmark of AD is the abnormal phosphorylation, mislocalization, and aggregation of tau^38,39, although tauopathy is also detected in non-AD neurodegenerative diseases including frontotemporal dementia (FTD) and progressive supranuclear palsy (PSP)⁴⁰. Indeed, tau phosphorylation is accelerated in iPSC-derived neurons from AD patients^13,41. Furthermore, APOE4 has been shown to exacerbate AD-related tau pathologies in cell and mouse models as well as human patients^10,42. Our results also found that AD status and APOE4 are independently associated with the increased levels of p-tau in iPSC-derived cerebral organoids. Given that neither APOE4 nor apoE amounts influenced Aβ levels in the iPSC-derived cerebral organoids, APOE4 may accelerate the p-tau accumulation in an Aβ-independent manner. Of note, apoE and p-tau levels show positive correlations in both RIPA and FA fractions, which is consistent with the results from mouse studies showing that Apoe deficiency ameliorates tauopathy⁴². In addition, Wang et al. has shown that p-tau levels were increased in APOE-deficient iPSC-derived neurons when treated with lysates of iPSC-derived neurons with APOE4¹⁹. Thus, although further studies are necessary, APOE4 likely induces a gain-of-toxic effect on p-tau accumulation in iPSC-derived cerebral organoids.

Intriguingly, we found marked increase in cellular apoptosis and reductions of synaptic proteins in mature organoids from AD patients compared to healthy subject-derived cerebral organoids. In supporting our findings, iPSC-derived neurons from a sporadic AD patient have been reported for their hyper-excitable calcium signaling phenotype, increased cytotoxicity and apoptosis, and reduced neurite length in comparison to those from healthy individuals³⁶. Wang et al. also found expression of apoE4 was associated with increased GABAergic neuron degeneration¹⁹. More supportively, the gain-of-toxic-function of apoE4 could be ameliorated when apoE4 was converted to apoE3 by gene editing, which is highly consistent with what we found in the cerebral organoid model¹⁹. In contrast, at the early neuronal differentiation stage, increased mature neuronal markers and synaptic proteins were detected in AD-derived organoids. Given that similar observations were reported in iPSC-derived neurons from AD subjects⁴³, specific factors facilitating excess neuronal maturation in the early stage may be commonly preserved in iPSCs from AD patients, resulting in enhanced neurodegeneration in the late stage. While APOE4 likely leads to early maturation of neurons²⁸ partially consistent with our results, we demonstrated that APOE4 further exacerbates the effects of AD status on apoptosis induction and synaptic loss in cerebral organoids at week 12. Thus, it is tempting to speculate that accelerated neuronal differentiation/maturation is synergistically induced by AD-related factors and APOE4 in the pre-symptomatic stage, and that the mechanistic exhaustion followed by neuronal dysfunction may contribute to disease development, although further studies are needed.

Our transcriptomic analysis demonstrated that an enrichment of gene sets involving RNA metabolisms is predominantly affected by AD status in the iPSC-derived cerebral organoids; while APOE4 did not have a strong effect on the profiles. RNA metabolism process is highly dynamic and requires a complex interplay among RNA-binding proteins (RBPs)⁴⁴. Thus, the disruption of RBP homeostasis is frequently involved in the pathology and genetics of neurodegenerative diseases including FTD and amyotrophic lateral sclerosis (ALS), where aggregation-prone RBPs often co-localize with stress granules^44,45,46. Interestingly, we also found a significant increase in stress granules in organoids from AD patients. While various RBPs have been shown to deposit in tau inclusions^47,48,49, internalized tau sensitizes cells to stress by promoting formation and stability of stress granules⁴⁹. Thus, enhanced p-tau accumulation may trigger transcriptional changes and stress granule formation in organoids from AD patients. There were also reductions of mRNAs in genes related to DNA repair (ERCC4 and DGCR8)^50,51 and anti-apoptosis (PNO1 and HSPA4)^52,53 in the organoids from AD patients. However, further studies are necessary to determine if those are causatively or consequently involved in neurodegeneration.

Although our results demonstrate that iPSC-derived cerebral organoid is a useful modeling system to investigate AD- and APOE4-related phenotypes, we should state several limitations. First, the core of cerebral organoids shows necrosis-like changes after week 12, likely due to the lack of vascular systems maintaining proper supplies of nutrition and gas exchange as well as immune cells to eliminate cell debris and toxic molecules. To address this caveat, emerging technology has further developed cerebral organoid models in which vascular cells^19,54 or microglia are incorporated^28,55,56. Thus, in future studies, we plan to utilize the iPSC-derived cerebral organoid system with neuroimmune and/or neurovascular system to establish a more advanced platform for AD research. Second, there was heterogeneity in organoid size and growth rate even though identical iPSC lines were used. We also noticed that maturation of astrocytes and the detection of Aβ or apoE in cerebral organoids may vary depending on differentiation protocols^28,33,57. For example, Lin et al.²⁸ found lower soluble apoE levels in cerebral organoids with APOE4 compared to those with APOE3, which is different from our results. The differences in the medium used for differentiation and the time points of medium change may potentially contribute to the different results under the stress conditions induced by APOE4 or AD status. The technical differences and the intrinsic variability using iPSC-derived cerebral organoid models should be considered whenever comparing results from different studies. Third, because of the limited availability of patient-derived iPSC lines, we could not match sex among our study groups in our study. Nonetheless, we did adjust sex as a variable in our analysis. Further optimization and standardization in the iPSC-derived cerebral organoid modeling system should be able to overcome these weaknesses.

In conclusion, our study established a true 3-D human cerebral organoid system to address AD pathogenesis. We successfully recapitulated AD-related pathologies related to ATN classification in our model system. APOE4 predominantly aggravates p-tau accumulation, while AD status is associated with higher levels of Aβ and p-tau, apoptosis, synaptic loss and increased stress granule formation. Importantly, APOE4 synergistically accelerates apoptosis and stress granule formation with AD status in our 3D model. Thus, exploring molecular mechanisms in the interaction between APOE4- and AD-related pathways should provide clues as to how APOE4 vastly increases AD risk. Because isogenic conversion of APOE4 into APOE3 reverses much of the AD-related phenotypes in cerebral organoids from AD, APOE4 might be a promising therapeutic target for AD.

Methods

Generation of iPSCs from human skin fibroblasts

Human skin biopsies from normal individuals and AD patients with APOE ε3/ε3 or ε4/ε4 genotype were obtained from Mayo Clinic patients under IRB protocols with patient consent for research, which was approved by the Mayo Clinic Institutional Review Board. APOE genotype was confirmed by Sanger sequencing using DNA samples from fibroblast lines. Cells were cultured in fibroblast medium containing 10% fetal bovine serum (FBS) (Gemini Bio-Products). The iPSCs were generated by electroporation of three episomal vectors into the fibroblasts using the NHDF nucleofector kit (Lonza)^23,58. Three micrograms of expression plasmid mixtures were electroporated into 6 × 10⁵ fibroblasts with 100 μl transfection reagents. After transfection, fibroblasts were plated onto a 100 mm Matrigel (Corning) coated dish. After 5 days of culture, the fibroblast medium was replaced with TeSR-E7 complete medium (Stemcell Technologies) and changed every day. IPSC colonies were isolated and expanded after 3–4 weeks in culture. The iPSC colonies were passaged using Dispase (Stemcell Technologies) and subjected to rock inhibitor Y27632 (Sigma-Aldrich) treatment for the first 24 h.

Trilineage differentiation of human iPSCs

Three germs layer differentiation was used to confirm the pluripotency of iPSCs using the STEMdiff Trilineage Differentiation kit (Stemcell Technologies) according to manufacturer’s instructions with some modification⁵⁹. When cells were ~70% confluent, iPSCs were passaged using Accutase (Stemcell Technologies) and plated into an AggreWell™800 plate (Stemcell Technologies) to form embryonic bodies (EB) in mTeSR1 medium for 2 days. EBs were washed out from the AggreWell plate and transferred to 6-well non-tissue culture treated plates (Corning) in specific differentiation mediums for each lineage. EBs were subjected to differentiation into mesoderm and endoderm lineages for 5 days or an ectoderm lineage for 7 days, then seeded onto Matrigel-coated plates for further analysis. Differentiation was assessed by immunostaining for germ layer-specific markers (Endoderm: SOX17; Mesoderm: Brachyury; Ectoderm: Nestin/Sox2).

Cerebral organoid culture

Cerebral organoids were generated according to manufacturer’s instructions of the commercial STEMdiff™ Cerebral Organoid Kit (Stemcell Technologies). On day 0, any pre-differentiated cells in the hiPSCs culture were removed by scraping under a microscope. Human iPSC colonies were dissociated into single-cell suspension with Accutase. In total, 9000 cells were then plated into each well of a U-bottom ultra-low-attachment 96-well plate in EB formation media (medium A) supplemented with 10 μM Y27632. An extra 100 μl EB medium was added on day 2 and day 4, respectively. EBs were moved to 24-well low attachment plates in neural induction medium (medium B) for another 3–5 days. EBs were further embedded into 20 μl of matrigel and cultured in neural expansion medium (medium C + D) for 3 days in 6-well low attachment plates for organoid formation. In the final stage, the organoids were transferred to 10-cm dishes with neural culture medium (medium E), and moved to an orbital shaker for further culture. Medium E was replaced with neuronal maturation medium after 4 weeks, which was composed of DMEM/F12 + Neurobasal Medium (1:1) supplemented with N2, B27, BDNF (20 ng/ml), GDNF (20 ng/ml), ascorbic acid (200 μM), and dbcAMP (100 nM) (Sigma-Aldrich)⁶⁰. Cerebral organoids were collected at different time points for immunostaining and other biochemical analysis. We differentiated 20 lines of iPSCs into cerebral organoids within a close time frame, and used cerebral organoids with normal size range (typically 2–4 mm of diameter) for our experiments.

Immunostaining

At pre-defined time points, cerebral organoids were fixed in 4% paraformaldehyde for 30 min then washed with PBS three times. After fixation, organoids were dehydrated with 30% sucrose in PBS at 4 °C. Organoids were then embedded with optical cutting temperature (OCT) compound (VWR) and frozen on dry ice. Frozen tissue was sectioned at 30 μm using a cryostat and collected on ultra-frosted glass microscope slides. Sections were stored at −20 °C. For immunostaining, sections were permeabilized in 0.25% Triton X-100 and blocked with blocking buffer containing 4% normal donkey serum, 2% BSA, and 1 M glycine in PBS. Sections were then incubated with primary antibodies in blocking buffer overnight at 4 °C. The information of primary antibodies and their dilutions used in this study are as follows: Nanog (Cell Signaling, 4903, 1:300), TRA-1-60 (Abcam, ab16288, 1:300), SSEA4 (Abcam, ab16287, 1:300), Sox17 (Abcam, ab84990, 1:300), Brachyury (R&D, AF2085, 1:300), Nestin (Abcam, ab18102, 1:500), Sox2 (Abcam, ab97959, 1:500), TUJ1 (Abcam, ab78078, 1:1000), TUJ1 (Sigma, T2200, 1:1000), CTIP2 (Abcam, ab18465, 1:100), SATB2 (Abcam, ab34735, 1:100), GFAP (Millipore, MAB360, 1:300), cleaved Caspase-3 (Cell Signaling Technology, 9579, 1:300), AT8 (Thermo Fisher Scientific, MN1020, 1:300), and G3BP (BD, 611126, 1:300). After washing three times with PBS, samples were incubated with fluorescently conjugated secondary antibodies (Alexa Fluor 488 and 594 conjugates, Invitrogen, 1:500) for 2 h at room temperature and washed three times with PBS before mounting with the glass coverslip. Fluorescent signals were detected by fluorescence microscopy (model IX71 Invert, Olympus), confocal laser scanning fluorescent microscopy (model LSM510 Invert, Carl Zeiss), and Keyence fluorescence microscopy (model BZ-X, Keyence). The fluorescence signal was quantified using ImageJ software.

Tissue processing

Cerebral organoids were harvested at pre-defined time points and lysed with RIPA Lysis and Extraction Buffer supplemented with Protease and Phosphatase Inhibitor Cocktails for Cell Lysis (Roche). Lysed samples were sonicated and incubated for 60 min on ice. Samples were centrifuged in an ultracentrifuge (Beckman–Coulter) at 100,000 × g, for 1 h at 4 °C. Supernatants (soluble fraction) were collected and the pellet was resuspended in 50 μl of 70% formic acid (FA), sonicated, and neutralized with 2.5 M Tris buffer (pH 8.5) (insoluble fraction). Total protein concentration in the soluble fraction was determined using a Pierce BCA Protein Assay Kit.

Western blotting

Samples in soluble fraction were loaded into a 4–20% sodium dodecyl sulfate-polyacrylamide gel (Bio-Rad), and transferred to PVDF Immobilon FL membranes (Millipore). After blocking with 5% non-fat milk in PBS, membranes were blotted overnight with primary antibodies in 5% non-fat milk containing 0.01% Tween-20, and then probed with LI-COR IRDye secondary antibodies or horseradish peroxidase-conjugated secondary antibody, detected by SuperSignal West Femto Chemiluminescent Substrate (Pierce). The information of primary antibodies and their dilutions used in this study are as follows: Cleaved Caspase-3 (Cell Signaling Technology, 9661, 1:1000), Caspase-3 (Cell Signaling Technology, 9662, 1:1000), PSD95 (Abcam, ab2723, 1:1000), Synaptophysin (Abcam, ab8049, 1:1000), AT8 (Thermo Fisher Scientific, MN1020, 1:1000), PHF1 (Abcam, ab184951, 1:1000), Tau5 (Millipore, 577801, 1:1000), APP (Thermo Fisher Scientific, 14-9749-82, 1:1000), GFAP (Millipore, MAB360, 1:1000), EEA1 (Cell Signaling Technology, 2411, 1:1000), LAMP1 (Cell Signaling Technology, 9091, 1:1000), G3BP (BD, 611126, 1:1000), ERCC4 (Fitzgerald, 10R-4026, 1:1000), POLR3A (Abcam, ab96328, 1:1000), and HSPA4 (Cell Signaling Technology, 3303, 1:1000). All uncropped blot images with molecular size markers were available in the Source data file.

ELISA quantification

The Aβ40, Aβ42, sAPPα, sAPPβ, CTF-β, and phospho-tau levels in both RIPA and formic acid fractions were measured using the Human β-Amyloid (1–40) ELISA Kit (Thermo Fisher, KHB3481), Human β-Amyloid (1–42) ELISA Kit (Thermo Fisher, KHB3441), Human sAPPα Assay Kit (IBL, 27734), Human sAPPβ Assay Kit (IBL, 27732), Human CTF-β ELISA kit (IBL, 27776), and the Tau (Phospho) [pS396] Human ELISA Kit (Thermo Fisher, KHB7031) according to the manufacturers’ instructions. For human apoE ELISA, 96-well plates were coated overnight with an apoE antibody (WUE4) in carbonate buffer at 4 °C. The plates were blocked with 1% milk in PBS, and washed three times with PBS. Recombinant apoE3 and apoE4 (Fitzgerald) were used as standards for the ELISA. Samples were diluted and incubated at 4 °C overnight. The plates were washed and incubated with biotin-conjugated goat anti-apoE antibody (Meridian Life Science) for 2 h at room temperature. After incubation with Horseradish Peroxidase Avidin D (Vector Laboratories) for 90 min at room temperature, the plate was developed by adding tetramethylbenzidine Super Slow substrate (Sigma). The reaction was stopped and read at 450 nm with a microplate reader (Biotek). Results were normalized to total protein concentration of the cell lysate.

RT-qPCR

RNA was extracted via the Trizol/chloroform method, followed by DNase and cleanup using the RNase-Free DNase Set and the RNeasy Mini Kit (QIAGEN). The quantity and quality of all RNA samples for RNA-seq was determined by the Agilent 2100 Bioanalyzer using the Agilent RNA 6000 Nano Chip (Agilent Technologies, CA). The cDNA was prepared with the iScript cDNA synthesis kit (Bio-Rad). Real-time qPCR was conducted with Universal SYBR Green Supermix (Bio-Rad) using an iCycler thermocycler (Bio-Rad). The 2exp (−ΔΔCt) method was used to determine the relative expression of each gene with ACTB gene coding β-actin as a reference. The primers used to amplify target genes by RT-qPCR are as follows: ACTB F (5′-CTGGCACCACACCTTCTACAATG-3′) and R (5′-AATGTCACGCACGATTTCCCGC-3′), MAP2 F (5′-CAGGTGGCGGACGTGTGAAAATTGAGAGTG-3′) and R (5′-CACGCTGGA TCTGCCTGGGGACTGTG-3′), CTIP2 F (5′-GTTGTGCAAATGTAGCTGGAA-3′) and R (5′-GAAGATGACCACCTGCTCTC-3′), SATB2 F (5′-CCTTACGCAGAATCTCAGACAA-3′) and R (5′-CCAGATATCTACCAGCAAGTCAG-3′), GFAP F (5′-CTGTTGCCAGAGATG GAGGTT-3′) and R (5′-TCATCGCTCAGGAGGTCCTT-3′), ERCC4 F (5′-TTTGTGAGG AAACTGTATCTG TGG-3′) and R (5′-GTCTGTATAGCAAGCATGGTAGG-3′), DGCR8 (PPH17455A, QIAGEN), POLR3A F (5′-TCTGGAGACCTGTAGGGACAA-3′) and R (5′-CTGGCTCACCAATGCTCT-3′), CLP1 (PPH10372A, QIAGEN), HSPA4 F (5′-GCAACA GCAGCAGACACCAGC-3′) and R (5′-GCCTTCTTGGCTTGGGGTGGT-3′), PNO1 (PPH15310A, QIAGEN), DND1 F (5′-CTCCACAGGCACCCTGAATG-3′) and R (5′-GGTGCCATAGGTCCCTGTCC-3′), NR2F6 F (5′-TTTAGTCGATTTCACGCGGAC-3′) and R (5′-ATCTGTTCAGACGGGTACTCG-3′), VPS9D1 F (5′-CTCAGACGTCCCCAGGAA CTG-3′) and R (5′-GCCAGACAAGGACAGCTCGTTC-3′), CTH F (5′-GAAGACCTACT GGAAGAT-3′) and R (5′-GGAATACTAGCTGTGACT-3′), CDHR3 F (5′-TGTGAAGGA TGAGGTTGGTG-3′) and R (5′-TCCAGGGTTTGCTCTTTCTAC-3′), CHIC1 (PPH09017A, QIAGEN), ASB3 F (5′-CATACTTATTTCATCGGGTGC-3′) and R (5′-GGTAACTGCCA ACTGTCCTC-3′), ARRDC3 F (5′-ATCCCAGTGTGATGTGACGA-3′) and R (5′-TTTGCAA CAGAATCGGAAAA-3′), TUBB3 F (5′-CCTCCGTGTAGTGACCCTT-3′) and R (5′-GGCCTTTGGACATCTCTTCAG-3′).

RNA-seq, quality control, and normalization

Twenty mRNA samples were sequenced at Mayo Clinic using Illumina HiSeq 4000. Reads were mapped to the human reference genome hg38. Raw gene read counts and sequencing quality control were generated using the Mayo Clinic RNA-Seq analytic pipeline: MAP-RSeq Version 3.0⁶¹. Conditional quantile normalization (CQN) was performed on the raw gene counts to correct for gene length differences, GC bias, global technical variations, and to obtain similar quantile-by-quantile distributions of gene expression levels across samples⁶². Based on the bi-modal distribution of the CQN-normalized and log2-transformed reads per kb per million (RPKM) gene expression values, genes with an average of log2 RPKM > = 3 in at least one group were considered expressed. Using this selection threshold, 18,291 genes were included in the downstream analysis.

Differential gene expression and pathway analysis

Differential gene expression analyses were performed using the Partek Genomics Suite (Partek Inc., St. Louis, MO). Gene expression between AD patient- and healthy subject-derived cerebral organoids, between APOE genotypes, and the interaction between AD disease status and APOE genotype were calculated using analyses of variance models (ANOVA). Pathway analyses of differentially expressed genes, defined by p < 0.05 and |fold change | ≥1.5, were performed using MetaCore pathway analysis (© MetaCore (Feb 2020) of Clarivate Analytics. All rights reserved).

DEGs cell-type distribution analysis of RNA-seq data

Cell proportion was estimated using marker genes described in BRETIGEA⁶³, a published reference dataset, and the CIBERSORT program⁶⁴. First, the top 50 marker genes of neuron, astrocyte, and oligodendrocyte cells were obtained from BRETIGEA (version 1.0.0)⁶³, among which 30 neuronal markers, 27 astrocytic markers, and 19 oligodendrocyte markers are also present in both reference dataset and our iPSC dataset. Next, expression levels of these markers in sorted neurons, astrocytes, and oligodendrocytes were obtained from the reference dataset of Zhang et al.⁶⁵. Finally, CIBERSORT analytical tool (version 1.0.1)⁶⁴ was applied to estimate cellular composition in the iPSC dataset using the selected maker gene expression in reference dataset and in the iPSC dataset for neurons, astrocytes, and oligodendrocytes. To estimate in which cell-type genes might be regulated, we applied CellCODE R package⁶⁶ to assess the interaction between group variable and surrogate proportion variables of each cell type. Because oligodendrocytes compose negligible proportion of the cell population in all but one sample, oligodendrocyte was not considered in the following analysis. Specifically, getAllSPV function in CellCODE was applied to refine the top 50 markers of neuron and astrocyte obtained from BRETIGEA, and to obtain surrogate variables of neuron and astrocyte populations using the remaining markers through singular value decomposition. Next, cellPopT function was applied to calculate the t statistics of the interaction term between group and the surrogate variables. Cell type with the largest t statistics is the designated cell type.

WGCNA

To identify groups of genes that are correlated with AD disease status or with APOE genotype, we performed weighted gene co-expression network analysis (WGCNA)⁶⁷ using the log2-transformed, CQN-normalized gene expression values. We used the soft power of 16, hybrid dynamic tree cutting, a minimum module size of 50 genes, and a minimum height for merging modules at 0.3 to build a signed hybrid co-expression networks. Each gene module was summarized by the first principal component of the scaled module expression profiles (module eigengene). Each module was assigned a unique color identifier, and genes that did not fulfill these criteria for any of the modules were assigned to the gray module. To assess the correlation of modules to AD disease status and APOE genotype, we defined healthy subject-derived cerebral organoids as 0 and AD patient-derived cerebral organoids as 1; and defined the APOE3 genotype as 0, and the APOE4 genotype as 1. Modules were annotated using R package anRichment (https://horvath.genetics.ucla.edu/html/CoexpressionNetwork/GeneAnnotation/). The connection among the top hub genes in the yellow module was visualized using VisANT⁶⁸.

Statistical analyses and reproducibility

ANCOVA for APOE4, AD status, sex, sampling age and source of iPSCs was performed to determine the interaction effect between APOE4 and AD status on each continuous variable for western blotting, ELISA, or RT-qPCR from four groups of cerebral organoids, followed by Tukey–Kramer tests to compare between selected two groups using JMP software version 15.0. For batch difference and isogenic cerebral organoids comparison, the Mann–Whitney U test was performed to determine the significance using GraphPad Prism version 8.0. Spearman correlation analysis was used to analyze the association between proteins using GraphPad Prism. Experiments were repeated in two independent differentiation batches. Data were presented as mean ± SEM. A p value of <0.05 was considered statistically significant.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability

Full scans of the gels and blots are available in Source data file. All relevant data are available from the corresponding author upon reasonable request. All the other data supporting the findings of this study are available within the article and its supplementary information files and from the corresponding author upon reasonable request. A reporting summary for this article is available as a Supplementary Information file. The RNA-seq data are available via the AD Knowledge Portal [https://adknowledgeportal.synapse.org]. The AD Knowledge Portal is a platform for accessing data, analyses, and tools generated by the Accelerating Medicines Partnership (AMP-AD) Target Discovery Program and other National Institute on Aging (NIA)-supported programs to enable open-science practices and accelerate translational learning. The data, analyses, and tools are shared early in the research cycle without a publication embargo on secondary use. Data are available for general research use according to the following requirements for data access and data attribution [https://adknowledgeportal.synapse.org/DataAccess/Instructions]. For access to content described in this manuscript see [https://doi.org/10.7303/syn22307008]. Source data are provided with this paper.

Change history

05 May 2021
A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-23081-4

References

Selkoe, D. J. Toward a remembrance of things past: deciphering Alzheimer disease. Harvey Lect. 99, 23–45 (2003).
PubMed Google Scholar
Jack, C. R. Jr. et al. Hypothetical model of dynamic biomarkers of the Alzheimer’s pathological cascade. Lancet Neurol. 9, 119–128 (2010).
Article CAS PubMed PubMed Central Google Scholar
DeTure, M. A. & Dickson, D. W. The neuropathological diagnosis of Alzheimer’s disease. Mol. Neurodegener. 14, 32 (2019).
Article PubMed PubMed Central Google Scholar
Corder, E. H. et al. Gene dose of apolipoprotein E type 4 allele and the risk of Alzheimer’s disease in late onset families. Science 261, 921–923 (1993).
Article CAS PubMed ADS Google Scholar
Kim, J., Basak, J. M. & Holtzman, D. M. The role of apolipoprotein E in Alzheimer’s disease. Neuron 63, 287–303 (2009).
Article CAS PubMed PubMed Central Google Scholar
Liu, C. C., Liu, C. C., Kanekiyo, T., Xu, H. & Bu, G. Apolipoprotein E and Alzheimer disease: risk, mechanisms and therapy. Nat. Rev. Neurol. 9, 106–118 (2013).
Article CAS PubMed PubMed Central Google Scholar
Yu, J. T., Tan, L. & Hardy, J. Apolipoprotein E in Alzheimer’s disease: an update. Annu Rev. Neurosci. 37, 79–100 (2014).
Article CAS PubMed Google Scholar
Najm, R., Jones, E. A. & Huang, Y. Apolipoprotein E4, inhibitory network dysfunction, and Alzheimer’s disease. Mol. Neurodegener. 14, 24 (2019).
Article PubMed PubMed Central Google Scholar
Huynh, T. V., Davis, A. A., Ulrich, J. D., Holtzman, D. M. & Apolipoprotein, E. and Alzheimer’s disease: the influence of apolipoprotein E on amyloid-beta and other amyloidogenic proteins. J. Lipid Res. 58, 824–836 (2017).
Article CAS PubMed PubMed Central Google Scholar
Huang, Y. & Mucke, L. Alzheimer mechanisms and therapeutic strategies. Cell 148, 1204–1222 (2012).
Article CAS PubMed PubMed Central Google Scholar
Golde, T. E., Schneider, L. S. & Koo, E. H. Anti-abeta therapeutics in Alzheimer’s disease: the need for a paradigm shift. Neuron 69, 203–213 (2011).
Article CAS PubMed PubMed Central Google Scholar
Ashe, K. H. & Zahs, K. R. Probing the biology of Alzheimer’s disease in mice. Neuron 66, 631–645 (2010).
Article CAS PubMed PubMed Central Google Scholar
Israel, M. A. et al. Probing sporadic and familial Alzheimer’s disease using induced pluripotent stem cells. Nature 482, 216–220 (2012).
Article CAS PubMed PubMed Central ADS Google Scholar
Yagi, T. et al. Modeling familial Alzheimer’s disease with induced pluripotent stem cells. Hum. Mol. Genet. 20, 4530–4539 (2011).
Article CAS PubMed Google Scholar
Kondo, T. et al. Modeling Alzheimer’s disease with iPSCs reveals stress phenotypes associated with intracellular Abeta and differential drug responsiveness. Cell Stem Cell 12, 487–496 (2013).
Article CAS PubMed Google Scholar
Liu, Q. et al. Effect of potent gamma-secretase modulator in human neurons derived from multiple presenilin 1-induced pluripotent stem cell mutant carriers. JAMA Neurol. 71, 1481–1489 (2014).
Article PubMed PubMed Central Google Scholar
Centeno, E. G. Z., Cimarosti, H. & Bithell, A. 2D versus 3D human induced pluripotent stem cell-derived cultures for neurodegenerative disease modelling. Mol. Neurodegener. 13, 27 (2018).
Article PubMed PubMed Central CAS Google Scholar
Muratore, C. R. et al. The familial Alzheimer’s disease APPV717I mutation alters APP processing and Tau expression in iPSC-derived neurons. Hum. Mol. Genet. 23, 3523–3536 (2014).
Article CAS PubMed PubMed Central Google Scholar
Wang, C. et al. Gain of toxic apolipoprotein E4 effects in human iPSC-derived neurons is ameliorated by a small-molecule structure corrector. Nat. Med. 24, 647–657 (2018).
Article CAS PubMed PubMed Central Google Scholar
Lancaster, M. A. et al. Cerebral organoids model human brain development and microcephaly. Nature 501, 373–379 (2013).
Article CAS PubMed ADS Google Scholar
Lancaster, M. A. & Knoblich, J. A. Generation of cerebral organoids from human pluripotent stem cells. Nat. Protoc. 9, 2329–2340 (2014).
Article CAS PubMed PubMed Central Google Scholar
Renner, M. et al. Self-organized developmental patterning and differentiation in cerebral organoids. EMBO J. 36, 1316–1329 (2017).
Article CAS PubMed PubMed Central Google Scholar
Okita, K. et al. A more efficient method to generate integration-free human iPS cells. Nat. Methods 8, 409–412 (2011).
Article CAS PubMed Google Scholar
Zhao, J. et al. APOE epsilon4/epsilon4 diminishes neurotrophic function of human iPSC-derived astrocytes. Hum. Mol. Genet. 26, 2690–2700 (2017).
Article CAS PubMed PubMed Central Google Scholar
Brookhouser, N., Zhang, P., Caselli, R., Kim, J. J. & Brafman, D. A. Generation and characterization of two human induced pluripotent stem cell (hiPSC) lines homozygous for the Apolipoprotein e4 (APOE4) risk variant-Alzheimer’s disease (ASUi005-A) and healthy non-demented control (ASUi006-A). Stem Cell Res. 32, 145–149 (2018).
Article CAS PubMed PubMed Central Google Scholar
Brookhouser, N., Zhang, P., Caselli, R., Kim, J. J. & Brafman, D. A. Generation and characterization of human induced pluripotent stem cell (hiPSC) lines from an Alzheimer’s disease (ASUi001-A) and non-demented control (ASUi002-A) patient homozygous for the Apolipoprotein e4 (APOE4) risk variant. Stem Cell Res. 24, 160–163 (2017).
Article CAS PubMed Google Scholar
Yamazaki, Y. et al. Apolipoprotein E and Alzheimer disease: pathobiology and targeting strategies. Nat. Rev. Neurol. 15, 501–518 (2019).
Lin, Y. T. et al. APOE4 causes widespread molecular and cellular alterations associated with Alzheimer’s disease phenotypes in human iPSC-derived brain cell types. Neuron 98, 1141–1154.e1147 (2018).
Article CAS PubMed PubMed Central Google Scholar
Niewidok, B. et al. Single-molecule imaging reveals dynamic biphasic partition of RNA-binding proteins in stress granules. J. Cell Biol. 217, 1303–1318 (2018).
Article CAS PubMed PubMed Central Google Scholar
Van Treeck, B. et al. RNA self-assembly contributes to stress granule formation and defining the stress granule transcriptome. Proc. Natl Acad. Sci. USA 115, 2734–2739 (2018).
Article PubMed PubMed Central CAS Google Scholar
Gaugler, J. et al. 2019 Alzheimer’s disease facts and figures. Alzheimer’s & Dementia 15, 321–387 (2019).
Jack, C. R. Jr. et al. NIA-AA Research Framework: toward a biological definition of Alzheimer’s disease. Alzheimers Dement 14, 535–562 (2018).
Article PubMed PubMed Central Google Scholar
Quadrato, G. et al. Cell diversity and network dynamics in photosensitive human brain organoids. Nature 545, 48–53 (2017).
Article CAS PubMed PubMed Central ADS Google Scholar
Choi, S. H. et al. A three-dimensional human neural cell culture model of Alzheimer’s disease. Nature 515, 274–278 (2014).
Article CAS PubMed PubMed Central ADS Google Scholar
Gonzalez, C. et al. Modeling amyloid beta and tau pathology in human cerebral organoids. Mol. Psychiatry 23, 2363–2374 (2018).
Article CAS PubMed PubMed Central Google Scholar
Balez, R. et al. Neuroprotective effects of apigenin against inflammation, neuronal excitability and apoptosis in an induced pluripotent stem cell model of Alzheimer’s disease. Sci. Rep. 6, 31450 (2016).
Article CAS PubMed PubMed Central ADS Google Scholar
Raja, W. K. et al. Self-organizing 3D human neural tissue derived from induced pluripotent stem cells recapitulate Alzheimer’s disease phenotypes. PLoS ONE 11, e0161969 (2016).
Article PubMed PubMed Central CAS Google Scholar
Alonso, A. C., Zaidi, T., Grundke-Iqbal, I. & Iqbal, K. Role of abnormally phosphorylated tau in the breakdown of microtubules in Alzheimer disease. Proc. Natl Acad. Sci. USA 91, 5562–5566 (1994).
Article CAS PubMed PubMed Central ADS Google Scholar
Mi, K. & Johnson, G. V. The role of tau phosphorylation in the pathogenesis of Alzheimer’s disease. Curr. Alzheimer Res 3, 449–463 (2006).
Article CAS PubMed Google Scholar
Meeter, L. H., Kaat, L. D., Rohrer, J. D. & van Swieten, J. C. Imaging and fluid biomarkers in frontotemporal dementia. Nat. Rev. Neurol. 13, 406–419 (2017).
Article CAS PubMed Google Scholar
Ochalek, A. et al. Neurons derived from sporadic Alzheimer’s disease iPSCs reveal elevated TAU hyperphosphorylation, increased amyloid levels, and GSK3B activation. Alzheimers Res. Ther. 9, 90 (2017).
Article PubMed PubMed Central CAS Google Scholar
Shi, Y. et al. ApoE4 markedly exacerbates tau-mediated neurodegeneration in a mouse model of tauopathy. Nature 549, 523–527 (2017).
Article PubMed PubMed Central ADS CAS Google Scholar
Meyer, K. et al. REST and neural gene network dysregulation in iPSC models of Alzheimer’s disease. Cell Rep. 26, 1112–1127.e1119 (2019).
Article CAS PubMed PubMed Central Google Scholar
Nussbacher, J. K., Tabet, R., Yeo, G. W. & Lagier-Tourenne, C. Disruption of RNA metabolism in neurological diseases and emerging therapeutic interventions. Neuron 102, 294–320 (2019).
Article CAS PubMed PubMed Central Google Scholar
Liu, E. Y., Cali, C. P. & Lee, E. B. RNA metabolism in neurodegenerative disease. Dis. Model Mech. 10, 509–518 (2017).
Article CAS PubMed PubMed Central Google Scholar
Johnson, E. C. B. et al. Deep proteomic network analysis of Alzheimer’s disease brain reveals alterations in RNA binding proteins and RNA splicing associated with disease. Mol. Neurodegener. 13, 52 (2018).
Article CAS PubMed PubMed Central Google Scholar
Vanderweyde, T. et al. Contrasting pathology of the stress granule proteins TIA-1 and G3BP in tauopathies. J. Neurosci. 32, 8270–8283 (2012).
Article CAS PubMed PubMed Central Google Scholar
Maziuk, B. F. et al. RNA binding proteins co-localize with small tau inclusions in tauopathy. Acta Neuropathol. Commun. 6, 71 (2018).
Article PubMed PubMed Central CAS Google Scholar
Brunello, C. A., Yan, X. & Huttunen, H. J. Internalized Tau sensitizes cells to stress by promoting formation and stability of stress granules. Sci. Rep. 6, 30498 (2016).
Article CAS PubMed PubMed Central ADS Google Scholar
Brookman, K. W. et al. ERCC4 (XPF) encodes a human nucleotide excision repair protein with eukaryotic recombination homologs. Mol. Cell Biol. 16, 6553–6562 (1996).
Article CAS PubMed PubMed Central Google Scholar
Calses, P. C. et al. DGCR8 mediates repair of UV-induced DNA damage independently of RNA processing. Cell Rep. 19, 162–174 (2017).
Article CAS PubMed PubMed Central ADS Google Scholar
Shen, A. et al. EBF1-mediated upregulation of ribosome assembly factor PNO1 contributes to cancer progression by negatively regulating the p53 signaling pathway. Cancer Res. 79, 2257–2270 (2019).
Article CAS PubMed Google Scholar
Zuo, D., Subjeck, J. & Wang, X. Y. Unfolding the role of large heat shock proteins: new insights and therapeutic implications. Front. Immunol. 7, 75 (2016).
Article PubMed PubMed Central CAS Google Scholar
Mansour, A. A. et al. An in vivo model of functional and vascularized human brain organoids. Nat. Biotechnol. 36, 432–441 (2018).
Article CAS PubMed PubMed Central Google Scholar
Abreu, C. M. et al. Microglia increase inflammatory responses in iPSC-derived human brainspheres. Front. Microbiol. 9, 2766 (2018).
Article PubMed PubMed Central Google Scholar
Ormel, P. R. et al. Microglia innately develop within cerebral organoids. Nat. Commun. 9, 4167 (2018).
Article PubMed PubMed Central ADS CAS Google Scholar
Linkous, A. et al. Modeling patient-derived glioblastoma with cerebral organoids. Cell Rep. 26, 3203–3211.e3205 (2019).
Article CAS PubMed PubMed Central Google Scholar
Wren, M. C. et al. Frontotemporal dementia-associated N279K tau mutant disrupts subcellular vesicle trafficking and induces cellular stress in iPSC-derived neural stem cells. Mol. Neurodegener. 10, 46 (2015).
Article PubMed PubMed Central CAS Google Scholar
Awe, J. P. et al. Putative immunogenicity expression profiling using human pluripotent stem cells and derivatives. Stem Cells Transl. Med. 4, 136–145 (2015).
Article CAS PubMed PubMed Central Google Scholar
Nimsanor, N. et al. Induced pluripotent stem cells (iPSCs) derived from a symptomatic carrier of a S305I mutation in the microtubule-associated protein tau (MAPT)-gene causing frontotemporal dementia. Stem Cell Res. 17, 564–567 (2016).
Article CAS PubMed Google Scholar
Kalari, K. R. et al. MAP-RSeq: Mayo analysis pipeline for RNA sequencing. BMC Bioinforma. 15, 224 (2014).
Article CAS Google Scholar
Hansen, K. D., Irizarry, R. A. & Wu, Z. Removing technical variability in RNA-seq data using conditional quantile normalization. Biostatistics 13, 204–216 (2012).
Article PubMed PubMed Central MATH Google Scholar
McKenzie, A. T. et al. Brain cell type specific gene expression and co-expression network architectures. Sci. Rep. 8, 8868 (2018).
Article PubMed PubMed Central ADS CAS Google Scholar
Newman, A. M. et al. Robust enumeration of cell subsets from tissue expression profiles. Nat. Methods 12, 453–457 (2015).
Article CAS PubMed PubMed Central Google Scholar
Zhang, Y. et al. Purification and characterization of progenitor and mature human astrocytes reveals transcriptional and functional differences with mouse. Neuron 89, 37–53 (2016).
Article CAS PubMed Google Scholar
Chikina, M., Zaslavsky, E. & Sealfon, S. C. CellCODE: a robust latent variable approach to differential expression analysis for heterogeneous cell populations. Bioinformatics 31, 1584–1591 (2015).
Article CAS PubMed PubMed Central Google Scholar
Langfelder, P. & Horvath, S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinforma. 9, 559 (2008).
Article CAS Google Scholar
Hu, Z. et al. VisANT 4.0: integrative network platform to connect genes, drugs, diseases and therapies. Nucleic Acids Res. 41, W225–W231 (2013).
Article PubMed PubMed Central Google Scholar

Download references

Acknowledgements

We are grateful to Dr. Yadong Huang and Chengzhong Wang for providing the APOE isogenic lines. This work was partially supported by Mayo Clinic Center for Regenerative Medicine, Neuroregeneration Lab. Lab Funding: This work was supported by NIH grants RF1AG051504, R37AG027924, RF1AG046205, RF1AG057181, R01AG066395, P01NS074969, and P30AG062677 (to G.B.), Younkin Fellowship and Alzheimer’s Association Research Fellowship 2018-AARF-592302 (to J.Z.), and R01AG061796 (to N.E.-T.). Dr. Zbigniew K. Wszolek is partially supported by the Mayo Clinic Center for Regenerative Medicine, the gifts from The Sol Goldman Charitable Trust, the Donald G. and Jodi P. Heeringa Family, the Haworth Family Professorship in Neurodegenerative Diseases fund, and by Albertson Parkinson’s Research Foundation.

Author information

Authors and Affiliations

Department of Neuroscience, Mayo Clinic, Jacksonville, FL, 32224, USA
Jing Zhao, Yuan Fu, Yu Yamazaki, Mary D. Davis, Chia-Chen Liu, Wenyan Lu, Kai Chen, Lin Jia, Yuka A. Martens, Lucy Job, Francis Shue, Thanh Thanh Nguyen, Steven G. Younkin, Nilüfer Ertekin-Taner, Takahisa Kanekiyo & Guojun Bu
Center for Regenerative Medicine, Neuroregeneration Lab, Mayo Clinic, Jacksonville, FL, 32224, USA
Jing Zhao, Mary D. Davis, Wenyan Lu, Lucy Job, Takahisa Kanekiyo & Guojun Bu
Department of Health Sciences Research, Mayo Clinic, Jacksonville, FL, 32224, USA
Yingxue Ren, Xue Wang, Yesesri Cherukuri & Yan W. Asmann
Department of Neurology, Mayo Clinic, Jacksonville, FL, 32224, USA
Neill R. Graff-Radford, Zbigniew K. Wszolek & Nilüfer Ertekin-Taner
School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, 85287, USA
David A. Brafman

Authors

Jing Zhao
View author publications
You can also search for this author in PubMed Google Scholar
Yuan Fu
View author publications
You can also search for this author in PubMed Google Scholar
Yu Yamazaki
View author publications
You can also search for this author in PubMed Google Scholar
Yingxue Ren
View author publications
You can also search for this author in PubMed Google Scholar
Mary D. Davis
View author publications
You can also search for this author in PubMed Google Scholar
Chia-Chen Liu
View author publications
You can also search for this author in PubMed Google Scholar
Wenyan Lu
View author publications
You can also search for this author in PubMed Google Scholar
Xue Wang
View author publications
You can also search for this author in PubMed Google Scholar
Kai Chen
View author publications
You can also search for this author in PubMed Google Scholar
Yesesri Cherukuri
View author publications
You can also search for this author in PubMed Google Scholar
Lin Jia
View author publications
You can also search for this author in PubMed Google Scholar
Yuka A. Martens
View author publications
You can also search for this author in PubMed Google Scholar
Lucy Job
View author publications
You can also search for this author in PubMed Google Scholar
Francis Shue
View author publications
You can also search for this author in PubMed Google Scholar
Thanh Thanh Nguyen
View author publications
You can also search for this author in PubMed Google Scholar
Steven G. Younkin
View author publications
You can also search for this author in PubMed Google Scholar
Neill R. Graff-Radford
View author publications
You can also search for this author in PubMed Google Scholar
Zbigniew K. Wszolek
View author publications
You can also search for this author in PubMed Google Scholar
David A. Brafman
View author publications
You can also search for this author in PubMed Google Scholar
Yan W. Asmann
View author publications
You can also search for this author in PubMed Google Scholar
Nilüfer Ertekin-Taner
View author publications
You can also search for this author in PubMed Google Scholar
Takahisa Kanekiyo
View author publications
You can also search for this author in PubMed Google Scholar
Guojun Bu
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

J.Z., C.L., F.S., N.E.T., T.K., and G.B. conceived and designed the project, and wrote the paper. M.D., L.J., S.G.Y., N.R.G.R., Z.W., and D.B. helped with collecting human skin biopsies and generating iPSC lines. J.Z., Y.F., Y.Y., W.L., Y.M., K.C., L.J., and T.N. executed the experiments and analyzed the data. Y.R., X.W., Y.C., and Y.A. performed analysis for RNA-sequencing data.

Corresponding author

Correspondence to Guojun Bu.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature Communications thanks Tracy Young-Pearse and the other anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary Information

Peer Review File

Description of Additional Supplementary Files

Supplementary Data 1

Reporting Summary

Source data

Source Data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Zhao, J., Fu, Y., Yamazaki, Y. et al. APOE4 exacerbates synapse loss and neurodegeneration in Alzheimer’s disease patient iPSC-derived cerebral organoids. Nat Commun 11, 5540 (2020). https://doi.org/10.1038/s41467-020-19264-0

Download citation

Received: 27 December 2019
Accepted: 02 October 2020
Published: 02 November 2020
DOI: https://doi.org/10.1038/s41467-020-19264-0

This article is cited by

Hotspots and trends of microglia in Alzheimer's disease: a bibliometric analysis during 2000–2022
- Lijie Zhang
- Qiuru Yao
- Qing Zeng
European Journal of Medical Research (2024)
Synaptic vesicle glycoprotein 2 A in serum is an ideal biomarker for early diagnosis of Alzheimer’s disease
- Xiaoling Wang
- Xiaomin Zhang
- Peichang Wang
Alzheimer's Research & Therapy (2024)
Neuropathogenesis-on-chips for neurodegenerative diseases
- Sarnai Amartumur
- Huong Nguyen
- Chaejeong Heo
Nature Communications (2024)
hiPSC-based models to decipher the contribution of human astrocytes to Alzheimer’s disease and potential therapeutics
- Julia TCW
- Amaia M. Arranz
Molecular Neurodegeneration (2023)
APOE deficiency impacts neural differentiation and cholesterol biosynthesis in human iPSC-derived cerebral organoids
- Jing Zhao
- Tadafumi C. Ikezu
- Guojun Bu
Stem Cell Research & Therapy (2023)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.