Abstract
During development, cortical neurons are produced in a temporally regulated sequence from apical progenitors, directly or indirectly, through the production of intermediate basal progenitors. The balance between these major progenitor types is critical for the production of the proper number and types of neurons, and it is thus important to decipher the cellular and molecular cues controlling this equilibrium. Here we address the role of a cell cycle regulator, the CDC25B phosphatase, in this process. We show that, in the developing mouse neocortex of both sex, deleting CDC25B in apical progenitors leads to a transient increase in the production of TBR1+ neurons at the expense of TBR2+ basal progenitors. This phenotype is associated with lengthening of the G2 phase of the cell cycle, the total cell cycle length being unaffected. Using in utero electroporation and cortical slice cultures, we demonstrate that the defect in TBR2+ basal progenitor production requires interaction with CDK1 and is because of the G2 phase lengthening in CDC25B mutants. Together, this study identifies a new role for CDC25B and G2 phase length in direct versus indirect neurogenesis at early stages of cortical development.
SIGNIFICANCE STATEMENT This study is the first analysis of the function of CDC25B, a G2/M regulator, in the developing neocortex. We show that removing CDC25B function leads to a transient increase in neuronal differentiation at early stages, occurring simultaneously with a decrease in basal intermediate progenitors (bIPs). Conversely, a CDC25B gain of function promotes production of bIPs, and this is directly related to CDC25B's ability to regulate CDK1 activity. This imbalance of neuron/progenitor production is linked to a G2 phase lengthening in apical progenitors; and using pharmacological treatments on cortical slice cultures, we show that shortening the G2 phase is sufficient to enhance bIP production. Our results reveal the importance of G2 phase length regulation for neural progenitor fate determination.
Introduction
During mammalian brain development, neural progenitor cells (NPCs) undergo precise sequential transitions that will define the number and type of neurons that will compose the neocortex (Kawaguchi, 2019). It is therefore important that the proliferation and differentiation rates of these progenitors be tightly regulated during neurogenesis, to avoid dramatic cortical defects, such as microcephaly in humans (Barbelanne and Tsang, 2014).
At the onset of corticogenesis, NPCs initially divide symmetrically to expand the stem cell pool, before giving rise to apical dividing progenitors (Govindan and Jabaudon, 2017). Apical radial glial cells (aRGs) first divide symmetrically and then shift to an asymmetric self-renewing division, producing an aRG and either a neuron (direct neurogenesis) or a basal intermediate progenitor (bIP), that divide once to produce two neurons (indirect neurogenesis) (Govindan and Jabaudon, 2017; Kawaguchi, 2019). Another type of basal progenitor, basal radial glia (bRG, also called outer radial glia), has emerged during evolution and is also present in very small proportions in the subventricular zone (SVZ) and intermediate zone of rodents (Namba and Huttner, 2017; Ostrem et al. 2017). All these progenitors produce sequentially the different types of projection neurons that colonize the six layers of the adult neocortex (Kawaguchi, 2019). The temporal transitions from one type of progenitor to another are therefore key parameters in controlling the size and functionality of the adult neocortex (Wilsch-Bräuninger et al., 2016).
Cell cycle kinetics, and in particular cell cycle phase duration, is of major importance for controlling these transitions (Dehay and Kennedy, 2007; Agius et al., 2015). Several studies showed that G1 lengthening is an important positive regulator of neuronal differentiation, G1 phase lengthening being associated with the transition from aRGs to bIPs during corticogenesis (Lange et al., 2009; Pilaz et al., 2009; Arai et al., 2011). In addition, it has been shown that proliferating and neurogenic progenitors display different S phase duration, the latter having a shorter S phase (Arai et al., 2011). Mitosis lengthening is also linked to an increase in neurogenic divisions (Pilaz et al., 2016). However, even if correlations have been made in different organs between the duration of G2 phase and cell fate (Locker et al., 2006; Agathocleous et al., 2007; Gonzales and Liang, 2015; Gonzales et al., 2015), the importance of controlling G2 length during neurogenesis remains poorly documented.
A few years ago, we showed that the CDC25B phosphatase, a master regulator of mitosis entry, favors neurogenic division at the expense of proliferative ones, and that part of this function is independent of its interaction with CDK1 (Benazeraf et al., 2006; Peco et al., 2012; Bonnet et al., 2018). CDC25B is a well-known G2/M regulator, belonging to the dual specificity CDC25 phosphatase family that activates CDK1-cyclinB complexes and thus promotes mitosis entry (Boutros et al., 2007). It is also recruited to the mother centrosome and is involved in the centrosome duplication cycle, as well as in microtubule nucleation (Boutros et al., 2007, 2011, 2013; Boutros and Ducommun, 2008). Our previous work shows that it allows progenitor maturation in the developing spinal cord, but it remains an open question whether this phosphatase can be considered a general player in neuronal maturation and whether regulation of G2 phase duration has a role to play in this process.
In this study, we investigate the role of CDC25B in corticogenesis. Using a conditional CDC25B KO mouse line, we show that CDC25B loss of function in apical progenitors leads to a transient imbalance in neuronal and bIP production, accompanied by a significant lengthening of the G2 phase. We provide evidence that this change in G2 duration in apical progenitors causes the switch in fate, producing neurons instead of bIPs. This study sheds light on an important role of CDC25B and G2 phase modulation for controlling the fate of apical progenitors, thus impacting the switch between direct to indirect neurogenesis.
Materials and Methods
Mice
Experiments were performed in accordance with European Community guidelines regarding care and use of animals, in agreement with the Ministère de l'Enseignement Supérieur et de la Recherche (number C3155511, reference 01024.01) and CNRS recommendations. Cdc25b+/−; Nestin-Cre, and Cdc25bflox/flox lines have been described previously (Tronche et al., 1999; Bonnet et al., 2018). These lines were crossed to obtain conditional mutant embryos Cdc25flox/−; nestin-Cre (cKO) and control embryos Cdc25bflox/–, Cdc25bflox/+ and Cdc25b flox/+; nestin-Cre (Control) of either sex. Tis21-GFP line (Tis21−/tm2(Gfp)Wb) has been described previously (Haubensak et al., 2004). This line was crossed to obtain conditional mutant embryos Cdc25bflox/–; nestin-Cre; Tis21-GFP that express GFP under the control of the Tis21 promotor. Embryonic (E) day 0.5 was assumed to start at midday of the day of the vaginal plug.
ISH and immunohistochemistry
ISH was performed on 60 µm vibratome coronal sections of mouse brain embryos as described previously (Bonnet et al., 2018). For immunohistochemistry, embryos were fixed in 4% PFA at 4°C overnight, washed in PBS, and embedded in paraffin using an STP 120 Spin Tissue Processor. Brains were sliced in 6-14 µm transversal sections using a HM 355 Microtom. Sections were directly mounted on coated histologic slides. Sections were rehydrated in successive baths of histoclear (twice), 100% ethanol (twice), 95% ethanol, 70% ethanol, 50% ethanol, and deionized H20 for 5 min each. Sections were heated in 1% citrate in H2O at ∼100°C for 45 min, permeabilized with PBS-0.5% Tween + 3% BSA + 1% horse serum and incubated with primary antibody overnight at 4°C. The following primary antibodies were used against: P-H3 (Millipore Biotechnology, 1:1000, rabbit), BrdU (G3G4, 1:400, mouse), active caspase 3 (BD Biosciences, 1:100, rabbit), GFP (Abcam, 1:1000, chick), CCND1 (Fisher Scientific, 1:200, rabbit), SATB2 (Abcam, 1:100, mouse), TBR1 (Abcam, 1:100, rabbit), PAX6 (Covance, 1:300, rabbit), PAX6 (MBL; 1:200, rabbit), TBR2 (eBioscience, 1:50, rat), RFP (SICGEN, 1:100, goat), CTIP2 (Abcam, 1:200, rat), and CUX1 (Santa Cruz Biotechnology, 1:100). Sections were then washed in PBS and incubated with secondary antibodies for 2 h at room temperature. We counterstained the sections with DAPI (1:1000, Invitrogen) to identify nuclei. For BrdU detection, we treated sections with 2N HCl and then with 0.1 m Na4B4O7 before incubation with primary antibodies.
Cell proliferation, cell cycle, and birth dating analysis/5-ethynyl-2′-deoxyuridine (EdU)/BrdU experiments
We determined G2 phase length using the percentage of labeled mitoses (PLM) paradigm (Quastler and Sherman, 1959). We injected pregnant mice intraperitoneally with 100 µl of EdU (Click-iT EdU AlexaFluor-647 Imaging Kit, Invitrogen) at 1 mg/ml. Embryos were harvested 120, 180, 240, 300, 360, or 480 min after EdU administration, and we quantified the percentage of P-H3+EdU+ colabeled nuclei with increasing times of exposure to EdU. The progression of this percentage is proportional to G2 phase duration. To determine cell proliferation, we injected pregnant mice intraperitoneally with 1 mg/ml of EdU. Embryos were harvested 60 min later. Total cell cycle length and S phase length were determined as described previously (Martynoga et al., 2005). We injected pregnant mice intraperitoneally with 100 µl of EdU 10 mg/ml, followed by 100 µl of BrdU (Sigma) at 10 mg/ml 90 min later. Embryos were collected 30 min later. We calculated Tc and Ts from counts of EdU+BrdU− and EdU+BrdU+ cells. To determine generation of neurons by pulse-chase experiments (birth dating), we injected pregnant mice intraperitoneally with EdU at E16.5, and we harvested embryos at birth (P0). We deduced G1 phase length by removing the average length of G2 + S + M phases from the average total cell cycle length in control and CDC25BcKO littermates.
Flow cytometry analysis
Flow cytometry analysis was conducted as described by Jungas et al. (2020). Cortical hemispheres were collected in PBS-1% FCS and cells mechanically dissociated through a 40 µm nylon mesh (Clearline). We suspended cells in 1.2 ml ice-cold PBS and fixed them by slowly adding 3 ml cold 100% ethanol while vortexing to obtain a 70% final concentration. The incubation procedure with primary and secondary antibodies was as described for immunofluorescence, with centrifugation at 2200 rpm for 10 min at each step. Finally, cells were resuspended in 500 µl PI/RNase solution for 30 min in the dark. Acquisition was performed using cytoflexS with a minimum of 1000 cells with the Cytexpert software. We excluded debris and doublets using morphometric parameters (FSC/SSC) and pulse area parameters from PI emission, respectively. To set the threshold of PAX6-, TBR2-, CTIP2-, SATB2-, and GFP-positive cells, control samples cells were incubated with secondary antibody and PI only. aRG progenitor DNA content was estimated by analyzing PI emission intensity of PAX6-positive cells only.
Organotypic brain slice culture and drug treatment
Coronal slices of E13.5 control and CDC25BcKO embryos were sectioned at 250 µm on a vibratome (Leica Microsystems). We deposited coronal slices on Millicell Cell Culture insert (0.4 µm, Millipore) in culture medium DMEM-F12, 1× penicilin/streptavidin. Slices remained incubated in 5% CO2 at 38°C for 24 h. For drug treatment, brain slices were incubated in culture medium containing 1 μm PD0166285 (Selleckchem).
To determine cell proliferation and the PLM paradigm, 250 μm EdU was dropped off directly on each slice, for the duration required, depending on the analysis (see above cell proliferation section). Brain slices were fixed 1 h at room temperature in 4% PFA, and sectioned at 60 µm with a vibratome before immunostaining.
DNA constructs and in utero electroporation
In utero electroporation experiments were performed in cortices of E13.5 embryos inside the pregnant mother as described previously (Fawal et al., 2018). Embryos were harvested 28 h later. Gain-of-function experiments were performed using a vector expressing two human Cdc25B forms, hCDC25B3 (Cdc25B WT) or hCDC25B3delCDK (Cdc25BdelCDK), under the control of a cis regulatory element of the mouse Cdc25B called pccRE (Körner et al., 2001; Bonnet et al., 2018). A control vector was generated with the β-Gal gene downstream of the pccRE; 1.5 mg/ml of either plasmid expressing Cdc25b or control plasmid was coinjected with 1 µg/µl pCIG-GFP and Fast Green (Sigma) in the lateral ventricle of embryos neocortex, manually using a beveled and calibrated glass micropipette. For electroporation, five 50 ms pulses of 35 V with a 950 ms interval were delivered across the uterus with two 3 mm electrode paddles positioned on either side of the head. Before surgery, mice were subcutaneously injected with 100 µl of morphine at 2.5 mg/kg. Mice were anesthetized with 4% isoflurane. After surgery, mice were subcutaneously injected with 100 µl of 2 mg/ml metacam.
Imaging and statistical analysis
All images of immunostained sections were acquired using a SP8 Leica confocal microscope. For each experiment, we analyzed at least three independent litters and three different slides per embryo. The quantifications were performed in the lateral half of the intermediary cortex, except for the G2 length estimation, for which the P-H3/EdU count was performed in the entire hemisection. For each experiment, we verified that the cell density was the same between Ctl and cKO. Quantitative data are expressed as mean ± SEM. They are represented as scatter plots with a bar, each point showing the value for one embryo, or as a violin plot, when the number of measures was large enough (see Fig. 5). Table 1 was established from the values reported in Figure 6. Statistical analyses were performed using Graph Pad Prism 9 or the RStudio software. Significance was assessed by performing the Mann–Whitney t test, the Wilcoxon test, or using the mixed effect model followed by an ANOVA test as described previously (Bonnet et al., 2018).
Cell cycle parameters in Ctl and cKO aRGs
Results
Cdc25b is expressed in apical progenitors
To determine whether Cdc25b is expressed in all cortical progenitors or in a subset of them, we performed ISH on brain coronal sections at key stages of corticogenesis (Fig. 1A). Cdc25b is detected in neuroepithelial cells from E10.5. At E11.5, Cdc25b is expressed in the ventricular zone (VZ) of the cortex, consisting of proliferating progenitors. At E14.5 and E16.5, it still appears restricted to the VZ that contains PAX6+ apical progenitors (aRGs) and TBR2+ newborn basal progenitors (bIPs) (Fig. 1A). To determine whether it is also expressed in the SVZ, we performed immunostaining for TBR2 following ISH of Cdc25b transcripts at E14.5. We observed that Cdc25b is not expressed in the SVZ, which contains TBR2+ bIPs, and that the TBR2+ cells in the VZ (newborn bIPs) (Arai et al., 2011) do not overlap with Cdc25b-expressing cells either (Fig. 1B). To confirm this observation, we analyzed published bulk RNA-Seq data from E14.5 samples (Florio et al., 2015) and scRNA-Seq data from E12 to E15 cortices (Telley et al., 2019). This data mining confirmed that Cdc25b is highly expressed in aRGs, barely in bIPs, and is not expressed in bRGs and neurons (Fig. 1C,D). Overall, these analyses indicate that, in the developing cortex, Cdc25b is predominantly expressed in aRGs, progressively extinguished in bIPs, and absent in neurons.
Cdc25B expression pattern during corticogenesis. Cdc25B is expressed in aRG progenitors (PAX6+ in the VZ) and progressively downregulated in basal progenitors (TBR2+). A, Cdc25B ISH on mouse embryo coronal sections at E10.5, E11.5, E14.5, and E16.5. B, PAX6 and TBR2 immunostaining and costaining of Cdc25b (ISH) and TBR2 expression (immunofluorescence) on E14.5 coronal brain sections. C, D, Quantification of CDC25B transcripts in progenitors and neurons using bulk RNA-Seq data from Florio et al. (2015) (C) and single-cell RNA-Seq from Telley et al. (2019) using the gene browser at http://genebrowser.unige.ch/telagirdon/ (D).
Cdc25b deletion transiently leads to increased neuronal differentiation and decreased basal progenitor production
Having determined that Cdc25b expression is restricted to the apical progenitors, we next investigated its function during the early phases of corticogenesis in Cdc25bfl/-; NestinCre (Cdc25bcKO or cKO) embryos, in which Cdc25b expression is specifically extinguished in NPCs from the onset of neurogenesis (Bonnet et al., 2018). In these embryos, deletion is complete in the cortex from E10.5 (Madisen et al., 2010); and at E14.5, Cdc25b expression is abolished in the whole brain (Fig. 2A). We first measured cortical hemisphere size at the postnatal stage (P0) and did not observe any gross size defect in Cdc25bcKO newborns (Fig. 2B,C). We then checked the production of different types of neurons over time. TBR1-, CTIP2-, and SATB2-expressing neurons are sequentially produced between E12.5 and E16.5 (Molyneaux et al., 2007; Vasistha et al., 2015). At E13.5, the majority of neurons express TBR1. In Cdc25bcKO, the number of TBR1+ neurons is transiently increased in the most lateral part of the cortical plate at E13.5 (Fig. 2D,E). This increase is no longer significant at E14.5 (Fig. 2F,G). We then checked CTIP2+ and SATB2+ neurons at E14.5 and found that their production occurs normally at this stage in Cdc25bcKO (Fig. 2H–J). At E15.5 and E16.5, no difference is observed in the number of TBR1+, CTIP2+, or SATB2+ neurons (Fig. 2K–M). These results indicate that, in Cdc25bcKO, neuronal production is transiently increased at E13.5 but has returned to normal at later stages.
Cdc25B loss of function triggers a transient increase of neuronal differentiation and a defect of neuronal specification. A, Cdc25B ISH at E14.5 on control (Ctl) and Cdc25b cKO embryo coronal sections. B, C, Cortical area (white dotted line) measurement in Ctl and cKO embryos at P0. C, Each point represents the cortex area in arbitrary unit (ua) for one embryo. Mann–Whitney t test and mixed model. Ctl, n = 29; cKO, n = 10. D, F, TBR1 immunostaining on coronal sections at E13.5 (D) and E14.5 (F). Nuclei are stained with DAPI (blue nuclei). E, G, Quantification of TBR1 neurons in the cortical plate (red nuclei) at E13.5 and E14.5. Each point is the mean value of 3 sections/embryo. Mixed model. E13.5, n = 9 or 10. *p = 0.019. E14.5, n = 6–8. H–J, Quantification of CTIP2 (I) and SATB2 neurons (J) in the cortical plate at E14.5. Each point is the mean value of 3 sections/embryo. Mixed model. CTIP2, n = 9 or 10; SATB2, n = 16. K, TBR1 (red nuclei) and SATB2 (green nuclei) immunostaining on coronal sections at E16.5; nuclei are stained with DAPI (blue nuclei). L, Quantification of SATB2+ and TBR1+ neurons in Ctl and cKO embryos at E16.5. Each point is the mean value of 3 sections/embryo. Ctl, n = 17; cKO, n = 11. Mixed model. M, Dynamics of neuronal production in Ctl and cKO quantified by FACS of SATB2 and CTIP2 neurons from E14.5 to E16.5. Each point represents the mean number of cells of several cortices: E14.5, Ctl, n = 6; cKO, n = 5; E15.5, Ctl, n = 6; cKO, n = 3; E16.5, Ctl, n = 6; cKO, n = 2. N, CUX1, CTIP2, and SATB2 immunostaining on P0 coronal brain sections of Ctl embryos. O–Q, Quantification of CUX1 (O), CTIP2 (P), and SATB2 (Q) cells in Ctl and cKO embryos at P0. Ctl, n = 9; cKO, n = 7. Mixed model. *p = 0.043. Each point is the mean value of 3 sections/embryo. R, TBR1 (gray nuclei) and SATB2 (green nuclei) immunostaining on coronal sections at P0. S, Total number of neurons in layers II-VI (SATB2+TBR1 neurons) at P0 in 100 μm wide on coronal sections from R. T, Quantification of TBR1+ neurons in layer VI in 100 μm wide on coronal sections from R. Each point represents the mean value of 3 sections/embryo. Ctl, n = 10; cKO, n = 8, Mixed model. CP, Cortical plate; PZ, progenitor zone. Scale bars, 20 µm.
To examine whether this transient increase of neuron production at E13.5 induces a change in neuron specification, we quantified the number of neurons in several layers at P0 using CTIP2 (layer V), SATB2, CUX1 (LII-III), and TBR1 (LVI) (Fig. 2N–T). Only the CUX1+ neuron number was slightly reduced in layers II-III (Fig. 2O). As CUX1+ neurons are mainly produced from TBR2+ basal progenitors (Vasistha et al., 2015), this suggests that, in cKO, the increased production of TBR1 neurons at E13.5 is made at the expense of basal progenitor production.
The production of basal progenitors is impaired in cKO
To assess whether basal progenitor production is impaired in cKOs or whether the increased number of TBR1+ neurons is because of premature differentiation of bIPs, we analyzed the dynamics of the production of the different progenitor types in Cdc25bcKO embryos. We examined aRG distribution (PAX6+ cells in the VZ and SVZ), bIP cell production (TBR2+ newborn cells in the VZ), and bIP accumulation (TBR2+ cells that have migrated in the SVZ), in the lateral part of the hemicortex between E12.5 and E16.5 (Fig. 3). The number of aRGs is not significantly affected at any of the four stages analyzed in Cdc25bcKO embryos (Fig. 3B,E,H,K), and they are correctly located in the VZ (Fig. 3A,D,G,J). bIP production starts to be impaired from E12.5, being significantly reduced by E13.5 (Fig. 3G,I), and this decrease leads to a deficit of bIP in the SVZ at E16.5 (Fig. 3J,L). To consolidate these data, we also quantified the total number of PAX6+ and TBR2+ cells in control (Ctl) and Cdc25bcKO cortices using cytometry (Jungas et al., 2020). While the total number of PAX6+ cells is not modified, the number of TBR2+ cells is significantly reduced in Cdc25bcKO at E14.5, confirming that bIP but not aRG production is affected in Cdc25bcKO (Fig. 3M). We verified that this reduction was not because of apoptosis, by quantifying activated Caspase3+ cells in the progenitor zone at E13.5 and E16.5, and we observed no difference between Ctl and Cdc25bcKO cortices (Fig. 4A,B). In addition, we quantified proliferation of bIPs (TBR2+ cells in VZ and SVZ) at E14.5 and found that proliferative and mitotic indexes were normal, ruling out that the reduced bIP number was because of a proliferation defect in these cells (Fig. 4C,D). Together, these results indicate that the reduced number of basal progenitors in Cdc25bcKO is most likely because of a reduction in their production from aRGs. As the number of aRGs is not modified, this suggests that, at E13.5, in the absence of CDC25B, the progeny of aRGs is changed from giving one aRG + one bIP to one aRG + one neuron. If this scenario is true, the balance between proliferative versus neurogenic division should not change in Cdc25bcKO, as all three types of division mentioned above are neurogenic divisions. To assess this, we generated Cdc25bcKO mice carrying a Tis21-GFP allele to distinguish aRGs committed to a neuronal or bIP fate (TIS21+) from self-expanding aRGs (TIS21–) (Haubensak et al., 2004; Attardo et al., 2008). We did not detect any change in the number of TIS21-GFP cells among aRGs at E13.5 in Cdc25bcKO compared with the Ctl (Fig. 4F,G), confirming that CDC25B does not change the mode of division of aRGs from proliferative to neurogenic but rather influences the fate of neurogenic divisions.
Cdc25b loss of function impairs basal progenitor production. A, D, G, J, PAX6 and TBR2 immunostaining on E12.5 (A), E13.5 (D), E14.5 (G), and E16.5 (J) coronal brain sections, in Ctl and cKO embryos. B, C, E, F, H, I, K, L, Quantification of PAX6+ cells in the VZ (aRGs) and TBR2+ cells in the VZ and SVZ (bIPs) at E12.5 (B,C), n = 5; E13.5 (E,F), n = 7 or 8; E14.5 (H,I), n = 9-14, E16.5 (K,L), n = 11-18 in Ctl and cKO embryos. Each point is the mean value of 3 sections/embryo. Mann–Whitney t test and mixed model. I, *p = 0.05. L. *p = 0.025, mixed model. M, Flow cytometry analysis of the proportion of PAX6 + (aRGs + bRGs) and TBR2+ (bIPs) singulates cells for Ctl and cKO E14.5 cortices. Each point represents the percentage of cells for one embryo. Ctl, n = 13; cKO, n = 12. Mann–Whitney. *p = 0.03. Scale bars, 20 µm.
The reduction of bIP cells in cKOs is not because of increased cell death, reduced proliferation rate, or a change of the mode of division of aRGs. A, Immunostaining of CAPS3+ cells on E16.5 coronal sections. B, Quantification of activated CASP3+ cells on frontal brain sections of Ctl and cKO embryos at E13.5 and E16.5. Each point is the mean value of counting in 3 sections/embryo. E13.5: Ctl, n = 7; cKO, n = 4; E16.5: Ctl, n = 9; cKO, n = 7. Mann–Whitney. C, Tbr2 and EdU labeling on coronal section after 1H EdU pulse just before harvesting embryos at E14.5 D, Quantification of EdU+ cells among total Tbr2+ cells (bIPs) in Ctl and cKO embryos at E14.5. Ctl, n = 12; cKO, n = 8. Mann–Whitney. E, Quantification of P-H3+TBR2+ cells (nonapical mitotic bIPs) in E14.5 Ctl and cKO embryos. Ctl, n = 11; cKO, n = 9. Mixed model. F, PAX6 Immunostaining on brain coronal sections at E13.5 in TIS21-GFP+ Ctl or cKO embryos. G, Quantification of PAX6+ TIS21-GFP+ cells (differentiating aRGs) among PAX6+ cells. Each point is the mean value of 3 sections/embryo. Mixed model. Scale bars, 20 µm.
Overall, these analyses of neuron and progenitor production at several developmental stages indicate a transient increase in neuron production, accompanied by a transient decrease in bIP production in Cdc25bcKO embryos, with no change in aRG numbers. This suggests that in Cdc25bcKO, aRGs transiently produce neurons instead of bIPS, and that CDC25B controls the tempo of switch between direct versus indirect neurogenesis.
The control of basal progenitor production by CDC25B requires CDK interaction
We previously showed that CDC25B can act independently of CDK interaction and G2 phase regulation in the developing spinal cord (Bonnet et al., 2018). To test whether its role in cortical neurogenesis is dependent or independent of CDK interaction, we compared the effect of CDC25B gain of function with that of a mutated form of CDC25B, CDC25BdelCDK, which is no longer able to interact with CDK1 to regulate G2 phase length (Bonnet et al., 2018). We used in utero electroporation at E13.5 and analyzed the fate of electroporated cells 28 h later (Fig. 5A). Compared with a GFP+ control vector, overexpressing CDC25B increases the number of TBR2+ cells without changing the number of aRGs (TBR2– cells in the VZ), indicating that CDC25B is sufficient to promote bIP production (Fig. 5B–D). Unlike WT CDC25B, ectopic expression of CDC25BdelCDK does not alter the number of TBR2+ cells (Fig. 5B,C). These results clearly show that CDC25B is sufficient to enhance bIP production and that this effect requires an interaction with CDK, suggesting that this control is exerted via a modification of the duration of the G2 phase.
Cdc25b regulation of bIP production is CDK1-dependent. A, Schematic representation and timeline of in utero electroporation experiments. B, GFP and TBR2 immunostaining on E14.5 brain coronal sections after in utero electroporation of a pmp-LacZ (CTL), pmp-CDC25BWT (CDC25BWT), or pmp-CDC25BdelCDK (CDC25BΔCDK) expression vectors at E13.5. C, Quantification of GFP+ cells, which are TBR2+ (bIPs) per section of CTL, CDC25BWT, or CDC25BΔCDK electroporated embryos. **p = 0.007; ***p < 0.0004; Kruskal–Wallis. D, Quantification of GFP+ cells that are TBR2– in the VZ (aRGs) after electroporation Kruskal–Wallis. *p = 0.003. C, D, Each point represents the mean value for 3 z sections/embryo. Ctl, n = 33; CDC25BWT, n = 29; CDC25BΔCDK, n = 39. IZ, Intermediate zone. Scale bars, 20 µm.
Cdc25b loss of function severely increases the length of the G2 phase without affecting the total duration of the cell cycle
Knowing the role of CDC25B in the control of G2/M cell cycle phase transition, and as bIP production involves CDC25B/CDK interaction, the defect in bIP production could be because of a defect in cell cycle parameters in aRGs of Cdc25bcKO embryos. We first measured the duration of the G2 phase using the PLM paradigm (Quastler and Sherman, 1959). Embryos were harvested 2, 3, 4 h (Fig. 6A) or 6 h (not shown) after EdU (5-ethynyl-2′-deoxyuridine) administration, and we quantified the percentage of P-H3+(Phospho-histone 3) EdU+/P-H3+ with increasing exposure times to EdU (Fig. 6A,B). We found that this percentage was consistently lower in aRGs from Cdc25bcKO embryos.
Cdc25b removal significantly lengthens the G2 phase of aRG progenitors. A, P-H3 labeling after 2, 3, and 4 h of EdU pulse in Ctl and cKO embryos at E14.5. B, Quantification of the proportion of apical (aRG) EdU+ P-H3+/P-H3+ cells (A, arrowheads). TG2 corresponds to the time when 50% of P-H3+ cells are EdU+ (indicated by the dotted lines). Mann–Whitney t test. n = 2-10. *p < 0.05. **p < 0.01. ***p < 0.0001. C, EdU+ cell quantification after 2 h EdU pulse, in 30 µm high bin from apical (bin 1) to basal side (bin 5), in Ctl and cKO embryos at E14.5. Mann–Whitney t test. n = 6-13. ***p < 0.0001. D, Quantification of apical P-H3+ cells in Ctl and cKO E14.5 cortices. Mann–Whitney and mixed model. n = 9-12. *p = 0.023. E, Quantification of EdU+ cells among VZ TBR2– (aRGs) cells after 1 h EdU pulse in Ctl and cKO E14.5 embryos. Mixed model. n = 12. F, EdU/BrdU double-labeling experiments to determine total cell cycle length (TC) and S phase duration (TS) of VZ PAX6+ cells (aRG), which were calculated from counts of EdU+BrdU− and EdU+BrdU+ cells as described previously (Martynoga et al., 2005). Bar plots represent the average length of measured TC and TS in aRGs of Ctl and cKO E14.5 embryos. Mann–Whitney t test and mixed model. n = 9-17. G, CCND1/TBR2 double immunostaining to determine the proportion of progenitor in G1 phase. Quantification of CCND1+ cells among VZ TBR2– (aRGs) cells in Ctl and cKO E14.5 cortices. Mixed model. n = 9. *p = 0.023. In bar plots, each point is the mean value of 3 sections/embryo. Scale bars, 20 µm.
Because the progression of this percentage is proportional to the duration of the G2 phase, we were able to extract an average length of the G2 phase. Loss of function of Cdc25b in aRGs led to a significant increase in the duration of the G2 phase at E14.5, from 2 h 12 in controls to 4 h 05 in Cdc25bcKO already detected at E12.5 (Fig. 6B; Table 1). In addition, we observed a strong delay in apical migration of EdU-labeled cycling progenitors (Fig. 6C), resulting in a reduction in the number of apical mitosis (P-H3+ cells; Fig. 6D), coherent with the G2 phase lengthening. We then determined whether this lengthening of G2 phase duration in progenitors affects their proliferation rate. We collected embryos after 1 h of EdU administration at E14.5 and quantified the EdU+ cells among the aRGs (VZ, TBR2– cells; Fig. 6E). This percentage reflects the proportion of aRGs in S phase, and we did not observe any difference between Ctl and Cdc25bcKO embryos (Fig. 6E). As the rate of proliferating cells was unaffected, we next investigated whether the lengthening of the G2 phase caused a general slowdown of the total cell cycle length in aRGs. We performed dual EdU and BrdU labeling to estimate total cell cycle length (Tc) and S phase (TS) durations, as described by Martynoga et al. (2005). As shown in Figure 6F, neither the total cell cycle length nor the TS phase is significantly modified in aRGs of Cdc25bcKO embryos. From this result, we conclude that a shortening of the G1 phase duration may be associated to the significant lengthening of the G2 phase. To confirm this hypothesis, we directly assessed the length of the G1 phase by counting the number of cells expressing CCND1 (specifically expressed during the G1 phase) among the aRGs. We observed a reduction in the number of CCND1+ labeled nuclei in Cdc25bcKO aRGs (Fig. 6G), indicating that the G1 phase is shorter in these cells.
In conclusion, the removal of CDC25B in aRGs leads to a substantial lengthening of the G2 phase (85% length increase; Table 1), as well as a shortening of the G1, resulting in a normal cell cycle length and proliferation rate. If these length changes are considered in proportion to their average length, then the most considerable change we observed is in the duration of the G2 phase (Table 1).
G2 phase shortening drives an increase in bIP production
As the duration of the G2 phase is almost doubled in Cdc25bcKO, we next investigated whether this lengthening could play a role in controlling the rate of bIP production. To this end, we turned to cortical slice cultures and applied pharmacological treatments to modulate the duration of the G2 phase. First, we treated cortical slices of E13.5 WT embryos with PD0166285, a Wee1/Myt1 inhibitor in conditions that do not induce toxicity (Fig. 7A) (Baffet et al., 2015). As Wee1 and Myt1 prevent mitosis entry through negatively regulating CDK1 activity (Potapova et al., 2009), their inhibition leads to an increase in CDK1 activity and thus to a shortening of the G2 phase. We first determined that PD0166285 treatment at 1 µm has an effect on G2 length in our cultures. For that, we quantified the number of EdU+/P-H3+ cells among the P-H3+ cells, following 2 h EdU incorporation just before culture arrest. The percentage of EdU+ cells was identical between untreated and treated controls (∼22% in both conditions, not shown). This indicates that the rate of proliferation is not significantly altered by pharmacological treatment and that the cortical slices are still in good physiological state at the end of the culture. We confirmed that this treatment effectively reduces the time spent in G2 phase by considering the percentage of EdU+P-H3+/P-H3+ cells in untreated and treated samples (increased 1.5× in treated cortices; Fig. 7B,C). We analyzed the impact of this treatment on the duration of the G1 and S phases and on the total duration of the cell cycle. We found that S phase is significantly reduced by PD0166285 treatment, whereas G1 phase duration and total cell cycle length are not significantly modified (Fig. 7D,E). We then counted the number of TBR2+ cells in treated and untreated slices and observed that pharmacological shortening of the G2 phase leads to a doubling of the number of bIPs in the VZ (Fig. 7F–H), leading to a global increase of bIP cells in the progenitor zone. We then investigated whether this treatment would be sufficient to correct the defect of bIP production in Cdc25bcKO cortices. We therefore applied the same experimental protocol as above to Cdc25bcKO and control littermates. As observed on WT embryos, treatment with 1 μm PD0166285 leads to a marked increase in EdU+P-H3+/P-H3+ cells in both control and Cdc25bcKO-treated slices, indicating a shortening of the G2 phase, with the EdU+P-H3+/P-H3+ ratio in treated Cdc25bcKO slices now at the level of untreated control littermates (Fig. 7I). We then checked the number of TBR2+ cells in the VZ and observed a significant increase of these cells in treated versus untreated Cdc25bcKO slices (1.9-fold increase) and treated versus untreated control littermates (1.5-fold increase) (Fig. 7J,K). This indicates that shortening the G2 length in Cdc25bcKO is sufficient to restore efficient bIP production.
G2 phase shortening triggers bIP production. A, Schematic representation of embryo brain slice culture experiments. B, P-H3 labeling after 2 h of EdU pulse in Ctl and cKO embryos at the end of the culture. C, Percentage of EDU+H3P+/H3P+ cells following 2 h EdU pulse, in untreated and PD0166285-treated WT brain slices. Wilcoxon test. n = 14. ***p < 0.0001. D, Quantification of CCND1+ TBR2– cells (aRG in G1 phase) in WT and WT+PD0166285 brain slices. Wilcoxon test. n = 11-16. E, Total cell cycle length (TC) and S phase duration (TS) of VZ TBR2– cells (aRGs) in WT and WT+PD0166285 brain slices. Wilcoxon test. n = 13. **p < 0.01. F, TBR2 immunostaining on WT and WT+PD0166285 brain slices G, Quantification of bIPs (TBR2+ cells) and aRGs (TBR2– cells) in the VZ of WT and WT+PD0166285 brain slices. Wilcoxon test. n = 17. **p < 0.01. H, Quantification of the number of TBR2+ cells in the SVZ and VZ+SVZ in the same sections as in Figure 6F, G (WT and WT+PD0166285 brain slices). Wilcoxon test. *p = 0.02. **p = 0.007. I, Percentage of EDU+H3P+/H3P+ cells following a 2 h EdU pulse in untreated or PD0166285-treated, Ctl or cKO brain slices. Wilcoxon test. n = 9. **p < 0.01. J, TBR2 immunostaining on untreated or PD0166285-treated Ctl or cKO brain slices. K, Quantification of TBR2+ cells (bIPs) in the VZ of untreated or PD0166285-treated, Ctl or cKO brain slices. Wilcoxon test. n = 9. *p < 0.05. **p < 0.01. Each point represents the mean value of 3 sections/embryo. Scale bars, 20 µm.
Overall, these results indicate that the duration of the G2 phase is an important parameter when controlling bIP production during early phases of cortical development.
Discussion
In this study, we investigate the function of CDC25B during corticogenesis. We show that removing CDC25B function leads to a transient increase in neuron numbers at early stages, accompanied by a decrease in the number of intermediate basal progenitors. This imbalance of neuron/progenitor number is because of a modification of the G2 phase length in apical progenitors. We propose a model in which CDC25B expression in aRGs controls the switch from direct to indirect neurogenesis through modulation of G2 phase length.
In the developing neocortex, neurons can be produced either directly from asymmetric division of apical progenitors or indirectly from symmetric division of intermediate basal progenitors (Govindan and Jabaudon, 2017). Controlling the balance between these two modes of production is important for the integrity and size of the cortex, but it also has an impact on the type of neurons that are produced (Haubensak et al., 2004; Attardo et al., 2008; Govindan and Jabaudon, 2017). In Cdc25bcKO, we observe a specific increase in TBR1+ neurons at E13.5, concomitant with a decrease in bIP number. Conversely, a CDC25B gain of function at E13.5 leads to an increase in bIPs in the VZ. We can exclude that the increase in neuron production is because of an increase in aRG production, as the number of aRG is not altered in Cdc25bcKO. Furthermore, if this scenario were true, one would expect that there would also be more bIPs produced. This alteration in the number of bIPs and neurons is also not because of an overproduction of neurons from bIPs. Indeed, the increase in neuron production in Cdc25bcKO embryos precedes the decrease in bIPs of the SVZ; and if bIPs differentiated into neurons prematurely, we should detect a decrease in bIP cells in the SVZ as early as E13.5, which is not the case. Since the number of aRGs is normal and the ratio of proliferative to neurogenic division is unchanged in Cdc25bcKO, our results indicate a change in the progeny of asymmetrically dividing aRGs, producing neurons instead of bIPs (Fig. 8). This phenotype is transient, as at E16.5 we no longer observe a defect in bIP number in the VZ, and the total number of neurons is unchanged in Cdc25bcKO. This could indicate either a compensation mechanism (Freret-Hodara et al., 2017) or, alternatively, that CDC25B is only required during early corticogenesis stages, at the time of the direct to indirect neurogenesis switch. Hence, we propose that, in Cdc25bcKO, the switch from direct to indirect neurogenesis is delayed, leading to a bIP/neuron production imbalance that is progressively compensated for (Fig. 8). This delay and the resulting decrease in the number of TBR2+ progenitors lead to a slight reduction of the CUX1+ upper layer neuron number at birth. This is consistent with the fact that these neurons are mainly produced from TBR2+ bIPs (Vasistha et al., 2015).
Model for CDC25B function during cortical neurogenesis. A, In CDC25BcKO, apical progenitor asymmetric divisions produce a neuron instead of a basal progenitor. B, In a normal situation, deep layer neurons are produced mainly between E12.5 and E13.5, directly from APs (direct neurogenesis). At E13.5, the production of bIPs begins in the VZ, which later divides in the SVZ to give rise to neurons (indirect neurogenesis). Indirect neurogenesis gradually becomes the main mechanism of neuron production. C, In CDC25BcKO, the ratio between direct and indirect neurogenesis is disturbed. Direct neurogenesis is maintained for a longer period of time, leading to an early production of deep layer neurons at the expense of bIP production. This imbalance is progressively compensated to reach normal neuronal production around E16.5. Because of this delay in the direct to indirect neurogenesis switch, upper layer neuron number is slightly reduced at P0.
In a previous study, we have shown that, in the developing spinal cord, CDC25B promotes progenitor maturation as it increases the number of progenitors performing neurogenic divisions at the expense of those that self-renew (Bonnet et al., 2018). Other studies have also made indirect links between defects in the balance between proliferative versus neurogenic division during cortical neurogenesis and CDC25B misregulation (Gruber et al., 2011; Wu et al., 2014). We did not see any change in the proliferative versus neurogenic mode of division in CDC25B mutant cortices. Our previous observations made in the chicken neural tube indicate that part of CDC25B function during neurogenesis is independent of its function on the cell cycle (Bonnet et al., 2018). In the neocortex, the phenotype we describe requires interaction with CDK1; but in both tissues, CDC25B promotes progenitor maturation.
Another important result is that this neurogenesis defect is linked to a significant modification of the G2 phase length. First, in Cdc25bcKO, the most important change in cell cycle parameters in aRGs is a substantial lengthening of the G2 length. Second, the CDC25BdelCDK gain of function that cannot regulate G2 phase length (Bonnet et al., 2018) does not phenocopy the WT form of CDC25B that increases bIP production, indicating that the defect in bIP production is directly related to the ability of CDC25B to regulate CDK1 activity. Third, we provide evidence that shortening the G2 phase in WT embryos using a pharmacological treatment is sufficient to stimulate bIP production. Fourth, this treatment is sufficient to restore an efficient bIP production in Cdc25bcKO. In these experiments, either the G1 phase length or the S phase length was slightly reduced. However, the changes in these phases are well below the levels that have been described as leading to alterations in aRG progenitor mode of division and fate (Lange et al., 2009; Pilaz et al., 2009; Arai et al., 2011). Furthermore, the only parameter that varies dramatically and significantly in both Cdc25bcKO and PD0166285 treatment is the length of the G2 phase. This argues for a role of G2 phase duration, rather than a role of duration in these experimental conditions. In Cdc25bcKO, the G2 lengthening is associated with a perturbation of the apical nuclear migration, as BrdU+ aRGs nuclei take more time to reach the apical surface to perform mitosis. This perturbation is likely because apical migration and G2 phase processing are linked, and cell cycle progression is a prerequisite to normal apical migration (Ueno et al., 2006; Baffet et al., 2015).
What is the relevance of this finding concerning the link between G2 phase length and bIP production in a normal context? Most studies in which GAP phase length has been measured indicate that, while G1 phase length increases during tissue maturation, G2 phase length appears to be more stable, ∼2 h in aRGs and bIPs, and does not change significantly during embryonic cortex maturation (Pilaz et al., 2009; Arai et al., 2011; Hasenpusch-Theil et al., 2018). However, these studies were performed at the level of cell population, which could have masked some heterogeneity in the duration of the G2 phase between different progenitors. Live imaging analysis performed in the developing spinal cord indicates that the length of G2 phase is more heterogeneous than previously described, ranging from 40 to 175 min (Molina et al., 2022). In the developing neocortex, data provided using real-time analysis of proliferating cell nuclear antigen-expressing apical progenitors show that the proportion of aRGs in G2 phase varies between E12 and E16 (Fousse et al., 2019). A recent study from Hagey et al. (2020) also emphasizes that cortical progenitors are not equivalent regarding G2/M phase length. Indeed, by performing single-cell RNA-Seq at different stages (from E9.5 to E18.5), they identify two temporally distinct cortical differentiation trajectories. At E11.5, progenitors can be separated into two groups: one group of lineage-committed progenitors, giving rise to deep layer neurons; and one group of uncommitted multipotent aRGs that will give rise to later born lineages. The group of lineage-committed progenitors spends more time in G1 and expresses G1-associated cyclin, whereas the uncommitted ones are more prone to be found in G2/M and express high levels of M phase associated genes, including CDK1. Interestingly, Cdc25B is found among the overrepresented genes in the upper layer path, compared with the deep layer 1 (Hagey et al., 2020, their supplementary data 5) (Cdc25B in upper-layer trajectory gene sets, p-adjust = 4.04E-15). Moreover, data mining (Telley et al., 2019; Hagey et al., 2020) indicates that Cdc25B transcripts have the same spatial and temporal profile than genes specifically expressed in the uncommitted E11.5 aRGs (not shown). Hence, CDC25B activity, by accelerating mitosis entry, could contribute to the maintenance of the uncommitted progenitor pool at early stages, to ensure correct production of upper layer neurons later on. Different G2 length duration, by changing the exposure time to instructive signaling pathways such as the Slit-Robo and/or the Notch signaling pathways (Murciano et al., 2002; Cárdenas et al., 2018; Fousse et al., 2019; Hagey et al., 2020), could lead to different fate decisions after mitosis.
In conclusion, our study sheds light on a new role for CDC25B and G2 phase length in controlling the fate of apical progenitors during early steps of corticogenesis, thereby controlling the switch from direct to indirect neuronal differentiation. The apparition of indirect neurogenesis is considered a milestone of mammalian cortical evolution leading to cortical expansion (Borrell, 2019). Expression of CDC25B in aRGs, by favoring bIP production, would be a way of supporting this evolutive trait. Future investigation will be necessary to understand the molecular mechanisms that support this function.
Footnotes
Work in F.P.'s laboratory was supported by Center National de la Recherche Scientifique, Université P. Sabatier, and Agence Nationale de la Recherche ANR-19-CE16-0006-01. M.R. was supported by Ministère de L'Enseignement Supérieur et de la Recherche. We thank Valerie Lobjois for fruitful discussions; and Myriam Roussigné, Elise Cau, and Mohamad Fawal for critical reading of the manuscript.
The authors declare no competing financial interests.
- Correspondence should be addressed to Sophie Bel-Vialar at sophie.vialar{at}univ-tlse3.fr