Interstitial fluid osmolarity modulates the action of differential tissue surface tension in progenitor cell segregation during gastrulation.

The segregation of different cell types into distinct tissues is a fundamental process in metazoan development. Differences in cell adhesion and cortex tension are commonly thought to drive cell sorting by regulating tissue surface tension (TST). However, the role that differential TST plays in cell segregation within the developing embryo is as yet unclear. Here, we have analyzed the role of differential TST for germ layer progenitor cell segregation during zebrafish gastrulation. Contrary to previous observations that differential TST drives germ layer progenitor cell segregation in vitro, we show that germ layers display indistinguishable TST within the gastrulating embryo, arguing against differential TST driving germ layer progenitor cell segregation in vivo We further show that the osmolarity of the interstitial fluid (IF) is an important factor that influences germ layer TST in vivo, and that lower osmolarity of the IF compared with standard cell culture medium can explain why germ layers display differential TST in culture but not in vivo Finally, we show that directed migration of mesendoderm progenitors is required for germ layer progenitor cell segregation and germ layer formation.

Elongated Cells Drive Morphogenesis in a Surface-Wrapped Finite-Element Model of Germband Retraction.

During Drosophila embryogenesis, the germband first extends to curl around the posterior end of the embryo and then retracts back; however, retraction is not simply the reversal of extension. At a tissue level, extension is coincident with ventral furrow formation, and at a cellular level, extension occurs via convergent cell neighbor exchanges in the germband, whereas retraction involves only changes in cell shape. To understand how cell shapes, tissue organization, and cellular forces drive germband retraction, we investigate this process using a whole-embryo, surface-wrapped cellular finite-element model. This model represents two key epithelial tissues-amnioserosa and germband-as adjacent sheets of two-dimensional cellular finite elements that are wrapped around an ellipsoidal three-dimensional approximation of an embryo. The model reproduces the detailed kinematics of in vivo retraction by fitting just one free model parameter, the tension along germband cell interfaces; all other cellular forces are constrained to follow ratios inferred from experimental observations. With no additional parameter adjustments, the model also reproduces quantitative assessments of mechanical stress using laser dissection and failures of retraction when amnioserosa cells are removed via mutations or microsurgery. Surprisingly, retraction in the model is robust to changes in cellular force values but is critically dependent on starting from a configuration with highly elongated amnioserosa cells. Their extreme cellular elongation is established during the prior process of germband extension and is then used to drive retraction. The amnioserosa is the one tissue whose cellular morphogenesis is reversed from germband extension to retraction, and this reversal coordinates the forces needed to retract the germband back to its pre-extension position and shape. In this case, cellular force strengths are less important than the carefully established cell shapes that direct them. VIDEO ABSTRACT.

Modeling cell elongation during germ band retraction: cell autonomy versus applied anisotropic stress.

The morphogenetic process of germ band retraction in Drosophila embryos involves coordinated movements of two epithelial tissues - germ band and amnioserosa. The germ band shortens along its rostral-caudal or head-to-tail axis, widens along its perpendicular dorsal-ventral axis, and uncurls from an initial 'U' shape. The amnioserosa mechanically assists this process by pulling on the crook of the U-shaped germ band. The amnioserosa may also provide biochemical signals that drive germ band cells to change shape in a mechanically autonomous fashion. Here, we use a finite-element model to investigate how these two contributions reshape the germ band. We do so by modeling the response to laser-induced wounds in each of the germ band's spatially distinct segments (T1-T3, A1-A9) during the middle of retraction when segments T1-A3 form the ventral arm of the 'U', A4-A7 form its crook, and A8-A9 complete the dorsal arm. We explore these responses under a range of externally applied stresses and internal anisotropy of cell edge tensions - akin to a planar cell polarity that can drive elongation of cells in a direction parallel to the minimum edge tension - and identify regions of parameter space (edge-tension anisotropy versus stress anisotropy) that best match previous experiments for each germ band segment. All but three germ band segments are best fit when the applied stress anisotropy and the edge-tension anisotropy work against one another - i.e., when the isolated effects would elongate cells in perpendicular directions. Segments in the crook of the germ band (A4-A7) have cells that elongate in the direction of maximum external stress, i.e., external stress anisotropy is dominant. In most other segments, the dominant factor is internal edge-tension anisotropy. These results are consistent with models in which the amnioserosa pulls on the crook of the germ band to mechanically assist retraction. In addition, they suggest a mechanical cue for edge-tension anisotropy whereby cells do not globally orient their internal elongation axis towards the amnioserosa, but instead orient this axis perpendicular to the local principal stress direction.

