Shape Transformations of Epithelial Shells

Local Mechanical Modulation-Driven Evagination in Invaginated Epithelia

Local cells can actively create reverse bending (evagination) in invaginated epithelia, which plays a crucial role in the formation of elaborate organisms. However, the precise physical mechanism driving the evagination remains elusive. Here, we present a three-dimensional vertex model, incorporating the intrinsic cell polarity, to explore the complex morphogenesis induced by local mechanical modulations. We find that invaginated tissues can spontaneously generate local reverse bending due to the shift of the apicobasal polarity. Their exact shapes can be analytically determined by the local apicobasal differential tension and the internal stress. Our continuum theory exhibits three regions in a phase diagram controlled by these two parameters, showing curvature transitions from ordered to disordered states. Additionally, we delve into epithelial curvature transition induced by the nucleus repositioning, revealing its active contribution to the apicobasal force generation. The uncovered mechanical principles could potentially guide more studies on epithelial folding in diverse systems.

Graph topological transformations in space-filling cell aggregates

Cell rearrangements are fundamental mechanisms driving large-scale deformations of living tissues. In three-dimensional (3D) space-filling cell aggregates, cells rearrange through local topological transitions of the network of cell-cell interfaces, which is most conveniently described by the vertex model. Since these transitions are not yet mathematically properly formulated, the 3D vertex model is generally difficult to implement. The few existing implementations rely on highly customized and complex software-engineering solutions, which cannot be transparently delineated and are thus mostly non-reproducible. To solve this outstanding problem, we propose a reformulation of the vertex model. Our approach, called Graph Vertex Model (GVM), is based on storing the topology of the cell network into a knowledge graph with a particular data structure that allows performing cell-rearrangement events by simple graph transformations. Importantly, when these same transformations are applied to a two-dimensional (2D) polygonal cell aggregate, they reduce to a well-known T1 transition, thereby generalizing cell-rearrangements in 2D and 3D space-filling packings. This result suggests that the GVM's graph data structure may be the most natural representation of cell aggregates and tissues. We also develop a Python package that implements GVM, relying on a graph-database-management framework Neo4j. We use this package to characterize an order-disorder transition in 3D cell aggregates, driven by active noise and we find aggregates undergoing efficient ordering close to the transition point. In all, our work showcases knowledge graphs as particularly suitable data models for structured storage, analysis, and manipulation of tissue data.

Mechanical control of neural plate folding by apical domain alteration

Vertebrate neural tube closure is associated with complex changes in cell shape and behavior, however, the relative contribution of these processes to tissue folding is not well understood. At the onset of Xenopus neural tube folding, we observed alternation of apically constricted and apically expanded cells. This apical domain heterogeneity was accompanied by biased cell orientation along the anteroposterior axis, especially at neural plate hinges, and required planar cell polarity signaling. Vertex models suggested that dispersed isotropically constricting cells can cause the elongation of adjacent cells. Consistently, in ectoderm, cell-autonomous apical constriction was accompanied by neighbor expansion. Thus, a subset of isotropically constricting cells may initiate neural plate bending, whereas a 'tug-of-war' contest between the force-generating and responding cells reduces its shrinking along the body axis. This mechanism is an alternative to anisotropic shrinking of cell junctions that are perpendicular to the body axis. We propose that apical domain changes reflect planar polarity-dependent mechanical forces operating during neural folding.

Single and two-cells shape analysis from energy functionals for three-dimensional vertex models

Vertex models have been extensively used for simulating the evolution of multicellular systems, and have given rise to important global properties concerning their macroscopic rheology or jamming transitions. These models are based on the definition of an energy functional, which fully determines the cellular response and conclusions. While two-dimensional vertex models have been widely employed, three-dimensional models are far more scarce, mainly due to the large amount of configurations that they may adopt and the complex geometrical transitions they undergo. We here investigate the shape of single and two-cells configurations as a function of the energy terms, and we study the dependence of the final shape on the model parameters: namely the exponent of the term penalising cell-cell adhesion and surface contractility. In single cell analysis, we deduce analytically the radius and limit values of the contractility for linear and quadratic surface energy terms, in 2D and 3D. In two-cells systems, symmetrical and asymmetrical, we deduce the evolution of the aspect ratio and the relative radius. While in functionals with linear surface terms yield the same aspect ratio in 2D and 3D, the configurations when using quadratic surface terms are distinct. We relate our results with well-known solutions from capillarity theory, and verify our analytical findings with a three-dimensional vertex model.

Cell-Based Modeling of Tissue Developing in the Scaffold Pores of Varying Cross-Sections

In this work, we present a mathematical model of cell growth in the pores of a perfusion bioreactor through which a nutrient solution is pumped. We have developed a 2-D vertex model that allows us to reproduce the microscopic dynamics of the microenvironment of cells and describe the occupation of the pore space with cells. In this model, each cell is represented by a polygon; the number of vertices and shapes may change over time. The model includes mitotic cell division and intercalation. We study the impact of two factors on cell growth. On the one hand, we consider a channel of variable cross-section, which models a scaffold with a porosity gradient. On the other hand, a cluster of cells grows under the influence of a nutrient solution flow, which establishes a non-uniform distribution of shear stresses in the pore space. We present the results of numerical simulation of the tissue growth in a wavy channel. The model allows us to obtain complete microscopic information that includes the dynamics of intracellular pressure, the local elastic energy, and the characteristics of cell populations. As we showed, in a functional-graded scaffold, the distribution of the shear stresses in the pore space has a complicated structure, which implies the possibility of controlling the growth zones by varying the pore geometry.

