Theoretical tool bridging cell polarities with development of robust morphologies

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

Abstract

Despite continual renewal and damages, a multicellular organism is able to maintain its complex morphology. How is this stability compatible with the complexity and diversity of living forms? Looking for answers at protein level may be limiting as diverging protein sequences can result in similar morphologies. Inspired by the progressive role of apical-basal and planar cell polarity in development, we propose that stability, complexity, and diversity are emergent properties in populations of proliferating polarized cells. We support our hypothesis by a theoretical approach, developed to effectively capture both types of polar cell adhesions. When applied to specific cases of development – gastrulation and the origins of folds and tubes – our theoretical tool suggests experimentally testable predictions pointing to the strength of polar adhesion, restricted directions of cell polarities, and the rate of cell proliferation to be major determinants of morphological diversity and stability.

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

eLife digest

Cells have the power to organise themselves to form complex and stable structures, whether it is to create a fully shaped baby from a single egg, or to allow adult salamanders to grow a new limb after losing a leg. This ability has been scrutinised at many different levels. For example, researchers have looked at the chemical messages exchanged by cells, or they have recorded the different shapes an embryo goes through during development. However, it is still difficult to reconcile the information from these approaches into a description that makes sense at multiple scales.

When an embryo develops, sheets of cells fold and unfold to create complex 3D shapes, like the tubes that make our lungs. Moulding sheets into tubes relies on interactions between cells that are not the same in all directions. In fact, two types of asymmetry (or polarity) guide these interactions. Apical-basal polarity runs across a sheet of cells, which means that the top surface of the sheet differs from the bottom. Planar cell polarity runs along the sheet and distinguishes one end from the other. For instance, apical-basal polarity marks the inner and outer surfaces of our skin, while planar cell polarity controls the direction in which our hair grows.

Nissen et al. set out to investigate how these polarities help cells in an embryo organise themselves to form complicated folds and tubes. To do this, simple mathematical representations of both apical-basal and planar cell polarities were designed. The representations were then combined to create computer simulations of groups of cells as these divide and interact with each other.

Simulations of ‘cells’ with only apical-basal polarity were able to generate different shapes in the ‘tissues’ produced, including many found in living organisms. External conditions, such as how cells were arranged to start with, determined the resulting shape. With both apical-basal and planar cell polarities, the simulations reproduced an important change that occurs during early development. They also replicated how the tubes that transport nutrients and oxygen form.

These results show that simple properties of individual cells, such as polarities, can produce different shapes in developing tissues and organs, without the need for a complicated overarching program. Abnormal changes in cell polarity are also associated with diseases such as cancer. The mathematical model developed by Nissen et al. could therefore be a useful tool to study these events.

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

Introduction

Multicellular organisms are amazing in their ability to maintain complex morphology in face of continuous cell renewal and damages. Adult salamander can regenerate entire limbs (Eguchi et al., 2011), and, during development, some regions can maintain patterning when moved to different parts of an embryo or if the size is varied (Lyons et al., 2012). Given the vast complexity and diversity of living shapes, how can we reconcile the robustness to perturbations with flexibility to diversify? While undoubtedly the end result is encoded in the DNA and protein networks, looking for an answer at this level is challenging. Examples of phenotypic plasticity (Libby and Rainey, 2011), convergent evolution, and contrasting rates of morphological and protein evolution (Cherry et al., 1979) show that morphological similarity may not couple to the protein sequence similarity (Gould, 1970). Inspired by the unfolding of morphological complexity in development, we propose that cellular polarity may be the key for reconciling complexity, robustness, and diversity of organismal morphologies.

During early development, the increase in morphological complexity coincides with the progressive polarization of cells – first apical-basal (AB) polarity and then planar cell polarity (PCP) (Müller and Bossinger, 2003; Roignot et al., 2013; Andrew and Ewald, 2010; Li and Bowerman, 2010). This theme is ubiquitous across vertebrates and invertebrates (Figure 1): starting from the single fertilized egg cell, first the morula formed by non-polarized cells turns into the blastula – a hollow sphere of cells with AB polarity. Then, as cells acquire additional PCP, primary head-tail axis forms and elongates during gastrulation and neurulation (Loh et al., 2016). Because of the optical transparency, these stages are particularly prominent in sea urchin. At the morula stage, a lumen in the center is formed and is gradually expanding as cells proliferate and rearrange into the hollow sphere. Next, during gastrulation, a group of cells invaginate and rearrange into a tube that narrows and elongates primarily by cell rearrangement and convergent extension movements (Martik and McClay, 2017). The tube then merges with the sphere at the side opposite to invagination, and as a result, the sphere transforms into a torus. Emerging data suggest that PCP drives both invagination and tube elongation (Nishimura et al., 2012; Croce et al., 2006; Long et al., 2015) – a recurring theme across species.

Figure 1 with 1 supplement see all

Download asset Open asset

Two symmetry-breaking events, gain of apical-basal (AB) polarity and planar cell polarity (PCP), on cellular level coincide with the appearance of a rich set of morphologies.

Starting from an aggregate of non-polarized cells (globular symmetry), individual cells can gain AB polarity and form one or multiple lumens (spherical symmetry). Additional, gain of PCP allows for tube formation (axial symmetry). Complex morphologies can be formed by combining cells with none, one, or two polarities. In Figure 1—figure supplement 1, we schematically illustrate how existing models capture different elements of development.

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

Mutations in PCP pathways produce shorter and wider tubes (Ochoa-Espinosa et al., 2012; Saburi et al., 2008; Kunimoto et al., 2017), somites (Song et al., 2010), and embryos (Gong et al., 2004). Formation and elongation of the tubes can proceed without cell division and cell death by cells rearranging along the tube’s axis, termed convergent extension (CE) (Andrew and Ewald, 2010; Tanimizu et al., 2009; Martik and McClay, 2017). While the importance of PCP in gastrulation and tubulogenesis is well established (Andrew and Ewald, 2010; Tanimizu et al., 2009; Martik and McClay, 2017; Kunimoto et al., 2017), it is unclear how polarity may control tube morphology.

The bulk-lumens-folds/tubes transition seen across animal species in early embryogenesis, is also a key feature of the later organ formation. The early stages of organogenesis in the liver, kidney, brain, gut, and pancreas are apparently so robust, that they can be recapitulated in vitro, allowing for advanced quantification and manipulation (Little, 2017). The case of pancreatic organoids is interesting as it illustrates an increase of morphological complexity from spheres to folds. Cells in in vitro pancreatic organoids first grow as a bulk and later acquire AB polarity and develop lumens. Depending on the growth conditions, organoids develop into a hollow sphere or acquire a complex folded shape (Greggio et al., 2013). It is currently unknown what drives the transition from sphere to folded state; the two possible hypotheses are rapid proliferation or physical pressure from growing into a stiff matrigel.

Is the apparent link between cellular polarity and morphological complexity accidental? Or, could it be that morphological transitions, stability, and diversity are emergent features in a population of proliferating polarized cells? If true, can we identify what drives the transition from lumens to folds and tubes? Why are these stable? Can we predict what controls fold depth, and tube length and width? To answer these questions, we lack a unified approach that could bridge polar interactions between single cells to the global features emerging on the scale of thousands of cells in 3D.

Starting with D’Arcy Thompson’s seminal contribution (Thompson, 1942), quantitative models aided in understanding specific morphogenetic events (Figure 1—figure supplement 1). Among these are invagination (Odell et al., 1981; Rauzi et al., 2015; Polyakov et al., 2014; Hočevar Brezavšček et al., 2012), primitive streak formation (Newman, 2008), convergent extension (Collinet et al., 2015; Belmonte et al., 2016), epithelial folding (Buske et al., 2012; Osterfield et al., 2013; Monier et al., 2015; Murisic et al., 2015), emergence of global PCP alignment from local cell–cell coupling (Amonlirdviman et al., 2005; Le Garrec et al., 2006; Burak and Shraiman, 2009), origins of tubulogenesis (Engelberg et al., 2008), and recently statistical properties of branching morphogenesis (Hannezo et al., 2017). However, they are often on either of the two ends of the spectra: those modeling single cells explicitly, often rely on vertex-based approaches and are limited to dozens of cells (Alt et al., 2017; Misra et al., 2016; Aigouy et al., 2010; Le Garrec et al., 2006). To capture the large features spanning thousands of cells, one typically turns to elastic models where AB polarity is implicit and epithelia is presented as a 2D elastic sheet (Hannezo et al., 2014; Etournay et al., 2015; Hufnagel et al., 2007; Nagai and Honda, 2009; Aliee et al., 2012; Nagai and Honda, 2001).

We developed a theoretical approach that, with only a few parameters, bridges cellular and organ scales by integrating both types of polarity. A main difference to earlier approaches is that a cell’s movement is coupled to how its AB polarity and PCP are oriented relative to each other and relative to neighbor cell polarities. In other words, in our approach, the adhesion strength between neighbor cells is modulated by the orientation of their polarities. We find that polarity enables complex shapes robust to noise but sensitive to changes in initial and boundary constraints, thus supporting that morphological stability and diversity are emergent properties of polarized cell populations. Lumens, folds, and stable tubes emerge as a result of energy minimization. We make testable predictions on morphological transitions in pancreatic organoids, tubulogenesis, and sea urchin gastrulation. Our approach illustrates the evolutionary flexibility in the regulatory proteins and networks, and suggests that despite differences in proteins between organisms, the same core principles may apply.

Model

There are three key elements that allow us to bridge the scale from cellular level to macroscopic stable morphologies.

(1) Cells are approximated by point particles

Cell–cell adhesion is modeled by repulsive and attractive forces acting between cell centers. This allows a substantially gain in computation time compared to vertex-based models where cell surface adhesion is explicitly considered (Alt et al., 2017). The potential for pairwise interaction between two interacting neighbors, i and j, separated by distance r_ij is

V_{i j} = e^{- r_{i j}} - S e^{- r_{i j} / β},

where the first term corresponds to repulsion, and the second term to attraction (see Figure 2A). For a pair of non-polar cells, the strength of attraction S = 1. β >1 is the parameter that sets how much longer the attraction range is compared to repulsion. We set β = 5 throughout the paper, but our results and conclusions are consistent for smaller β. The main results are also not sensitive to the exact choice of the potential, thus for example the higher power in the exponential,

V_{i j} = e^{- {(r_{i j})}^{4}} - S e^{- {(r_{i j} / β)}^{4}},

give qualitatively similar results (see Figure 2—figure supplement 1A). The potential energy of a cell is the sum of pairwise neighbor interactions

V_{i} = \sum_{j}^{} V_{i j} .

Figure 2 with 4 supplements see all

Download asset Open asset

Cells are modeled as interacting particles with a polarity-dependent potential.

