|Mathematical modeling supports substantial mouse neural progenitor cell death.|
|Jump to Full Text|
|PMID: 19602274 Owner: NLM Status: MEDLINE|
|BACKGROUND: Existing quantitative models of mouse cerebral cortical development are not fully constrained by experimental data. RESULTS: Here, we use simple difference equations to model neural progenitor cell fate decisions, incorporating intermediate progenitor cells and initially low rates of neural progenitor cell death. Also, we conduct a sensitivity analysis to investigate possible uncertainty in the fraction of cells that divide, differentiate, and die at each cell cycle. CONCLUSION: We demonstrate that uniformly low-level neural progenitor cell death, as concluded in previous models, is incompatible with normal mouse cortical development. Levels of neural progenitor cell death up to and exceeding 50% are compatible with normal cortical development and may operate to prevent forebrain overgrowth as observed following cell death attenuation, as occurs in caspase 3-null mutant mice.|
|Michael J McConnell; Hugh R MacMillan; Jerold Chun|
Related Documents :
|17640634 - Temporal and spatial differences in expression of trkb isoforms in rat retina during co...
15204064 - Survival and migration of transplanted neural stem cell-derived dopamine cells in the b...
8674864 - Neural and hormonal regulation of growth of corpora allata in the cockroach, diploptera...
3307484 - Neuroglial and neuroblastoma cell lines are capable of metabolizing ethanol via an alco...
1671664 - Functional hepatocyte heterogeneity in glutamate, aspartate and alpha-ketoglutarate upt...
8239494 - Activation of a calcium magnesium independent endonuclease in human leukemic cell apopt...
|Type: Journal Article; Research Support, N.I.H., Extramural; Research Support, U.S. Gov't, Non-P.H.S. Date: 2009-07-14|
|Title: Neural development Volume: 4 ISSN: 1749-8104 ISO Abbreviation: - Publication Date: 2009|
|Created Date: 2009-08-21 Completed Date: 2009-10-20 Revised Date: -|
Medline Journal Info:
|Nlm Unique ID: 101286574 Medline TA: Neural Dev Country: England|
|Languages: eng Pagination: 28 Citation Subset: IM|
|Crick-Jacobs Center for Theoretical and Computational Biology, Salk Institute for Biological Studies, La Jolla, CA 92037, USA. email@example.com|
|APA/MLA Format Download EndNote Download BibTex|
Cell Death / physiology
Cell Differentiation / physiology
Neurons / physiology*
Stem Cells / physiology*
|MH51699/MH/NIMH NIH HHS|
Journal ID (nlm-ta): Neural Dev
Publisher: BioMed Central
Copyright © 2009 McConnell et al; licensee BioMed Central Ltd.
open-access: This is an Open Access article distributed under the terms of the Creative Commons Attribution License (), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Received Day: 18 Month: 3 Year: 2009
Accepted Day: 14 Month: 7 Year: 2009
collection publication date: Year: 2009
Electronic publication date: Day: 14 Month: 7 Year: 2009
Volume: 4First Page: 28 Last Page: 28
Publisher Id: 1749-8104-4-28
PubMed Id: 19602274
|Mathematical modeling supports substantial mouse neural progenitor cell death|
|Michael J McConnell12||Email: firstname.lastname@example.org|
|Hugh R MacMillan3||Email: email@example.com|
|Jerold Chun2||Email: firstname.lastname@example.org|
1Crick-Jacobs Center for Theoretical and Computational Biology, Salk Institute for Biological Studies, La Jolla, CA 92037, USA
2Helen L Dorris Neuropsychiatric Disorder Institute and Department of Molecular Biology, The Scripps Research Institute, La Jolla, CA 92037, USA
3Department of Mathematical Sciences, Clemson University, Clemson, SC 29634, USA
Essentially every excitatory neuron in the cerebral cortex is born from a heterogeneous pool of mitotic cells (referred to collectively as neural progenitor/precursor cells (NPCs)) in the embryonic ventricular zone (VZ) [1-4]. During a 'neurogenic interval' in mouse – commencing at embryonic day (E)10 in the rostro-medial cortex and concluding at E18 in the caudo-lateral cortex [5-7] – the founding NPC population expands through proliferative divisions until it is exhausted by terminal differentiation and programmed cell death (PCD). NPC proliferation must be balanced with the operation of PCD to produce a sufficient, but not supernumerary, neuronal population. Understanding of the cellular and genetic mechanisms controlling the size of the cerebral cortex, among the most notable distinctions of the brain's evolution [8,9], could benefit from an accurate quantitative model of the fate decisions made by NPCs.
Initial models of mouse NPC fate decisions are insufficiently constrained because it was assumed that cell divisions outside of the VZ do not contribute cortical neurons ; however, it is now clear that non-VZ mitoses contribute significant numbers of cortical neurons [11-13]. Earlier models also lacked direct measurement of an important parameter, the founding NPC population size, instead relying on two related but less direct measurements : changes in the size of the VZ (from which the number of NPCs is extrapolated) and the fraction of cells emigrating from the VZ (from which VZ neuronal output is calculated). Direct measurement of the size of the founding NPC pool [15,16] provides a more accurate description than extrapolation and addresses a limitation of early models. These data are employed to constrain new models as pursued here.
A direct quantitative implication of non-VZ neuronal production is that additional PCD will be required to offset additional neuronal production. Reports of previous models contend that NPC death is an insignificant component of neurodevelopment [17-19] despite empirical data that are consistent with significant NPC death. Data supporting much higher levels of NPC death than proposed in prior models were first reported using a sensitive DNA end-labeling technique, 'in situ end-labeling plus' (ISEL+), and ligation-mediated PCR [20-24]. Additional support came from recent analyses of NPC progeny, marked using a genomically encoded lineage tracer; here, the progeny clone size was found to diminish markedly at E14 . Perhaps most compelling, deletion of numerous pro-cell death genes, including those encoding caspase 3 , caspase 9 , APAF1 , Bax, Bak , and Pten , as well as novel molecules like ephrins , all lead to brain overgrowth phenotypes. Conversely and consistently, null mutations in pro-survival genes (for example, those encoding Bcl-x , Survivin  and Mcl-1 ) lead to smaller brains. Given this large body of empirical evidence, new models should account for the more extensive operation of PCD during cortical neurogenesis.
Here we report new quantitative models that incorporate new data and are consistent with cortical PCD empirical evidence. In particular, with respect to total neuronal production, we demonstrate a clear requirement for substantial NPC death during mouse cerebral cortical development. These models further provide a quantitative explanation of neurodevelopmental cortical overgrowth phenotypes produced by PCD attenuation as observed in caspase 3-null mutant mice.
Using cumulative labeling of newly synthesized DNA and explicit counting of each labeled cell migrating from the VZ and subventricular zone, the fraction of terminally differentiating cells has been measured experimentally [14,18]. We denote these fractions qi, where i = 1, 2,..., 11 for the 11 cell cycles (CCs) estimated to occur during the neurogenic interval (Table 1). These data, qi, and the estimated initial population of NPCs (P0), allow us to determine, through mathematical modeling, what fractions of NPC death at each CC, di for I = 1, 2,..., 11, yield VZ output consistent with experimental counts of the total number of excitatory neurons in the mature mouse cerebral cortex. To this end, we require an estimate for the plausible range of total VZ neuronal output.
The total number of neurons in the adult mouse cerebral cortex has been estimated at between 1.0 × 107 and 1.6 × 107 (Table 1 and references therein), and we use this estimate to constrain simulations of VZ output. Approximately 15 to 30% of the estimated total neuronal population consists of inhibitory interneurons. In addition, the estimate of 1.0 to 1.6 × 107 neurons includes cell loss that is due to 30 to 50% post-mitotic (non-VZ) cell death (Table 1 and references therein). Taking these figures into consideration, we can determine a range for plausible VZ output as follows. First, a lower bound on VZ output corresponds to taking the low estimate of 1.0 × 107 neurons and assuming that 30% of these are interneurons and that there was only 30% post-mitotic death; that is, a lower bound on VZ output is ((1.0 × 107) × (1 – 0.3))/(1 – 0.3) = 1.0 × 107. Similarly, an upper bound on VZ output follows from taking 1.6 × 107 neurons and assuming that 15% are interneurons and that 50% post-mitotic cell death occurred; thus, a plausible upper bound on VZ output is ((1.6 × 107) × (1 – 0.15))/(1 – 0.5)) = 2.72 × 107.
Therefore, data suggest that an accurate quantitative model of mouse cerebral cortical neurogenesis should yield a plausible range of VZ output between 1.0 × 107 and 2.72 × 107. However, the published estimate of 140 progeny per founding NPC  corresponds to 7.7 × 107 NPC progeny produced by 5.5 × 105 founding NPCs. NPC death was presumed negligible when calculating this progeny-per-NPC estimate and it cannot be reconciled with the maximum plausible VZ output.
To address this apparent contradiction, we calculated VZ output at various levels of NPC death using two related models of daughter cell fate decision-making (Figure 1A,B; Materials and methods). Model DG1 is derived from the experimental observation that NPC cell death occurs during the G1 phase of the CC , whereas model DG2 is more similar to existing models [17,18]. In both, the residual proliferative population before the ithCC, denoted Pi-1, is doubled at mitosis. In model DG1 (Figure 1A), no cells die prior to cell division while in model DG2, death is imposed before each division. As a simple consequence, when di is uniformly constant, model DG1 with founding population P0 is equivalent to model DG1 with founding population P0/(1-di). Thus, without NPC death (di = 0 for all i), the same VZ output is calculated for each model: 8.4 × 107, or 153 progeny per NPC. Assuming di = 0.05 for all i, as measured using the less-sensitive technique terminal dUTP nick-end labeling (TUNEL) , VZ output for DG1 and DG2 are 5.8 × 107 and 5.5 × 107respectively, both more than twice the plausible upper bound on VZ output (2.72 × 107) established above.
Caviness and colleagues  calculated that a founding population of 2.5 × 105 NPCs was compatible with their model. The plausibility window when P0 = 2.5 × 105 accommodates between 5 and 17% NPC death for model DG1 and between 4 and 15% NPC death for model DG2 (data not shown). These calculations suggest that a significant reduction of the founding NPC pool could be consistent with lower levels of cell death during development; however, this explanation is not consistent with two independent measurements of the founding NPC pool size [15,16].
Our models allow us to view the range of plausible VZ output as a function of PCD over the course of 11 CCs, and we refer to this correspondence as the 'plausibility window.' For example, as shown in Figure 1C,D, a plausible range of VZ output is observed when 13% <di < 26% for all i. Because of the delayed depletion of the proliferative population, the corresponding range of plausible di values in model DG1 is broader than that calculated using model DG2 (that is, 15 to 26% versus 13 to 23%, respectively). Neither model DG1 or DG2 matches the experimental measurements of either 5% NPC death using TUNEL or 50% NPC death using ISEL+; however, these calculations do demonstrate that 5% NPC death is too low to calculate VZ output adequately in the normal mouse brain.
Measurements of CC duration together with the fraction of NPCs that differentiate at experimentally defined ages permits extrapolation of model parameter qi for each i . Sensitivity analysis  provides a means of determining the relationship between inherent uncertainty in these estimates of qi and uncertainty in simulated VZ output.
Using a Monte Carlo approach, we sampled each qi and di as normally distributed random variables with means and taken from [14,37] and variance given by and (Materials and methods). For example, for and ν1 = 0.5 we have qi ~N(0.42, 0.122). Sampling the 22-dimensional parameter space for qi and di, and then computing VZ output for each sample point in parameter space, allowed us to determine that ν1 = ν2 = 0.5 gives plausible VZ output for >90% of the model realizations (Figure 2A).
For ν1 = ν2 = 0.5, we computed the first-order sensitivities of VZ output to variation in each qi and di (Materials and methods). The first-order sensitivities and are measures of how uncertainty in each qi and di drives uncertainty in VZ output. For example, when is close to zero for a given qi, significant variation in VZ output can still occur despite fixing this qi. We find that death rates in the first few CCs and differentiation rates in the middle CCs can account for the most significant portion of model variability; yet even the largest and are only approximately 10% (Table 2). As such, variation in VZ output is hardly attributable to variation in any one qi or di. An obvious explanation is that, given exponential growth conditions, outlying values of either qi or di in early CCs are easily compensated by outliers in later CCs. Interestingly, significant variation in qi and di from the mean values and can still yield plausible VZ output (Figure 2B).
NPC mitoses occur on the ventricular surface, but 'non-surface' mitoses are also observed in the developing cortex. Initially, these cell divisions were erroneously considered non-neuronogenic , but proliferative intermediate progenitor cells (IPCs) – daughter cells of NPCs that have migrated to the subventricular zone and intermediate zones, and are immunoreactive for the transcription factor Tbr2 – are now known to contribute additional excitatory neurons to the mouse cerebral cortex [11-13].
To determine how IPCs may affect the level of NPC death required for plausible VZ output, IPC subpopulation levels were extrapolated from the literature (Table 1 and references therein) and these constraints were incorporated into both models. We take a fraction of the Q cell population after the ith division – that is, we denote qi'Qi as the size of an IPC subpopulation – and these cells undergo an additional CC, doubling their contribution to the total Q population (Figure 3A,B). In initial modeling, we assume that additional IPC CCs and IPC death either do not occur or at least offset each other, and we include those IPCs generated after the last NPC CC in the VZ output calculated for CC11. With the inclusion of IPCs, VZ output without NPC death increases to 1.1 × 109 in model DG1 and more than doubles to 2.1 × 109 in model DG2 (Figure 3C,D). With 5% NPC death, model DG1 calculates VZ output to be 7.3 × 108 and model DG2 calculates 1.3 × 109 (Figure 3C,D). Since the qi' fraction is largest in the last three CCs, IPCs have a relatively small effect on the plausibility window; it expands slightly to accommodate 17 to 29% NPC death in model DG1 and 20 to 29% NPC death in model DG2. Importantly, a larger IPC population (from a larger IPC fraction or additional IPC CCs without IPC death) could increase, but not decrease, the requirement for NPC death.
Cortical development with 5% NPC death at early CCs requires more than 50% NPC death at later CCs
A primary objection to the idea of substantial NPC death is the intuition that ≥ 50% NPC death throughout neurodevelopment would preclude expansion of the NPC pool. However, the level of NPC death need not be constant at each CC. For example, ISEL+ labels approximately 5% of cells at E10 and the percentage increases significantly thereafter [20,35]. Moreover, NPC lineage analyses using a genomically encoded marker found that clone size increased during early development but then diminished after E14 , suggesting significant cell death with corpse elimination at later, but not earlier, CCs. These empirical observations are also consistent with increased model sensitivity to the levels of death at earlier CCs (Table 2). Together, these data suggest that early 'expansion' CCs occur; initially, low PCD levels – for example, setting di = 0.05 for 1 ≤ i ≤ 4 provides for four expansion CCs – further constrain VZ output. Here we calculate the corresponding amount of NPC death required during later CCs for plausible VZ output.
The fourth neurogenic CC takes place on E12, a time when significant ISEL+ labeling was observed . VZ output with three, four, or five expansion CCs amplifies differences between models DG1 and DG2. In model DG1 with three expansion CCs, the plausibility window permits 21 to 42% NPC death; four expansion CCs permit 26 to 56% NPC death; and five expansion CCs permit 36 to >80% NPC death (Figure 4A,C,E). Using model DG2, the plausibility window is less broad and accommodates less NPC death. After three expansion CCs the DG2 plausibility window permits 17 to 34% NPC death; after four expansion CCs it permits 20 to 41% NPC death; and after five expansion CCs it permits 25 to 54% NPC death (Figure 4B,D,F).
We note here that inclusion of four or five expansion CCs makes model DG1 and model DG2 readily distinguishable from one another. This distinction is a consequence of the differential involvement of cell death when expansion CCs begin and end, and illustrates some additional sensitivity to the size of the founding NPC population. As described above, during the first CC in model DG2 the founding population is reduced before mitosis; therefore, for uniformly constant cell death, model DG2 with initial population P0 is equivalent to model DG1 with initial population d1P0. Although this distinction leads to only subtle differences in most simulations, when considering the relationship between expansion CCs and plausible levels of NPC death, model DG1 with four or five expansion CCs illustrates additional potential for biological variation.
Early expansion CCs and IPCs together are most compatible with experimental data and further support high levels of NPC death
Both expansion CCs and IPC CCs occur during normal mouse brain development. When NPC population kinetics include expansion CCs (as shown in Figure 4) and IPCs (as shown in Figure 3), the plausibility window is broadest. For model DG1 with IPC CCs, three expansion CCs permit 26 to 46% NPC death and four permit 32 to 60% NPC death (Figure 5A,C). After five expansion CCs, VZ output already exceeds the minimum calculated estimate of cortical neurons (data not shown), with = 43% NPC death required during the remaining CCs to restrain VZ output (Figure 5E). As above (Figure 3), model DG2 with IPC CCs yields a lower and narrower range of NPC death levels: 28 to 41% with three expansion CCs, 32 to 49% with four expansion CCs, and 38 to 61% with five expansion CCs (Figure 5B,D,F).
Cell death amongst NPCs is a prominent feature of neurogenesis in other regions of the nervous system (for example, retina; reviewed in [38,39]); however, the amount of NPC death during mouse cerebral cortical development is debated [40,41]. Employing published cell counts to estimate ranges of neuronal population size in the mouse cerebral cortex, we calculated a plausible range of VZ output for normal mouse neurodevelopment. Models DG1 and DG2 use simple difference equations to capture essential features of previous probabilistic models. Model DG1 with three to five early expansion CCs, a pool of IPCs, and NPC death near 50% incorporates the most experimentally observed constraints and is consistent with NPC death levels as observed in ISEL+ analyses [22,42,43].
Apart from initial modeling by Takahashi, Caviness, and colleagues (referenced throughout), Gohlke, Faustman, and colleagues [17,44] have used a Kolmogorov forward equation to compute the probability distribution for mouse VZ output explicitly, as a continuous-time Markov chain. In contrast, we estimate this same probability distribution by repeatedly simulating a deterministic difference equation with random perturbations taken from assumed distributions. This allows us to illustrate how significant variation in qi and di fractions, over the course of 11 CCs, can occur and still yield plausible VZ output (Figure 2B). In population measurements using stereological counting, variation at the level of subpopulations of NPCs would probably go unnoticed, yet this may be an important feature of cerebral cortical development. While such variation is implicit in the Kolmogorov forward equation approach, Monte Carlo allows us to view sample trajectories that illustrate this variation.
The Gohlke models also observe that the original Takahashi models lead to VZ output that exceeds experimental counts of cortical neurons by at least threefold . However, these authors use two additional parameters (in addition to low-level NPC death) to reduce cortical neuron production to plausible levels: a diminished growth fraction insofar as not all VZ cells are progenitor cells and a clearance time for dying NPCs. Many NPCs co-label for bromodeoxyuridine and ISEL+ [20,24,42]; therefore, Gohlke and colleagues have likely recast some additional NPC death as a 'diminished' growth fraction. These parameters have the effect of reducing the proliferative population from which any given q fraction is taken. Consequently, the overproduction observed using the Takahashi model is limited and plausible levels of VZ output are obtained from the Gohlke mouse model [17,44].
The Gohlke mouse model is reportedly compatible with levels of NPC death up to 24% . This value is in good agreement with 13 to 26% NPC death calculated using models DG1 and DG2 (Figure 1). Despite this, Gohlke and colleagues [17,45] report subsequent model analysis using NPC death rates at or near 0% NPC death. The contrasting higher levels of NPC death required in Gohlke models of primate cortical neurogenesis , relative to murine cortical neurogenesis, may simply reflect an underestimation of mouse NPC death.
A strong implication of all modeling experiments presented here is that TUNEL significantly underestimates NPC death during mouse cerebral cortical development. This reflects technical differences in sensitivity between the two procedures, with ISEL+ being approximately ten times more sensitive than the originally reported TUNEL technique [21,22]. ISEL+ detects more dying cells not only amongst NPCs, but also in other tissues like the thymus and small intestinal villus . Consistent with a tenfold reduced sensitivity relative to ISEL+, TUNEL detects as few as one-tenth of the dying NPCs (5% versus 50%).
In order to reconcile the approximately 5% NPC death proposed by prior models  with the forebrain overgrowth phenotype observed in caspase 3-deficient mice, it has been proposed that NPC death may occur normally in a small population of neuroepithelial stem cells or radial glia at an early age (<E12) . The claim is that additional survival of a few such cells could underlie forebrain overgrowth in caspase 3-deficient mice because each of these individual cells might ultimately give rise to many neurons (approximately 140 according to estimates from concurrent models). However, this notion is inconsistent with experimental data demonstrating marked overgrowth of NPCs in caspase 3-null embryos by E12, before large numbers of NPC progeny emigrate to the cortex [26,42].
Our models are consistent with changes in NPC and neuronal populations that have been observed following experimental attenuation of PCD. For example, a 20% reduction in NPC death from the midpoint of any plausibility window corresponds to VZ output that exceeds the upper bound for that plausibility window. This is strikingly consistent with a 30% reduction in ISEL+ labeling of NPCs observed in caspase 3-deficient mice , where exceptional forebrain overgrowth is observed. We suggest that an in vivo correlate of 'exceeding the plausibility window' is forebrain overgrowth.
The breadth of the plausibility window provides a comparative measure of model robustness with respect to viable levels of NPC death. In some scenarios the slope of VZ output is steep (Figure 1D), and a 10% difference in NPC death leads to marked differences in VZ output, while a similar 10% change has little impact on VZ output in models where the slope is less steep (for example, Figure 4E). Notably, early expansion CCs (Figure 4), rather than the size of the founding NPC population or additional IPC progeny (compare Figures 2 and 3), significantly extend the range of NPC death that is compatible with plausible VZ output. Although modeled during the first three to five CCs here because of experimental evidence for pre-E12 expansion CCs [20,25], transient low-level NPC death could theoretically operate during any CC. It is tempting to speculate that high levels of NPC death at later CCs follow from low levels of NPC death during early expansion CCs. A similar process has been reported in embryonic stem cells earlier in development [47,48] and may be related to a permissive decatenation checkpoint observed in NPCs .
One might predict that other perturbations leading to transient high levels of NPC death at early CCs could lead to low level NPC death at later CCs. In this scenario, developmental accommodation of atypical NPC death might occur at the level of stem cell niches [49,50], providing sufficient, but not necessarily 'normal,' VZ output. Perhaps local control of NPC death could insulate cerebral cortical development against genetic differences and chemical or environmental insults. Given genetic diversity among NPCs, produced in part by chromosomal aneuploidy  and retrotransposition , differences in NPC death amongst individuals suggests selection and/or survival mechanisms that influence the mosaic composition of an individual's cerebral cortex.
Models DG1 and DG2 resolve discrepancies existing between previous models and experimental data; furthermore, these models provide a quantitative account for qualitative differences observed during PCD-attenuated cerebral cortical development. This theoretical framework should motivate additional experimental investigation of expansion CCs and reinterpretation of other cortical development phenotypes that measured PCD among NPCs using only TUNEL staining.
The order in which NPC death, differentiation, and proliferation are imposed alters VZ output, so it is natural to consider two related models. For the first model, DG1, the number of generated neurons at the ith CC is Qi = 2qiPi-1, where Pi-1 is given by Pi = (1 - di)(1 - qi)2Pi-1. Here, (1 - qi)2Pi-1 is the number of non-emigrating daughter cells after the ith division and di(1 - qi)2Pi-1 of these die, leaving Pi(1 - di)(1 - qi)2Pi-1 to divide again. Alternatively, imposing death before imposing differentiation gives Qi = qi(1 - di) 2Pi-1 while, again, Pi(1 - di)(1 - qi)2Pi-1 many NPCs go through the ith CC. We refer to this model as DG2. In each of these two related models, the total VZ output over 11 NPC CCs is . Viewing VZ = VZ(q, d) as a function of all qi and di – that is, q = (qi) and d = (di) are 11-dimensional vectors – allows us to explore how total VZ output is sensitive to variation in qi and di.
Defining and as the mean differentiation and death rates, we generate random values q = (qi) and d = (di) using a normal distribution. As variances we use and so that when either qi or di are close to 0 or 1, sampling rarely provides values that are nonbiological; that is, less than 0 or greater than 1. To indicate this distribution, we write and . Then, following Saltelli et al. , we define first-order sensitivities:
where Var [VZ(q, d)] is the variance in VZ output that arises when all qi and di are randomized (for example, 10,000 runs) and where Var [E[VZ(q, d|qk)]], for example, is the variance in expected VZ output that arises when, for some fixed k, all qi with i ≠ k are allowed to vary. That is, to compute Var [E[VZ(q, d|qk)]], qk is fixed repeatedly while for all i ≠ k and for all i = 1, 2,..., 11. Then, the average (expected) VZ output is taken, and this entire process is repeated (for example, approximately 100 times) to determine the variance of such an average.
Calculations were performed using MatLab version 7.7.0 (Mathworks, Natick, MA, USA) as detailed in the text. Figures were prepared using Illustrator and Photoshop (Adobe Systems Inc., San Jose, CA, USA). The MatLab scripts are available as additional files 1 and 2.
CC: cell cycle; E: embryonic day; IPC: intermediate progenitor cell; ISEL+: in situ end-labeling plus; NPC: neural progenitor/precursor cell; PCD: programmed cell death; TUNEL: terminal dUTP nick-end labeling; VZ: ventricular zone.
Click here for additional data file (1749-8104-4-28-S1.zip)
Additional file 2
Click here for additional data file (1749-8104-4-28-S2.zip)
We are very grateful to Drs S Rehen, D Kaushal, M Kingsbury, A Yang, L Cai, L Chao, T Haydar, and P Yeh for helpful discussions. This work was supported by a NIGMS pharmacology training grant (MM), NSF/EPSCoR grant no. EPS 0447660 (HM), and MH51699 (JC).
|Gotz M,Huttner WB. The cell biology of neurogenesisNat Rev Mol Cell Biol 2005;6:777–788. [pmid: 16314867] [doi: 10.1038/nrm1739]|
|Capela A,Temple S. LeX is expressed by principle progenitor cells in the embryonic nervous system, is secreted into their environment and binds Wnt-1Dev Biol 2006;291:300–313. [pmid: 16458284] [doi: 10.1016/j.ydbio.2005.12.030]|
|Fishell G,Kriegstein A. Cortical development: new conceptsNeuron 2005;46:361–362. [pmid: 15882631] [doi: 10.1016/j.neuron.2005.04.016]|
|Bystron I,Blakemore C,Rakic P. Development of the human cerebral cortex: Boulder Committee revisitedNat Rev Neurosci 2008;9:110–122. [pmid: 18209730] [doi: 10.1038/nrn2252]|
|Angevine JB,Sidman RL. Autoradiographic study of cell migration during histogenesis of the cerebral cortex of the mouseNature 1961;192:766–768. [pmid: 17533671] [doi: 10.1038/192766b0]|
|Takahashi T,Nowakowski RS,Caviness VS Jr. The cell cycle of the pseudostratified ventricular epithelium of the embryonic murine cerebral wallJ Neurosci 1995;15:6046–6057. [pmid: 7666188]|
|Bayar SA,Altman J. Neocortical Development. 1991New York: Raven Press Ltd;|
|Rakic P. A small step for the cell, a giant leap for mankind: a hypothesis of neocortical expansion during evolutionTrends Neurosci 1995;18:383–388. [pmid: 7482803] [doi: 10.1016/0166-2236(95)93934-P]|
|Kriegstein A,Noctor S,Martinez-Cerdeno V. Patterns of neural stem and progenitor cell division may underlie evolutionary cortical expansionNat Rev Neurosci 2006;7:883–890. [pmid: 17033683] [doi: 10.1038/nrn2008]|
|Takahashi T,Nowakowski RS,Caviness VS Jr. Early ontogeny of the secondary proliferative population of the embryonic murine cerebral wallJ Neurosci 1995;15:6058–6068. [pmid: 7666189]|
|Haubensak W,Attardo A,Denk W,Huttner WB. Neurons arise in the basal neuroepithelium of the early mammalian telencephalon: a major site of neurogenesisProc Natl Acad Sci USA 2004;101:3196–3201. [pmid: 14963232] [doi: 10.1073/pnas.0308600100]|
|Noctor SC,Martinez-Cerdeno V,Ivic L,Kriegstein AR. Cortical neurons arise in symmetric and asymmetric division zones and migrate through specific phasesNat Neurosci 2004;7:136–144. [pmid: 14703572] [doi: 10.1038/nn1172]|
|Englund C,Fink A,Lau C,Pham D,Daza RA,Bulfone A,Kowalczyk T,Hevner RF. Pax6, Tbr2, and Tbr1 are expressed sequentially by radial glia, intermediate progenitor cells, and postmitotic neurons in developing neocortexJ Neurosci 2005;25:247–251. [pmid: 15634788] [doi: 10.1523/JNEUROSCI.2899-04.2005]|
|Takahashi T,Nowakowski RS,Caviness VS Jr. The leaving or Q fraction of the murine cerebral proliferative epithelium: a general model of neocortical neuronogenesisJ Neurosci 1996;16:6183–6196. [pmid: 8815900]|
|Haydar TF,Nowakowski RS,Yarowsky PJ,Krueger BK. Role of founder cell deficit and delayed neuronogenesis in microencephaly of the trisomy 16 mouseJ Neurosci 2000;20:4156–4164. [pmid: 10818151]|
|Vaccarino FM,Schwartz ML,Raballo R,Rhee J,Lyn-Cook R. Fibroblast growth factor signaling regulates growth and morphogenesis at multiple steps during brain developmentCurr Top Dev Biol 1999;46:179–200. [pmid: 10417880] [doi: 10.1016/S0070-2153(08)60329-4]|
|Gohlke JM,Griffith WC,Faustman EM. The role of cell death during neocortical neurogenesis and synaptogenesis: implications from a computational model for the rat and mouseBrain Res Dev Brain Res 2004;151:43–54. [pmid: 15246691] [doi: 10.1016/j.devbrainres.2004.03.020]|
|Caviness VS Jr,Takahashi T,Nowakowski RS. Neuronogenesis and the early events of neocortical histogenesisResults Probl Cell Differ 2000;30:107–143. [pmid: 10857187]|
|Cai L,Hayes NL,Takahashi T,Caviness VS Jr,Nowakowski RS. Size distribution of retrovirally marked lineages matches prediction from population measurements of cell cycle behaviorJ Neurosci Res 2002;69:731–744. [pmid: 12205666] [doi: 10.1002/jnr.10398]|
|Blaschke AJ,Staley K,Chun J. Widespread programmed cell death in proliferative and postmitotic regions of the fetal cerebral cortexDevelopment 1996;122:1165–1174. [pmid: 8620843]|
|Yung YC,Kennedy G,Chun J. Identification of neural programmed cell death through the detection of DNA fragmentation in situ and by PCRCurr Protoc Neurosci 2009;Chapter 3 [pmid: 19575470]|
|Chun J,Blaschke AJ. Identification of neural programmed cell death through the detection of DNA fragmentation in situ and by PCRCurr Protoc Neurosci 2001;Chapter 3 [pmid: 18428472]|
|Staley K,Blaschke AJ,Chun J. Apoptotic DNA fragmentation is detected by a semi-quantitative ligation-mediated PCR of blunt DNA endsCell Death Differ 1997;4:66–75. [pmid: 16465212] [doi: 10.1038/sj.cdd.4400207]|
|Blaschke AJ,Weiner JA,Chun J. Programmed cell death is a universal feature of embryonic and postnatal neuroproliferative regions throughout the central nervous systemJ Comp Neurol 1998;396:39–50. [pmid: 9623886] [doi: 10.1002/(SICI)1096-9861(19980622)396:1<39::AID-CNE4>3.0.CO;2-J]|
|Wilkie AL,Jordan SA,Sharpe JA,Price DJ,Jackson IJ. Widespread tangential dispersion and extensive cell death during early neurogenesis in the mouse neocortexDev Biol 2004;267:109–118. [pmid: 14975720] [doi: 10.1016/j.ydbio.2003.11.006]|
|Kuida K,Zheng TS,Na S,Kuan C,Yang D,Karasuyama H,Rakic P,Flavell RA. Decreased apoptosis in the brain and premature lethality in CPP32-deficient miceNature 1996;384:368–372. [pmid: 8934524] [doi: 10.1038/384368a0]|
|Kuida K,Haydar TF,Kuan CY,Gu Y,Taya C,Karasuyama H,Su MS,Rakic P,Flavell RA. Reduced apoptosis and cytochrome c-mediated caspase activation in mice lacking caspase 9Cell 1998;94:325–337. [pmid: 9708735] [doi: 10.1016/S0092-8674(00)81476-2]|
|Yoshida H,Kong YY,Yoshida R,Elia AJ,Hakem A,Hakem R,Penninger JM,Mak TW. Apaf1 is required for mitochondrial pathways of apoptosis and brain developmentCell 1998;94:739–750. [pmid: 9753321] [doi: 10.1016/S0092-8674(00)81733-X]|
|Lindsten T,Golden JA,Zong WX,Minarcik J,Harris MH,Thompson CB. The proapoptotic activities of Bax and Bak limit the size of the neural stem cell poolJ Neurosci 2003;23:11112–11119. [pmid: 14657169]|
|Groszer M,Erickson R,Scripture-Adams DD,Lesche R,Trumpp A,Zack JA,Kornblum HI,Liu X,Wu H. Negative regulation of neural stem/progenitor cell proliferation by the Pten tumor suppressor gene in vivoScience 2001;294:2186–2189. [pmid: 11691952] [doi: 10.1126/science.1065518]|
|Depaepe V,Suarez-Gonzalez N,Dufour A,Passante L,Gorski JA,Jones KR,Ledent C,Vanderhaeghen P. Ephrin signalling controls brain size by regulating apoptosis of neural progenitorsNature 2005;435:1244–1250. [pmid: 15902206] [doi: 10.1038/nature03651]|
|Motoyama N,Wang F,Roth KA,Sawa H,Nakayama K,Negishi I,Senju S,Zhang Q,Fujii S,et al. Massive cell death of immature hematopoietic cells and neurons in Bcl-x-deficient miceScience 1995;267:1506–1510. [pmid: 7878471] [doi: 10.1126/science.7878471]|
|Jiang Y,de Bruin A,Caldas H,Fangusaro J,Hayes J,Conway EM,Robinson ML,Altura RA. Essential role for survivin in early brain developmentJ Neurosci 2005;25:6962–6970. [pmid: 16049172] [doi: 10.1523/JNEUROSCI.1446-05.2005]|
|Arbour N,Vanderluit JL,Le Grand JN,Jahani-Asl A,Ruzhynsky VA,Cheung EC,Kelly MA,MacKenzie AE,Park DS,Opferman JT,Slack RS. Mcl-1 is a key regulator of apoptosis during CNS development and after DNA damageJ Neurosci 2008;28:6068–6078. [pmid: 18550749] [doi: 10.1523/JNEUROSCI.4940-07.2008]|
|Thomaidou D,Mione MC,Cavanagh JF,Parnavelas JG. Apoptosis and its relation to the cell cycle in the developing cerebral cortexJ Neurosci 1997;17:1075–1085. [pmid: 8994062]|
|Saltelli A,Tarantola S,Campolongo F,Ratto M. Sensitivity Analysis in Practice: a Guide to Assessing Scientific Models. 2004Hoboken, NJ: Wiley;|
|Miyama S,Takahashi T,Nowakowski RS,Caviness VS Jr. A gradient in the duration of the G1 phase in the murine neocortical proliferative epitheliumCereb Cortex 1997;7:678–689. [pmid: 9373022] [doi: 10.1093/cercor/7.7.678]|
|de la Rosa EJ,de Pablo F. Cell death in early neural development: beyond the neurotrophic theoryTrends Neurosci 2000;23:454–458. [pmid: 11006461] [doi: 10.1016/S0166-2236(00)01628-3]|
|Lossi L,Merighi A. In vivo cellular and molecular mechanisms of neuronal apoptosis in the mammalian CNSProg Neurobiol 2003;69:287–312. [pmid: 12787572] [doi: 10.1016/S0301-0082(03)00051-0]|
|Gilmore EC,Nowakowski RS,Caviness VS Jr,Herrup K. Cell birth, cell death, cell diversity and DNA breaks: how do they all fit together?Trends Neurosci 2000;23:100–105. [pmid: 10675909] [doi: 10.1016/S0166-2236(99)01503-9]|
|Chun J. Cell death, DNA breaks and possible rearrangements: an alternative viewTrends Neurosci 2000;23:407–409. [pmid: 10941186] [doi: 10.1016/S0166-2236(00)01612-X]|
|Pompeiano M,Blaschke AJ,Flavell RA,Srinivasan A,Chun J. Decreased apoptosis in proliferative and postmitotic regions of the Caspase 3-deficient embryonic central nervous systemJ Comp Neurol 2000;423:1–12. [pmid: 10861532] [doi: 10.1002/1096-9861(20000717)423:1<1::AID-CNE1>3.0.CO;2-S]|
|Pompeiano M,Hvala M,Chun J. Onset of apoptotic DNA fragmentation can precede cell elimination by days in the small intestinal villusCell Death Differ 1998;5:702–709. [pmid: 10200526] [doi: 10.1038/sj.cdd.4400403]|
|Gohlke JM,Griffith WC,Bartell SM,Lewandowski TA,Faustman EM. A computational model for neocortical neuronogenesis predicts ethanol-induced neocortical neuron number deficitsDev Neurosci 2002;24:467–477. [pmid: 12697984] [doi: 10.1159/000069357]|
|Gohlke JM,Griffith WC,Faustman EM. Computational models of neocortical neuronogenesis and programmed cell death in the developing mouse, monkey, and humanCereb Cortex 2007;17:2433–2442. [pmid: 17204816] [doi: 10.1093/cercor/bhl151]|
|Haydar TF,Kuan CY,Flavell RA,Rakic P. The role of cell death in regulating the size and shape of the mammalian forebrainCereb Cortex 1999;9:621–626. [pmid: 10498280] [doi: 10.1093/cercor/9.6.621]|
|Mantel C,Guo Y,Lee MR,Kim MK,Han MK,Shibayama H,Fukuda S,Yoder MC,Pelus LM,Kim KS,Broxmeyer HE. Checkpoint-apoptosis uncoupling in human and mouse embryonic stem cells: a source of karyotpic instabilityBlood 2007;109:4518–4527. [pmid: 17289813] [doi: 10.1182/blood-2006-10-054247]|
|Damelin M,Sun YE,Sodja VB,Bestor TH. Decatenation checkpoint deficiency in stem and progenitor cellsCancer Cell 2005;8:479–484. [pmid: 16338661] [doi: 10.1016/j.ccr.2005.11.004]|
|Shen Q,Goderie SK,Jin L,Karanth N,Sun Y,Abramova N,Vincent P,Pumiglia K,Temple S. Endothelial cells stimulate self-renewal and expand neurogenesis of neural stem cellsScience 2004;304:1338–1340. [pmid: 15060285] [doi: 10.1126/science.1095505]|
|Jiao JW,Feldheim DA,Chen DF. Ephrins as negative regulators of adult neurogenesis in diverse regions of the central nervous systemProc Natl Acad Sci USA 2008;105:8778–8783. [pmid: 18562299] [doi: 10.1073/pnas.0708861105]|
|Rehen SK,McConnell MJ,Kaushal D,Kingsbury MA,Yang AH,Chun J. Chromosomal variation in neurons of the developing and adult mammalian nervous systemProc Natl Acad Sci USA 2001;98:13361–13366. [pmid: 11698687] [doi: 10.1073/pnas.231487398]|
|Muotri AR,Chu VT,Marchetto MC,Deng W,Moran JV,Gage FH. Somatic mosaicism in neuronal precursor cells mediated by L1 retrotranspositionNature 2005;435:903–910. [pmid: 15959507] [doi: 10.1038/nature03663]|
|Cobos I,Calcagnotto ME,Vilaythong AJ,Thwin MT,Noebels JL,Baraban SC,Rubenstein JL. Mice lacking Dlx1 show subtype-specific loss of interneurons, reduced inhibition and epilepsyNat Neurosci 2005;8:1059–1068. [pmid: 16007083] [doi: 10.1038/nn1499]|
|Sarkisian MR,Frenkel M,Li W,Oborski JA,LoTurco JJ. Altered interneuron development in the cerebral cortex of the flathead mutantCereb Cortex 2001;11:734–743. [pmid: 11459763] [doi: 10.1093/cercor/11.8.734]|
|Verney C,Takahashi T,Bhide PG,Nowakowski RS,Caviness VS Jr. Independent controls for neocortical neuron production and histogenetic cell deathDev Neurosci 2000;22:125–138. [pmid: 10657705] [doi: 10.1159/000017434]|
Model parameters and constraints
|Founding NPC population||P0||5 to 6 × 105||[15,16]|
|NPC CC number||i||1 to 11|||
|Fraction of NPCs that become newly||qi||q1 = 0.005; q2 = 0.04;||[14,37]|
|post-mitotic neurons||q3 = 0.09;q4 = 0.14;|
|q5 = 0.21;q6 = 0.31;|
|q7 = 0.42; q8 = 0.54;|
|q9 = 0.69;q10 = 0.84;|
|q11 = 1.0|
|Fraction of NPCs that die at ith CC||di||5% to >50%||[12,13,17]|
|After ith CC, number of new:|
|Post-mitotic cells||QI||Model specific,|
|Dying daughter cells||DI|
|Total neurons produced from NPCs||1 × 107 to 2.72 × 107||See text|
|Cortical neurons||10 to 16 × 106||[53,54]|
|Cortical interneurons||15 to 30% of cortical neurons||[20,24,35]|
|Post-natal cell death||30 to 50% of cortical neurons||[20,24,55]|
|Fraction of IPCs at each CC||q'i||q'1 = 0; q'2 = 0; q'3 = 0;||[10,12,13]|
|q'4 = 0; q'5 = 0.12;|
|q'6 = 0.12; q'7 = 0.16;|
|q'8 = 0.16; q'9 = 0.54;|
|q'10 = 0.54; q'11 = 0.54|
Representative sensitivity analysis
Previous Document: Relationships among neurocognition, symptoms and functioning in patients with schizophrenia: a path-...
Next Document: Three different Plasmodium species show similar patterns of clinical tolerance of malaria infection.