RESUMEN

Cross-sections of cell shapes in a tissue monolayer typically resemble a tiling of convex polygons. Yet, examples exist where the polygons are not convex with curved cell-cell interfaces, as seen in the adaxial epidermis. To date, two-dimensional vertex models predicting the structure and mechanics of cell monolayers have been mostly limited to convex polygons. To overcome this limitation, we introduce a framework to study curvy cell-cell interfaces at the subcellular scale within vertex models by using a parametrized curve between vertices that is expanded in a Fourier series and whose coefficients represent additional degrees of freedom. This extension to non-convex polygons allows for cells with the same shape index, or dimensionless perimeter, to be, for example, either elongated or globular with lobes. In the presence of applied, anisotropic stresses, we find that local, subcellular curvature or buckling can be energetically more favourable than larger scale deformations involving groups of cells. Inspired by recent experiments, we also find that local, subcellular curvature at cell-cell interfaces emerges in a group of cells in response to the swelling of additional cells surrounding the group. Our framework, therefore, can account for a wider array of multicellular responses to constraints in the tissue environment.

Asunto(s)

RESUMEN

The epithelial adaptations to mechanical stress are facilitated by molecular and tissue-scale changes that include the strengthening of junctions, cytoskeletal reorganization, and cell-proliferation-mediated changes in tissue rheology. However, the role of cell size in controlling these properties remains underexplored. Our experiments in the zebrafish embryonic epidermis, guided by theoretical estimations, reveal a link between epithelial mechanics and cell size, demonstrating that an increase in cell size compromises the tissue fracture strength and compliance. We show that an increase in E-cadherin levels in the proliferation-deficient epidermis restores epidermal compliance but not the fracture strength, which is largely regulated by Ezrin-an apical membrane-cytoskeleton crosslinker. We show that Ezrin fortifies the epithelium in a cell-size-dependent manner by countering non-muscle myosin-II-mediated contractility. This work uncovers the importance of cell size maintenance in regulating the mechanical properties of the epithelium and fostering protection against future mechanical stresses.

Asunto(s)

RESUMEN

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

Asunto(s)

RESUMEN

Most of normal proliferative epithelia of plants and metazoans are topologically invariant and characterized by similar cell distributions according to the number of cell neighbors (DCNs). Here we study peculiarities of these distributions and explain why the DCN obtained from the location of intercellular boundaries and that based on the Voronoi tessellation with nodes located on cell nuclei may differ from each other. As we demonstrate, special microdomains where four or more intercellular boundaries converge are topologically charged. Using this fact, we deduce a new equation describing the topological balance of the DCNs. The developed theory is applied for a series of microphotographs of non-tumoral epithelial cells of the human cervix (HCerEpiC) to improve the image processing near the edges of microphotographs and reveal the topological invariance of the examined monolayers. Special contact microdomains may be present in epithelia of various natures, however, considering the well-known vertex model of epithelium, we show that such contacts are absent in the usual solid-like state of the model and appear only in the liquid-like cancer state. Also, we discuss a possible biological role of special contacts in context of proliferative epithelium dynamics and tissue morphogenesis.

Asunto(s)

RESUMEN

In this work, we present a mathematical model of cell growth in the pores of a perfusion bioreactor through which a nutrient solution is pumped. We have developed a 2-D vertex model that allows us to reproduce the microscopic dynamics of the microenvironment of cells and describe the occupation of the pore space with cells. In this model, each cell is represented by a polygon; the number of vertices and shapes may change over time. The model includes mitotic cell division and intercalation. We study the impact of two factors on cell growth. On the one hand, we consider a channel of variable cross-section, which models a scaffold with a porosity gradient. On the other hand, a cluster of cells grows under the influence of a nutrient solution flow, which establishes a non-uniform distribution of shear stresses in the pore space. We present the results of numerical simulation of the tissue growth in a wavy channel. The model allows us to obtain complete microscopic information that includes the dynamics of intracellular pressure, the local elastic energy, and the characteristics of cell populations. As we showed, in a functional-graded scaffold, the distribution of the shear stresses in the pore space has a complicated structure, which implies the possibility of controlling the growth zones by varying the pore geometry.