(A) Potential between two interacting cells with apical-basal polarity (see Equation 6). Cells repulse when polarities are antiparallel (top/green part) and attract when they are parallel (orange/bottom part). (**B–C**) Two cells interact only if no other cells block the line of sight between them. (B) Cell i and j do not interact if ij’s midpoint (black dot) is inside of the Voronoi diagram for cell k (shaded in grey). (C) Cell i and j interact because cell k is further away than the distance ij/2 and ij’s midpoint therefore lie outside of cell k’s Voronoi diagram. In the related Figure 2—figure supplement 1A–D, we test the sensitivity of our model to the details of the potential and neighborhood assignments. In Figure 2—figure supplement 2, we relate changes in cell shapes to the model components, and in Figure 2—figure supplement 3 (Figure 2—video 1), we illustrate how altering polarity affects the dynamics of the systems with two and six cells.

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

(2) Cells interact with (a subset of Voronoi) neighbors

Interacting neighbors of cell i are selected from a subset of cells sharing a Voronoi surface. The subset is limited to the nearest neighbors j which are closest to the midpoint between i and j (Figure 2B–C). This constrain effectively corrects for the finite volume associated with point particles and assures that two cells will not interact if the line of sight between their centers is separated by a surface of a third cell. Without this constraint, the macroscopic morphologies collapse. However, our results are robust to replacing the line of sight constraint with full Voronoi and a cut-off distance for attraction force (Figure 2—figure supplement 1B).

(3) Cell–cell adhesion depends on the orientation of polarity

To capture directional adhesion, we set the strength of attraction, S, to be dependent on the relative orientation of the polarities in each of the cells. We assume that AB polarity and PCP are orthogonal, and that the polarities of one cell align with the polarities of its neighbor cells. Mathematically, we introduce unit vectors, ${\hat{p}}_{i}$ and ${\hat{p}}_{j}$ representing AB polarity, and ${\hat{q}}_{i}$ and ${\hat{q}}_{j}$ representing PCP for cell i and j, respectively. We set

S = λ_{1} S_{1} + λ_{2} S_{2} + λ_{3} S_{3}

where λ₁, λ₂, and λ₃ are the strengths of the different polarity terms. We require that

λ_{1} + λ_{2} + λ_{3} = 1,

to satisfy the constraint that perfectly aligned cells always have a steady state distance of 2 cell radii (for β = 5). To capture that in an epithelial sheet, AB polarities align parallel to each other and tight adherens junctions form in the plane perpendicular to AB polarity, we introduce the quadruple product

S_{1} = ({\hat{p}}_{i} \times {\hat{r}}_{i j}) \cdot ({\hat{p}}_{j} \times {\hat{r}}_{i j}) .

This makes two interacting cells with AB polarity maximally attracted (S₁ = 1) if the two apical sides are next to each other. On the other hand, if apical side of one cell is next to the basal side of another cell, the two cells will be maximally repulsing (S₁ = −1), see Figure 2A.

In case of planar polarization, we define

S_{2} = ({\hat{p}}_{i} \times {\hat{q}}_{i}) \cdot ({\hat{p}}_{j} \times {\hat{q}}_{j}),

which makes the attraction maximal if the PCP of two cells are parallel to each other and perpendicular to their AB polarities. In addition, we assume that similarly to AB polarity, two cells with PCP are maximally attracted if their PCP are parallel and cells have the same kind of pole (e.g. Vangl-enriched) next to each other,

S_{3} = ({\hat{q}}_{i} \times {\hat{r}}_{i j}) \cdot ({\hat{q}}_{j} \times {\hat{r}}_{i j}) .

We later show that this assumption makes neighbor exchange on a sheet possible and results in CE. However, unlike with tight junctions, preferred directional adhesion with PCP is not as well established. Cells adhere to each other by membrane proteins assembled in adherens junctions just below the apical surface. Both proteins regulating adherens junctions, for example Smash (Beati et al., 2018), as well as adherence proteins forming adherens junctions can be planar polarized, for example Bazooka, E-cadherins (Simões et al., 2010; Tamada et al., 2017; Levayer and Lecuit, 2013; Warrington et al., 2013; Aigouy and Le Bivic, 2016). These indirectly support our assumption of anisotropic, planar polarized adhesion.

The motion of the cells and their polarities are calculated assuming overdamped dynamics

\frac{{d \bar{r}}_{i}}{d t} = - \frac{d V_{i}}{d {\bar{r}}_{i}} + η,

\frac{{d \bar{p}}_{i}}{d t} = - \frac{d V_{i}}{d {\bar{p}}_{i}} + η,

and

\frac{{d \bar{q}}_{i}}{d t} = - \frac{d V_{i}}{d {\bar{q}}_{i}} + η,

where the ${\bar{p}}_{i}$ and ${\bar{q}}_{i}$ differentiation takes into account the rotation of polarity vectors, and η is a random uncorrelated Gaussian noise. In practice, we implemented the model in a MATLAB script (available in the Materials and methods section), where we use the Euler method. We perform the differentiation along the polarity by differentiating along all three cartesian coordinates (see Model details in the Materials and methods section). After each time step, we normalize the updated polarity vectors. The above differentiation does not include the change in partners when neighborhood changes. This is treated as a non-equilibrium step where potential energy can increase (Equation 3). Biologically, this is similar to cells spending biochemical energy as they rearrange their neighborhood.

The point particle approximation has been utilized earlier for modeling non-polar cell adhesion in early blastocyst (Krupinski et al., 2011), slug formation in amebae (Dallon and Othmer, 2004), and PCP organization in primitive streak formation (Newman, 2008). The main novelty of our approach is the dynamical coupling of cell positions and polarity orientations (Equation 6–8).

Model implications

One of the implications of the coupling between position and polarity is that in a sheet of cells, turning AB polarity in one cell will cause a force on its neighbors. In case of two cells (Figure 2—figure supplement 3 and Figure 2—video 1), the pair relaxes the imposed stress by rotating both the polarity and their positions. In biological terms, turning AB polarity in one cell (e.g. by apical constriction, illustrated in Figure 2—figure supplement 2) of an epithelial sheet will induce bending of the sheet as is the case with bottle cells in invagination.

The present formulation of PCP has several implications. First, we restrict the effects of PCP to directed (anisotropic) cell–cell adhesion and do not consider its other possible roles, in for example asymmetric cell differentiation, thus primarily focusing on its role on CE. Second, in our current formulation, AB polarity and PCP influence each other’s orientation on equal footing (Equation 7). PCP, however, is typically constrained to the apical plane and thus is expected not to influence the orientation of AB polarity. Disabling PCP’s effect on AB polarity (see Materials and methods) does not influence our main results on tube formation and gastrulation. However, the symmetry in polarities is appealing for its simplicity and is indirectly supported by the following experimental observations: First, cells can acquire PCP without AB polarity present (Baer et al., 2009; Zorn and Wells, 2009). Second, proteins required for AB polarity can be planar-polarized (Warrington et al., 2013; Aigouy and Le Bivic, 2016; Beati et al., 2018; Choi and Sokol, 2009; Dollar et al., 2005; Kaplan and Tolwinski, 2010). Third, changes in cell shapes during invagination (e.g. sliding of adherens junctions and formation of bottle cells) are regulated by PCP in neural tube closure (Ossipova et al., 2014; Nishimura et al., 2012; Kinoshita et al., 2008), gastrulation in C. elegans (Lee et al., 2006), sea urchin (Croce et al., 2006), and Xenopus (Choi and Sokol, 2009). These changes in cell shape effectively reorient AB polarity (Figure 2—figure supplement 2).

Results

We have recently introduced effective representation of AB polarity, and showed that it is sufficient for capturing spherical trophectoderm in the early blastocyst (Nissen et al., 2017). Expanding on that work, we here explore how AB polarity supports diverse yet stable and complex morphologies.

Stable complex shapes emerge from randomly polarized cell aggregates

Adult organismal shapes are stable over long time, maintaining sizes and relative positions of lumens and folds, despite continual local damages and cell renewal. To test if cellular polarity could enable such stability in time and to random local perturbations, we first performed a series of tests with AB polarized cells (Figure 3 and Figure 3—video 1).

Figure 3 with 3 supplements see all

Download asset Open asset

Development of 8000 cells from a compact aggregate starting at time 0.

(A) Cells are assigned random apical-basal polarity directions and attract each other through polar interactions (see Equation 6). (**A–D**) Cross-section of the system at different time points with red and blue marking two opposite sides of the polar cells. Cells closest to the viewer are marked red/blue, whereas cells furthest away are yellow/white. (E) Full system at the time point shown in (D). (F) Development of the number of neighbors per cell (red) and the energy per cell (blue), as defined by the potential between neighbor cells in Figure 2. Dark colors show the mean over all cells while light-shaded regions show the cell–cell variations. The yellow dot marks the energy for a hollow sphere with the same number of cells. See Figure 3—video 1 for full time series. In Figure 2—figure supplement 1E–G and Figure 3—figure supplement 1, we study how the final morphology depends on noise. In Figure 3—figure supplement 2, we show how the outer surface self-seals, and that the shape is maintained when cells divide.

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

When starting a bulk of cells with AB polarities pointing randomly, an initial rapid expansion (Figure 3A–C) stabilizes into a complex morphology of interconnected channels (Figure 3C–E). The shape remains unchanged for at least 10 times longer than the initial expansion (Figure 3C–E). The stability of the shape is illustrated by the time evolution of the average energy per cell (Figure 3F) that after an initial fast drop converges to a constant value. As expected, this value is higher than the energy of a hollow sphere (yellow dot in Figure 3F) – a configuration obtained if we start with radially, instead of random, polarized cells and preserve radial polarization at all times. The observed behavior is not sensitive to the shape of the potential (Figure 2—figure supplement 1A) but is sensitive to how the neighborhood is defined (Figure 2—figure supplement 1B–D). Rerunning the simulation in Figure 3 with different initial conditions results in a different stable shape (Figure 2—figure supplement 1E and Figure 3—figure supplement 1).

The macroscopic features of the shapes are robust to noise (Figure 2—figure supplement 1E–F and Figure 3—figure supplement 1). While the shapes emerging under high and low noise are not identical, the relative position and sizes of the majority of channels and lumens are preserved. The changes caused by noise stem from perturbations during initial expansion stage. If the same level of noise is applied after the system reached stable state, after time t = 10000, noise does not cause any major macroscopic changes (Figure 2—figure supplement 1G). The obtained shapes have self-sealing features, as an initial cut and unwrapping of a section of a surface refolds and seals back into the original morphology (Figure 3—figure supplement 2A–C). Furthermore, the shapes (Figure 3D–E) are also robust to overall growth (Figure 3—figure supplement 2D–F) retaining the same macroscopic features, just scaled to a larger size. Robustness to noise and cell proliferation further support the link between polarity and stability of morphologies, for example organ shapes, as they expand from infant to adult.

