Cerebellar folding is initiated by mechanical constraints on a fluid-like layer without a cellular pre-pattern

Abstract
eLife digest
Introduction
Results
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

Models based in differential expansion of elastic material, axonal constraints, directed growth, or multi-phasic combinations have been proposed to explain brain folding. However, the cellular and physical processes present during folding have not been defined. We used the murine cerebellum to challenge folding models with in vivo data. We show that at folding initiation differential expansion is created by the outer layer of proliferating progenitors expanding faster than the core. However, the stiffness differential, compressive forces, and emergent thickness variations required by elastic material models are not present. We find that folding occurs without an obvious cellular pre-pattern, that the outer layer expansion is uniform and fluid-like, and that the cerebellum is under radial and circumferential constraints. Lastly, we find that a multi-phase model incorporating differential expansion of a fluid outer layer and radial and circumferential constraints approximates the in vivo shape evolution observed during initiation of cerebellar folding.

https://doi.org/10.7554/eLife.45019.001

eLife digest

The human brain has a characteristic pattern of ridges and grooves that make up its folded shape. Folds in the outer layer of the brain, known as the cortex, increase the surface area and make more space for cells to connect and form complex circuitries. Different models have been put forward to explain how these folds form during development. Examples include tension from cells pulling areas of the cortex together, or layers of the cortex growing at different rates, causing the cortex to buckle and create folds. Discriminating between these different models requires biological information about the cells and tissue of the brain at the start of the folding process. However, it has been difficult to extract this information when considering the development of the human brain in three dimensions.

Lawton et al. have overcome these difficulties by using a part of the mouse brain called the cerebellum as a simpler system. As in humans, the mouse cerebellum is a densely folded structure, sitting underneath the brain, that plays a major role in regulating movements, as well as cognition. The symmetrical structure of the mouse cerebellum means it can be analyzed in two dimensions, making it easier to track the mechanics of folding.

By applying the extracted biological data onto a mathematical model, Lawton et al. showed folding was driven by a combination of previously unknown features. For instance, that cells in the outer layer of the cerebellum grow faster than cells in the center, with cells growing uniformly across the outer layer. Other features include the fluid-like composition of the outer layer, which allows cells to move freely and regularly change position, and tensions surrounding the cerebellum mechanically straining its growth. Notably, the pattern of cells and tissue fibers in the cortex had no influence over these mechanical properties and provided no pre-indication of where the sites of folding would occur. The data collected deviates from other models, and has led Lawton et al. to propose a new explanation for how the brain folds, incorporating these newly found features.

Problems with brain folding during human development can lead to debilitating conditions. Applying this new model to folding disorders of the human brain could help scientists to understand how these folding defects arise.

https://doi.org/10.7554/eLife.45019.002

Introduction

Recent work to elucidate the mechanics of brain folding has primarily focused on the human cerebral cortex and involved models of directed growth, axonal tension, or differential expansion of elastic materials that generate compressive forces to drive mechanical instabilities leading to folding (Tallinen et al., 2014; Ronan et al., 2014; Bayly et al., 2013; Xu et al., 2010; Hohlfeld and Mahadevan, 2011; Bayly et al., 2014; Lejeune et al., 2016; Karzbrun et al., 2018). Current elastic material models are able to create three-dimensional shapes strikingly similar to the final folds seen in the adult human cortex (Tallinen et al., 2016). A recent multi-phase model (Engstrom et al., 2018) that includes elastic and fluid-like layers, differential expansion and radial constraints takes into consideration that multiple factors could lead to folding in the developing brain. However, the cell and tissue level mechanics actually present at the initiation of folding have not been considered or defined, as technological limitations are significant in animals with a folded cerebrum.

The murine cerebellum has a simple alignment of 8–10 stereotypical folds along the anterior-posterior axis. Combined with the genetic tools available in mouse this allows for precise developmental interrogation to identify and analyze the in vivo cellular and tissue level behaviors driving growth and folding. The developing cerebellum is distinct from the cerebral cortex, as it has a temporary external granule cell layer (EGL) of proliferating granule cell precursors that cover the surface and generate growth primarily in the anterior-posterior (AP) direction (Leto et al., 2016; Legué et al., 2015; Legué et al., 2016). During development a thickening occurs in the EGL at the base of each forming fissure, termed anchoring center (AC) (Sudarov and Joyner, 2007), whereas in the adult cerebellum the inner granule cell layer (IGL), generated by the EGL during the first two weeks of life, is thinnest at the ACs. Previous work on cerebellar folding utilized a tri-layer elastic model incorporating the EGL, the adjacent molecular layer, rich in axons and dendrites, and the IGL (Lejeune et al., 2016). However, neither the molecular layer nor the IGL are present when folding is initiated in the embryo. Therefore we argue that a bilayer system consisting of the EGL and underlying core, is a more appropriate approximation for cerebellar folding.

Here we show that cerebellar folding emerges from differential expansion between an un-patterned, rapidly expanding EGL and an underlying core. Additionally, we demonstrate that the measured stiffness differential, compressive forces, and the thickness variation in the EGL are all inconsistent with traditional elastic wrinkling models driven by differential growth. Furthermore, we demonstrate that the expansion of the EGL is uniform, and fluid-like, and that the cerebellum is under radial and circumferential constraints when folding initiates. Lastly, we constrain the recent multi-phase model with our in vivo data and find we can capture the temporal shape evolution seen during mouse cerebellum folding initiation. The implications of our findings for human cerebral cortex folding are discussed.

Results

Tissue level mechanics drive folding

It is well known that differentially expanding bilayer systems can wrinkle to relax building stress (Richman et al., 1975; Nelson, 2016; Hannezo et al., 2012; Shyer et al., 2013; Wiggs et al., 1997). We reasoned that in the cerebellum the EGL could behave as a quickly expanding outer layer and its attachment to a more slowly growing core could generate forces that result in a wrinkling-like phenotype. To test whether the cerebellum has differential expansion between the two layers, we measured the expansion of the EGL and the core during the time of initiation of folding from midline sagittal sections (Figure 1a–d). Unlike the cerebral cortex, the unfolded murine cerebellum is a simple cylinder-like structure elongated in the medio-lateral axis (Figure 1e) (Szulc et al., 2015). All folds in the medial cerebellum (vermis) are aligned in the same axis allowing 2-D measurements to estimate expansion in the anterior-posterior axis of the vermis. Therefore the length of the surface of the EGL was used as a measure of the cerebellum surface area and the area of the core as an approximation of cerebellum volume (Figure 1d), and measurements were made each day from embryonic day 16.5 (E16.5) through postnatal day 0 (P0). In cross-section the unfolded cerebellum approximates a semicircle, therefore we reasoned that if the cerebellum were to remain unfolded then the ratio of expansion between the length of the EGL and the area of the core should approximate the ratio of the circumference of a semi-circle to its area. Of significance, we found that at E16.5 and E17.5 the ratios of growth between the EGL and core closely approximated the expansion of a semi-circle. However, at E18.5 and P0 the expansion rate of the EGL was greater than the rate of core expansion (Figure 1f). Thus we uncovered that the cerebellum does indeed go through a phase of differential expansion. We next determined whether differential expansion correlates with when folding occurs by calculating a folding index (the convex curvature of the EGL divided by the length of the EGL) at each stage (Mota and Herculano-Houzel, 2015). Indeed, we found that the cerebellum remains unfolded during the initial proportional expansion between the EGL and core and only folds when the differential expansion is initiated (Figure 1g). These results provide quantitative evidence that cerebellar folding involves tissue level mechanical forces arising from differential expansion.

Figure 1

Download asset Open asset

Initiation of cerebellar folding correlates with timing of differential expansion.

(**A–D**), H and E stained midline sagittal sections of FVB/N mice at the indicated embryonic (E) and postnatal days (P). Anterior to the left. Stars: ACs. Brackets: anterior/posterior ends of the EGL. Black line and red line in (D): EGL and convex length, respectively. Shaded area: core. (E) manganese enhanced magnetic resonance imaging of P11 cerebellum outlined in white adapted from Szulc et al. (2015). Anterior to the top. Vermis indicated by dotted yellow lines. (F), At E16.5 and E17.5 expansion of EGL length and cerebellar area fit the proportional expansion of a semi-circle (curve). At E18.5 and P0 EGL expansion is greater than core area growth creating differential expansion. (G), Folding index [1 - (convex length/EGL length) x 100] reveals folding initiates during differential expansion. Scale bars: 200 μm.

https://doi.org/10.7554/eLife.45019.003

In vivo data contradict elastic bilayer models

Since there is differential expansion between the EGL and the core and as this type of expansion is the driver of elastic bilayer models we tested whether the properties of cerebellar tissue are consistent with the requirements and predictions of such models. Briefly, the initial resulting wrinkling instability defines the distance between folds as the initial sinusoidal undulations increase in amplitude to ultimately turn into lobules. The folding wavelength depends on the thickness of the external layer (EGL) and the ratio of the stiffness of the two layers (EGL/core). In particular, for a planar geometry, with the stiffness of the external layer defined as E_o, the stiffness of the core as E_i, and the thickness of the external layer denoted as t, the folding wavelength λ is given by Allen (1969)

