Learning a conserved mechanism for early neuroectoderm morphogenesis.

Morphogenesis is the process whereby the body of an organism develops its target shape. The morphogen BMP is known to play a conserved role across bilaterian organisms in determining the dorsoventral (DV) axis. Yet, how BMP governs the spatio-temporal dynamics of cytoskeletal proteins driving morphogenetic flow remains an open question. Here, we use machine learning to mine a morphodynamic atlas of Drosophila development, and construct a mathematical model capable of predicting the coupled dynamics of myosin, E-cadherin, and morphogenetic flow. Mutant analysis shows that BMP sets the initial condition of this dynamical system according to the following signaling cascade: BMP establishes DV pair-rule-gene patterns that set-up an E-cadherin gradient which in turn creates a myosin gradient in the opposite direction through mechanochemical feedbacks. Using neural tube organoids, we argue that BMP, and the signaling cascade it triggers, prime the conserved dynamics of neuroectoderm morphogenesis from fly to humans.

YAP dysregulation triggers hypertrophy by CCN2 secretion and TGFß uptake in human pluripotent stem cell-derived cardiomyocytes.

Hypertrophy Cardiomyopathy (HCM) is the most prevalent hereditary cardiovascular disease - affecting >1:500 individuals. Advanced forms of HCM clinically present with hypercontractility, hypertrophy and fibrosis. Several single-point mutations in b-myosin heavy chain (MYH7) have been associated with HCM and increased contractility at the organ level. Different MYH7 mutations have resulted in increased, decreased, or unchanged force production at the molecular level. Yet, how these molecular kinetics link to cell and tissue pathogenesis remains unclear. The Hippo Pathway, specifically its effector molecule YAP, has been demonstrated to be reactivated in pathological hypertrophic growth. We hypothesized that changes in force production (intrinsically or extrinsically) directly alter the homeostatic mechano-signaling of the Hippo pathway through changes in stresses on the nucleus. Using human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), we asked whether homeostatic mechanical signaling through the canonical growth regulator, YAP, is altered 1) by changes in the biomechanics of HCM mutant cardiomyocytes and 2) by alterations in the mechanical environment. We use genetically edited hiPSC-CM with point mutations in MYH7 associated with HCM, and their matched controls, combined with micropatterned traction force microscopy substrates to confirm the hypercontractile phenotype in MYH7 mutants. We next modulate contractility in healthy and disease hiPSC-CMs by treatment with positive and negative inotropic drugs and demonstrate a correlative relationship between contractility and YAP activity. We further demonstrate the activation of YAP in both HCM mutants and healthy hiPSC-CMs treated with contractility modulators is through enhanced nuclear deformation. We conclude that the overactivation of YAP, possibly initiated and driven by hypercontractility, correlates with excessive CCN2 secretion (connective tissue growth factor), enhancing cardiac fibroblast/myofibroblast transition and production of known hypertrophic signaling molecule TGFß. Our study suggests YAP being an indirect player in the initiation of hypertrophic growth and fibrosis in HCM. Our results provide new insights into HCM progression and bring forth a testbed for therapeutic options in treating HCM.

Active cell divisions generate fourfold orientationally ordered phase in living tissue.

Morphogenesis, the process through which genes generate form, establishes tissue-scale order as a template for constructing the complex shapes of the body plan. The extensive growth required to build these ordered substrates is fuelled by cell proliferation, which, naively, should destroy order. Understanding how active morphogenetic mechanisms couple cellular and mechanical processes to generate order-rather than annihilate it-remains an outstanding question in animal development. We show that cell divisions are the primary drivers of tissue flow, leading to a fourfold orientationally ordered phase. Waves of anisotropic cell proliferation propagate across the embryo with precise patterning. Defects introduced into the nascent lattice by cell divisions are moved out of the tissue bulk towards the boundary by subsequent divisions. Specific cell proliferation rates and orientations enable cell divisions to organize rather than fluidize the tissue. We observe this using live imaging and tissue cartography to analyse the dynamics of fourfold tissue ordering in the trunk segmental ectoderm of the crustacean Parhyale hawaiensis beginning 72 h after egg lay. The result is a robust, active mechanism for generating global orientational order in a non-equilibrium system that sets the stage for the subsequent development of shape and form.

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.

Septins regulate border cell surface geometry, shape, and motility downstream of Rho in Drosophila.