The final shapes are robust to noise but sensitive to initial and boundary conditions

The orientation of polarities in a subpopulation of cells may be set by the environment that the cells are embedded in, for example signaling molecules deposited into extracellular matrix can influence orientation of the AB polarity (Overeem et al., 2015) or signals from neighboring cells of a different type can orient PCP (Chu and Sokol, 2016). We will refer to these constraints as boundary conditions.

To investigate sensitivity to boundary conditions, we consider three cases where polarities are fixed at all times and point either radially out from the center of mass (Figure 4A), radially out from a central axis (Figure 4B), or pointing away from a central plane (Figure 4C). As anticipated, the difference in symmetries of boundary conditions results in a sphere, a cylinder, or two parallel planes (Figure 4—video 1). At the same time, in these symmetric cases, the differences in initial conditions but without imposed boundary conditions are not sufficient to generate different structures; they all converge to the nested ‘Russian doll’-like hollow spheres (Figure 4D). In development, this highlights the importance of the neighboring tissues for defining boundary conditions.

Figure 4 with 1 supplement see all

Download asset Open asset

Different morphologies can be obtained by varying boundary conditions (Figure 4—video 1).

(A) A hollow sphere emerges if polarities are fixed and initially point radially out from the center of mass. (B) A hollow tube is obtained if polarities point radially out from a central axis. (C) Two flat planes pointing in opposite directions are obtained if polarities point away from a central plane. (D) For all three initial conditions (**A–C**), if the polarities are allowed to change dynamically and the noise is high (η = 10⁰ compared to η = 10⁻¹ in A–C), the resulting shape consists of three nested ‘Russian doll’-like hollow spheres that will never merge due to opposing polarities. In contrast to the random initial condition in Figure 3, the initial conditions in (D) are symmetric.

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

Our results thus support the idea that polar adhesion enables stable and robust macroscopic shapes. The closest biological parallels would be the complex luminal morphologies emerging in reaggregation experiments on for example Hydra (Seybold et al., 2016) or in vitro culture of purjunkie brain cells (Muguruma et al., 2015). Together with our simulations, these experiments highlight how stable and complex morphologies can develop in non-proliferating populations from cell rearrangements alone.

Folding by pressure or rapid proliferation result in different fold-morphologies

Transitions from spheres to folded shapes are ubiquitous in development. Folds are an important part of in vivo organ development, and the composition of cell types in the folded organoids is closer to that in real organs (Greggio et al., 2013). To date, it is unclear what drives the transition from spheres to folded lumens. One possibility is that it is driven by the mechanical properties of the matrigel that effectively may place the growing organoid under pressure. Alternatively, data from 3D brain organoids suggests that the rapid cell proliferation leads to the emergence of surface folding (Li et al., 2017).

The simplicity of our tool allows to explore both of these scenarios. To model dividing cells, we pick a random cell from the entire population and introduce a new daughter cell with inherited polarity direction placed in a random location a half cell radius away from the mother cell. This event introduces dynamic perturbation by locally increasing cell density and requires some time to relax back to equilibrium. If proliferation is slow, and the time between two cell divisions anywhere in the system is longer than relaxation time of the whole system (the time it takes to reach equilibrium), the system approaches global equilibrium and will expand as a sphere. However, if proliferation is increased, the system will be pushed out of equilibrium and folds will emerge (Figure 5A, Materials and methods). In more quantitative terms, our forces are such that a single cell can move up to 0.2 cell diameter per time unit. For cells that divide every 1000 time units, the transition to non-equilibrium buckling happens when the system has grown to about 5000 cells (Figure 5—video 1).

Figure 5 with 2 supplements see all

Download asset Open asset

The number of complex folds in a growing organoid depends on the generation time and the pressure from the surrounding medium (Figure 5—video 1).

(A) Number of local minima as a function of 1/(generation time), *t_G*⁻¹. In silico organoids grow from 200 cells up to 8000, 12,000, or 16,000 cells with different generation times and no outer pressure. (B) Number of local minima as a function of pressure, P. In silico organoids grow to the same size with the same 1/(generation time), *t_G*⁻¹ = 1.4⋅10⁻⁴ but different outer pressure. The images illustrate the 16,000 cells stage. Blue dots mark the average, while light shaded regions show the SEM based on triplicates. See also Figure 5—figure supplement 1 for additional measurements on the differences between rapid growth and pressure.

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

As cells divide faster, our simulations predict a transition from a smooth spherical shell to an increasingly folded structure with multiple pronounced folds, in line with the observation of brain organoids proliferating at different rates (Li et al., 2017). In comparison with the model for cortical convolutions by Tallinen et al., 2016 in which folding is a result of expanding cortical sheet adhered to the non-expanding white matter core, our mechanism does not require a bulk core. Instead folds emerge in a fast-expanding sheet when the growth is faster than the global relaxation to dynamical equilibrium.

While we find that the external pressure is not necessary for folding, pressure alone can also drive folding (Figure 5B, Materials and methods). However, this scenario contradicts the observation that pancreatic organoids can grow as spheres or folded morphologies in gels with the same stiffness but different media composition (Greggio et al., 2013).

In principle, both scenarios may contribute to folding, but visually the fold morphologies are different. To differentiate between the two, we have quantified the final folded structures in terms of their local minima (Figure 5, Materials and methods, see also Figure 5—video 1). Our simulations predict that in the pressure-driven case, the number of local minima will reach an upper limit as organoids increase in size (Figure 5B). In the case of out-of-equilibrium proliferation, new folds can continue forming as organoids grow (Figure 5A). Increased proliferation causes more and shallower folds. These folds are different than obtained with pressure which causes fewer but deeper minima. Quantitatively, both the depth and the horizontal extension of the folds are about double as large with pressure than with growth-induced folding (Figure 5—figure supplement 1).

PCP enables convergent extension and robust tubulogenesis

Despite the numerous evidences supporting the role of PCP in tubulogenesis, it remains unresolved whether oriented cell division or the extent of CE controls tube length and width (Karner et al., 2009; Carroll and Yu, 2012). It is also debated if it is important for the tubes to maintain regular shape, or if it is only important for tube initiation and growth (Kunimoto et al., 2017).

The simplicity of our approach allows us to address these questions by introducing cell–cell interactions through PCP. This term favors front-rear cell alignment in the interaction potential with only two additional parameters: the strength of the orientational constraint of AB polarity with respect to PCP, λ₂, and the strength of PCP, λ₃ (see Equations 4–8). For simplicity, we focus on the stability (ability to maintain regular diameter over time) and tube morphogenesis in systems without cell division.

Inducing PCP in a spherical lumen leads to two significant events. First, independent of initial orientation, after some transient time PCPs becomes globally ordered and point in direction parallel to an emerging equator that self-organizes around the sphere (inset in Figure 6). This arrangement has the lowest energy. Second, cells start intercalating along the axis perpendicular to PCP orientation, gradually elongating the lumen (Figure 6—video 1). During intercalations, cells exchange their neighbors through T1-like transitions as reported experimentally, Figure 6—figure supplement 3 (Nishimura et al., 2012; Sanchez-Corrales et al., 2018). The intercalations along the axis continue until the force balance between AB polarity and PCP is restored at a new equilibrium. Thus, our model predicts that the strength of PCP (λ₃) relative to AB polarity (λ₁) determines the width and the length of the tube (Figure 6). We obtain similar results if we constrain PCP to always remain in the apical plane and thus does not allow PCP to reorient AB polarity (Figure 6—figure supplement 1, Materials and methods). Note, that this result is very different in nature from the tube presented in Figure 4B as both AB polarity and PCP can now reorient in each cell at any time.

Figure 6 with 4 supplements see all

Download asset Open asset

The length and width of tubes are set by the strength of planar cell polarity (PCP, λ₃).

For each value of λ₃, we initialize 1000 cells on a hollow sphere with PCP whirling around an internal axis (PCP orientation marked by cyan arrows in the top-left inset). Semi-major axis (dark blue) and semi-minor axis (light blue) are measured at the final stage (Materials and methods). Images show the final state. Throughout the figure, λ₂ = 0.5 and λ₁ = 1 - λ₂ - λ₃. The animated evolution from sphere to tube is shown in Figure 6—video 1. See also Figure 6—figure supplement 1 where we show that tubes also form when we disable the direct influence of PCP on apical-basal polarity, and Figure 6—figure supplement 2 where we vary the degree of PCP along the axis of the tube. In Figure 6—figure supplement 3, we show that cell intercalations result in experimentally reported T1 neighbor exchanges during convergent extension.

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

These results support the observations that stable tubes can emerge without cell proliferation. In addition, when first the tube is formed, loss of PCP does not lead to cyst formation as recently shown by Kunimoto et al. (2017). However, localized cysts could result if the lumen is initialized with varying strength of PCP along the axis (Figure 6—figure supplement 2).

Two polarities are sufficient to explain major features of sea urchin gastrulation

Currently invagination in neurulation and gastrulation is understood and quantitatively modeled as a process driven by changes in cell shapes or the mechanical properties of cells with AB polarity (Rauzi et al., 2013; Tamulonis et al., 2011; Misra et al., 2016; Hočevar Brezavšček et al., 2012). This process is often assumed to be driven by apical constriction and decoupled from the eventual tube formation and elongation. However, emerging data suggests that PCP drives both invagination and tube elongation (Nishimura et al., 2012; Croce et al., 2006; Long et al., 2015) likely because apical constriction is controlled by PCP (Ossipova et al., 2014; Nishimura et al., 2012; Croce et al., 2006) and cell intercalations, similar to those in CE, contribute to invagination (Nishimura et al., 2012; Rembold et al., 2006; Sanchez-Corrales et al., 2018).

To probe the limits of our approach, we investigated if AB polarity, PCP, and boundary conditions reminiscent of posterior organizer (Loh et al., 2016) are sufficient to recapitulate the main stages of sea urchin gastrulation: invagination, tube formation, elongation (by CE), and finally merging of the gastrula tube with the pole opposite to invagination site.

The current understanding is that sea urchin gastrulation consists of primary invagination, driven by swelling of the inner layer of extracellular matrix beneath invaginating cells (Lane et al., 1993), and formation of a ring of bottle cells due to apical constriction (Kimberly and Hardin, 1998), and secondary invagination where tube elongates due to CE (Lyons et al., 2012). PCP is necessary for both invagination, possibly through its effect on apical constriction in bottle cells (Nishimura et al., 2012), and tube extension (Croce et al., 2006).

