Computer simulation of wound closure in epithelial tissues: cell-basal-lamina adhesion

Free and Interfacial Boundaries in Individual-Based Models of Multicellular Biological systems

Coordination of cell behaviour is key to a myriad of biological processes including tissue morphogenesis, wound healing, and tumour growth. As such, individual-based computational models, which explicitly describe inter-cellular interactions, are commonly used to model collective cell dynamics. However, when using individual-based models, it is unclear how descriptions of cell boundaries affect overall population dynamics. In order to investigate this we define three cell boundary descriptions of varying complexities for each of three widely used off-lattice individual-based models: overlapping spheres, Voronoi tessellation, and vertex models. We apply our models to multiple biological scenarios to investigate how cell boundary description can influence tissue-scale behaviour. We find that the Voronoi tessellation model is most sensitive to changes in the cell boundary description with basic models being inappropriate in many cases. The timescale of tissue evolution when using an overlapping spheres model is coupled to the boundary description. The vertex model is demonstrated to be the most stable to changes in boundary description, though still exhibits timescale sensitivity. When using individual-based computational models one should carefully consider how cell boundaries are defined. To inform future work, we provide an exploration of common individual-based models and cell boundary descriptions in frequently studied biological scenarios and discuss their benefits and disadvantages.

Interacting active surfaces: A model for three-dimensional cell aggregates

We introduce a modelling and simulation framework for cell aggregates in three dimensions based on interacting active surfaces. Cell mechanics is captured by a physical description of the acto-myosin cortex that includes cortical flows, viscous forces, active tensions, and bending moments. Cells interact with each other via short-range forces capturing the effect of adhesion molecules. We discretise the model equations using a finite element method, and provide a parallel implementation in C++. We discuss examples of application of this framework to small and medium-sized aggregates: we consider the shape and dynamics of a cell doublet, a planar cell sheet, and a growing cell aggregate. This framework opens the door to the systematic exploration of the cell to tissue-scale mechanics of cell aggregates, which plays a key role in the morphogenesis of embryos and organoids.

Two types of critical cell density for mechanical elimination of abnormal cell clusters from epithelial tissue

Recent technological advances in high-resolution imaging and artificial modulation of genetic functions at different times and regions have enabled direct observations of the formation and elimination of abnormal cell populations. A recent trend in cell competition research is the incorporation of cell mechanics. In different tissues and species, abnormal cells developing in epithelial tissues are mechanically eliminated by cell contraction via actomyosin accumulation at the interface between normal and abnormal cells. This mechanical cell elimination process has attracted attention as a potential universal defense mechanism. Here, we theoretically examined the conditions for mechanical elimination of growing abnormal cell populations. Simulations and mathematical analyses using a vertex dynamics model revealed two types of critical cell density associated with mechanical elimination of abnormal cell clusters. One is a subtype of homeostatic density, in which the frequencies of spontaneous mechanical cell elimination and proliferation are balanced, even if no explicit dependence of proliferation or apoptosis on the cell density is assumed. This density is related to the mechanical stability of a single cell. The other is density related to mechanical stability as a cell population under external pressure. Both density types are determined by tissue mechanical properties. In solid tissues, the former type is reached first as the intensity of interfacial contraction increases, and it functions as a critical density. On the other hand, the latter type becomes critical when tissues are highly fluid. The derived analytical solution explicitly reveals the dependence of critical contractile force and density on different parameters. We also found a negative correlation between the proliferation rate of abnormal cells and the likelihood of the abnormal cell population expanding by escaping elimination. This is counterintuitive because in the context of cell competition, fast-growing cell populations generally win. These findings provide new insight into, and interpretation of, the results from experimental studies.

High-resolution imaging techniques have revealed that abnormal cells developing in epithelial tissues are mechanically eliminated via contraction at the interface between the abnormal cells and normal surrounding cells. This phenomenon is seen in various species and tissues and thus is regarded as a primitive defense system against precancerous cells common to all animals. For comprehensive understanding of this potential defense system, we derived mathematical conditions to achieve mechanical elimination of growing abnormal cell populations. We identified two characteristic cell density types associated with successful mechanical elimination of abnormal cell clusters. Both are determined by tissue physical properties, and the smaller of the two functions as a critical density above which abnormal cell populations cannot exist. We also found a counterintuitive phenomenon in which slower proliferation of abnormal cells promotes their growth as a population. Our results will help elucidate the mechanisms of intrinsic tissue defenses against cancer from the perspective of cell mechanics.

