Unified quantitative characterization of epithelial tissue development

The Geometric Basis of Epithelial Convergent Extension

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

Choice of friction coefficient deeply affects tissue behaviour in stochastic epithelial vertex models

To understand the mechanisms that coordinate the formation of biological tissues, the use of numerical implementations is necessary. The complexity of such models involves many assumptions and parameter choices that result in unpredictable consequences, obstructing the comparison with experimental data. Here, we focus on vertex models, a family of spatial models used extensively to simulate the dynamics of epithelial tissues. Usually, in the literature, the choice of the friction coefficient is not addressed using quasi-static deformation arguments that generally do not apply to realistic scenarios. In this manuscript, we discuss the role that the choice of friction coefficient has on the relaxation times and consequently in the conditions of cell cycle progression and division. We explore the effects that these changes have on the morphology, growth rate and topological transitions of the tissue dynamics. These results provide a deeper understanding of the role that an accurate mechanical description plays in the use of vertex models as inference tools. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Balancing competing effects of tissue growth and cytoskeletal regulation during Drosophila wing disc development

How a developing organ robustly coordinates the cellular mechanics and growth to reach a final size and shape remains poorly understood. Through iterations between experiments and model simulations that include a mechanistic description of interkinetic nuclear migration, we show that the local curvature, height, and nuclear positioning of cells in the Drosophila wing imaginal disc are defined by the concurrent patterning of actomyosin contractility, cell-ECM adhesion, ECM stiffness, and interfacial membrane tension. We show that increasing cell proliferation via different growth-promoting pathways results in two distinct phenotypes. Triggering proliferation through insulin signaling increases basal curvature, but an increase in growth through Dpp signaling and Myc causes tissue flattening. These distinct phenotypic outcomes arise from differences in how each growth pathway regulates the cellular cytoskeleton, including contractility and cell-ECM adhesion. The coupled regulation of proliferation and cytoskeletal regulators is a general strategy to meet the multiple context-dependent criteria defining tissue morphogenesis.

Live imaging basement membrane assembly under the pupal notum epithelium

Basement membranes are sheet-like extracellular matrices containing Collagen IV, and they are conserved across the animal kingdom. Basement membranes usually line the basal surfaces of epithelia, where they contribute to structure, maintenance, and signaling. Although adult epithelia contact basement membranes, in early embryos the epithelia contact basement membranes only after basement membranes are assembled in embryogenesis. In Drosophila , the pupal notum epithelium is a useful model for live imaging epithelial cell behaviors, yet it is unclear when the basement membrane assembles in the pupa, as pupae are undergoing metamorphosis, similar to embryogenesis. To characterize the basement membrane in the pupal notum, we used spinning disk fluorescent microscopy to visualize Collagen IV subunit Vkg-GFP and adherens junction protein p120ctnRFP. Bright punctae of Vkg-GFP were observed in the X-Y plane, possibly representing Vkg-containing cells. We found that a thin continuous Vkg-containing basement membrane was evident at 14 h APF, which became more enriched with Vkg-GFP over the next 6 h, indicating the basement membrane is still assembling during that time. Live imaging of the pupal notum during this time could provide insight into formation, assembly, and repair of the basement membranes.

Wounding increases nuclear ploidy in wound-proximal epidermal cells of the Drosophila pupal notum

After injury, tissues must replace cell mass and genome copy number. The mitotic cycle is one mechanism for replacement, but non-mitotic strategies have been observed in quiescent tissues to restore tissue ploidy after wounding. Here we report that nuclei of the mitotically capable Drosophila pupal notum enlarged following nearby laser ablation. Measuring DNA content, we determined that nuclei within 100 µm of a laser-wound increased their ploidy to ~8C, consistent with one extra S-phase. These data indicate non-mitotic repair strategies are not exclusively utilized by quiescent tissues and may be an underexplored wound repair strategy in mitotic tissues.

Patterned apoptosis has an instructive role for local growth and tissue shape regulation in a fast-growing epithelium

What regulates organ size and shape remains one fundamental mystery of modern biology. Research in this area has primarily focused on deciphering the regulation in time and space of growth and cell division, while the contribution of cell death has been overall neglected. This includes studies of the Drosophila wing, one of the best-characterized systems for the study of growth and patterning, undergoing massive growth during larval stage and important morphogenetic remodeling during pupal stage. So far, it has been assumed that cell death was relatively neglectable in this tissue both during larval stage and pupal stage, and as a result, the pattern of growth was usually attributed to the distribution of cell division. Here, using systematic mapping and registration combined with quantitative assessment of clone size and disappearance as well as live imaging, we outline a persistent pattern of cell death and clone elimination emerging in the larval wing disc and persisting during pupal wing morphogenesis. Local variation of cell death is associated with local variation of clone size, pointing to an impact of cell death on local growth that is not fully compensated by proliferation. Using morphometric analyses of adult wing shape and genetic perturbations, we provide evidence that patterned death locally and globally affects adult wing shape and size. This study describes a roadmap for precise assessment of the contribution of cell death to tissue shape and outlines an important instructive role of cell death in modulating quantitatively local growth and morphogenesis of a fast-growing tissue.

Enabling neighbour labelling: using synthetic biology to explore how cells influence their neighbours

Cell-cell interactions are central to development, but exploring how a change in any given cell relates to changes in the neighbour of that cell can be technically challenging. Here, we review recent developments in synthetic biology and image analysis that are helping overcome this problem. We highlight the opportunities presented by these advances and discuss opportunities and limitations in applying them to developmental model systems.

Core PCP mutations affect short-time mechanical properties but not tissue morphogenesis in the Drosophila pupal wing

