The brain is one of the most fascinating organs in the human body. Its complex structure on both micro- and macroscopic scales closely correlates with the unique cognitive abilities of humans. Cortical folding is one of the most important features of the human brain. Still, compared to other mammals, the human brain is neither the largest nor the most folded brain. However, relative to its size, it has the largest number of cortical neurons that connect with billions of neuronal synapses (Herculano-Houzel, 2009). This fact attracted the attention of neuroscientists over the past few years to explore the source of these cells and how they develop in the early stages of brain development.

The number of brain cells is determined in utero through the proliferation process. Previous studies on different lissencephalic species, such as mice, have shown that cell division in the brain is confined to a small region near the cerebral ventricles (Hansen et al., 2010). However, in gyrencephalic species, this seems to be different. Recent findings show that the human brain, for example, is characterized by two proliferation zones with two different types of progenitor cells. Both zones produce neurons that later migrate towards the outer brain surface and form the cortex Lui et al., 2011; Pebworth et al., 2021.

In rodents, progenitor cells around the ventricular zone (VZ) generate intermediate progenitor cells as their daughters, which accumulate above the VZ and form a new layer called the subventricular zone (Noctor et al., 2002). In humans, there is an additional outer layer of the subventricular zone, often referred to as outer subventricular zone (OSVZ) (Hansen et al., 2010; Lui et al., 2011; Noctor et al., 2007). This zone was first discovered in the monkey brain by Colette Dehay and her colleagues (Smart et al., 2002), and confirmed in the human brain by several following studies (Huttner and Kosodo, 2005). The OSVZ seems to play a significant role in the proliferation process and affects the size and complexity of the human cortex. The evidence for this allegation is the wave of cortical neurogenesis that coincides with the cell division in the OSVZ (Lukaszewicz et al., 2005). At the macroscopic scale, the high proliferation in the OSVZ coincides with a significant tangential expansion of the cortical layers. The latter is an essential factor for the formation of cortical folds (Reillo et al., 2011). Still, it remains unknown, how exactly this proliferation process in the OSVZ affects gyrification of the forming cortex. Different approaches have been used to understand the relation between cellular mechanisms at the microscopic scale and cortical development at the macroscopic scale. Genetic analyses and experimental studies using cell culture models and brain organoids have given first valuable insights concerning the source of cells and their behavior (Hansen et al., 2010). Here, we intend to complement these studies by using a numerical approach to bridge the scales from the behavior of different progenitor cell types at the cell scale to the emergence of cortical folds at the tissue or organ scale.

From a mechanics point of view, forces that are generated due to cellular processes may act as a link to understand the underlying mechanisms behind cortical folding (Budday et al., 2015b). Many previous studies tried to explain normal and abnormal cortical folding either from a purely biological or mechanical perspective (Tallinen et al., 2014; Razavi et al., 2015). However, it will not be possible to capture the folding mechanism without considering both perspectives at the same time (de Rooij and Kuhl, 2018; Zarzor et al., 2021; Wang et al., 2022). In other words, to fully understand the physiological and pathological mechanisms underlying cortical folding in the developing human brain, we need to study the coupling between cellular processes and mechanical forces – to eventually assess how disruption of cellular processes affect the folding pattern and lead to malformations of cortical development (Guerrini et al., 2008; Blumcke et al., 2021; Llinares-Benadero and Borrell, 2019).

To fill this knowledge gap, we establish a two-field computational model that accounts for both proliferating zones in the human brain, the VZ and OSVZ. The first field in the model describes the growth and deformation of brain tissue based on the theory of finite growth (Rodriguez et al., 1994; Göktepe et al., 2010). The second field describes the cellular processes occurring during human brain development, where we use an advection-diffusion equation to mimic the proliferation and migration in the subcortex and neuronal connectivity in the cortex (de Rooij and Kuhl, 2018; Zarzor et al., 2021). We add two source terms to consider the division in different zones. We validate our model through a comparison of the simulation results with histologically stained sections of the human fetal brain and address unresolved questions regarding the role of the OSVZ for cortical folding.

The central cellular unit that plays a critical role in essential processes of brain development is a type of progenitor cells called radial glial cells. In the early stage of brain development, when neurogenesis begins, neuroepithelial cells transform into radial glial cells (Noctor et al., 2007). Around gestational week 5, these cells locate near cerebral ventricles, in the VZ, where they undergo interkinetic nuclear migration. The associated symmetric division behavior leads to a significant increase in the number of radial glial cells and results in both increased thickness and surface area of the VZ (Blows, 2003; Fish et al., 2008; Bystron et al., 2008). Subsequently, the cells switch to an asymmetric division behavior and generate intermediate progenitor cells (Noctor et al., 2004). The latter migrate to the subventricular zone, where they proliferate and produce neurons, as illustrated in Figure 1 (Noctor et al., 2007; Pebworth et al., 2021). Ultimately, the majority of cortical neurons are produced by intermediate progenitor cells (Lui et al., 2011; Libé-Philippot and Vanderhaeghen, 2021).

Figure 1

Download asset Open asset

According to the radial unit hypothesis proposed by Pasko Rakic over 30 years ago, the radial glial cell fibers organize the migration process, which starts around gestational week 6 (Nonaka-Kinoshita et al., 2013). He postulated that these fibers form a scaffold to guide neurons during their migration from the proliferating zones to their final destination in the cortex, which forms the outer brain surface (Rakic, 1988; Lui et al., 2011). In gyrencephalic species, those fibers have a characteristic fan-like distribution (Borrell and Götz, 2014; Nonaka-Kinoshita et al., 2013). However, it is still under debate whether this unique distribution is the cause or rather the result of cortical folding. Our recent computational analyses support the latter, i.e., that it is the result of cortical folding (Zarzor et al., 2021). The migration process synchronizes with a radial expansion of all brain layers. Still, the VZ does not expand remarkably as the intermediate progenitor cells move outwards to the subventricular zone. The migrated neurons finally organize themselves in the six-layered cortex in an inside-out sequence, where the early-born neurons occupy the inner layers (Gilmore and Herrup, 1997). Until gestational week 23, the outer brain surface is still smooth, although the bottom four layers of the cortex are already filled with neurons (Shinmyo et al., 2017). The first folds appear between gestational weeks 20 and 28, as the cortical layer significantly expands tangentially (Budday et al., 2015b; Habas et al., 2012). Importantly, at around gestational week 25, the cortical neuronal connectivity emerges and comes along with the horizontal elongation of neuronal dendrites (Takahashi et al., 2012).

