Effect of cellular rearrangement time delays on the rheology of vertex models for confluent tissues

Strain softening and hysteresis arising from 3D multicellular dynamics during long-term large deformation

Living tissues exhibit complex mechanical properties, including viscoelastic and elastoplastic responses, that are crucial for regulating cell behaviors and tissue deformations. Despite their significance, the intricate properties of three-dimensional (3D) cell constructs are not well understood and are inadequately implemented in biomaterial engineering. To address this gap, we developed a numerical method to analyze the dynamic properties of cell constructs using a 3D vertex model framework. By focusing on 3D tissues composed of confluent homogeneous cells, we characterized their properties in response to various deformation magnitudes and time scales. Stress relaxation tests revealed that large deformations initially induced relaxation in the shapes of individual cells. This process is amplified by subsequent transient cell rearrangements, homogenizing cell shapes and leading to tissue fluidization. Additionally, dynamic viscoelastic analyses showed that tissues exhibited strain softening and hysteresis during large deformations. Interestingly, this strain softening originates from multicellular structures independent of cell rearrangement, while hysteresis arises from cell rearrangement. Moreover, tissues exhibit elastoplastic responses over the long term, which are well represented by the Ramberg-Osgood model. These findings highlight the characteristic properties of cell constructs emerging from their structures and rearrangements, especially during long-term large deformations. The developed method offers a new approach to uncover the dynamic nature of 3D tissue mechanics and could serve as a technical foundation for exploring tissue mechanics and advancing biomaterial engineering.

Contact area and tissue growth dynamics shape synthetic juxtacrine signaling patterns

Cell-cell communication through direct contact, or juxtacrine signaling, is important in development, disease, and many areas of physiology. Synthetic forms of juxtacrine signaling can be precisely controlled and operate orthogonally to native processes, making them a powerful reductionist tool with which to address fundamental questions in cell-cell communication in vivo. Here, we investigate how cell-cell contact length and tissue growth dynamics affect juxtacrine signal responses through implementing a custom synthetic gene circuit in Drosophila wing imaginal discs alongside mathematical modeling to determine synthetic Notch (synNotch) activation patterns. We find that the area of contact between cells largely determines the extent of synNotch activation, leading to the prediction that the shape of the interface between signal-sending and signal-receiving cells will impact the magnitude of the synNotch response. Notably, synNotch outputs form a graded spatial profile that extends several cell diameters from the signal source, providing evidence that the response to juxtacrine signals can persist in cells as they proliferate away from source cells, or that cells remain able to communicate directly over several cell diameters. Our model suggests that the former mechanism may be sufficient, since it predicts graded outputs without diffusion or long-range cell-cell communication. Overall, we identify that cell-cell contact area together with output synthesis and decay rates likely govern the pattern of synNotch outputs in both space and time during tissue growth, insights that may have broader implications for juxtacrine signaling in general.

Regulation of epithelial cell jamming transition by cytoskeleton and cell-cell interactions

Multicellular systems, such as epithelial cell collectives, undergo transitions similar to those in inert physical systems like sand piles and foams. To remodel or maintain tissue organization during development or disease, these collectives transition between fluid-like and solid-like states, undergoing jamming or unjamming transitions. While these transitions share principles with physical systems, understanding their regulation and implications in cell biology is challenging. Although cell jamming and unjamming follow physics principles described by the jamming diagram, they are fundamentally biological processes. In this review, we explore how cellular processes and interactions regulate jamming and unjamming transitions. We begin with an overview of how these transitions control tissue remodeling in epithelial model systems and describe recent findings of the physical principles governing tissue solidification and fluidization. We then explore the mechanistic pathways that modulate the jamming phase diagram axes, focusing on the regulation of cell fluctuations and geometric compatibility. Drawing upon seminal works in cell biology, we discuss the roles of cytoskeleton and cell-cell adhesion in controlling cell motility and geometry. This comprehensive view illustrates the molecular control of cell jamming and unjamming, crucial for tissue remodeling in various biological contexts.

Interplay of geometry and mechanics in epithelial wound healing