Septins self-assemble into polymers that bind and deform membranes in vitro and regulate diverse cell behaviors in vivo. How their in vitro properties relate to their in vivo functions is under active investigation. Here, we uncover requirements for septins in detachment and motility of border cell clusters in the Drosophila ovary. Septins and myosin colocalize dynamically at the cluster periphery and share phenotypes but, surprisingly, do not impact each other. Instead, Rho independently regulates myosin activity and septin localization. Active Rho recruits septins to membranes, whereas inactive Rho sequesters septins in the cytoplasm. Mathematical analyses identify how manipulating septin expression levels alters cluster surface texture and shape. This study shows that the level of septin expression differentially regulates surface properties at different scales. This work suggests that downstream of Rho, septins tune surface deformability while myosin controls contractility, the combination of which governs cluster shape and movement.

Geometric control of myosin II orientation during axis elongation.

The actomyosin cytoskeleton is a crucial driver of morphogenesis. Yet how the behavior of large-scale cytoskeletal patterns in deforming tissues emerges from the interplay of geometry, genetics, and mechanics remains incompletely understood. Convergent extension in Drosophila melanogaster embryos provides the opportunity to establish a quantitative understanding of the dynamics of anisotropic non-muscle myosin II. Cell-scale analysis of protein localization in fixed embryos suggests that gene expression patterns govern myosin anisotropy via complex rules. However, technical limitations have impeded quantitative and dynamic studies of this process at the whole embryo level, leaving the role of geometry open. Here, we combine in toto live imaging with quantitative analysis of molecular dynamics to characterize the distribution of myosin anisotropy and the corresponding genetic patterning. We found pair rule gene expression continuously deformed, flowing with the tissue frame. In contrast, myosin anisotropy orientation remained approximately static and was only weakly deflected from the stationary dorsal-ventral axis of the embryo. We propose that myosin is recruited by a geometrically defined static source, potentially related to the embryo-scale epithelial tension, and account for transient deflections by cytoskeletal turnover and junction reorientation by flow. With only one parameter, this model quantitatively accounts for the time course of myosin anisotropy orientation in wild-type, twist, and even-skipped embryos, as well as embryos with perturbed egg geometry. Geometric patterning of the cytoskeleton suggests a simple physical strategy to ensure a robust flow and formation of shape.

Learning a conserved mechanism for early neuroectoderm morphogenesis.

Morphogenesis is the process whereby the body of an organism develops its target shape. The morphogen BMP is known to play a conserved role across bilaterian organisms in determining the dorsoventral (DV) axis. Yet, how BMP governs the spatio-temporal dynamics of cytoskeletal proteins driving morphogenetic flow remains an open question. Here, we use machine learning to mine a morphodynamic atlas of Drosophila development, and construct a mathematical model capable of predicting the coupled dynamics of myosin, E-cadherin, and morphogenetic flow. Mutant analysis shows that BMP sets the initial condition of this dynamical system according to the following signaling cascade: BMP establishes DV pair-rule-gene patterns that set-up an E-cadherin gradient which in turn creates a myosin gradient in the opposite direction through mechanochemical feedbacks. Using neural tube organoids, we argue that BMP, and the signaling cascade it triggers, prime the conserved dynamics of neuroectoderm morphogenesis from fly to humans.

Patterned mechanical feedback establishes a global myosin gradient.

Morphogenesis, the coordinated execution of developmental programs that shape embryos, raises many fundamental questions at the interface between physics and biology. In particular, how the dynamics of active cytoskeletal processes are coordinated across the surface of entire embryos to generate global cell flows is poorly understood. Two distinct regulatory principles have been identified: genetic programs and dynamic response to mechanical stimuli. Despite progress, disentangling these two contributions remains challenging. Here, we combine in toto light sheet microscopy with genetic and optogenetic perturbations of tissue mechanics to examine theoretically predicted dynamic recruitment of non-muscle myosin II to cell junctions during Drosophila embryogenesis. We find dynamic recruitment has a long-range impact on global myosin configuration, and the rate of junction deformation sets the rate of myosin recruitment. Mathematical modeling and high frequency analysis reveal myosin fluctuations on junctions around a mean value set by mechanical feedback. Our model accounts for the early establishment of the global myosin pattern at 80% fidelity. Taken together our results indicate spatially modulated mechanical feedback as a key regulatory input in the establishment of long-range gradients of cytoskeletal configurations and global tissue flow patterns.