Motivated by these observations, we set boundary condition such that PCP of the invaginating cells are oriented around the anterior-posterior (top-bottom) axis, and are always in the apical plane. This constrain on PCP orientation allows for CE. While this particular configuration is not documented, it is consistent with observed effects of WNT orienting PCP within the apical plane (Humphries and Mlodzik, 2018). Second, we simulate the combined effect of bottle cells (Figure 2—figure supplement 2) and bending by swelling of the extracellular matrix by applying an external force, F, on AB polarity (see Materials and methods). This force gradually reorients AB polarity away from the anterior-posterior axis, thus leading to bending of the epithelial sheet. The effect is maximal for cells closest to the anterior-posterior axis. This external force is a phenomenological description aiming at capturing the observed effects of how change in AB polarity results in tissue bending and does not aim at capturing mechanisms driving the reorientation of AB polarity.

As a result, cells start to rearrange, the bottom flattens (Figure 7B) and bends inward (Figure 7C). Subsequently, the CE-driven by PCP causes the invaginated cells to rearrange, tube elongates, and merges with the top of the sphere (Figure 7F and Figure 7—video 1). In line with experimental observations, the tube elongates due to cells moving in the tube (Martik and McClay, 2017).

Figure 7 with 2 supplements see all

Download asset Open asset

External constraints on apical-basal (AB) polarity and planar cell polarity (PCP) can initiate invagination and drive gastrulation in sea urchin.

(A) The lower third of the cells in a blastula with AB polarity (apical is blue–white, basal is red–orange) pointing radially out acquire PCP (cyan–green) in apical plane pointing around the anterior-posterior (top-bottom) axis (as in the inset to Figure 6). (B) Flattening of the blastula and (C) invagination occur due to external force reorienting AB polarity (Materials and methods). (**D–E**) Tube elongation is due to PCP-driven convergent extension and (F) merging with the top of the blastula happens when the tube approaches the top. Throughout the simulation, λ₁ = 0.5, λ₂ = 0.4, and λ₃ = 0.1 for the lower cells while the top cells have λ₁ = 1 and λ₂ = λ₃ = 0. For full time dynamics see Figure 7—video 1. In Figure 7—figure supplement 1, we consider alternative scenarios of sea urchin gastrulation and and neurulation.

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

The contribution of extracellular matrix swelling is specific to sea urchin invagination. As bottle cells alone are sufficient to drive invagination in other systems, we have tested a cell-intrinsic scenario where AB polarities of bottle cell neighbors prefer to be tilted towards each other (by e.g. modifying the potential in Equation 6 to capture AB polarities as in Figure 2—figure supplement 2, unpublished results). This also resulted in successful invagination.

Several observations suggest that radial intercalation movements towards the center of invagination may drive tissue bending (Sanchez-Corrales et al., 2018; Panousopoulou and Green, 2016; Rembold et al., 2006; Chung et al., 2017). When tested in our model, PCP alone (omitting external or cell-intrinsic reorientation of AB polarity) results in exvagination and tube elongating outside of the sphere (this is more energetically favorable than extending inwards). A similar, exogastrulating, phenotype was observed in a PCP mutant (Long et al., 2015). While the similarity may be accidental, it is possible that this mutation abrogates PCP-driven apical constriction (Ossipova et al., 2014; Nishimura et al., 2012) and since cellular apices face outside, abrogating apical constriction eliminates the bias in direction of tube formation. We have also tested the consequences of apical constriction (reorienting AB polarity) in the model without PCP and while this was sufficient for invagination, the cavity remained spherically symmetric and failed to form a tube.

Discussion

Despite the stunning diversity and complexity of morphologies, the same concepts seem to emerge across organismal development. One of them is the link between local, cellular, and global, organs/whole organism, symmetry breaking. We know, from experimental and, to a lesser degree, theoretical work that cellular polarity is essential for forming axis, complex folded sheets, and interconnected tubes (Figure 1—figure supplement 1). What we do not know is why are these shapes so stable, and where do the differences between species and organs come from. To understand the differences, we typically compare genes or gene regulatory networks, thus limiting our understanding to analog processes within a few related species.

Phenomenological description bridges cell polarity to macroscopic morphologies

To address the origins of morphological diversity and stability across species and organs, we focused on a phenomenological description of polarized cell–cell interactions. This allowed us to bridge local single-cell symmetry breaking events to global changes in morphologies spanning tens of thousands of interacting cells. With this tool at hand, we find that with only a few parameters, we can recapitulate the two global symmetry breaking events: formation of epithelial sheets and folds by cells with AB polarity, and emergence of global axial symmetry (tubes) among cells with PCP.

Remarkably, our results show that interactions among AB polarized cells lead to stable morphologies, that after initial relaxation remain indefinitely in their final configuration. The morphologies are robust to noise, growth, and local damage (Figure 3—figure supplement 2). These results may explain how organs and embryos preserve their architecture while growing. Polar cell–cell interactions not only provide clue to the morphological stability, but also point to a simple explanation to the origin of the diversity. We find that the exact morphological details are defined by initial conditions, for example initial positions and orientation of polarities, and boundary conditions, for example polarities restricted to certain direction for a fraction of the cells. It is thus tempting to speculate that diverse shapes do not require multiple interacting morphogen gradients, but can be a result of differences in initial and/or boundary conditions: as for example presence of yolk cells at start and boundary constraints by vitelline membrane (Schierenberg and Junkersdorf, 1992; Wu et al., 2010).

The diversity of shapes and forms is further enriched by a second symmetry breaking event, PCP, oriented perpendicular to AB polarity. Within our phenomenological framework, addition of PCP component is simple, and requires only two additional parameters: one favoring perpendicular orientation of AB polarity and PCP within a cell, and another, favoring parallel PCP alignments between neighbor cells. These constraints are the coarse-grained representation of the well-established experimental and computational results on intracellular symmetry breaking events and global ordering of planar polarities mediated by cell–cell coupling (Le Garrec et al., 2006; Amonlirdviman et al., 2005; Wang et al., 2006). The first constraint allowed formation of axial symmetry and in combination with AB polarity, stable tubes, with length and diameter remaining constant with time. The second constraint resulted in cell rearrangements and intercalations consistent with the cell-autonomous CE typically associated with PCP. The patterns of neighbor exchanges during cell-rearrangements are in line with the ubiquitous T1 exchanges through formation and resolution of four cell vertex, Figure 6—figure supplement 3 (Nishimura et al., 2012; Sanchez-Corrales et al., 2018). The mechanism of the CE in our model is in line with the results by ‘filopodia tension model’ where elongated structures of many cells emerge from local cell–cell interactions in a direction defined by PCP (Belmonte et al., 2016). The presented formulations of our model captures only some of the known events contributing to CE and does not include PCP-driven changes in cell shape mediated by for example apical constriction (Nishimura et al., 2012) or contribution of external forces (Lye and Sanson, 2011).

Combining AB polarization and a local induction of PCP in a subpopulation of cells was sufficient to obtain main stages of sea urchin gastrulation: invagination, tube formation, and elongation through CE as well as merging of the tube with the animal pole at the top of the blastula. It is important to notice that the model for gastrulation uses an external force (Equation 25) acting on AB polarity, and the model for neurulation uses an imposed external constraint on PCP. This is not fully satisfying, and suggests an extension of the model to capture tissue bending induced by changes in cell shapes at the edge of the region that invaginates.

The existing in silico models treat invagination and CE-driven tube elongation as independent processes (Figure 1—figure supplement 1). Recent data, however, suggests that multiple mechanisms (intracellular apical constriction, intercellular directed cell division and cell intercalations, and supracellular actomyosin cables) act simultaneously and contribute to both invagination and tubulogenesis (Chung et al., 2017; Nishimura et al., 2012; Ossipova et al., 2014). Within our approach apical constriction (modeled as reorientation of AB polarity) and CE can act in parallel (Figure 7). It will be interesting to parallel recent experimental work (Chung et al., 2017) and computationally investigate how a combination of intra-, inter- and supracellular mechanisms contribute to the robustness of tubulogenesis, and to what extent the model can capture the range of observed phenotypes.

It has been proposed that the above mechanisms may all be coordinated by PCP (Nishimura et al., 2012). Besides the reported molecular links, a simple logic suggests that these mechanisms cannot be isotropic as in this case the initial bending will result in spherical structures. Thus, apical constriction, cell intercalations, and actomyosin cables have to be anisotropic (planar polarized) in directions consistent with the eventual tube orientation. This anisotropy is reported for both ‘wrapping’ tubes forming parallel to the epithelial plane, for example neurulation (Nishimura et al., 2012), and ‘budding’ tubes forming orthogonally to the epithelial plane, for example salivary glands (Chung et al., 2017; Sanchez-Corrales et al., 2018).

As organizing signals such as WNT can induce and orient PCP (Chu and Sokol, 2016) within the apical plane, we asked if it is in principle possible to design PCP constraints (not limited to the apical plane) that would result in ‘wrapping’ and ‘budding’. First, pointing PCP out of the plane was sufficient for a sheet of cells to bend (Figure 7—figure supplement 1). This is because in our formulation, PCP can drive reorientation of AB polarity and that in its turn is able to bend the sheet (Figure 7). Both ‘budding’ and ‘wrapping’ were qualitatively captured by the model when the axial and radial anisotropy were set by constraining orientation of PCP for cells within a circle (‘budding’ in sea urchin) and two stripes of cells (mimicking hinge points in neurulation ‘wrapping’). While CE was needed for proper tube forming in sea urchin example, the tube formed without CE in neurulation (Figure 7—figure supplement 1).

Thus, the simulations suggest that ‘wrapping’ in neurulation and gastrulation in Drosophila (Figure 7—figure supplement 1) vs. ‘budding’ in sea urchin and organogenesis (Andrew and Ewald, 2010; Zegers, 2014) may be outcomes of different constraints imposed on PCP. While it is intriguing to speculate that PCP may be oriented out of epithelial plane directly by organizing signals, this may also be an indirect effect of a sequence of intermediate steps. As the organizing signals not only induce and orient PCP but also drive apical constriction (effectively reorienting AB polarity) the PCP may be gradually oriented out of the original plane of epithelium by the following sequence of events: PCP → apical constriction → tilt in AB polarity → tilt in apical plane → PCP out of original epithelial plane. This less precise but simpler interpretation highlights the fact that PCP may drive many of the alternative mechanisms of tubulogenesis and shifts the focus from the differences in mechanisms driving tubulogenesis to the differences in boundary conditions – a set of constraints imposed on cell polarities by neighboring tissue (e.g. notochord in neurulation and organizer in gastrulation).

Testable predictions