Rho activation drives luminal collapse and eversion in epithelial acini

Epithelial cells lining a gland and cells grown in a soft extracellular matrix polarize with apical proteins exposed to the lumen and basal proteins in contact with the extracellular matrix. Alterations to polarity, including an apical-out polarity, occur in human cancers. Although some aberrant polarity states may result from altered protein trafficking, recent observations of an extraordinary tissue-level inside-out unfolding suggest an alternative pathway for altered polarity. Because mechanical alterations are common in human cancer, including an upregulation of RhoA-mediated actomyosin tension in acinar epithelia, we explored whether perturbing mechanical homeostasis could cause apical-out eversion. Acinar eversion was robustly induced by direct activation of RhoA in normal and tumor epithelial acini, or indirect activation of RhoA through blockage of β1-integrins, disruption of the LINC complex, oncogenic Ras activation, or Rac1 inhibition. Furthermore, laser ablation of a portion of the untreated acinus was sufficient to induce eversion. Analyses of acini revealed high curvature and low phosphorylated myosin in the apical cell surfaces relative to the basal surfaces. A vertex-based mathematical model that balances tension at cell-cell interfaces revealed a fivefold greater basal cell surface tension relative to the apical cell surface tension. The model suggests that the difference in surface energy between the apical and basal surfaces is the driving force for acinar eversion. Our findings raise the possibility that a loss of mechanical homeostasis may cause apical-out polarity states in human cancers.

Stability of asymmetric cell division: A deformable cell model of cytokinesis applied to C. elegans

Cell division during early embryogenesis is linked to key morphogenic events such as embryo symmetry breaking and tissue patterning. It is thought that the physical surrounding of cells together with cell intrinsic cues act as a mechanical "mold," guiding cell division to ensure these events are robust. To quantify how cell division is affected by the mechanical and geometrical environment, we present a novel computational mechanical model of cytokinesis, the final phase of cell division. Simulations with the model reproduced experimentally observed furrow dynamics and describe the volume ratio of daughter cells in asymmetric cell divisions, based on the position and orientation of the mitotic spindle. For dividing cells in geometrically confined environments, we show how the orientation of confinement relative to the division axis modulates the volume ratio in asymmetric cell division. Further, we quantified how cortex viscosity and surface tension determine the shape of a dividing cell and govern bubble-instabilities in asymmetric cell division. Finally, we simulated the formation of the three body axes via sequential (a)symmetric divisions up until the six-cell stage of early C. elegans development, which proceeds within the confines of an eggshell. We demonstrate how model input parameters spindle position and orientation provide sufficient information to reliably predict the volume ratio of daughter cells during the cleavage phase of development. However, for egg geometries perturbed by compression, the model predicts that a change in confinement alone is insufficient to explain experimentally observed differences in cell volume. This points to an effect of the compression on the spindle positioning mechanism. Additionally, the model predicts that confinement stabilizes asymmetric cell divisions against bubble-instabilities.

Mechanical control of neural plate folding by apical domain alteration

Vertebrate neural tube closure is associated with complex changes in cell shape and behavior, however, the relative contribution of these processes to tissue folding is not well understood. In this study, we evaluated morphology of the superficial cell layer in the Xenopus neural plate. At the stages corresponding to the onset of tissue folding, we observed the alternation of cells with apically constricting and apically expanding apical domains. The cells had a biased orientation along the anteroposterior (AP) axis. This apical domain heterogeneity required planar cell polarity (PCP) signaling and was especially pronounced at neural plate hinges. Vertex model simulations suggested that spatially dispersed isotropically constricting cells cause the elongation of their non-constricting counterparts along the AP axis. Consistent with this hypothesis, cell-autonomous induction of apical constriction in Xenopus ectoderm cells was accompanied by the expansion of adjacent non-constricting cells. Our observations indicate that a subset of isotropically constricting cells can initiate neural plate bending, whereas a 'tug-of-war' contest between the force-generating and responding cells reduces its shrinking along the AP axis. This mechanism is an alternative to anisotropic shrinking of cell junctions that are perpendicular to the body axis. We propose that neural folding relies on PCP-dependent transduction of mechanical signals between neuroepithelial cells.

Impact of environmental asymmetry on epithelial morphogenesis

Epithelial folding is a universal biological phenomenon in morphogenesis, typical examples being brain gyri, villi of the intestinal tract, and imaginal discs in invertebrates. During epithelial morphogenesis, the physical constraints imposed by the surrounding microenvironment on epithelial tissue play critical roles in folding morphology. In this study, we focused on the asymmetry of the environmental constraints sandwiching the epithelial sheet and introduced the degree of asymmetry, which indicates whether the basal or apical side of the epithelium is closer to the constraint wall. Then, we investigated the relationship between the degree of asymmetry and epithelial folding morphology using three-dimensional vertex simulations. The results show that the folding patterns of the epithelial sheets change from spot patterns to labyrinth patterns and then to hole patterns as the degree of asymmetry changes. Furthermore, we examined the pattern formation in terms of the equation of out-of-plane displacement of the sheet derived from the mechanical energy functional.

