Force networks, torque balance and Airy stress in the planar vertex model of a confluent epithelium

A Geometric Tension Dynamics Model of Epithelial Convergent Extension

Epithelial tissue elongation by convergent extension is a key motif of animal morphogenesis. On a coarse scale, cell motion resembles laminar fluid flow; yet in contrast to a fluid, epithelial cells adhere to each other and maintain the tissue layer under actively generated internal tension. To resolve this apparent paradox, we formulate a model in which tissue flow in the tension-dominated regime occurs through adiabatic remodeling of force balance in the network of adherens junctions. We propose that the slow dynamics within the manifold of force-balanced configurations is driven by positive feedback on myosin-generated cytoskeletal tension. Shifting force balance within a tension network causes active cell rearrangements (T1 transitions) resulting in net tissue deformation oriented by initial tension anisotropy. Strikingly, we find that the total extent of tissue deformation depends on the initial cellular packing order. T1s degrade this order so that tissue flow is self-limiting. We explain these findings by showing that coordination of T1s depends on coherence in local tension configurations, quantified by a geometric order parameter in tension space. Our model reproduces the salient tissue- and cell-scale features of germ band elongation during Drosophila gastrulation, in particular the slowdown of tissue flow after approximately twofold elongation concomitant with a loss of order in tension configurations. This suggests local cell geometry contains morphogenetic information and yields experimentally testable predictions. Furthermore, defining biologically controlled active tension dynamics on the manifold of force-balanced states may provide a general approach to the description of morphogenetic flow.

Couple stresses and discrete potentials in the vertex model of cellular monolayers

The vertex model is widely used to simulate the mechanical properties of confluent epithelia and other multicellular tissues. This inherently discrete framework allows a Cauchy stress to be attributed to each cell, and its symmetric component has been widely reported, at least for planar monolayers. Here, we consider the stress attributed to the neighbourhood of each tricellular junction, evaluating in particular its leading-order antisymmetric component and the associated couple stresses, which characterise the degree to which individual cells experience (and resist) in-plane bending deformations. We develop discrete potential theory for localised monolayers having disordered internal structure and use this to derive the analogues of Airy and Mindlin stress functions. These scalar potentials typically have broad-banded spectra, highlighting the contributions of small-scale defects and boundary layers to global stress patterns. An affine approximation attributes couple stresses to pressure differences between cells sharing a trijunction, but simulations indicate an additional role for non-affine deformations.

How dynamic prestress governs the shape of living systems, from the subcellular to tissue scale

Cells and tissues change shape both to carry out their function and during pathology. In most cases, these deformations are driven from within the systems themselves. This is permitted by a range of molecular actors, such as active crosslinkers and ion pumps, whose activity is biologically controlled in space and time. The resulting stresses are propagated within complex and dynamical architectures like networks or cell aggregates. From a mechanical point of view, these effects can be seen as the generation of prestress or prestrain, resulting from either a contractile or growth activity. In this review, we present this concept of prestress and the theoretical tools available to conceptualize the statics and dynamics of living systems. We then describe a range of phenomena where prestress controls shape changes in biopolymer networks (especially the actomyosin cytoskeleton and fibrous tissues) and cellularized tissues. Despite the diversity of scale and organization, we demonstrate that these phenomena stem from a limited number of spatial distributions of prestress, which can be categorized as heterogeneous, anisotropic or differential. We suggest that in addition to growth and contraction, a third type of prestress-topological prestress-can result from active processes altering the microstructure of tissue.

Derivation of continuum models from discrete models of mechanical forces in cell populations

In certain discrete models of populations of biological cells, the mechanical forces between the cells are center based or vertex based on the microscopic level where each cell is individually represented. The cells are circular or spherical in a center based model and polygonal or polyhedral in a vertex based model. On a higher, macroscopic level, the time evolution of the density of the cells is described by partial differential equations (PDEs). We derive relations between the modelling on the micro and macro levels in one, two, and three dimensions by regarding the micro model as a discretization of a PDE for conservation of mass on the macro level. The forces in the micro model correspond on the macro level to a gradient of the pressure scaled by quantities depending on the cell geometry. The two levels of modelling are compared in numerical experiments in one and two dimensions.