While the processes summarized above are common among mammals, the human brain has some specific features that play a significant role in increasing the number of cortical neurons, and enhancing the complexity of cortical folds (Libé-Philippot and Vanderhaeghen, 2021). At the beginning of the second trimester, around gestational week 11, the original radial glial cells switch from producing intermediate progenitor cells to producing a special kind of cells that is found in all gyrencephalic species but is enriched in the human brain. The newly generated cells are similar to the radial glial cells in terms of shape and function, but unlike the original radial glial cells, they migrate to the outer layer of the subventricular zone, referred to as OSVZ, after they are born. Therefore, they are called outer radial glial cells (ORGCs) (Lui et al., 2011; Fietz et al., 2010; Hansen et al., 2010; Reillo et al., 2011; Nonaka-Kinoshita et al., 2013). We would like to note that some literature refers to this type of cells as basal radial glial cells. While the original radial glial cells have a bipolar morphology with two processes – one extending to the cerebral ventricle and one to the outer cortical surface – the ORGCs have a distinct unipolar structure with only a single process extending to the outer cortical surface (Hansen et al., 2010; Betizeau et al., 2013; Reillo et al., 2011; Nonaka-Kinoshita et al., 2013).

The OSVZ shows a significantly more pronounced radial expansion compared to the inner subventricular zone and VZ between gestational weeks 11.5 and 32. The immediate reason causing this difference is the characteristic division behavior of ORGCs: they translocate rapidly in radial direction before they divide, which scientists have referred to as ‘mitotic small translocation (MST)’ (Fietz et al., 2010). Importantly, the MST behavior pushes the boundary of the OSVZ outward, which increases its capacity to produce new neurons. The intermediate progenitor cells have enough space to undergo multiple rounds of division before producing neurons, which increases the overall number of generated neurons (Kriegstein et al., 2006; Lui et al., 2011).

The ORGCs, like radial glial cells, play an important role in the proliferation process: they divide symmetrically and asymmetrically to produce further ORGCs and intermediate progenitor cells (Libé-Philippot and Vanderhaeghen, 2021). Intermediate progenitor cells divide to generate a pair of neurons (Lui et al., 2011). According to previous studies, 40% of produced neurons are generated by ORGCs at gestational week 13, but this ratio increases to 60% by gestational week 14, and exceeds 75% by gestational week 15.5. Then, after gestational week 17, the ORGCs become the only source of cortical neurons in the upper cortical layers (Hansen et al., 2010). Besides their role in increasing the number of neurons, the ORGCs generate additional scaffolds that elongate to the outer brain surface and serve as paths for neuronal migration (Llinares-Benadero and Borrell, 2019; Nonaka-Kinoshita et al., 2013). Compared to other mammals, the neurogenesis of the human cortex can thus be divided into two main stages. The first stage is characterized by migration along a continuous scaffold consisting of radial glial cell fibers, which run from the ventricular surface to the outer cortical layer around gestational week 15. During the second stage, the migration path switches to a discontinuous form. After gestational week 17, the radial glial cell fibers run from the ventricular surface to the inner subventricular zone, while ORGC fibers run from the OSVZ to the outer brain surface, as indicated in Figure 1 (Nowakowski et al., 2016). Consequently, migrating neurons follow a sinuous path through numerous radial fibers before they reach its final location in the cortex (Lui et al., 2011). The additional scaffolds formed by ORGC fibers are not only important for neuronal migration, but also for the tangential expansion of the cortex. Previous studies show that a reduced number of ORGCs lead to reduced tangential expansion. In these cases, the cortex is less folded or even lissencephalic (Poluch and Juliano, 2015). In contrast, increasing the number of ORGCs leads to more excessive folding (Florio et al., 2017; Borrell, 2018). Still, it remains unknown whether these effects are a result of the specific proliferation behavior of ORGCs or the associated scaffold of ORGC fibers. What is known, though, is that the existence of ORGCs is a necessary but not sufficient condition for cortical folding (Llinares-Benadero and Borrell, 2019).

Around gestational week 28, the migrating neurons occupy the first (top) cortical layer, while the migration and proliferation processes come to an end. After finishing their role during the neurogenesis stage, radial glial cells and ORGCs switch to produce different types of glial cells, e.g., astrocytes and oligodendrocytes (Schmechel and Rakic, 1979). Also, radial glial cells may later convert into ependymal cells that locate around the cerebral ventricles. The oligodendrocytes form myelin sheaths around neuronal axons, wherefore the subcortical layer gains its characteristic white color.

To numerically study the effect of the VZ and the OSVZ on the resulting folding pattern, we simulate human brain development by using the finite element method. The influence of various factors on the emergence of cortical folds can be best shown on a simple two-dimensional (2D) quarter-circular geometry (Darayi et al., 2022), as illustrated in Figure 2A. In addition, we also investigate the folding evolution on a simplified half-sphere three-dimensional (3D) geometry. In the following, we introduce the main equations describing the coupling between cellular mechanisms in different proliferating zones and cortical folding, which we solve numerically.

Figure 2

Download asset Open asset