Thickness of epithelia on wavy substrates: measurements and continuous models

We measured the thickness of MDCK epithelia grown on substrates with a sinusoidal profile. We show that while at long wavelength the profile of the epithelium follows that of the substrate, at short wavelengths cells are thicker in valleys than on ridges. This is reminiscent of the so-called «healing length in the case of a thin liquid film wetting a rough solid substrate. We explore the ability of continuum mechanics models to account for these observations. Modeling the epithelium as a thin liquid film, with surface tension, does not fully account for the measurements. Neither does modeling the epithelium as a thin incompressible elastic film. On the contrary, the addition of an apical active stress gives satisfactory agreement with measurements, with one fitting parameter, the ratio between the active stress and the elastic modulus.

Transformation of Amorphous Nanobowls to Crystalline Ellipsoids Induced by Trans-Cis Isomerization of Azobenzene

The stimuli-responsive transition of nanostructures from amorphous to crystalline state is of high interest in polymer science, but is still challenging. Herein, we demonstrate the transformation of amorphous nanobowls to crystalline ellipsoids triggered by UV induced trans-cis isomerization, using an azobenzene-containing amphiphilic homopolymer (PAzoAA) as building block. The amide bond and azobenzene pendants are introduced to the side chain of PAzoAA to afford hydrogen bonding and π-π interaction, which promotes the formation of nanobowls rather than spherical nanostructures. Upon exposed to UV irradiation, trans-cis isomerization of azobenzene pendants occurs, leading to the increase of hydrophilicity and destruction of π-π interaction, further resulting in the disassembly of the nanobowls. Then the PAzoAA re-assembles to form crystalline ellipsoids instead of amorphous nanostructures when recovered at 70°C without UV light. We further confirm that the high incubation temperature after UV irradiation is critical for the cis-trans transformation and the high mobility of the polymer chains to facilitate the regular rearrangement of azobenzene pendants. Overall, we propose a facile method to achieve the transformation of amorphous nanobowls to crystalline ellipsoids, which may bring new insight into preparation of crystalline nanoparticles using amorphous precursors. This article is protected by copyright. All rights reserved.

Appreciating the role of cell shape changes in the mechanobiology of epithelial tissues

The wide range of epithelial cell shapes reveals the complexity and diversity of the intracellular mechanisms that serve to construct their morphology and regulate their functions. Using mechanosensitive steps, epithelial cells can sense a variety of different mechanochemical stimuli and adapt their behavior by reshaping their morphology. These changes of cell shape rely on a structural reorganization in space and time that generates modifications of the tensional state and activates biochemical cascades. Recent studies have started to unveil how the cell shape maintenance is involved in mechanical homeostatic tasks to sustain epithelial tissue folding, identity, and self-renewal. Here, we review relevant works that integrated mechanobiology to elucidate some of the core principles of how cell shape may be conveyed into spatial information to guide collective processes such as epithelial morphogenesis. Among many other parameters, we show that the regulation of the cell shape can be understood as the result of the interplay between two counteracting mechanisms: actomyosin contractility and intercellular adhesions, and that both do not act independently but are functionally integrated to operate on molecular, cellular, and tissue scales. We highlight the role of cadherin-based adhesions in force-sensing and mechanotransduction, and we report recent developments that exploit physics of liquid crystals to connect cell shape changes to orientational order in cell aggregates. Finally, we emphasize that the further intermingling of different disciplines to develop new mechanobiology assays will lead the way toward a unified picture of the contribution of cell shape to the pathophysiological behavior of epithelial tissues.

Morphologies of compressed active epithelial monolayers

Using a three-dimensional active vertex model, we numerically study the shapes of strained unsupported epithelial monolayers subject to active junctional noise due to stochastic binding and unbinding of myosin. We find that while uniaxial, biaxial, and isotropic in-plane compressive strains do lead to the formation of longitudinal, herringbone pattern, and labyrinthine folds, respectively, the villus morphology characteristic of, e.g., the small intestine appears only if junctional tension fluctuations are strong enough to fluidize the tissue. Moreover, the fluidized epithelium features villi even in the absence of compressive strain provided that the apico-basal differential surface tension is large enough. We analyze several details of the different epithelial forms including the role of strain rate and the modulation of tissue thickness across folds. Our results show that even unsupported, non-patterned epithelia can form nontrivial morphologies.

A three-dimensional vertex model forDrosophilasalivary gland invagination

During epithelial morphogenesis, force generation at the cellular level not only causes cell deformation, but may also produce coordinated cell movement and rearrangement on the tissue level. In this paper, we use a novel three-dimensional vertex model to explore the roles of cellular forces during the formation of the salivary gland in theDrosophilaembryo. Representing the placode as an epithelial sheet of initially columnar cells, we focus on the spatial and temporal patterning of contractile forces due to three actomyosin pools: the apicomedial actomyosin in the pit of the placode, junctional actomyosin arcs outside the pit, and a supracellular actomyosin cable along the circumference of the placode. In anin silico'wild type' model, these pools are activated at different times according to experimental data. To identify the role of each myosin pool, we have also simulated variousin silico'mutants' in which only one or two of the myosin pools are activated. We find that the apicomedial myosin initiates a small dimple in the pit, but this is not essential for the overall invagination of the placode. The myosin arcs are the main driver of invagination and are responsible for the internalization of the apical surface. The circumferential actomyosin cable acts to constrict the opening of the developing tube, and is responsible for forming a properly shaped lumen. Cell intercalation tends to facilitate the invagination, but the geometric constraints of our model only allow a small number of intercalations, and their effect is minor. The placode invagination predicted by the model is in general agreement with experimental observations. It confirms some features of the current 'belt-and-braces' model for the process, and provides new insights on the separate roles of the various myosin pools and their spatio-temporal coordination.