Visceral organ morphogenesis via calcium-patterned muscle constrictions.

Organ architecture is often composed of multiple laminar tissues arranged in concentric layers. During morphogenesis, the initial geometry of visceral organs undergoes a sequence of folding, adopting a complex shape that is vital for function. Genetic signals are known to impact form, yet the dynamic and mechanical interplay of tissue layers giving rise to organs' complex shapes remains elusive. Here, we trace the dynamics and mechanical interactions of a developing visceral organ across tissue layers, from subcellular to organ scale in vivo. Combining deep tissue light-sheet microscopy for in toto live visualization with a novel computational framework for multilayer analysis of evolving complex shapes, we find a dynamic mechanism for organ folding using the embryonic midgut of Drosophila as a model visceral organ. Hox genes, known regulators of organ shape, control the emergence of high-frequency calcium pulses. Spatiotemporally patterned calcium pulses trigger muscle contractions via myosin light chain kinase. Muscle contractions, in turn, induce cell shape change in the adjacent tissue layer. This cell shape change collectively drives a convergent extension pattern. Through tissue incompressibility and initial organ geometry, this in-plane shape change is linked to out-of-plane organ folding. Our analysis follows tissue dynamics during organ shape change in vivo, tracing organ-scale folding to a high-frequency molecular mechanism. These findings offer a mechanical route for gene expression to induce organ shape change: genetic patterning in one layer triggers a physical process in the adjacent layer - revealing post-translational mechanisms that govern shape change.

Human neural tube morphogenesis in vitro by geometric constraints.

Understanding human organ formation is a scientific challenge with far-reaching medical implications1,2. Three-dimensional stem-cell cultures have provided insights into human cell differentiation3,4. However, current approaches use scaffold-free stem-cell aggregates, which develop non-reproducible tissue shapes and variable cell-fate patterns. This limits their capacity to recapitulate organ formation. Here we present a chip-based culture system that enables self-organization of micropatterned stem cells into precise three-dimensional cell-fate patterns and organ shapes. We use this system to recreate neural tube folding from human stem cells in a dish. Upon neural induction5,6, neural ectoderm folds into a millimetre-long neural tube covered with non-neural ectoderm. Folding occurs at 90% fidelity, and anatomically resembles the developing human neural tube. We find that neural and non-neural ectoderm are necessary and sufficient for folding morphogenesis. We identify two mechanisms drive folding: (1) apical contraction of neural ectoderm, and (2) basal adhesion mediated via extracellular matrix synthesis by non-neural ectoderm. Targeting these two mechanisms using drugs leads to morphological defects similar to neural tube defects. Finally, we show that neural tissue width determines neural tube shape, suggesting that morphology along the anterior-posterior axis depends on neural ectoderm geometry in addition to molecular gradients7. Our approach provides a new route to the study of human organ morphogenesis in health and disease.

Tissue topography steers migrating Drosophila border cells.

Moving cells can sense and respond to physical features of the microenvironment; however, in vivo, the significance of tissue topography is mostly unknown. Here, we used Drosophila border cells, an established model for in vivo cell migration, to study how chemical and physical information influences path selection. Although chemical cues were thought to be sufficient, live imaging, genetics, modeling, and simulations show that microtopography is also important. Chemoattractants promote predominantly posterior movement, whereas tissue architecture presents orthogonal information, a path of least resistance concentrated near the center of the egg chamber. E-cadherin supplies a permissive haptotactic cue. Our results provide insight into how cells integrate and prioritize topographical, adhesive, and chemoattractant cues to choose one path among many.

Dynamic morphoskeletons in development.

Morphogenetic flows in developmental biology are characterized by the coordinated motion of thousands of cells that organize into tissues, naturally raising the question of how this collective organization arises. Using only the kinematics of tissue deformation, which naturally integrates local and global mechanisms along cell paths, we identify the dynamic morphoskeletons behind morphogenesis, i.e., the evolving centerpieces of multicellular trajectory patterns. These features are model- and parameter-free, frame-invariant, and robust to measurement errors and can be computed from unfiltered cell-velocity data. We reveal the spatial attractors and repellers of the embryo by quantifying its Lagrangian deformation, information that is inaccessible to simple trajectory inspection or Eulerian methods that are local and typically frame-dependent. Computing these dynamic morphoskeletons in wild-type and mutant chick and fly embryos, we find that they capture the early footprint of known morphogenetic features, reveal new ones, and quantitatively distinguish between different phenotypes.