λ = 2 π t {(\frac{1}{3} \frac{E_{o}}{E_{i}})}^{1 / 3} .

If the length of the system is $l$ , then the number of folds is inversely proportional to the thickness of the EGL

n = \frac{l}{λ} \propto \frac{l}{t} {(\frac{E_{i}}{E_{o}})}^{1 / 3} .

We explored a standard elastic bilayer model in a circular geometry using the observed ratio of thickness of the EGL to radius of the cerebellum near the onset of shape change (E16.5) and invoked a neo-Hookean elastic solid for both layers (Zhao and Zhao, 2017). The resulting shape change was studied as a function of the ratio of the layer stiffness values (Figure 2a). We found that to produce the observed number of folds (three in the semi-circular cerebellum and six in the circular model) at initiation of folding through wrinkling based models constrained by our measurements of the embryonic cerebellum, a large stiffness ratio was required of around 50. To map the stiffness contrast in the cerebellum we used scanning acoustic microscopy (SAM) to measure the bulk modulus of the cerebellum daily from E16.5 to P18.5 (Figure 2b–c, Figure 2—figure supplement 1) using established methods (Rohrbach et al., 2015; Rohrbach et al., 2018; Rohrbach and Mamou, 2018). For small deformations, the instantaneous bulk modulus should linearly relate to the stiffness and, therefore, the ratio of the instantaneous bulk moduli should scale similarly to the ratio of stiffnesses (assuming the same Poisson’s ratio for the EGL and for the core, neither of which have been directly measured). While this qualitative approach and SAM tissue preparation protocols may not be able to produce the absolute values of the elastic properties of the tissues, it can give a reasonable indication of the relative stiffnesses of different parts in the cerebellum. Using this estimation, we found that the EGL has a slightly higher instantaneous bulk modulus than the core at all stages measured. Unsurprisingly, the ratio (~1.05:1) was not close to being sufficient to produce a folding wavelength similar to that in the cerebellum (Figure 2d). Consistent with our finding, small modulus contrasts have been reported for other brain regions with multiple loading modes, such as shear, compression, and tension (Xu et al., 2010; Lejeune et al., 2016; Budday et al., 2017). Elastic material models with graded growth profiles have been developed that predict folding of cerebral cortex without a large stiffness differential (Tallinen et al., 2014). However, these models are still bound by other measurable requirements as discussed below.

Figure 2 with 3 supplements see all

Download asset Open asset

Measured tissue stiffness, stress, and shape at folding initiation are inconsistent with wrinkling models.

(A), Inducing the correct number of folds through a wrinkling model requires a stiffness differential between the layers of 50 fold (μf/μs = 50, g = 1.05). (**B,C**) Acoustic mapping of cerebellar slices show a slightly stiffer EGL than core at each stage (anova df = 5; p=1.0e⁻⁴, F = 13.59), but not the required differential. Stars indicate statistical differences. (D) Wrinkling simulations constrained by developmental data produce wavelengths inconsistent with the embryonic mouse cerebellum. (E) Elastic simulations predict the EGL remains closed after cutting. (**F,G**) Images of a live cerebellar slice before and after cutting, and images from time lapse movie, (H) show the EGL opens, revealing circumferential tension along the EGL. Red and yellow dots: cut edges. Lines: relaxation angle. (I) Staining of nuclei with EGL outlined and lines used to measure thickness. (**J–L**) Normalized EGL thickness (thickness/mean thickness) at the ACs increases during folding initiation (anova E16.5 df = 29, p=8.2e⁻²⁰, F = 12.59. E17.5 df = 29, p=3.4e⁻¹¹⁶, F = 62.78, E18.5 df = 57, p=6.8e⁻⁶⁷, F = 13.28). At E16.5 only brains with visible ACs were included. Error bars: S.D.

https://doi.org/10.7554/eLife.45019.004

Elastic bi-layer wrinkling models predict compressive forces in the outer layer. Simulations performed of cuts through the outer layer and into the inner layer predict that upon relaxation the outer layer should not open (Figure 2e). We tested whether this prediction reflects the biology using surgical dissection blades to make radial cuts across the meninges, EGL, and into the core of live E16.5 tissue slices. Time-lapse imaging revealed that, in contrast to the prediction, the EGL opens as well as part of the underlying cut in the core (Figure 2f–h, Figure 2—figure supplement 2a–c, and Video 1). This result indicates there is circumferential tension within the outer layers of the cerebellum. This finding also rules out the elastic models with graded growth profiles as they predict compressive forces in the outer region as well.

Video 1

Download asset

posterframe for video — Live slice cutting and relaxation reveals circumferential tension along the EGL.

Time-lapse video shows relaxation of live tissue slice after cutting radially through the EGL and into the underlying core. Images were acquired every 10 s for 10 min. The time-lapse was started moments after the tissue was collected in frame after the cut. The slice shown in the video is the same as in Figure 2f–h.

https://doi.org/10.7554/eLife.45019.008

The elastic bi-layer model requires the EGL to be thinnest at the base of each AC, which are the lowest parts of the cerebellar surface. Thus, the EGL should have an ‘in-phase’ thickness variation. Without this feature, a purely elastic model – bi-layer based or even graded growth profile based – cannot be in mechanical equilibrium (in the quasistatic limit) (Engstrom et al., 2018). However, we previously reported that the embryonic EGL is thickest in the ACs when folding initiates, that is it has an ‘out-of-phase’ thickness variation (Sudarov and Joyner, 2007). To validate this observation, we quantified the thickness variations in the EGL centered at the ACs present at E16.5–18.5. Not all cerebella have visible AC at E16.5. However in the subset that do and in the three ACs present at E17.5, the EGL was found to be 1.2–1.4 times thicker in the ACs than in the surrounding EGL (Figure 2i–l and Figure 2—figure supplement 3). Moreover, the thickness ratio increased to 1.7 times at E18.5 (Figure 2l). As described above, the final thickness variations of the IGL (as well as the molecular layer) of the cerebellar cortex are in-phase, just as the layers of the adult cerebral cortex. These results further show that traditional elastic wrinkling models cannot capture the initiation of cerebellum folding, and highlight the importance of making biological measurements at the time of folding rather than when it is complete.

Uniform outer layer expansion without a cellular pre-pattern

As elastic bi-layer models do not align with the biology of cerebellar folding, we looked for other drivers of morphometric changes. Since the EGL drives the majority of cerebellar growth (Leto et al., 2016; Legué et al., 2015; Legué et al., 2016), we first tested whether regional differences in EGL proliferation rates are present that could influence the folding pattern of the cerebellum. Proliferation rates (S phase index) were measured in the EGL during folding initiation (E16.5 and E17.5) in the inbred FVB/N strain to reduce variation between samples. First we asked if the regions that will give rise to distinct sets of lobules have different rates of proliferation that could contribute to the larger and smaller sizes that the lobules ultimately attain. We focused on the anterior cerebellum that divides into a larger region with lobules 1–3 (L123) and smaller region (L45), as well as the central area that comprises lobules 6–8 (L678) of the cerebellum (Figure 3a–b). The more posterior cerebellum does not consistently fold at this stage, thus measurements were not included. Interestingly, we found that the proliferation rates were similar in the three regions at E16.5 (Figure 3c). The EGL proliferation rate at E17.5 in L678 was slightly reduced compared to the L123 region, but no other differences were found (Figure 3d). Thus proliferation is uniform just before initiation of folding and the small difference found during folding does not correlate with lobule size. This result indicates that lobule size is not determined by modulating the levels of proliferation at the onset of folding. Rather, lobule size could be set by both the timing of invagination, and the distance between ACs as granule cell precursors in one lobule do not cross the surrounding ACs to contribute to an adjacent lobule (Legué et al., 2015).

Figure 3 with 1 supplement see all

Download asset Open asset

At folding initiation the EGL has uniform proliferation, cell size, and cell shape.

(**A,B**) Low and high power images of immuno-histochemical (IHC) staining of sagittal cerebellar sections to measure proliferation in the lobules (L) indicated at 25 μm windows surrounding the ACs (stars). Scale bar: 200 μm. (**C,D**) EGL proliferation rates are shown before and during the onset of invagination (Two-way anova: df = 2. (C) p=0.10, F = 4.36 (D) p=0.03, F = 10.31). (**E,F**) Proliferation rates are shown in the AC and in the surrounding EGL showing uniformity (Two-way anova: df = 18. (E) p=0.03, F = 2.15 (F) p=2.1e⁻³, F = 3.06). (G) Section of Atoh1-CreER/+; R26^MTMG/+ E16.5 cerebellum showing masked labeled cells. Scale bar: 20 μm. (**H,I**) Cell shape (sphericity) measurements before and during folding (anova df = 3. (H) p=0.34, F = 1.13 i p=0.61, F = 0.61). (**J,K**) Cell size measurements before and during folding (anova df = 3. (J) p=3.6e⁻³, F = 4.75, (K) p=0.85, F = 0.26). Stars indicate statistical differences. Error bars: S.D.