Wound healing is a complex biological process critical for maintaining an organism's structural integrity and tissue repair following an infection or injury. Recent studies have unveiled the mechanisms involving the coordination of biochemical and mechanical responses in the tissue in wound healing. In this article, we focus on the healing property of an epithelial tissue as a material while the effects of biological mechanisms such as cell proliferation, tissue intercalation, cellular migration, cell crawling, and filopodia protrusion is minimal. We present a mathematical framework that predicts the fate of a wounded tissue based on the wound's geometrical features and the tissue's mechanical properties. Precisely, adapting the vertex model of tissue mechanics, we predict whether a wound of a specific size in an epithelial monolayer characterized by certain levels of actomyosin contractility and cell-cell adhesion will heal (i.e., close), shrink in size, or rupture the tissue further. Moreover, we show how tissue-mediated mechanisms such as purse-string tension at the wound boundary facilitate wound healing. Finally, we validate the predictions of our model by designing an experimental setup that enables us to create wounds of specific sizes in kidney epithelial cells (MDCK) monolayers. Altogether, this work sets up a basis for interpreting the interplay of mechanical and geometrical features of a tissue in the process of wound healing.

Contact area and tissue growth dynamics shape synthetic juxtacrine signaling patterns

Cell-cell communication through direct contact, or juxtacrine signaling, is important in development, disease, and many areas of physiology. Synthetic forms of juxtacrine signaling can be precisely controlled and operate orthogonally to native processes, making them a powerful reductionist tool with which to address fundamental questions in cell-cell communication in vivo. Here we investigate how cell-cell contact length and tissue growth dynamics affect juxtacrine signal responses through implementing a custom synthetic gene circuit in Drosophila wing imaginal discs alongside mathematical modeling to determine synthetic Notch (synNotch) activation patterns. We find that the area of contact between cells largely determines the extent of synNotch activation, leading to the prediction that the shape of the interface between signal-sending and signal-receiving cells will impact the magnitude of the synNotch response. Notably, synNotch outputs form a graded spatial profile that extends several cell diameters from the signal source, providing evidence that the response to juxtacrine signals can persist in cells as they proliferate away from source cells, or that cells remain able to communicate directly over several cell diameters. Our model suggests the former mechanism may be sufficient, since it predicts graded outputs without diffusion or long-range cell-cell communication. Overall, we identify that cell-cell contact area together with output synthesis and decay rates likely govern the pattern of synNotch outputs in both space and time during tissue growth, insights that may have broader implications for juxtacrine signaling in general.

General, open-source vertex modeling in biological applications using Tissue Forge

Vertex models are a widespread approach for describing the biophysics and behaviors of multicellular systems, especially of epithelial tissues. Vertex models describe a wide variety of developmental scenarios and behaviors like cell rearrangement and tissue folding. Often, these models are implemented as single-use or closed-source software, which inhibits reproducibility and decreases accessibility for researchers with limited proficiency in software development and numerical methods. We developed a physics-based vertex model methodology in Tissue Forge, an open-source, particle-based modeling and simulation environment. Our methodology describes the properties and processes of vertex model objects on the basis of vertices, which allows integration of vertex modeling with the particle-based formalism of Tissue Forge, enabling an environment for developing mixed-method models of multicellular systems. Our methodology in Tissue Forge inherits all features provided by Tissue Forge, delivering open-source, extensible vertex modeling with interactive simulation, real-time simulation visualization and model sharing in the C, C++ and Python programming languages and a Jupyter Notebook. Demonstrations show a vertex model of cell sorting and a mixed-method model of cell migration combining vertex- and particle-based models. Our methodology provides accessible vertex modeling for a broad range of scientific disciplines, and we welcome community-developed contributions to our open-source software implementation.

Tension Remodeling Regulates Topological Transitions in Epithelial Tissues

