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.

A mechanical wave travels along a genetic guide to drive the formation of an epithelial furrow during Drosophila gastrulation.

Epithelial furrowing is a fundamental morphogenetic process during gastrulation, neurulation, and body shaping. A furrow often results from a fold that propagates along a line. How fold formation and propagation are controlled and driven is poorly understood. To shed light on this, we study the formation of the cephalic furrow, a fold that runs along the embryo dorsal-ventral axis during Drosophila gastrulation and the developmental role of which is still unknown. We provide evidence of its function and show that epithelial furrowing is initiated by a group of cells. This cellular cluster works as a pacemaker, triggering a bidirectional morphogenetic wave powered by actomyosin contractions and sustained by de novo medial apex-to-apex cell adhesion. The pacemaker's Cartesian position is under the crossed control of the anterior-posterior and dorsal-ventral gene patterning systems. Thus, furrow formation is driven by a mechanical trigger wave that travels under the control of a multidimensional genetic guide.

Wrinkling Instability in Unsupported Epithelial Sheets.

We investigate the elasticity of an unsupported epithelial monolayer and we discover that unlike a thin solid plate, which wrinkles if geometrically incompatible with the underlying substrate, the epithelium may do so even in the absence of the substrate. From a cell-based model, we derive an exact elasticity theory and discover wrinkling driven by the differential apico-basal surface tension. Our theory is mapped onto that for supported plates by introducing a phantom substrate whose stiffness is finite beyond a critical differential tension. This suggests a new mechanism for an autonomous control of tissues over the length scale of their surface patterns.

Active Instability and Nonlinear Dynamics of Cell-Cell Junctions.

Active cell-junction remodeling is important for tissue morphogenesis, yet its underlying physics is not understood. We study a mechanical model that describes junctions as dynamic active force dipoles. Their instability can trigger cell intercalations by a critical collapse. Nonlinearities in tissue's elastic response can stabilize the collapse either by a limit cycle or condensation of junction lengths at cusps of the energy landscape. Furthermore, active junction networks undergo collective instability to drive active in-plane ordering or develop a limit cycle of collective oscillations, which extends over regions of the energy landscape corresponding to distinct network topologies.

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.

Elasticity, Stability, and Quasioscillations of Cell-Cell Junctions in Solid Confluent Epithelia.

Macroscopic properties and shapes of biological tissues depend on the remodeling of cell-cell junctions at the microscopic scale. We propose a theoretical framework that couples a vertex model of solid confluent tissues with the dynamics describing generation of local force dipoles in the junctional actomyosin. Depending on the myosin turnover rate, junctions either preserve stable length or collapse to initiate cell rearrangements. We find that noise can amplify and sustain transient oscillations to the fixed point, giving rise to quasiperiodic junctional dynamics. We also discover that junctional stability is affected by cell arrangements and junctional rest tensions, which may explain junctional collapse during convergence and extension in embryos.

Publisher Correction: Mechanics of a multilayer epithelium instruct tumour architecture and function.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Mechanics of a multilayer epithelium instruct tumour architecture and function.

Loss of normal tissue architecture is a hallmark of oncogenic transformation1. In developing organisms, tissues architectures are sculpted by mechanical forces during morphogenesis2. However, the origins and consequences of tissue architecture during tumorigenesis remain elusive. In skin, premalignant basal cell carcinomas form 'buds', while invasive squamous cell carcinomas initiate as 'folds'. Here, using computational modelling, genetic manipulations and biophysical measurements, we identify the biophysical underpinnings and biological consequences of these tumour architectures. Cell proliferation and actomyosin contractility dominate tissue architectures in monolayer, but not multilayer, epithelia. In stratified epidermis, meanwhile, softening and enhanced remodelling of the basement membrane promote tumour budding, while stiffening of the basement membrane promotes folding. Additional key forces stem from the stratification and differentiation of progenitor cells. Tumour-specific suprabasal stiffness gradients are generated as oncogenic lesions progress towards malignancy, which we computationally predict will alter extensile tensions on the tumour basement membrane. The pathophysiologic ramifications of this prediction are profound. Genetically decreasing the stiffness of basement membranes increases membrane tensions in silico and potentiates the progression of invasive squamous cell carcinomas in vivo. Our findings suggest that mechanical forces-exerted from above and below progenitors of multilayered epithelia-function to shape premalignant tumour architectures and influence tumour progression.

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.

Solid-fluid transition and cell sorting in epithelia with junctional tension fluctuations.

Tissues transition between solid-like and fluid-like behavior, which has major implications for morphogenesis and disease. These transitions can occur due to changes in the intrinsic shape of constituent cells and cell motility. We consider an alternative mechanism by studying tissues that explore the energy landscape through stochastic dynamics, driven by turnover of junctional molecular motors. To identify the solid-fluid transition, we start with single-component tissues and show that the mean cell-shape index uniquely describes the effective diffusion coefficient of cell movements, which becomes finite at the transition. We generalize our approach to two-component tissues, and explore cell-sorting dynamics both due to differential adhesion and due to differential degree of junctional fluctuations. We recover some known characteristic scaling relations describing the sorting kinetics, and discover some discrepancies from these relations in the case of differential-fluctuations-driven sorting. Finally, we show that differential fluctuations efficiently sort two solid-like tissues with a fluid intercompartmental boundary.

Cell-Size Pleomorphism Drives Aberrant Clone Dispersal in Proliferating Epithelia.