RESUMEN

Cellular systems are known to exhibit some of the fastest movements in biology, but little is known as to how single cells can dissipate this energy rapidly and adapt to such large accelerations without disrupting internal architecture. To address this, we investigate Spirostomum ambiguum-a giant cell (1-4 mm in length) well-known to exhibit ultrafast contractions (50% of body length) within 5 ms with a peak acceleration of 15[Formula: see text]. Utilizing transmitted electron microscopy and confocal imaging, we identify an association of rough endoplasmic reticulum (RER) and vacuoles throughout the cell-forming a contiguous fenestrated membrane architecture that topologically entangles these two organelles. A nearly uniform interorganelle spacing of 60 nm is observed between RER and vacuoles, closely packing the entire cell. Inspired by the entangled organelle structure, we study the mechanical properties of entangled deformable particles using a vertex-based model, with all simulation parameters matching 10 dimensionless numbers to ensure dynamic similarity. We demonstrate how entangled deformable particles respond to external loads by an increased viscosity against squeezing and help preserve spatial relationships. Because this enhanced damping arises from the entanglement of two networks incurring a strain-induced jamming transition at subcritical volume fractions, which is demonstrated through the spatial correlation of velocity direction, we term this phenomenon "topological damping." Our findings suggest a mechanical role of RER-vacuolar meshwork as a metamaterial capable of damping an ultrafast contraction event.

Asunto(s)

RESUMEN

The mesoderm invagination of the Drosophila embryo is known as an archetypal morphogenic process. To explore the roles of the active cellular forces and the regulation of these forces, we developed an integrated vertex model that combines the regulation of morphogen expression with cell movements and tissue mechanics. Our results suggest that a successful furrow formation requires an apical tension gradient, decreased basal tension, and increased lateral tension, which corresponds to apical constriction, basal expansion, and apicobasal shortening respectively. Our model also considers the mechanical feedback which leads to an ectopic twist expression with external compression as observed in experiments. Our model predicts that ectopic invagination could happen if an external compressive gradient is applied.

Asunto(s)

RESUMEN

Heart development begins with the formation of a tube as cardiac progenitors migrate from opposite sides of the embryo. Abnormal cardiac progenitor movements cause congenital heart defects. However, the mechanisms of cell migration during early heart development remain poorly understood. Using quantitative microscopy, we found that in Drosophila embryos, cardiac progenitors (cardioblasts) migrated through a sequence of forward and backward steps. Cardioblast steps were associated with oscillatory non-muscle myosin II waves that induced periodic shape changes and were necessary for timely heart tube formation. Mathematical modeling predicted that forward cardioblast migration required a stiff boundary at the trailing edge. Consistent with this, we found a supracellular actin cable at the trailing edge of the cardioblasts that limited the amplitude of the backward steps, thus biasing the direction of cell movement. Our results indicate that periodic shape changes coupled with a polarized actin cable produce asymmetrical forces that promote cardioblast migration.

Asunto(s)

RESUMEN

Convergence-extension in embryos is controlled by chemical and mechanical signalling. A key cellular process is the exchange of neighbours via T1 transitions. We propose and analyse a model with positive feedback between recruitment of myosin motors and mechanical tension in cell junctions. The model produces active T1 events, which act to elongate the tissue perpendicular to the main direction of tissue stress. Using an idealised tissue patch comprising several active cells embedded in a matrix of passive hexagonal cells, we identified an optimal range of mechanical stresses to trigger an active T1 event. We show that directed stresses also generate tension chains in a realistic patch made entirely of active cells of random shapes and leads to convergence-extension over a range of parameters. Our findings show that active intercalations can generate stress that activates T1 events in neighbouring cells, resulting in tension-dependent tissue reorganisation, in qualitative agreement with experiments on gastrulation in chick embryos.

