Mathematical modeling supports fate restriction in neurogenic progenitors of the embryonic ventral spinal cord

In the developing neural tube in chicken and mammals, neural stem cells proliferate and differentiate according to a stereotyped spatio-temporal pattern. Several actors have been identified in the control of this process, from tissue-scale morphogens patterning (Shh, BMP) to intrinsic determinants in neural progenitor cells. In a previous study (Bonnet et al. eLife 7, 2018), we have shown that the CDC25B phosphatase promotes the transition from proliferation to differentiation in a cell-cycle independent fashion. In this study, we set up a mathematical model linking progenitor modes of division to the dynamics of progenitors and differentiated populations. Here, we build on this previous model to propose a complete dynamical picture of this process. We start from the standard model in which progenitors are homogeneous and can perform any type of divisions (proliferative division yielding two progenitors, asymmetric neurogenic divisions yielding one progenitor and one neuron, and terminal symmetric divisions yielding two neurons). We constraint this model using published data about mode of divisions and population dynamics of progenitors/neurons at different developmental stages (Saade et al. Cell Reports 4, 2013), and check the effect of CDC25B gain of function in this context. Next, we explore the scenarios in which progenitors population is actually split into two different pools, one of which composed of cells that have lost the capacity to perform proliferative divisions (fate restriction). We show that one such scenario appears relevant and calls for further identification of the alternative role of CDC25B in such a fate restriction.