https://doi.org/10.7554/eLife.45019.009

Each AC is first detected as a regional inward thickening of the EGL (Sudarov and Joyner, 2007) (Figure 2i–l and Figure 2—figure supplement 3). We measured the proliferation of the EGL specifically within the forming AC regions to test whether altered proliferation rates could explain the thickenings and therefore the initiation of an AC. We found the rate of proliferation within each forming AC region at E16.5 and E17.5 was the same as in the surrounding EGL (Figure 3e,f), thus proliferation within all regions of the EGL at the initiation of folding is uniform. Furthermore, regional modulation of proliferation does not form or position the ACs.

At E18.5, after the initiation of folding, we found that the rate of proliferation was significantly lower in the L678 region compared with the L123 and L45 regions (Figure 3—figure supplement 1a). However, proliferation within the ACs at E18.5 remained uniform with the surrounding regions (Figure 3—figure supplement 1b). Since ACs compartmentalize the EGL, our results show that regional differences in proliferation rates arise in lobule regions after initiation of folding, which thus could be important for determining the ultimate size of the folds.

Changes in cell size and shape have been shown to induce morphological changes (Mammoto and Ingber, 2010; Harding et al., 2014; Stemple, 2005; He et al., 2014). To test if regionally specific regulation of cell shape or size directs folding, we fluorescently labeled cell membranes of scattered granule cell precursors (GCPs) in the EGL using genetics (Atoh1-CreER/+; R26^MTMG/+ mice injected with tamoxifen two days prior to analysis). We then segmented the cells in 3D and quantified their sphericity (Figure 3g). We discovered that GCPs in the EGL take on a large variation of shapes and sizes at E16.5 and E18.5. However, we found no difference in cell shape in the different lobule regions of the EGL or between the AC areas and the surrounding EGL at each age (Figure 3h,i). Cell size was uniform at both stages except for a slight reduction in L678 at E16.5 when compared with L123 and the AC regions. However, the size of cells is reduced at E18.5 compared to E16.5 (Figure 3j,k and Figure 3—figure supplement 1c). Thus, the proliferating GCPs that drive expansion of the EGL have both uniform proliferation rates and similar shapes and sizes across the lobule regions defined by the first three ACs at folding initiation.

Uniform fiber distribution and radial tension at folding initiation

The EGL is traversed by fibers of Bergmann glial and radial glial cells (Leung and Li, 2018; Yuasa, 1996; Yamada and Watanabe, 2002). We tested whether the fibers are distributed in patterns that could locally change the physical properties of the EGL and induce invaginations. Genetics was used to fluorescently label cell membranes of scattered glial cells (nestin-creER/+;R26^MTMG/+ mice injected with tamoxifen at E14.5) (Figure 4a). Fibers crossing the EGL at E16.5 were counted in sagittal slices and aligned relative to the ACs (Figure 4b). This analysis showed that the Bergmann glial and radial glial fibers are distributed evenly along the AP axis of the EGL, and therefore are not directing the positions where folding initiates based on an uneven regional distribution.

Figure 4

Download asset Open asset

The EGL has a uniform distribution of crossing fibers at folding initiation.

(A) E16.5 nestin-CreER/+; R26^MTMG/+ cerebellum section showing fluorescent labeling of radial and Bergmann glial fibers. Stars: AC. Dotted lines denote EGL. Scale bar: 50 μm. (B) Measurements of fiber density in the ACs compared to the surrounding EGL (anova df = 10; p=0.76, F = 0.66). Error bars: S.D. (**C–E**) Still images of a tissue cutting experiment to test for radial tension between the EGL and the VZ. Red and yellow dots: cut edges at t = 0, 10. Lines: relaxation angle.

https://doi.org/10.7554/eLife.45019.011

Tension based folding models suggest constraints from axons and other fibers could direct folding (Xu et al., 2010; Van Essen, 1997). Since the cerebellum is under circumferential tension, as demonstrated above, we examined evidence of radial tension between the EGL and the ventricular zone (VZ) at the initiation of folding. Cuts were made in live E16.5 tissue slices between the EGL and VZ running approximately parallel to them so that they cut across radial fibers in the anterior cerebellum (Figure 4c). As predicted, after cutting the tissue relaxed revealing tension directed radially within the cerebellum (Figure 4d,e and Figure 2—figure supplement 2 and Video 2). Interestingly, quantification of how the radial and horizontal cuts open revealed that only the horizontal cuts opened along the full length of the cut although they opened more slowly than radial cuts (Figure 2—figure supplement 2g–j), indicating different stress profiles in the two orientations.

Video 2

Download asset

Taken together, at the time of folding initiation the EGL, which is driving the differential expansion, is itself growing uniformly and the cerebellum is under both radial and circumferential constraints. Finally, there is no evidence of any pre-patterning in the EGL in either cellular behaviors or fiber distribution.

The EGL is fluid-like as cells undergo dynamic rearrangement

As the granule cells within the EGL have such varied shapes as shown above, we looked to see if the cells within the EGL were undergoing any rearrangement movements that may indicate fluid properties. A small, scattered fraction of nuclei in the EGL were fluorescently labeled (Atoh1-CreER/+; R26^ntdTom/+ injected with tamoxifen two days prior to imaging) and ex vivo slice-culture time-lapse imaging was performed for up to five hours. Tracking the cell positions through time revealed that granule cells within the EGL are highly motile within the EGL. Furthermore, there was no obvious directionality or collectivity to the movement. However, the dynamic motility resulted in the constant exchanging of nearest neighbors over the course of tens to hundreds of minutes and shows that at the timescale of folding the EGL is more fluid-like than a solid epithelial layer (Figure 5 and Videos 3 and 4).

Figure 5

Download asset Open asset

During folding initiation, cells within the EGL are motile and undergo rearrangement.

(A) Image of E16.5 live cerebellar slice (Atoh1-creER/+R26^Ai75/+) showing scattered labeling within the EGL. (B) Red dots indicate starting position, displacement arrows show final position of marked cells after 5 hr. White dashed line indicates outer edge of EGL. (C) Still images from time-lapse, inset above. Cells tracked and marked with colored spheres exchange nearest neighbors over a time-scale of tens of minutes. Scale bars are 20 μm.

https://doi.org/10.7554/eLife.45019.013

Video 3

Download asset

Video 4

Download asset

Multi-phase wrinkling model simulates cerebellar shape change during folding initiation

We recently developed a model for folding from differentially expanding bi-layer tissues that takes into account the out-of-phase thickness of the outer layer of several systems and possible contribution of radial mechanical constraints present in neurological tissue (Engstrom et al., 2018). We applied the model here to the initiation of cerebellar folding based on five primary assumptions. First, the core is an incompressible material (μ) as indicated by the bulk modulus measurements. Second the outer layer, that is the EGL, expands uniformly (k_t) as shown by the proliferation rate. Third, the EGL is assumed to be a fluid-like material as demonstrated by the live-imaging of neighbor exchanges. Fourth, there is an elastic component radially to the entire cerebellum (k_r), seen in the cutting and relaxation experiment and possibly mediated by radial glia. Fifth, the EGL is constrained towards a uniform thickness (β), possibly by Bergmann glia fibers spanning the EGL. Given the interplay between incompressible material, compressible fibrous material, and a proliferating non-elastic EGL, this model is multi-phase.

An energy functional parameterized by both the inner and outer boundary of the EGL and incorporating the above five assumptions into three dimensionless parameters (μ/k_r, k_r/k_t, k_t/β) is minimized to yield an equation for a driven harmonic oscillator resulting in sinusoidal shapes for both the inner and outer boundary of the EGL given an initial elliptical shape. In contrast with the elastic bilayer wrinkling model, EGL thickness oscillations are found to be out-of-phase with the surface height (radius) oscillations when 0 < μ/k_r <1. Additionally, the model predicts that the ratio of the measured surface height amplitude (A_r) and the EGL thickness amplitude (A_t) is given by

\frac{A_{r}}{A_{t}} = \frac{\frac{μ}{k_{r}}}{1 - \frac{μ}{k_{r}}},

which need not be $≫ 1$ as is typical of elastic bilayer wrinkling, and the number of initial folds at E16.5 is determined by

n = \sqrt{\frac{k_{t}}{β}} \sqrt{1 + \frac{\frac{μ}{k_{t}}}{1 - \frac{μ}{k_{r}}}}

Note that in contrast with elastic wrinkling, the number of initial folds does not depend on the thickness (a length scale) of the EGL, but only on material properties.

Given that our tissue cutting and relaxation experiment revealed circumferential tension in the cerebellum at folding initiation (Figure 2f–h, Video 1), we returned to the mathematics and found a previously unrealized geometric relationship in the circular limit of the model that in fact assumes circumferential tension in addition to the previously discussed radial tension given that the perimeter of a circle is determined by its radius.