How morphogenetic movements are robustly coordinated in space and time is a fundamental open question in biology. We study this question using the wing of Drosophila melanogaster, an epithelial tissue that undergoes large-scale tissue flows during pupal stages. Previously, we showed that pupal wing morphogenesis involves both cellular behaviors that allow relaxation of mechanical tissue stress, as well as cellular behaviors that appear to be actively patterned (Etournay et al., 2015). Here, we show that these active cellular behaviors are not guided by the core planar cell polarity (PCP) pathway, a conserved signaling system that guides tissue development in many other contexts. We find no significant phenotype on the cellular dynamics underlying pupal morphogenesis in mutants of core PCP. Furthermore, using laser ablation experiments, coupled with a rheological model to describe the dynamics of the response to laser ablation, we conclude that while core PCP mutations affect the fast timescale response to laser ablation they do not significantly affect overall tissue mechanics. In conclusion, our work shows that cellular dynamics and tissue shape changes during Drosophila pupal wing morphogenesis do not require core PCP as an orientational guiding cue.

An archetype and scaling of developmental tissue dynamics across species

Morphometric studies have revealed the existence of simple geometric relationships among various animal shapes. However, we have little knowledge of the mathematical principles behind the morphogenetic dynamics that form the organ/body shapes of different species. Here, we address this issue by focusing on limb morphogenesis in Gallus gallus domesticus (chicken) and Xenopus laevis (African clawed frog). To compare the deformation dynamics between tissues with different sizes/shapes as well as their developmental rates, we introduce a species-specific rescaled spatial coordinate and a common clock necessary for cross-species synchronization of developmental times. We find that tissue dynamics are well conserved across species under this spacetime coordinate system, at least from the early stages of development through the phase when basic digit patterning is established. For this developmental period, we also reveal that the tissue dynamics of both species are mapped with each other through a time-variant linear transformation in real physical space, from which hypotheses on a species-independent archetype of tissue dynamics and morphogenetic scaling are proposed.

Epithelial homeostasis: Cell size shapes cell fate

A new study shows that cell size, in conjunction with specific signaling pathways, controls apoptosis within developing tissues. Cells with smaller sizes and relatively smaller sizes compared to their neighbors exhibit an increased likelihood of undergoing apoptosis. These processes are regulated by the Hippo/YAP and Notch pathways, respectively.

Epithelial apoptotic pattern emerges from global and local regulation by cell apical area

Geometry is a fundamental attribute of biological systems, and it underlies cell and tissue dynamics. Cell geometry controls cell-cycle progression and mitosis and thus modulates tissue development and homeostasis. In sharp contrast and despite the extensive characterization of the genetic mechanisms of caspase activation, we know little about whether and how cell geometry controls apoptosis commitment in developing tissues. Here, we combined multiscale time-lapse microscopy of developing Drosophila epithelium, quantitative characterization of cell behaviors, and genetic and mechanical perturbations to determine how apoptosis is controlled during epithelial tissue development. We found that early in cell lives and well before extrusion, apoptosis commitment is linked to two distinct geometric features: a small apical area compared with other cells within the tissue and a small relative apical area with respect to the immediate neighboring cells. We showed that these global and local geometric characteristics are sufficient to recapitulate the tissue-scale apoptotic pattern. Furthermore, we established that the coupling between these two geometric features and apoptotic cells is dependent on the Hippo/YAP and Notch pathways. Overall, by exploring the links between cell geometry and apoptosis commitment, our work provides important insights into the spatial regulation of cell death in tissues and improves our understanding of the mechanisms that control cell number and tissue size.

Wound-Induced Syncytia Outpace Mononucleate Neighbors during Drosophila Wound Repair

All organisms have evolved to respond to injury. Cell behaviors like proliferation, migration, and invasion replace missing cells and close wounds. However, the role of other wound-induced cell behaviors is not understood, including the formation of syncytia (multinucleated cells). Wound-induced epithelial syncytia were first reported around puncture wounds in post-mitotic Drosophila epidermal tissues, but have more recently been reported in mitotically competent tissues such as the Drosophila pupal epidermis and zebrafish epicardium. The presence of wound-induced syncytia in mitotically active tissues suggests that syncytia offer adaptive benefits, but it is unknown what those benefits are. Here, we use in vivo live imaging to analyze wound-induced syncytia in mitotically competent Drosophila pupae. We find that almost half the epithelial cells near a wound fuse to form large syncytia. These syncytia use several routes to speed wound repair: they outpace diploid cells to complete wound closure; they reduce cell intercalation during wound closure; and they pool the resources of their component cells to concentrate them toward the wound. In addition to wound healing, these properties of syncytia are likely to contribute to their roles in development and pathology.

DISSECT is a tool to segment and explore cell and tissue mechanics in highly deformed 3D epithelia

Understanding morphogenesis strongly relies on the characterization of tissue topology and mechanical properties deduced from imaging data. The development of new imaging techniques offers the possibility to go beyond the analysis of mostly flat surfaces and image and analyze complex tissue organization in depth. An important bottleneck in this field is the need to analyze imaging datasets and extract quantifications not only of cell and tissue morphology but also of the cytoskeletal network's organization in an automatized way. Here, we describe a method, called DISSECT, for DisPerSE (Discrete Persistent Structure Extractor)-based Segmentation and Exploration of Cells and Tissues, that offers the opportunity to extract automatically, in strongly deformed epithelia, a precise characterization of the spatial organization of a given cytoskeletal network combined with morphological quantifications in highly remodeled three-dimensional (3D) epithelial tissues. We believe that this method, applied here to Drosophila tissues, will be of general interest in the expanding field of morphogenesis and tissue biomechanics.

Constitutive model for the rheology of biological tissue

The rheology of biological tissue is key to processes such as embryo development, wound healing, and cancer metastasis. Vertex models of confluent tissue monolayers have uncovered a spontaneous liquid-solid transition tuned by cell shape; and a shear-induced solidification transition of an initially liquidlike tissue. Alongside this jamming/unjamming behavior, biological tissue also displays an inherent viscoelasticity, with a slow time and rate-dependent mechanics. With this motivation, we combine simulations and continuum theory to examine the rheology of the vertex model in nonlinear shear across a full range of shear rates from quastistatic to fast, elucidating its nonlinear stress-strain curves after the inception of shear of finite rate, and its steady state flow curves of stress as a function of strain rate. We formulate a rheological constitutive model that couples cell shape to flow and captures both the tissue solid-liquid transition and its rich linear and nonlinear rheology.

Apical size and deltaA expression predict adult neural stem cell decisions along lineage progression