In addition to our conceptual findings, we propose three testable predictions. First, we predict that two potential mechanisms behind the emergence of folds in pancreatic organoids – matrigel resistance and rapid, out-of-equilibrium, cell proliferation – will result in distinct morphologies. Our results suggest that in case of rapid proliferation, the growing structure will develop many shallow folds close to the surface which later tend to deepen. In contrast, external pressure causes fewer but deeper and longer folds (Figure 5—figure supplement 1). And further, as organoids grow in size, the number of folds will reach an upper limit when under pressure, however, in case of rapid proliferation, the number of folds will keep growing (Figure 5). Visual inspection of published morphologies seems to support the out-of-equilibrium growth (Greggio et al., 2013; Li et al., 2017). To assess if the growth is out-of-equilibrium in 3D organoids, one can quantify the distributions of cell shapes (Cerruti et al., 2013). Our model thus predicts that quantitative counting of folds and measurements of the fold depth and length relative to the size of the growing structure may discriminate between the alternative hypothesis. Quantification of the folds can be done in in vitro organoids by either phase or confocal fluorescence microscopy of whole-mount immunostained samples (Greggio et al., 2013; Li et al., 2017). The fold depth and length can be quantified with the same approach as applied to simulated shapes (Materials and methods) in binarized images of the 3D organoid surfaces. As oganoids cannot be cultured without gel supporting 3D growth, it will be necessary to vary both gel stiffness and generation time to uncouple their respective contribution to the folding. The work by Little, 2017 shows that in brain organoids generation time can be both slowed down and speeded up by either genetic manipulation or by adding small molecule inhibitors of pathways regulating cell proliferation. Unfortunately, changing stiffness of matrigel also changes its biochemical composition and may affect cell proliferation and differentiation. One will have to turn to synthetic hydrogels, where it is now possible to uncouple mechanical and biochemical clues (Gjorevski et al., 2016). To illustrate possible applications of our approach, we have only focused on two out of several possible mechanisms that may contribute to folding. The other likely alternative is that folding may result from differences in biomechanical interactions or generation times characteristic to the different cell types. These alternative scenarios are straightforward to consider in our model and will be an exciting venue to explore when more quantitative data on differences in organoid morphologies is available.

Our second prediction is that in case of tubes formed by non-proliferating cells, the length and width of the tubes are controlled by the relative strength of AB polarity and PCP. This result calls for quantification of adhesion proteins along the AB polarity and PCP axes. In PCP mutants with shorter and wider tubes, one would expect less planar polarization in adherens junctions and actomyosin, for example larger spread compared to wild type in their orientation quantified relative to the tube axis (Nishimura et al., 2012). Alternatively, the balance between AB polarity and PCP can be altered by weakening AB polarity, for example mutating tight junction proteins should result in longer tubes. A similar phenotype has already been reported for Drosophila tracheal tube (Laprise et al., 2010). With the recent advancements in in vitro systems of tubulogenesis, allowing for easy genetic manipulations and more amenable for quantitative imaging, it may in principle be possible to relate the extent of planar anisotropy in PCP mutants and strength of AB adhesion in tight junction mutants with tube length and diameter. The existing coupling between PCP and AB polarity may, however, make it challenging to tweak one polarity at a time.

Our third, and probably most challenging to test, prediction is on the conditions differentiating between tubes forming perpendicular (e.g. sea urchin gastrulation) or parallel (as in Drosophila gastrulation or neurulation) to the plane of epithelium. We predict that the outcome will be defined by the orientation of PCP in the invaginating region and the geometry of the boundaries (circular for budding and axial for wrapping) set by for example WNT organizing signals. Recent development in imaging localization of PCP complexes in single cells (Wu et al., 2013; Chu and Sokol, 2016; Minegishi et al., 2017; Habib et al., 2013) allows monitoring localization of PCP complexes, and thus PCP orientation, in individual cells. By placing WNT-soaked (Habib et al., 2013) beads or WNT-secreting cells (Chu and Sokol, 2016) one can vary PCP orientation in the cells at the epithelial boundary facing WNT and test for the direction of the epithelial bending and possibly tube formation.

Our approach is by purpose phenomenological and by its nature cannot make predictions about specific molecular details. In all cases, we do not see our simulations as finalized predictions, but rather as pointing in the most promising direction for further exploration of these complex developmental processes. Our setup easily allows for changes as we learn more. The proposed tool should be used in close collaboration with gained experimental knowledge on initial conditions, cell generation times and differentiation processes where polarities play a central role.

Our results open for a series of biological generalizations both in development and diseases. On one hand, we now may be able to explain and unify the apparently very distinct morphological transitions during gastrulation in flies, frogs, fish, mice, and humans by accounting for different initial and boundary conditions. Our model suggests how a moderate change in expression of polarities during some critical evolutionary stages could lead to widely different final morphologies. Thereby, development driven by cell–cell polarity interactions could provide major morphological transitions from local and transient modulations in polarity.

On the other hand, it becomes possible to think of gastrulation, neurulation, tubulogenesis, and organogenesis as the same class of phenomena, where the orientation of the tube is guided by local organizers, and lengths/widths of the tubes are determined by the relative strength of AB polarity and PCP. At the same time, there is an emerging view that wound healing and cancer are local perturbations – for example local loss of cells, dysregulation of cell polarities (Martin-Belmonte and Perez-Moreno, 2012), proliferation, or autonomously induced organizing signals – of otherwise conceptually the same developmental processes (Humphries and Mlodzik, 2018). The power of our model is that it allows to address these hypotheses through predictive models for the dynamics of many cells that interact through combinations of AB polarity and PCP.

Materials and methods

Throughout the paper, we use the Euler method to integrate the ordinary differential equations stated in the Model section with dt set to 0.1 or 0.2. Lower dt values will give qualitatively similar results but with increased simulation time. Higher dt values, will result in a collapse of the presented morphologies. The time unit is arbitrary, and the same throughout the paper.

The noise parameter η = 10⁻⁴ where nothing else is stated. Lower noise values will give more smooth simulations, while η on the order of 10⁰ will result in collapsing shapes (see also Figure 2—figure supplement 1E–F).

Model details

Request a detailed protocol

In our model, we use the following potential to describe the pairwise interaction between cells

V_{i j} = e^{- r_{i j}} - S e^{- r_{i j} / β},

where r_ij is the center–center distance between cell i and cell j, and S is the polarity factor

S = λ_{1} S_{1} + λ_{2} S_{2} + λ_{3} S_{3}

Here, λ₁, λ₂, and λ₃ are the strengths of the respective polarity terms which are given as

S_{1} = ({\hat{p}}_{i} \times {\hat{r}}_{i j}) \cdot ({\hat{p}}_{j} \times {\hat{r}}_{i j}),

S_{2} = ({\hat{p}}_{i} \times {\hat{q}}_{i}) \cdot ({\hat{p}}_{j} \times {\hat{q}}_{j}),

S_{3} = ({\hat{q}}_{i} \times {\hat{r}}_{i j}) \cdot ({\hat{q}}_{j} \times {\hat{r}}_{i j}) .

The unit vectors ${\hat{p}}_{i}, {\hat{p}}_{j}$ and ${\hat{q}}_{i}, {\hat{q}}_{j}$ represent the apical-basal polarity and planar cell polarity of cell i and j. Throughout the paper, β is a constant which we set to 5. In order to use the Euler method, we need the gradient of V_ij differentiated with respect to position, ${\bar{r}}_{i}$ , and the two polarities, ${\bar{p}}_{i}$ and ${\bar{q}}_{i}$ :

\frac{d V_{i j}}{{d \bar{r}}_{i}} = e^{- r_{i j} / β} {γ {\hat{r}}_{i j} - \frac{λ_{1}}{r_{i j}} [({\hat{r}}_{i j} \cdot {\hat{p}}_{j}) {\hat{p}}_{i} + ({\hat{r}}_{i j} \cdot {\hat{p}}_{i}) {\hat{p}}_{j}] - \frac{λ_{3}}{r_{i j}} [({\hat{r}}_{i j} \cdot {\hat{q}}_{j}) {\hat{q}}_{i} + ({\hat{r}}_{i j} \cdot {\hat{q}}_{i}) {\hat{q}}_{j}]},

\frac{d V_{i j}}{{d \bar{p}}_{i}} = e^{- r_{i j} / β} {λ_{1} [S_{1} {\hat{p}}_{i} - {\hat{p}}_{j} + ({\hat{r}}_{i j} \cdot {\hat{p}}_{j}) {\hat{r}}_{i j}] + λ_{2} [S_{2} {\hat{p}}_{i} - ({\hat{q}}_{i} \cdot {\hat{q}}_{j}) {\hat{p}}_{j} + ({\hat{q}}_{i} \cdot {\hat{p}}_{j}) {\hat{q}}_{j}]},

\frac{d V_{i j}}{d {\bar{q}}_{i}} = e^{- r_{i j} / β} {λ_{2} [S_{2} {\hat{q}}_{i} - ({\hat{p}}_{i} \cdot {\hat{p}}_{j}) {\hat{q}}_{j} + ({\hat{p}}_{i} \cdot {\hat{q}}_{j}) {\hat{p}}_{j}] + λ_{3} [S_{3} {\hat{q}}_{i} - {\hat{q}}_{j} + ({\hat{r}}_{i j} \cdot {\hat{q}}_{j}) {\hat{r}}_{i j}]} .

In order to derive Equations 17, 18, and 19, we have used the following:

\frac{d}{d {\bar{r}}_{i}} e^{- r_{i j} / β} = \frac{1}{β} e^{- r_{i j} / β} {\hat{r}}_{i j},

γ = e^{- r_{i j} (β - 1) / β} - \frac{S}{β} + \frac{2}{r_{i j}} [λ_{1} ({\hat{r}}_{i j} \cdot {\hat{p}}_{i}) ({\hat{r}}_{i j} \cdot {\hat{p}}_{j}) + λ_{3} ({\hat{r}}_{i j} \cdot {\hat{q}}_{i}) ({\hat{r}}_{i j} \cdot {\hat{q}}_{j})],

\frac{d S_{1}}{d {\bar{r}}_{i}} = \frac{1}{r_{i j}} [({\hat{r}}_{i j} \cdot {\hat{p}}_{i}) {\hat{p}}_{j} + ({\hat{r}}_{i j} \cdot {\hat{p}}_{j}) {\hat{p}}_{i}] - 2 [({\hat{r}}_{i j} \cdot {\hat{p}}_{i}) ({\hat{r}}_{i j} \cdot {\hat{p}}_{j}) {\hat{r}}_{i j}],

\frac{d S_{1}}{d {\bar{p}}_{i}} = \frac{1}{p_{i}} [{\hat{p}}_{j} - ({\hat{r}}_{i j} \cdot {\hat{p}}_{j}) {\hat{r}}_{i j} - S_{1} {\hat{p}}_{i}],

where p_i is the length of the polarity of cell i which is equal to one at all times.

Generation time in growing organoids

Request a detailed protocol