To rigorously test the shape prediction of the model, we first constrained 3 of the 5 parameters for a circular version of the model by using both the thickness amplitude, and average thickness of the EGL, as measured at E16.5, and the number of initial folds. Secondly, the parameter μ/k_r (denoted as ε) was assumed to scale linearly with time. Together, this allowed for the generation of shape predictions at later developmental stages (E17.5 and E18.5) from the E16.5 starting approximation. Solving the analytical model as constrained by our measured embryonic data we found that it closely approximates the phase and amplitude behavior of EGL thickness and radius oscillations from E16.5 through E18.5. (Figure 6a–c). However, the model is not able to produce self-contacting folds or hierarchical folding, both of which are seen in the cerebellum at later stages.

Figure 6

Download asset Open asset

A multiphase model with radial and circumferential constraints and liquid-fibrous EGL composition approximates evolution of cerebellar shape.

(A) Schematic of multiphase model showing types of tension. (B) Thickness variations that arise concomitant with folding approximate those seen in the cerebellum. (C) Shape factor analysis: model for semicircle (red), multi-phase model shape prediction (blue) and actual shape measured from sections (black). Assumed linear relationship between ε and time, ε(T) = 0.3(T −15.5).

https://doi.org/10.7554/eLife.45019.016

Hierarchical folding involves differential growth

The cerebellum has hierarchical folding in which the initial folds become subdivided. Given that ACs hold their position during development and compartmentalize granule cells within lobules of the EGL (Legué et al., 2015) we reasoned that the ACs could be acting as mechanical boundaries enabling similar mechanics to drive the secondary folding. To test this possibility we measured the expansion of the EGL and the core of the individual lobule regions from E18.5 to P3. We found that indeed in the lobule regions that undergo folding there is a temporal correlation between when the onset of sub-folding and differential expansion occur (Figure 7a–d). In contrast, the region (L45) that does not fold during the same time period has a different, more rectangular shape, and the ratio of EGL growth to core growth is proportional for a rectangle during the time measured (Figure 7).

Figure 7

Download asset Open asset

Differential expansion correlates with progressive subfolding of initial folds.

(**A,B**) H and E stained midline sagittal sections of FVB/N cerebella at P1 and P3 with three lobule regions highlighted in red (L123), green (**L45**), and blue (L678). (C) Expansion of lobule length and lobule area for L123 and L678 approximate the proportional expansion of a semi-circle (curve) in both regions at E18.5 and P0. After P0 the EGL expansion in both regions increases more than the underlying area creating differential expansion. (D) Folding initiates during regional differential expansion. (E) The expansion in length and area of L45 is proportional to a columnar shape (curve and inset figure) from E18.5 to P3. (F) the L45 region remains unfolded through P3. Scale bars: 200 μm.

https://doi.org/10.7554/eLife.45019.017

Discussion

Here we have provided experimental evidence that cerebellar folding emerges without obvious pre-patterning. Additionally, the outer layer has fluid-like properties and expands uniformly, and the growth creates a differential expansion between the outer layer and the core. Thus, traditional morphometric cellular behaviors such as changes in cell shape, size and proliferation do not direct where cerebellar folding initiates. Furthermore, our developmental interrogation revealed tissue moduli, mechanical constraints, and emergent thickness variations in the EGL that are fundamentally inconsistent with traditional elastic bilayer wrinkling models. Therefore our results call for a new understanding of brain folding.

By applying a multi-phase model constrained by our measured data we were able to capture the correct shape variations and number of folds at the onset of folding. Our new framework accounts for: the rapidly expanding fluid EGL, whose thickness is proposed to be regulated by Bergmann glial fibers, the slower growing incompressible core, and fibrous material in the form of glial fibers and possibly axons as well as the meninges that potentially provide radial and circumferential tension (Figure 8). This multi-phase model of folding makes many new predictions. One such prediction is that adjusting the amount of tension spanning the cerebellum will change the degree of folding. Indeed, alterations of the cells that likely create tension-based forces could explain the dramatically disrupted folding seen in mouse mutants in which radial glia do not produce Bergmann glia (Li et al., 2014). Without Bergmann glia, the EGL would be expected to not form a layer with regular thickness and it should be more sensitive to variations in radial glial tension. Consistent with this prediction, mutants without Bergmann glia have more localized and less regular folds (Li et al., 2014) Our combination of experimental studies and modeling thus provide new insights into cerebellar folding, including an underappreciated role for tension.

Figure 8

Download asset Open asset

Uniform cell behaviors within a fluid-like EGL create differential expansion based folding approximated by a multiphase model with radial and circumferential tension.

(A) Cartoon of E16.5 cerebellum showing the EGL (dark purple) overlying the in incompressible core (light purple) with fibers (lines) of radial glia (closed ovals) and Bergmann glia (open ovals) crossing the cerebellum and the EGL. (B) Map of stress within the cerebellum at the onset of foliation. (C) Schematics showing that an AC is first detected as a regional inward thickening of the EGL (left). The constraining tensions shape the fluid-like EGL such that the EGL becomes thicker at the AC (middle). As radial glia transition to Bergmann glia, modeling predicts a reduction in radial tension (right). (D) Since ACs hold their position in space, and compartmentalize the cells within the EGL, we propose that they behave as mechanical boundaries allowing local domains of differential expansion to arise and progressive folding to occur.

https://doi.org/10.7554/eLife.45019.018

Under the new framework revealed by our measurements made in the developing mouse cerebellum, to approximate the observed shape changes in the murine cerebellum from E17.5 to E18.5 the ratio of the core stiffness over the radial tension must increase. Yet, the measured bulk modulus of the core shows no increase during development. Therefore a second prediction is that radial tension must decrease during development. While the cerebellum is crossed by many fibers at folding initiation, radial glial fibers are an attractive candidate to mediate this change in radial tension (Sillitoe and Joyner, 2007; Rahimi-Balaei et al., 2015). First, they span from the VZ to the surface of the cerebellum at E16.5. Additionally, during folding initiation the radial glia undergo a transition into Bergmann glia where they release their basal connection to the VZ and the cell body migrates towards the surface (Mota and Herculano-Houzel, 2015). This transition could lead to a reduction in the global radial tension and thus would be consistent with our model prediction.

The mechanics underlying hierarchical folding remain an open challenge. However, our developmental data may provide a way forward. As ACs maintain their spatial positions, and as they compartmentalize granule cells within the EGL into the lobule regions (Legué et al., 2015), we propose that they create fixed mechanical boundaries that divide the cerebellum into self-similar domains. These domains, with their similar physical properties to the initial unfolded cerebellum, can then undergo additional folding. Furthermore since ACs compartmentalize granule cells within the lobule regions, once separated the lobule regions can develop distinct characteristics, like the observed differential proliferation rates at E18.5. We speculate, therefore, that the folding patterns seen across cerebella in different species evolved by adjustment of global as well as regional levels of differential expansion and tension which ultimately mold the functionality of the cerebellum.

Finally it is interesting to note the similarities and differences between the developing cerebellum and the cerebral cortex. Radial glia span the entire cerebral cortex just as in the cerebellum (Götz et al., 2002). Furthermore, species with folded cerebrums have evolved outer radial glial cells for which the cell body leaves the ventricular zone to become positioned near the surface while retaining fibers anchored on the surface, similar to Bergmann glia in the cerebellum (Leung and Li, 2018; Reillo et al., 2011). While we have emphasized the notion of tension via glial fibers in the developing cerebellum, axonal tension has been discussed in the context of shaping the developing cerebrum (Van Essen, 1997). Tissue cutting in the cerebral cortex of ferrets has revealed a similar tension pattern during folding as we found in the cerebellum (Xu et al., 2010). We therefore submit that our work calls for a revival of the notion of how tension affects the shape of the developing cerebrum.

Unlike the cerebellum, the cerebral cortex is not divided into a simple bilayer system. However, outer radial glial cells proliferate, much like the GCPs of the EGL, to drive the expansion of the outer regions of the cerebral cortex around the time of initiation of folding (Hansen et al., 2010; Heng et al., 2017; Nowakowski et al., 2016). Moving the zone of proliferation out from the VZ gives more space for the increased proliferation required in folding systems. The cerebellum, housing 80% of the neurons in the human brain may be an extreme example requiring the region of proliferation to be completely on the outer surface (Andersen et al., 1992). Constraining models of folding of different brain regions with developmental data will bring about a more accurate quantitative understanding of the shaping of the developing brain.

Materials and methods

Key resources table