We validate our computational model by comparing the simulation results with histologically stained sections of the human fetal brain. For details on the corresponding preparation and staining, we refer to Zarzor et al., 2021. The sections belong to human fetuses aborted at gestational weeks 17, 24, 30, and 34. After defining the areas representative of the simulation domain introduced in the previous section (see Figure 2A) for each gestational week, we detect and assess the cell density using the software Qupath. This study was approved by the ethics review board of the Friedrich-Alexander-Universität Erlangen-Nürnberg with the reference number 22-209-Bp, and all procedures were conducted in accordance with the Declaration of Helsinki. In the following, we summarize the steps we followed to determine the cell density in human fetal brain sections. Since this analysis is intended to mainly serve as a means of comparison between different gestational weeks, a more detailed analysis using individual markers was not necessary at this stage.

1. Identify regions of interest to which the analysis should be applied, as shown in Figure 3A.

2. Pre-process the annotated area to ensure better cell detection through image color transfer, contrast ratio modification, and stain vector settings.

3. Automatically detect the cells in the relevant area by using the ‘positive cell detection’ command in Qupath, which distinguishes between cell types according to the staining. Here, the negative cells (blue) are mostly glial cells or extracellular matrix components, while positive cells (red) are mostly neurons or progenitor cells, as shown in Figure 3B. While this distinction between cell types might not be fully accurate, the results are satisfactory for our specific application.

4. Count the nearby detections in a small circle with an arbitrarily chosen radius $ϵ=100\mathrm{\mu }\mathrm{m}$ around each detected positive cell (i) to determine the corresponding cell density by dividing the number of nearby detections by the circular area, as demonstrated in Figure 3C.

5. Visualize the cell density, as illustrated in Figure 3A.

Figure 3

Download asset Open asset

Figure 3 shows the cell density distribution around gestational week 17 and demonstrates the densely packed VZ due to the high proliferation rate of radial glial cells, with about ${\displaystyle 13,600\text{cells}/\mathrm{m}\mathrm{m}{}^2}$. In addition, we can locate the other zones introduced in Figure 1. The OSVZ shows a higher cell density (approximately ${\displaystyle 9550\text{cell}/\mathrm{m}\mathrm{m}{}^2}$) than both the inner subventricular and the intermediate zone. This zone is composed of several types of cells, the original intermediate progenitor cells that migrated from the VZ, migrated neurons, ORGCs that produce more intermediate progenitor cells, and newborn neurons that are produced through the intermediate progenitor cells’ asymmetric division. Our analysis of the histologically stained sections of the human fetal brain (see Figure 3) also shows that the OSVZ is 1.5 times thicker than the inner subventricular zone, even though it did not emerge at an earlier stage of development (the inner subventricular zone emerges around gestational week 7 and the OSVZ around gestational week 11). This implies that the thickness of the OSVZ increases with time. The intermediate zone is characterized by a low cell density with about ${\displaystyle 1600\text{cells}/\mathrm{m}\mathrm{m}{}^2}$. Still, it is a transit area for the migrating neurons. The higher cell density in the inner subventricular zone ${\displaystyle 4800\text{cells}/\mathrm{m}\mathrm{m}{}^2}$ corroborates the presence of another type of cell besides migrating neurons, i.e., intermediate progenitor cells. The migration process in gestational week 17 is still ongoing – the cortex is not yet fully developed and filled with neurons with only about ${\displaystyle 12,000\text{cells}/\mathrm{m}\mathrm{m}{}^2}$. Therefore, at this stage of development, the VZ still has the highest cell density in the brain. Figure 4 compares the cell density distribution along a line from the ventricular surface to the outer cortical surface in the human fetal brain at gestational week 17 (indicated by line L in Figure 3) with different simulation results. For better comparability and to avoid differences in the dimensions between the histologically stained sections of the human fetal brain and the simulation domain, we normalize the domain’s radius according to the extension from the ventricular to the outer cortical surface in the stained sections. In addition, we normalize the cell density with respect to its maximum value in the cortex. According to our previous work (Zarzor et al., 2021), we include a varying stiffness in the cortex during human brain development in our simulations, adopt a stiffness ratio of 3, and a division rate in the VZ of 120.

Figure 4

Download asset Open asset

The simulation results for an initial division rate in the OSVZ of $G_{\text{osvz}}=10$ well capture the trends observed in the histologically stained sections of the human fetal brain. The cell density shows a first local peak representing the VZ for a normalized radius of approximately 0.05. It then gradually decreases to reach its first local minimum in the inner subventricular zone for a normalized radius of 0.3. This effect is less pronounced in the simulations than in the actual human fetal brain. The curves start to rise again in the OSVZ to reach the second peak for a normalized radius of 0.4. Again, the simulation results capture this peak quite accurately, with a difference of only 3%. The second local minimum represents the intermediate zone for a normalized radius of 0.7, while the third peak represents the cortex. The simulation results for $G_{\text{osvz}}=20$ and 30 result in a higher cell density in the OSVZ than in the actual human fetal brain. In contrast, the curve for $G_{\text{osvz}}=0$ (i.e., without including the effect of the OSVZ) shows a significantly decreased cell density in the OSVZ and intermediate zone.

In this section, we apply our computational model to answer some of the major questions regarding the role of the OSVZ for cortical folding during human brain development. As mentioned above, we consider both cases, constant and varying cortical stiffness. Before addressing the individual questions, we would like to first highlight some general features of the mechanical instability problem as the underlying mechanism of cortical folding. Figure 5 shows an exemplary temporal course of both sulcus depth and folding evolution – defined as the ratio between the outer perimeter at time step $t$ and the initial perimeter (see Figure 6) – for the varying stiffness case with $G_{\text{vz}}=120$ and $G_{\text{osvz}}=20$. The outer brain maintains its smooth surface until the first instability point, where compressive stresses in the cortex induced by cell-density-driven cortical growth reach a critical value, the cortex starts to fold and sulci deepen rapidly and uniformly. At the instability point, the curve of the folding evolution shows a kink, as the outer brain surface now expands more rapidly. After a second instability point, where a secondary instability occurs and we see a pitchfork-like bifurcation in Figure 5, only every second sulcus continues to deepen, while those in between become shallower again. We refer to the resulting folding pattern as period doubling pattern. For a detailed discussion of secondary instabilities, period doubling, and even period tripling patterns and their role in cortical folding, we refer to our previous works in Budday et al., 2015a and Budday et al., 2015c.