Cell neighbor exchanges play a critical role in regulating tissue fluidity during epithelial morphogenesis and repair. In vivo, these neighbor exchanges are often hindered by the formation of transiently stable fourfold vertices, which can develop into complex multicellular rosettes where five or more cell junctions meet. Despite their importance, the mechanical origins of multicellular rosettes have remained elusive, and current cellular models lack the ability to explain their formation and maintenance. Here we present a dynamic vertex model of epithelial tissues with strain-dependent tension remodeling and mechanical memory dissipation. We show that an increase in cell junction tension upon contraction and reduction in tension upon extension can stabilize higher-order vertices, temporarily stalling cell rearrangements. On the other hand, inducing mechanical memory dissipation via relaxation of junction strain and stress promotes the resolution of higher-order vertices, facilitating cell neighbor exchanges. We demonstrate that by tuning the rates of tension remodeling and mechanical memory dissipation, we can control topological transitions and tissue material properties, recapitulating complex cellular topologies seen in developing organisms.

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.

General, Open-Source Vertex Modeling in Biological Applications Using Tissue Forge

Vertex models are a widespread approach for describing the biophysics and behaviors of multicellular systems, especially of epithelial tissues. Vertex models describe a wide variety of developmental scenarios and behaviors like cell rearrangement and tissue folding. Often, these models are implemented as single-use or closed-source software, which inhibits reproducibility and decreases accessibility for researchers with limited proficiency in software development and numerical methods. We developed a physics-based vertex model methodology in Tissue Forge, an open-source, particle-based modeling and simulation environment. Our methodology describes the properties and processes of vertex model objects on the basis of vertices, which allows integration of vertex modeling with the particle-based formalism of Tissue Forge, enabling an environment for developing mixed-method models of multicellular systems. Our methodology in Tissue Forge inherits all features provided by Tissue Forge, delivering opensource, extensible vertex modeling with interactive simulation, real-time simulation visualization and model sharing in the C , C + + and Python programming languages and a Jupyter Notebook. Demonstrations show a vertex model of cell sorting and a mixed-method model of cell migration combining vertex- and particle-based models. Our methodology provides accessible vertex modeling for a broad range of scientific disciplines, and we welcome community-developed contributions to our open-source software implementation.

Cell cycle dynamics control fluidity of the developing mouse neuroepithelium

As developing tissues grow in size and undergo morphogenetic changes, their material properties may be altered. Such changes result from tension dynamics at cell contacts or cellular jamming. Yet, in many cases, the cellular mechanisms controlling the physical state of growing tissues are unclear. We found that at early developmental stages, the epithelium in the developing mouse spinal cord maintains both high junctional tension and high fluidity. This is achieved via a mechanism in which interkinetic nuclear movements generate cell area dynamics that drive extensive cell rearrangements. Over time, the cell proliferation rate declines, effectively solidifying the tissue. Thus, unlike well-studied jamming transitions, the solidification uncovered here resembles a glass transition that depends on the dynamical stresses generated by proliferation and differentiation. Our finding that the fluidity of developing epithelia is linked to interkinetic nuclear movements and the dynamics of growth is likely to be relevant to multiple developing tissues.

Powering morphogenesis: multiscale challenges at the interface of cell adhesion and the cytoskeleton

Among the defining features of the animal kingdom is the ability of cells to change shape and move. This underlies embryonic and postembryonic development, tissue homeostasis, regeneration, and wound healing. Cell shape change and motility require linkage of the cell's force-generating machinery to the plasma membrane at cell-cell and cell-extracellular matrix junctions. Connections of the actomyosin cytoskeleton to cell-cell adherens junctions need to be both resilient and dynamic, preventing tissue disruption during the dramatic events of embryonic morphogenesis. In the past decade, new insights radically altered the earlier simple paradigm that suggested simple linear linkage via the cadherin-catenin complex as the molecular mechanism of junction-cytoskeleton interaction. In this Perspective we provide a brief overview of our current state of knowledge and then focus on selected examples highlighting what we view as the major unanswered questions in our field and the approaches that offer exciting new insights at multiple scales from atomic structure to tissue mechanics.

All-in-one rheometry and nonlinear rheology of multicellular aggregates