Reagent type (species) or resource	Designation	Source or reference	Identifiers	Additional information
Genetic reagent (M.Musculus)	FVB/N	Jax Mice, The Jackson Laboratory	Jackson Labs Stock number: 001800
Genetic reagent (M.Musculus)	Atoh1-CreER	Machold and Fishell, 2005	Jackson Labs Stock number: 007684
Genetic reagent (M.Musculus)	Nestin-CreER	Imayoshi et al., 2006
Genetic reagent (M.Musculus)	Rosa26^MTMG	Muzumdar et al., 2007	Jackson labs Stock number: 007676
Genetic reagent (M.Musculus)	Rosa26^Ai75	Daigle et al., 2018	Jackson labs Stock number: 025106
Antibody	mouse monoclonal anti-P27	BD Pharmingen	610241	Dilution 1:500
Antibody	rabbit polyclonal anti-GFP	Life Technologies	A11122	Dilution 1:500
Antibody	rat monoclonal anti-GFP	Nacalai Tesque	04404–84	Dilution 1:500
Commercial Assay or Kit	EDU	Invitrogen	C10340

Animals

All experiments were performed following protocols approved by Memorial Sloan Kettering Cancer Center’s Institutional Animal Care and Use Committee. The inbred FVB/N stain was used for all proliferation rate, area, length, and expansion rate measurements. Atoh1-CreER (Machold and Fishell, 2005), Nestin-CreER (Imayoshi et al., 2006), Rosa26^MTMG (Muzumdar et al., 2007), Rosa26^Ai75 (Daigle et al., 2018) were used to quantify cell shape and size as well as fiber distribution and were maintained on the outbred Swiss Webster background. The Swiss Webster strain was used for scanning acoustic microscopy. Both sexes were used for the analysis. Animals were kept on a 12 hr light/dark cycle and food and water were supplied ad libitum. All experiments were performed following protocols approved by Memorial Sloan Kettering Cancer Center’s Institutional Animal Care and Use Committee.

The appearance of a vaginal plug set noon of the day as Embryonic day 0.5 (E0.5). All animals were collected within two hours of noon on the day of collection. Tamoxifen (Tm, Sigma-Aldrich) was dissolved in corn oil (Sigma-Aldrich) at 20 mg/mL. Pregnant females carrying litters with Atoh1-CreER/+;R26^MTMG/MTMG or NestinCER/+;R26^MTMG/MTMG embryos were given one 20 μg/g dose of TM via subcutaneous injection two days prior to analysis. 25 μg/g of 5-ethynyl-2-deoxyruidine (EDU; Invitrogen) was administered via subcutaneous injection one hour prior to collection.

Tissue processing, immunohistochemistry, and imaging

Request a detailed protocol

For embryonic stages heads were fixed in 4% paraformaldehyde overnight at 4°C. For postnatal animals, the brain was dissected out first before fixation. Tissues were stored in 30% sucrose. For all proliferation, area, length, and thickness measurements brains were embedded in optimal cutting temperature (OCT) compound. Parasagittal sections were collect with a Leica cryostat (CM3050s) at 10 μm.

Prior to IHC, EdU was detected using a commercial kit (Invitrogen, C10340). Following EdU reaction the following primary antibodies were used either overnight at 4°C or 4 hr at room temperature: mouse anti-P27 (BD Pharmingen, 610241), rabbit anti-GFP (Life Technologies, A11122), rat anti-GFP (Nacalai Tesque, 04404–84). All antibodies were diluted to 1:500 in 2% milk (American Bioanalytical) and 0.2% Triton X-100 (Fisher Scientific). Alexa Fluor secondary antibodies (1:500; Invitrogen) were used: Alexa Fluor 488 donkey anti-rabbit, A21206, Alexa Fluor 488 donkey anti-rat, A21208, Alexa Fluor 488 donkey anti-mouse, A21202, Alexa Fluor 647 donkey anti-mouse, A31571. EdU was detected using a commercial kit (Invitrogen, C10340).

For cell size, shape and fiber density analysis 60 μm parasagittal sections were collected on a Leica vibratome (VT100S). Primary and secondary antibodies were diluted 1:500 in 2% milk and incubated overnight at 4°C.

For scanning acoustic microscopy brains were processed for paraffin embedding and parasagittal sections of 10 μm thick were collected on a microtome (Leica RM2255). Structured illumination and confocal Imaging was done with Zeiss Observer Z.1 with Apatome or Zeiss LSM 880 respectively.

Quantification of proliferation, length, area, folding index and thickness

Request a detailed protocol

Measurements for all analysis were taken from the three most midline sagittal sections and averaged. The most midline section was determined by dividing the distance in half between the lateral edges where the third ventrical and the mesencephalic vesicle are no longer connected. Quantifications were made using Imaris (Bitplane) and MATLAB (Mathworks) software.

EGL Proliferation rate was calculated as EDU+/(Dapi+;P27-) cells. All cells were counted within the lobule region to the midpoint of the Anchoring Centers. For proliferation measurements through the ACs and the surrounding EGL at E16.5 and E17.5 a 50 μm window measured from the outer surface of the EGL was centered at the AC. The measuring window was centered at every 25 μm anterior and posterior to the EGL for a total distance of 250 μm anterior and posterior to the AC. At E18.5, when the AC is fully formed, everything proximal to the centroid of the cerebellum under the midpoint of the AC was counted as the AC. Non-overlapping regions of 50 μm also were measured in either direction for a total of 200 μM anterior and posterior to the AC. Proliferation was measured in three cerebella at E16.5 and E17.5 and in four cerebella at E18.5.

EGL length was measured from the outer surface of the EGL following the curvature of the EGL. Cerebellar area was calculated as the area within the outer surface of the EGL and the ventricular zone. A short strait edge was made perpendicular to the ventricular zone to close the area back upon to the anterior end of the EGL. The convex curvature of the cerebellum was measured by following only the positive curvature of the EGL. The folding index was determined as FI = 1 - (Positive curvature/EGL length). Data collected for E16.5, E17.5, E18.5 and P0 came from 6,8,7, and 9 cerebella respectively.

EGL thickness was measured by defining the outer and inner curvature of the EGL. The shortest distance lines were drawn to the outer curvature from discrete points distributed at every 12.5 μm along the inner curvature of the EGL. Nine ACs and surrounding regions from five cerebella were quantified at E16.5 and 13 ACs from five cerebellar were analyzed at E17.5. At E18.5 six ACs from two cerebellar were quantified.

Quantification of cell shape

Request a detailed protocol

Midline sections were imaged with a Zeiss LSM 880. Serial images were taken to cover the entire EGL of lobule regions L123, L45, and L678 and the ACs. Manual cell masks were created with Imaris software defining the curvature at every z-slice. Every cell that was completely included in the imaging window and that was distinguishable from surrounding cells was counted to reduce sampling bias. Cells from three brains were measured at each stage for a total of 131 at E16.5 and 201 at E18.5. Shape was defined via sphericity, which is the surface area of a sphere having the same volume as the cell of interest divided by the surface area of the cell of interest.

Quantification of fibers within the EGL

Request a detailed protocol

Midline sections were imaged with a Zeiss LSM 880. Image tiling was used to cover the EGL. Using Python, a fourth or fifth order polynomial was fitted to the outer edge of the EGL in each image, and five scan lines were positioned at 12.2 μm intervals beneath the surface, and parallel to it. A bin width of 50 μm as measured along the polynomial contour was centered at the AC. Bins of equal distance were extended both anteriorly and posteriorly. Staining intensity was counted along each scan line at every z-slice of the confocal stack. Each image was normalized to the mean intensity and smoothed with a Gaussian filter. Peak counting was done using minimum and maximum filters, keeping neighborhood size and threshold parameters constant for all images. The results from the five scan lines were averaged.

Tissue cutting

Request a detailed protocol

Live cerebella of E16.5 FVB/N mice were collected in dissection buffer as previously described (Wojcinski et al., 2017) and embedded in low-melting point agarose (Invitrogen). Sagittal slices at a thickness of 250 μm were collected. Slices were removed from the agarose and place in petri-dishes coated with Poly(2-hydroxyethyl methacrylate)(Sigma-Aldrich). Tissue cuts (eight horizontal, 10 radial) were made with a 30° Premier Edge stab knife (Oasis Medical). Slices were allowed to relax for 10 min. Time-lapse images were acquired on a Leica MZ75 dissection scope.

Live imaging analysis

Request a detailed protocol

Live cerebella of E16.5 Atoh1-CreER/+; R26Ai75/+mice were collected and slices of a thickness of 250 μm were cultured on Millicell cell culture inserts (Millipore) in glass bottom plates (Matek) as previously described (Wojcinski et al., 2017). Image stacks were acquired on a Zeiss LSM 880 at intervals of around 3.5 min for up to 5 hr. Cell positions were tracked using Imaris (Bitplane) software. Three time-lapses were analzyed.

Scanning acoustic microscopy

Request a detailed protocol

Mechanical tissue properties were analyzed using a 250 MHz Scanning Acoustic Microscope (SAM), described previously (Rohrbach et al., 2015; Rohrbach et al., 2018; Rohrbach and Mamou, 2018). Briefly, 12 μm paraffin sections of mouse embryonic brains were de-parafinized, hydrated in de-ionized water and raster scanned (2 µm steps in both direction) on the SAM to acquire radio-frequency (RF) ultrasound data. At each scan location, signal processing was performed to compute the amplitude, sample thickness, speed of sound, acoustic impedance, attenuation, bulk modulus, and mass density (Rohrbach and Mamou, 2018). Two-dimensional maps of tissues properties were formed using the values obtained at each scan location. Bulk modulus was computed from the product of the acoustic impedance and the speed of sound. Co-registered histology and SAM amplitude images were used to identify regions-of-interest (ROIs) corresponding to the EGL layer and underlying core of the cerebellum in each sample. Bulk modulus was analyzed as a measure of tissue stiffness: ROI measurements were acquired from 3 sections from three embryos at each developmental stage.