The complex three-dimensional organization of epithelial tissues

Understanding the cellular organization of tissues is key to developmental biology. In order to deal with this complex problem, researchers have taken advantage of reductionist approaches to reveal fundamental morphogenetic mechanisms and quantitative laws. For epithelia, their two-dimensional representation as polygonal tessellations has proved successful for understanding tissue organization. Yet, epithelial tissues bend and fold to shape organs in three dimensions. In this context, epithelial cells are too often simplified as prismatic blocks with a limited plasticity. However, there is increasing evidence that a realistic approach, even from a reductionist perspective, must include apico-basal intercalations (i.e. scutoidal cell shapes) for explaining epithelial organization convincingly. Here, we present an historical perspective about the tissue organization problem. Specifically, we analyze past and recent breakthroughs, and discuss how and why simplified, but realistic, in silico models require scutoidal features to address key morphogenetic events.

Collective cell mechanics of epithelial shells with organoid-like morphologies

The study of organoids, artificially grown cell aggregates with the functionality and small-scale anatomy of real organs, is one of the most active areas of research in biology and biophysics, yet the basic physical origins of their different morphologies remain poorly understood. Here, we propose a mechanistic theory of epithelial shells which resemble small-organoid morphologies. Using a 3D surface tension-based vertex model, we reproduce the characteristic shapes from branched and budded to invaginated structures. We find that the formation of branched morphologies relies strongly on junctional activity, enabling temporary aggregations of topological defects in cell packing. To elucidate our numerical results, we develop an effective elasticity theory, which allows one to estimate the apico-basal polarity from the tissue-scale modulation of cell height. Our work provides a generic interpretation of the observed epithelial shell morphologies, highlighting the role of physical factors such as differential surface tension, cell rearrangements, and tissue growth.

The physical origins of man-made tissue morphologies with organ-like microscopic anatomy and functionality remain poorly understood. Here, authors propose a mechanistic theory of these structures, employing a 3D vertex model to reproduce the characteristic morphologies such as branched shapes.

Development of a New 3D Hybrid Model for Epithelia Morphogenesis

Many epithelial developmental processes like cell migration and spreading, cell sorting, or T1 transitions can be described as planar deformations. As such, they can be studied using two-dimensional tools and vertex models that can properly predict collective dynamics. However, many other epithelial shape changes are characterized by out-of-plane mechanics and three-dimensional effects, such as bending, cell extrusion, delamination, or invagination. Furthermore, during planar cell dynamics or tissue repair in monolayers, spatial intercalation between the apical and basal sides has even been detected. Motivated by this lack of symmetry with respect to the midsurface, we here present a 3D hybrid model that allows us to model differential contractility at the apical, basal or lateral sides. We use the model to study the effects on wound closure of solely apical or lateral contractile contributions and show that an apical purse-string can be sufficient for full closure when it is accompanied by volume preservation.

Curling of epithelial monolayers reveals coupling between active bending and tissue tension

Epithelial monolayers are two-dimensional cell sheets which compartmentalize the body and organs of multicellular organisms. Their morphogenesis during development or pathology results from patterned endogenous and exogenous forces and their interplay with tissue mechanical properties. In particular, bending of epithelia is thought to result from active torques generated by the polarization of myosin motors along their apicobasal axis. However, the contribution of these out-of-plane forces to morphogenesis remains challenging to evaluate because of the lack of direct mechanical measurement. Here we use epithelial curling to characterize the out-of-plane mechanics of epithelial monolayers. We find that curls of high curvature form spontaneously at the free edge of epithelial monolayers devoid of substrate in vivo and in vitro. Curling originates from an enrichment of myosin in the basal domain that generates an active spontaneous curvature. By measuring the force necessary to flatten curls, we can then estimate the active torques and the bending modulus of the tissue. Finally, we show that the extent of curling is controlled by the interplay between in-plane and out-of-plane stresses in the monolayer. Such mechanical coupling emphasizes a possible role for in-plane stresses in shaping epithelia during morphogenesis.

A Mechanical Instability in Planar Epithelial Monolayers Leads to Cell Extrusion