Figure 5

Download asset Open asset

Figure 6

Download asset Open asset

How does cell proliferation in the OSVZ affect cortical folding patterns?

Many previous studies have tried to address this point (Hansen et al., 2010). However, through purely experimental approaches, it is difficult to answer this question. In our model, we can apply different values of the division rate in the OSVZ ($G_{\text{osvz}}$) to show how ORGC proliferation affects cortical folding. Figure 7 shows the folding patterns emerging at gestational week 36 for both varying and constant stiffness cases and different values of the initial division rate in the OSVZ $G_{\text{osvz}}$. In the case of constant stiffness, the effect is marginal. In the case of varying stiffness, there is a more noticeable change in the folding patterns. In general, the distance between neighboring sulci decreases with increasing $G_{\text{osvz}}$, as marked in Figure 7. For the displayed cases, the distance decreases from d = $8.796\mathrm{mm}$ for $G_{\text{osvz}}=0$ to d = $8.67\mathrm{mm}$ for $G_{\text{osvz}}=10$ and d = $8.2\mathrm{mm}$ for $G_{\text{osvz}}=20$. Interestingly, the cortical thickness and effective stiffness ratio at the first instability point (denoted by w in Figure 5) are the same for all these cases. Therefore, we attribute the observed differences to the faster increase in the cell density and thus cortical growth, cortical stiffness, and the effective stiffness after the instability has been initiated. In addition, we observe period doubling patterns emerge (Budday et al., 2015a; Budday et al., 2015c), which are most pronounced at a value of $G_{\text{osvz}}=20$. This indicates that the proliferation in the OSVZ enhances secondary mechanical instabilities and leads to more complex folding patterns earlier.

Figure 7

Download asset Open asset

Figure 8 demonstrates that the observed trends also hold true when extending the model to 3D. For the case of varying stiffness with a stiffness ratio of 3, a growth ratio of 3, and an initial division rate in the VZ $G_{\text{vz}}=600$, the folding complexity increases with increasing initial division rate in the OSVZ $G_{\text{osvz}}$.

Figure 8

Download asset Open asset

Besides the direct relation between the proliferation in the OSVZ and the folding morphology, there are indirect effects and other aspects concerning ORGC proliferation, which will be discussed in more detail in the following sections.

How does cell proliferation in the OSVZ affect the cell density and folding evolution?

After we have assessed the effect of cell proliferation in the OSVZ on the final folding pattern, we investigate its effect on the evolution of both cell density and folding morphology. Figure 9 shows the temporal evolution of the maximum cell density in the domain and the folding evolution between time steps 260 and 380 for different initial division rates in the OSVZ $G_{\text{osvz}}$ and a constant division rate in the VZ $G_{\text{vz}}=120$. Again, we consider both cases, constant and varying cortical stiffness. Increasing the initial division rate in the OSVZ leads to a significant increase in both the cell density and folding evolution. Consequently, for the case of $G_{\text{osvz}}=30$, the cell density reaches the highest value of ${\displaystyle 1100\mathrm{m}{\mathrm{m}}^{-2}}$, corresponding to a folding evolution of 1.25. For the case of $G_{\text{osvz}}=0$, in contrast, the cell density does not even exceed the migration threshold. Comparing the cases $G_{\text{osvz}}=20$ and $G_{\text{osvz}}=30$ shows that the instability occurs earlier with increasing initial division rate in the OSVZ $G_{\text{osvz}}$. Thus, the ORGC proliferation decreases the time required to reach the final folding pattern. In general, the differences between the results for the constant and varying cortical stiffness case are minor. However, the curves for the varying cortical stiffness case rise faster than for the constant case. We note that we could not generate results for the initial division rate in the OSVZ $G_{\text{osvz}}=30$ and the constant cortical stiffness case due to numerical issues.

Figure 9

Download asset Open asset

Which proliferation zone is more influential for fetal human brain development?

To answer the question which one of the proliferation zones (VZ or OSVZ) is more important for cellular brain development and cortical folding, we have implemented four sets of division rates. The first set assumes a high initial division rate in the VZ $G_{\text{vz}}$ and a low initial division rate in the OSVZ $G_{\text{osvz}}$. For the following parameter sets, we gradually decrease $G_{\text{vz}}$, while we increase $G_{\text{osvz}}$. Figure 10 shows the corresponding results for the maximum cell density in the domain and the folding evolution between time steps 250 and 450 for both cases, constant and varying cortical stiffness. Unexpectedly, the cell density for the set ($G_{\text{vz}}=30,G_{\text{osvz}}=30$) rises faster than the one for the set ($G_{\text{vz}}=120,G_{\text{osvz}}=0$). Indeed, not only the cell density is affected by increasing the division rate in the OSVZ (at the expense of decreasing it in the VZ), but also convolutions (cortical folds) appear earlier, as illustrated by the curve for the folding evolution (Figure 10, right). Our results thus indicate an unproportionally strong effect of ORGC proliferation in the OSVZ on cellular brain development and cortical folding. We attribute this observation to the larger volume occupied by the OSVZ compared to the VZ. Concerning the difference between the results for the constant and varying cortical stiffness case, the curves for a varying stiffness rise faster than for a constant stiffness.

Figure 10

Download asset Open asset

How does the MST behavior of ORGCs affect brain development?