Activation of Topological Defects Induces a Brittle-to-Ductile Transition in Epithelial Monolayers

Epithelial monolayers are subjected to various mechanical forces, such as stretching, shearing, and compression. Thus, its mechanical response to external loadings is essential for its biological functions. However, the mechanism of the fracture failure of the epithelial monolayer remains poorly understood. Here, by introducing a new type of topological transition, i.e., detach transition or T4 transition, we develop a modified cellular vertex model to investigate the rupture of the cell monolayer. Interestingly, we find a brittle-to-ductile transition in epithelial monolayers, which is controlled by the mechanical properties of single cells and cell-cell contacts. We reveal that the external loadings can activate cell rearrangement in ductile cell monolayers. The plastic deformation results from the nucleation and propagation of "pentagon-heptagon defects" in analogy with the topological defects commonly seen in 2D materials. By using a simplified four-cell model, we further demonstrate that the brittle-to-ductile transition is induced by the competition between cell rearrangement and cell detachment. Our work provides a new theoretical framework to study the rupture of living tissues and may have important implications for many other biological processes, such as wound healing and tissue morphogenesis.

Large-scale demixing in a binary mixture of cells with rigidity disparity in biological tissues

Physical demixing on large scales of embryonic cell populations is fundamental to metazoan development, but whether a rigidity disparity alone is sufficient to driving large-scale demixing in a binary mixture of cell tissues is still an open question. To answer this question, we study mixing and demixing in a binary mixture of rigidity disparity cell tissues without heterotypic interactions using the Voronoi-based cellular model. Under suitable system parameters, the solid-like cells in the mixture can aggregate into a large cluster and the large-scale demixing occurs, which addresses that a rigidity disparity alone is sufficient to drive large-scale demixing. Remarkably, there exists an optimal temperature or rigidity disparity at which the binary mixture can be separated to the maximum extent. The necessary condition for the separation of mixtures is that the two types of cells are solid-like and liquid-like, respectively. The observation of robust demixing on large scales suggests that the sorting of progenitor cells may occur very early in the development process before robust heterotypic interfacial tensions are established. Our findings are relevant to understanding the mechanisms that drive cell sorting in confluent tissues.

A mesoscale mechanical model of cellular interactions

Computational models of cell mechanics allow the precise interrogation of cell shape change. These morphological changes are required for cells to survive in diverse tissue environments. Here, we present a mesoscale mechanical model of cell-substrate interactions using the level set method based on experimentally measured parameters. By implementing a viscoelastic mechanical equivalent circuit, we accurately model whole-cell deformations that are important for a variety of cellular processes. To effectively model shape changes as a cell interacts with a substrate, we have included receptor-mediated adhesion, which is governed by catch-slip bond behavior. The effect of adhesion was explored by subjecting cells to a variety of different substrates including flat, curved, and deformable surfaces. Finally, we increased the accuracy of our simulations by including a deformable nucleus in our cells. This model sets the foundation for further exploration into computational analyses of multicellular interactions.

Cell-substrate mechanics guide collective cell migration through intercellular adhesion: a dynamic finite element cellular model

During the process of tissue formation and regeneration, cells migrate collectively while remaining connected through intercellular adhesions. However, the roles of cell-substrate and cell-cell mechanical interactions in regulating collective cell migration are still unclear. In this study, we employ a newly developed finite element cellular model to study collective cell migration by exploring the effects of mechanical feedback between cell and substrate and mechanical signal transmission between adjacent cells. Our viscoelastic model of cells consists many triangular elements and is of high resolution. Cadherin adhesion between cells is modeled explicitly as linear springs at subcellular level. In addition, we incorporate a mechano-chemical feedback loop between cell-substrate mechanics and Rac-mediated cell protrusion. Our model can reproduce a number of experimentally observed patterns of collective cell migration during wound healing, including cell migration persistence, separation distance between cell pairs and migration direction. Moreover, we demonstrate that cell protrusion determined by the cell-substrate mechanics plays an important role in guiding persistent and oriented collective cell migration. Furthermore, this guidance cue can be maintained and transmitted to submarginal cells of long distance through intercellular adhesions. Our study illustrates that our finite element cellular model can be employed to study broad problems of complex tissue in dynamic changes at subcellular level.