Finite element simulations

Request a detailed protocol

The wrinkle of a circular bilayer structure in Figure 3a was simulated with commercial software ABAQUS. Both film and substrate were modeled as incompressible neo-Hookean materials. The ratio between shear moduli of the film and substrate was 50 and the initial radius of the simulated structure was 16 times that of the film thickness. The differential growth of the EGL and core was modeled by an isotropic expansion of the film in the bilayer structure.

To test the elastic wrinkling model, we conducted finite element (FE) simulations for bilayer structures with a film bonded on a substrate, which represents the EGL layer and core structure, respectively. The structures were assumed to be under 2D plane strain deformation to mimic the quasi-2D nature of cerebellum wrinkles. Neo-Hookean model was adopted to describe the elastic properties of both film and substrate, whose strain energy can be expressed as

U = \frac{1}{2} μ (I_{1} - 3)

where $μ$ is the shear modulus and $I_{1}$ represents the first invariant of the right Cauchy-Green strain tensor. The Poisson’s ratios for the film and substrate were set to be 0.5, based on experimental observations that the bulk modulus of EGL and core are in the order of GPa, much larger than the shear modulus of soft tissues (~ kPa).

We carried out FE simulations through commercial software ABAQUS. A second order 6 node hybrid element (CPE6MH) was utilized to discretize the film and substrate. Very fine FE meshes were used to make sure the results independent of mesh size. To incorporate differential growth in real EGL layer and core, an isotropic growth deformation tension was applied to the modeled film by decoupling the deformation tenor $F$ into elastic deformation part $A$ and growth part $G$ .

F = A ∙ G

For simplicity, we assume the growth part is isotropic and controlled by a scalar variable $g$