We have introduced a specific simulation parameter to capture the MST behavior of ORGCs, which now allows us to study the effect of the MST factor on cortical folding. For this parameter study, we limit ourselves to the varying cortical stiffness case with an initial division rate in the VZ $G_{\text{vz}}=30$ and an initial division rate in the OSVZ $G_{\text{osvz}}=30$. Figure 11 shows the temporal evolution of the maximum cell density and the folding evolution between time steps 250 and 430 for different values of the MST factor. Our simulations indicate that with increasing MST factor, the value of the maximum cell density in the domain increases exponentially. Since the MST factor incorporates the expansion of the OSVZ with time, these results are consistent with the observation in the previous section, showing an increasing spatial expansion of the OSVZ.

Figure 11

Download asset Open asset

Is the OSVZ affected by cortical folding?

After we have discussed the effect of the OSVZ on the formation of cortical folds, we will now discuss the opposite – the effect of cortical folding on the OSVZ. Although the OSVZ has a constant thickness in our model throughout the entire domain in the intermediate stage of human brain development (before cortical folds emerge), as demonstrated in red in Figure 12, left, the simulations show that this quickly changes after the first folds start to emerge. The thickness of the OSVZ starts to show spatial variations and becomes thicker beneath gyri and thinner beneath sulci, as shown in Figure 12, right. This result is consistent with what was previously observed experimentally in the human brain (Kostović et al., 2002). While it is to date not clear whether this phenomenon is rather the cause or the result of cortical folding, our study clearly indicates that the (mechanics-driven) process of cortical folding is sufficient to induce OSVZ variations. The forces generated by the mechanical growth problem not only fold the cortical layer but also lead to undulations in the deeper zones. Still, the deeper proliferating zones (i.e., inner subventricular zone and VZ) remain equally smooth as the ventricular surface.

Figure 12

Download asset Open asset

Are cortical folds affected by regional proliferation variations in the OSVZ?

Previous studies have emphasized the existence of variations in the ORGC proliferation rate in different brain regions (Hansen et al., 2010). Some of those have suggested that the proliferation rate is higher beneath gyri than beneath sulci (Borrell, 2018). Similarly, in the ferret brain, where a region close in structure to the primate’s OSVZ was found, this region shows a unique mosaic-like structure (Fietz et al., 2010; Reillo and Borrell, 2012). In this section, we aim to assess the effect of regional proliferation variations in the OSVZ on the emerging cortical folding pattern. We discuss two different heterogeneous patterns here, but have included more variations online through our user interface on GitHub, as described in the Data availability section. In the first case, the OSVZ division rate gradually decreases along the circumferential direction. In the second case, the division rate varies in a more random pattern. Figures 13 and 14 show how cortical folds develop in both cases for the varying cortical stiffness case, an initial division rate in the VZ $G_{\text{vz}}=120$, and an initial division rate in the OSVZ $G_{\text{osvz}}=20$. As expected, the evolving folding patterns slightly differ. In both cases, the first folds appear where the cell proliferation rate is highest. As expected, those regions also show a higher cell density in the cortex than regions nearby. However, both cases lead to final patterns with similar distances between sulci and folding complexity (one period doubling pattern). In addition, gyri and sulci are distributed equally – regardless of the division rate. Therefore, we may conclude that inhomogeneous cell proliferation in the OSVZ controls the location of first gyri and sulci but does not necessarily affect the distance between sulci (also referred to as folding wavelength) and the overall complexity of the emerging folding pattern. This agrees well with our previous finding that the characteristic wavelength of folding remains relatively stable for inhomogeneous cortical growth patterns (Budday and Steinmann, 2018). The simulation results are also consistent with the previously found remarkable surface expansion above the regions with higher proliferation in the OSVZ (Llinares-Benadero and Borrell, 2019).

Figure 13

Download asset Open asset

Figure 14

Download asset Open asset

In this work, we have made use of a computational model for brain growth to provide insights into the role of the OSVZ – a unique additional proliferating zone in humans – during cortical folding in the developing brain. By using computational tools, we have addressed different open questions with the aim of valuably supplementing classical experimental approaches. The advantage of the computational model is that we can, on the one hand, clearly isolate purely physical mechanisms. On the other hand, it is possible to systematically study and understand the effect of individual parameters, which is much more difficult through in vitro and in vivo experiments. Our simulations have demonstrated how the proliferation in the OSVZ enhances the complexity of cortical folding patterns in 2D and 3D. Our results show that the existence of the OSVZ particularly triggers the emergence of secondary mechanical instabilities leading to more complex folding patterns. Furthermore, the proliferation of ORGCs reduces the time required to induce the mechanical instability and thus cortical folding. Interestingly, our simulation results suggest that the generated mechanical forces not only ‘fold’ the cortex but also deeper subcortical zones including the OSVZ, which becomes thicker beneath gyri and thinner beneath sulci as a result of cortical folding. Consequently, our analyses suggest that the purely mechanics-driven process of cortical folding is sufficient to induce regional differences in the thickness of the OSVZ. Finally, our simulations reveal that inhomogeneous cell proliferation patterns in the OSVZ can control the location of first gyri and sulci but do not necessarily affect the distance between sulci and the overall complexity of the emerging folding pattern.

In conclusion, our physics-based computational modeling approach has allowed us to systematically assess the role of the OSVZ during human brain development and its effect on cortical folding – the classical hallmark of the human cortex at the organ scale. In the future, the computational framework can be used to not only better understand physiological brain development but also pathological processes – especially those involving abnormal cortical folding patterns. The computational model is able to shed new light on the interplay between the multiple processes at different scales and can help identify the main controlling parameters. However, it can only complement, not substitute, sophisticated experimental approaches that are still needed to answer questions, e.g., regarding the functional difference between progenitor cell types and corresponding lineage decisions.

