Integrating morphogenesis with underlying mechanics and cell biology

Multiscale analysis of architecture, cell size and the cell cortex reveals cortical F-actin density and composition are major contributors to mechanical properties during convergent extension

The large-scale movements that construct complex three-dimensional tissues during development are governed by universal physical principles. Fine-grained control of both mechanical properties and force production is crucial to the successful placement of tissues and shaping of organs. Embryos of the frog Xenopus laevis provide a dramatic example of these physical processes, as dorsal tissues increase in Young's modulus by six-fold to 80 Pascal over 8 h as germ layers and the central nervous system are formed. These physical changes coincide with emergence of complex anatomical structures, rounds of cell division, and cytoskeletal remodeling. To understand the contribution of these diverse structures, we adopt the cellular solids model to relate bulk stiffness of a solid foam to the unit size of individual cells, their microstructural organization, and their material properties. Our results indicate that large-scale tissue architecture and cell size are not likely to influence the bulk mechanical properties of early embryonic or progenitor tissues but that F-actin cortical density and composition of the F-actin cortex play major roles in regulating the physical mechanics of embryonic multicellular tissues.

Large, long range tensile forces drive convergence during Xenopus blastopore closure and body axis elongation

Indirect evidence suggests that blastopore closure during gastrulation of anamniotes, including amphibians such as Xenopus laevis, depends on circumblastoporal convergence forces generated by the marginal zone (MZ), but direct evidence is lacking. We show that explanted MZs generate tensile convergence forces up to 1.5 μN during gastrulation and over 4 μN thereafter. These forces are generated by convergent thickening (CT) until the midgastrula and increasingly by convergent extension (CE) thereafter. Explants from ventralized embryos, which lack tissues expressing CE but close their blastopores, produce up to 2 μN of tensile force, showing that CT alone generates forces sufficient to close the blastopore. Uniaxial tensile stress relaxation assays show stiffening of mesodermal and ectodermal tissues around the onset of neurulation, potentially enhancing long-range transmission of convergence forces. These results illuminate the mechanobiology of early vertebrate morphogenic mechanisms, aid interpretation of phenotypes, and give insight into the evolution of blastopore closure mechanisms.

Tissue-specific roles of Fgfr2 in development of the external genitalia

Congenital anomalies frequently occur in organs that undergo tubulogenesis. Hypospadias is a urethral tube defect defined by mislocalized, oversized, or multiple openings of the penile urethra. Deletion of Fgfr2 or its ligand Fgf10 results in severe hypospadias in mice, in which the entire urethral plate is open along the ventral side of the penis. In the genital tubercle, the embryonic precursor of the penis and clitoris, Fgfr2 is expressed in two epithelial populations: the endodermally derived urethral epithelium and the ectodermally derived surface epithelium. Here, we investigate the tissue-specific roles of Fgfr2 in external genital development by generating conditional deletions of Fgfr2 in each of these cell types. Conditional deletion of Fgfr2 results in two distinct phenotypes: endodermal Fgfr2 deletion causes mild hypospadias and inhibits maturation of a complex urethral epithelium, whereas loss of ectodermal Fgfr2 results in severe hypospadias and absence of the ventral prepuce. Although these cell type-specific mutants exhibit distinctive genital anomalies, cellular analysis reveals that Fgfr2 regulates epithelial maturation and cell cycle progression in the urethral endoderm and in the surface ectoderm. The unexpected finding that ectodermal deletion of Fgfr2 results in the most severe hypospadias highlights a major role for Fgfr2 in the developing genital surface epithelium, where epithelial maturation is required for maintenance of a closed urethral tube. These results demonstrate that urethral tubulogenesis, prepuce morphogenesis, and sexually dimorphic patterning of the lower urethra are controlled by discrete regions of Fgfr2 activity.

CellFIT: a cellular force-inference toolkit using curvilinear cell boundaries

