- 1Instituto de Biomedicina and Departamento de Biología Molecular, Universidad de León, León, Spain
- 2Instituto de Desarrollo Ganadero y Sanidad Animal and Departamento de Producción Animal, Universidad de León, León, Spain
- 3Instituto de Neurociencias, Consejo Superior de Investigaciones Científicas and Universidad Miguel Hernández, Alicante, Spain
- 4Instituto de Biomedicina and Departamento de Producción Animal, Universidad de León, León, Spain
Introduction: Neurogenesis is tightly regulated in space and time, ensuring the correct development and organization of the central nervous system. Critical regulators of brain development and morphogenesis in mice include two members of the p53 family: p53 and p73. However, dissecting the in vivo functions of these factors and their various isoforms in brain development is challenging due to their pleiotropic effects. Understanding their role, particularly in neurogenesis and brain morphogenesis, requires innovative experimental approaches.
Methods: To address these challenges, we developed an efficient and highly reproducible protocol to generate mouse brain organoids from pluripotent stem cells. These organoids contain neural progenitors and neurons that self-organize into rosette-like structures resembling the ventricular zone of the embryonic forebrain. Using this model, we generated organoids from p73-deficient mouse cells to investigate the roles of p73 and its isoforms (TA and DNp73) during brain development.
Results and Discussion: Organoids derived from p73-deficient cells exhibited increased neuronal apoptosis and reduced neural progenitor proliferation, linked to compensatory activation of p53. This closely mirrors previous in vivo observations, confirming that p73 plays a pivotal role in brain development. Further dissection of p73 isoforms function revealed a dual role of p73 in regulating brain morphogenesis, whereby TAp73 controls transcriptional programs essential for the establishment of the neurogenic niche structure, while DNp73 is responsible for the precise and timely regulation of neural cell fate. These findings highlight the distinct roles of p73 isoforms in maintaining the balance of neural progenitor cell biology, providing a new understanding of how p73 regulates brain morphogenesis.
1 Introduction
The mammalian brain is an extremely complex structure that relies heavily on the cerebral cortex for higher cognitive function (Borrell, 2019). The development of the cerebral cortex begins with the amplification of neuroepithelial stem cells, which then transition into apical Radial Glial Cells (aRGCs), initiating neurogenesis. aRGCs are the primary type of cortical progenitor cells, constituting the primary germinal layer: the Ventricular Zone (VZ). In rodents, aRGCs primarily undergo symmetric proliferative divisions at the onset of neurogenesis, leading to the expansion of their pool. As development proceeds, aRGC self-amplification decreases with the concomitant increase in asymmetric divisions, giving rise to secondary, basal progenitor cells (indirect neurogenesis) and neurons (direct neurogenesis) (Miyata et al., 2004; Takahashi et al., 1995; Takahashi et al., 1996). Basal progenitors migrate from the apical surface to the basal border of the VZ, forming the Subventricular Zone (SVZ), a secondary germinal zone. In rodents, the SVZ is a significant neurogenic niche in the adult brain, where predominantly excitatory neurons for all cortical layers are generated (Miyata et al., 2004; Noctor et al., 2004; Noctor et al., 2008; Taverna et al., 2014). Towards the end of neurogenesis, most aRGCs differentiate into glial or ependymal cells, depleting the neural stem cell (NSC) pool and marking the end of neurogenesis (Villalba et al., 2021). The neurogenic capacity and homeostasis of the rodents SVZ is dependent on the integrity of the ependymal cells and their capacity to assemble into the unique structures that define this region: the pinwheels (Kuo et al., 2006; Paez-Gonzalez et al., 2011).
The p53 family of transcription factors is involved in several developmental processes (Maeso-Alonso et al., 2021; Van Nostrand et al., 2017), including the regulation of morphogenesis and function of the central nervous system (CNS) (Agostini et al., 2010; Fujitani et al., 2010; Gonzalez-Cano et al., 2010; Hooper et al., 2006; Mendrysa et al., 2011; Moreau et al., 2023; Talos et al., 2010; Xiong et al., 2020). Among the p53 family members, p73 is considered a key regulator of CNS development and function. Indeed, p73 knockout (p73KO) mice exhibit a plethora of severe developmental defects (Marques et al., 2019; Nemajerova and Moll, 2019; Rufini et al., 2011) including anomalies in the organization and homeostasis of the neurogenic niche of the SVZ (Gonzalez-Cano et al., 2016). Global p73KO mice display hippocampal dysgenesis characterized by a missing or truncated dentate gyrus, cortical thinning, loss of sympathetic and Cajal-Retzius neurons and enlarged ventricles with hydrocephalus (Medina-Bolivar et al., 2014; Meyer et al., 2004; Yang et al., 2000). Moreover, p73 deficiency results in deficient ependyma that fails to assemble into neurogenic pinwheels and consequently to organize a functional SVZ niche (Gonzalez-Cano et al., 2016). In mice, p73 is necessary to maintain self-renewal and proliferation, and to inhibit senescence of NSCs or neural progenitor cells (NPCs), therefore maintaining the neurogenic pool (Agostini et al., 2010; Fujitani et al., 2010; Gonzalez-Cano et al., 2010; Talos et al., 2010). However, despite strong in vivo and in vitro evidence supporting the relevance of p73 in neural development, the specific cellular and molecular mechanisms of action remain unknown.
The Trp73 gene gives rise to two different isoforms, TAp73 and DNp73, which are typically expressed in a tissue specific manner (Grespi et al., 2012). TAp73 is transcriptionally competent, while the N-terminal truncated DNp73 acts as dominant-negative of both p53 and TAp73; however, DNp73 has its own unique functions (Marques-Garcia et al., 2009; Niemantsverdriet et al., 2012). Consequently, the activity and function of p73 result from a context-dependent balance between its isoforms. The brain phenotypes observed in p73 isoform-specific knockout mice are milder compared to those in global p73KO mice. DNp73 knockout (DNp73KO) mice exhibit neuronal loss (Tissir et al., 2009; Wilhelm et al., 2010), while TAp73 knockout (TAp73KO) mice display hippocampal dysgenesis (Fujitani et al., 2010; Tomasini et al., 2008). Notably, some phenotypic features present in p73KO mice do not manifest in any of the isoform-specific knockouts. These findings suggest that complex interactions between the p73 isoforms play a crucial role in brain development and morphogenesis, which are yet to be fully understood.
It is well established that inherent differences exist between the in vivo models and the traditional 2D in vitro systems (Gopalakrishnan, 2019), highlighting the need for multi-system strategies to generate comprehensive datasets and dissect the molecular and cellular function of p73 in neural development and brain morphogenesis. Brain organoids, defined as complex 3D structures that develop from stem cells through a self-organization process, have emerged as powerful in vitro tools for studying brain cellular morphogenesis, structure and microenvironment (Pasca et al., 2022; Xu et al., 2021). In this work, we present an improved, robust and cost-effective model of mouse brain organoids (mBOs). Leveraging this protocol, we have generated mBOs to study the role of p53 and p73 isoforms in mouse neurodevelopment.
Our results reveal specific functions of p73 isoforms, with TAp73 governing transcriptional programs related to cell adhesion, extracellular matrix organization and epithelial cell morphogenesis, which are essential for establishing the archetypical proliferative ventricular-like zones in mBOs, resembling the mouse neurogenic niches. On the other hand, DNp73 regulates the early neuroectodermal fate switch and is thus crucial for the precise regulation of neural cell fate, as well as for preserving the undifferentiated and proliferative state of the NSC pool.
2 Materials and methods
2.1 Cell culture
Four mouse induced pluripotent stem cell (miPSC) lines were used (2 clones per line): (i) wild-type (WT), (ii) p73KO, (iii) p53 knockout (p53KO) and (iv) p73KO/p53KO (DKO)-miPSCs, whose generation and culture have been described previously (Martin-Lopez et al., 2017). Briefly, miPSCs were cultured on mouse embryonic fibroblast feeder cells (5 × 104 cells/cm2) in Dulbecco’s Modified Eagle Medium (DMEM, Cat# D5671, Sigma-Aldrich, MO, United States) supplemented with 15% fetal bovine serum (FBS, Cat# 17479633, Thermo Fisher Scientific, MA, United States), 2 mM L-glutamine (Cat# 7513, Sigma-Aldrich, MO, United States), 1 mM sodium pyruvate (Cat# 11360039, Thermo Fisher Scientific, MA, United States), 1 mM nonessential amino acids (Cat# 11140035, Thermo Fisher Scientific, MA, United States), 0.1 mM β-mercaptoethanol (Cat# M3148, Sigma-Aldrich, MO, United States) and 1,000 U/mL Leukemia Inhibitory Factor (LIF, Cat# ESG1107, Millipore, MA, United States).
Also, the following mouse embryonic stem cell (mESC) lines were used: (i) parental E14TG2α (Hooper et al., 1987) were kindly provided by Dr. Jim McWhir (former Researcher at the Roslin Institute, Edinburgh, Scotland, United Kingdom); (ii) E14-TAp73KO (3 clones) and (iii) E14-DNp73KO cells (3 clones), which were generated in our lab (Lopez-Ferreras et al., 2021). These cells were cultured on 0.1% gelatin-coated plates (105 cells/cm2) in Glasgow Minimum Essential Medium (GMEM, Cat# G-5154, Sigma-Aldrich, MO, United States) supplemented with 10% FBS, 2 mM L-glutamine, 1 mM sodium pyruvate, 1 mM nonessential amino acids, 0.1 mM β-mercaptoethanol and 500 U/mL LIF.
2.2 Mouse brain organoid generation
For mBOs generation, mPSCs were dissociated to single cells using 0.25% trypsin-EGTA and 8 × 103 cells were seeded and quickly aggregated in each well of 96-well ultra-low attachment plates (Cat# 7007, Costar, ME, United States), in differentiation medium (100 μL/well) containing GMEM supplemented with 10% KnockOut Serum Replacement (KSR, Cat# 10828010, Gibco, MA, United States), 2 mM GlutaMAX (Cat# 35050038, Gibco, MA, United States), 1 mM sodium pyruvate, 0.1 mM nonessential amino acids and 0.1 mM β-mercaptoethanol. A partial medium change (half of the volume) was performed on day 4, and 3 days later, cell aggregates were transferred to a 10-cm bacterial-grade dish in N2 medium: DMEM/F12 (Cat# D8437, Sigma-Aldrich, MO, United States) supplemented with 1:100 N2 (Cat# 17502–048, Gibco, MA, United States) and GlutaMAX. The medium was changed every 3 days until day 14. The day on which mPSCs were seeded in ultra-low attachment plates to differentiate is referred to as differentiation day 0.
The initial optimization of the protocol was performed using the E14TG2⍺ mESC line. We conducted several rounds of experiments, introducing stepwise modifications. These changes were systematically compared to the efficiency of the reference protocol developed by Eiraku et al. (2008), allowing us to iteratively refine and enhance our method. After optimization, we successfully generated mBOs across seven independent batches of mESCs. To further validate the robustness of our protocol, we extended its application to miPSCs, where we conducted experiments with four additional batches.
2.3 Tissue processing
mBOs were fixed in 4% paraformaldehyde (PFA) containing 6% sucrose at 4°C for 30 min. After fixation, mBOs were washed with PBS, cryoprotected overnight with 30% sucrose and stored at 4°C in PBS. Then, mBOs were embedded in PolyFreeze Tissue Freezing Medium (Cat# SHH0026, Sigma-Aldrich, MO, United States) and frozen in Peel-A-Way™ embedding molds (Cat# 18986-1, Polysciences, PA, United States) using dry ice. mBOs were finally cryosectioned at 8 μm using a Leica cryostat (Cat# CM 1950, Leica, Germany) and the obtained sections were collected in SuperFrost Plus™ Adhesion slides (Cat# J1800AMNZ, Epredia, NH, United States) and stored at −80°C.
2.4 Immunofluorescence
mBO cryosections were permeabilized in PBS containing 0.25% Triton X-100 and blocked in 10% of horse serum and 2% bovine serum albumin (BSA) for 2 h at room temperature (RT). mBO slices were incubated with primary antibodies overnight at 4°C in blocking solution, followed by the appropriate fluorophore-conjugated secondary antibodies incubation for 2 h. Nuclei were counterstained with DAPI.
Primary antibodies used were: anti-SOX2 (1:200, goat polyclonal, Cat# AF 2018, R&D Systems, MN, United States); anti-bIII tubulin (1:1,000, mouse monoclonal, Cat# MMS-435P, Covance, NJ, United States); anti-N-Cadherin (1:200, rabbit polyclonal, Cat# 13116, Cell Signaling, MA, United States); anti-NeuN (1:200, mouse monoclonal, Cat# MAB377, Merck Millipore, MA, United States); anti-phospho Histone H3 (pHH3) (1:500, mouse monoclonal, Cat# 9706, Cell Signaling); anti-Cleaved Caspase-3 (CC3) (1:200, rabbit polyclonal, Cat# 9664, Cell Signaling). The secondary antibodies were: Alexa568 anti-goat (Cat# A11057, Invitrogen, MA, United States); Alexa647 anti-mouse (Cat# A31571, Invitrogen); Alexa488 anti-rabbit (Cat# A21206, Invitrogen); Alexa647 anti-mouse IgG2a (Cat# A21241, Invitrogen); Alexa594 anti-mouse IgG1 (Cat# A21125, Invitrogen); all diluted 1:1,000.
2.5 Image acquisition, analysis and quantification
Confocal microscopy images (8 bits) were obtained in a Zeiss LSM800 Confocal Laser Scanning Microscope (Carl Zeiss Microscopy GmbH, Germany) using 25 × Plan-Apo/0.8 numerical aperture or 63 × Plan-Apo/1.4 numerical aperture oil objectives at RT. In some cases, a 0.5 × digital zoom was used. Confocal Z-stack images were acquired, and stacks were z-projected to maximum intensity. Images were processed with the ZEN blue software (Carl Zeiss Microscopy GmbH). All images were further analyzed using ImageJ/Fiji software (MD, United States). Three batches of mESCs-derived organoids and another three batches generated from miPSCs were used for quantifications. For each genotype, five different mBO cryosections per clone were analyzed across the three independent batches resulting in a total number of, at least, 15 mBOs per condition.
For the quantification of the number of ventricles, 25 × 0.5 Z-stack confocal images of mBO slices were employed to manually count the number of ventricles in each mBO cryosection based on the presence of SOX2+ neural rosettes with a defined N-CAD+ lumen.
For the quantification of NPCs, 63 × Z-stack confocal images were used and the number of SOX2+ and the total number of cells (DAPI-stained nuclei) were manually counted. Also, NeuN+ cells were manually quantified to calculate the percentage of NeuN+ cells.
The percentage of proliferating or apoptotic cells was calculated by analyzing 63 × Z-stack pictures to manually count the number of pHH3+ or CC3+ cells and the total number of cells (DAPI-stained nuclei). To quantify the percentage of proliferating progenitors and apoptotic progenitors or neurons, colocalization studies were performed using 63 × Z-stack confocal images. Images of the pHH3 or the CC3 channel were binarized using an automatic threshold. A selection was created (ROI) and applied to the SOX2 or TUJ1 channel. Then, the “Clear Outside” plug-in was applied to eliminate the signal without double staining. Finally, cells co-expressing pHH3/SOX2, CC3/SOX2 or CC3/TUJ1 were counted using ImageJ/Fiji.
2.6 RNA sequencing and transcriptome data analysis
Bulk RNA from WT-, TA- and DNp73KO-mBOs was isolated at the indicated time points (see figure legends), using RNeasy mini Kit (Cat# 74106, QIAGEN, Hilden, Germany). RNA concentration and quality were determined using a Nanodrop ND-100 (Thermo Fisher scientific, Waltham, MA, United States). For each genotype and timepoint, RNA samples from three biological replicates were sent to Novogene Company Limited (Cambridge, England, United Kingdom) for RNA-sequencing (RNA-seq). The library preparations were sequenced on an Illumina platform (Illumina Inc., CA, United States) generating paired-end reads with 40 million reads per direction. The quality of sequencing was validated with FASTQC software (Simon, 2010). The alignment and quantification of the reads on the reference sequences were performed using Salmon software (Patro et al., 2017). The differential expression analysis was conducted using the R statistical analysis program, with the DESeq2 package (Love et al., 2014). Mapping of the reads to the reference genome, normalization to obtain the FPKMs (fragments per kilobase of transcript per million mapped reads) and DESeq2 analysis were performed at the core service Bioinformática USAL (Universidad de Salamanca, Spain). The FPKM raw data is available at OPEN SCAYLE: https://open.scayle.es/dataset/alonso-olivares-et-al-2024. Functional annotation analysis on differentially expressed genes (DEGs) was carried out using The Database for Annotation, Visualization and Integrated Discovery (DAVID) (Huang da et al., 2009; Sherman et al., 2022). Principal component analysis (PCA) was performed using Clustvis web tool (Metsalu and Vilo, 2015). Venn diagrams were created using BioVenn (Hulsen et al., 2008).
2.7 Gene expression analysis by reverse transcription and quantitative PCR (qRT-PCR)
cDNA synthesis from RNA samples was carried out using the High Capacity RNA-to-cDNA kit (Cat# 4387406, Applied Biosystems, Waltham, MA, United States). Gene expression was analyzed by qRT-PCR using FastStart Universal SYBR Green Master (Cat# 4913850001, Roche, Basel, CH) in a StepOnePlusTM Real-Time PCR System (Applied Biosystems). Primer sequences for target genes used are included in Supplementary Table S1. The comparative threshold cycle method was used to quantify relative mRNA expression.
2.8 Statistical analysis
Statistical analyses were performed using GraphPad Prism software version 10.2.3 for Windows (GraphPad Software, Boston, MA, United States). Values were expressed as mean ± standard error of the mean (SEM). Data distribution was analyzed in all cases. When comparing two groups, an unpaired Student’s t-test was used. For comparisons between multiple groups, statistical differences were determined using either parametric one-way ANOVA followed by Tukey’s correction or nonparametric Kruskal–Wallis tests with Dunn’s correction. Differences were considered significant when p < 0.05. Comparisons with WT mBOs are indicated by asterisks (*p < 0.05, **p < 0.01, ***p < 0.001), whereas hash (#) denote comparisons with p73KO organoids and dagger (†) is used to show comparisons with p53KO. ns: non-significant.
3 Results
3.1 Optimization of a rapid and robust system to generate mBOs from mouse pluripotent stem cells (mPSCs)
Several protocols to generate mBOs have been described (Eiraku et al., 2008; Nasu et al., 2012; Watanabe et al., 2005), but they are not highly reproducible and require some adaptation to each cell line. We sought to establish an improved protocol to consistently generate high-quality, reproducible mBOs. We first used a previous method (Eiraku et al., 2008) to generate mBOs from WT-mPSCs, both mESC and miPSC lines. With this protocol, some mESCs differentiated into the neural lineage and gave rise to TUJ1+ neurons, but neural progenitors (SOX2+) failed to self-organize and form N-CAD+ ventricle-like structures (Supplementary Figure S1). We then adjusted some parameters to optimize the initial size of the organoids and improve cell viability: we increased the number of seeded cells from 3,000 to 8,000 cells/well and reduced the frequency of culture medium changing during the process (Figure 1A). We also modified the composition of the initial differentiation medium, eliminating Dickkopf Wnt Signaling Pathway Inhibitor 1 (DKK-1) and left-right determination factor 1 (LEFTY-1). DKK-1 inactivates Wnt/β-catenin signaling, which is necessary for NSC proliferation (Chenn and Walsh, 2002). On the other hand, LEFTY-1 inhibits Nodal signaling through pSMAD2/3 (Kim et al., 2014), so the presence of LEFTY-1 in the medium may prevent differentiation. Based on this reasoning, and in agreement with (Eiraku and Sasai, 2011), we eliminated both factors from the initial differentiation medium.
Figure 1. Optimization of a robust system to generate mBOs from mouse pluripotent stem cells. (A) Schematic overview of the optimized protocol for generating mBOs from mPSCs. (B, B’) Representative confocal images of mBOs immunostained for SOX2 (neural progenitor marker), TUJ1 (neuronal marker) and N-CAD (apical location in cells forming the ventricles) at day 14 of the differentiation process from the mESC line E14TG2α (B) and WT-iPSCs (B’); white arrows indicate ventricle-like structures. The lower panels show a higher magnification image of the squared regions in the upper panels. Scale bars: upper panel: 250 μm; lower panel: 50 μm. (C) Functional annotation of DEGs upregulated between day 14 and day 3 WT-mBOs; FDR: false discovery rate. (D) Heatmap profiles showing RNA-seq expression Z-scores for representative gene set modules of DEGs among WT-mBOs at different time points. Representative gene names within the indicated functional module are shown. (E) Time course analysis of gene expression (FPKM values from RNA-seq) for genes related to the indicated categories: stem cells, intermediate progenitor cells, neurons, cortical plate, caudal brain and cerebellum.
To determine whether the modified protocol could generate bona fide mBOs, we analyzed the expression and localization of key markers: SOX2 (a neural stem cell marker), TUJ1 (an early neuronal marker), and neural-cadherin N-CAD (a cell adhesion molecule). These markers are essential for assessing the presence and organization of VZ-like neural rosettes, a crucial indicator of proper brain organoid formation and development (Eichmuller and Knoblich, 2022). Indeed, rosette-like structures represent the in vitro counterpart of the neural tube, the structure of the early nervous system in vivo (Roll et al., 2022). We observed that WT mESC-derived mBOs contained SOX2+ neural progenitors at day 14, which were clearly organized in rosette-like structures surrounded by TUJ1+ neurons (Figure 1B). Rosettes were structurally similar to the VZ of the embryonic forebrain, with an internal lumen nicely decorated with N-CAD corresponding to the apical domain of neuroepithelial cells. This protocol was also efficient at generating mBOs from WT miPSCs (Figure 1B’), again with rosettes structurally similar to telencephalic VZ, thus confirming that the implemented modifications resulted in a robust protocol to generate mBOs from mPSCs.
To substantiate the quality of the mBOs generated with the improved protocol, we performed bulk transcriptomic analyses on mBOs derived from WT E14TG2α mESCs, analyzing them at different time points during development. Analysis of differentially expressed genes (DEGs) identified 3,083 genes that were significantly upregulated from day 3 to day 14 (padj <0.01; log2FC > 1), while 914 genes were downregulated (padj <0.01; log2FC < −1). Functional annotation analysis of the upregulated DEGs confirmed that our culture conditions induced transcriptional programs related to neurogenesis, highlighting Gene Ontology (GO) terms related to nervous system development, axon guidance, neuron fate specification, and forebrain neuron differentiation (Figure 1C). The transcriptomic data also showed progressive upregulation of pan neuronal markers such as Map2, NeuN, Ncam1 and Dcx, as well as genes related to axon guidance (Dpysl4 and Dpysl5, Robo2) and neurite outgrowth (Ncan) during the process of mBO formation (Figure 1D; Supplementary Table S2). Upregulation of some of these DEGs associated with neuronal identity was further validated by qRT-PCR (Supplementary Figure S2). This upregulation was particularly striking in genes related to the organization of the cortical plate, such as Foxg1 and Lhx2, whereas markers of the caudal brain (Gbx2 and Lmx1a) or cerebellum (Pcp2) remained low (Figure 1E). Concomitantly, pluripotency genes such as Pou5f1, Zscan10 or Lin28A were downregulated from day 3 to day 14 (Figures 1D, E; Supplementary Figure S2; Supplementary Table S2), together with genes associated with cell cycle progression and DNA synthesis, indicating a decrease in stem and proliferative cells. Some of the most prominent signaling pathways involved in neural induction and development include Notch, Sonic Hedgehog and Wnt (Manzari-Tavakoli et al., 2022). As expected, these signaling pathways were activated during the differentiation process and generation of mBOs, as shown by the upregulation of genes directly involved in signal transduction (Shh and Notch1), as well as some of their target genes (Emx2, Ascl1 or Neurogenin), involved in neural development (Figures 1C–E).
Altogether, this confirmed that our optimized protocol promoted neural development and differentiation, leading to the generation and organization of mature cortical neurons from mPSCs. Therefore, we have fine-tuned and improved a protocol to generate mBOs by implementing changes in the initial cell number and the cocktail of factors used to induce neural commitment. This resulted in a rapid, affordable and highly reproducible protocol for the generation of mBOs from both induced and embryonic mPSCs.
3.2 p73 and p53 regulate neural progenitor proliferation and neuronal apoptosis in mBOs
p53 and p73 are known cell cycle and cell death regulators during brain development (Agostini et al., 2018). Thus, we decided to corroborate these facts in the mBO model and to exploit our improved protocol to delve into the role of p53 and p73 in several aspects of mouse cortical brain development. We used previously generated p73KO-miPSCs, together with p53KO- and DKO-miPSCs (Martin-Lopez et al., 2017). Additionally, we made use of isoform specific mESCs (TAp73KO- and DNp73KO-mESCs) generated by CRISPR/Cas9 genome editing (Lopez-Ferreras et al., 2021). First, we assessed the capacity of these mBOs to recapitulate some of the effects of p53 and p73 deficiency on neural cell proliferation and/or cell death described in vivo (Gonzalez-Cano et al., 2010; Meletis et al., 2006). Similarly to what was observed in neurospheres assays (Meletis et al., 2006; Armesilla-Diaz et al., 2009), histological analysis of p53KO-mBOs after 14 days in differentiation conditions revealed a significant increase in the percentage of proliferative cells (positive for the mitosis marker pHH3) compared to controls (WT-mBOs) (Figure 2A). On the contrary, the absence of p73 led to a significant reduction in the percentage of proliferating cells. Interestingly, elimination of p73 in the context of p53-deficiency (DKO-mBOs), restored the proliferation index to levels close to WT-mBOs, suggesting that p73 is essential to sustain the high proliferation index unleashed by the absence of p53 within the mBO (Figure 2A).
Figure 2. p73 and p53 regulate proliferation of neural progenitors and neuronal apoptosis within mBO in opposite but coordinated ways: TAp73 supports neural progenitor proliferation whereas DNp73 prevents p53-dependent neuronal apoptosis. All images are representative immunostaining confocal microscopy images of mBO cryosections at day 14 from miPSCs (A–D) or mESCs (E–H) with the indicated genotypes. (A, E) Immunostaining and quantification of proliferating cells (pHH3+). (B, F) Immunostaining and quantification of proliferating neural progenitors (pHH3+/SOX2+). (C, G) Immunostaining and quantification of apoptotic cells (CC3+). (D, H) Immunostaining and quantification of apoptotic neural progenitors (CC3+/SOX2+) or neurons (CC3+/TUJ1+). Scale bar: 50 μm. Data were collected from at least 30 mBOs from 2 different miPSC-clones or 3 different mESC-clones across 3 independent experiments and are represented as mean ± SEM. Statistical differences using Kruskal–Wallis test are indicated as: *, #, †, p-value <0.05; **, ##, ††, p-value <0.01; ***, ###, †††, p-value <0.001.
To identify if the reduced proliferation observed in p73KO-mBOs affected their overall cellular populations, or if it was restricted to a specific cell type, we quantified the percentage of proliferating SOX2+ NPCs (pHH3+/SOX2+, Figure 2B). In WT-mBOs, the percentage of proliferating NPCs was around 2%, while this percentage drastically decreased to 0.45% in p73KO-mBOs (Figure 2B). Again, p53KO-mBOs doubled the percentage of SOX2+ proliferating cells, while DKO-mBOs had NPC proliferation rates comparable to their WT counterparts (Figure 2B). It is important to keep in mind that the decrease of proliferating cells in the p73KO could be related to a compensatory p53-activation leading the NPCs to exit the cell cycle, to senesce or to die.
Next, we analyzed apoptosis by quantification of CC3+ cells. In WT-mBOs, nearly 10% of cells were CC3+ (Figure 2C). Interestingly, while the absence of p53 in p53KO-mBOs significantly reduced, but not obliterated, this basal level of apoptosis, lack of p73 increased significantly the percentage of apoptotic cells in p73KO-mBOs compared to the WT-mBOs (Figure 2C). This effect was reduced to WT levels upon p53 loss in DKO-mBOs (Figure 2C), confirming a role of p73 in preventing p53-dependent cell death. However, it is noteworthy that a significant percentage of apoptotic cell death remains in double mutants, compared to mBOs lacking only p53, suggesting that p73 deficiency elicits a certain level of p53-independent apoptosis.
To determine whether p73 function is necessary to prevent overall apoptosis within mBOs or if its role is specific to a particular cell type, we performed double staining for SOX2/CC3 and TUJ1/CC3. As shown in Figure 2D, the lack of p73 did not significantly affect the frequency of apoptotic NPCs (SOX2+/CC3+ cells) but led to a sharp increase in the number of apoptotic neurons (TUJ1+/CC3+ cells).
Next, we sought to discriminate which p73 isoform was responsible for maintaining the proliferative pool of NPCs by analyzing cell proliferation and apoptosis in the p73 isoform-specific knockout-mBOs. In WT-mBOs, over 1% of the cells were pHH3+, a percentage severely reduced in mBOs lacking either TA- or DNp73 (Figure 2E), indicating that both isoforms are required to maintain the mitotic rate within the mBOs. Intriguingly, TAp73 deficiency led to a significant decrease in the number of SOX2+ NPCs, while DNp73 deficiency did not have a significant effect (Figure 2F). Together, these results indicate that TAp73 is the isoform necessary to sustain proliferation and maintain the pool of NPCs, whereas DNp73 is dispensable for the latter function, but essential for sustaining the proliferative capacity of other cell types in the mBO. The analysis of CC3+ cells revealed that while the loss of TAp73 had no effect on the overall apoptosis rate compared to WT-mBOs, DNp73-deficiency significantly increased the number of apoptotic cells (Figure 2G). Moreover, the lack of DNp73 primarily affected early immature neurons, but not SOX2+ progenitors (Figure 2H, right and left panels, respectively).
Collectively, our mBO model confirmed that p73 plays a key role in neuronal survival in vitro, consistent with previous observations in vivo (Tissir et al., 2009), where DNp73 plays a pro-survival role in discrete neuron types in a p53-dependent manner (Tissir et al., 2009; Osterburg and Dotsch, 2022; Pozniak et al., 2002; Pozniak et al., 2000). Moreover, our results evidenced that both TA- and DNp73 are necessary to maintain a functional pool of NPCs within the mBO, but they operate differently depending on the cell type: TAp73 maintains the proliferating NPC pool, while DNp73 sustains general proliferation and prevents neuronal apoptosis.
3.3 p53 and p73 are required for the self-assembly of mBOs
We next examined the effect of total loss of p73 and/or p53, as well as the elimination of the specific p73 isoforms, in the cellular and molecular events of mouse brain morphogenesis, beyond the already described impact in cell proliferation and apoptosis. mBOs at day 14 displayed obvious morphological differences between the four genotypes (Figure 3A). WT-mBOs exhibited multiple buds with translucent lumens resembling neuroepithelia (Figure 3A, yellow arrows) whereas the lack of p73, p53 or both (DKO) resulted in mBOs with less visible neuroepithelial structures, especially in the case of p73KO-mBOs which had almost smooth borders (Figure 3A, dashed arrows). This suggested that deficiencies in these p53 family members interfere with the correct development and morphogenesis of mBOs. A more detailed analysis of these epithelia-like buds revealed that WT-mBOs exhibited rosette-like structures formed by SOX2+ NPCs with an N-CAD+ ventricular surface (Figure 3B, arrows), surrounded by TUJ1+ neurons on the organoid surface, resembling the cortical plate in developing embryos (Roll et al., 2022). In contrast, the existing ventricles in mBOs lacking p73, p53, or both, were ill-defined, with an overall disrupted pattern (Figure 3B, dashed arrows). This was particularly striking in p73KO-mBOs, where neither N-CAD expression, nor defined ventricles, were observed (dashed arrows). Accordingly, while WT-mBOs presented an average of three defined neural ventricles per organoid cryosection, the mutant mBOs, and the p73KO-mBOs in particular, displayed a reduced number of ventricles (Figure 3B’). These observations suggest defects in cell-cell adhesion and cellular self-assembly, highlighting the essential requirement of Trp73 for proper NPC organization and self-assembly into defined VZ-like structures.
Figure 3. TAp73 is necessary for the correct organization of neural progenitors into ventricle-like structures, while DNp73 is dispensable. (A, C) Phase contrast images of mBOs derived from different miPSCs (A) or mESCs (C) genotypes. Orange arrows indicate buds with transparent lumens resembling neuroepithelium. Scale bar: 250 μm. (B, D) Representative confocal images of mBOs immunostained for SOX2 (neural progenitor marker), TUJ1 (neuronal marker) and N-CAD (apical location in cells forming the ventricles) at day 14 of the differentiation process. Dashed-squared areas on the left panels are magnified on the right panels; white arrows indicate well-formed ventricular structures; white dashed arrows indicate the absence of defined ventricular structures. Scale bars: 250 μm (left panel) and 50 μm (right panel). (B’, D’) Quantification of the number of ventricles per cryosection in different mBOs. The number of ventricles per cryosection is represented as mean ± SEM. Data were collected from at least 20 mBOs of the indicated clones from 3 independent experiments. Statistical differences using Kruskal-Wallis test are represented as: * p-value <0.05; *** p-value <0.001.
3.4 TAp73 is necessary for the organization of NPCs into VZ-like structures
Having identified that Trp73 is fundamental for mBO organization and the formation of VZ-like structures, we aimed to disclose the specific roles of the different p73 isoforms in this process. While WT- and DNp73KO-mBOs presented translucent neuroepithelial buds (Figure 3C, yellow arrows), mBOs lacking TAp73 were opaque, indicating defective formation of optically translucent and radially organized neuroectoderm. Besides, the neuroepithelial structures were less numerous and conspicuous in TAp73KO-mBOs (Figure 3C, dashed arrows).
Immunofluorescence analysis of mBO cryosections showed that mBOs generated from WT-mESCs exhibited NPCs organized into spherical rosettes with defined N-CAD+ ventricles, similar to the ones generated from WT miPSCs (Figure 3D, arrows). In contrast, NPCs in TAp73KO-mBOs were rather scattered and loose and did not form neural rosettes nor N-CAD+ ventricles (Figure 3D, dashed arrows), suggesting defects in NPC cell-cell adhesion, polarization and self-organization. In contrast, DNp73KO-mBOs showed well defined ventricles (Figure 3D, arrows), which were even more numerous than in WT-mBOs (Figure 3D’). This increased abundance of rosettes might be due to elevated TAp73 levels detected in these cells (Lopez-Ferreras et al., 2021), which could enhance their ability to establish correct cell polarity and form cell-to-cell interactions. In this differentiation context, this may favor the formation of VZ-like structures.
Considering all the above, the defective organization observed in p73KO-mBOs appears to be primarily due to the absence of the TAp73 isoform, seemingly essential for the organization of NPCs into rosette-like structures, and hence for the proper formation of ventricles within the mBO.
3.5 p73 isoforms, particularly DNp73, play a prominent role in the differentiation of NPCs
In the course of brain morphogenesis, NSCs differentiate to progressively acquire a specific neuronal identity and to eventually establish the complex structures of the brain. The processes of fate commitment and spatial organization are tightly regulated and, in fact, closely coupled. Thus, we explored whether the observed defects in neural structure organization within the mBOs were related to irregularities in the process of cellular decision-making and cell fate determination. In a typical differentiation process, SOX2+ progenitors give rise to TUJ1+ neurons which further differentiate into mature, post-mitotic neurons expressing the nuclear marker NeuN. Hence, we analyzed the expression of SOX2 and NeuN in mBOs at day 14 by immunofluorescence (Figure 4). In WT-mBOs, nearly 40% of cells were SOX2+ neural progenitors and 14% were NeuN+ neurons (Figures 4A, A’). The number of SOX2+ progenitors in mBOs lacking p53 was not significantly affected, while they exhibited a minor, but significant, increase (up to 23%) in the percentage of NeuN+ neurons (Figure 4A’), indicating that lack of p53 slightly affected the differentiation programs involved in neuronal maturation rather than the generation of neurons from SOX2+ progenitors. This is in agreement with previous studies highlighting the role of p53 in neuronal differentiation and maturation (Armesilla-Diaz et al., 2009; Forsberg et al., 2013; Marin Navarro et al., 2020). In contrast, mBOs lacking p73 showed significantly fewer SOX2+ neural progenitors and a noticeable significant increase in mature NeuN + cells (from 14% to 65%) (Figure 4A’). This indicates that the absence of p73 is associated with deregulation of neural cell fate, leading to a premature neuronal differentiation in these organoids, supporting previous data obtained using neurosphere assays (Fujitani et al., 2010; Gonzalez-Cano et al., 2010). Interestingly, the elimination of p53 in this context (DKO), decreased significantly, but not completely, this effect, suggesting a functional interaction between p73 and p53 in regulating cell fate and neural differentiation.
Figure 4. DNp73 absence leads to a depletion in the pool of neuronal precursors with a concomitant premature neuronal differentiation in a p53-dependent manner (A, A’) Representative confocal images (A) and quantification (A’) of mBOs generated from miPSCs at day 14, immunostained for SOX2 (neural progenitor marker), and NeuN (marker of postmitotic neurons). Dashed-squared areas on the left panel are magnified on the right panel. (B, B’) Representative confocal images (B) and quantification (B’) of mBOs generated from mESCs at day 14, immunostained for SOX2 (neural progenitor marker), and NeuN (marker of postmitotic neurons). Dashed-squared areas on the left panel are magnified on the right panel. Scale bars: 250 μm (left panel) and 50 μm (right panel). For the bar graphs in (A’) and (B’), data were collected from at least 30 mBOs of 2 different miPSC-clones or 20 mBOs from 3 different mESC-clones across 3 independent experiments and are represented as mean ± SEM. Statistical differences using Kruskal–Wallis test are represented as: ***, ### p-value <0.001.
When analyzing the contribution of each p73 isoform, we observed that TAp73 deficiency did not significantly affect the number of SOX2+ progenitors nor the percentage of NeuN + neurons compared to WT-mBO (Figures 4B, B’). In contrast, the absence of DNp73 led to a notable decrease in the number of SOX2+ progenitors, accompanied by a highly significant increase in the number of mature neurons (from 8% in WT-mBOs to 50% in DNp73-mBOs). In DNp73-mBOs, these neurons were detected throughout the organoid rather than mainly in the cortical border, as was the case in WT- and TAp73KO-mBOs (Figure 4B, white arrows). This is in line with the phenotype observed in p73KO-mBOs and confirms that accelerated neuronal differentiation in these organoids is caused by the absence of DNp73. These results demonstrate that DNp73 has a remarkable effect in orchestrating neuronal differentiation from neural progenitors and may specifically act as a negative regulator of neuronal cell fate.
3.6 TAp73 and DNp73 exert their function through the regulation of distinct transcriptional programs
To further pinpoint the transcriptional changes occurring in the p73 isoform-deficient mBOs, we conducted global transcriptomics analyses using bulk RNA sequencing at days 7 and 14 (Supplementary Table S3). PCA analysis showed that clustering of biological replicates from WT, TA- and DN-p73KO organoids improved over time, with a clear separation between days 7 and 14 (Supplementary Figure S3A). Moreover, the major source of variation was linked to the mBOs genotype (with principal component PC1 accounting for almost 40% of the variance). The specific loss of either TA- or DN-p73 resulted in altered transcriptional profiles (Supplementary Figure S3B), with the effect being more pronounced at later stages. Moreover, DNp73 deficiency had a greater impact on the number of DEGs, particularly at day 14: 2,835 genes were upregulated (padj <0.01; log2FC > 1) and 2,350 genes became downregulated (padj <0.01; log2FC < −1) compared to WT-mBOs.
TAp73 is a bona fide transcriptional activator, while DNp73 acts mostly, but not always (Toh et al., 2008), as a repressor. Thus, we focused on genes that were downregulated in the absence of TAp73 and upregulated upon DNp73 deficiency. Functional annotation analyses of the DEG lists revealed that downregulated genes in TAp73KO-mBOs were associated with cell adhesion, extracellular matrix organization, animal organ morphogenesis (e.g., heart tube development, mammary gland development, etc.) and axis specification (Figure 5A). This reinforces the role of TAp73 in regulating a transcriptional module related to tissue morphogenesis (Maeso-Alonso et al., 2021) and may explain the observed inability to form archetypal structures, such as proper VZ-like structures, in the absence of TAp73. Notably, genes associated with functions previously linked to TAp73, such as angiogenesis, stem cell differentiation, and cytoskeleton regulation (Gonzalez-Cano et al., 2010; Fernandez-Alonso et al., 2015; Fuertes-Alvarez et al., 2018), were also enriched. On the other hand, upregulated genes in DNp73KO-mBOs were highly significantly associated to biological terms such as nervous system development, synaptic transmission, axonogenesis or axon guidance (Figure 5B), correlating with the phenotype of premature neuronal differentiation observed in DNp73KO-mBOs. Additionally, GO terms linked to DNA replication, cell cycle progression, cell division and stem cell population maintenance were enriched among downregulated DEGs in DNp73KO-mBOs (Supplementary Figure S3C), indicating that DNp73 sustains the capacity of NPCs to proliferate.
Figure 5. TA- and DNp73 regulate distinct transcriptional programs related to brain development and nervous system function. (A, B) DEGs with a p-adj <0.05 were used for GO analysis using DAVID Software. (A) Functional annotation of downregulated DEGs in TAp73KO-mBOs compared to WT after 14 days in culture. (B) Functional annotation clustering of upregulated DEGs in DNp73KO-mBO compared to WT after 14 days in culture. (C) Venn diagrams illustrating the percentage of transcripts related to neural stem cells (NSCs), neural progenitor cells (NPCs) and mature neuron profiles which are upregulated (red) or downregulated (purple) in TAp73KO or DNp73KO-mBOs at day 14. (D) Heatmaps of DEGs related to different neuronal identities in TAp73 and DNp73KO-mBOs at day 14. (E) Heatmaps of DEGs related to metabolism, axon growth, synaptic functionality, signal transduction and oxidative metabolism in TAp73KO and DNp73KO-mBOs at day 14.
To evaluate how the elimination of TA- or DN-p73 affected neural cell fate determination, we intersected the up- and downregulated DEGs at different stages (compared to their WT counterparts), with gene sets that have been previously linked to NSCs, NPCs and mature neurons (Bifari et al., 2020; Cahoy et al., 2008; Ciarpella et al., 2021). As shown in Figure 5C, none of the genes related to the NPCs- or the mature neurons-profiles were downregulated in DNp73KO-mBOs, and only a few were repressed in TAp73KO. It called our attention that three NSC-associated genes were commonly downregulated in DNp73KO and TAp73KO-mBOs both at day 7 and 14: NADPH-Dependent Carbonyl Reductase 3 (Cbr3), Thiosulfate Sulfurtransferase (Tst) and ETS Variant Transcription Factor 4 (Etv4). The downregulation of Etv4 is of particular interest, since it works as a mechanical transducer that drives spatiotemporal lineage specification and its inactivation in human embryonic stem cell epithelia derepresses the neuroectoderm fate (Yang et al., 2024). On the other hand, DNp73KO-mBOs remarkably upregulated genes related to the neural fate when compared to the WT counterparts, whereas the lack of TAp73 had a milder outcome. Indeed, elimination of TAp73 had almost no impact at early stages (day 7) (Figure 5C). On the contrary, the effect in DNp73KO-mBOs was striking at early stages, with DNp73 deficiency resulting in the upregulation of 60% of the genes included in the NPC profile and a 37% overlap with genes related to mature neurons at day 7. This suggests that DNp73 absence accelerates neuronal lineage induction.
The importance of DNp73 repression, and TAp73 to a minor extent, in establishing neuronal identity is further illustrated by the significant upregulation of individual markers related to various neuron types (Figure 5D). For instance, isoform-deficient mBOs upregulate genes that have been shown to be important for the function of interneurons, such as Dlx1, Scn1a or Erbb4 (Ciarpella et al., 2021). In addition, DNp73KO-mBOs robustly overexpress genes related to excitatory and inhibitory synapses. The enrichment in glutamatergic synapses was evidenced, for example, by the upregulation of the glutamate transporter Slc17a6 (VGLUT2), as well as the overexpression of NMDA and AMPA receptors (Grin2b, and Gria1a/Gria2, respectively). On the other hand, the presence of GABAergic synapses was indirectly assessed by the increased expression of GABA transporter (Slc32a1) and gephyrin (Gphn). Conversely, markers indicative of other neuronal lineages (e.g., serotonergic) were not consistently expressed. Collectively, our results indicate the presence of a diversely enriched population of neurons in p73 deficient organoids. Finally, we examined a previously described gene set of mature neuron markers (Ciarpella et al., 2021). The expression levels of genes related to: i) metabolism and axon growth, ii) synaptic functionality and vesicle trafficking and iii) signal transduction and neurotransmission regulation were increased in both genotypes, most notably in DNp73KO-mBOs, indicating an enhanced neuronal maturation (Figure 5E). However, genes related to beta-oxidation pathways at peroxisomal (Acox1, Ech1) or mitochondrial level (Acadsb, Echdc2, and Slc25a20) (Watanabe et al., 2007) were not significantly upregulated, suggesting that at the latest analyzed time point the metabolic switch to oxidative metabolism did not accompany the neuronal maturation observed in those organoids.
Brain cortex development is a highly orchestrated sequence of events that requires precise signaling inputs. Shh signaling is vital for the formation and maturation of cortical layers, as it regulates the timing of progenitor cell differentiation and promotes their proliferation (Cai et al., 2023; Saade et al., 2017). These processes affect the thickness and cell composition of each cortical layer (Belmonte-Mateos and Pujades, 2021; Nowakowski et al., 2017). The Notch signaling pathway is pivotal in maintaining NPC populations and regulating their differentiation into specialized neurons and glia (Gaiano and Fishell, 2002), thereby influencing the overall size and complexity of the cortical region within the organoid. Additionally, the timing of Wnt activation is particularly crucial for establishing neural identity and suppressing unwanted cell fates (Rosso and Inestrosa, 2013), ensuring that cortical brain organoids develop with the appropriate cellular composition and organization. Supporting the role of DNp73 in regulating neural cell fate, we observed significant overexpression of some key genes involved in these pathways in DNp73KO-mBOs at day 7 (Figure 6A). This was evidenced by the upregulation of key genes involved in the activation and transduction of these pathways (Notch1, Shh, and Wnt5a), as well as their target genes (Emx2, Ascl-1, and Neurog1) (Figures 6A, B). Additionally, the induction of the neural lineage requires dual SMAD inhibition (Chambers et al., 2009). In accordance, DNp73KO-mBOs also exhibited downregulation of Bmp-4 and Nodal, which promote the activation of Smad1/5/8 and Smad2/3 signaling, respectively, suggesting that the absence of DNp73KO leads to dual-SMAD inhibition.
Figure 6. Lack of DNp73 results in the enhanced activation of the Shh, Wnt, and Notch signaling pathways, along with an early and strong repression of BMP and Activin/Nodal signaling. (A) Bar plot showing the expression of a series of transcripts corresponding to the indicated signaling pathways. The log2 Fold Change of DNp73 mBO at day 14 compared to WT mBOs at day 14 is shown. (B) Time course analysis of gene expression corresponding to genes related to the indicated signaling pathways. (C) Expression levels of the indicated transcripts were analyzed in WT and DNp73KO-mBOs at days 3, 7, and 14 of the optimized protocol using quantitative PCR (qPCR). The data represent the mean ± SEM of four replicates, derived from two independent experiments. Statistical significance was assessed using Student’s t-test and is denoted as follows: * p < 0.05; ** p < 0.01; *** p < 0.001.
Gene expression dynamics in DNp73KO-mBOs from day 3 to day 14, as indicted by our transcriptomic data and validated by qRT-PCR analysis, reflected the activation of the Shh, Wnt, and Notch signaling pathways (Figures 6B, C). Regarding the Activin/Nodal signaling, it is known that this cascade regulates Nodal and Lefty1 gene expression and, that while Lefty1 blocks Nodal signaling (Chen and Shen, 2004), Lefty1 expression pattern follows essentially that of Nodal in early stages of development (Besser, 2004; Zhang et al., 2019). Thus, in agreement with a repression of this signaling cascade, we detected a significant reduction of Lefty1 expression in DNp73KO-mBOs from days 7–14 (Figures 6B, C). This repression of Nodal signaling could account for the switch towards a predominant neuroectodermal fate observed in the absence of DNp73 at early stages (Figure 5C), as previously reported in ESCs (Vallier et al., 2004; Wang et al., 2017), pointing to DNp73 as a node regulator of neuronal fate. Altogether, these characteristic molecular signatures highlight the importance of p73, particularly DNp73, for a precise and timely regulation of the onset of neural cell fate commitment and shed new light on the specific role of each p73 isoform, and how they tightly coordinate the maintenance of a pool of proliferating neural progenitors with the timely generation of neurons, and their maturation.
4 Discussion
Generating cellular models that accurately reflect the complexity and uniqueness of the mammalian brain is a major challenge in cell biology. In this work, we refined a 3D in vitro system to generate mBOs that recapitulate the main structural and developmental characteristics of the mouse cerebral cortex. We introduced several modifications into the protocol previously defined by Eiraku and coworkers (Eiraku et al., 2008), aimed at optimizing the initial number of cells to constitute the aggregates and increase their survival. These changes resulted in a rapid, affordable and robust method that efficiently induced the differentiation of pluripotent stem cells, mESCs and miPSCs, into neural progenitors and neurons, which constituted distinctive 3D VZ-like neurogenic structures. As it has been suggested before (Eiraku et al., 2008; Watanabe et al., 2005; Eiraku and Sasai, 2011; Watanabe et al., 2007; Gaspard et al., 2008), the onset of neural differentiation is predominantly dependent on cell-autonomous mechanisms rather than on inductive signaling. This illustrates why the process of mBO generation is exquisitely influenced by minor changes in the initial cell density and basal culture conditions (Eiraku et al., 2008; Watanabe et al., 2005; Wataya et al., 2008).
Murine cerebral organoids have also been generated from NSC (Ciarpella et al., 2021). The protocol developed by Ciarpella and colleagues efficiently generates brain organoids composed of mature and functional neurons after 32 days in culture. Although these authors reported a very interesting tool for screening purposes, the fact that they used mouse embryos as a source of NSCs, instead of mESCs or miPSCs, makes the protocol more laborious, limits the throughput and still relies on the use of animals compared to ours. In addition, our protocol is shorter (14 days) and requires less exogenous factors, which determines a faster and cheaper production of neuron-bearing organoids that can be useful to study the early developmental stages of mouse brain development. Therefore, despite previous attempts of generating 3D forebrain structures from mPSCs (Eiraku et al., 2008) and more recent achievements in obtaining these structures from NSCs (Ciarpella et al., 2021), the protocol described here reports the efficient production of mBOs from PSCs in a fast, affordable and robust manner.
Applying our optimized protocol to analyze the function of p53 family members in mouse brain morphogenesis, we have been able not only to confirm some biological findings obtained from other in vitro and in vivo models (Agostini et al., 2010; Gonzalez-Cano et al., 2010; Talos et al., 2010), but also to provide a more comprehensive reading of the distinct roles played by p73 isoforms in this particular context. Our data showed that p73 deficiency results in a dramatic reduction in NPCs, likely due to a combination of premature neuronal differentiation and decreased proliferation capacity of neural progenitors. This supports a model in which the progenitor population fails to properly expand, exits the cell cycle prematurely, and initiates neuronal differentiation. Previous in vivo studies from our group have shown that p73 deficiency halts the transition of RGCs into ependymal cells (Gonzalez-Cano et al., 2016). RGCs function as neurogenic progenitors during brain development, and therefore, the previously defined role of p73 in RGCs is consistent with our in vitro findings regarding neural progenitor differentiation.
During development, neurogenic progenitors must polarize and establish cell-cell adhesions to organize in ventricular structures. Disruptions in these functions in vivo can lead to structural defects in the neurogenic niches, including the ependymal cell layer that lines and defines the VZ (Gonzalez-Cano et al., 2016). In addition, p73 has been implicated in regulating planar cell polarity (PCP) in different contexts, including the neurogenic niche of the SVZ of the brain (Gonzalez-Cano et al., 2016), as well as in the proper establishment of adherens and tight junctions in different cell types (Maeso-Alonso et al., 2021; Lopez-Ferreras et al., 2021; Holembowski et al., 2014). This endorses the results obtained in mBOs, in which we report that p73 deficiency hinders the formation of ventricular structures early during development.
The differential contribution of the p73 isoforms to the regulation of neurogenesis and brain development has been a longstanding question, as the phenotypes observed in the TA- and DN-specific knockout mice are intriguing. Trp73KO mice, that lack all p73 isoforms, are runty and exhibit very severe brain anomalies such as hydrocephalus and hippocampal dysgenesis (Yang et al., 2000). TAp73KO mice partially recapitulate this phenotype and present hippocampal dysgenesis too, whereas DNp73KO show a mild phenotype characterized by neurodegeneration (Tissir et al., 2009; Wilhelm et al., 2010). To distinguish the molecular function of the p73 isoforms, we took advantage of the mPSCs previously generated in our group, in which TA- and DNp73 isoforms were specifically ablated (Lopez-Ferreras et al., 2021), and generated mBOs. We found that TAp73 is essential for precise structural organization of neural progenitors. Also, our transcriptome analysis shows that TAp73 regulates a transcriptional module involved in cell adhesion, extracellular matrix organization and morphogenesis, which mechanistically underlies the establishment of the structures which are found disrupted in TAp73KO-mBOs. On the other hand, DNp73KO-mBOs exhibit premature neuronal differentiation, which seems counterintuitive with the reported role of DNp73 as a pivotal regulator of neuronal survival (Tissir et al., 2009). However, the neurons observed in mBOs are highly apoptotic, which is indeed consistent with the previously in vivo observation.
The increased apoptosis observed in p73KO-mBOs is compensated after p53 loss, indicating that neuronal apoptosis resulting from p73 elimination is p53-dependent. In mBOs, the absence of p53 led to a subtle disorganization, even in the absence of p73 (DKO-mBOs), as well as a slight accelerated neuronal differentiation, in accordance with previous data shown in human p53 knockdown organoids (Marin Navarro et al., 2020). All these results show that p53 family members play important and largely non-redundant roles in mouse brain development.
Mechanistically, the functional interaction between p73 isoforms and p53 regulation has been well studied. Like p53, TAp73 can transactivate target genes that regulate apoptosis and senescence (Rozenberg et al., 2021). On the other hand, DNp73 acts as a dominant negative for TAp73 (and p53) through a direct competition for the promoter or by formation of inactive hetero-oligomeric complexes, showing an antiapoptotic effect (Di et al., 2013). Moreover, both TAp73 and p53 activate the DNp73 promoter creating a negative feedback loop (Rozenberg et al., 2021). Thus, DNp73 is activated after DNA damage in a p53-dependent manner to regulate p53-induced cell cycle arrest (Vossio et al., 2002). However, the pro-survival role of DNp73 can only partially explain the neural phenotypes, both in the various p73 knockout models (Niklison-Chirou et al., 2017) and in our mBOs.
Based on our proliferation, apoptosis and transcriptomic data, we hypothesize that DNp73 is essential to sustain proliferating neural progenitor cells, but it is also a regulator of early neuronal fate acquisition during development. The precise patterning of brain organoids relies heavily on the strength and duration of exposure to diffusible morphogens such as Shh and Wnts, as well as the intricate crosstalk among different signaling pathways such as Notch (Zhao and Haddad, 2024). In our model, lack of DNp73 results in the enhanced activation of the Shh, Wnt, and Notch signaling cascades, along with an early and strong repression of BMP and Activin/Nodal signaling. The functional interplay of p73 isoforms with these signaling pathways, especially DNp73, is not yet completely understood. Several reports have linked p73 function to the Notch signaling pathway since the loss of total p73 leads to the transcriptional dysregulation of NOTCH-1, NOTCH-2, HES-5, JAG2, HEY-2, and DELTEX (Niklison-Chirou et al., 2016), being HEY2 a direct TAp73 transcriptional target (Fujitani et al., 2010). Moreover, TAp73 isoforms, but not DNp73, are capable of directly binding the Notch-1 intracellular domain (N1ICD) and antagonize its transcriptional activity (Hooper et al., 2006). However, additional studies are required to unravel how DNp73 physiologically regulates these factors.
On the other hand, Wnt3 and its receptor Fzd1 are direct p53, TAp63 and TAp73 target genes in ESCs (Wang et al., 2017). Thus, the elevated TAp73 levels detected in DNp73KO cells (Lopez-Ferreras et al., 2021) could be responsible for the enhanced Wnt signaling. However, whether DNp73 acts directly over this signaling cascade remains to be clarified. The interaction between p73 and the BMP pathway is better understood. P73 is a positive modulator of the BMP circuit; DNp73 expression enhances BMP4-induced Smad1/5/8 activation and directly represses the activation of the SMAD6 promoter (Martin-Lopez et al., 2017). Thus, DNp73 depletion could hinder BMP signaling cascade during mBO formation. Moreover, similarly to our observation in mBOs, during ESC differentiation, depletion of p53 and p73 inhibited Smad2/3/Nodal signaling, and the p53/p73-depleted cells adopted a neuroectoderm fate (Wang et al., 2017). This change into an alternative neuroectoderm fate is consistent with Nodal inhibition during development (Vallier et al., 2004). Thus, our data support the idea of DNp73 as a regulator of the early neuroectoderm fate switch.
Overall, we have optimized a model to efficiently generate mouse cortical brain organoids. In fact, we have utilized this in vitro culture system to gain further insight into the specific function of p53 and p73 isoforms in brain morphogenesis, laying the foundation of a new understanding of p53 and p73 in neurodevelopment and validating the reliability of this 3D model as a useful tool for biomedical and neurodevelopmental biology research. Moreover, we propose a dual role of p73 regulating brain morphogenesis, whereby TAp73 governs transcriptional programs essential for the establishment of the neurogenic niche structure, and DNp73 is central for a precise and timely regulation of neural cell fate.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://open.scayle.es/dataset/alonso-olivares-et-al-2024, Alonso-olivares-et-al-2024.
Ethics statement
The animal study was approved by the appropriate institutional committee (Comité de Etica, Universidad de León, 10/06/2011; OEBA-ULE-018-2017; ETICA-ULE-053-2022). All miPSC lines used in this study were previously generated in our lab by reprogramming of mouse embryonic fibroblasts from different genotypes (Martin-Lopez et al., 2017). Breeding of wild-type mice (WT), Trp73 knockout (p73KO), Trp53 knockouts (p53KO) and the double mutants (DKO) was carried out in agreement with European (Council Directive 2010/63/UE) and Spanish regulations (RD 53/2013) on the protection of experimental animals. E14TG2α mESCs were kindly provided by Dr. Jim McWhir (former Researcher at the Roslin Institute, Edinburgh, Scotland, United Kingdom). E14TG2α-TAp73KO and E14TG2α-DNp73KO cells were previously generated in our lab (Lopez-Ferreras et al., 2021). The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
HA-O: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Writing–original draft, Writing–review and editing, Visualization. MMM: Conceptualization, Data curation, Formal Analysis, Supervision, Writing–original draft, Writing–review and editing, Funding acquisition, Project administration, Visualization. AP-C: Writing–review and editing, Methodology. LL-F: Methodology, Writing–review and editing, Investigation. NM-G: Writing–review and editing, Methodology, Investigation. AV-J: Methodology, Writing–review and editing, Investigation. VB: Conceptualization, Writing–review and editing, Funding acquisition, Supervision, Writing–original draft. MCM: Conceptualization, Funding acquisition, Resources, Supervision, Writing–original draft, Writing–review and editing, Formal Analysis, Visualization. RF-A: Conceptualization, Supervision, Writing–original draft, Data curation, Formal Analysis, Writing–review and editing, Visualization.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by Grants PID 2019-105169RB-I00 funded by MCIN/AEI/10.13039/501100011033 and LE022P20 from Junta de Castilla y León (cofinanced with FEDER funds) (to MCM), and by European Research Council (101118729) to VB, who also acknowledges financial support from the AEI through the “Severo Ochoa” Program for Centers of Excellence in R&D (CEX2021-001165-S funded by MCIN/AEI/10.13039/501100011033). Fundación Asociación Española Contra el Cáncer grant (PRYGN21941BORR). RF-A was supported by a María Zambrano fellowship from AEI. HA-O and NM-G are supported by a predoctoral scholarship from the Asociación Española Contra el Cáncer (AECC), and AP-C was supported by a Fundación Tatiana Pérez de Guzmán el Bueno fellowship.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2024.1464932/full#supplementary-material
References
Agostini, M., Niklison-Chirou, M. V., Annicchiarico-Petruzzelli, M. M., Grelli, S., Di Daniele, N., Pestlikis, I., et al. (2018). p73 regulates primary cortical neuron metabolism: a global metabolic profile. Mol. Neurobiol. 55 (4), 3237–3250. doi:10.1007/s12035-017-0517-3
Agostini, M., Tucci, P., Chen, H., Knight, R. A., Bano, D., Nicotera, P., et al. (2010). p73 regulates maintenance of neural stem cell. Biochem. Biophys. Res. Commun. 403 (1), 13–17. doi:10.1016/j.bbrc.2010.10.087
Armesilla-Diaz, A., Bragado, P., Del Valle, I., Cuevas, E., Lazaro, I., Martin, C., et al. (2009). p53 regulates the self-renewal and differentiation of neural precursors. Neuroscience 158 (4), 1378–1389. doi:10.1016/j.neuroscience.2008.10.052
Belmonte-Mateos, C., and Pujades, C. (2021). From cell states to cell fates: how cell proliferation and neuronal differentiation are coordinated during embryonic development. Front. Neurosci. 15, 781160. doi:10.3389/fnins.2021.781160
Besser, D. (2004). Expression of nodal, lefty-a, and lefty-B in undifferentiated human embryonic stem cells requires activation of Smad2/3. J. Biol. Chem. 279 (43), 45076–45084. doi:10.1074/jbc.M404979200
Bifari, F., Dolci, S., Bottani, E., Pino, A., Di Chio, M., Zorzin, S., et al. (2020). Complete neural stem cell (NSC) neuronal differentiation requires a branched chain amino acids-induced persistent metabolic shift towards energy metabolism. Pharmacol. Res. 158, 104863. doi:10.1016/j.phrs.2020.104863
Borrell, V. (2019). Recent advances in understanding neocortical development. F1000Res 8. doi:10.12688/f1000research.20332.1
Cahoy, J. D., Emery, B., Kaushal, A., Foo, L. C., Zamanian, J. L., Christopherson, K. S., et al. (2008). A transcriptome database for astrocytes, neurons, and oligodendrocytes: a new resource for understanding brain development and function. J. Neurosci. 28 (1), 264–278. doi:10.1523/JNEUROSCI.4178-07.2008
Cai, E., Barba, M. G., and Ge, X. (2023). Hedgehog signaling in cortical development. Cells 13 (1), 21. doi:10.3390/cells13010021
Chambers, S. M., Fasano, C. A., Papapetrou, E. P., Tomishima, M., Sadelain, M., and Studer, L. (2009). Highly efficient neural conversion of human ES and iPS cells by dual inhibition of SMAD signaling. Nat. Biotechnol. 27 (3), 275–280. doi:10.1038/nbt.1529
Chen, C., and Shen, M. M. (2004). Two modes by which Lefty proteins inhibit nodal signaling. Curr. Biol. 14 (7), 618–624. doi:10.1016/j.cub.2004.02.042
Chenn, A., and Walsh, C. A. (2002). Regulation of cerebral cortical size by control of cell cycle exit in neural precursors. Science. 297 (5580), 365–369. doi:10.1126/science.1074192
Ciarpella, F., Zamfir, R. G., Campanelli, A., Ren, E., Pedrotti, G., Bottani, E., et al. (2021). Murine cerebral organoids develop network of functional neurons and hippocampal brain region identity. iScience 24 (12), 103438. doi:10.1016/j.isci.2021.103438
Di, C., Yang, L., Zhang, H., Ma, X., Zhang, X., Sun, C., et al. (2013). Mechanisms, function and clinical applications of DNp73. Cell Cycle 12 (12), 1861–1867. doi:10.4161/cc.24967
Eichmuller, O. L., and Knoblich, J. A. (2022). Human cerebral organoids - a new tool for clinical neurology research. Nat. Rev. Neurol. 18 (11), 661–680. doi:10.1038/s41582-022-00723-9
Eiraku, M., and Sasai, Y. (2011). Mouse embryonic stem cell culture for generation of three-dimensional retinal and cortical tissues. Nat. Protoc. 7 (1), 69–79. doi:10.1038/nprot.2011.429
Eiraku, M., Watanabe, K., Matsuo-Takasaki, M., Kawada, M., Yonemura, S., Matsumura, M., et al. (2008). Self-organized formation of polarized cortical tissues from ESCs and its active manipulation by extrinsic signals. Cell Stem Cell 3 (5), 519–532. doi:10.1016/j.stem.2008.09.002
Fernandez-Alonso, R., Martin-Lopez, M., Gonzalez-Cano, L., Garcia, S., Castrillo, F., Diez-Prieto, I., et al. (2015). p73 is required for endothelial cell differentiation, migration and the formation of vascular networks regulating VEGF and TGFβ signaling. Cell Death Differ. 22 (8), 1287–1299. doi:10.1038/cdd.2014.214
Forsberg, K., Wuttke, A., Quadrato, G., Chumakov, P. M., Wizenmann, A., and Di Giovanni, S. (2013). The tumor suppressor p53 fine-tunes reactive oxygen species levels and neurogenesis via PI3 kinase signaling. J. Neurosci. 33 (36), 14318–14330. doi:10.1523/JNEUROSCI.1056-13.2013
Fuertes-Alvarez, S., Maeso-Alonso, L., Villoch-Fernandez, J., Wildung, M., Martin-Lopez, M., Marshall, C., et al. (2018). p73 regulates ependymal planar cell polarity by modulating actin and microtubule cytoskeleton. Cell Death Dis. 9 (12), 1183. doi:10.1038/s41419-018-1205-6
Fujitani, M., Cancino, G. I., Dugani, C. B., Weaver, I. C., Gauthier-Fisher, A., Paquin, A., et al. (2010). TAp73 acts via the bHLH Hey2 to promote long-term maintenance of neural precursors. Curr. Biol. 20 (22), 2058–2065. doi:10.1016/j.cub.2010.10.029
Gaiano, N., and Fishell, G. (2002). The role of notch in promoting glial and neural stem cell fates. Annu. Rev. Neurosci. 25, 471–490. doi:10.1146/annurev.neuro.25.030702.130823
Gaspard, N., Bouschet, T., Hourez, R., Dimidschstein, J., Naeije, G., van den Ameele, J., et al. (2008). An intrinsic mechanism of corticogenesis from embryonic stem cells. Nature 455 (7211), 351–357. doi:10.1038/nature07287
Gonzalez-Cano, L., Fuertes-Alvarez, S., Robledinos-Anton, N., Bizy, A., Villena-Cortes, A., Farinas, I., et al. (2016). p73 is required for ependymal cell maturation and neurogenic SVZ cytoarchitecture. Dev. Neurobiol. 76 (7), 730–747. doi:10.1002/dneu.22356
Gonzalez-Cano, L., Herreros-Villanueva, M., Fernandez-Alonso, R., Ayuso-Sacido, A., Meyer, G., Garcia-Verdugo, J. M., et al. (2010). p73 deficiency results in impaired self renewal and premature neuronal differentiation of mouse neural progenitors independently of p53. Cell Death Dis. 1 (12), e109. doi:10.1038/cddis.2010.87
Gopalakrishnan, J. (2019). The emergence of stem cell-based brain organoids: trends and challenges. Bioessays 41 (8), e1900011. doi:10.1002/bies.201900011
Grespi, F., Amelio, I., Tucci, P., Annicchiarico-Petruzzelli, M., and Melino, G. (2012). Tissue-specific expression of p73 C-terminal isoforms in mice. Cell Cycle 11 (23), 4474–4483. doi:10.4161/cc.22787
Holembowski, L., Kramer, D., Riedel, D., Sordella, R., Nemajerova, A., Dobbelstein, M., et al. (2014). TAp73 is essential for germ cell adhesion and maturation in testis. J. Cell Biol. 204 (7), 1173–1190. doi:10.1083/jcb.201306066
Hooper, C., Tavassoli, M., Chapple, J. P., Uwanogho, D., Goodyear, R., Melino, G., et al. (2006). TAp73 isoforms antagonize Notch signalling in SH-SY5Y neuroblastomas and in primary neurones. J. Neurochem. 99 (3), 989–999. doi:10.1111/j.1471-4159.2006.04142.x
Hooper, M., Hardy, K., Handyside, A., Hunter, S., and Monk, M. (1987). HPRT-deficient (Lesch-Nyhan) mouse embryos derived from germline colonization by cultured cells. Nature 326 (6110), 292–295. doi:10.1038/326292a0
Huang da, W., Sherman, B. T., and Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4 (1), 44–57. doi:10.1038/nprot.2008.211
Hulsen, T., de Vlieg, J., and Alkema, W. (2008). BioVenn - a web application for the comparison and visualization of biological lists using area-proportional Venn diagrams. BMC Genomics 9, 488. doi:10.1186/1471-2164-9-488
Kim, D. K., Cha, Y., Ahn, H. J., Kim, G., and Park, K. S. (2014). Lefty1 and lefty2 control the balance between self-renewal and pluripotent differentiation of mouse embryonic stem cells. Stem Cells Dev. 23 (5), 457–466. doi:10.1089/scd.2013.0220
Kuo, C. T., Mirzadeh, Z., Soriano-Navarro, M., Rasin, M., Wang, D., Shen, J., et al. (2006). Postnatal deletion of Numb/Numblike reveals repair and remodeling capacity in the subventricular neurogenic niche. Cell 127 (6), 1253–1264. doi:10.1016/j.cell.2006.10.041
Lopez-Ferreras, L., Martinez-Garcia, N., Maeso-Alonso, L., Martin-Lopez, M., Diez-Matilla, A., Villoch-Fernandez, J., et al. (2021). Deciphering the nature of Trp73 isoforms in mouse embryonic stem cell models: generation of isoform-specific deficient cell lines using the CRISPR/Cas9 gene editing system. Cancers (Basel) 13 (13), 3182. doi:10.3390/cancers13133182
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15 (12), 550. doi:10.1186/s13059-014-0550-8
Maeso-Alonso, L., Lopez-Ferreras, L., Marques, M. M., and Marin, M. C. (2021). p73 as a tissue architect. Front. Cell Dev. Biol. 9, 716957. doi:10.3389/fcell.2021.716957
Manzari-Tavakoli, A., Babajani, A., Farjoo, M. H., Hajinasrollah, M., Bahrami, S., and Niknejad, H. (2022). The cross-talks among bone morphogenetic protein (BMP) signaling and other prominent pathways involved in neural differentiation. Front. Mol. Neurosci. 15, 827275. doi:10.3389/fnmol.2022.827275
Marin Navarro, A., Pronk, R. J., van der Geest, A. T., Oliynyk, G., Nordgren, A., Arsenian-Henriksson, M., et al. (2020). p53 controls genomic stability and temporal differentiation of human neural stem cells and affects neural organization in human brain organoids. Cell Death Dis. 11 (1), 52. doi:10.1038/s41419-019-2208-7
Marques, M. M., Villoch-Fernandez, J., Maeso-Alonso, L., Fuertes-Alvarez, S., and Marin, M. C. (2019). The Trp73 mutant mice: a ciliopathy model that uncouples ciliogenesis from planar cell polarity. Front. Genet. 10, 154. doi:10.3389/fgene.2019.00154
Marques-Garcia, F., Ferrandiz, N., Fernandez-Alonso, R., Gonzalez-Cano, L., Herreros-Villanueva, M., Rosa-Garrido, M., et al. (2009). p73 plays a role in erythroid differentiation through GATA1 induction. J. Biol. Chem. 284 (32), 21139–21156. doi:10.1074/jbc.M109.026849
Martin-Lopez, M., Maeso-Alonso, L., Fuertes-Alvarez, S., Balboa, D., Rodriguez-Cortez, V., Weltner, J., et al. (2017). p73 is required for appropriate BMP-induced mesenchymal-to-epithelial transition during somatic cell reprogramming. Cell Death Dis. 8 (9), e3034. doi:10.1038/cddis.2017.432
Medina-Bolivar, C., Gonzalez-Arnay, E., Talos, F., Gonzalez-Gomez, M., Moll, U. M., and Meyer, G. (2014). Cortical hypoplasia and ventriculomegaly of p73-deficient mice: developmental and adult analysis. J. Comp. Neurol. 522 (11), 2663–2679. doi:10.1002/cne.23556
Meletis, K., Wirta, V., Hede, S. M., Nister, M., Lundeberg, J., and Frisen, J. (2006). p53 suppresses the self-renewal of adult neural stem cells. Development 133 (2), 363–369. doi:10.1242/dev.02208
Mendrysa, S. M., Ghassemifar, S., and Malek, R. (2011). p53 in the CNS: perspectives on development, stem cells, and cancer. Stem Cells, Cancer. Genes. Cancer 2 (4), 431–442. doi:10.1177/1947601911409736
Metsalu, T., and Vilo, J. (2015). ClustVis: a web tool for visualizing clustering of multivariate data using Principal Component Analysis and heatmap. Nucleic Acids Res. 43 (W1), W566–W570. doi:10.1093/nar/gkv468
Meyer, G., Cabrera, S. A., Perez Garcia, C. G., Martinez Millan, L., Walker, N., and Caput, D. (2004). Developmental roles of p73 in Cajal-Retzius cells and cortical patterning. J. Neurosci. 24 (44), 9878–9887. doi:10.1523/JNEUROSCI.3060-04.2004
Miyata, T., Kawaguchi, A., Saito, K., Kawano, M., Muto, T., and Ogawa, M. (2004). Asymmetric production of surface-dividing and non-surface-dividing cortical progenitor cells. Development 131 (13), 3133–3145. doi:10.1242/dev.01173
Moreau, M. X., Saillour, Y., Elorriaga, V., Bouloudi, B., Delberghe, E., Deutsch Guerrero, T., et al. (2023). Repurposing of the multiciliation gene regulatory network in fate specification of Cajal-Retzius neurons. Dev. Cell 58 (15), 1365–1382.e6. doi:10.1016/j.devcel.2023.05.011
Nasu, M., Takata, N., Danjo, T., Sakaguchi, H., Kadoshima, T., Futaki, S., et al. (2012). Robust formation and maintenance of continuous stratified cortical neuroepithelium by laminin-containing matrix in mouse ES cell culture. PLoS One 7 (12), e53024. doi:10.1371/journal.pone.0053024
Nemajerova, A., and Moll, U. M. (2019). Tissue-specific roles of p73 in development and homeostasis. J. Cell Sci. 132 (19), jcs233338. doi:10.1242/jcs.233338
Niemantsverdriet, M., Nagle, P., Chiu, R. K., Langendijk, J. A., Kampinga, H. H., and Coppes, R. P. (2012). ΔNp73 enhances promoter activity of TGF-β induced genes. PLoS One 7 (12), e50815. doi:10.1371/journal.pone.0050815
Niklison-Chirou, M. V., Erngren, I., Engskog, M., Haglof, J., Picard, D., Remke, M., et al. (2017). TAp73 is a marker of glutamine addiction in medulloblastoma. Genes. Dev. 31 (17), 1738–1753. doi:10.1101/gad.302349.117
Niklison-Chirou, M. V., Killick, R., Knight, R. A., Nicotera, P., Melino, G., and Agostini, M. (2016). How does p73 cause neuronal defects? Mol. Neurobiol. 53 (7), 4509–4520. doi:10.1007/s12035-015-9381-1
Noctor, S. C., Martinez-Cerdeno, V., Ivic, L., and Kriegstein, A. R. (2004). Cortical neurons arise in symmetric and asymmetric division zones and migrate through specific phases. Nat. Neurosci. 7 (2), 136–144. doi:10.1038/nn1172
Noctor, S. C., Martinez-Cerdeno, V., and Kriegstein, A. R. (2008). Distinct behaviors of neural stem and progenitor cells underlie cortical neurogenesis. J. Comp. Neurol. 508 (1), 28–44. doi:10.1002/cne.21669
Nowakowski, T. J., Bhaduri, A., Pollen, A. A., Alvarado, B., Mostajo-Radji, M. A., Di Lullo, E., et al. (2017). Spatiotemporal gene expression trajectories reveal developmental hierarchies of the human cortex. Science 358 (6368), 1318–1323. doi:10.1126/science.aap8809
Osterburg, C., and Dotsch, V. (2022). Structural diversity of p63 and p73 isoforms. Cell Death Differ. 29 (5), 921–937. doi:10.1038/s41418-022-00975-4
Paez-Gonzalez, P., Abdi, K., Luciano, D., Liu, Y., Soriano-Navarro, M., Rawlins, E., et al. (2011). Ank3-dependent SVZ niche assembly is required for the continued production of new neurons. Neuron 71 (1), 61–75. doi:10.1016/j.neuron.2011.05.029
Pasca, S. P., Arlotta, P., Bateup, H. S., Camp, J. G., Cappello, S., Gage, F. H., et al. (2022). A nomenclature consensus for nervous system organoids and assembloids. Nature 609 (7929), 907–910. doi:10.1038/s41586-022-05219-6
Patro, R., Duggal, G., Love, M. I., Irizarry, R. A., and Kingsford, C. (2017). Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14 (4), 417–419. doi:10.1038/nmeth.4197
Pozniak, C. D., Barnabe-Heider, F., Rymar, V. V., Lee, A. F., Sadikot, A. F., and Miller, F. D. (2002). p73 is required for survival and maintenance of CNS neurons. J. Neurosci. 22 (22), 9800–9809. doi:10.1523/JNEUROSCI.22-22-09800.2002
Pozniak, C. D., Radinovic, S., Yang, A., McKeon, F., Kaplan, D. R., and Miller, F. D. (2000). An anti-apoptotic role for the p53 family member, p73, during developmental neuron death. Science 289 (5477), 304–306. doi:10.1126/science.289.5477.304
Roll, L., Lessmann, K., Brustle, O., and Faissner, A. (2022). Cerebral organoids maintain the expression of neural stem cell-associated glycoepitopes and extracellular matrix. Cells 11 (5), 760. doi:10.3390/cells11050760
Rosso, S. B., and Inestrosa, N. C. (2013). WNT signaling in neuronal maturation and synaptogenesis. Front. Cell Neurosci. 7, 103. doi:10.3389/fncel.2013.00103
Rozenberg, J. M., Zvereva, S., Dalina, A., Blatov, I., Zubarev, I., Luppov, D., et al. (2021). The p53 family member p73 in the regulation of cell stress response. Biol. Direct 16 (1), 23. doi:10.1186/s13062-021-00307-5
Rufini, A., Agostini, M., Grespi, F., Tomasini, R., Sayan, B. S., Niklison-Chirou, M. V., et al. (2011). p73 in Cancer. Genes. Cancer 2 (4), 491–502. doi:10.1177/1947601911408890
Saade, M., Gonzalez-Gobartt, E., Escalona, R., Usieto, S., and Marti, E. (2017). Shh-mediated centrosomal recruitment of PKA promotes symmetric proliferative neuroepithelial cell division. Nat. Cell Biol. 19 (5), 493–503. doi:10.1038/ncb3512
Sherman, B. T., Hao, M., Qiu, J., Jiao, X., Baseler, M. W., Lane, H. C., et al. (2022). DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 50 (W1), W216–W221. doi:10.1093/nar/gkac194
Simon, A. (2010). Bioinformatics. Available at: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/.
Takahashi, T., Nowakowski, R. S., and Caviness, V. S. (1995). The cell cycle of the pseudostratified ventricular epithelium of the embryonic murine cerebral wall. J. Neurosci. 15 (9), 6046–6057. doi:10.1523/JNEUROSCI.15-09-06046.1995
Takahashi, T., Nowakowski, R. S., and Caviness, V. S. (1996). The leaving or Q fraction of the murine cerebral proliferative epithelium: a general model of neocortical neuronogenesis. J. Neurosci. 16 (19), 6183–6196. doi:10.1523/JNEUROSCI.16-19-06183.1996
Talos, F., Abraham, A., Vaseva, A. V., Holembowski, L., Tsirka, S. E., Scheel, A., et al. (2010). p73 is an essential regulator of neural stem cell maintenance in embryonal and adult CNS neurogenesis. Cell Death Differ. 17 (12), 1816–1829. doi:10.1038/cdd.2010.131
Taverna, E., Gotz, M., and Huttner, W. B. (2014). The cell biology of neurogenesis: toward an understanding of the development and evolution of the neocortex. Annu. Rev. Cell Dev. Biol. 30, 465–502. doi:10.1146/annurev-cellbio-101011-155801
Tissir, F., Ravni, A., Achouri, Y., Riethmacher, D., Meyer, G., and Goffinet, A. M. (2009). DeltaNp73 regulates neuronal survival in vivo. Proc. Natl. Acad. Sci. U. S. A. 106 (39), 16871–16876. doi:10.1073/pnas.0903191106
Toh, W. H., Logette, E., Corcos, L., and Sabapathy, K. (2008). TAp73beta and DNp73beta activate the expression of the pro-survival caspase-2S. Nucleic Acids Res. 36 (13), 4498–4509. doi:10.1093/nar/gkn414
Tomasini, R., Tsuchihara, K., Wilhelm, M., Fujitani, M., Rufini, A., Cheung, C. C., et al. (2008). TAp73 knockout shows genomic instability with infertility and tumor suppressor functions. Genes. Dev. 22 (19), 2677–2691. doi:10.1101/gad.1695308
Vallier, L., Reynolds, D., and Pedersen, R. A. (2004). Nodal inhibits differentiation of human embryonic stem cells along the neuroectodermal default pathway. Dev. Biol. 275 (2), 403–421. doi:10.1016/j.ydbio.2004.08.031
Van Nostrand, J. L., Bowen, M. E., Vogel, H., Barna, M., and Attardi, L. D. (2017). The p53 family members have distinct roles during mammalian embryonic development. Cell Death Differ. 24 (4), 575–579. doi:10.1038/cdd.2016.128
Villalba, A., Gotz, M., and Borrell, V. (2021). The regulation of cortical neurogenesis. Curr. Top. Dev. Biol. 142, 1–66. doi:10.1016/bs.ctdb.2020.10.003
Vossio, S., Palescandolo, E., Pediconi, N., Moretti, F., Balsano, C., Levrero, M., et al. (2002). DN-p73 is activated after DNA damage in a p53-dependent manner to regulate p53-induced cell cycle arrest. Oncogene 21 (23), 3796–3803. doi:10.1038/sj.onc.1205465
Wang, Q., Zou, Y., Nowotschin, S., Kim, S. Y., Li, Q. V., Soh, C. L., et al. (2017). The p53 family coordinates Wnt and nodal inputs in mesendodermal differentiation of embryonic stem cells. Cell Stem Cell 20 (1), 70–86. doi:10.1016/j.stem.2016.10.002
Watanabe, K., Kamiya, D., Nishiyama, A., Katayama, T., Nozaki, S., Kawasaki, H., et al. (2005). Directed differentiation of telencephalic precursors from embryonic stem cells. Nat. Neurosci. 8 (3), 288–296. doi:10.1038/nn1402
Watanabe, K., Nakamura, M., Okano, H., and Toyama, Y. (2007). Establishment of three-dimensional culture of neural stem/progenitor cells in collagen Type-1 Gel. Restor. Neurol. Neurosci. 25 (2), 109–117.
Wataya, T., Ando, S., Muguruma, K., Ikeda, H., Watanabe, K., Eiraku, M., et al. (2008). Minimization of exogenous signals in ES cell culture induces rostral hypothalamic differentiation. Proc. Natl. Acad. Sci. U. S. A. 105 (33), 11796–11801. doi:10.1073/pnas.0803078105
Wilhelm, M. T., Rufini, A., Wetzel, M. K., Tsuchihara, K., Inoue, S., Tomasini, R., et al. (2010). Isoform-specific p73 knockout mice reveal a novel role for delta Np73 in the DNA damage response pathway. Genes. Dev. 24 (6), 549–560. doi:10.1101/gad.1873910
Xiong, Y., Zhang, Y., Xiong, S., and Williams-Villalobo, A. E. (2020). A glance of p53 functions in brain development, neural stem cells, and brain cancer. Brain Cancer. Biol. (Basel). 9 (9), 285. doi:10.3390/biology9090285
Xu, X., Li, L., Luo, L., Shu, L., Si, X., Chen, Z., et al. (2021). Opportunities and challenges of glioma organoids. Cell Commun. Signal 19 (1), 102. doi:10.1186/s12964-021-00777-0
Yang, A., Walker, N., Bronson, R., Kaghad, M., Oosterwegel, M., Bonnin, J., et al. (2000). p73-deficient mice have neurological, pheromonal and inflammatory defects but lack spontaneous tumours. Nature 404 (6773), 99–103. doi:10.1038/35003607
Yang, S., Golkaram, M., Oh, S., Oh, Y., Cho, Y., Yoe, J., et al. (2024). ETV4 is a mechanical transducer linking cell crowding dynamics to lineage specification. Nat. Cell Biol. 26 (6), 903–916. doi:10.1038/s41556-024-01415-w
Zhang, H., Chen, S., Shang, C., Wu, X., Wang, Y., and Li, G. (2019). Interplay between Lefty and Nodal signaling is essential for the organizer and axial formation in amphioxus embryos. Dev. Biol. 456 (1), 63–73. doi:10.1016/j.ydbio.2019.08.006
Keywords: mouse brain organoids, brain development, brain morphogenesis, p53-family, p73, TAp73, DNp73, neurogenesis
Citation: Alonso-Olivares H, Marques MM, Prieto-Colomina A, López-Ferreras L, Martínez-García N, Vázquez-Jiménez A, Borrell V, Marin MC and Fernandez-Alonso R (2024) Mouse cortical organoids reveal key functions of p73 isoforms: TAp73 governs the establishment of the archetypical ventricular-like zones while DNp73 is central in the regulation of neural cell fate. Front. Cell Dev. Biol. 12:1464932. doi: 10.3389/fcell.2024.1464932
Received: 15 July 2024; Accepted: 04 September 2024;
Published: 23 September 2024.
Edited by:
Finosh Thankam, Western University of Health Sciences, United StatesReviewed by:
Megha Agarwal, Stanford University, United StatesXiaolei Li, University of Pennsylvania, United States
Copyright © 2024 Alonso-Olivares, Marques, Prieto-Colomina, López-Ferreras, Martínez-García, Vázquez-Jiménez, Borrell, Marin and Fernandez-Alonso. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Maria C. Marin, carmen.marin@unileon.es; Rosalia Fernandez-Alonso, rfera@unileon.es