In cell extrusion, a cell embedded in an epithelial monolayer loses its apical or basal surface and is subsequently squeezed out of the monolayer by neighboring cells. Cell extrusions occur during apoptosis, epithelial-mesenchymal transition, or precancerous cell invasion. They play important roles in embryogenesis, homeostasis, carcinogenesis, and many other biological processes. Although many of the molecular factors involved in cell extrusion are known, little is known about the mechanical basis of cell extrusion. We used a three-dimensional (3D) vertex model to investigate the mechanical stability of cells arranged in a monolayer with 3D foam geometry. We found that when the cells composing the monolayer have homogeneous mechanical properties, cells are extruded from the monolayer when the symmetry of the 3D geometry is broken because of an increase in cell density or a decrease in the number of topological neighbors around single cells. Those results suggest that mechanical instability inherent in the 3D foam geometry of epithelial monolayers is sufficient to drive epithelial cell extrusion. In the situation in which cells in the monolayer actively generate contractile or adhesive forces under the control of intrinsic genetic programs, the forces act to break the symmetry of the monolayer, leading to cell extrusion that is directed to the apical or basal side of the monolayer by the balance of contractile and adhesive forces on the apical and basal sides. Although our analyses are based on a simple mechanical model, our results are in accordance with observations of epithelial monolayers in vivo and consistently explain cell extrusions under a wide range of physiological and pathophysiological conditions. Our results illustrate the importance of a mechanical understanding of cell extrusion and provide a basis by which to link molecular regulation to physical processes.

Viscoelastic voyages - Biophysical perspectives on cell intercalation during Drosophila gastrulation

Developmental processes are driven by a combination of cytoplasmic, cortical, and surface-associated forces. However, teasing apart the contributions of these forces and how a viscoelastic cell responds has long been a key question in developmental biology. Recent advances in applying biophysical approaches to these questions is leading to a fundamentally new understanding of morphogenesis. In this review, we discuss how computational analysis of experimental findings and in silico modeling of Drosophila gastrulation processes has led to a deeper comprehension of the physical principles at work in the early embryo. We also summarize many of the emerging methodologies that permit biophysical analysis as well as those that provide direct and indirect measurements of force directions and magnitudes. Finally, we examine the multiple frameworks that have been used to model tissue and cellular behaviors.

Feedback from Tissue Mechanics Self-Organizes Efficient Outgrowth of Plant Organ

Plant organ outgrowth superficially appears like the continuous mechanical deformation of a sheet of cells. Yet, how precisely cells as individual mechanical entities can act to morph a tissue reliably and efficiently into three dimensions during outgrowth is still puzzling, especially when cells are tightly connected as in plant tissue. In plants, the mechanics of cells within a tissue is particularly well-defined because individual cell growth is essentially the mechanical yielding of the cell wall in response to internal turgor pressure. Cell-wall stiffness is controlled by biological signaling, which is impacted by stresses, and hence, cell growth is observed to respond to mechanical stresses building up within a tissue. What is the role of the mechanical feedback during morphing of tissue in three dimensions? Here, we develop a three-dimensional vertex model to investigate tissue mechanics at the onset of organ outgrowth at the tip of a plant shoot. We find that organ height is primarily governed by the ratio of growth rates of faster-growing cells initiating the organ versus slower-growing cells surrounding them. Remarkably, the outgrowth rate is higher when cell growth responds to the tissue-wide mechanical stresses. Our quantitative analysis of simulation data shows that tissue mechanical feedback on cell growth can act via a twofold mechanism. First, the feedback guides patterns of cellular growth. Second, the feedback modifies the stress patterns on the cells, consequently amplifying and propagating growth anisotropies. This mechanism may allow plants to grow organs efficiently out of the meristem by reorganizing the cellular growth rather than inflating growth rates.

Mechanobiology of cells and cell systems, such as organoids

Organoids are in vitro 3D self-organizing tissues that mimic embryogenesis. Organoid research is advancing at a tremendous pace, since it offers great opportunities for disease modeling, drug development and screening, personalized medicine, as well as understanding organogenesis. Mechanobiology of organoids is an unexplored area, which can shed light to several unexplained aspects of self-organization behavior in organogenesis. It is becoming evident that collective cell behavior is distinctly different from individual cells' conduct against certain stimulants. Inherently consisting of higher number of degrees of freedom for cell motility and more complex cell-to-cell and cell-to-extracellular matrix behavior, understanding mechanotransduction in organoids is even more challenging compared with cell communities in 2D culture conditions. Yet, deciphering mechanobiology of organoids can help us understand effects of mechanical cues in health and disease, and translate findings of basic research toward clinical diagnosis and therapy.

An energy landscape approach to understanding variety and robustness in tissue morphogenesis

During morphogenesis in development, multicellular tissues deform by mechanical forces induced by spatiotemporally regulated cellular activities, such as cell proliferation and constriction. Various morphologies are formed because of various spatiotemporal combinations and sequences of multicellular activities. Despite its potential to variations, morphogenesis is a surprisingly robust process, in which qualitatively similar morphologies are reproducibly formed even under spatiotemporal fluctuation of multicellular activities. To understand these essential characteristics of tissue morphogenesis, which involves the coexistence of various morphologies and robustness of the morphogenetic process, in this study, we propose a novel approach to capture the overall view of morphogenesis from mechanical viewpoints. This approach will enable visualization of the energy landscape, which includes morphogenetic processes induced by admissible histories of cellular activities. This approach was applied to investigate the morphogenesis of a sheet-like tissue with curvature, where it deformed to a concave or convex morphology depending on the history of growth and constriction. Qualitatively different morphologies were produced by bifurcation of the valley in the energy landscape. The depth and steepness of the valley near the stable states represented the degree of robustness to fluctuations of multicellular activities. Furthermore, as a realistic example, we showed an application of this approach to luminal folding observed in the initial stage of intestinal villus formation. This approach will be helpful to understand the mechanism of how various morphologies are formed and how tissues reproducibly achieve specific morphologies.