The maintenance of neural stem cells (NSCs) in the adult brain depends on their activation frequency and division mode. Using long-term intravital imaging of NSCs in the zebrafish adult telencephalon, we reveal that apical surface area and expression of the Notch ligand DeltaA predict these NSC decisions. deltaA-negative NSCs constitute a bona fide self-renewing NSC pool and systematically engage in asymmetric divisions generating a self-renewing deltaAneg daughter, which regains the size and behavior of its mother, and a neurogenic deltaApos daughter, eventually engaged in neuronal production following further quiescence-division phases. Pharmacological and genetic manipulations of Notch, DeltaA, and apical size further show that the prediction of activation frequency by apical size and the asymmetric divisions of deltaAneg NSCs are functionally independent of Notch. These results provide dynamic qualitative and quantitative readouts of NSC lineage progression in vivo and support a hierarchical organization of NSCs in differently fated subpopulations.

Robust statistical properties of T1 transitions in a multi-phase field model of cell monolayers

Large-scale tissue deformation which is fundamental to tissue development hinges on local cellular rearrangements, such as T1 transitions. In the realm of the multi-phase field model, we analyse the statistical and dynamical properties of T1 transitions in a confluent monolayer. We identify an energy profile that is robust to changes in several model parameters. It is characterized by an asymmetric profile with a fast increase in energy before the T1 transition and a sudden drop after the T1 transition, followed by a slow relaxation. The latter being a signature of the fluidity of the cell monolayer. We show that T1 transitions are sources of localised large deformation of the cells undergoing the neighbour exchange, and they induce other T1 transitions in the nearby cells leading to a chaining of events that propagate local cell deformation to large scale tissue flows.

Zasp52 strengthens whole embryo tissue integrity through supracellular actomyosin networks

During morphogenesis, large-scale changes of tissue primordia are coordinated across an embryo. In Drosophila, several tissue primordia and embryonic regions are bordered or encircled by supracellular actomyosin cables, junctional actomyosin enrichments networked between many neighbouring cells. We show that the single Drosophila Alp/Enigma-family protein Zasp52, which is most prominently found in Z-discs of muscles, is a component of many supracellular actomyosin structures during embryogenesis, including the ventral midline and the boundary of the salivary gland placode. We reveal that Zasp52 contains within its central coiled-coil region a type of actin-binding motif usually found in CapZbeta proteins, and this domain displays actin-binding activity. Using endogenously-tagged lines, we identify that Zasp52 interacts with junctional components, including APC2, Polychaetoid and Sidekick, and actomyosin regulators. Analysis of zasp52 mutant embryos reveals that the severity of the embryonic defects observed scales inversely with the amount of functional protein left. Large tissue deformations occur where actomyosin cables are found during embryogenesis, and in vivo and in silico analyses suggest a model whereby supracellular Zasp52-containing cables aid to insulate morphogenetic changes from one another.

Systematic analysis of RhoGEF/GAP localizations uncovers regulators of mechanosensing and junction formation during epithelial cell division

Cell proliferation is central to epithelial tissue development, repair, and homeostasis. During cell division, small RhoGTPases control both actomyosin dynamics and cell-cell junction remodeling to faithfully segregate the genome while maintaining tissue polarity and integrity. To decipher the mechanisms of RhoGTPase spatiotemporal regulation during epithelial cell division, we generated a transgenic fluorescently tagged library for the 48 Drosophila Rho guanine exchange factors (RhoGEFs) and GTPase-activating proteins (GAPs), and we systematically characterized their endogenous distributions by time-lapse microscopy. Therefore, we unveiled candidate regulators of the interplay between actomyosin and junctional dynamics during epithelial cell division. Building on these findings, we established that the conserved RhoGEF Cysts and RhoGEF4 play sequential and distinct roles to couple cytokinesis with de novo junction formation. During ring contraction, Cysts via Rho1 participates in the neighbor mechanosensing response, promoting daughter-daughter cell membrane juxtaposition in preparation to de novo junction formation. Subsequently and upon midbody formation, RhoGEF4 via Rac acts in the dividing cell to ensure the withdrawal of the neighboring cell membranes, thus controlling de novo junction length and cell-cell arrangements upon cytokinesis. Altogether, our findings delineate how the RhoGTPases Rho and Rac are locally and temporally activated during epithelial cytokinesis, highlighting the RhoGEF/GAP library as a key resource to understand the broad range of biological processes regulated by RhoGTPases.

Curvature gradient drives polarized tissue flow in the Drosophila embryo

Tissue flow during morphogenesis is commonly driven by local constriction of cell cortices, which is caused by the activation of actomyosin contractility. This can lead to long-range flows due to tissue viscosity. However, in the absence of cell-intrinsic polarized forces or polarity in forces external to the tissue, these flows must be symmetric and centered around the region of contraction. Polarized tissue flows have been previously demonstrated to arise from the coupling of such contractile flows to points of increased friction or adhesion to external structures. However, we show with experiments and modeling that the onset of polarized tissue flow in early Drosophila morphogenesis occurs independent of adhesion and is instead driven by a geometric coupling of apical actomyosin contractility to tissue curvature. Particularly, the onset of polarized flow is driven by a mismatch between the position of apical myosin activation and the position of peak curvature at the posterior pole of the embryo. Our work demonstrates how genetic and geometric information inherited from the mother interact to create polarized flow during embryo morphogenesis.

Attachment and detachment of cortical myosin regulates cell junction exchange during cell rearrangement in the Drosophila wing epithelium

Epithelial cells remodel cell adhesion and change their neighbors to shape a tissue. This cellular rearrangement proceeds in three steps: the shrinkage of a junction, exchange of junctions, and elongation of the newly generated junction. Herein, by combining live imaging and physical modeling, we showed that the formation of myosin-II (myo-II) cables around the cell vertices underlies the exchange of junctions in the Drosophila wing epithelium. The local and transient detachment of myo-II from the cell cortex is regulated by the LIM domain-containing protein Jub and the tricellular septate junction protein M6. Moreover, we found that M6 shifts to the adherens junction plane on jub RNAi and that Jub is persistently retained at reconnecting junctions in m6 RNAi cells. This interplay between Jub and M6 can depend on the junction length and thereby couples the detachment of cortical myo-II cables and the shrinkage/elongation of the junction during cell rearrangement. Furthermore, we developed a mechanical model based on the wetting theory and clarified how the physical properties of myo-II cables are integrated with the junction geometry to induce the transition between the attached and detached states and support the unidirectionality of cell rearrangement. Collectively, this study elucidates the orchestration of geometry, mechanics, and signaling for exchanging junctions.