Asunto(s)

RESUMEN

Biological tissues acquire reproducible shapes during development through dynamic cell behaviors. Most of these behaviors involve the remodeling of cell-cell contacts. During epithelial morphogenesis, contractile actomyosin networks remodel cell-cell contacts by shrinking and extending junctions between lateral cell surfaces. However, actomyosin networks not only generate mechanical stresses but also respond to them, confounding our understanding of how mechanical stresses remodel cell-cell contacts. Here, we develop a two-point optical manipulation method to impose different stress patterns on cell-cell contacts in the early epithelium of the Drosophila embryo. The technique allows us to produce junction extension and shrinkage through different push and pull manipulations at the edges of junctions. We use these observations to expand classical vertex-based models of tissue mechanics, incorporating negative and positive mechanosensitive feedback depending on the type of remodeling. In particular, we show that Myosin-II activity responds to junction strain rate and facilitates full junction shrinkage. Altogether our work provides insight into how stress produces efficient deformation of cell-cell contacts in vivo and identifies unanticipated mechanosensitive features of their remodeling.

Asunto(s)

RESUMEN

Internal organs such as the heart demonstrate apparent left-right (LR) asymmetric morphology and positioning. Cellular chirality and associated LR biased mechanical behavior such as cell migration have been attributed to LR symmetry breaking during embryonic development. Mathematical models have shown that chiral directional migration can be driven by cellular intrinsic torque. Tissue jamming state (i.e., solid-like vs fluid-like state) strongly regulates collective migratory behavior, but how it might affect chiral morphogenesis is still unknown. Here, we develop a cell vertex model to study the role of tissue rigidity or jamming state on chiral morphogenesis of the cells on a patterned ring-shaped tissue, simulating a previously reported experimental setup for measuring cell chirality. We simulate chirality as torsional forces acting on cell vertices. As expected, the cells undergo bidirectional migration at the opposing (inner and outer) boundaries of the ring-shaped tissue. We discover that more fluid-like tissues (unjammed) demonstrate a stronger chiral cell alignment and elongation than more solid-like (jammed) tissues and maintain a bigger difference in migration velocity between opposing tissue boundaries. Finally, we find that fluid-like tissues undergo more cell-neighbor exchange events. This study reveals that chiral torque is sufficient to achieve a biased cellular alignment as seen in vitro. It further sheds light on the mechanical regulation of chiral morphogenesis of tissues and reveals a role of cell density-independent tissue rigidity in this process.

Asunto(s)

RESUMEN

The vertex model is widely used to simulate the mechanical properties of confluent epithelia and other multicellular tissues. This inherently discrete framework allows a Cauchy stress to be attributed to each cell, and its symmetric component has been widely reported, at least for planar monolayers. Here, we consider the stress attributed to the neighbourhood of each tricellular junction, evaluating in particular its leading-order antisymmetric component and the associated couple stresses, which characterise the degree to which individual cells experience (and resist) in-plane bending deformations. We develop discrete potential theory for localised monolayers having disordered internal structure and use this to derive the analogues of Airy and Mindlin stress functions. These scalar potentials typically have broad-banded spectra, highlighting the contributions of small-scale defects and boundary layers to global stress patterns. An affine approximation attributes couple stresses to pressure differences between cells sharing a trijunction, but simulations indicate an additional role for non-affine deformations.

Asunto(s)

RESUMEN

