Deconstructing gastrulation at single-cell resolution

Multicell-Fold: geometric learning in folding multicellular life

During developmental processes such as embryogenesis, how a group of cells fold into specific structures, is a central question in biology. However, it remains a major challenge to understand and predict the behavior of every cell within the living tissue over time during such intricate processes. Here we present a geometric deep-learning model that can accurately capture the highly convoluted interactions among cells. We demonstrate that multicellular data can be represented with both granular and foam-like physical pictures through a unified graph data structure, considering both cellular interactions and cell junction networks. Using this model, we achieve interpretable 4-D morphological sequence alignment, and predicting cell rearrangements before they occur at single-cell resolution. Furthermore, using neural activation map and ablation studies, we demonstrate cell geometries and cell junction networks together regulate morphogenesis at single-cell precision. This approach offers a pathway toward a unified dynamic atlas for a variety of developmental processes.

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.

The Geometric Basis of Epithelial Convergent Extension

Shape changes of epithelia during animal development, such as convergent extension, are achieved through concerted mechanical activity of individual cells. While much is known about the corresponding large scale tissue flow and its genetic drivers, fundamental questions regarding local control of contractile activity on cellular scale and its embryo-scale coordination remain open. To address these questions, we develop a quantitative, model-based analysis framework to relate cell geometry to local tension in recently obtained timelapse imaging data of gastrulating Drosophila embryos. This analysis provides a systematic decomposition of cell shape changes and T1-rearrangements into internally driven, active, and externally driven, passive, contributions. Our analysis provides evidence that germ band extension is driven by active T1 processes that self-organize through positive feedback acting on tensions. More generally, our findings suggest that epithelial convergent extension results from controlled transformation of internal force balance geometry which combines the effects of bottom-up local self-organization with the top-down, embryo-scale regulation by gene expression.

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.

EGFR-dependent actomyosin patterning coordinates morphogenetic movements between tissues

The movements that give rise to the body's structure are powered by cell shape changes and rearrangements that are coordinated at supracellular scales. How such cellular coordination arises and integrates different morphogenetic programs is unclear. Using quantitative imaging, we found a complex pattern of adherens junction (AJ) levels in the ectoderm prior to gastrulation onset in Drosophila. AJ intensity exhibited a double-sided gradient, with peaks at the dorsal midline and ventral neuroectoderm. We show that this dorsal-ventral AJ pattern is regulated by epidermal growth factor (EGF) signaling and that this signal is required for ectoderm cell movement during mesoderm invagination and axis extension. We identify AJ levels and junctional actomyosin as downstream effectors of EGFR signaling. Overall, our study demonstrates a mechanism of coordination between tissue folding and convergent extension that facilitates embryo-wide gastrulation movements.

The dynamics and biophysics of shape formation: Common themes in plant and animal morphogenesis

The emergence of tissue form in multicellular organisms results from the complex interplay between genetics and physics. In both plants and animals, cells must act in concert to pattern their behaviors. Our understanding of the factors sculpting multicellular form has increased dramatically in the past few decades. From this work, common themes have emerged that connect plant and animal morphogenesis-an exciting connection that solidifies our understanding of the developmental basis of multicellular life. In this review, we will discuss the themes and the underlying principles that connect plant and animal morphogenesis, including the coordination of gene expression, signaling, growth, contraction, and mechanical and geometric feedback.

Maintaining symmetry during body axis elongation

Bilateral symmetry defines much of the animal kingdom and is crucial for numerous functions of bilaterian organisms. Genetic approaches have discovered highly conserved patterning networks that establish bilateral symmetry in early embryos,1 but how this symmetry is maintained throughout subsequent morphogenetic events remains largely unknown.2 Here we show that the terminal patterning system-which relies on Ras/ERK signaling through activation of the Torso receptor by its ligand Trunk3-is critical for preserving bilateral symmetry during Drosophila body axis elongation, a process driven by cell rearrangements in the two identical lateral regions of the embryo and specified by the dorsal-ventral and anterior-posterior patterning systems.4 We demonstrate that fluctuating asymmetries in this rapid convergent-extension process are attenuated in normal embryos over time, possibly through noise-dissipating forces from the posterior midgut invagination and movement. However, when Torso signaling is attenuated via mutation of Trunk or RNAi directed against downstream Ras/ERK pathway components, body axis elongation results in a characteristic corkscrew phenotype,5 which reflects dramatic reorganization of global tissue flow and is incompatible with viability. Our results reveal a new function downstream of the Drosophila terminal patterning system in potentially active control of bilateral symmetry and should motivate systematic search for similar symmetry-preserving regulatory mechanisms in other bilaterians.