Resolving morphogenesis into quantifiable cell behaviours

Morphogenesis is extremely diverse, but its systematic quantification to determine the physical mechanisms that produce different phenotypes is possible by quantifying the underlying cell behaviours. These are limited and definable: they consist of cell proliferation, orientation of cell division, cell rearrangement, directional matrix production, cell addition/subtraction and cell size/shape change. Although minor variations in these categories are possible, in sum they capture all possible morphogenetic behaviours. This article summarises these processes, discusses their measurement, and highlights some salient examples.

aPKC regulates apical constriction to prevent tissue rupture in the Drosophila follicular epithelium

Apical-basal polarity is an essential epithelial trait controlled by the evolutionarily conserved PAR-aPKC polarity network. Dysregulation of polarity proteins disrupts tissue organization during development and in disease, but the underlying mechanisms are unclear due to the broad implications of polarity loss. Here, we uncover how Drosophila aPKC maintains epithelial architecture by directly observing tissue disorganization after fast optogenetic inactivation in living adult flies and ovaries cultured ex vivo. We show that fast aPKC perturbation in the proliferative follicular epithelium produces large epithelial gaps that result from increased apical constriction, rather than loss of apical-basal polarity. Accordingly, we can modulate the incidence of epithelial gaps by increasing and decreasing actomyosin-driven contractility. We traced the origin of these large epithelial gaps to tissue rupture next to dividing cells. Live imaging shows that aPKC perturbation induces apical constriction in non-mitotic cells within minutes, producing pulling forces that ultimately detach dividing and neighboring cells. We further demonstrate that epithelial rupture requires a global increase of apical constriction, as it is prevented by the presence of non-constricting cells. Conversely, a global induction of apical tension through light-induced recruitment of RhoGEF2 to the apical side is sufficient to produce tissue rupture. Hence, our work reveals that the roles of aPKC in polarity and actomyosin regulation are separable and provides the first in vivo evidence that excessive tissue stress can break the epithelial barrier during proliferation.

Tissue flow regulates planar cell polarity independently of the Frizzled core pathway

Planar cell polarity (PCP) regulates the orientation of external structures. A core group of proteins that includes Frizzled forms the heart of the PCP regulatory system. Other PCP mechanisms that are independent of the core group likely exist, but their underlying mechanisms are elusive. Here, we show that tissue flow is a mechanism governing core group-independent PCP on the Drosophila notum. Loss of core group function only slightly affects bristle orientation in the adult central notum. This near-normal PCP results from tissue flow-mediated rescue of random bristle orientation during the pupal stage. Manipulation studies suggest that tissue flow can orient bristles in the opposite direction to the flow. This process is independent of the core group and implies that the apical extracellular matrix functions like a "comb" to align bristles. Our results reveal the significance of cooperation between tissue dynamics and extracellular substances in PCP establishment.

Extracting multiple surfaces from 3D microscopy images in complex biological tissues with the Zellige software tool

Background

Efficient tools allowing the extraction of 2D surfaces from 3D-microscopy data are essential for studies aiming to decipher the complex cellular choreography through which epithelium morphogenesis takes place during development. Most existing methods allow for the extraction of a single and smooth manifold of sufficiently high signal intensity and contrast, and usually fail when the surface of interest has a rough topography or when its localization is hampered by other surrounding structures of higher contrast. Multiple surface segmentation entails laborious manual annotations of the various surfaces separately.

Results

As automating this task is critical in studies involving tissue-tissue or tissue-matrix interaction, we developed the Zellige software, which allows the extraction of a non-prescribed number of surfaces of varying inclination, contrast, and texture from a 3D image. The tool requires the adjustment of a small set of control parameters, for which we provide an intuitive interface implemented as a Fiji plugin.

Conclusions

As a proof of principle of the versatility of Zellige, we demonstrate its performance and robustness on synthetic images and on four different types of biological samples, covering a wide range of biological contexts.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01378-0.

Instabilities and Geometry of Growing Tissues

We present a covariant continuum formulation of a generalized two-dimensional vertexlike model of epithelial tissues which describes tissues with different underlying geometries, and allows for an analytical macroscopic description. Using a geometrical approach and out-of-equilibrium statistical mechanics, we calculate both mechanical and dynamical instabilities of a tissue, and their dependences on various variables, including activity, and cell-shape heterogeneity (disorder). We show how both plastic cellular rearrangements and the tissue elastic response depend on the existence of mechanical residual stresses at the cellular level. Even freely growing tissues may exhibit a growth instability depending on the intrinsic proliferation rate. Our main result is an explicit calculation of the cell pressure in a homeostatic state of a confined growing tissue. We show that the homeostatic pressure can be negative and depends on the existence of mechanical residual stresses. This geometric model allows us to sort out elastic and plastic effects in a growing, flowing, tissue.

Image-based parameter inference for epithelial mechanics

Measuring mechanical parameters in tissues, such as the elastic modulus of cell-cell junctions, is essential to decipher the mechanical control of morphogenesis. However, their in vivo measurement is technically challenging. Here, we formulated an image-based statistical approach to estimate the mechanical parameters of epithelial cells. Candidate mechanical models are constructed based on force-cell shape correlations obtained from image data. Substitution of the model functions into force-balance equations at the cell vertex leads to an equation with respect to the parameters of the model, by which one can estimate the parameter values using a least-squares method. A test using synthetic data confirmed the accuracy of parameter estimation and model selection. By applying this method to Drosophila epithelial tissues, we found that the magnitude and orientation of feedback between the junction tension and shrinkage, which are determined by the spring constant of the junction, were correlated with the elevation of tension and myosin-II on shrinking junctions during cell rearrangement. Further, this method clarified how alterations in tissue polarity and stretching affect the anisotropy in tension parameters. Thus, our method provides a novel approach to uncovering the mechanisms governing epithelial morphogenesis.