Video force microscopy reveals the mechanics of ventral furrow invagination in Drosophila.

The absence of tools for mapping the forces that drive morphogenetic movements in embryos has impeded our understanding of animal development. Here we describe a unique approach, video force microscopy (VFM), that allows detailed, dynamic force maps to be produced from time-lapse images. The forces at work in an embryo are considered to be decomposed into active and passive elements, where active forces originate from contributions (e.g., actomyosin contraction) that do mechanical work to the system and passive ones (e.g., viscous cytoplasm) that dissipate energy. In the present analysis, the effects of all passive components are considered to be subsumed by an effective cytoplasmic viscosity, and the driving forces are resolved into equivalent forces along the edges of the polygonal boundaries into which the region of interest is divided. Advanced mathematical inverse methods are used to determine these driving forces. When applied to multiphoton sections of wild-type and mutant Drosophila melanogaster embryos, VFM is able to calculate the equivalent driving forces acting along individual cell edges and to do so with subminute temporal resolution. In the wild type, forces along the apical surface of the presumptive mesoderm are found to be large and to vary parabolically with time and angular position, whereas forces along the basal surface of the ectoderm, for example, are found to be smaller and nearly uniform with position. VFM shows that in mutants with reduced junction integrity and myosin II activity, the driving forces are reduced, thus accounting for ventral furrow failure.

Digital tracking algorithm reveals the influence of structural irregularities on joint movements in the human cervical spine.

BACKGROUND: Disc height loss and osteophytes change the local mechanical environment in the spine; while previous research has examined kinematic dysfunction under degenerative change, none has looked at the influence of disc height loss and osteophytes throughout movement. METHODS: Twenty patients with pain related to the head, neck or shoulders were imaged via videofluoroscopy as they underwent sagittal-plane flexion and extension. A clinician graded disc height loss and osteophytes as "severe/moderate", "mild", or "none". A novel tracking algorithm quantified motions of each vertebra. This information was used to calculate intervertebral angular and shear displacements. The digital algorithm made it practical to track individual vertebrae in multiple patients through hundreds of images without bias. FINDINGS: Cases without height loss/osteophytes had a consistent increase in intervertebral angular displacement from C2/C3 to C5/C6, like that of healthy individuals, and mild height losses did not produce aberrations that were systematic or necessarily discernable. However, joints with moderate to severe disc height loss and osteophytes exhibited reduced range of motion compared to adjacent unaffected joints in that patient and corresponding joints in patients without structural irregularities. INTERPRETATION: Digitally-obtained motion histories of individual joints allowed anatomical joint changes to be linked with changes in joint movement patterns. Specifically, disc height loss and osteophytes were found to influence cervical spine movement in the sagittal plane, reducing angular motions at affected joints by approximately 10% between those with and without height loss and osteophytes. Further, these joint changes were associated with perturbed intervertebral angular and shear movements.

Coordination of Receptor Tyrosine Kinase Signaling and Interfacial Tension Dynamics Drives Radial Intercalation and Tube Elongation.

We sought to understand how cells collectively elongate epithelial tubes. We first used 3D culture and biosensor imaging to demonstrate that epithelial cells enrich Ras activity, phosphatidylinositol (3,4,5)-trisphosphate (PIP3), and F-actin to their leading edges during migration within tissues. PIP3 enrichment coincided with, and could enrich despite inhibition of, F-actin dynamics, revealing a conserved migratory logic compared with single cells. We discovered that migratory cells can intercalate into the basal tissue surface and contribute to tube elongation. We then connected molecular activities to subcellular mechanics using force inference analysis. Migration and transient intercalation required specific and similar anterior-posterior ratios of interfacial tension. Permanent intercalations were distinguished by their capture at the boundary through time-varying tension dynamics. Finally, we integrated our experimental and computational data to generate a finite element model of tube elongation. Our model revealed that intercalation, interfacial tension dynamics, and high basal stress are together sufficient for mammary morphogenesis.

A new cell-based FE model for the mechanics of embryonic epithelia.