The datasets and code generated and/or analyzed during the current study are available on GitHub with a user interface to change model parameters and run simulations: https://github.com/SaeedZarzor/BFSimulator (copy archived at swh:1:rev:5fd9298c69a446b5ff21fff2ecaff5ef8a0ac085).

1. 1. Betizeau M
   2. Cortay V
   3. Patti D
   4. Pfister S
   5. Gautier E
   6. Bellemin-Ménard A
   7. Afanassieff M
   8. Huissoud C
   9. Douglas RJ
   10. Kennedy H
   11. Dehay C

   (2013) Precursor diversity and complexity of lineage relationships in the outer subventricular zone of the primate

   Neuron 80:442–457.

   https://doi.org/10.1016/j.neuron.2013.09.032

   • PubMed
   • Google Scholar
2. 1. Blows WT

   (2003)

   Child brain development

   Nursing Times 99:28–31.

   • PubMed
   • Google Scholar
3. 1. Blumcke I
   2. Budday S
   3. Poduri A
   4. Lal D
   5. Kobow K
   6. Baulac S

   (2021) Neocortical development and epilepsy: insights from focal cortical dysplasia and brain tumours

   The Lancet. Neurology 20:943–955.

   https://doi.org/10.1016/S1474-4422(21)00265-9

   • PubMed
   • Google Scholar
4. 1. Borrell V
   2. Götz M

   (2014) Role of radial glial cells in cerebral cortex folding

   Current Opinion in Neurobiology 27:39–46.

   https://doi.org/10.1016/j.conb.2014.02.007

   • PubMed
   • Google Scholar
5. 1. Borrell V

   (2018) How cells fold the cerebral cortex

   The Journal of Neuroscience 38:776–783.

   https://doi.org/10.1523/JNEUROSCI.1106-17.2017

   • PubMed
   • Google Scholar
6. 1. Budday S
   2. Steinmann P
   3. Kuhl E

   (2014) The role of mechanics during brain development

   Journal of the Mechanics and Physics of Solids 72:75–92.

   https://doi.org/10.1016/j.jmps.2014.07.010

   • PubMed
   • Google Scholar
7. 1. Budday S
   2. Kuhl E
   3. Hutchinson JW

   (2015a) Period-doubling and period-tripling in growing bilayered systems

   Philosophical Magazine 95:3208–3224.

   https://doi.org/10.1080/14786435.2015.1014443

   • PubMed
   • Google Scholar
8. 1. Budday S
   2. Steinmann P
   3. Kuhl E

   (2015b) Physical biology of human brain development

   Frontiers in Cellular Neuroscience 9:257.

   https://doi.org/10.3389/fncel.2015.00257

   • PubMed
   • Google Scholar
9. 1. Budday S
   2. Steinmann P
   3. Kuhl E

   (2015c) Secondary instabilities modulate cortical complexity in the mammalian brain

   Philosophical Magazine 95:3244–3256.

   https://doi.org/10.1080/14786435.2015.1024184

   • PubMed
   • Google Scholar
10. 1. Budday S
    2. Steinmann P

    (2018) On the influence of inhomogeneous stiffness and growth on mechanical instabilities in the developing brain

    International Journal of Solids and Structures 132–133:31–41.

    https://doi.org/10.1016/j.ijsolstr.2017.08.010

    • Google Scholar
11. 1. Budday S
    2. Ovaert TC
    3. Holzapfel GA
    4. Steinmann P
    5. Kuhl E

    (2020) Fifty shades of brain: a review on the mechanical testing and modeling of brain tissue

    Archives of Computational Methods in Engineering 27:1187–1230.

    https://doi.org/10.1007/s11831-019-09352-w

    • Google Scholar
12. 1. Bystron I
    2. Blakemore C
    3. Rakic P

    (2008) Development of the human cerebral cortex: Boulder Committee revisited

    Nature Reviews. Neuroscience 9:110–122.

    https://doi.org/10.1038/nrn2252

    • PubMed
    • Google Scholar
13. 1. Darayi M
    2. Hoffman ME
    3. Sayut J
    4. Wang S
    5. Demirci N
    6. Consolini J
    7. Holland MA

    (2022) Computational models of cortical folding: a review of common approaches

    Journal of Biomechanics 139:110851.

    https://doi.org/10.1016/j.jbiomech.2021.110851

    • PubMed
    • Google Scholar
14. 1. de Rooij R
    2. Kuhl E

    (2018) A physical multifield model predicts the development of volume and structure in the human brain

    Journal of the Mechanics and Physics of Solids 112:563–576.

    https://doi.org/10.1016/j.jmps.2017.12.011

    • Google Scholar
15. 1. Fietz SA
    2. Kelava I
    3. Vogt J
    4. Wilsch-Bräuninger M
    5. Stenzel D
    6. Fish JL
    7. Corbeil D
    8. Riehn A
    9. Distler W
    10. Nitsch R
    11. Huttner WB

    (2010) OSVZ progenitors of human and ferret neocortex are epithelial-like and expand by integrin signaling

    Nature Neuroscience 13:690–699.

    https://doi.org/10.1038/nn.2553

    • PubMed
    • Google Scholar
16. 1. Fish JL
    2. Dehay C
    3. Kennedy H
    4. Huttner WB

    (2008) Making bigger brains-the evolution of neural-progenitor-cell division

    Journal of Cell Science 121:2783–2793.

    https://doi.org/10.1242/jcs.023465

    • PubMed
    • Google Scholar
17. 1. Florio M
    2. Borrell V
    3. Huttner WB

    (2017) Human-Specific genomic signatures of neocortical expansion

    Current Opinion in Neurobiology 42:33–44.

    https://doi.org/10.1016/j.conb.2016.11.004

    • PubMed
    • Google Scholar