Topological structure and dynamics of three-dimensional active nematics.

Topological structures are effective descriptors of the nonequilibrium dynamics of diverse many-body systems. For example, motile, point-like topological defects capture the salient features of two-dimensional active liquid crystals composed of energy-consuming anisotropic units. We dispersed force-generating microtubule bundles in a passive colloidal liquid crystal to form a three-dimensional active nematic. Light-sheet microscopy revealed the temporal evolution of the millimeter-scale structure of these active nematics with single-bundle resolution. The primary topological excitations are extended, charge-neutral disclination loops that undergo complex dynamics and recombination events. Our work suggests a framework for analyzing the nonequilibrium dynamics of bulk anisotropic systems as diverse as driven complex fluids, active metamaterials, biological tissues, and collections of robots or organisms.

A variational method for image-based inference of internal stress in epithelial tissues.

Cellular mechanics drives epithelial morphogenesis, the process wherein cells collectively rearrange to produce tissue-scale deformations that determine organismal shape. However, quantitative understanding of tissue mechanics is impaired by the difficulty of direct measurement of stress in-vivo. This difficulty has spurred the development of image-based inference algorithms that estimate stress from snapshots of epithelial geometry. Such methods are challenged by sensitivity to measurement error and thus require accurate geometric segmentation for practical use. We overcome this difficulty by introducing a novel approach - the Variational Method of Stress Inference (VMSI) - which exploits the fundamental duality between stress and geometry at equilibrium of discrete mechanical networks that model confluent cellular layers. We approximate the apical geometry of an epithelial tissue by a 2D tiling with Circular Arc Polygons (CAP) in which arcs represent intercellular interfaces defined by the balance of local line tension and pressure differentials between adjacent cells. The mechanical equilibrium of such networks imposes extensive local constraints on CAP geometry. These constraints provide the foundation of VMSI which, starting with images of epithelial monolayers, simultaneously approximates both tissue geometry and internal forces, subject to the constraint of equilibrium. We find VMSI to be more robust than previous methods. Specifically, the VMSI performance is validated by the comparison of the predicted cellular and mesoscopic scale stress with the measured myosin II patterns during early Drosophila embryogenesis. VMSI prediction of mesoscopic stress tensor correlates at the 80% level with the measured myosin distribution and reveals that most of the myosin activity in that case is involved in a static internal force balance within the epithelial layer. In addition to insight into cell mechanics, this study provides a practical method for non-destructive estimation of stress in live epithelial tissue.

Extracellular matrix stiffness cues junctional remodeling for 3D tissue elongation.

Organs are sculpted by extracellular as well as cell-intrinsic forces, but how collective cell dynamics are orchestrated in response to environmental cues is poorly understood. Here we apply advanced image analysis to reveal extracellular matrix-responsive cell behaviors that drive elongation of the Drosophila follicle, a model system in which basement membrane stiffness instructs three-dimensional tissue morphogenesis. Through in toto morphometric analyses of wild type and round egg mutants, we find that neither changes in average cell shape nor oriented cell division are required for appropriate organ shape. Instead, a major element is the reorientation of elongated cells at the follicle anterior. Polarized reorientation is regulated by mechanical cues from the basement membrane, which are transduced by the Src tyrosine kinase to alter junctional E-cadherin trafficking. This mechanosensitive cellular behavior represents a conserved mechanism that can elongate edgeless tubular epithelia in a process distinct from those that elongate bounded, planar epithelia.

EGFR signaling coordinates patterning with cell survival during Drosophila epidermal development.

Extensive apoptosis is often seen in patterning mutants, suggesting that tissues can detect and eliminate potentially harmful mis-specified cells. Here, we show that the pattern of apoptosis in the embryonic epidermis of Drosophila is not a response to fate mis-specification but can instead be explained by the limiting availability of prosurvival signaling molecules released from locations determined by patterning information. In wild-type embryos, the segmentation cascade elicits the segmental production of several epidermal growth factor receptor (EGFR) ligands, including the transforming growth factor Spitz (TGFα), and the neuregulin, Vein. This leads to an undulating pattern of signaling activity, which prevents expression of the proapoptotic gene head involution defective (hid) throughout the epidermis. In segmentation mutants, where specific peaks of EGFR ligands fail to form, gaps in signaling activity appear, leading to coincident hid up-regulation and subsequent cell death. These data provide a mechanistic understanding of how cell survival, and thus appropriate tissue size, is made contingent on correct patterning.