Two types of critical cell density for mechanical elimination of abnormal cell clusters from epithelial tissue

Recent technological advances in high-resolution imaging and artificial modulation of genetic functions at different times and regions have enabled direct observations of the formation and elimination of abnormal cell populations. A recent trend in cell competition research is the incorporation of cell mechanics. In different tissues and species, abnormal cells developing in epithelial tissues are mechanically eliminated by cell contraction via actomyosin accumulation at the interface between normal and abnormal cells. This mechanical cell elimination process has attracted attention as a potential universal defense mechanism. Here, we theoretically examined the conditions for mechanical elimination of growing abnormal cell populations. Simulations and mathematical analyses using a vertex dynamics model revealed two types of critical cell density associated with mechanical elimination of abnormal cell clusters. One is a subtype of homeostatic density, in which the frequencies of spontaneous mechanical cell elimination and proliferation are balanced, even if no explicit dependence of proliferation or apoptosis on the cell density is assumed. This density is related to the mechanical stability of a single cell. The other is density related to mechanical stability as a cell population under external pressure. Both density types are determined by tissue mechanical properties. In solid tissues, the former type is reached first as the intensity of interfacial contraction increases, and it functions as a critical density. On the other hand, the latter type becomes critical when tissues are highly fluid. The derived analytical solution explicitly reveals the dependence of critical contractile force and density on different parameters. We also found a negative correlation between the proliferation rate of abnormal cells and the likelihood of the abnormal cell population expanding by escaping elimination. This is counterintuitive because in the context of cell competition, fast-growing cell populations generally win. These findings provide new insight into, and interpretation of, the results from experimental studies.

High-resolution imaging techniques have revealed that abnormal cells developing in epithelial tissues are mechanically eliminated via contraction at the interface between the abnormal cells and normal surrounding cells. This phenomenon is seen in various species and tissues and thus is regarded as a primitive defense system against precancerous cells common to all animals. For comprehensive understanding of this potential defense system, we derived mathematical conditions to achieve mechanical elimination of growing abnormal cell populations. We identified two characteristic cell density types associated with successful mechanical elimination of abnormal cell clusters. Both are determined by tissue physical properties, and the smaller of the two functions as a critical density above which abnormal cell populations cannot exist. We also found a counterintuitive phenomenon in which slower proliferation of abnormal cells promotes their growth as a population. Our results will help elucidate the mechanisms of intrinsic tissue defenses against cancer from the perspective of cell mechanics.

Left-right symmetry of zebrafish embryos requires somite surface tension

The body axis of vertebrate embryos is periodically segmented into bilaterally symmetric pairs of somites^1,2. The anteroposterior length of somites, their position and left-right symmetry are thought to be molecularly determined before somite morphogenesis^3,4. Here we show that, in zebrafish embryos, initial somite anteroposterior lengths and positions are imprecise and, consequently, many somite pairs form left-right asymmetrically. Notably, these imprecisions are not left unchecked and we find that anteroposterior lengths adjust within an hour after somite formation, thereby increasing morphological symmetry. We find that anteroposterior length adjustments result entirely from changes in somite shape without change in somite volume, with changes in anteroposterior length being compensated by corresponding changes in mediolateral length. The anteroposterior adjustment mechanism is facilitated by somite surface tension, which we show by comparing in vivo experiments and in vitro single-somite explant cultures using a mechanical model. Length adjustment is inhibited by perturbation of molecules involved in surface tension, such as integrin and fibronectin. By contrast, the adjustment mechanism is unaffected by perturbations to the segmentation clock, therefore revealing a distinct process that influences morphological segment lengths. We propose that tissue surface tension provides a general mechanism to adjust shapes and ensure precision and symmetry of tissues in developing embryos.

Shear-Driven Solidification and Nonlinear Elasticity in Epithelial Tissues

Biological processes, from morphogenesis to tumor invasion, spontaneously generate shear stresses inside living tissue. The mechanisms that govern the transmission of mechanical forces in epithelia and the collective response of the tissue to bulk shear deformations remain, however, poorly understood. Using a minimal cell-based computational model, we investigate the constitutive relation of confluent tissues under simple shear deformation. We show that an initially undeformed fluidlike tissue acquires finite rigidity above a critical applied strain. This is akin to the shear-driven rigidity observed in other soft matter systems. Interestingly, shear-driven rigidity can be understood by a critical scaling analysis in the vicinity of the second order critical point that governs the liquid-solid transition of the undeformed system. We further show that a solidlike tissue responds linearly only to small strains and but then switches to a nonlinear response at larger stains, with substantial stiffening. Finally, we propose a mean-field formulation for cells under shear that offers a simple physical explanation of shear-driven rigidity and nonlinear response in a tissue.

ECM degradation in the Drosophila abdominal epidermis initiates tissue growth that ceases with rapid cell-cycle exit

During development, multicellular organisms undergo stereotypical patterns of tissue growth in space and time. How developmental growth is orchestrated remains unclear, largely due to the difficulty of observing and quantitating this process in a living organism. Drosophila histoblast nests are small clusters of progenitor epithelial cells that undergo extensive growth to give rise to the adult abdominal epidermis and are amenable to live imaging. Our quantitative analysis of histoblast proliferation and tissue mechanics reveals that tissue growth is driven by cell divisions initiated through basal extracellular matrix degradation by matrix metalloproteases secreted by the neighboring larval epidermal cells. Laser ablations and computational simulations show that tissue mechanical tension does not decrease as the histoblasts fill the abdominal epidermal surface. During tissue growth, the histoblasts display oscillatory cell division rates until growth termination occurs through the rapid emergence of G0/G1 arrested cells, rather than a gradual increase in cell-cycle time as observed in other systems such as the Drosophila wing and mouse postnatal epidermis. Different developing tissues can therefore achieve their final size using distinct growth termination strategies. Thus, adult abdominal epidermal development is characterized by changes in the tissue microenvironment and a rapid exit from the cell cycle.

Active T1 transitions in cellular networks