Mechanical forces play a key role in a wide range of biological processes, from embryogenesis to cancer metastasis, and there is considerable interest in the intuitive question, "Can cellular forces be inferred from cell shapes?" Although several groups have posited affirmative answers to this stimulating question, nagging issues remained regarding equation structure, solution uniqueness and noise sensitivity. Here we show that the mechanical and mathematical factors behind these issues can be resolved by using curved cell edges rather than straight ones. We present a new package of force-inference equations and assessment tools and denote this new package CellFIT, the Cellular Force Inference Toolkit. In this approach, cells in an image are segmented and equilibrium equations are constructed for each triple junction based solely on edge tensions and the limiting angles at which edges approach each junction. The resulting system of tension equations is generally overdetermined. As a result, solutions can be obtained even when a modest number of edges need to be removed from the analysis due to short length, poor definition, image clarity or other factors. Solving these equations yields a set of relative edge tensions whose scaling must be determined from data external to the image. In cases where intracellular pressures are also of interest, Laplace equations are constructed to relate the edge tensions, curvatures and cellular pressure differences. That system is also generally overdetermined and its solution yields a set of pressures whose offset requires reference to the surrounding medium, an open wound, or information external to the image. We show that condition numbers, residual analyses and standard errors can provide confidence information about the inferred forces and pressures. Application of CellFIT to several live and fixed biological tissues reveals considerable force variability within a cell population, significant differences between populations and elevated tensions along heterotypic boundaries.

Combined hydrogels that switch human pluripotent stem cells from self-renewal to differentiation

The ability of materials to define the architecture and microenvironment experienced by cells provides new opportunities to direct the fate of human pluripotent stem cells (HPSCs) [Robinton DA, Daley GQ (2012) Nature 481(7381):295-305]. However, the conditions required for self-renewal vs. differentiation of HPSCs are different, and a single system that efficiently achieves both outcomes is not available [Giobbe GG, et al. (2012) Biotechnol Bioeng 109(12):3119-3132]. We have addressed this dual need by developing a hydrogel-based material that uses ionic de-cross-linking to remove a self-renewal permissive hydrogel (alginate) and switch to a differentiation-permissive microenvironment (collagen). Adjusting the timing of this switch can preferentially steer the HPSC differentiation to mimic lineage commitment during gastrulation to ectoderm (early switch) or mesoderm/endoderm (late switch). As an exemplar differentiated cell type, we showed that directing early lineage specification using this single system can promote cardiogenesis with increased gene expression in high-density cell populations. This work will facilitate regenerative medicine by allowing in situ HPSC expansion to be coupled with early lineage specification within defined tissue geometries.

Calcium wave propagation in networks of endothelial cells: model-based theoretical and experimental study

In this paper, we present a combined theoretical and experimental study of the propagation of calcium signals in multicellular structures composed of human endothelial cells. We consider multicellular structures composed of a single chain of cells as well as a chain of cells with a side branch, namely a “T” structure. In the experiments, we investigate the result of applying mechano-stimulation to induce signaling in the form of calcium waves along the chain and the effect of single and dual stimulation of the multicellular structure. The experimental results provide evidence of an effect of architecture on the propagation of calcium waves. Simulations based on a model of calcium-induced calcium release and cell-to-cell diffusion through gap junctions shows that the propagation of calcium waves is dependent upon the competition between intracellular calcium regulation and architecture-dependent intercellular diffusion.

Calcium wave signal has been found in a wide variety of cell types. Over the last years, a large number of calcium experiments have shown that calcium signal is not only an intracellular regulator but is also able to be transmitted to surrounding cells as intercellular signal. This paper focuses on the development of an approach with complementary integration of theoretical and experimental methods for studying the multi-level interactions in multicellular architectures and their effect on collective cell dynamic behavior. We describe new types of higher-order (across structure) behaviors arising from lower-order (within cells) phenomena, and make predictions concerning the mechanisms underlying the dynamics of multicellular biological systems. The theoretical approach describes numerically the dynamics of non-linear behavior of calcium-based signaling in model networks of cells. Microengineered, geometrically constrained networks of human umbilical vein endothelial cells (HUVEC) serve as platforms to arbitrate the theoretical predictions in terms of the effect of network topology on the spatiotemporal characteristics of emerging calcium signals.