As epithelial tissues develop, groups of cells related by descent tend to associate in clonal populations rather than dispersing within the cell layer. While this is frequently assumed to be a result of differential adhesion, precise mechanisms controlling clonal cohesiveness remain unknown. Here we employ computational simulations to modulate epithelial cell size in silico and show that junctions between small cells frequently collapse, resulting in clone-cell dispersal among larger neighbors. Consistent with similar dynamics in vivo, we further demonstrate that mosaic disruption of Drosophila Tor generates small cells and results in aberrant clone dispersal in developing wing disc epithelia. We propose a geometric basis for this phenomenon, supported in part by the observation that soap-foam cells exhibit similar size-dependent junctional rearrangements. Combined, these results establish a link between cell-size pleomorphism and the control of epithelial cell packing, with potential implications for understanding tumor cell dispersal in human disease.

Self-Similar Dynamics of Nuclear Packing in the Early Drosophila Embryo.

Embryonic development starts with cleavages, a rapid sequence of reductive divisions that result in an exponential increase of cell number without changing the overall size of the embryo. In Drosophila, the final four rounds of cleavages occur at the surface of the embryo and give rise to ∼6000 nuclei under a common plasma membrane. We use live imaging to study the dynamics of this process and to characterize the emergent nuclear packing in this system. We show that the characteristic length scale of the internuclear interaction scales with the density, which allows the densifying embryo to sustain the level of structural order at progressively smaller length scales. This is different from nonliving materials, which typically undergo disorder-order transition upon compression. To explain this dynamics, we use a particle-based model that accounts for density-dependent nuclear interactions and synchronous divisions. We reproduce the pair statistics of the disordered packings observed in embryos and recover the scaling relation between the characteristic length scale and the density both in real and reciprocal space. This result reveals how the embryo can robustly preserve the nuclear-packing structure while being densified. In addition to providing quantitative description of self-similar dynamics of nuclear packings, this model generates dynamic meshes for the computational analysis of pattern formation and tissue morphogenesis.

Nuclear (Bio)physics in the Embryo.

Deneke et al. (2019) discover that dynamic interactions of cell cycle and actomyosin contractility systems synchronize nuclear cleavages, generating a cytoplasmic flow that results in a spatially uniform distribution of zygotic nuclei in the early Drosophila embryo. This work underscores the importance of self-organizing mechanisms before the onset of zygotic transcription.

Metabolic Regulation of Developmental Cell Cycles and Zygotic Transcription.

The thirteen nuclear cleavages that give rise to the Drosophila blastoderm are some of the fastest known cell cycles [1]. Surprisingly, the fertilized egg is provided with at most one-third of the dNTPs needed to complete the thirteen rounds of DNA replication [2]. The rest must be synthesized by the embryo, concurrent with cleavage divisions. What is the reason for the limited supply of DNA building blocks? We propose that frugal control of dNTP synthesis contributes to the well-characterized deceleration of the cleavage cycles and is needed for robust accumulation of zygotic gene products. In support of this model, we demonstrate that when the levels of dNTPs are abnormally high, nuclear cleavages fail to sufficiently decelerate, the levels of zygotic transcription are dramatically reduced, and the embryo catastrophically fails early in gastrulation. Our work reveals a direct connection between metabolism, the cell cycle, and zygotic transcription.

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.

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.

Quantitative Morphology of Epithelial Folds.

The shape of spatially modulated epithelial morphologies such as villi and crypts is usually associated with the epithelium-stroma area mismatch leading to buckling. We propose an alternative mechanical model based on intraepithelial stresses generated by differential tensions of apical, lateral, and basal sides of cells as well as on the elasticity of the basement membrane. We use it to theoretically study longitudinal folds in simple epithelia and we identify four types of corrugated morphologies: compact, invaginated, evaginated, and wavy. The obtained tissue contours and thickness profiles are compared to epithelial folds observed in invertebrates and vertebrates, and for most samples, the agreement is within the estimated experimental error. Our model establishes the groove-crest modulation of tissue thickness as a morphometric parameter that can, together with the curvature profile, be used to estimate the relative differential apicobasal tension in the epithelium.

Theory of epithelial elasticity.

We propose an elastic theory of epithelial monolayers based on a two-dimensional discrete model of dropletlike cells characterized by differential surface tensions of their apical, basal, and lateral sides. We show that the effective tissue bending modulus depends on the apicobasal differential tension and changes sign at the transition from the flat to the fold morphology. We discuss three mechanisms that stabilize the finite-wavelength fold structures: Physical constraint on cell geometry, hard-core interaction between non-neighboring cells, and bending elasticity of the basement membrane. We show that the thickness of the monolayer changes along the waveform and thus needs to be considered as a variable rather than a parameter. Next we show that the coupling between the curvature and the thickness is governed by the apicobasal polarity and that the amplitude of thickness modulation along the waveform is proportional to the apicobasal differential tension. This suggests that intracellular stresses can be measured indirectly by observing easily measurable morphometric parameters. We also study the mechanics of three-dimensional structures with cylindrical symmetry.

Embryo-scale tissue mechanics during Drosophila gastrulation movements.

Morphogenesis of an organism requires the development of its parts to be coordinated in time and space. While past studies concentrated on defined cell populations, a synthetic view of the coordination of these events in a whole organism is needed for a full understanding. Drosophila gastrulation begins with the embryo forming a ventral furrow, which is eventually internalized. It is not understood how the rest of the embryo participates in this process. Here we use multiview selective plane illumination microscopy coupled with infrared laser manipulation and mutant analysis to dissect embryo-scale cell interactions during early gastrulation. Lateral cells have a denser medial-apical actomyosin network and shift ventrally as a compact cohort, whereas dorsal cells become stretched. We show that the behaviour of these cells affects furrow internalization. A computational model predicts different mechanical properties associated with tissue behaviour: lateral cells are stiff, whereas dorsal cells are soft. Experimental analysis confirms these properties in vivo.