In amorphous solids as in tissues, neighbor exchanges can relax local stresses and allow the material to flow. In this paper, we use an anisotropic vertex model to study T1 rearrangements in polygonal cellular networks. We consider two different physical realizations of the active anisotropic stresses: (i) anisotropic bond tension and (ii) anisotropic cell stress. Interestingly, the two types of active stress lead to patterns of relative orientation of T1 transitions and cell elongation that are different. Our work suggests that these two realizations of anisotropic active stresses can be observed in vivo. We describe and explain these results through the lens of a continuum description of the tissue as an anisotropic active material. We furthermore discuss the energetics of the dynamic tissue and express the energy balance in terms of internal elastic energy, mechanical work, chemical work and heat. This allows us to define active T1 transitions that can perform mechanical work while consuming chemical energy.

The wing imaginal disc

The Drosophila wing imaginal disc is a tissue of undifferentiated cells that are precursors of the wing and most of the notum of the adult fly. The wing disc first forms during embryogenesis from a cluster of ∼30 cells located in the second thoracic segment, which invaginate to form a sac-like structure. They undergo extensive proliferation during larval stages to form a mature larval wing disc of ∼35,000 cells. During this time, distinct cell fates are assigned to different regions, and the wing disc develops a complex morphology. Finally, during pupal stages the wing disc undergoes morphogenetic processes and then differentiates to form the adult wing and notum. While the bulk of the wing disc comprises epithelial cells, it also includes neurons and glia, and is associated with tracheal cells and muscle precursor cells. The relative simplicity and accessibility of the wing disc, combined with the wealth of genetic tools available in Drosophila, have combined to make it a premier system for identifying genes and deciphering systems that play crucial roles in animal development. Studies in wing imaginal discs have made key contributions to many areas of biology, including tissue patterning, signal transduction, growth control, regeneration, planar cell polarity, morphogenesis, and tissue mechanics.

Live 3D imaging and mapping of shear stresses within tissues using incompressible elastic beads

To investigate the role of mechanical constraints in morphogenesis and development, we have developed a pipeline of techniques based on incompressible elastic sensors. These techniques combine the advantages of incompressible liquid droplets, which have been used as precise in situ shear stress sensors, and of elastic compressible beads, which are easier to tune and to use. Droplets of a polydimethylsiloxane mix, made fluorescent through specific covalent binding to a rhodamin dye, are produced by a microfluidics device. The elastomer rigidity after polymerization is adjusted to the tissue rigidity. Its mechanical properties are carefully calibrated in situ, for a sensor embedded in a cell aggregate submitted to uniaxial compression. The local shear stress tensor is retrieved from the sensor shape, accurately reconstructed through an active contour method. In vitro, within cell aggregates, and in vivo, in the prechordal plate of the zebrafish embryo during gastrulation, our pipeline of techniques demonstrates its efficiency to directly measure the three dimensional shear stress repartition within a tissue.

Multi-view confocal microscopy enables multiple organ and whole organism live-imaging

Understanding how development is coordinated in multiple tissues and gives rise to fully functional organs or whole organisms necessitates microscopy tools. Over the last decade numerous advances have been made in live-imaging, enabling high resolution imaging of whole organisms at cellular resolution. Yet, these advances mainly rely on mounting the specimen in agarose or aqueous solutions, precluding imaging of organisms whose oxygen uptake depends on ventilation. Here, we implemented a multi-view multi-scale microscopy strategy based on confocal spinning disk microscopy, called Multi-View confocal microScopy (MuViScopy). MuViScopy enables live-imaging of multiple organs with cellular resolution using sample rotation and confocal imaging without the need of sample embedding. We illustrate the capacity of MuViScopy by live-imaging Drosophila melanogaster pupal development throughout metamorphosis, highlighting how internal organs are formed and multiple organ development is coordinated. We foresee that MuViScopy will open the path to better understand developmental processes at the whole organism scale in living systems that require gas exchange by ventilation.

Rapid and robust optogenetic control of gene expression in Drosophila

Deciphering gene function requires the ability to control gene expression in space and time. Binary systems such as the Gal4/UAS provide a powerful means to modulate gene expression and to induce loss or gain of function. This is best exemplified in Drosophila, where the Gal4/UAS system has been critical to discover conserved mechanisms in development, physiology, neurobiology, and metabolism, to cite a few. Here we describe a transgenic light-inducible Gal4/UAS system (ShineGal4/UAS) based on Magnet photoswitches. We show that it allows efficient, rapid, and robust activation of UAS-driven transgenes in different tissues and at various developmental stages in Drosophila. Furthermore, we illustrate how ShineGal4 enables the generation of gain and loss-of-function phenotypes at animal, organ, and cellular levels. Thanks to the large repertoire of UAS-driven transgenes, ShineGal4 enriches the Drosophila genetic toolkit by allowing in vivo control of gene expression with high temporal and spatial resolutions.

From heterogeneous morphogenetic fields to homogeneous regions as a step towards understanding complex tissue dynamics

Within developing tissues, cell proliferation, cell motility and other cell behaviors vary spatially, and this variability gives a complexity to the morphogenesis. Recently, novel formalisms have been developed to quantify tissue deformation and underlying cellular processes. A major challenge for the study of morphogenesis now is to objectively define tissue sub-regions exhibiting different dynamics. Here, we propose a method to automatically divide a tissue into regions where the local deformation rate is homogeneous. This was achieved by several steps including image segmentation, clustering and region boundary smoothing. We illustrate the use of the pipeline using a large dataset obtained during the metamorphosis of the Drosophila pupal notum. We also adapt it to determine regions in which the time evolution of the local deformation rate is homogeneous. Finally, we generalize its use to find homogeneous regions for cellular processes such as cell division, cell rearrangement, or cell size and shape changes. We also illustrate it on wing blade morphogenesis. This pipeline will contribute substantially to the analysis of complex tissue shaping, and the biochemical and biomechanical regulations driving tissue morphogenesis.

Configurational fingerprints of multicellular living systems