Tissues are generally subjected to external stresses, a potential stimulus for their differentiation or remodeling. While single-cell rheology has been extensively studied leading to controversial results about nonlinear response, mechanical tissue behavior under external stress is still poorly understood, in particular, the way individual cell properties translate at the tissue level. Herein, using magnetic cells we were able to form perfectly monitored cellular aggregates (magnetic molding) and to deform them under controlled applied stresses over a wide range of timescales and amplitudes (magnetic rheometer). We explore the rheology of these minimal tissue models using both standard assays (creep and oscillatory response) as well as an innovative broad spectrum solicitation coupled with inference analysis thus being able to determine in a single experiment the best rheological model. We find that multicellular aggregates exhibit a power-law response with nonlinearities leading to tissue stiffening at high stress. Moreover, we reveal the contribution of intracellular (actin network) and intercellular components (cell-cell adhesions) in this aggregate rheology.

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.

Non-monotonic fluidization generated by fluctuating edge tensions in confluent tissues

In development and homeostasis, multi-cellular systems exhibit spatial and temporal heterogeneity in their biochemical and mechanical properties. Nevertheless, it remains unclear how spatiotemporally heterogeneous forces affect the dynamical and mechanical properties of confluent tissue. To address this question, we study the dynamical behavior of the two-dimensional cellular vertex model for epithelial monolayers in the presence of fluctuating cell-cell interfacial tensions, which is a biologically relevant source of mechanical spatiotemporal heterogeneity. In particular, we investigate the effects of the amplitude and persistence time of fluctuating tension on the tissue dynamics. We unexpectedly find that the long-time diffusion constant describing cell rearrangements depends non-monotonically on the persistence time, while it increases monotonically as the amplitude increases. Our analysis indicates that at low and intermediate persistence times tension fluctuations drive motion of vertices and promote cell rearrangements, while at the highest persistence times the tension in the network evolves so slowly that rearrangements become rare.

Adhesion-regulated junction slippage controls cell intercalation dynamics in an Apposed-Cortex Adhesion Model

Cell intercalation is a key cell behaviour of morphogenesis and wound healing, where local cell neighbour exchanges can cause dramatic tissue deformations such as body axis extension. Substantial experimental work has identified the key molecular players facilitating intercalation, but there remains a lack of consensus and understanding of their physical roles. Existing biophysical models that represent cell-cell contacts with single edges cannot study cell neighbour exchange as a continuous process, where neighbouring cell cortices must uncouple. Here, we develop an Apposed-Cortex Adhesion Model (ACAM) to understand active cell intercalation behaviours in the context of a 2D epithelial tissue. The junctional actomyosin cortex of every cell is modelled as a continuous viscoelastic rope-loop, explicitly representing cortices facing each other at bicellular junctions and the adhesion molecules that couple them. The model parameters relate directly to the properties of the key subcellular players that drive dynamics, providing a multi-scale understanding of cell behaviours. We show that active cell neighbour exchanges can be driven by purely junctional mechanisms. Active contractility and cortical turnover in a single bicellular junction are sufficient to shrink and remove a junction. Next, a new, orthogonal junction extends passively. The ACAM reveals how the turnover of adhesion molecules regulates tension transmission and junction deformation rates by controlling slippage between apposed cell cortices. The model additionally predicts that rosettes, which form when a vertex becomes common to many cells, are more likely to occur in actively intercalating tissues with strong friction from adhesion molecules.

Sculpting with stem cells: how models of embryo development take shape

During embryogenesis, organisms acquire their shape given boundary conditions that impose geometrical, mechanical and biochemical constraints. A detailed integrative understanding how these morphogenetic information modules pattern and shape the mammalian embryo is still lacking, mostly owing to the inaccessibility of the embryo in vivo for direct observation and manipulation. These impediments are circumvented by the developmental engineering of embryo-like structures (stembryos) from pluripotent stem cells that are easy to access, track, manipulate and scale. Here, we explain how unlocking distinct levels of embryo-like architecture through controlled modulations of the cellular environment enables the identification of minimal sets of mechanical and biochemical inputs necessary to pattern and shape the mammalian embryo. We detail how this can be complemented with precise measurements and manipulations of tissue biochemistry, mechanics and geometry across spatial and temporal scales to provide insights into the mechanochemical feedback loops governing embryo morphogenesis. Finally, we discuss how, even in the absence of active manipulations, stembryos display intrinsic phenotypic variability that can be leveraged to define the constraints that ensure reproducible morphogenesis in vivo.