G = g [\begin{matrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{matrix}]

where $g > 1$ represents a faster growth in EGL than the core. To trigger instabilities in numerical simulations, random perturbations (e.g., White Gaussian noise with 0.001t mean magnitude) were applied to the nodal positions at the top surface of the film and the interface between the film and substrate.

To qualitatively understand the cut experiments we ran a FE simulation of a pre-cut circular bilayer structure and then assigned swelling strain to the film. This neglected the dynamical process in the real cut experiments and only focused on the final equilibrium of the cerebellum after long time relaxation. All the simulation parameters were the same as those in the wrinkling simulation. The initial cut length a is equal to 8 t. The minimum in-plane principal stress corresponds to the hoop stress in the film.

Details of multi-phase model as applied to initiation of cerebellar folds

Request a detailed protocol

For a full treatment of the mathematics please see Engstrom et al. (2018).

We, formulated a two-dimensional model based on the parameters of a midsagittal section of the cerebellum. The distance of the outer edge of the EGL and, hence, the outer edge of the cerebellum from the center of the cerebellum was defined as r(θ) with θ as the angular coordinate. We assumed that r(θ) was single-valued. The thickness of the EGL was defined as t(θ). See model schematic in.

Taking into account the four assumptions discussed in the main text, we constructed the following energy functional to be minimized

E [r, t, \frac{d t}{d θ}] = \int d θ {k_{r} (r - r_{0})^{2} - k_{t} (t - t_{0})^{2} + β ({\frac{d t}{d θ}}^{2})},

with k_r as the stiffness modulus (a spring constant in one-dimension) of the radial glial fibers and the pial surface contained in the meninges surrounding the cerebellum since the cerebellar radius is proportional to its perimeter, r₀ as the preferred radius of the cerebellum, k_t denoting a growth potential due to cell proliferation, t₀ as thickness of the EGL (cortex), and, β quantified the mechanical resistance to changing the thickness of the EGL. Given our first assumption of an incompressible cerebellar core, we imposed the constraint

\frac{1}{2} \int d θ {(r - t)}^{2} = A_{0},

with A₀ as a preferred cerebellar area. We applied the variational principle to minimize the energy functional subject to the core constraint, that is

δ (E - μ \int d θ {(r - t)}^{2}) = 0,

where $μ$ is a Lagrange multiplier. Assuming the preferred radius of the cerebellum is constant and the thickness of the EGL/cortex is also constant, then the preferred cerebellar shape was a circle and the EGL an annulus.

The variational analysis yielded the following equation of shape for t(θ);

\frac{d^{2} t}{d θ^{2}} + q^{2} t (θ) = \frac{k_{t}}{β} (t_{0} + \frac{\frac{μ r_{0}}{k_{t}}}{1 - \frac{μ}{k_{r}}}),

with $q^{2} = \frac{k_{t}}{β} (1 + \frac{\frac{μ}{k_{t}}}{1 - \frac{μ}{k_{r}}}) .$ The solution to the equation of shape was

t (θ) = A_{t} \sin (q θ + ϕ) + C_{1} (r_{0}, t_{0}, k_{r}, k_{t}, μ),

with C₁ independent of θ and $A_{t} = \sqrt{2} (1 - \frac{μ}{k_{r}}) \sqrt{\frac{A_{0}}{π} - C_{2} (r_{0}, t_{0}, k_{r}, k_{t}, μ)}$ such that $A_{0} > π C_{2} .$ There was an additional equation of shape for r(θ) from the variational principle that depended on t(θ) and so was determined

r (θ) = - \frac{\frac{μ}{k_{r}}}{1 - \frac{μ}{k_{r}}} A_{t} \sin (q θ + ϕ) + C_{3} (r_{0}, k_{r}, μ) .

We used the measured data at E16.5 to set the parameters to make predictions for the shape of both the EGL and core (and so the relationship between the two) at later times. Because we are primarily interested in shape changes, rather than size changes, a nondimensionalized model solution was used, that is we chose units where r₀ = 1. This reduces the total number of parameters specifying the model to five dimensionless parameters. Plots assumed a circular preferred shape, and with other parameters as follows: $ϵ = μ / k_{r}$ is shown in Figure 6b,c, $c = k_{r} / k_{t} = 0.06 / ϵ$ , $A_{t} / r_{0} = ϵ / 9.6$ , $t_{0} / r_{0} = ϵ / 4.8$ , and $q = 6$ . Note that for $ϵ = 0.3$ , these parameters are numerically consistent with our E16.5 measurements: $A_{t} / t_{0} = 0.5$ and $r_{0} / t_{0} = 16$ , as well as the observed number of invaginations in the half circle: $q / 2 = 3$ . All of these parameters are either constant or depend on the time-like parameter $ϵ$ . One of these dependencies has a functional form that is physically justifiable ( $A_{t} ~ ϵ$ ), another has a form that is biologically justifiable ( $c ~ 1 / ϵ$ ), owing to the decrease in the number of radial glia over time.

We defined a dimensionless 'shape factor' as half of the perimeter divided by the square root of half of the area as appropriate for a semi-circle. To compare the model’s predictive deviation of this quantity form the semi-circular value we assumed a linear relationship between $ϵ$ and time T measured in embryonic days: $ϵ$ (T) = 0.3(T-15.5).

Statistical analyses

Request a detailed protocol

Statistical analyses were performed using Matlab software. Significance was determined at p<0.05. Two-way ANOVA was used for proliferation analysis as two variables were tracked, mouse and region. Cell shape, volume, fiber distribution, EGL thickness and bulk modulus were run under a standard ANOVA. After ANOVA analysis a multiple comparison was run with Tukey’s honestly significant difference criterion. F-test for variance and two-tailed student’s paired t-test were used for slice cutting and relaxation quantifications. The degrees of freedom, where appropriate, and P values are given in the figure legends. All error bars are standard deviations. No statistical methods were used to predetermine the sample sizes. We used sample sizes aligned with the standard in the field. No randomization was used nor was data collection or analysis performed blind.

Data availability

All data is presented in the manuscript and supporting files.

References

Book
1. Allen HG
(1969)
Analysis and Design of Structural Sandwich Panels (1st ed)

Oxford, New York: Pergamon Press.
- Google Scholar
(1992) A quantitative study of the human cerebellum with unbiased stereological techniques
The Journal of Comparative Neurology 326:549–560.

https://doi.org/10.1002/cne.903260405
- PubMed
- Google Scholar
1. Bayly PV
2. Okamoto RJ
3. Xu G
4. Shi Y
5. Taber LA
(2013) A cortical folding model incorporating stress-dependent growth explains gyral wavelengths and stress patterns in the developing brain
Physical Biology 10:016005.

https://doi.org/10.1088/1478-3975/10/1/016005
- PubMed
- Google Scholar
(2014) Mechanical forces in cerebral cortical folding: a review of measurements and models
Journal of the Mechanical Behavior of Biomedical Materials 29:568–581.

https://doi.org/10.1016/j.jmbbm.2013.02.018
- PubMed
- Google Scholar
1. Budday S
2. Sommer G
3. Haybaeck J
4. Steinmann P
5. Holzapfel GA
6. Kuhl E
(2017) Rheological characterization of human brain tissue
Acta Biomaterialia 60:315–329.

https://doi.org/10.1016/j.actbio.2017.06.024
- PubMed
- Google Scholar
1. Daigle TL
2. Madisen L
3. Hage TA
4. Valley MT
5. Knoblich U
6. Larsen RS
7. Takeno MM
8. Huang L
9. Gu H
10. Larsen R
11. Mills M
12. Bosma-Moody A
13. Siverts LA
14. Walker M
15. Graybuck LT
16. Yao Z
17. Fong O
18. Nguyen TN
19. Garren E
20. Lenz GH
21. Chavarha M
22. Pendergraft J
23. Harrington J
24. Hirokawa KE
25. Harris JA
26. Nicovich PR
27. McGraw MJ
28. Ollerenshaw DR
29. Smith KA
30. Baker CA
31. Ting JT
32. Sunkin SM
33. Lecoq J
34. Lin MZ
35. Boyden ES
36. Murphy GJ
37. da Costa NM
38. Waters J
39. Li L
40. Tasic B
41. Zeng H
(2018) A suite of transgenic driver and reporter mouse lines with enhanced Brain-Cell-Type targeting and functionality
Cell 174:465–480.

https://doi.org/10.1016/j.cell.2018.06.035
- PubMed
- Google Scholar
(2018) Buckling without bending: a new paradigm in morphogenesis
Physical Review X 8:.

https://doi.org/10.1103/PhysRevX.8.041053
- Google Scholar
(2002) Radial glial cells as neuronal precursors: a new perspective on the correlation of morphology and lineage restriction in the developing cerebral cortex of mice
Brain Research Bulletin 57:777–788.

https://doi.org/10.1016/S0361-9230(01)00777-8
- PubMed
- Google Scholar
(2012) Mechanical instabilities of biological tubes
Physical Review Letters 109:018101.

https://doi.org/10.1103/PhysRevLett.109.018101
- PubMed
- Google Scholar
(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
(2014) The roles and regulation of multicellular rosette structures during morphogenesis
Development 141:2549–2558.

https://doi.org/10.1242/dev.101444
- PubMed
- Google Scholar
(2014) Apical constriction drives tissue-scale hydrodynamic flow to mediate cell elongation
Nature 508:392–396.

https://doi.org/10.1038/nature13070
- PubMed
- Google Scholar
1. Heng X
2. Guo Q
3. Leung AW
4. Li JY
(2017) Analogous mechanism regulating formation of neocortical basal radial glia and cerebellar bergmann glia
eLife 6:e23253.

https://doi.org/10.7554/eLife.23253
- PubMed
- Google Scholar
1. Hohlfeld E
2. Mahadevan L
(2011) Unfolding the sulcus
Physical Review Letters 106:105702.

https://doi.org/10.1103/PhysRevLett.106.105702
- PubMed
- Google Scholar
(2006) Temporal regulation of cre recombinase activity in neural stem cells
Genesis 44:233–238.

https://doi.org/10.1002/dvg.20212
- PubMed
- Google Scholar
(2018) Human brain organoids on a chip reveal the physics of folding
Nature Physics 14:515–522.

https://doi.org/10.1038/s41567-018-0046-7
- PubMed
- Google Scholar
(2015) Clonal analysis reveals granule cell behaviors and compartmentalization that determine the folded morphology of the cerebellum
Development 142:1661–1671.

https://doi.org/10.1242/dev.120287
- PubMed
- Google Scholar
(2016) Differential timing of granule cell production during cerebellum development underlies generation of the foliation pattern
Neural Development 11:17.

https://doi.org/10.1186/s13064-016-0072-z
- PubMed
- Google Scholar
1. Lejeune E
2. Javili A
3. Weickenmeier J
4. Kuhl E
5. Linder C
(2016) Tri-layer wrinkling as a mechanism for anchoring center initiation in the developing cerebellum
Soft Matter 12:5613–5620.

https://doi.org/10.1039/C6SM00526H
- PubMed
- Google Scholar
1. Leto K
2. Arancillo M
3. Becker EB
4. Buffo A
5. Chiang C
6. Ding B
7. Dobyns WB
8. Dusart I
9. Haldipur P
10. Hatten ME
11. Hoshino M
12. Joyner AL
13. Kano M
14. Kilpatrick DL
15. Koibuchi N
16. Marino S
17. Martinez S
18. Millen KJ
19. Millner TO
20. Miyata T
21. Parmigiani E
22. Schilling K
23. Sekerková G
24. Sillitoe RV
25. Sotelo C
26. Uesaka N
27. Wefers A
28. Wingate RJ
29. Hawkes R
(2016) Consensus paper: cerebellar development
The Cerebellum 15:789–828.

https://doi.org/10.1007/s12311-015-0724-2
- PubMed
- Google Scholar
1. Leung AW
2. Li JYH
(2018) The molecular pathway regulating Bergmann Glia and folia generation in the cerebellum
The Cerebellum 17:42–48.

https://doi.org/10.1007/s12311-017-0904-3
- PubMed
- Google Scholar
1. Li K
2. Leung AW
3. Guo Q
4. Yang W
5. Li JY
(2014) Shp2-dependent ERK signaling is essential for induction of Bergmann Glia and foliation of the cerebellum
The Journal of Neuroscience 34:922–931.

https://doi.org/10.1523/JNEUROSCI.3476-13.2014
- PubMed
- Google Scholar
1. Machold R
2. Fishell G
(2005) Math1 is expressed in temporally discrete pools of cerebellar rhombic-lip neural progenitors
Neuron 48:17–24.

https://doi.org/10.1016/j.neuron.2005.08.028
- PubMed
- Google Scholar
1. Mammoto T
2. Ingber DE
(2010) Mechanical control of tissue and organ development
Development 137:1407–1420.

https://doi.org/10.1242/dev.024166
- PubMed
- Google Scholar
1. Mota B
2. Herculano-Houzel S
(2015) Cortical folding scales universally with surface area and thickness, not number of neurons
Science 349:74–77.

https://doi.org/10.1126/science.aaa9101
- PubMed
- Google Scholar
1. Muzumdar MD
2. Tasic B
3. Miyamichi K
4. Li L
5. Luo L
(2007) A global double-fluorescent cre reporter mouse
Genesis 45:593–605.

https://doi.org/10.1002/dvg.20335
- PubMed
- Google Scholar
1. Nelson CM
(2016) On buckling morphogenesis
Journal of Biomechanical Engineering 138:021005.

https://doi.org/10.1115/1.4032128
- PubMed
- Google Scholar
(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
(2015) Embryonic stages in cerebellar afferent development
Cerebellum & Ataxias 2:7.

https://doi.org/10.1186/s40673-015-0026-y
- PubMed
- Google Scholar
(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
(1975) Mechanical model of brain convolutional development
Science 189:18–21.

https://doi.org/10.1126/science.1135626
- PubMed
- Google Scholar
(2015) Fine-resolution maps of acoustic properties at 250 MHz of unstained fixed murine retinal layers
The Journal of the Acoustical Society of America 137:EL381–EL387.

https://doi.org/10.1121/1.4916790
- PubMed
- Google Scholar
1. Rohrbach D
2. Silverman RH
3. Chun D
4. Lloyd HO
5. Urs R
6. Mamou J
(2018) Improved High-Frequency ultrasound corneal biometric accuracy by Micrometer-Resolution Acoustic-Property maps of the cornea
Translational Vision Science & Technology 7:21.

https://doi.org/10.1167/tvst.7.2.21
- PubMed
- Google Scholar
1. Rohrbach D
2. Mamou J
(2018) Autoregressive signal processing applied to High-Frequency acoustic microscopy of soft tissues
IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control 65:2054–2072.

https://doi.org/10.1109/TUFFC.2018.2869876
- PubMed
- Google Scholar
1. Ronan L
2. Voets N
3. Rua C
4. Alexander-Bloch A
5. Hough M
6. Mackay C
7. Crow TJ
8. James A
9. Giedd JN
10. Fletcher PC
(2014) Differential tangential expansion as a mechanism for cortical gyrification
Cerebral Cortex 24:2219–2228.

https://doi.org/10.1093/cercor/bht082
- PubMed
- Google Scholar
1. Shyer AE
2. Tallinen T
3. Nerurkar NL
4. Wei Z
5. Gil ES
6. Kaplan DL
7. Tabin CJ
8. Mahadevan L
(2013) Villification: how the gut gets its villi
Science 342:212–218.

https://doi.org/10.1126/science.1238842
- PubMed
- Google Scholar
1. Sillitoe RV
2. Joyner AL
(2007) Morphology, molecular codes, and circuitry produce the three-dimensional complexity of the cerebellum
Annual Review of Cell and Developmental Biology 23:549–577.

https://doi.org/10.1146/annurev.cellbio.23.090506.123237
- PubMed
- Google Scholar
1. Stemple DL
(2005) Structure and function of the notochord: an essential organ for chordate development
Development 132:2503–2512.

https://doi.org/10.1242/dev.01812
- PubMed
- Google Scholar
1. Sudarov A
2. Joyner AL
(2007) Cerebellum morphogenesis: the foliation pattern is orchestrated by multi-cellular anchoring centers
Neural Development 2:26.

https://doi.org/10.1186/1749-8104-2-26
- PubMed
- Google Scholar
1. Szulc KU
2. Lerch JP
3. Nieman BJ
4. Bartelle BB
5. Friedel M
6. Suero-Abreu GA
7. Watson C
8. Joyner AL
9. Turnbull DH
(2015) 4d MEMRI atlas of neonatal FVB/N mouse brain development
NeuroImage 118:49–62.

https://doi.org/10.1016/j.neuroimage.2015.05.029
- PubMed
- Google Scholar
(2014) Gyrification from constrained cortical expansion
PNAS 111:12667–12672.

https://doi.org/10.1073/pnas.1406015111
- PubMed
- Google Scholar
(2016) On the growth and form of cortical convolutions
Nature Physics 12:588–593.

https://doi.org/10.1038/nphys3632
- Google Scholar
1. Van Essen DC
(1997) A tension-based theory of morphogenesis and compact wiring in the central nervous system
Nature 385:313–318.

https://doi.org/10.1038/385313a0
- PubMed
- Google Scholar
1. Wiggs BR
2. Hrousis CA
3. Drazen JM
4. Kamm RD
(1997) On the mechanism of mucosal folding in normal and asthmatic airways
Journal of Applied Physiology 83:1814–1821.

https://doi.org/10.1152/jappl.1997.83.6.1814
- PubMed
- Google Scholar
1. Wojcinski A
2. Lawton AK
3. Bayin NS
4. Lao Z
5. Stephen DN
6. Joyner AL
(2017) Cerebellar granule cell replenishment postinjury by adaptive reprogramming of nestin⁺ progenitors
Nature Neuroscience 20:1361–1370.

https://doi.org/10.1038/nn.4621
- PubMed
- Google Scholar
1. Xu G
2. Knutsen AK
3. Dikranian K
4. Kroenke CD
5. Bayly PV
6. Taber LA
(2010) Axons pull on the brain, but tension does not drive cortical folding
Journal of Biomechanical Engineering 132:071013.

https://doi.org/10.1115/1.4001683
- PubMed
- Google Scholar
1. Yamada K
2. Watanabe M
(2002) Cytodifferentiation of Bergmann Glia and its relationship with purkinje cells
Anatomical Science International 77:94–108.

https://doi.org/10.1046/j.0022-7722.2002.00021.x
- PubMed
- Google Scholar
1. Yuasa S
(1996) Bergmann glial development in the mouse cerebellum as revealed by tenascin expression
Anatomy and Embryology 194:223–234.

https://doi.org/10.1007/BF00187133
- PubMed
- Google Scholar
1. Zhao R
2. Zhao X
(2017) Multimodal surface instabilities in curved Film–Substrate Structures
Journal of Applied Mechanics 84:081001–081013.

https://doi.org/10.1115/1.4036940
- Google Scholar

Article and author information

Author details

Andrew K Lawton

Developmental Biology Program, Sloan Kettering Institute, New York, United States

Contribution
Conceptualization, Formal analysis, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing—original draft, Project administration, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-8633-6637
Tyler Engstrom

Department of Physics, Syracuse University, Syracuse, United States

Contribution
Conceptualization, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing—review and editing

Competing interests
No competing interests declared
Daniel Rohrbach

Lizzi Center for Biomedical Engineering, Riverside Research, New York, United States

Contribution
Formal analysis, Supervision, Investigation, Methodology, Writing—review and editing

Competing interests
No competing interests declared
Masaaki Omura
1. Lizzi Center for Biomedical Engineering, Riverside Research, New York, United States
2. Department of Radiology, Skirball Institute of Biomolecular Medicine, NYU School of Medicine, New York, United States
3. Graduate School of Science and Engineering, Chiba University, Chiba, Japan
Contribution
Formal analysis, Investigation

Competing interests
No competing interests declared
Daniel H Turnbull

Department of Radiology, Skirball Institute of Biomolecular Medicine, NYU School of Medicine, New York, United States

Contribution
Supervision, Project administration, Writing—review and editing

Competing interests
No competing interests declared
Jonathan Mamou

Lizzi Center for Biomedical Engineering, Riverside Research, New York, United States

Contribution
Resources, Supervision, Project administration, Writing—review and editing

Competing interests
No competing interests declared
Teng Zhang

Department of Mechanical & Aerospace Engineering, Syracuse University, Syracuse, United States

Contribution
Formal analysis, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing—review and editing

Competing interests
No competing interests declared
J M Schwarz

Department of Physics, Syracuse University, Syracuse, United States

Contribution
Conceptualization, Supervision, Funding acquisition, Investigation, Methodology, Project administration, Writing—review and editing

Competing interests
No competing interests declared
Alexandra L Joyner
1. Developmental Biology Program, Sloan Kettering Institute, New York, United States
2. Biochemistry, Cell and Molecular Biology Program, Weill Graduate School of Medical Sciences, Cornell University, New York, United States
Contribution
Conceptualization, Supervision, Funding acquisition, Writing—original draft, Project administration, Writing—review and editing

For correspondence
joynera@mskcc.org

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-7090-9605

Funding

National Institute of Neurological Disorders and Stroke (Postdoctoral Training Fellowship: F32NS086163)

Andrew K Lawton

National Institute of Biomedical Imaging and Bioengineering (R21EB016117)

Jonathan Mamou

National Science Foundation (Comet Cluster in the Extreme Science and Engineering Discovery Environment: TG-MSS170004)

Teng Zhang

National Science Foundation (NSF-DMR-CMMT 1507938)

J M Schwarz

National Science Foundation (NSF-PHY-PoLS 1607416)

J M Schwarz

National Institute of Mental Health (R37MH085726)

Alexandra L Joyner

National Cancer Institute (Cancer Center Support Grant: P30 CA008748-48)

Alexandra L Joyner

National Institute of Mental Health (R01NS092096)

Alexandra L Joyner

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We are grateful to Anna-Katerina Hadjantonakis and Songhai Shi for use of Imaris software. We thank Jennifer Zallen and Anna-Katerina Hadjantonakis for experimental advice, Nathanael Kim for help with the acoustic microscopy and Professor Tadashi Yamagushi for his support for MO’s visit to New York. We appreciate the discussions we have had with Alexandre Wojcinski and the entire Joyner Laboratory, and the administrative support from Cara Monaco.

Ethics

Animal experimentation: All experiments were performed following protocols approved by Memorial Sloan Kettering Cancer Center's Institutional Animal Care and Use Committee. Protocol number: 07-01-001 Mouse Developmental Genetics

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

2,787

views
362

downloads
30

citations

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

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Andrew K Lawton
Tyler Engstrom
Daniel Rohrbach
Masaaki Omura
Daniel H Turnbull
Jonathan Mamou
Teng Zhang
J M Schwarz
Alexandra L Joyner

(2019)

Cerebellar folding is initiated by mechanical constraints on a fluid-like layer without a cellular pre-pattern

eLife 8:e45019.

https://doi.org/10.7554/eLife.45019

Categories and tags

Research organism

Mouse

Share this article

Cite this article

Initiation of cerebellar folding correlates with timing of differential expansion.

Measured tissue stiffness, stress, and shape at folding initiation are inconsistent with wrinkling models.

Live slice cutting and relaxation reveals circumferential tension along the EGL.

At folding initiation the EGL has uniform proliferation, cell size, and cell shape.

The EGL has a uniform distribution of crossing fibers at folding initiation.

Live slice cutting and relaxation reveals radial tension between the EGL and The VZ.

During folding initiation, cells within the EGL are motile and undergo rearrangement.

Live slice imaging reveals fluid-like behavior in the EGL.

Live slice imaging reveals fluid-like behavior in the EGL.

A multiphase model with radial and circumferential constraints and liquid-fibrous EGL composition approximates evolution of cerebellar shape.

Differential expansion correlates with progressive subfolding of initial folds.

Uniform cell behaviors within a fluid-like EGL create differential expansion based folding approximated by a multiphase model with radial and circumferential tension.

Author details

Andrew K Lawton

Contribution

Competing interests

Tyler Engstrom

Contribution

Competing interests

Daniel Rohrbach

Contribution

Competing interests

Masaaki Omura

Contribution

Competing interests

Daniel H Turnbull

Contribution

Competing interests

Jonathan Mamou

Contribution

Competing interests

Teng Zhang

Contribution

Competing interests

J M Schwarz

Contribution

Competing interests

Alexandra L Joyner

Contribution

For correspondence

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Research organism

Further reading