Junctional and cytoplasmic contributions in wound healing

Wound healing is characterized by the re-epitheliation of a tissue through the activation of contractile forces concentrated mainly at the wound edge. While the formation of an actin purse string has been identified as one of the main mechanisms, far less is known about the effects of the viscoelastic properties of the surrounding cells, and the different contribution of the junctional and cytoplasmic contractilities. In this paper, we simulate the wound healing process, resorting to a hybrid vertex model that includes cell boundary and cytoplasmic contractilities explicitly, together with a differentiated viscoelastic rheology based on an adaptive rest-length. From experimental measurements of the recoil and closure phases of wounds in the Drosophila wing disc epithelium, we fit tissue viscoelastic properties. We then analyse in terms of closure rate and energy requirements the contributions of junctional and cytoplasmic contractilities. Our results suggest that reduction of junctional stiffness rather than cytoplasmic stiffness has a more pronounced effect on shortening closure times, and that intercalation rate has a minor effect on the stored energy, but contributes significantly to shortening the healing duration, mostly in the later stages.

Self-sustained planar intercalations due to mechanosignaling feedbacks lead to robust axis extension during morphogenesis

Tissue elongation is a necessary process in metazoans to implement their body plans that is not fully understood. Here we propose a mechanism based on the interplay between cellular mechanics and primordia patterning that results in self-sustained planar intercalations. Thus, we show that a location-dependent modulation of the mechanical properties of cells leads to robust axis extension. To illustrate the plausibility of this mechanism, we test it against different patterning models by means of computer simulations of tissues where we implemented mechano-signaling feedbacks. Our results suggest that robust elongation relies on a trade-off between cellular and tissue strains that is orchestrated through the cleavage orientation. In the particular context of axis extension in Turing-patterned tissues, we report that different directional cell activities cooperate synergetically to achieve elongation. Altogether, our findings help to understand how the axis extension phenomenon emerges from the dynamics of individual cells.

Mechanosensitive Junction Remodeling Promotes Robust Epithelial Morphogenesis

Morphogenesis of epithelial tissues requires tight spatiotemporal coordination of cell shape changes. In vivo, many tissue-scale shape changes are driven by pulsatile contractions of intercellular junctions, which are rectified to produce irreversible deformations. The functional role of this pulsatory ratchet and its mechanistic basis remain unknown. Here we combine theory and biophysical experiments to show that mechanosensitive tension remodeling of epithelial cell junctions promotes robust epithelial shape changes via ratcheting. Using optogenetic control of actomyosin contractility, we find that epithelial junctions show elastic behavior under low contractile stress, returning to their original lengths after contraction, but undergo irreversible deformation under higher magnitudes of contractile stress. Existing vertex-based models for the epithelium are unable to capture these results, with cell junctions displaying purely elastic or fluid-like behaviors, depending on the choice of model parameters. To describe the experimental results, we propose a modified vertex model with two essential ingredients for junction mechanics: thresholded tension remodeling and continuous strain relaxation. First, junctions must overcome a critical strain threshold to trigger tension remodeling, resulting in irreversible junction length changes. Second, there is a continuous relaxation of junctional strain that removes mechanical memory from the system. This enables pulsatile contractions to further remodel cell shape via mechanical ratcheting. Taken together, the combination of mechanosensitive tension remodeling and junctional strain relaxation provides a robust mechanism for large-scale morphogenesis.

DLITE Uses Cell-Cell Interface Movement to Better Infer Cell-Cell Tensions