In order to overcome a significant stiffening artefact associated with current finite element (FE) models for the mechanics of embryonic epithelia, two new FE formulations were developed. Cell-cell interfacial tensions gamma are represented by constant-force rod elements as in previous models. However, the viscosity of the cytoplasm with its embedded organelles and filament networks is modeled using viscous triangular elements, it is modeled using either radial and circumferential dashpots or an orthogonal dashpot system rather than the viscous triangular elements typical of previous two-dimensional FE models. The models are tested against tissue (epithelium) stretching because it gives rise to significant changes in cell shape and against cell sorting because it involves high rates of cell rearrangement. The orthogonal dashpot system is found to capture cell size and shape effects well, give the model cells characteristics that are consistent with those of real cells, provide high computational efficiency and hold promise for future three-dimensional analyses.

Inferring cellular forces from image stacks.

Although the importance of cellular forces to a wide range of embryogenesis and disease processes is widely recognized, measuring these forces is challenging, especially in three dimensions. Here, we introduce CellFIT-3D, a force inference technique that allows tension maps for three-dimensional cellular systems to be estimated from image stacks. Like its predecessors, video force microscopy and CellFIT, this cell mechanics technique assumes boundary-specific interfacial tensions to be the primary drivers, and it constructs force-balance equations based on triple junction (TJ) dihedral angles. The technique involves image processing, segmenting of cells, grouping of cell outlines, calculation of dihedral planes, averaging along three-dimensional TJs, and matrix equation assembly and solution. The equations tend to be strongly overdetermined, allowing indistinct TJs to be ignored and solution error estimates to be determined. Application to clean and noisy synthetic data generated using Surface Evolver gave tension errors of 1.6-7%, and analyses of eight-cell murine embryos gave estimated errors smaller than the 10% uncertainty of companion aspiration experiments. Other possible areas of application include morphogenesis, cancer metastasis and tissue engineering.This article is part of the themed issue 'Systems morphodynamics: understanding the development of tissue hardware'.

A videofluoroscopy-based tracking algorithm for quantifying the time course of human intervertebral displacements.

The motions of individual intervertebral joints can affect spine motion, injury risk, deterioration, pain, treatment strategies, and clinical outcomes. Since standard kinematic methods do not provide precise time-course details about individual vertebrae and intervertebral motions, information that could be useful for scientific advancement and clinical assessment, we developed an iterative template matching algorithm to obtain this data from videofluoroscopy images. To assess the bias of our approach, vertebrae in an intact porcine spine were tracked and compared to the motions of high-contrast markers. To estimate precision under clinical conditions, motions of three human cervical spines were tracked independently ten times and vertebral and intervertebral motions associated with individual trials were compared to corresponding averages. Both tests produced errors in intervertebral angular and shear displacements no greater than 0.4° and 0.055 mm, respectively. When applied to two patient cases, aberrant intervertebral motions in the cervical spine were typically found to correlate with patient-specific anatomical features such as disc height loss and osteophytes. The case studies suggest that intervertebral kinematic time-course data could have value in clinical assessments, lead to broader understanding of how specific anatomical features influence joint motions, and in due course inform clinical treatments.

Non-straight cell edges are important to invasion and engulfment as demonstrated by cell mechanics model.

Computational models of cell-cell mechanical interactions typically simulate sorting and certain other motions well, but as demands on these models continue to grow, discrepancies between the cell shapes, contact angles and behaviours they predict and those that occur in real cells have come under increased scrutiny. To investigate whether these discrepancies are a direct result of the straight cell-cell edges generally assumed in these models, we developed a finite element model that approximates cell boundaries using polylines with an arbitrary number of segments. We then compared the predictions of otherwise identical polyline and monoline (straight-edge) models in a variety of scenarios, including annealing, single- and multi-cell engulfment, sorting, and two forms of mixing--invasion and checkerboard pattern formation. Keeping cell-cell edges straight influences cell motion, cell shape, contact angle, and boundary length, especially in cases where one cell type is pulled between or around cells of a different type, as in engulfment or invasion. These differences arise because monoline cells have restricted deformation modes. Polyline cells do not face these restrictions, and with as few as three segments per edge yielded realistic edge shapes and contact angle errors one-tenth of those produced by monoline models, making them considerably more suitable for situations where angles and shapes matter, such as validation of cellular force-inference techniques. The findings suggest that non-straight cell edges are important both in modelling and in nature.

Practical aspects of the cellular force inference toolkit (CellFIT).