18. 1. Gilmore EC
    2. Herrup K

    (1997) Cortical development: layers of complexity

    Current Biology 7:R231–R234.

    https://doi.org/10.1016/s0960-9822(06)00108-4

    • PubMed
    • Google Scholar
19. 1. Göktepe S
    2. Abilez OJ
    3. Kuhl E

    (2010) A generic approach towards finite growth with examples of athlete’s heart, cardiac dilation, and cardiac wall thickening

    Journal of the Mechanics and Physics of Solids 58:1661–1680.

    https://doi.org/10.1016/j.jmps.2010.07.003

    • Google Scholar
20. 1. Guerrini R
    2. Dobyns WB
    3. Barkovich AJ

    (2008) Abnormal development of the human cerebral cortex: genetics, functional consequences and treatment options

    Trends in Neurosciences 31:154–162.

    https://doi.org/10.1016/j.tins.2007.12.004

    • PubMed
    • Google Scholar
21. 1. Habas PA
    2. Scott JA
    3. Roosta A
    4. Rajagopalan V
    5. Kim K
    6. Rousseau F
    7. Barkovich AJ
    8. Glenn OA
    9. Studholme C

    (2012) Early folding patterns and asymmetries of the normal human brain detected from in utero MRI

    Cerebral Cortex 22:13–25.

    https://doi.org/10.1093/cercor/bhr053

    • PubMed
    • Google Scholar
22. 1. Hansen DV
    2. Lui JH
    3. Parker PRL
    4. Kriegstein AR

    (2010) Neurogenic radial glia in the outer subventricular zone of human neocortex

    Nature 464:554–561.

    https://doi.org/10.1038/nature08845

    • PubMed
    • Google Scholar
23. 1. Herculano-Houzel S

    (2009) The human brain in numbers: a linearly scaled-up primate brain

    Frontiers in Human Neuroscience 3:31.

    https://doi.org/10.3389/neuro.09.031.2009

    • PubMed
    • Google Scholar
24. 1. Huttner WB
    2. Kosodo Y

    (2005) Symmetric versus asymmetric cell division during neurogenesis in the developing vertebrate central nervous system

    Current Opinion in Cell Biology 17:648–657.

    https://doi.org/10.1016/j.ceb.2005.10.005

    • PubMed
    • Google Scholar
25. 1. Kostović I
    2. Judas M
    3. Rados M
    4. Hrabac P

    (2002) Laminar organization of the human fetal cerebrum revealed by histochemical markers and magnetic resonance imaging

    Cerebral Cortex 12:536–544.

    https://doi.org/10.1093/cercor/12.5.536

    • PubMed
    • Google Scholar
26. 1. Kriegstein A
    2. Noctor S
    3. Martínez-Cerdeño V

    (2006) Patterns of neural stem and progenitor cell division may underlie evolutionary cortical expansion

    Nature Reviews. Neuroscience 7:883–890.

    https://doi.org/10.1038/nrn2008

    • PubMed
    • Google Scholar
27. 1. Libé-Philippot B
    2. Vanderhaeghen P

    (2021) Cellular and molecular mechanisms linking human cortical development and evolution

    Annual Review of Genetics 55:555–581.

    https://doi.org/10.1146/annurev-genet-071719-020705

    • PubMed
    • Google Scholar
28. 1. Llinares-Benadero C
    2. Borrell V

    (2019) Deconstructing cortical folding: genetic, cellular and mechanical determinants

    Nature Reviews. Neuroscience 20:161–176.

    https://doi.org/10.1038/s41583-018-0112-2

    • PubMed
    • Google Scholar
29. 1. Lui JH
    2. Hansen DV
    3. Kriegstein AR

    (2011) Development and evolution of the human neocortex

    Cell 146:18–36.

    https://doi.org/10.1016/j.cell.2011.06.030

    • PubMed
    • Google Scholar
30. 1. Lukaszewicz A
    2. Savatier P
    3. Cortay V
    4. Giroud P
    5. Huissoud C
    6. Berland M
    7. Kennedy H
    8. Dehay C

    (2005) G1 phase regulation, area-specific cell cycle control, and cytoarchitectonics in the primate cortex

    Neuron 47:353–364.

    https://doi.org/10.1016/j.neuron.2005.06.032

    • PubMed
    • Google Scholar
31. 1. Noctor SC
    2. Flint AC
    3. Weissman TA
    4. Wong WS
    5. Clinton BK
    6. Kriegstein AR

    (2002) Dividing precursor cells of the embryonic cortical ventricular zone have morphological and molecular characteristics of radial glia

    The Journal of Neuroscience 22:3161–3173.

    https://doi.org/10.1523/JNEUROSCI.22-08-03161.2002

    • PubMed
    • Google Scholar
32. 1. Noctor SC
    2. Martínez-Cerdeño V
    3. Ivic L
    4. Kriegstein AR

    (2004) Cortical neurons arise in symmetric and asymmetric division zones and migrate through specific phases

    Nature Neuroscience 7:136–144.

    https://doi.org/10.1038/nn1172

    • PubMed
    • Google Scholar
33. 1. Noctor SC
    2. Martínez-Cerdeño V
    3. Kriegstein AR

    (2007) Contribution of intermediate progenitor cells to cortical histogenesis

    Archives of Neurology 64:639–642.

    https://doi.org/10.1001/archneur.64.5.639

    • PubMed
    • Google Scholar
34. 1. Nonaka-Kinoshita M
    2. Reillo I
    3. Artegiani B
    4. Martínez-Martínez MÁ
    5. Nelson M
    6. Borrell V
    7. Calegari F

    (2013) Regulation of cerebral cortex size and folding by expansion of basal progenitors

    The EMBO Journal 32:1817–1828.

    https://doi.org/10.1038/emboj.2013.96

    • PubMed
    • Google Scholar