Growth orientations, rather than heterogeneous growth rates, dominate jaw joint morphogenesis in the larval zebrafish

In early limb embryogenesis, synovial joints acquire specific shapes which determine joint motion and function. The process by which the opposing cartilaginous joint surfaces are moulded into reciprocal and interlocking shapes, called joint morphogenesis, is one of the least understood aspects of joint formation and the cell-level dynamics underlying it are yet to be unravelled. In this research, we quantified key cellular dynamics involved in growth and morphogenesis of the zebrafish jaw joint and synthesised them in a predictive computational simulation of joint development. Cells in larval zebrafish jaw joints labelled with cartilage markers were tracked over a 48-h time window using confocal imaging. Changes in distance and angle between adjacent cell centroids resulting from cell rearrangement, volume expansion and extracellular matrix (ECM) deposition were measured and used to calculate the rate and direction of local tissue deformations. We observed spatially and temporally heterogeneous growth patterns with marked anisotropy over the developmental period assessed. There was notably elevated growth at the level of the retroarticular process of the Meckel's cartilage, a feature known to undergo pronounced shape changes during zebrafish development. Analysis of cell dynamics indicated a dominant role for cell volume expansion in growth, with minor influences from ECM volume increases and cell intercalation. Cell proliferation in the joint was minimal over the timeframe of interest. Synthesising the dynamic cell data into a finite element model of jaw joint development resulted in accurate shape predictions. Our biofidelic computational simulation demonstrated that zebrafish jaw joint growth can be reasonably approximated based on cell positional information over time, where cell positional information derives mainly from cell orientation and cell volume expansion. By modifying the input parameters of the simulation, we were able to assess the relative contributions of heterogeneous growth rates and of growth orientation. The use of uniform rather than heterogeneous growth rates only minorly impacted the shape predictions, whereas isotropic growth fields resulted in altered shape predictions. The simulation results suggest that growth anisotropy is the dominant influence on joint growth and morphogenesis. This study addresses the gap of the cellular processes underlying joint morphogenesis, with implications for understanding the aetiology of developmental joint disorders such as developmental dysplasia of the hip and arthrogryposis.

Morphogenetic forces planar polarize LGN/Pins in the embryonic head during Drosophila gastrulation

Spindle orientation is often achieved by a complex of Partner of Inscuteable (Pins)/LGN, Mushroom Body Defect (Mud)/Nuclear Mitotic Apparatus (NuMa), Gαi, and Dynein, which interacts with astral microtubules to rotate the spindle. Cortical Pins/LGN recruitment serves as a critical step in this process. Here, we identify Pins-mediated planar cell polarized divisions in several of the mitotic domains of the early Drosophila embryo. We found that neither planar cell polarity pathways nor planar polarized myosin localization determined division orientation; instead, our findings strongly suggest that Pins planar polarity and force generated from mesoderm invagination are important. Disrupting Pins polarity via overexpression of a myristoylated version of Pins caused randomized division angles. We found that disrupting forces through chemical inhibitors, depletion of an adherens junction protein, or blocking mesoderm invagination disrupted Pins planar polarity and spindle orientation. Furthermore, directional ablations that separated mesoderm from mitotic domains disrupted spindle orientation, suggesting that forces transmitted from mesoderm to mitotic domains can polarize Pins and orient division during gastrulation. To our knowledge, this is the first in vivo example where mechanical force has been shown to polarize Pins to mediate division orientation.