Cells cooperate as groups to achieve structure and function at the tissue level, during which specific material characteristics emerge. Analogous to phase transitions in classical physics, transformations in the material characteristics of multicellular assemblies are essential for a variety of vital processes including morphogenesis, wound healing, and cancer. In this work, we develop configurational fingerprints of particulate and multicellular assemblies and extract volumetric and shear order parameters based on this fingerprint to quantify the system disorder. Theoretically, these two parameters form a complete and unique pair of signatures for the structural disorder of a multicellular system. The evolution of these two order parameters offers a robust and experimentally accessible way to map the phase transitions in expanding cell monolayers and during embryogenesis and invasion of epithelial spheroids.

RASSF8-mediated transport of Echinoid via the exocyst promotes Drosophila wing elongation and epithelial ordering

Cell-cell junctions are dynamic structures that maintain cell cohesion and shape in epithelial tissues. During development, junctions undergo extensive rearrangements to drive the epithelial remodelling required for morphogenesis. This is particularly evident during axis elongation, where neighbour exchanges, cell-cell rearrangements and oriented cell divisions lead to large-scale alterations in tissue shape. Polarised vesicle trafficking of junctional components by the exocyst complex has been proposed to promote junctional rearrangements during epithelial remodelling, but the receptors that allow exocyst docking to the target membranes remain poorly understood. Here, we show that the adherens junction component Ras Association domain family 8 (RASSF8) is required for the epithelial re-ordering that occurs during Drosophila pupal wing proximo-distal elongation. We identify the exocyst component Sec15 as a RASSF8 interactor. Loss of RASSF8 elicits cytoplasmic accumulation of Sec15 and Rab11-containing vesicles. These vesicles also contain the nectin-like homophilic adhesion molecule Echinoid, the depletion of which phenocopies the wing elongation and epithelial packing defects observed in RASSF8 mutants. Thus, our results suggest that RASSF8 promotes exocyst-dependent docking of Echinoid-containing vesicles during morphogenesis.

Nonlinear rheology of cellular networks

Morphogenesis depends crucially on the complex rheological properties of cell tissues and on their ability to maintain mechanical integrity while rearranging at long times. In this paper, we study the rheology of polygonal cellular networks described by a vertex model in the presence of fluctuations. We use a triangulation method to decompose shear into cell shape changes and cell rearrangements. Considering the steady-state stress under constant shear, we observe nonlinear shear-thinning behavior at all magnitudes of the fluctuations, and an even stronger nonlinear regime at lower values of the fluctuations. We successfully capture this nonlinear rheology by a mean-field model that describes the tissue in terms of cell elongation and cell rearrangements. We furthermore introduce anisotropic active stresses in the vertex model and analyze their effect on rheology. We include this anisotropy in the mean-field model and show that it recapitulates the behavior observed in the simulations. Our work clarifies how tissue rheology is related to stochastic cell rearrangements and provides a simple biophysical model to describe biological tissues. Further, it highlights the importance of nonlinearities when discussing tissue mechanics.

Dynamic changes in epithelial cell packing during tissue morphogenesis

Cell packing - the spatial arrangement of cells - determines the shapes of organs. Recently, investigations of organ development in a variety of model organisms have uncovered cellular mechanisms that are used by epithelial tissues to change cell packing, and thereby their shapes, to generate functional architectures. Here, we review these cellular mechanisms across a wide variety of developmental processes in vertebrates and invertebrates and identify a set of common motifs in the morphogenesis toolbox that, in combination, appear to allow any change in tissue shape. We focus on tissue elongation, folding and invagination, and branching. We also highlight how these morphogenetic processes are achieved by cell-shape changes, cell rearrangements, and oriented cell division. Finally, we describe approaches that have the potential to engineer three-dimensional tissues for both basic science and translational purposes. This review provides a framework for future analyses of how tissues are shaped by the dynamics of epithelial cell packing.

A rapid segmentation method of cell boundary for developing embryos using machine learning with a personal computer

Cell segmentation is crucial in the study of morphogenesis in developing embryos, but it had been limited in its accuracy until machine learning methods for image segmentation like U-Net. However, these methods take too much time. In this study, we provide a rapid method for cell segmentation using machine learning with a personal computer, termed Cell Segmentator using Machine Learning (CSML). CSML took four seconds per image with a personal computer for segmentation on average, much less than time to obtain an image. We observed that F-value of segmentation by CSML was around 0.97, showing better performance than state-of-the-art methods like RACE and watershed in assessing the segmentation of Xenopus ectodermal cells. CSML also showed slightly better performance and faster than other machine learning-based methods such as U-Net. CSML required only one whole embryo image for training a Fully Convolutional Network classifier and only two parameters. To validate its accuracy, we compared CSML to other methods in assessing several indicators of cell shape. We also examined the generality of this approach by measuring its performance of segmentation of independent images. Our data demonstrate the superiority of CSML, and we expect this application to improve efficiency in cell shape studies.

In Vitro Measurements of Cellular Forces and their Importance in the Lung-From the Sub- to the Multicellular Scale

Throughout life, the body is subjected to various mechanical forces on the organ, tissue, and cellular level. Mechanical stimuli are essential for organ development and function. One organ whose function depends on the tightly connected interplay between mechanical cell properties, biochemical signaling, and external forces is the lung. However, altered mechanical properties or excessive mechanical forces can also drive the onset and progression of severe pulmonary diseases. Characterizing the mechanical properties and forces that affect cell and tissue function is therefore necessary for understanding physiological and pathophysiological mechanisms. In recent years, multiple methods have been developed for cellular force measurements at multiple length scales, from subcellular forces to measuring the collective behavior of heterogeneous cellular networks. In this short review, we give a brief overview of the mechanical forces at play on the cellular level in the lung. We then focus on the technological aspects of measuring cellular forces at many length scales. We describe tools with a subcellular resolution and elaborate measurement techniques for collective multicellular units. Many of the technologies described are by no means restricted to lung research and have already been applied successfully to cells from various other tissues. However, integrating the knowledge gained from these multi-scale measurements in a unifying framework is still a major future challenge.

Active flows and deformable surfaces in development