If we are to fully understand the reasons that cells and tissues move and acquire their distinctive geometries during processes such as embryogenesis and wound healing, we will need detailed maps of the forces involved. One of the best current prospects for obtaining this information is noninvasive force-from-images techniques such as CellFIT, the Cellular Force Inference Toolkit, whose various steps are discussed here. Like other current quasistatic approaches, this one assumes that cell shapes are produced by interactions between interfacial tensions and intracellular pressures. CellFIT, however, allows cells to have curvilinear boundaries, which can significantly improve inference accuracy and reduce noise sensitivity. The quality of a CellFIT analysis depends on how accurately the junction angles and edge curvatures are measured, and a software tool we describe facilitates determination and evaluation of this information. Special attention is required when edges are crenulated or significantly different in shape from a circular arc. Because the tension and pressure equations are overdetermined, a select number of edges can be removed from the analysis, and these might include edges that are poorly defined in the source image, too short to provide accurate angles or curvatures, or noncircular. The approach works well for aggregates with as many as 1000 cells, and introduced errors have significant effects on only a few adjacent cells. An understanding of these considerations will help CellFIT users to get the most out of this promising new technique.

Computer simulations of mitosis and interdependencies between mitosis orientation, cell shape and epithelia reshaping.

Finite element-based computer simulations are used to investigate mitosis and how mitosis, cell shape, and epithelium reshaping depend on each other. Frame- and cell-oriented patterns of mitosis with growing and non-growing daughter cells are considered. Previous simulations have shown that applied stresses or strains can reshape cells so that their long axes are aligned in the principal stretch direction. The simulations reported here show that this can produce global alignment of the mitosis cleavage planes. Other simulations reported here show that mitoses with suitably aligned cleavage planes can drive epithelium reshaping. Formulas that quantify these and other dependencies are derived. These formulas provide quantitative relationships against which current hypotheses regarding epithelia reshaping in real biological systems can be evaluated.

A computer model for reshaping of cells in epithelia due to in-plane deformation and annealing.

Although cell reshaping is fundamental to the mechanics of epithelia, technical barriers have prevented the methods of mechanics from being used to investigate it. These barriers have recently been overcome by the cell-based finite element formulation of Chen and Brodland. Here, parameters to describe the fabric of an epithelium in terms of cell shape and orientation and cell edge density are defined. Then, rectangular "patches" of model epithelia having various initial fabric parameters are generated and are either allowed to anneal or are subjected to one of several patterns of in-plane deformation. The simulations show that cell reshaping lags the deformation history, that it is allayed by cell rearrangement and that it causes the epithelium as a whole to exhibit viscoelastic mechanical properties. Equations to describe changes in cell shape due to annealing and in-plane deformation are presented.

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.

Forces driving epithelial wound healing.

A fundamental feature of multicellular organisms is their ability to self-repair wounds through the movement of epithelial cells into the damaged area. This collective cellular movement is commonly attributed to a combination of cell crawling and "purse-string" contraction of a supracellular actomyosin ring. Here we show by direct experimental measurement that these two mechanisms are insufficient to explain force patterns observed during wound closure. At early stages of the process, leading actin protrusions generate traction forces that point away from the wound, showing that wound closure is initially driven by cell crawling. At later stages, we observed unanticipated patterns of traction forces pointing towards the wound. Such patterns have strong force components that are both radial and tangential to the wound. We show that these force components arise from tensions transmitted by a heterogeneous actomyosin ring to the underlying substrate through focal adhesions. The structural and mechanical organization reported here provides cells with a mechanism to close the wound by cooperatively compressing the underlying substrate.

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.

Assessing the mechanical energy costs of various tissue reshaping mechanisms.

Early-stage embryos must reshape the tissues of which they are made into organs and other life-sustaining structures; and if non-mammalian embryos fail to complete these tasks before the energy provided by their yolk runs out, they die. The aim of this study is to use a cell-level computational model to investigate the energetic cost of a variety of mechanisms that can drive an in-plane reshaping pattern known as convergent extension--a motif in which a tissue narrows in one in-plane direction and expands in another. Mechanisms considered include oriented lamellipodia, directed mitosis, stress fibers, and anisotropic external tension. Both isolated patches of tissue and actively contracting tissues that deform adjacent passive areas are considered. The cell-level finite element model used here assumes that the cell membrane and its associated proteins generate a net tension γ along each cell-cell interface and that the cytoplasm and its embedded networks and structures have an effective viscosity µ. Work costs are based exclusively on mechanical considerations such as edge lengths and tensions, and because a traditional mechanical efficiency cannot be calculated, mechanisms are compared on the basis of the work they must do to the tissue to cause a specified rate of in-plane reshaping. Although the model contains a number of simplifications compared to real embryonic tissues, it is able to show that the work requirements for tissue reshaping by mitoses and by lamellipodia are of the same order. Lamellipodia are energetically most effective when their tensions are approximately twice as large as the interfacial tensions in the surrounding cells. The model also shows that stress fibers or other direct stretch or compression mechanisms are at least five times more efficient for tissue reshaping than are mitoses or lamellipodia and that the work needed to deform a typical cellular tissue is more than thirty times greater than if it did not contain cell boundaries. Collectively, these findings indicate that common tissue reshaping mechanisms have mechanical efficiencies of less than one percent and that mechanical efficiency is not the primary determinant of which mechanism(s) an embryo uses to reshape its tissues.