Cell shapes and connectivities evolve over time as the colony changes shape or embryos develop. Shapes of intercellular interfaces are closely coupled with the forces resulting from actomyosin interactions, membrane tension, or cell-cell adhesions. Although it is possible to computationally infer cell-cell forces from a mechanical model of collective cell behavior, doing so for temporally evolving forces in a manner robust to digitization difficulties is challenging. Here, we introduce a method for dynamic local intercellular tension estimation (DLITE) that infers such evolution in temporal force with less sensitivity to digitization ambiguities or errors. This method builds upon previous work on single time points (cellular force-inference toolkit). We validate our method using synthetic geometries. DLITE's inferred cell colony tension evolutions correlate better with ground truth for these synthetic geometries as compared to tension values inferred from methods that consider each time point in isolation. We introduce cell connectivity errors, angle estimate errors, connection mislocalization, and connection topological changes to synthetic data and show that DLITE has reduced sensitivity to these conditions. Finally, we apply DLITE to time series of human-induced pluripotent stem cell colonies with endogenously expressed GFP-tagged zonulae occludentes-1. We show that DLITE offers improved stability in the inference of cell-cell tensions and supports a correlation between the dynamics of cell-cell forces and colony rearrangement.

Critical contractility and cell size for mechanical cell elimination from epithelial tissue

A defense mechanism in epithelial tissue can mechanically eliminate abnormal cells by contracting the cell boundary or area via actomyosin activity. From numerical simulations of the vertex dynamics model and approximate analytical solutions based on deviations from the ground state, here we derived general conditions for mechanical cell elimination (MCE) occurring via cell contraction. In particular, we found that MCE is realized by saddle-node bifurcation in a wide parameter range, and that the size of the eliminated cell is almost constant at the bifurcation point, suggesting the existence of an intrinsic threshold cell area for MCE.

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.

Cooperation of dual modes of cell motility promotes epithelial stress relaxation to accelerate wound healing

Collective cell migration in cohesive units is vital for tissue morphogenesis, wound repair, and immune response. While the fundamental driving forces for collective cell motion stem from contractile and protrusive activities of individual cells, it remains unknown how their balance is optimized to maintain tissue cohesiveness and the fluidity for motion. Here we present a cell-based computational model for collective cell migration during wound healing that incorporates mechanochemical coupling of cell motion and adhesion kinetics with stochastic transformation of active motility forces. We show that a balance of protrusive motility and actomyosin contractility is optimized for accelerating the rate of wound repair, which is robust to variations in cell and substrate mechanical properties. This balance underlies rapid collective cell motion during wound healing, resulting from a tradeoff between tension mediated collective cell guidance and active stress relaxation in the tissue.

Simulation of Cell Patterning Triggered by Cell Death and Differential Adhesion in Drosophila Wing

The Drosophila wing exhibits a well-ordered cell pattern, especially along the posterior margin, where hair cells are arranged in a zigzag pattern in the lateral view. Based on an experimental result observed during metamorphosis of Drosophila, we considered that a pattern of initial cells autonomously develops to the zigzag pattern through cell differentiation, intercellular communication, and cell death (apoptosis) and performed computer simulations of a cell-based model of vertex dynamics for tissues. The model describes the epithelial tissue as a monolayer cell sheet of polyhedral cells. Their vertices move according to equations of motion, minimizing the sum total of the interfacial and elastic energies of cells. The interfacial energy densities between cells are introduced consistently with an ideal zigzag cell pattern, extracted from the experimental result. The apoptosis of cells is modeled by gradually reducing their equilibrium volume to zero and by assuming that the hair cells prohibit neighboring cells from undergoing apoptosis. Based on experimental observations, we also assumed wing elongation along the proximal-distal axis. Starting with an initial cell pattern similar to the micrograph experimentally obtained just before apoptosis, we carried out the simulations according to the model mentioned above and successfully reproduced the ideal zigzag cell pattern. This elucidates a physical mechanism of patterning triggered by cell apoptosis theoretically and exemplifies, to our knowledge, a new framework to study apoptosis-induced patterning. We conclude that the zigzag cell pattern is formed by an autonomous communicative process among the participant cells.

