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.

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.

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.

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.

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'.

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.

How computational models can help unlock biological systems.

With computation models playing an ever increasing role in the advancement of science, it is important that researchers understand what it means to model something; recognize the implications of the conceptual, mathematical and algorithmic steps of model construction; and comprehend what models can and cannot do. Here, we use examples to show that models can serve a wide variety of roles, including hypothesis testing, generating new insights, deepening understanding, suggesting and interpreting experiments, tracing chains of causation, doing sensitivity analyses, integrating knowledge, and inspiring new approaches. We show that models can bring together information of different kinds and do so across a range of length scales, as they do in multi-scale, multi-faceted embryogenesis models, some of which connect gene expression, the cytoskeleton, cell properties, tissue mechanics, morphogenetic movements and phenotypes. Models cannot replace experiments nor can they prove that particular mechanisms are at work in a given situation. But they can demonstrate whether or not a proposed mechanism is sufficient to produce an observed phenomenon. Although the examples in this article are taken primarily from the field of embryo mechanics, most of the arguments and discussion are applicable to any form of computational modelling.

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.

On the origins of the mitotic shift in proliferating cell layers.

BACKGROUND: During plant and animal development, monolayer cell sheets display a stereotyped distribution of polygonal cell shapes. In interphase cells these shapes range from quadrilaterals to decagons, with a robust average of six sides per cell. In contrast, the subset of cells in mitosis exhibits a distinct distribution with an average of seven sides. It remains unclear whether this 'mitotic shift' reflects a causal relationship between increased polygonal sidedness and increased division likelihood, or alternatively, a passive effect of local proliferation on cell shape. METHODS: We use a combination of probabilistic analysis and mathematical modeling to predict the geometry of mitotic polygonal cells in a proliferating cell layer. To test these predictions experimentally, we use Flp-Out stochastic labeling in the Drosophila wing disc to induce single cell clones, and confocal imaging to quantify the polygonal topologies of these clones as a function of cellular age. For a more generic test in an idealized cell layer, we model epithelial sheet proliferation in a finite element framework, which yields a computationally robust, emergent prediction of the mitotic cell shape distribution. RESULTS: Using both mathematical and experimental approaches, we show that the mitotic shift derives primarily from passive, non-autonomous effects of mitoses in neighboring cells on each cell's geometry over the course of the cell cycle. Computationally, we predict that interphase cells should passively gain sides over time, such that cells at more advanced stages of the cell cycle will tend to have a larger number of neighbors than those at earlier stages. Validating this prediction, experimental analysis of randomly labeled epithelial cells in the Drosophila wing disc demonstrates that labeled cells exhibit an age-dependent increase in polygonal sidedness. Reinforcing these data, finite element simulations of epithelial sheet proliferation demonstrate in a generic framework that passive side-gaining is sufficient to generate a mitotic shift. CONCLUSIONS: Taken together, our results strongly suggest that the mitotic shift reflects a time-dependent accumulation of shared cellular interfaces over the course of the cell cycle. These results uncover fundamental constraints on the relationship between cell shape and cell division that should be general in adherent, polarized cell layers.

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.

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.

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.

Measuring the modulus of silicone hydrogel contact lenses.

PURPOSE: The purpose of this study is to demonstrate a novel method for measuring the modulus of contact lenses in their as-received, variable-thickness form and to determine whether modulus varies with location within commercial lenses and whether it is dependent on lens geometry and temperature. METHODS: The thickness profiles of lenses having powers from -8 diopters (D) to +4 D were measured using a Rehder electronic thickness gauge. Strip-shaped specimens having a width of 5.5 mm were then cut from the lenses. Graphite particles were sprinkled on the specimen surface so that its motions could be tracked using digital image-correlation techniques. The specimens were mounted in a BioTester test system using BioRakes (rather than clamps) and stretched uniaxially until all parts of the lens between the attachment points had elongated by at least 10%. This procedure allowed local modulus values to be determined at 110 locations over the surface of each lens and any property variations within the lenses to be characterized. Tests were performed at 5, 23, and 37°C. RESULTS: Material modulus was found to be essentially constant within any given lens and was independent of the optical power of the lens. Young's Modulus values ranged from 0.3 to 1.9 MPa, depending on the lens manufacturer and product, and some lens materials showed a decrease in modulus with temperature. For the materials tested, those with lower water content had a tendency to exhibit higher moduli. CONCLUSIONS: Testing of the kind reported here is important for assessing the efficacy of current and proposed contact lens materials and designs, especially if such designs make use of variable properties to enhance function or fit.

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.

DRhoGEF2 regulates cellular tension and cell pulsations in the Amnioserosa during Drosophila dorsal closure.

Coordination of apical constriction in epithelial sheets is a fundamental process during embryogenesis. Here, we show that DRhoGEF2 is a key regulator of apical pulsation and constriction of amnioserosal cells during Drosophila dorsal closure. Amnioserosal cells mutant for DRhoGEF2 exhibit a consistent decrease in amnioserosa pulsations whereas overexpression of DRhoGEF2 in this tissue leads to an increase in the contraction time of pulsations. We probed the physical properties of the amnioserosa to show that the average tension in DRhoGEF2 mutant cells is lower than wild-type and that overexpression of DRhoGEF2 results in a tissue that is more solid-like than wild-type. We also observe that in the DRhoGEF2 overexpressing cells there is a dramatic increase of apical actomyosin coalescence that can contribute to the generation of more contractile forces, leading to amnioserosal cells with smaller apical surface than wild-type. Conversely, in DRhoGEF2 mutants, the apical actomyosin coalescence is impaired. These results identify DRhoGEF2 as an upstream regulator of the actomyosin contractile machinery that drives amnioserosa cells pulsations and apical constriction.

A framework for connecting gene expression to morphogenetic movements in embryos.

Although much has been learned about genetic networks, cell mechanics, and whole-embryo mechanics through experimental and computational studies, the challenge of connecting these separate bodies of knowledge into an integrated whole remains. Here, we offer a multiscale biochemical-mechanical framework from which such integration might proceed. We identify components of the framework for which quantitative descriptions are currently available, and use the framework to gain insight into convergent extension and gastrulation--crucial tissue movements that occur in early stage amphibian embryos.

Control of the mitotic cleavage plane by local epithelial topology.

For nearly 150 years, it has been recognized that cell shape strongly influences the orientation of the mitotic cleavage plane (e.g., Hofmeister, 1863). However, we still understand little about the complex interplay between cell shape and cleavage-plane orientation in epithelia, where polygonal cell geometries emerge from multiple factors, including cell packing, cell growth, and cell division itself. Here, using mechanical simulations, we show that the polygonal shapes of individual cells can systematically bias the long-axis orientations of their adjacent mitotic neighbors. Strikingly, analyses of both animal epithelia and plant epidermis confirm a robust and nearly identical correlation between local cell topology and cleavage-plane orientation in vivo. Using simple mathematics, we show that this effect derives from fundamental packing constraints. Our results suggest that local epithelial topology is a key determinant of cleavage-plane orientation, and that cleavage-plane bias may be a widespread property of polygonal cell sheets in plants and animals.

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.

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.