Dorsal closure is a process that occurs during embryogenesis of Drosophila melanogaster. During dorsal closure, the amnioserosa (AS), a one-cell thick epithelial tissue that fills the dorsal opening, shrinks as the lateral epidermis sheets converge and eventually merge. During this process, the aspect ratio of amnioserosa cells increases markedly. The standard 2-dimensional vertex model, which successfully describes tissue sheet mechanics in multiple contexts, would in this case predict that the tissue should fluidize via cell neighbor changes. Surprisingly, however, the amnioserosa remains an elastic solid with no such events. We here present a minimal extension to the vertex model that explains how the amnioserosa can achieve this unexpected behavior. We show that continuous shrink-age of the preferred cell perimeter and cell perimeter polydispersity lead to the retention of the solid state of the amnioserosa. Our model accurately captures measured cell shape and orientation changes and predicts non-monotonic junction tension that we confirm with laser ablation experiments.

RESUMEN

Leaf meristem is a cell proliferative zone present in the lateral organ primordia. In this study, we examined how cell proliferative zones in primordia of planar floral organs and polar auxin transport inhibitor (PATI)-treated leaf organs differ from those of non-treated foliage leaves of Arabidopsis thaliana, with a focus on the accumulation pattern of ANGUSTIFOLIA3 (AN3) protein, a key element for leaf meristem positioning. We found that PATI-induced leaf shape changes were correlated with cell division angle but not with meristem positioning/size or AN3 localisation. In contrast, different shapes between sepals and petals compared with foliage leaves were associated with both altered meristem position, due to altered AN3 expression patterns, and different distributions of cell division angles. A numerical simulation showed that meristem position majorly affected the final shape but biased cell division angles had a minor effect. Taken together, these results suggest that the unique shapes of different lateral organs depend on the position of the meristem in the case of floral organs and cell division angles in the case of leaf organs with different auxin flow.

Asunto(s)

RESUMEN

Remodeling of multicellular architecture is a critical developmental process for shaping the axis of a bilaterally symmetric animal body and involves coordinated cell-cell interactions and cell rearrangement. In arthropods, the early embryonic process that leads to the segmented body axis varies at the cellular and molecular levels depending on the species. Developmental studies using insect and spider model species have provided specific examples of these diversified mechanisms that regulate axis formation and segmentation in arthropod embryos. However, there are few theoretical models for how diversity in the early embryonic process occurred during evolution, in part because of a limited computational infrastructure. We developed a virtual spherical-shaped multicellular platform to reproduce body axis-forming processes. Each virtual cell behaves according to the cell vertex model, with the computational program organized in a hierarchical order from cells and tissues to whole embryos. Using an initial set of two different mechanical states for cell differentiation and global directional signals that are linked to the planar polarity of each cell, the virtual cell assembly exhibited morphogenetic processes similar to those observed in spider embryos. We found that the development of an elongating body axis is achieved through implementation of an interactive cell polarity parameter associated with edge tension at the cell-cell adhesion interface, with no local control of the cell division rate and direction. We also showed that modifying the settings can cause variation in morphogenetic processes. This platform also can embed a gene network that generates waves of gene expression in a virtual dynamic multicellular field. This study provides a computational platform for testing the development and evolution of animal body patterns.

RESUMEN

Pattern formation and morphogenesis of cell populations is essential for successful embryogenesis. Steinberg proposed the differential adhesion hypothesis, and differences in cell-cell adhesion and interfacial tension have proven to be critical for cell sorting. Standard theoretical models such as the vertex model consider not only cell-cell adhesion/tension but also area elasticity of apical cell surfaces and viscous friction forces. However, the potential contributions of the latter two parameters to pattern formation and morphogenesis remain to be determined. In this theoretical study, we analyzed the effect of both area elasticity and the coefficient of friction on pattern formation and morphogenesis. We assumed the presence of two cell populations, one population of which is surrounded by the other. Both populations were placed on the surface of a uniformly expanding environment analogous to growing embryos, in which friction forces are exerted between cell populations and their expanding environment. When the area elasticity or friction coefficient in the cell cluster was increased relative to that of the surrounding cell population, the cell cluster was elongated. In comparison with experimental observations, elongation of the notochord in mice is consistent with the hypothesis based on the difference in area elasticity but not the difference in friction coefficient. Because area elasticity is an index of cellular stiffness, we propose that differential cellular stiffness may contribute to tissue elongation within an expanding environment.