Dynamic cellular finite-element method for modelling large-scale cell migration and proliferation under the control of mechanical and biochemical cues: a study of re-epithelialization

Computational modelling of cells can reveal insight into the mechanisms of the important processes of tissue development. However, current cell models have limitations and are challenged to model detailed changes in cellular shapes and physical mechanics when thousands of migrating and interacting cells need to be modelled. Here we describe a novel dynamic cellular finite-element model (DyCelFEM), which accounts for changes in cellular shapes and mechanics. It also models the full range of cell motion, from movements of individual cells to collective cell migrations. The transmission of mechanical forces regulated by intercellular adhesions and their ruptures are also accounted for. Intra-cellular protein signalling networks controlling cell behaviours are embedded in individual cells. We employ DyCelFEM to examine specific effects of biochemical and mechanical cues in regulating cell migration and proliferation, and in controlling tissue patterning using a simplified re-epithelialization model of wound tissue. Our results suggest that biochemical cues are better at guiding cell migration with improved directionality and persistence, while mechanical cues are better at coordinating collective cell migration. Overall, DyCelFEM can be used to study developmental processes when a large population of migrating cells under mechanical and biochemical controls experience complex changes in cell shapes and mechanics.

Possible roles of mechanical cell elimination intrinsic to growing tissues from the perspective of tissue growth efficiency and homeostasis

Cell competition is a phenomenon originally described as the competition between cell populations with different genetic backgrounds; losing cells with lower fitness are eliminated. With the progress in identification of related molecules, some reports described the relevance of cell mechanics during elimination. Furthermore, recent live imaging studies have shown that even in tissues composed of genetically identical cells, a non-negligible number of cells are eliminated during growth. Thus, mechanical cell elimination (MCE) as a consequence of mechanical cellular interactions is an unavoidable event in growing tissues and a commonly observed phenomenon. Here, we studied MCE in a genetically-homogeneous tissue from the perspective of tissue growth efficiency and homeostasis. First, we propose two quantitative measures, cell and tissue fitness, to evaluate cellular competitiveness and tissue growth efficiency, respectively. By mechanical tissue simulation in a pure population where all cells have the same mechanical traits, we clarified the dependence of cell elimination rate or cell fitness on different mechanical/growth parameters. In particular, we found that geometrical (specifically, cell size) and mechanical (stress magnitude) heterogeneities are common determinants of the elimination rate. Based on these results, we propose possible mechanical feedback mechanisms that could improve tissue growth efficiency and density/stress homeostasis. Moreover, when cells with different mechanical traits are mixed (e.g., in the presence of phenotypic variation), we show that MCE could drive a drastic shift in cell trait distribution, thereby improving tissue growth efficiency through the selection of cellular traits, i.e. intra-tissue "evolution". Along with the improvement of growth efficiency, cell density, stress state, and phenotype (mechanical traits) were also shown to be homogenized through growth. More theoretically, we propose a mathematical model that approximates cell competition dynamics, by which the time evolution of tissue fitness and cellular trait distribution can be predicted without directly simulating a cell-based mechanical model.

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'.

Mechanical model of geometric cell and topological algorithm for cell dynamics from single-cell to formation of monolayered tissues with pattern

Geometric and mechanical properties of individual cells and interactions among neighboring cells are the basis of formation of tissue patterns. Understanding the complex interplay of cells is essential for gaining insight into embryogenesis, tissue development, and other emerging behavior. Here we describe a cell model and an efficient geometric algorithm for studying the dynamic process of tissue formation in 2D (e.g. epithelial tissues). Our approach improves upon previous methods by incorporating properties of individual cells as well as detailed description of the dynamic growth process, with all topological changes accounted for. Cell size, shape, and division plane orientation are modeled realistically. In addition, cell birth, cell growth, cell shrinkage, cell death, cell division, cell collision, and cell rearrangements are now fully accounted for. Different models of cell-cell interactions, such as lateral inhibition during the process of growth, can be studied in detail. Cellular pattern formation for monolayered tissues from arbitrary initial conditions, including that of a single cell, can also be studied in detail. Computational efficiency is achieved through the employment of a special data structure that ensures access to neighboring cells in constant time, without additional space requirement. We have successfully generated tissues consisting of more than 20,000 cells starting from 2 cells within 1 hour. We show that our model can be used to study embryogenesis, tissue fusion, and cell apoptosis. We give detailed study of the classical developmental process of bristle formation on the epidermis of D. melanogaster and the fundamental problem of homeostatic size control in epithelial tissues. Simulation results reveal significant roles of solubility of secreted factors in both the bristle formation and the homeostatic control of tissue size. Our method can be used to study broad problems in monolayered tissue formation. Our software is publicly available.