A Preferred Curvature-Based Continuum Mechanics Framework for Modeling Embryogenesis

Mechanics plays a key role in the development of higher organisms. However, understanding this relationship is complicated by the difficulty of modeling the link between local forces generated at the subcellular level and deformations observed at the tissue and whole-embryo levels. Here we propose an approach first developed for lipid bilayers and cell membranes, in which force-generation by cytoskeletal elements enters a continuum mechanics formulation for the full system in the form of local changes in preferred curvature. This allows us to express and solve the system using only tissue strains. Locations of preferred curvature are simply related to products of gene expression. A solution, in that context, means relaxing the system's mechanical energy to yield global morphogenetic predictions that accommodate a tendency toward the local preferred curvature, without a need to explicitly model force-generation mechanisms at the molecular level. Our computational framework, which we call SPHARM-MECH, extends a 3D spherical harmonics parameterization known as SPHARM to combine this level of abstraction with a sparse shape representation. The integration of these two principles allows computer simulations to be performed in three dimensions on highly complex shapes, gene expression patterns, and mechanical constraints. We demonstrate our approach by modeling mesoderm invagination in the fruit-fly embryo, where local forces generated by the acto-myosin meshwork in the region of the future mesoderm lead to formation of a ventral tissue fold. The process is accompanied by substantial changes in cell shape and long-range cell movements. Applying SPHARM-MECH to whole-embryo live imaging data acquired with light-sheet microscopy reveals significant correlation between calculated and observed tissue movements. Our analysis predicts the observed cell shape anisotropy on the ventral side of the embryo and suggests an active mechanical role of mesoderm invagination in supporting the onset of germ-band extension.

Mathematical Modeling of Tissue Folding and Asymmetric Tissue Flow during Epithelial Morphogenesis

Recent studies have revealed that intrinsic, individual cell behavior can provide the driving force for deforming a two-dimensional cell sheet to a three-dimensional tissue without the need for external regulatory elements. However, whether intrinsic, individual cell behavior could actually generate the force to induce tissue deformation was unclear, because there was no experimental method with which to verify it in vivo. In such cases, mathematical modeling can be effective for verifying whether a locally generated force can propagate through an entire tissue and induce deformation. Moreover, the mathematical model sometimes provides potential mechanistic insight beyond the information obtained from biological experimental results. Here, we present two examples of modeling tissue morphogenesis driven by cell deformation or cell interaction. In the first example, a mathematical study on tissue-autonomous folding based on a two-dimensional vertex model revealed that active modulations of cell mechanics along the basal–lateral surface, in addition to the apical side, can induce tissue-fold formation. In the second example, by applying a two-dimensional vertex model in an apical plane, a novel mechanism of tissue flow caused by asymmetric cell interactions was discovered, which explained the mechanics behind the collective cellular movement observed during epithelial morphogenesis.

Theoretical tool bridging cell polarities with development of robust morphologies

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.

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.

Syncytial germline architecture is actively maintained by contraction of an internal actomyosin corset

Syncytial architecture is an evolutionarily-conserved feature of the germline of many species and plays a crucial role in their fertility. However, the mechanism supporting syncytial organization is largely unknown. Here, we identify a corset-like actomyosin structure within the syncytial germline of Caenorhabditis elegans, surrounding the common rachis. Using laser microsurgery, we demonstrate that actomyosin contractility within this structure generates tension both in the plane of the rachis surface and perpendicular to it, opposing membrane tension. Genetic and pharmacological perturbations, as well as mathematical modeling, reveal a balance of forces within the gonad and show how changing the tension within the actomyosin corset impinges on syncytial germline structure, leading, in extreme cases, to sterility. Thus, our work highlights a unique tissue-level cytoskeletal structure, and explains the critical role of actomyosin contractility in the preservation of a functional germline.

Germline cells in many species are fused to form a syncytium but the mechanics behind the maintenance of these structures are poorly defined. Here, the authors propose an inner contractile actomyosin corset provides a supportive framework to maintain germline architecture in C. elegans.

Fluidization of epithelial sheets by active cell rearrangements

We theoretically explore fluidization of epithelial tissues by active T1 neighbor exchanges. We show that the geometry of cell-cell junctions encodes important information about the local features of the energy landscape, which we support by an elastic theory of T1 transformations. Using a 3D vertex model, we show that the degree of active noise driving forced cell rearrangements governs the stress-relaxation timescale of the tissue. We study tissue response to in-plane shear at different timescales. At short time, the tissue behaves as a solid, whereas its long-time fluid behavior can be associated with an effective viscosity which scales with the rate of active T1 transformations. Furthermore, we develop a coarse-grained theory, where we treat the tissue as an active fluid and confirm the results of the vertex model. The impact of cell rearrangements on tissue shape is illustrated by studying axial compression of an epithelial tube.

Scutoids are a geometrical solution to three-dimensional packing of epithelia

As animals develop, tissue bending contributes to shape the organs into complex three-dimensional structures. However, the architecture and packing of curved epithelia remains largely unknown. Here we show by means of mathematical modelling that cells in bent epithelia can undergo intercalations along the apico-basal axis. This phenomenon forces cells to have different neighbours in their basal and apical surfaces. As a consequence, epithelial cells adopt a novel shape that we term "scutoid". The detailed analysis of diverse tissues confirms that generation of apico-basal intercalations between cells is a common feature during morphogenesis. Using biophysical arguments, we propose that scutoids make possible the minimization of the tissue energy and stabilize three-dimensional packing. Hence, we conclude that scutoids are one of nature's solutions to achieve epithelial bending. Our findings pave the way to understand the three-dimensional organization of epithelial organs.