Global morphogenetic flow is accurately predicted by the spatial distribution of myosin motors.

During embryogenesis tissue layers undergo morphogenetic flow rearranging and folding into specific shapes. While developmental biology has identified key genes and local cellular processes, global coordination of tissue remodeling at the organ scale remains unclear. Here, we combine in toto light-sheet microscopy of the Drosophila embryo with quantitative analysis and physical modeling to relate cellular flow with the patterns of force generation during the gastrulation process. We find that the complex spatio-temporal flow pattern can be predicted from the measured meso-scale myosin density and anisotropy using a simple, effective viscous model of the tissue, achieving close to 90% accuracy with one time dependent and two constant parameters. Our analysis uncovers the importance of a) spatial modulation of myosin distribution on the scale of the embryo and b) the non-locality of its effect due to mechanical interaction of cells, demonstrating the need for the global perspective in the study of morphogenetic flow.

Identification of a neural crest stem cell niche by Spatial Genomic Analysis.

The neural crest is an embryonic population of multipotent stem cells that form numerous defining features of vertebrates. Due to lack of reliable techniques to perform transcriptional profiling in intact tissues, it remains controversial whether the neural crest is a heterogeneous or homogeneous population. By coupling multiplex single molecule fluorescence in situ hybridization with machine learning algorithm based cell segmentation, we examine expression of 35 genes at single cell  resolution in vivo. Unbiased hierarchical clustering reveals five spatially distinct subpopulations within the chick dorsal neural tube. Here we identify a neural crest stem cell niche that centers around the dorsal midline with high expression of neural crest genes, pluripotency factors, and lineage markers. Interestingly, neural and neural crest stem cells express distinct pluripotency signatures. This Spatial Genomic Analysis toolkit provides a straightforward approach to study quantitative multiplex gene expression in numerous biological systems, while offering insights into gene regulatory networks via synexpression analysis.

Active Tension Network model suggests an exotic mechanical state realized in epithelial tissues.

Mechanical interactions play a crucial role in epithelial morphogenesis, yet understanding the complex mechanisms through which stress and deformation affect cell behavior remains an open problem. Here we formulate and analyze the Active Tension Network (ATN) model, which assumes that the mechanical balance of cells within a tissue is dominated by cortical tension and introduces tension-dependent active remodeling of the cortex. We find that ATNs exhibit unusual mechanical properties. Specifically, an ATN behaves as a fluid at short times, but at long times supports external tension like a solid. Furthermore, an ATN has an extensively degenerate equilibrium mechanical state associated with a discrete conformal - "isogonal" - deformation of cells. The ATN model predicts a constraint on equilibrium cell geometries, which we demonstrate to approximately hold in certain epithelial tissues. We further show that isogonal modes are observed in the fruit y embryo, accounting for the striking variability of apical areas of ventral cells and helping understand the early phase of gastrulation. Living matter realizes new and exotic mechanical states, the study of which helps to understand biological phenomena.

Symmetry Breaking in an Edgeless Epithelium by Fat2-Regulated Microtubule Polarity.

Planar cell polarity (PCP) information is a critical determinant of organ morphogenesis. While PCP in bounded epithelial sheets is increasingly well understood, how PCP is organized in tubular and acinar tissues is not. Drosophila egg chambers (follicles) are an acinus-like "edgeless epithelium" and exhibit a continuous, circumferential PCP that does not depend on pathways active in bounded epithelia; this follicle PCP directs formation of an ellipsoid rather than a spherical egg. Here, we apply an imaging algorithm to "unroll" the entire 3D tissue surface and comprehensively analyze PCP onset. This approach traces chiral symmetry breaking to plus-end polarity of microtubules in the germarium, well before follicles form and rotate. PCP germarial microtubules provide chiral information that predicts the direction of whole-tissue rotation as soon as independent follicles form. Concordant microtubule polarity, but not microtubule alignment, requires the atypical cadherin Fat2, which acts at an early stage to translate plus-end bias into coordinated actin-mediated collective cell migration. Because microtubules are not required for PCP or migration after follicle rotation initiates, while dynamic actin and extracellular matrix are, polarized microtubules lie at the beginning of a handoff mechanism that passes early chiral PCP of the cytoskeleton to a supracellular planar polarized extracellular matrix and elongates the organ.