RESUMEN

Several models have been proposed to describe the dynamics of epithelial tissues undergoing morphogenetic changes driven by apical constriction pulses, which differ in where the constriction is applied, either at the perimeter or in the medial regions. To help discriminate between these models, we analyse the impact of where constriction is applied on the final geometry of the active contracted cell, using the two-dimensional vertex model. We find that medial activity, characterized by a reduction in the reference area, generates anisotropic cell shapes, whereas isotropic cell shapes are produced when the reference perimeter is reduced. When plasticity is included, sufficiently slow processes of medial contractile activity, compared with the characteristic times of elasticity and plasticity, cells can achieve less elongated shapes. Similarly, for perimeter activity, the highest level of contraction is achieved. Finally, we apply the model to describe the apical contractile pulses observed within the epithelial enveloping cell layer during the pre-epiboly of the annual killifish Austrolebias nigripinnis. The analysis of the cell shape changes allowed a global fit of all parameters of the vertex model, with the pulses being quantitatively captured using perimeter activity and area plasticity.

Asunto(s)

RESUMEN

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.

Asunto(s)

RESUMEN

The left-right organizer in zebrafish embryos, Kupffer's Vesicle (KV), is a simple organ that undergoes programmed asymmetric cell shape changes that are necessary to establish the left-right axis of the embryo. We use simulations and experiments to investigate whether 3D mechanical drag forces generated by the posteriorly-directed motion of the KV through the tailbud tissue are sufficient to drive such shape changes. We develop a fully 3D vertex-like (Voronoi) model for the tissue architecture, and demonstrate that the tissue can generate drag forces and drive cell shape changes. Furthermore, we find that tailbud tissue presents a shear-thinning, viscoelastic behavior consistent with those observed in published experiments. We then perform live imaging experiments and particle image velocimetry analysis to quantify the precise tissue velocity gradients around KV as a function of developmental time. We observe robust velocity gradients around the KV, indicating that mechanical drag forces must be exerted on the KV by the tailbud tissue. We demonstrate that experimentally observed velocity fields are consistent with the viscoelastic response seen in simulations. This work also suggests that 3D viscoelastic drag forces could be a generic mechanism for cell shape change in other biological processes.

Asunto(s)

RESUMEN

Organogenesis is the phase of embryonic development leading to the formation of fully functional organs. In the case of the thyroid, organogenesis starts from the endoderm and generates a multitude of closely packed independent spherical follicular units surrounded by a dense network of capillaries. Follicular organisation is unique and essential for thyroid function, i.e. thyroid hormone production. Previous in vivo studies showed that, besides their nutritive function, endothelial cells play a central role during thyroid gland morphogenesis. However, the precise mechanisms and biological parameters controlling the transformation of the multi-layered thyroid epithelial primordium into a multitude of single-layered follicles are mostly unknown. Animal studies used to improve understanding of organogenesis are costly and time-consuming, with recognised limitations. Here, we developed and used a 2-D vertex model of thyroid growth, angiogenesis and folliculogenesis, within the open-source Chaste framework. Our in silico model, based on in vivo images, correctly simulates the differential growth and proliferation of central and peripheral epithelial cells, as well as the morphogen-driven migration of endothelial cells, consistently with our experimental data. Our simulations further showed that reduced epithelial cell adhesion was critical to allow endothelial invasion and fission of the multi-layered epithelial mass. Finally, our model also allowed epithelial cell polarisation and follicular lumen formation by endothelial cell abundance and proximity. Our study illustrates how constant discussion between theoretical and experimental approaches can help us to better understand the roles of cellular movement, adhesion and polarisation during thyroid embryonic development. We anticipate that the use of in silico models like the one we describe can push forward the fields of developmental biology and regenerative medicine.

Asunto(s)