In Figure 5, the number of cells, N, at a given time, t, is N = 200 exp(ln(2) t/t_G) where t_G is the generation time. In these simulations, the AB potential (Equation 6) between cells is set to zero when the angle between their polarity is larger than π/2.

Modeling resistance from matrigel

Request a detailed protocol

In Figure 5, we model resistance from the matrigel by imposing a surface force pointing toward the center of mass. The potential of the pressure in the growth medium is given by

V_{M} = - \frac{P r^{2}}{2 r_{m a x}}

where P is the stiffness of the medium, r is distance from the center of mass, and r_max is the distance to the cell that is the furthest away from the center of mass. The resulting force will be constant in time at the periphery. Thus, all cells on a growing sphere will be exposed to a force of equal size. However, cells that end up deep inside a folded morphology will experience weaker resistance.

Quantification of the local minima

Request a detailed protocol

In Figure 5, the number of local minima is defined as the number of cells that do not have any neighbor cells that are closer to the center of mass than themselves, and at the same time have an average angle between their AB polarity and their neighbor cells displacement vector that is less than π/2.

Measuring the tube length and width

Request a detailed protocol

In Figure 6, the semi-minor and semi-major axes correspond to the half-width and half-length of the tubes, respectively. As cells on opposite sides of a tube have AB polarity pointing in opposite directions, we approximate the semi-major and semi-minor axes, by finding the half of the maximum and minimum distance between two cells with AB polarity pointing in opposite directions.

Modeling cells with different polarities

Request a detailed protocol

In the gastrulation simulation (Figure 7), each cell is assigned a specific value of polarity strengths (λ_1,i, λ_2,i, and λ_3,i). We define the mutual interaction strength between a pair, i and j, of cells in Equation 4 with different polarity strengths by setting λ₁ = mean(λ_1,i, λ_1,j), and λ₂ = mean(λ_2,i, λ_2,j) as well as λ₃ = mean(λ_3,i, λ_3,j). This choice makes sure that two neighbor cells interact with a force with equal magnitude but opposite sign. Furthermore, it makes sure that λ₁ + λ₂ + λ₃ = 1 holds for all cells.

Disabling PCP’s effect on AB polarity

Request a detailed protocol

In the model described by Equations 6–9, AB polarity and PCP can influence each other’s orientation. To constrain PCP to the apical plane and thus disable its influence on AB polarity, we set λ₂ = 0 when updating AB polarity in time (during the numerical integration of Equation 10).

Modeling bottle cells and apical constriction by an external force

Request a detailed protocol

The invagination in gastrulation is implemented by adding an external force, F, that act on the AB polarity in addition to our usual intrinsic forces from Equation 10

F = - k \hat{r} e^{(- x^{2} - y^{2}) / σ^{2}},

The two parameters are k which is the strength of the force (in Figure 7, k = 0.02), and σ which defines the decline of the gaussian force (in Figure 7, σ = 10). $\hat{r}$ is the unit vector pointing from x = y = 0 to the cells’ position, and x and y are the respective coordinates of the cells. This force is applied on the AB polarity and bends the orientation of the polarity away from a z-axis (the anterior-posterior axis).

MATLAB script

Request a detailed protocol

MATLAB script to generate and visualize data. MATLAB R2016b or newer is required together with the Statistics and Machine Learning Toolbox. In addition, the Parallel Computing Toolbox is required if the PAR parameter in the basic script is set to 1. The input folder contains initial conditions for three standardized systems. Bulk systems have neither apical-basal (AB) polarity nor planar cell polarity. Plane and shell systems have only AB polarity. ‘N’ in the file names gives the number of cells in the system. The initial polarity directions can be modified on line 4–5 in the basic.m file. Inside this file, it is also possible to set the degree of noise (η), the size of the time steps (dt), and the relation between the polarity strengths (λ₁, λ₂, and λ₃). The parameter inc is used to speed up the simulations by only applying the neighborhood function to the nearest 100 neighbors. Generated data is saved in the output folder, and the visualization script is in a separate folder.

Data availability

All data generated or analyzed during this study are included in the manuscript and supporting files. MATLAB code to reproduce or generate new data is added as a supplementary zip file together with a MATLAB script to visualize the data.

References

1. Aigouy B
2. Farhadifar R
3. Staple DB
4. Sagner A
5. Röper JC
6. Jülicher F
7. Eaton S
(2010) Cell flow reorients the axis of planar polarity in the wing epithelium of Drosophila
Cell 142:773–786.

https://doi.org/10.1016/j.cell.2010.07.042
- PubMed
- Google Scholar
1. Aigouy B
2. Le Bivic A
(2016) The PCP pathway regulates Baz planar distribution in epithelial cells
Scientific Reports 6:33420.

https://doi.org/10.1038/srep33420
- PubMed
- Google Scholar
1. Aliee M
2. Röper JC
3. Landsberg KP
4. Pentzold C
5. Widmann TJ
6. Jülicher F
7. Dahmann C
(2012) Physical mechanisms shaping the Drosophila dorsoventral compartment boundary
Current Biology 22:967–976.

https://doi.org/10.1016/j.cub.2012.03.070
- PubMed
- Google Scholar
(2017) Vertex models: from cell mechanics to tissue morphogenesis
Philosophical Transactions of the Royal Society B: Biological Sciences 372:20150520.

https://doi.org/10.1098/rstb.2015.0520
- Google Scholar
1. Amonlirdviman K
2. Khare NA
3. Tree DR
4. Chen WS
5. Axelrod JD
6. Tomlin CJ
(2005) Mathematical modeling of planar cell polarity to understand domineering nonautonomy
Science 307:423–426.

https://doi.org/10.1126/science.1105471
- PubMed
- Google Scholar
1. Andrew DJ
2. Ewald AJ
(2010) Morphogenesis of epithelial tubes: Insights into tube formation, elongation, and elaboration
Developmental Biology 341:34–55.

https://doi.org/10.1016/j.ydbio.2009.09.024
- PubMed
- Google Scholar
(2009) Cellular and molecular mechanisms underlying the formation of biological tubes
Current Topics in Developmental Biology 89:137–162.

https://doi.org/10.1016/S0070-2153(09)89006-6
- PubMed
- Google Scholar
1. Beati H
2. Peek I
3. Hordowska P
4. Honemann-Capito M
5. Glashauser J
6. Renschler FA
7. Kakanj P
8. Ramrath A
9. Leptin M
10. Luschnig S
11. Wiesner S
12. Wodarz A
(2018) The adherens junction-associated LIM domain protein Smallish regulates epithelial morphogenesis
The Journal of Cell Biology 217:1079–1095.

https://doi.org/10.1083/jcb.201610098
- PubMed
- Google Scholar
(2016) Filopodial-Tension Model of Convergent-Extension of Tissues
PLOS Computational Biology 12:e1004952.

https://doi.org/10.1371/journal.pcbi.1004952
- PubMed
- Google Scholar
1. Burak Y
2. Shraiman BI
(2009) Order and stochastic dynamics in Drosophila planar cell polarity
PLoS Computational Biology 5:e1000628.

https://doi.org/10.1371/journal.pcbi.1000628
- PubMed
- Google Scholar
1. Buske P
2. Przybilla J
3. Loeffler M
4. Sachs N
5. Sato T
6. Clevers H
7. Galle J
(2012) On the biomechanics of stem cell niche formation in the gut--modelling growing organoids
FEBS Journal 279:3475–3487.

https://doi.org/10.1111/j.1742-4658.2012.08646.x
- PubMed
- Google Scholar
Book
1. Carroll TJ
2. Yu J
(2012)
Chapter Eight - The Kidney and Planar Cell Polarity.”

In: Yang Y. ingzi, editors. Current Topics in Developmental Biology. Academic Press. pp. 185–212.
- Google Scholar
1. Cerruti B
2. Puliafito A
3. Shewan AM
4. Yu W
5. Combes AN
6. Little MH
7. Chianale F
8. Primo L
9. Serini G
10. Mostov KE
11. Celani A
12. Gamba A
(2013) Polarity, cell division, and out-of-equilibrium dynamics control the growth of epithelial structures
The Journal of Cell Biology 203:359–372.

https://doi.org/10.1083/jcb.201305044
- PubMed
- Google Scholar
1. Cherry LM
2. Case SM
3. Kunkel JG
4. Wilson AC
(1979) Comparisons of frogs, humans, and chimpanzees
Science 204:435.

https://doi.org/10.1126/science.204.4391.435
- PubMed
- Google Scholar
1. Choi SC
2. Sokol SY
(2009) The involvement of lethal giant larvae and Wnt signaling in bottle cell formation in Xenopus embryos
Developmental Biology 336:68–75.

https://doi.org/10.1016/j.ydbio.2009.09.033
- PubMed
- Google Scholar
1. Chu CW
2. Sokol SY
(2016) Wnt proteins can direct planar cell polarity in vertebrate ectoderm
eLife 5:e16463.

https://doi.org/10.7554/eLife.16463
- PubMed
- Google Scholar
1. Chung S
2. Kim S
3. Andrew DJ
(2017) Uncoupling apical constriction from tissue invagination
eLife 6:e22235.

https://doi.org/10.7554/eLife.22235
- PubMed
- Google Scholar
1. Collinet C
2. Rauzi M
3. Lenne PF
4. Lecuit T
(2015) Local and tissue-scale forces drive oriented junction growth during tissue extension
Nature Cell Biology 17:1247–1258.

https://doi.org/10.1038/ncb3226
- PubMed
- Google Scholar
1. Croce J
2. Duloquin L
3. Lhomond G
4. McClay DR
5. Gache C
(2006) Frizzled5/8 is required in secondary mesenchyme cells to initiate archenteron invagination during sea urchin development
Development 133:547–557.

https://doi.org/10.1242/dev.02218
- PubMed
- Google Scholar
1. Dallon JC
2. Othmer HG
(2004) How cellular movement determines the collective force generated by the Dictyostelium discoideum slug
Journal of Theoretical Biology 231:203–222.

https://doi.org/10.1016/j.jtbi.2004.06.015
- PubMed
- Google Scholar
1. Dollar GL
2. Weber U
3. Mlodzik M
4. Sokol SY
(2005) Regulation of Lethal giant larvae by Dishevelled
Nature 437:1376–1380.

https://doi.org/10.1038/nature04116
- PubMed
- Google Scholar
1. Eguchi G
2. Eguchi Y
3. Nakamura K
4. Yadav MC
5. Millán JL
6. Tsonis PA
(2011) Regenerative capacity in newts is not altered by repeated regeneration and ageing
Nature Communications 2:.

https://doi.org/10.1038/ncomms1389
- Google Scholar
Conference
(2008)
In Silico Simulation of Epithelial Cell Tubulogenesis