Topological graph description of multicellular dynamics based on vertex model

Vertex models are generally powerful tools for exploring biological insights into multicellular dynamics. In these models, a multicellular structure is represented by a network, which is dynamically rearranged using topological operations. Remarkably, the topological dynamics of the network are important in guaranteeing the results from the models and their biological implications. However, it remains unclear whether the entire pattern of multicellular topological dynamics can be accurately expressed by a set of operators in the models. Surprisingly, vertex models have been empirically used for several decades without any mathematical verification. In this study, we propose a rigorous two-/three-dimensional (2D/3D) vertex model to describe multicellular topological dynamics. To do this, we classify several types of vertex models from a graph-theoretic perspective. Based on the classification, mathematical analyses reveal several conditions that enable us to apply the operators accurately without topological errors. Under these conditions, the operators can completely express the entire pattern of multicellular topological dynamics. From these results, we newly propose rigorous 2D/3D vertex models that can be applied to general multicellular dynamics, and we clarify several points to verify the results obtained from previous models.

Geometric constraints alter cell arrangements within curved epithelial tissues

Organ and tissue formation are complex three-dimensional processes involving cell division, growth, migration, and rearrangement, all of which occur within physically constrained regions. However, analyzing such processes in three dimensions in vivo is challenging. Here, we focus on the process of cellularization in the anterior pole of the early Drosophila embryo to explore how cells compete for space under geometric constraints. Using microfluidics combined with fluorescence microscopy, we extract quantitative information on the three-dimensional epithelial cell morphology. We observed a cellular membrane rearrangement in which cells exchange neighbors along the apical-basal axis. Such apical-to-basal neighbor exchanges were observed more frequently in the anterior pole than in the embryo trunk. Furthermore, cells within the anterior pole skewed toward the trunk along their long axis relative to the embryo surface, with maximum skew on the ventral side. We constructed a vertex model for cells in a curved environment. We could reproduce the observed cellular skew in both wild-type embryos and embryos with distorted morphology. Further, such modeling showed that cell rearrangements were more likely in ellipsoidal, compared with cylindrical, geometry. Overall, we demonstrate that geometric constraints can influence three-dimensional cell morphology and packing within epithelial tissues.

Complex structures from patterned cell sheets

The formation of three-dimensional structures from patterned epithelial sheets plays a key role in tissue morphogenesis. An important class of morphogenetic mechanisms relies on the spatio-temporal control of apical cell contractility, which can result in the localized bending of cell sheets and in-plane cell rearrangements. We have recently proposed a modified vertex model that can be used to systematically explore the connection between the two-dimensional patterns of cell properties and the emerging three-dimensional structures. Here we review the proposed modelling framework and illustrate it through the computational analysis of the vertex model that captures the salient features of the formation of the dorsal appendages during Drosophila oogenesis.This article is part of the themed issue 'Systems morphodynamics: understanding the development of tissue hardware'.

The physics of organoids: a biophysical approach to understanding organogenesis

Organoids representing a diversity of tissues have recently been created, bridging the gap between cell culture and experiments performed in vivo Being small and amenable to continuous monitoring, they offer the opportunity to scrutinize the dynamics of organ development, including the exciting prospect of observing aspects of human embryo development live. From a physicist's perspective, their ability to self-organize - to differentiate and organize cells in space - calls for the identification of the simple rules that underlie this capacity. Organoids provide tractable conditions to investigate the effects of the growth environment, including its molecular composition and mechanical properties, along with the initial conditions such as cell number and type(s). From a theoretical standpoint, different types of in silico modeling can complement the measurements performed in organoids to understand the role of chemical diffusion, contact signaling, differential cell adhesion and mechanical controls. Here, we discuss what it means to take a biophysical approach to understanding organogenesis in vitro and how we might expect such approaches to develop in the future.

Epithelial Folding Driven by Apical or Basal-Lateral Modulation: Geometric Features, Mechanical Inference, and Boundary Effects

During embryonic development, epithelial sheets fold into complex structures required for tissue and organ functions. Although substantial efforts have been devoted to identifying molecular mechanisms underlying epithelial folding, far less is understood about how forces deform individual cells to sculpt the overall sheet morphology. Here we describe a simple and general theoretical model for the autonomous folding of monolayered epithelial sheets. We show that active modulation of intracellular mechanics along the basal-lateral as well as the apical surfaces is capable of inducing fold formation in the absence of buckling instability. Apical modulation sculpts epithelia into shallow and V-shaped folds, whereas basal-lateral modulation generates deep and U-shaped folds. These characteristic tissue shapes remain unchanged when subject to mechanical perturbations from the surroundings, illustrating that the autonomous folding is robust against environmental variabilities. At the cellular scale, how cells change shape depends on their initial aspect ratios and the modulation mechanisms. Such cell deformation characteristics are verified via experimental measurements for a canonical folding process driven by apical modulation, indicating that our theory could be used to infer the underlying folding mechanisms based on experimental data. The mechanical principles revealed in our model could potentially guide future studies on epithelial folding in diverse systems.