Computational models for mechanics of morphogenesis

In the developing embryo, tissues differentiate, deform, and move in an orchestrated manner to generate various biological shapes driven by the complex interplay between genetic, epigenetic, and environmental factors. Mechanics plays a key role in regulating and controlling morphogenesis, and quantitative models help us understand how various mechanical forces combine to shape the embryo. Models allow for the quantitative, unbiased testing of physical mechanisms, and when used appropriately, can motivate new experimentaldirections. This knowledge benefits biomedical researchers who aim to prevent and treat congenital malformations, as well as engineers working to create replacement tissues in the laboratory. In this review, we first give an overview of fundamental mechanical theories for morphogenesis, and then focus on models for specific processes, including pattern formation, gastrulation, neurulation, organogenesis, and wound healing. The role of mechanical feedback in development is also discussed. Finally, some perspectives aregiven on the emerging challenges in morphomechanics and mechanobiology.

The mechanics of metastasis: insights from a computational model

Although it may seem obvious that mechanical forces are required to drive metastatic cell movements, understanding of the mechanical aspects of metastasis has lagged far behind genetic and biochemical knowledge. The goal of this study is to learn about the mechanics of metastasis using a cell-based finite element model that proved useful for advancing knowledge about the forces that drive embryonic cell and tissue movements. Metastasis, the predominant cause of cancer-related deaths, involves a series of mechanical events in which one or more cells dissociate from a primary tumour, migrate through normal tissue, traverse in and out of a multi-layer circulatory system vessel and resettle. The present work focuses on the dissemination steps, from dissociation to circulation. The model shows that certain surface tension relationships must be satisfied for cancerous cells to dissociate from a primary tumour and that these equations are analogous to those that govern dissociation of embryonic cells. For a dissociated cell to then migrate by invadopodium extension and contraction and exhibit the shapes seen in experiments, the invadopodium must generate a contraction equal to approximately twice that produced by the interfacial tension associated with surrounding cells. Intravasation through the wall of a vessel is governed by relationships akin to those in the previous two steps, while release from the vessel wall is governed by equations that involve surface and interfacial tensions. The model raises a number of potential research questions. It also identifies how specific mechanical properties and the sub-cellular structural components that give rise to them might be changed so as to thwart particular metastatic steps and thereby block the spread of cancer.

Independent migration of cell populations in the early gastrulation of the amphipod crustacean Parhyale hawaiensis

Cells are the principal component of tissues and can drive morphogenesis through dynamic changes in structure and interaction. During gastrulation, the primary morphogenetic event of early development, cells change shape, exchange neighbors, and migrate long distances to establish cell layers that will form the tissues of the adult animal. Outside of Drosophila, little is known about how changes in cell behavior might drive gastrulation among arthropods. Here, we focus on three cell populations that form two aggregations during early gastrulation in the crustacean Parhyale hawaiensis. Using cytoskeletal markers and lineage tracing we observe bottle cells in anterior and visceral mesoderm precursors as gastrulation commences, and find that both Cytochalasin D, an inhibitor of actin polymerization, and ROCKOUT, an inhibitor of Rho-kinase activity, prevent gastrulation. Furthermore, by ablating specific cells, we show that each of the three populations acts independently during gastrulation, confirming previous hypotheses that cell behavior during Parhyale gastrulation relies on intrinsic signals instead of an inductive mechanism.

Mechanics of head fold formation: investigating tissue-level forces during early development