30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. pp. 1036–1039.
- Google Scholar
1. Etournay R
2. Popović M
3. Merkel M
4. Nandi A
5. Blasse C
6. Aigouy B
7. Brandl H
8. Myers G
9. Salbreux G
10. Jülicher F
11. Eaton S
(2015) Interplay of cell dynamics and epithelial tension during morphogenesis of the Drosophila pupal wing
eLife 4:e07090.

https://doi.org/10.7554/eLife.07090
- PubMed
- Google Scholar
(2016) Designer matrices for intestinal stem cell and organoid culture
Nature 539:560–564.

https://doi.org/10.1038/nature20168
- PubMed
- Google Scholar
1. Gong Y
2. Mo C
3. Fraser SE
(2004) Planar cell polarity signalling controls cell division orientation during zebrafish gastrulation
Nature 430:689–693.

https://doi.org/10.1038/nature02796
- PubMed
- Google Scholar
1. Gould SJ
(1970) Evolutionary paleontology and the science of form
Earth-Science Reviews 6:77–119.

https://doi.org/10.1016/0012-8252(70)90027-9
- Google Scholar
1. Gould SJ
2. Eldredge N
(1993) Punctuated equilibrium comes of age
Nature 366:223–227.

https://doi.org/10.1038/366223a0
- PubMed
- Google Scholar
(2013) Artificial three-dimensional niches deconstruct pancreas development in vitro
Development 140:4452–4462.

https://doi.org/10.1242/dev.096628
- PubMed
- Google Scholar
1. Habib SJ
2. Chen BC
3. Tsai FC
4. Anastassiadis K
5. Meyer T
6. Betzig E
7. Nusse R
(2013) A localized Wnt signal orients asymmetric stem cell division in vitro
Science 339:1445–1448.

https://doi.org/10.1126/science.1231077
- PubMed
- Google Scholar
(2014) Theory of epithelial sheet morphology in three dimensions
PNAS 111:27–32.

https://doi.org/10.1073/pnas.1312076111
- PubMed
- Google Scholar
1. Hannezo E
2. Scheele C
3. Moad M
4. Drogo N
5. Heer R
6. Sampogna RV
7. van Rheenen J
8. Simons BD
(2017) A Unifying Theory of Branching Morphogenesis
Cell 171:242–255.

https://doi.org/10.1016/j.cell.2017.08.026
- PubMed
- Google Scholar
(2012) A model of epithelial invagination driven by collective mechanics of identical cells
Biophysical Journal 103:1069–1077.

https://doi.org/10.1016/j.bpj.2012.07.018
- PubMed
- Google Scholar
(2007) On the mechanism of wing size determination in fly development
PNAS 104:3835–3840.

https://doi.org/10.1073/pnas.0607134104
- PubMed
- Google Scholar
1. Humphries AC
2. Mlodzik M
(2018) From instruction to output: Wnt/PCP signaling in development and cancer
Current Opinion in Cell Biology 51:110–116.

https://doi.org/10.1016/j.ceb.2017.12.005
- PubMed
- Google Scholar
1. Kaplan NA
2. Tolwinski NS
(2010) Spatially defined Dsh-Lgl interaction contributes to directional tissue morphogenesis
Journal of Cell Science 123:3157–3165.

https://doi.org/10.1242/jcs.069898
- PubMed
- Google Scholar
(2009) Wnt9b signaling regulates planar cell polarity and kidney tubule morphogenesis
Nature Genetics 41:793–799.

https://doi.org/10.1038/ng.400
- PubMed
- Google Scholar
1. Kimberly EL
2. Hardin J
(1998) Bottle cells are required for the initiation of primary invagination in the sea urchin embryo
Developmental Biology 204:235–250.

https://doi.org/10.1006/dbio.1998.9075
- PubMed
- Google Scholar
(2008) Apical accumulation of Rho in the neural plate is important for neural plate cell shape change and neural tube formation
Molecular Biology of the Cell 19:2289–2299.

https://doi.org/10.1091/mbc.e07-12-1286
- PubMed
- Google Scholar
(2011) Simulating the mammalian blastocyst--molecular and mechanical interactions pattern the embryo
PLoS Computational Biology 7:e1001128.

https://doi.org/10.1371/journal.pcbi.1001128
- PubMed
- Google Scholar
(2017) Disruption of Core Planar Cell Polarity Signaling Regulates Renal Tubule Morphogenesis but Is Not Cystogenic
Current Biology 27:3120–3131.

https://doi.org/10.1016/j.cub.2017.09.011
- PubMed
- Google Scholar
1. Lane MC
2. Koehl MA
3. Wilt F
4. Keller R
(1993)
A Role for Regulated Secretion of Apical Extracellular Matrix during Epithelial Invagination in the Sea Urchin

Development 117:1049–1060.
- PubMed
- Google Scholar
1. Laprise P
2. Paul SM
3. Boulanger J
4. Robbins RM
5. Beitel GJ
6. Tepass U
(2010) Epithelial polarity proteins regulate Drosophila tracheal tube size in parallel to the luminal matrix pathway
Current Biology 20:55–61.

https://doi.org/10.1016/j.cub.2009.11.017
- PubMed
- Google Scholar
(2006) Establishment and maintenance of planar epithelial cell polarity by asymmetric cadherin bridges: a computer model
Developmental Dynamics 235:235–246.

https://doi.org/10.1002/dvdy.20617
- PubMed
- Google Scholar
(2006) Wnt/Frizzled signaling controls C. elegans gastrulation by activating actomyosin contractility
Current Biology 16:1986–1997.

https://doi.org/10.1016/j.cub.2006.08.090
- PubMed
- Google Scholar
1. Levayer R
2. Lecuit T
(2013) Oscillation and polarity of E-cadherin asymmetries control actomyosin flow patterns during morphogenesis
Developmental Cell 26:162–175.

https://doi.org/10.1016/j.devcel.2013.06.020
- PubMed
- Google Scholar
1. Li R
2. Bowerman B
(2010) Symmetry breaking in biology
Cold Spring Harbor Perspectives in Biology 2:a003475.

https://doi.org/10.1101/cshperspect.a003475
- PubMed
- Google Scholar
1. Li Y
2. Muffat J
3. Omer A
4. Bosch I
5. Lancaster MA
6. Sur M
7. Gehrke L
8. Knoblich JA
9. Jaenisch R
(2017) Induction of Expansion and Folding in Human Cerebral Organoids
Cell Stem Cell 20:385–396.

https://doi.org/10.1016/j.stem.2016.11.017
- PubMed
- Google Scholar
1. Libby E
2. Rainey PB
(2011) Exclusion rules, bottlenecks and the evolution of stochastic phenotype switching
Proceedings of the Royal Society B: Biological Sciences 278:3574–3583.

https://doi.org/10.1098/rspb.2011.0146
- Google Scholar
1. Little MH
(2017) Organoids: a Special Issue
Development 144:935–937.

https://doi.org/10.1242/dev.150292
- PubMed
- Google Scholar
(2016) Generating Cellular Diversity and Spatial Form: Wnt Signaling and the Evolution of Multicellular Animals
Developmental Cell 38:643–655.

https://doi.org/10.1016/j.devcel.2016.08.011
- PubMed
- Google Scholar
1. Long JT
2. Irwin L
3. Enomoto AC
4. Grow Z
5. Ranck J
6. Peeler MT
(2015) Jun N-terminal kinase activity is required for invagination but not differentiation of the sea urchin archenteron
Genesis 53:762–769.

https://doi.org/10.1002/dvg.22898
- PubMed
- Google Scholar
Book
1. Lye CM
2. Sanson B
(2011)
"Chapter Five - Tension and Epithelial Morphogenesis in Drosophila Early Embryos"

In: Labouesse Michel, editors. Current Topics in Developmental Biology, 95. Academic Press. pp. 145–187.
- Google Scholar
(2012) Morphogenesis in sea urchin embryos: linking cellular events to gene regulatory network states
Wiley Interdisciplinary Reviews: Developmental Biology 1:231–252.

https://doi.org/10.1002/wdev.18
- PubMed
- Google Scholar
1. Martik ML
2. McClay DR
(2017) New insights from a high-resolution look at gastrulation in the sea urchin, Lytechinus variegatus
Mechanisms of Development 148:3–10.

https://doi.org/10.1016/j.mod.2017.06.005
- PubMed
- Google Scholar
1. Martin-Belmonte F
2. Perez-Moreno M
(2012) Epithelial cell polarity, stem cells and cancer
Nature Reviews Cancer 12:23–38.

https://doi.org/10.1038/nrc3169
- Google Scholar
1. Minegishi K
2. Hashimoto M
3. Ajima R
4. Takaoka K
5. Shinohara K
6. Ikawa Y
7. Nishimura H
8. McMahon AP
9. Willert K
10. Okada Y
11. Sasaki H
12. Shi D
13. Fujimori T
14. Ohtsuka T
15. Igarashi Y
16. Yamaguchi TP
17. Shimono A
18. Shiratori H
19. Hamada H
(2017) A Wnt5 Activity Asymmetry and Intercellular Signaling via PCP Proteins Polarize Node Cells for Left-Right Symmetry Breaking
Developmental Cell 40:439–452.

https://doi.org/10.1016/j.devcel.2017.02.010
- PubMed
- Google Scholar
(2016) Shape Transformations of Epithelial Shells
Biophysical Journal 110:1670–1678.

https://doi.org/10.1016/j.bpj.2016.03.009
- PubMed
- Google Scholar
1. Monier B
2. Gettings M
3. Gay G
4. Mangeat T
5. Schott S
6. Guarner A
7. Suzanne M
(2015) Apico-basal forces exerted by apoptotic cells drive epithelium folding
Nature 518:245–248.

https://doi.org/10.1038/nature14152
- PubMed
- Google Scholar
(2015) Self-organization of polarized cerebellar tissue in 3D culture of human pluripotent stem cells
Cell Reports 10:537–550.

https://doi.org/10.1016/j.celrep.2014.12.051
- PubMed
- Google Scholar
1. Müller HA
2. Bossinger O
(2003) Molecular networks controlling epithelial cell polarity in development
Mechanisms of Development 120:1231–1256.

https://doi.org/10.1016/j.mod.2003.06.001
- PubMed
- Google Scholar
(2015) From discrete to continuum models of three-dimensional deformations in epithelial sheets
Biophysical Journal 109:154–163.

https://doi.org/10.1016/j.bpj.2015.05.019
- PubMed
- Google Scholar
1. Nagai T
2. Honda H
(2001) A dynamic cell model for the formation of epithelial tissues
Philosophical Magazine B 81:699–719.

https://doi.org/10.1080/13642810108205772
- Google Scholar
1. Nagai T
2. Honda H
(2009) Computer simulation of wound closure in epithelial tissues: Cell–basal-lamina adhesion
Physical Review E 80:.