Vertex models of epithelial morphogenesis

The dynamic behavior of epithelial cell sheets plays a central role during numerous developmental processes. Genetic and imaging studies of epithelial morphogenesis in a wide range of organisms have led to increasingly detailed mechanisms of cell sheet dynamics. Computational models offer a useful means by which to investigate and test these mechanisms, and have played a key role in the study of cell-cell interactions. A variety of modeling approaches can be used to simulate the balance of forces within an epithelial sheet. Vertex models are a class of such models that consider cells as individual objects, approximated by two-dimensional polygons representing cellular interfaces, in which each vertex moves in response to forces due to growth, interfacial tension, and pressure within each cell. Vertex models are used to study cellular processes within epithelia, including cell motility, adhesion, mitosis, and delamination. This review summarizes how vertex models have been used to provide insight into developmental processes and highlights current challenges in this area, including progressing these models from two to three dimensions and developing new tools for model validation.

Modelling wound closure in an epithelial cell sheet using the cellular Potts model

We use a two-dimensional cellular Potts model to represent the behavior of an epithelial cell layer and describe its dynamics in response to a microscopic wound.

Multiscale Modeling of Cellular Epigenetic States: Stochasticity in Molecular Networks, Chromatin Folding in Cell Nuclei, and Tissue Pattern Formation of Cells

Genome sequences provide the overall genetic blueprint of cells, but cells possessing the same genome can exhibit diverse phenotypes. There is a multitude of mechanisms controlling cellular epigenetic states and that dictate the behavior of cells. Among these, networks of interacting molecules, often under stochastic control, depending on the specific wirings of molecular components and the physiological conditions, can have a different landscape of cellular states. In addition, chromosome folding in three-dimensional space provides another important control mechanism for selective activation and repression of gene expression. Fully differentiated cells with different properties grow, divide, and interact through mechanical forces and communicate through signal transduction, resulting in the formation of complex tissue patterns. Developing quantitative models to study these multi-scale phenomena and to identify opportunities for improving human health requires development of theoretical models, algorithms, and computational tools. Here we review recent progress made in these important directions.

Cell models lead to understanding of multi-cellular morphogenesis consisting of successive self-construction of cells

Morphogenesis of multi-cellular organisms occurs through cell behaviours within a cell aggregate. Cell behaviours have been described using cell models involving equations of motion for cells. Cells in cell models construct shapes of the cell aggregate by themselves. Here, a history of cell models, the cell centre model and the vertex cell model, which we have constructed, are described. Furthermore, the application of these cell models is explained in detail. These cell models have been applied to transformation of cell aggregates to become spherical, formation of mammalian blastocysts and cell intercalation in elongating tissues. These are all elemental processes of morphogenesis and take place in succession during the whole developmental process. A chain of successive elemental processes leads to morphogenesis. Finally, we highlight that cell models are indispensable to understand the process whereby genes direct biological shapes.

Canalization-based vein formation in a growing leaf