A biomechanical analysis of ventral furrow formation in the Drosophila melanogaster embryo.

The article provides a biomechanical analysis of ventral furrow formation in the Drosophila melanogaster embryo. Ventral furrow formation is the first large-scale morphogenetic movement in the fly embryo. It involves deformation of a uniform cellular monolayer formed following cellularisation, and has therefore long been used as a simple system in which to explore the role of mechanics in force generation. Here we use a quantitative framework to carry out a systematic perturbation analysis to determine the role of each of the active forces observed. The analysis confirms that ventral furrow invagination arises from a combination of apical constriction and apical-basal shortening forces in the mesoderm, together with a combination of ectodermal forces. We show that the mesodermal forces are crucial for invagination: the loss of apical constriction leads to a loss of the furrow, while the mesodermal radial shortening forces are the primary cause of the internalisation of the future mesoderm as the furrow rises. Ectodermal forces play a minor but significant role in furrow formation: without ectodermal forces the furrow is slower to form, does not close properly and has an aberrant morphology. Nevertheless, despite changes in the active mesodermal and ectodermal forces lead to changes in the timing and extent of furrow, invagination is eventually achieved in most cases, implying that the system is robust to perturbation and therefore over-determined.

Novel lap test determines the mechanics of delamination between annular lamellae of the intervertebral disc.

Delamination between lamellae of the annulus fibrosus is a crucial stage of intervertebral disc herniation, and to better understand the mechanics of the delamination process, a novel lap test was devised. Specimens consisting of two adjacent, naturally bonded lamellae were obtained from the cervical region of frozen porcine spines. They were cut into specimens nominally 3.5mm wide by 7 mm long and tabs of the deep and superficial layers were removed from opposite ends of the specimens so that a 4.5-5.0mm long intact interface remained between the lamellae. Specimens were mounted in a BioTester tensile instrument using BioRake attachments having 5 sharpened points side-by-side, and they were strained at 2%/s. Force-time curves were obtained and, using tracking software, a detailed map was made of the time course of the displacements within the specimens. Extensibility of the lamellae themselves was found to substantially complicate interpretation of the data. The experiments, together with mathematical analyses and finite element models, show that much of the shear load is transferred between lamellae at the ends of the bonded region, a finding of clinical importance. The inter-lamellae bond was found to have a peak strength of 0.30 ± 0.05 N/mm of specimen width (not to be confused with lap length), and the remarkable ability to carry substantial load even when lamellae had displaced up to 10mm relative to each other.

Identifying same-cell contours in image stacks: a key step in making 3D reconstructions.

Identification of contours belonging to the same cell is a crucial step in the analysis of confocal stacks and other image sets in which cell outlines are visible, and it is central to the making of 3D cell reconstructions. When the cells are close packed, the contour grouping problem is more complex than that found in medical imaging, for example, because there are multiple regions of interest, the regions are not separable from each other by an identifiable background and regions cannot be distinguished by intensity differences. Here, we present an algorithm that uses three primary metrics-overlap of contour areas in adjacent images, co-linearity of the centroids of these areas across three images in a stack, and cell taper-to assign cells to groups. Decreasing thresholds are used to successively assign contours whose membership is less obvious. In a final step, remaining contours are assigned to existing groups by setting all thresholds to zero and groups having strong hour-glass shapes are partitioned. When applied to synthetic data from isotropic model aggregates, a curved model epithelium in which the long axes of the cells lie at all possible angles to the transection plane, and a confocal image stack, algorithm assignments were between 97 and 100% accurate in sets having at least four contours per cell. The algorithm is not particularly sensitive to the thresholds used, and a single set of parameters was used for all of the tests. The algorithm, which could be extended to time-lapse data, solves a key problem in the translation of image data into cell information.