Contractile actin belt and mesh structures provide the opposite dependence of epithelial stiffness on the spontaneous curvature of constituent cells

Actomyosin generates contractile forces within cells, which have a crucial role in determining the macroscopic mechanical properties of epithelial tissues. Importantly, actin cytoskeleton, which propagates actomyosin contractile forces, forms several characteristic structures in a 3D intracellular space, such as a circumferential actin belt lining adherence junctions and an actin mesh beneath the apical membrane. However, little is known about how epithelial mechanical property depends on the intracellular contractile structures. We performed computational simulations using a 3D vertex model, and demonstrated the longitudinal tensile test of an epithelial tube, whose inside and outside are defined as the apical and basal surfaces, respectively. As a result, these subcellular structures provide the contrary dependence of epithelial stiffness and fracture force on the spontaneous curvature of constituent cells; the epithelial stiffness increases with increasing the spontaneous curvature in the case of belt, meanwhile it decreases in the case of mesh. This qualitative difference emerges from the different anisotropic deformability of apical cell surfaces; while belt preserves isotropic apical cell shapes, mesh does not. Moreover, the difference in the anisotropic deformability determines the frequency of cell rearrangements, which in turn effectively decrease the tube stiffness. These results illustrate the importance of the intracellular contractile structures, which may be regulated to optimize mechanical functions of individual epithelial tissues.

Vertex models: from cell mechanics to tissue morphogenesis

Tissue morphogenesis requires the collective, coordinated motion and deformation of a large number of cells. Vertex model simulations for tissue mechanics have been developed to bridge the scales between force generation at the cellular level and tissue deformation and flows. We review here various formulations of vertex models that have been proposed for describing tissues in two and three dimensions. We discuss a generic formulation using a virtual work differential, and we review applications of vertex models to biological morphogenetic processes. We also highlight recent efforts to obtain continuum theories of tissue mechanics, which are effective, coarse-grained descriptions of vertex models.This article is part of the themed issue 'Systems morphodynamics: understanding the development of tissue hardware'.

Generating tissue topology through remodeling of cell-cell adhesions

During tissue morphogenesis, cellular rearrangements give rise to a large variety of three-dimensional structures. Final tissue architecture varies greatly across organs, and many develop to include combinations of folds, tubes, and branched networks. To achieve these different tissue geometries, constituent cells must follow different programs that dictate changes in shape and/or migratory behavior. One essential component of these changes is the remodeling of cell-cell adhesions. Invasive migratory behavior and separation between tissues require localized breakdown of cadherin-mediated adhesions. Conversely, tissue folding and fusion require the formation and reinforcement of cell-cell adhesions. Cell-cell adhesion plays a critical role in tissue morphogenesis; its manipulation may therefore prove to be invaluable in generating complex topologies ex vivo. Recapitulating these shapes in engineered tissues would enable a better understanding of how these processes occur in vivo, and may lead to improved design of organs for clinical applications. In this review, we discuss work investigating the formation of folds, tubes, and branched networks with an emphasis on known or possible roles for cell-cell adhesion. We then examine recently developed tools that could be adapted to manipulate cell-cell adhesion in engineered tissues.

Contact inhibition of locomotion and mechanical cross-talk between cell-cell and cell-substrate adhesion determine the pattern of junctional tension in epithelial cell aggregates

A computational approach is used to analyze the biomechanics of epithelial cells based on their capacity to adhere to one another and to the substrate and exhibit contact inhibition of locomotion. This approach reproduces emergent properties of epithelial cell aggregates and makes predictions for experimental validation.

We used a computational approach to analyze the biomechanics of epithelial cell aggregates—islands, stripes, or entire monolayers—that combines both vertex and contact-inhibition-of-locomotion models to include cell–cell and cell–substrate adhesion. Examination of the distribution of cell protrusions (adhesion to the substrate) in the model predicted high-order profiles of cell organization that agree with those previously seen experimentally. Cells acquired an asymmetric distribution of basal protrusions, traction forces, and apical aspect ratios that decreased when moving from the edge to the island center. Our in silico analysis also showed that tension on cell–cell junctions and apical stress is not homogeneous across the island. Instead, these parameters are higher at the island center and scale up with island size, which we confirmed experimentally using laser ablation assays and immunofluorescence. Without formally being a three-dimensional model, our approach has the minimal elements necessary to reproduce the distribution of cellular forces and mechanical cross-talk, as well as the distribution of principal stress in cells within epithelial cell aggregates. By making experimentally testable predictions, our approach can aid in mechanical analysis of epithelial tissues, especially when local changes in cell–cell and/or cell–substrate adhesion drive collective cell behavior.

The dynamic mechanical properties of cellularised aggregates

Cellularised materials are composed of cells interfaced through specialised intercellular junctions that link the cytoskeleton of one cell to that of its neighbours allowing for transmission of forces. Cellularised materials are common in early development and adult tissues where they can be found in the form of cell sheets, cysts, or amorphous aggregates and in pathophysiological conditions such as cancerous tumours. Given the growing realisation that forces can regulate cell physiology and developmental processes, understanding how cellularised materials deform under mechanical stress or dissipate stress appear as key biological questions. In this review, we will discuss the dynamic mechanical properties of cellularised materials devoid of extracellular matrix.