35. 1. Nowakowski TJ
    2. Pollen AA
    3. Sandoval-Espinosa C
    4. Kriegstein AR

    (2016) Transformation of the radial glia scaffold demarcates two stages of human cerebral cortex development

    Neuron 91:1219–1227.

    https://doi.org/10.1016/j.neuron.2016.09.005

    • PubMed
    • Google Scholar
36. 1. Pebworth MP
    2. Ross J
    3. Andrews M
    4. Bhaduri A
    5. Kriegstein AR

    (2021) Human intermediate progenitor diversity during cortical development

    PNAS 118:e2019415118.

    https://doi.org/10.1073/pnas.2019415118

    • PubMed
    • Google Scholar
37. 1. Poluch S
    2. Juliano SL

    (2015) Fine-tuning of neurogenesis is essential for the evolutionary expansion of the cerebral cortex

    Cerebral Cortex 25:346–364.

    https://doi.org/10.1093/cercor/bht232

    • PubMed
    • Google Scholar
38. 1. Rakic P

    (1988) Specification of cerebral cortical areas

    Science 241:170–176.

    https://doi.org/10.1126/science.3291116

    • PubMed
    • Google Scholar
39. 1. Razavi MJ
    2. Zhang T
    3. Li X
    4. Liu T
    5. Wang X

    (2015) Role of mechanical factors in cortical folding development

    Physical Review. E, Statistical, Nonlinear, and Soft Matter Physics 92:032701.

    https://doi.org/10.1103/PhysRevE.92.032701

    • PubMed
    • Google Scholar
40. 1. Reillo I
    2. de Juan Romero C
    3. García-Cabezas MÁ
    4. Borrell V

    (2011) A role for intermediate radial glia in the tangential expansion of the mammalian cerebral cortex

    Cerebral Cortex 21:1674–1694.

    https://doi.org/10.1093/cercor/bhq238

    • PubMed
    • Google Scholar
41. 1. Reillo I
    2. Borrell V

    (2012) Germinal zones in the developing cerebral cortex of ferret: ontogeny, cell cycle kinetics, and diversity of progenitors

    Cerebral Cortex 22:2039–2054.

    https://doi.org/10.1093/cercor/bhr284

    • PubMed
    • Google Scholar
42. 1. Rodriguez EK
    2. Hoger A
    3. McCulloch AD

    (1994) Stress-Dependent finite growth in soft elastic tissues

    Journal of Biomechanics 27:455–467.

    https://doi.org/10.1016/0021-9290(94)90021-3

    • PubMed
    • Google Scholar
43. 1. Schmechel DE
    2. Rakic P

    (1979) A Golgi study of radial glial cells in developing monkey telencephalon: morphogenesis and transformation into astrocytes

    Anatomy and Embryology 156:115–152.

    https://doi.org/10.1007/BF00300010

    • PubMed
    • Google Scholar
44. 1. Shinmyo Y
    2. Terashita Y
    3. Dinh Duong TA
    4. Horiike T
    5. Kawasumi M
    6. Hosomichi K
    7. Tajima A
    8. Kawasaki H

    (2017) Folding of the cerebral cortex requires Cdk5 in upper-layer neurons in gyrencephalic mammals

    Cell Reports 20:2131–2143.

    https://doi.org/10.1016/j.celrep.2017.08.024

    • PubMed
    • Google Scholar
45. 1. Smart IHM
    2. Dehay C
    3. Giroud P
    4. Berland M
    5. Kennedy H

    (2002) Unique morphological features of the proliferative zones and postmitotic compartments of the neural epithelium giving rise to striate and extrastriate cortex in the monkey

    Cerebral Cortex 12:37–53.

    https://doi.org/10.1093/cercor/12.1.37

    • PubMed
    • Google Scholar
46. 1. Takahashi E
    2. Folkerth RD
    3. Galaburda AM
    4. Grant PE

    (2012) Emerging cerebral connectivity in the human fetal brain: an mr tractography study

    Cerebral Cortex 22:455–464.

    https://doi.org/10.1093/cercor/bhr126

    • PubMed
    • Google Scholar
47. 1. Tallinen T
    2. Chung JY
    3. Biggins JS
    4. Mahadevan L

    (2014) Gyrification from constrained cortical expansion

    PNAS 111:12667–12672.

    https://doi.org/10.1073/pnas.1406015111

    • PubMed
    • Google Scholar
48. 1. Wang S
    2. Saito K
    3. Kawasaki H
    4. Holland MA

    (2022) Orchestrated neuronal migration and cortical folding: a computational and experimental study

    PLOS Computational Biology 18:e1010190.

    https://doi.org/10.1371/journal.pcbi.1010190

    • PubMed
    • Google Scholar
49. 1. Zarzor MS
    2. Kaessmair S
    3. Steinmann P
    4. Blümcke I
    5. Budday S

    (2021) A two-field computational model couples cellular brain development with cortical folding

    Brain Multiphysics 2:100025.

    https://doi.org/10.1016/j.brain.2021.100025

    • Google Scholar

Author details

1. Mohammad Saeed Zarzor

   Friedrich-Alexander-Universität Erlangen-Nürnberg, Institute of Applied Mechanics, Erlangen, Germany

   Contribution

   Conceptualization, Software, Validation, Investigation, Visualization, Methodology, Writing - original draft

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-3005-6115

2. Ingmar Blumcke

   University Hospitals Erlangen, Institute of Neuropathology, Erlangen, Germany

   Contribution

   Resources, Validation, Writing - review and editing

   Competing interests

   No competing interests declared

3. Silvia Budday

   Friedrich-Alexander-Universität Erlangen-Nürnberg, Institute of Applied Mechanics, Erlangen, Germany

   Contribution

   Conceptualization, Supervision, Funding acquisition, Methodology, Writing - original draft, Project administration, Writing - review and editing

   For correspondence

   silvia.budday@fau.de

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-7072-8174

• 1,382

  views

• 139

  downloads

• 9

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.