Vein formation is an important process in plant leaf development. The phytohormone auxin is known as the most important molecule for the control of venation patterning; and the canalization model, in which cells experiencing higher auxin flux differentiate into specific cells for auxin transportation, is widely accepted. To date, several mathematical models based on the canalization hypothesis have been proposed that have succeeded in reproducing vein patterns similar to those observed in actual leaves. However, most previous studies focused on patterning in fixed domains, and, in a few exceptional studies, limited tissue growth - such as cell proliferation at leaf margins and small deformations without large changes in cell number - were dealt with. Considering that, in actual leaf development, venation patterning occurs in an exponentially growing tissue, whether the canalization hypothesis still applies is an important issue to be addressed. In this study, we first show through a pilot simulation that the coupling of chemical dynamics for canalization and tissue growth as independent models cannot reproduce normal venation patterning. We then examine conditions sufficient for achieving normal patterning in a growing leaf by introducing various constraints on chemical dynamics, tissue growth, and cell mechanics; in doing so, we found that auxin flux- or differentiation-dependent modification of the cell cycle and elasticity of cell edges are essential. The predictions given by our simulation study will serve as guideposts in experiments aimed at finding the key factors for achieving normal venation patterning in developing plant leaves.

A vertex model of Drosophila ventral furrow formation

Ventral furrow formation in Drosophila is an outstanding model system to study the mechanisms involved in large-scale tissue rearrangements. Ventral cells accumulate myosin at their apical sides and, while being tightly coupled to each other via apical adherens junctions, execute actomyosin contractions that lead to reduction of their apical cell surface. Thereby, a band of constricted cells along the ventral epithelium emerges which will form a tissue indentation along the ventral midline (the ventral furrow). Here we adopt a 2D vertex model to simulate ventral furrow formation in a surface view allowing easy comparison with confocal live-recordings. We show that in order to reproduce furrow morphology seen in vivo, a gradient of contractility must be assumed in the ventral epithelium which renders cells more contractile the closer they lie to the ventral midline. The model predicts previous experimental findings, such as the gain of eccentric morphology of constricting cells and an incremental fashion of apical cell area reduction. Analysis of the model suggests that this incremental area reduction is caused by the dynamical interplay of cell elasticity and stochastic contractility as well as by the opposing forces from contracting neighbour cells. We underpin results from the model through in vivo analysis of ventral furrow formation in wildtype and twi mutant embryos. Our results show that ventral furrow formation can be accomplished as a “tug-of-war” between stochastically contracting, mechanically coupled cells and may require less rigorous regulation than previously thought.

Computational and experimental study of the mechanics of embryonic wound healing

Wounds in the embryo show a remarkable ability to heal quickly without leaving a scar. Previous studies have found that an actomyosin ring (purse string) forms around the wound perimeter and contracts to close the wound over the course of several dozens of minutes. Here, we report experiments that reveal an even faster mechanism which remarkably closes wounds by more than 50% within the first 30s. Circular and elliptical wounds (~100μm in size) were made in the blastoderm of early chick embryos and allowed to heal, with wound area and shape characterized as functions of time. The closure rate displayed a biphasic behavior, with rapid constriction lasting about a minute, followed by a period of more gradual closure to complete healing. Fluorescent staining suggests that both healing phases are driven by actomyosin contraction, with relatively rapid contraction of fibers at cell borders within a relatively thick ring of tissue (several cells wide) around the wound followed by slower contraction of a thin supracellular actomyosin ring along the margin, consistent with a purse string mechanism. Finite-element modeling showed that this idea is biophysically plausible, with relatively isotropic contraction within the thick ring giving way to tangential contraction in the thin ring. In addition, consistent with experimental results, simulated elliptical wounds heal with little change in aspect ratio, and decreased membrane tension can cause these wounds to open briefly before going on to heal. These results provide new insight into the healing mechanism in embryonic epithelia.

Relaxation-expansion model for self-driven retinal morphogenesis: a hypothesis from the perspective of biosystems dynamics at the multi-cellular level

The generation of complex organ structures such as the eye requires the intricate orchestration of multiple cellular interactions. In this paper, early retinal development is discussed with respect to the structure formation of the optic cup. Although recent studies have elucidated molecular mechanisms of retinal differentiation, little is known about how the unique shape of the optic cup is determined. A recent report has demonstrated that optic-cup morphogenesis spontaneously occurs in three-dimensional stem-cell culture without external forces, indicating a latent intrinsic order to generate the structure. Based on this self-organizing phenomenon, we introduce the "relaxation-expansion" model to mechanically interpret the tissue dynamics that enable the spontaneous invagination of the neural retina. This model involves three consecutive local rules (relaxation, apical constriction, and expansion), and its computer simulation recapitulates the optic-cup morphogenesis in silico.