During its earliest stages, the avian embryo is approximately planar. Through a complex series of folds, this flat geometry is transformed into the intricate three-dimensional structure of the developing organism. Formation of the head fold (HF) is the first step in this cascading sequence of out-of-plane tissue folds. The HF establishes the anterior extent of the embryo and initiates heart, foregut and brain development. Here, we use a combination of computational modeling and experiments to determine the physical forces that drive HF formation. Using chick embryos cultured ex ovo, we measured: (1) changes in tissue morphology in living embryos using optical coherence tomography (OCT); (2) morphogenetic strains (deformations) through the tracking of tissue labels; and (3) regional tissue stresses using changes in the geometry of circular wounds punched through the blastoderm. To determine the physical mechanisms that generate the HF, we created a three-dimensional computational model of the early embryo, consisting of pseudoelastic plates representing the blastoderm and vitelline membrane. Based on previous experimental findings, we simulated the following morphogenetic mechanisms: (1) convergent extension in the neural plate (NP); (2) cell wedging along the anterior NP border; and (3) autonomous in-plane deformations outside the NP. Our numerical predictions agree relatively well with the observed morphology, as well as with our measured stress and strain distributions. The model also predicts the abnormal tissue geometries produced when development is mechanically perturbed. Taken together, the results suggest that the proposed morphogenetic mechanisms provide the main tissue-level forces that drive HF formation.

Calcium wave propagation in chains of endothelial cells with nonlinear reaction dynamics: Green's function approach

We present a Green's function-based perturbative approach to solving nonlinear reaction-diffusion problems in networks of endothelial cells. We focus on a single component (Ca2+), piecewise nonlinear model of endoplasmic calcium dynamics and trans-membrane diffusion. The decoupling between nonlinear reaction dynamics and the linear diffusion enables the calculation of the diffusion part of the Green's function for network of cells with nontrivial topologies. We verify analytically and then numerically that our approach leads to the known transition from propagation of calcium front to failure of propagation when the diffusion rate is varied relative to the reaction rates. We then derive the Green's function for a semi-infinite chain of cells with various boundary conditions. We show that the calcium dynamics of cells in the vicinity of the end of the semi-infinite chain is strongly dependent on the boundary conditions. The behavior of the semi-infinite chain with absorbing boundary conditions, a simple model of a multicellular structure with an end in contact with the extracellular matrix, suggests behavioral differentiation between cells at the end and cells embedded within the chain.

Biological evolution based on nonrandom variability regulated by the organism

A hypothetical mechanism for rapid and nonrandom emergence of evolutionary adaptations is proposed. It is supposed that some transcription factors and transcription regulators that are able to cross membranes can leave the cells of their origin and move within the organism using a specialized transport system when individual development occurs under conditions extreme for the given species. This system, in particular, connects soma with germline. The supply of germline cells with unusual transcription regulators changes the balance of their nuclear regulatory RNAs, thus initiating RNA-dependent epigenetic modifications of the germline genome and therefore changes in phenotypes of the progeny. It is highly probable that some of these phenotypes are adaptive and lay the basis for the origin of the next biological species. The proposed mechanism can serve as a basis for a new theory of the origin of species.

Mechanical determinants of epithelium thickness in early-stage embryos

To address the fundamental question, "How is the thickness of an embryonic epithelium determined?" equations of force equilibrium are derived and used to construct a theory applicable to embryonic epithelia and other labile monolayers, including confluent cultured cells. Under typical physiological conditions, the underlying mathematical relationships can be reworked into one simple algebraic equation. The theory shows that sheet thickness is determined by mechanical interactions between microfilaments, microtubules, cell adhesion systems and in-plane loading. It also explains the perplexing negative in-plane Poisson's ratio observed in some uniaxial tensile tests of embryonic epithelia. Thickness anomalies can affect phenotype, and the equations make it possible to calculate the degree to which gene expression, signalling pathways or other factors must alter the forces at work in order to produce a specified change in sheet thicknesses. Key findings are confirmed using a cell-based finite element model.