https://doi.org/10.1103/PhysRevE.80.061903
- Google Scholar
1. Newman TJ
(2008) Grid-free models of multicellular systems, with an application to large-scale vortices accompanying primitive streak formation
Current Topics in Developmental Biology 81:157–182.

https://doi.org/10.1016/S0070-2153(07)81005-2
- PubMed
- Google Scholar
1. Nikolopoulou E
2. Galea GL
3. Rolo A
4. Greene ND
5. Copp AJ
(2017) Neural tube closure: cellular, molecular and biomechanical mechanisms
Development 144:552–566.

https://doi.org/10.1242/dev.145904
- PubMed
- Google Scholar
(2012) Planar cell polarity links axes of spatial dynamics in neural-tube closure
Cell 149:1084–1097.

https://doi.org/10.1016/j.cell.2012.04.021
- PubMed
- Google Scholar
1. Nissen SB
2. Perera M
3. Gonzalez JM
4. Morgani SM
5. Jensen MH
6. Sneppen K
7. Brickman JM
8. Trusina A
(2017) Four simple rules that are sufficient to generate the mammalian blastocyst
PLOS Biology 15:e2000737.

https://doi.org/10.1371/journal.pbio.2000737
- PubMed
- Google Scholar
(2012) Tubulogenesis: Src42A goes to great lengths in tube elongation
Current Biology 22:R446–R449.

https://doi.org/10.1016/j.cub.2012.04.033
- PubMed
- Google Scholar
1. Odell GM
2. Oster G
3. Alberch P
4. Burnside B
(1981)
The mechanical basis of morphogenesis. I. Epithelial folding and invagination

Developmental Biology 85:446–462.
- PubMed
- Google Scholar
1. Ossipova O
2. Kim K
3. Lake BB
4. Itoh K
5. Ioannou A
6. Sokol SY
(2014) Role of Rab11 in planar cell polarity and apical constriction during vertebrate neural tube closure
Nature Communications 5:3734.

https://doi.org/10.1038/ncomms4734
- PubMed
- Google Scholar
(2013) Three-dimensional epithelial morphogenesis in the developing Drosophila egg
Developmental Cell 24:400–410.

https://doi.org/10.1016/j.devcel.2013.01.017
- PubMed
- Google Scholar
(2015) Mechanisms of apical-basal axis orientation and epithelial lumen positioning
Trends in Cell Biology 25:476–485.

https://doi.org/10.1016/j.tcb.2015.04.002
- PubMed
- Google Scholar
1. Panousopoulou E
2. Green JB
(2016) Invagination of ectodermal placodes is driven by cell intercalation-mediated contraction of the suprabasal tissue canopy
PLOS Biology 14:e1002405.

https://doi.org/10.1371/journal.pbio.1002405
- PubMed
- Google Scholar
1. Polyakov O
2. He B
3. Swan M
4. Shaevitz JW
5. Kaschube M
6. Wieschaus E
(2014) Passive mechanical forces control cell-shape change during Drosophila ventral furrow formation
Biophysical Journal 107:998–1010.

https://doi.org/10.1016/j.bpj.2014.07.013
- PubMed
- Google Scholar
(2013) Physical models of mesoderm invagination in Drosophila embryo
Biophysical Journal 105:3–10.

https://doi.org/10.1016/j.bpj.2013.05.039
- PubMed
- Google Scholar
1. Rauzi M
2. Krzic U
3. Saunders TE
4. Krajnc M
5. Ziherl P
6. Hufnagel L
7. Leptin M
(2015) Embryo-scale tissue mechanics during Drosophila gastrulation movements
Nature Communications 6:.

https://doi.org/10.1038/ncomms9677
- Google Scholar
(2006) Individual cell migration serves as the driving force for optic vesicle evagination
Science 313:1130–1134.

https://doi.org/10.1126/science.1127144
- PubMed
- Google Scholar
(2013) Polarity in mammalian epithelial morphogenesis
Cold Spring Harbor Perspectives in Biology 5:a013789.

https://doi.org/10.1101/cshperspect.a013789
- PubMed
- Google Scholar
1. Saburi S
2. Hester I
3. Fischer E
4. Pontoglio M
5. Eremina V
6. Gessler M
7. Quaggin SE
8. Harrison R
9. Mount R
10. McNeill H
(2008) Loss of Fat4 disrupts PCP signaling and oriented cell division and leads to cystic kidney disease
Nature Genetics 40:1010–1015.

https://doi.org/10.1038/ng.179
- PubMed
- Google Scholar
(2018) Radially patterned cell behaviours during tube budding from an epithelium
eLife 7:e35717.

https://doi.org/10.7554/eLife.35717
- PubMed
- Google Scholar
1. Schierenberg E
2. Junkersdorf B
(1992) The role of eggshell and underlying vitelline membrane for normal pattern formation in the early C. elegans embryo
Roux's Archives of Developmental Biology 202:10–16.

https://doi.org/10.1007/BF00364592
- PubMed
- Google Scholar
(2016) Sequential development of apical-basal and planar polarities in aggregating epitheliomuscular cells of Hydra
Developmental Biology 412:148–159.

https://doi.org/10.1016/j.ydbio.2016.02.022
- PubMed
- Google Scholar
(2010) Rho-Kinase Directs Bazooka/Par-3 Planar Polarity during Drosophila Axis Elongation
Developmental Cell 19:377–388.

https://doi.org/10.1016/j.devcel.2010.08.011
- PubMed
- Google Scholar
1. Song H
2. Hu J
3. Chen W
4. Elliott G
5. Andre P
6. Gao B
7. Yang Y
(2010) Planar cell polarity breaks bilateral symmetry by controlling ciliary positioning
Nature 466:378–382.

https://doi.org/10.1038/nature09129
- PubMed
- Google Scholar
(2016)
On the growth and form of cortical convolutions

Nature Physics 12588:588–593.
- Google Scholar
(2017) Abl Regulates Planar Polarized Junctional Dynamics through β-Catenin Tyrosine Phosphorylation
Developmental Cell 41:330.

https://doi.org/10.1016/j.devcel.2017.04.018
- PubMed
- Google Scholar
1. Tamulonis C
2. Postma M
3. Marlow HQ
4. Magie CR
5. de Jong J
6. Kaandorp J
(2011) A cell-based model of Nematostella vectensis gastrulation including bottle cell formation, invagination and zippering
Developmental Biology 351:217–228.

https://doi.org/10.1016/j.ydbio.2010.10.017
- PubMed
- Google Scholar
(2009) Liver progenitor cells fold up a cell monolayer into a double-layered structure during tubular morphogenesis
Molecular Biology of the Cell 20:2486–2494.

https://doi.org/10.1091/mbc.e08-02-0177
- PubMed
- Google Scholar
Book
1. Thompson DW
(1942)
On Growth and Form

Cambridge: Cambridge University Press.
- Google Scholar
1. Wang Y
2. Badea T
3. Nathans J
(2006) Order from disorder: Self-organization in mammalian hair patterning
PNAS 103:19800–19805.

https://doi.org/10.1073/pnas.0609712104
- PubMed
- Google Scholar
(2013) The Frizzled-dependent planar polarity pathway locally promotes E-cadherin turnover via recruitment of RhoGEF2
Development 140:1045–1054.

https://doi.org/10.1242/dev.088724
- PubMed
- Google Scholar
(2010) Drosophila vitelline membrane assembly: a critical role for an evolutionarily conserved cysteine in the "VM domain" of sV23
Developmental Biology 347:360–368.

https://doi.org/10.1016/j.ydbio.2010.08.037
- PubMed
- Google Scholar
(2013) Wg and Wnt4 provide long-range directional input to planar cell polarity orientation in Drosophila
Nature Cell Biology 15:1045–1055.

https://doi.org/10.1038/ncb2806
- PubMed
- Google Scholar
1. Zegers MM
(2014) 3D in vitro cell culture models of tube formation
Seminars in Cell & Developmental Biology 31:132–140.

https://doi.org/10.1016/j.semcdb.2014.02.016
- PubMed
- Google Scholar
1. Zorn AM
2. Wells JM
(2009) Vertebrate endoderm development and organ formation
Annual Review of Cell and Developmental Biology 25:221–251.

https://doi.org/10.1146/annurev.cellbio.042308.113344
- PubMed
- Google Scholar

Article and author information

Author details

Silas Boye Nissen
1. Center for Models of Life, Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark
2. StemPhys, Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark
Contribution
Conceptualization, Resources, Data curation, Software, Formal analysis, Supervision, 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-0002-9473-4755
Steven Rønhild

Center for Models of Life, Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark

Contribution
Software, Writing—review and editing

Competing interests
No competing interests declared
Ala Trusina
1. Center for Models of Life, Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark
2. StemPhys, Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark
Contribution
Conceptualization, Resources, Data curation, Software, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing—original draft, Project administration, Writing—review and editing

For correspondence
trusina@nbi.ku.dk

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-1945-454X
Kim Sneppen

Center for Models of Life, Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark

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

For correspondence
sneppen@nbi.dk

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-9820-3567

Funding

Danmarks Grundforskningsfond (DNRF116)

Silas Boye Nissen
Ala Trusina

Seventh Framework Programme (FP/2007/2013/ERC no. 740704)

Kim Sneppen

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

Acknowledgements

We thank Anne Grapin-Botton for insightful discussions on organoids and Sigurd Carlsen for discussion on tube formation. This research has received funding from the Danish National Research Foundation (grant number: DNRF116) and the European Research Council under the European Union's Seventh Framework Programme (FP/2007 2013)/ERC Grant Agreement number 740704.

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

6,963

views
707

downloads
25

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)

Silas Boye Nissen
Steven Rønhild
Ala Trusina
Kim Sneppen

(2018)

Theoretical tool bridging cell polarities with development of robust morphologies

eLife 7:e38407.

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

Share this article

Cite this article

Two symmetry-breaking events, gain of apical-basal (AB) polarity and planar cell polarity (PCP), on cellular level coincide with the appearance of a rich set of morphologies.

Cells are modeled as interacting particles with a polarity-dependent potential.

Development of 8000 cells from a compact aggregate starting at time 0.

Different morphologies can be obtained by varying boundary conditions (Figure 4—video 1).

The number of complex folds in a growing organoid depends on the generation time and the pressure from the surrounding medium (Figure 5—video 1).

The length and width of tubes are set by the strength of planar cell polarity (PCP, λ3).

External constraints on apical-basal (AB) polarity and planar cell polarity (PCP) can initiate invagination and drive gastrulation in sea urchin.

Author details

Silas Boye Nissen

Contribution

Competing interests

Steven Rønhild

Contribution

Competing interests

Ala Trusina

Contribution

For correspondence

Competing interests

Kim Sneppen

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

Further reading

The length and width of tubes are set by the strength of planar cell polarity (PCP, λ₃).