Dynamics and mechanical stability of the developing dorsoventral organizer of the wing imaginal disc

Shaping the primordia during development relies on forces and mechanisms able to control cell segregation. In the imaginal discs of Drosophila the cellular populations that will give rise to the dorsal and ventral parts on the wing blade are segregated and do not intermingle. A cellular population that becomes specified by the boundary of the dorsal and ventral cellular domains, the so-called organizer, controls this process. In this paper we study the dynamics and stability of the dorsal-ventral organizer of the wing imaginal disc of Drosophila as cell proliferation advances. Our approach is based on a vertex model to perform in silico experiments that are fully dynamical and take into account the available experimental data such as: cell packing properties, orientation of the cellular divisions, response upon membrane ablation, and robustness to mechanical perturbations induced by fast growing clones. Our results shed light on the complex interplay between the cytoskeleton mechanics, the cell cycle, the cell growth, and the cellular interactions in order to shape the dorsal-ventral organizer as a robust source of positional information and a lineage controller. Specifically, we elucidate the necessary and sufficient ingredients that enforce its functionality: distinctive mechanical properties, including increased tension, longer cell cycle duration, and a cleavage criterion that satisfies the Hertwig rule. Our results provide novel insights into the developmental mechanisms that drive the dynamics of the DV organizer and set a definition of the so-called Notch fence model in quantitative terms.

During development, tissues are shaped in order to form organs with specific functionalities. This process relies on mechanisms that control cell segregation and migration. These concepts are beautifully illustrated by the morphogenetic events that the imaginal discs of Drosophila undergo during metamorphosis. In particular, the cellular populations that will give rise to the dorsal (D) and ventral (V) parts on the wing blade are segregated and do not intermingle. The so-called organizer, a cellular population that becomes specified by the boundary of the D and V cellular domains, is responsible for this. Yet, how does the DV organizer robustly deal with the cellular growth in order to prevent cell mixing? Moreover, how can the organizer be conveniently scaled as the tissue grows? Herein we address these questions using a computational approach that takes into account the available experimental data. Thus, our study unveils the elements that are necessary and sufficient for understanding in a quantitative and predictive manner the dynamics, structure, and stability of a robust growing DV organizer: distinctive mechanical properties of cells, differences in cell cycle duration, and a well-defined cleavage criterion.

Self-organizing optic-cup morphogenesis in three-dimensional culture

Balanced organogenesis requires the orchestration of multiple cellular interactions to create the collective cell behaviours that progressively shape developing tissues. It is currently unclear how individual, localized parts are able to coordinate with each other to develop a whole organ shape. Here we report the dynamic, autonomous formation of the optic cup (retinal primordium) structure from a three-dimensional culture of mouse embryonic stem cell aggregates. Embryonic-stem-cell-derived retinal epithelium spontaneously formed hemispherical epithelial vesicles that became patterned along their proximal-distal axis. Whereas the proximal portion differentiated into mechanically rigid pigment epithelium, the flexible distal portion progressively folded inward to form a shape reminiscent of the embryonic optic cup, exhibited interkinetic nuclear migration and generated stratified neural retinal tissue, as seen in vivo. We demonstrate that optic-cup morphogenesis in this simple cell culture depends on an intrinsic self-organizing program involving stepwise and domain-specific regulation of local epithelial properties.

Mechanics and remodelling of cell packings in epithelia

Epithelia are sheets of cells that are dynamically remodelled by cell division and cell death during development. Here we describe the cell shapes and packings as networks of polygons: stable and stationary network configurations obey force balance and are represented as local minima of a potential function. We characterize the physical properties of this vertex model, including the set of ground states, and the energetics of topological rearrangements. We furthermore discuss a quasistatic description of cell division that allows us to study the mechanics and dynamics of tissue remodelling during growth. The biophysics of cells and their rearrangements can account for the morphology of cell packings observed in experiments.