We review progress in active hydrodynamic descriptions of flowing media on curved and deformable manifolds: the state-of-the-art in continuum descriptions of single-layers of epithelial and/or other tissues during development. First, after a brief overview of activity, flows and hydrodynamic descriptions, we highlight the generic challenge of identifying the dependence on dynamical variables of so-called active kinetic coefficients- active counterparts to dissipative Onsager coefficients. We go on to describe some of the subtleties concerning how curvature and active flows interact, and the issues that arise when surfaces are deformable. We finish with a broad discussion around the utility of such theories in developmental biology. This includes limitations to analytical techniques, challenges associated with numerical integration, fitting-to-data and inference, and potential tools for the future, such as discrete differential geometry.

The cell in the age of the genomic revolution: Cell Regulatory Networks

Over the last few years an intense activity in the areas of advanced microscopy and quantitative cell biology has put the focus on the morphogenetic events that shape embryos. The interest in these processes is taking place against the backdrop of genomic studies, particularly of global patterns of gene expression at the level of single cells, which cannot fully account for the way cells build tissues and organs. Here we discuss the need to integrate the activity of genes with that of cells and propose the need to develop a framework, based on cellular processes and cell interactions, that parallels that which has been created for gene activity in the form of Gene Regulatory Networks (GRNs). We begin to do this by suggesting elements for building Cell Regulatory Networks (CRNs). In the same manner that GRNs create schedules of gene expression that result in the emergence of cell fates over time, CRNs create tissues and organs i.e. space. We also suggest how GRNs and CRNs might interact in the building of embryos through feedback loops involving mechanics and tissue tectonics.

Mechanical Patterning in Animal Morphogenesis

Morphogenesis is one of the most remarkable examples of biological pattern formation. Despite substantial progress in the field, we still do not understand the organizational principles responsible for the robust convergence of the morphogenesis process across scales to form viable organisms under variable conditions. Achieving large-scale coordination requires feedback between mechanical and biochemical processes, spanning all levels of organization and relating the emerging patterns with the mechanisms driving their formation. In this review, we highlight the role of mechanics in the patterning process, emphasizing the active and synergistic manner in which mechanical processes participate in developmental patterning rather than merely following a program set by biochemical signals. We discuss the value of applying a coarse-grained approach toward understanding this complex interplay, which considers the large-scale dynamics and feedback as well as complementing the reductionist approach focused on molecular detail. A central challenge in this approach is identifying relevant coarse-grained variables and developing effective theories that can serve as a basis for an integrated framework for understanding this remarkable pattern-formation process. Expected final online publication date for the Annual Review of Cell and Developmental Biology, Volume 37 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

LocalZProjector and DeProj: a toolbox for local 2D projection and accurate morphometrics of large 3D microscopy images

Background

Quantitative imaging of epithelial tissues requires bioimage analysis tools that are widely applicable and accurate. In the case of imaging 3D tissues, a common preprocessing step consists of projecting the acquired 3D volume on a 2D plane mapping the tissue surface. While segmenting the tissue cells is amenable on 2D projections, it is still very difficult and cumbersome in 3D. However, for many specimen and models used in developmental and cell biology, the complex content of the image volume surrounding the epithelium in a tissue often reduces the visibility of the biological object in the projection, compromising its subsequent analysis. In addition, the projection may distort the geometry of the tissue and can lead to strong artifacts in the morphology measurement.

Results

Here we introduce a user-friendly toolbox built to robustly project epithelia on their 2D surface from 3D volumes and to produce accurate morphology measurement corrected for the projection distortion, even for very curved tissues. Our toolbox is built upon two components. LocalZProjector is a configurable Fiji plugin that generates 2D projections and height-maps from potentially large 3D stacks (larger than 40 GB per time-point) by only incorporating signal of the planes with local highest variance/mean intensity, despite a possibly complex image content. DeProj is a MATLAB tool that generates correct morphology measurements by combining the height-map output (such as the one offered by LocalZProjector) and the results of a cell segmentation on the 2D projection, hence effectively deprojecting the 2D segmentation in 3D. In this paper, we demonstrate their effectiveness over a wide range of different biological samples. We then compare its performance and accuracy against similar existing tools.

Conclusions

We find that LocalZProjector performs well even in situations where the volume to project also contains unwanted signal in other layers. We show that it can process large images without a pre-processing step. We study the impact of geometrical distortions on morphological measurements induced by the projection. We measured very large distortions which are then corrected by DeProj, providing accurate outputs.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12915-021-01037-w).

Leaf Shape Diversity: From Genetic Modules to Computational Models

Plant leaves display considerable variation in shape. Here, we introduce key aspects of leaf development, focusing on the morphogenetic basis of leaf shape diversity. We discuss the importance of the genetic control of the amount, duration, and direction of cellular growth for the emergence of leaf form. We highlight how the combined use of live imaging and computational frameworks can help conceptualize how regulated cellular growth is translated into different leaf shapes. In particular, we focus on the morphogenetic differences between simple and complex leaves and how carnivorous plants form three-dimensional insect traps. We discuss how evolution has shaped leaf diversity in the case of complex leaves, by tinkering with organ-wide growth and local growth repression, and in carnivorous plants, by modifying the relative growth of the lower and upper sides of the leaf primordium to create insect-digesting traps.

Modeling and live imaging of mechanical instabilities in the zebrafish aorta during hematopoiesis

All blood cells originate from hematopoietic stem/progenitor cells (HSPCs). HSPCs are formed from endothelial cells (ECs) of the dorsal aorta (DA), via endothelial-to-hematopoietic transition (EHT). The zebrafish is a primary model organism to study the process in vivo. While the role of mechanical stress in controlling gene expression promoting cell differentiation is actively investigated, mechanisms driving shape changes of the DA and individual ECs remain poorly understood. We address this problem by developing a new DA micromechanical model and applying it to experimental data on zebrafish morphogenesis. The model considers the DA as an isotropic tubular membrane subjected to hydrostatic blood pressure and axial stress. The DA evolution is described as a movement in the dimensionless controlling parameters space: normalized hydrostatic pressure and axial stress. We argue that HSPC production is accompanied by two mechanical instabilities arising in the system due to the plane stress in the DA walls and show how a complex interplay between mechanical forces in the system drives the emerging morphological changes.