Generating active T1 transitions through mechanochemical feedback
Abstract
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.
Editor's evaluation
This theoretical investigation provides important findings on the role of active mechanical feedback on tissue remodelling. The authors present convincing evidence that mechanically enforced myosin recruitment at cell-cell junctions can lead to tissue expansion in the direction perpendicular to an externally applied uniaxial mechanical stress. The relevance of the proposed mechanism for convergence–extension systems requires more investigation through comparison with experimental data.
https://doi.org/10.7554/eLife.79862.sa0Introduction
Embryonic development involves complex tissue dynamics, including rearrangements and shape changes of the cells. This is particularly evident during gastrulation where the presumptive ectoderm, mesoderm, and endoderm take up their correct positions in the embryo (Wolpert et al., 2015). Key cellular processes that underlie tissue formation and morphogenesis during gastrulation are cell division, differentiation, and cell movement. Directed cell intercalation is a major mechanism driving large-scale tissue shape changes both in epithelial and mesenchymal tissues (Huebner and Wallingford, 2018). The narrowing and lengthening of epithelial tissues resulting from such intercalations, known as convergent extension (Keller et al., 2000), underlie germband extension in Drosophila (Bertet et al., 2004; Blankenship et al., 2006), as well as primitive streak formation in the chick embryo (Voiculescu et al., 2007; Rozbicki et al., 2015). In the latter, cell intercalations facilitate coordinated movements of hundreds of thousands of cells in two counter-rotating millimetre-scale cell flows that drive the formation of the primitive streak at the site where the flows meet (Rozbicki et al., 2015; Saadaoui et al., 2020; Serrano Nájera and Weijer, 2020). Unlike cell migration (Alert and Trepat, 2020), which typically involves a significant contribution from crawling against a substrate such as the extracellular matrix, during intercalation, cells pull against each other in order to exchange their neighbours (Huebner and Wallingford, 2018). This is a complex, active process that requires a carefully coordinated shrinking and subsequent expansion of cell–cell interfaces, known as the T1 transition (Kong et al., 2017).
The widely studied morphological process of germband extension in Drosophila involves directed cell intercalations in the ventral ectoderm, where during the associated T1 events dorsal–ventral (DV)-oriented junctions shrink and new junctions are generated in anterior–posterior (AP) direction. The contraction of DV-oriented junctions has been shown to correlate with increased accumulation of apical myosin II in these junctions that can form supercellular cables in aligned junctions (Bertet et al., 2004; Blankenship et al., 2006). The extension of new junctions has been associated with the activity of medial myosin (Rauzi et al., 2010; Collinet et al., 2015). Laser ablation experiments in the ventral ectoderm have shown that the myosin-rich DV-oriented junctions are under higher tension than AP junctions (Rauzi et al., 2008; Fernandez-Gonzalez et al., 2009; Collinet et al., 2015), and aspiration and optical tweezing and optogenetic experiments have shown that myosin can be recruited to junctions in response to increased tension in these junctions, demonstrating the existence of mechanical feedback (Fernandez-Gonzalez et al., 2009; Clément et al., 2017; Gustafson et al., 2022). This is in agreement with observations in the Drosophila wing disk when it has recently been shown that myosin accumulates on apical junctions in response to mechanical stretching of the disk (Duda et al., 2019). It has recently been shown that the distribution of apical myosin can accurately predict the observed tissue flow patterns during germband extension (Streichan et al., 2018). However, an as yet unresolved question is by which mechanism the anisotropic distribution of myosin cables is initially generated. It is thought to depend on the family of Toll receptors under the control of pair rule genes (Paré et al., 2014) and their interactions with an adhesion G protein0coupled receptor (Lavalou et al., 2021) that could generate asymmetries in cells and can possibly signal to Rho-kinase and myosin; however, the precise molecular details remain to be resolved. Recently, a strong correlation between the DV junctional strain rate gradient and the junctional myosin recruitment rate gradient, both high at the ventral side, has been observed (Gustafson et al., 2022). This has led to the renewed suggestion that the myosin anisotropy may arise in response to extrinsic forces, such as those generated by the mesoderm invagination of the ventral furrow (Butler et al., 2009), a process that starts somewhat before germband extension and in combination with other extrinsic events such as posterior hindgut invagination, and other geometric constraints could drive germband extension (Collinet et al., 2015; Gehrels et al., 2023). In this scenario, the AP and DV patterning system could be involved in setting the level of mechanical feedback.
Experiments in the chick embryo showed that directed intercalations of mesendoderm precursors located in a sickle-shaped region in the posterior of the epiblast drive the tissue flows underlying the formation of the primitive streak (Voiculescu et al., 2007; Rozbicki et al., 2015). This sickle-shaped mesendoderm precursor region contracts along its long axis towards the AP midline of the embryo and extends along this midline in anterior direction to form the primitive streak. Measurements of the directions of intercalations show that they are aligned along the long axis of the mesendoderm sickle in the direction of the contraction of the sickle and correlate well with the direction and magnitude of the anisotropic strain rate component (Rozbicki et al., 2015; Chuai et al., 2023). These directed intercalations are mediated by super cellular myosin cables in aligned junctions in the direction of intercalation (Rozbicki et al., 2015; Saadaoui et al., 2020; Serrano Nájera and Weijer, 2020). The onset of tissue motion starts near the central midline of the sickle and then rapidly extends outwards to more lateral regions, suggesting the existence of an outward-propagating signal (Rozbicki et al., 2015). Furthermore, the intercalating cells go through a characteristic elongation in the direction of intercalation around the time of the onset of motion, which disappears when the epiblast tissue flows pickup speed. These observations, coupled with the fact that the chick embryo epiblast contains more than 60,000 cells at the onset of gastrulation requiring coordination of cell intercalation over large distances, led us to suggest that long-range mechanical signals coordinate the intercalations in the large scale (Rozbicki et al., 2015; Serrano Nájera and Weijer, 2020). More specifically, Rozbicki et al., 2015; Serrano Nájera and Weijer, 2020 proposed that local contraction of a junction would, through an increase in tension in aligned junctions of neighbouring cells, activate a mechanical feedback process in those junctions, in turn, resulting in their contraction. This process could explain the formation of the observed myosin cables in response to tension and result in coordinated and directed intercalations. So far, it has not yet been possible to test tension-dependent myosin recruitment directly in chick embryos since there has not been a live indicator of myosin activity. Currently, it is only possible to observe active myosin using a phospho-myosin light chain antibody in fixed embryos. The orientation and alignment of the myosin cables, however, correlate well with the anisotropic component of the strain rate tensor (Rozbicki et al., 2015; Chuai et al., 2023), making it likely that a tension-dependent recruitment of myosin occurs in chick embryos.
To develop a cell-level model of the convergence–extension process, it is necessary to understand how externally applied and internally generated mechanical stresses couple to the signalling pathways that regulate the cell’s mechanical response. One, therefore, needs to understand the feedback between mechanical stress anisotropy and the anisotropy of the distribution of force-generating molecular motors in the cell, that is, how it emerges and is propagated and coordinated over large distances. The aim of this study is to formulate and analyse a model for cell intercalations that includes explicit mechanochemical coupling and does not require initial chemical prepatterning. The initial symmetry breaking is driven by anisotropic mechanical stresses rather than anisotropic distribution of signalling molecules. The focus of this work is on the mechanism of the T1 transition that occurs perpendicular to the direction of the maximum principal mechanical stress. We refer to such T1 transition that generate stress as active, which is different from T1 transitions that relieve stress by intercalating in a direction perpendicular to stresses generated by surrounding tissues. In the present model, this stress is assumed to be anisotropic and externally applied, while in an embryo it is produced by the tissue surrounding the region of interest, for example, the sickle-shaped region in the posterior of the chick embryo that develops into the primitive streak (Serrano Nájera and Weijer, 2020). Unlike passive T1 events that are local plastic rearrangements that relieve the applied stresses as, for example, in foams (Weaire and Hutzler, 2001), active T1 transitions studied here require the cell to induce junction contractions via self-amplifying generation of tension. The key ingredient of the model is, therefore, a feedback mechanism between the kinetics of the force-producing molecules, here assumed to be myosin, and mechanical tension in cell junctions. One of the simplest forms in which this feedback could be implemented is through the formation of well-documented catch bounds between myosin heads and actin filaments, which is highly relevant for proper muscle function (Veigel et al., 2003). In this catch-bond mechanism, the dissociation rate of myosin from actin filaments is a decaying exponential function of tension, where increased tension results in dissociation rate of myosin from the actin-cytoskeleton to be a simple exponential function of tension where increased tension results in a lower dissociation rate of myosin. Assuming that the association rate is not tension sensitive, this process will result in a net tension-dependent myosin accumulation.
Here we construct a model that generically provides a mechanism for active T1 events that underlie convergent extension flows such as those observed during primitive streak formation in the chick embryo (Rozbicki et al., 2015). Our analysis indicates that the viscoelastic nature of the cell–cell junctions is essential for an active T1 event, in agreement with studies on ratcheting during junction contractions (Clément et al., 2017; Staddon et al., 2019). In addition, for the active T1 transition to be possible, there must be a separation of elastic (), viscoelastic remodelling (), and motor turnover timescales (), with .
We first analyse the mechanochemical feedback in the case of a single junction, which describes the key ingredients of the proposed mechanism but avoids complications associated with cell rearrangements. The analysis then proceeds to the two-dimensional case, implemented as an extension of the vertex model (Farhadifar et al., 2007; Fletcher et al., 2014), where cell rearrangement is not only possible but leads to shape changes at the tissue scale. We find active T1 transitions and convergence–extension flows over a broad region of externally applied stresses and relaxation timescales, confirming that the proposed mechanism is robust.
Results
Single-junction model
To understand the mechanism that couples the kinetics of myosin motors to the local mechanical tension and leads to the activation of contractility in cell–cell junctions, we first analyse a model of a single junction. The single-junction model thus provides insight into the conditions under which the junction length can contract to zero and trigger a T1 transition. In this model, for simplicity, the junction is assumed to be surrounded by a tissue which provides an elastic, tension-generating background against which it actively contracts. Guided by experiments, there are three key ingredients of the single-junction model. (1) The junction is viscoelastic, as established by pull-release optical tweezer experiments on cell–cell junctions in Drosophila and chick embryos (Clément et al., 2017; Ferro et al., 2020). This means that the junction is able to remove imposed tension by remodelling itself. (2) The junction can generate tension via the action of myosin motor minifilaments that slide actin filaments against each other. (3) There is self-amplifying feedback due to the exponentially decreasing unbinding rate of myosin motors with tension in the junction (Veigel et al., 2003; Kovács et al., 2007; Figure 1A and Figure 1—figure supplement 1A and B). Furthermore, epithelial tissues, and early-stage embryos, in particular, are under mechanical tension as revealed by tissue and cell junction cutting and tweezing experiments (Clément et al., 2017; Ferro et al., 2020). This tension generates an elastic background against which the junction contracts and expands. We refer to that background elasticity as the elastic barrier since it is assumed that it acts to prevent junction remodelling. In the full model, it additionally captures the yield stress of the underlying material that inhibits T1 transitions (Bi et al., 2015). Details of the model are discussed in ‘Materials and methods’, with parameter values given in Table 1.

An active junction.
(A) An external pulling force of magnitude induces tension in a cell–cell junction of length , which consists of passive viscoelastic and active components. The passive component consists of a Maxwell element with stiffness and viscous relaxation time and a harmonic spring of stiffness and rest length connected to it in parallel. The active component is due to myosin motors (green and blue dots) with concentration that act to contract cortical actin filaments (red lines), exerting a force of magnitude . Myosin motors bind to the actin cortex with association rate and unbind with a tension-dependent dissociation rate . (B) Heatmap plot of the contraction force . For , the junction contraction rate is , where is the friction coefficient with the surrounding medium. The mechanochemical feedback loop is contractile in the top-right quadrant where , , and . Negative values of correspond to an extending junction. (C) Junction length vs. time for , , (black dot in B) for increasing values of the elastic barrier . An active corresponds to reaching . Increasing slows down contractions, until, for , the equilibrium length and no T1 is possible. Inset: myosin dynamics for the same set of junctions; the horizontal dashed line indicates , , , and . Length is measured in units of , time in units of , and force in units of .
Values of the parameters in the single-junction model.
Units: length (), time (), force ().
Base | ||
---|---|---|
Parameter | Description | |
Spring constant | ||
Barrier rest length | ||
Friction with substrate | ||
Model | ||
Parameter | Description | Value range |
Barrier spring constant | ||
Applied external tension | ||
Myosin activity | ||
Viscoelastic time | ||
Myosin time | ||
Myosin reference level | 0.5 | |
Threshold tension | ||
Slope of vs. at | ||
Tension-independent myosin dissociation | 1 |
In order to describe the three key ingredients that characterise a cell–cell junction, while taking into account the effects of the elastic background, we adopt a minimal description, that is, we assume that the junction consists of four elements connected in parallel: (1) a Maxwell element with stiffness and viscous relaxation timescale , which models the viscoelastic character of the junction; (2) an elastic spring with spring constant and rest length , which represents the elastic background, that is, the elastic barrier; (3) an active element that models the contribution of the cytoskeleton by generating active tension , where is the activity, is the ratio of the number of myosin motors bound to the junction and the maximum possible number of bound motors, and m0 is the reference value of ; and (4) a dashpot with dissipation rate , which models dissipation with the surrounding medium. The first two elements form a standard linear solid (SLS) element (Figure 1—figure supplement 1A). The presence of m0 in the active element is necessary to account for the possibility that active contractions of the surrounding tissue are stronger than those in the junction, which would result in it expanding. Furthermore, since and m0 measure the relative numbers of bound motors, the expression for the active force does not include the junction length (for further discussion, see ‘Materials and methods’). Under these assumptions, the dynamics of a junction with length and the rest length l0 is
where is the junction tension, and is external tension, which generates stresses in the junction. The feedback loop between the concentration of bound myosin motors and the mechanical tension is captured by the equation for the myosin dynamics, which incorporates a tension-independent myosin binding rate , and an unbinding rate that decreases with tension as . The third equation in Equation 1, therefore, describes a catch-bond-type mechanism (Dembo et al., 1988; Veigel et al., 2003; Thomas et al., 2008; Prezhdo and Pereverzev, 2009) for myosin kinetics. We remark that it is possible to construct different models for tension-dependent myosin kinetics, for example, with the motor binding rate being tension-dependent. This is, however, not expected to lead to qualitative differences. The key for the mechanism of the active T1 transition studied here is that there is increased contractile activity in response to external mechanical tension, without the need to specify its precise molecular origin. The reason, however, to consider tension-dependent unbinding is because the nonmuscle myosin II, which is the primary molecular motor that drives contractions in early embryonic tissues, is known to have actin association rates that are tension-independent, but exhibits tension-dependent dissociation (Kovács et al., 2007). At steady state , and the equilibrium myosin is a sigmoid function of tension. , therefore, sets the threshold that separates low and high levels of attached myosin motors. and k0 are constants with choices of their values discussed in ‘Materials and methods’.
The first two equations in (1) can be combined as , where . The intersection of nullclines and defines the fixed points of the dynamics, (Figure 1—figure supplement 1C). Experiments of Clément et al., 2017 showed that prolonged pulling forced the junction to remodel and retain the elongated shape. Therefore, the relevant regime consistent with observations in real tissues is where the viscoelastic remodelling and myosin association timescales are longer than the elastic relaxation timescale, that is, for . For , there is a unique stable fixed point that determines the long-time dynamics of the junction and the length of the junction continues to change at a constant rate (Figure 1C, the curve). For a junction with no external load, therefore, a fixed point with corresponds to steady contraction, while a fixed point with corresponds to expansion.
In the presence of a finite elastic barrier of height , the junction behaves as an elastic solid in the long-time limit, and it is in mechanical equilibrium with . This corresponds to a single steady-state solution of Equations 1 at , indicated by the fixed point GB (Figure 1—figure supplement 1C). The corresponding steady-state value of myosin,
is independent of , reflecting the fact that in mechanical equilibrium external tension is balanced by the tension generated by myosin motors. The condition for a T1 transition to occur is, therefore, that the active tension due to myosin motors is sufficiently strong to shrink the junction to . The mechanical equilibrium condition at that point allows one to compute the magnitude of the contraction force the junction generates, or equivalently, the maximum barrier height, that a contracting junction can overcome. Figure 1B shows isolines of , where positive values of correspond to junctions that can contract down to a T1 in the presence of a load, while negative values of correspond to junctions that cannot. Above a threshold in and , the junction is able to gradually generate sufficiently large contraction forces required to overcome the elastic barrier and shrink down. Conversely, for the junction expands, which is the appropriate regime for elongation after a T1. Therefore, there is positive feedback between mechanical tension and activity, which results from the assumption that the myosin association rate is independent of tension while the dissociation rate decays exponentially with it. The isoline (Figure 1B, thick black curve), separates the contracting and the expanding regimes, and it corresponds to a critical threshold for a T1 transition,
where is given by Equation 2. Figure 1C shows the junction dynamics in the presence of barriers of different heights for , with (Figure 1B, black dot). The junction shrinks to a point for , while for larger values of the contraction stops at finite . The initial elongation of the junction (Figure 1C) is due to our choice of the initial value of and reflects the fact that it takes for the active contractile machinery to kick in.
We conclude by emphasising the contractile response to applied tension of the single-junction model that will be at the heart of convergence–extension mechanism discussed below. That is, applying small external forces will lead to an expanding junction. Increasing the force, however, leads to the junction shrinking due to the increase of bound myosin motors. Once initiated, the process of activation can continue spontaneously. To showcase this, we applied a pulling force corresponding to (neither contracting nor expanding) to a chain of active junctions connected in series (Figure 1—figure supplement 2A). An initial myosin pulse applied to the central junction causes it to contract. The contraction of the central junction activates contractions of the neighbouring junctions, leading to shrinking of the entire chain (Figure 1—figure supplement 2B). Finally, the contraction rate strongly depends on the timescale of viscoelastic relaxation in the junction. Both effects should be measurable experimentally.
Vertex model with active junctions
The single-junction model serves as a building block for a model of the entire epithelial tissue. A natural way to proceed is to extend an existing model for tissue mechanics. Setting aside apicobasal polarity, which affects cell intercalations in real tissues (Huebner and Wallingford, 2018), it can be assumed that the mechanical properties of the epithelial tissue arise chiefly from the apical junction cortex, an approximation that is able to qualitatively capture many aspects of tissue mechanics (Fletcher et al., 2014; Murisic et al., 2015).
We, therefore, model the mechanical response of the tissue with the vertex model (Farhadifar et al., 2007; Fletcher et al., 2014). The appeal of the vertex model is that it is able to capture both fluid and solid behaviours of epithelial tissues, including some aspects of the viscoelastic rheology (Tong et al., 2021; Tong et al., 2022). In the vertex model, the transition between the fluid state where cells can easily intercalate and tissue scale flows is possible, and the solid phase where intercalations are arrested is controlled by a single dimensionless geometric parameter, (Staple et al., 2010; Bi et al., 2015; Park et al., 2015; Bi et al., 2016).
Here we extend the vertex model to include activity via the mechanochemical coupling introduced for the single junction and define the vertex model with active junctions. Details of the model are given in ‘Materials and methods’. Each junction, shared by two cells denoted as 1 and 2, is augmented by two active elements that model active contractions by the actomyosin cortex on either side. The tension in the junction is then
where is the passive contribution from the standard vertex model energy given in Equation 11. The junction rest length is, however, not constant but subject to viscoelastic relaxation with . It is important to note that depends only on the cell perimeter and not on junction length and, therefore, does not allow for anisotropic mechanical response to an applied anisotropic tension. This term is however responsible for the yield stress of the vertex model when attempting T1 transitions (Bi et al., 2015) and corresponds to the barrier term in the single-junction model. The additional spring term of the Maxwell element in Equation 4 is, therefore, crucial for generating anisotropic tension. The myosin dynamics of each active element is coupled to a conserved myosin pool of size for each cell that determines the association rate of myosin motors to the junctions. As in the single-junction model, the myosin dissociation rate is modulated by mechanical tension, leading to the following equation for myosin kinetics:
Here, is the total amount of activated myosin bound to the junctions of cell , and we have included a noise component with zero mean and variance to model stochastic binding and unbinding of myosin. The overdamped dynamics of vertices is determined by force balance between friction, elastic forces due to deformations of the passive vertex model, and active forces due to pairs of active elements acting along the junctions connected to the vertex, that is. , where is the position of vertex , is the energy of the passive vertex model (Equation 9), and the active term derives from the active elements introduced in Equation 4.
Single active T1 transition in a hexagonal patch
We begin by discussing a single active T1 event (Figure 2 and Figure 3). The unit of time is set by the elastic timescale , the unit of length by the side of a regular hexagon , and the unit of force by , where is the perimeter modulus of the passive vertex model introduced in Equation 9. Values of the parameters used in simulations are listed in Table 2.

An active T1 transition event.
(A) Top panel: the mechanical anisotropy in the initial state is produced by applying pulling forces (green arrows) in the horizontal direction to the left and right boundaries. Bottom panel: the final state after the active T1 shows a clear convergence–extension deformation (red arrows). Cells are coloured by type: passive (light grey), buffer (medium grey), and active (dark grey). Junctions are coloured by junctional myosin. (B) Time sequence of the active T1 transition measured from the moment activity and viscoelasticity were switched on. Cells in the top row are coloured by type and junctions are coloured by tension. Cells in the bottom row are coloured by activated myosin , and junctions are coloured by myosin. Parameters: , , (active), (buffer), (passive), , , , , , , , , with () cells in the horizontal (vertical) direction. Units: length (), time (), force ().

Junction dynamics during the active T1 transition shown in Figure 2.
(A) Definition of central (red/blue – junction that disappears/appears), inner (orange), and outer (blue-green) shoulder junctions through the T1 transition. (B) Central junction: myosin, (green; two curves for myosin on two sides of the junction), tension, (yellow), junction length, (black), and rest length l0 (purple) vs. time. The vertical line indicates the T1 transition, at which point junctional myosin is redistributed according to the rules outlined in Figure 2—figure supplement 1B. (C) Same as in panel (A) but averaged over four shoulder junctions, with variance indicated as shade. (D) Polar histogram of the orientation of the first T1 event measured with respect to the pulling direction, from simulations. Blue (red) indicates appearing (disappearing) junctions. Parameters are and and as in Figure 2.
Values of the parameters used in the vertex model with active junctions.
Units: length (), time (), force ().
Base | ||
---|---|---|
Parameter | Description | |
Hexagonal cell edge length | ||
Perimeter modulus | ||
Spring constant | ||
Friction with substrate | ||
Model | ||
Parameter | Description | Value range |
Area modulus | ||
Target cell area | ||
Target cell perimeter | ||
Pulling force | ||
Myosin activity | ||
Viscoelastic time | ||
Myosin time | ||
Threshold tension | ||
Slope of vs. at | ||
Myosin reference level | 0.5 | |
Total cell myosin | 6 | |
Variance of myosin fluctuations | 1 | |
Tension-independent myosin dissociation | 0.1 |
To understand the dynamics of a single active T1, we first studied a regular lattice of hexagonal cells that are passive except for an inclusion of four central active cells surrounded by a buffer ring at half activity (Figure 2A). The buffer cells were used to prevent distortions associated with large differences in activity between cells. Mechanical anisotropy was created by applying forces of equal magnitude and opposite direction perpendicular to the left and right boundaries to induce anisotropic mechanical stresses in the tissue. Furthermore, both activity and viscoelastic relaxation were switched off until mechanical equilibrium was reached (see ‘Materials and methods’). The initial state is mechanically anisotropic (Figure 2B, top left), with differential tension in horizontal (h) vs. shoulder (s) junctions (Figure 3A), with being significantly larger than . Near the equilibrium point , and for , , the dynamics of the model tissue closely resembles the dynamics of the single-junction model. It is easy to show that , and mechanical anisotropy leads to anisotropic distribution of myosin, that is, (Figure 2B, bottom left). There is a range of applied pulling forces that, therefore, produce tensions in the system. In this regime, for a suitable (Equation 3) horizontal junctions are contractile, while shoulder junctions are extensile. Here, the factor of 1/2 is due to both active elements of a junction acting in parallel.
In steady state, both activity and viscoelasticity were switched on (Figure 2—videos 1–6), with timescales chosen such that , which is the biologically relevant regime. Figure 2B, second column, shows the system after . For the central horizontal active junction, contractility has been triggered, and the junction is steadily contracting at high myosin and against high tension (Figure 3B).
The active T1 transition is reached at , when the central junction shrinks to a point and a four-vertex is created (Figure 2B, third column). If the sum of the forces on the four-vertex is favourable for it to split (Spencer et al., 2017) in the vertical direction, a T1 event occurs, accompanied by myosin redistribution (see ‘Materials and methods’). The newly created vertical junction expands at intermediate values of myosin and tension (Figure 2B, fourth column, and Figure 3B). In contrast, aided by the anisotropic redistribution of myosin, the shoulder junctions are now strongly polarised and begin contracting, with higher myosin levels on the side of the junction belonging to the expanding pair of cells (Figure 3C). This expansion phase is followed by several secondary T1 events, the first of which typically occurs at a shoulder junction (Figure 2B, fifth column). During this propagation phase, there is a strong mechanical anisotropy in the direction of applied pulling forces. Together, the central T1 and the subsequent T1 events lead to substantial convergence–extension flow as can be seen qualitatively in the shape of the region formed by the 14 central and buffer cells (Figure 2, medium and dark grey).
Timescales of active T1 events and local convergence–extension strain
We proceed to quantify T1 transitions using the method introduced by Graner et al., 2008 (for a brief summary, see ‘Materials and methods’). First, using the topological tensor, (Equation 19), we measured the time and orientation of the T1 transition along the direction of the central junction of the active region (hereafter the ‘central T1’). Figure 4A shows the probability of a central T1 as a function of and , with other parameters held constant at the same values as in Figure 2. The probability was computed from simulations with different realisations of the myosin noise as the fraction of simulations where the first observed T1 was along the central active junction rather than elsewhere in the system. The probability of any T1 in the system was also measured, with the red contour in Figure 4A corresponding to 50% of realisations having a T1.

Existence and timescales of T1 transitions in the vertex model with active junctions as a function of and , averaged over simulations with different realisations of the myosin noise.
(A) Probability of a central T1 transition. The red line is the 50% probability contour of any T1 occurring in the simulation. (B) Typical timescale for the T1 transition to occur, measured as the length of the contraction phase. The other parameters are the same as in Figure 2. (C) Magnitude of the convergence–extension deformation as a function of and characterised by measuring induced by the T1 transition.
These results show that there is an absolute lower threshold, , for any form of T1 to occur. This is qualitatively consistent with both sides of the junction acting as two parallel instances of the single-junction model. Second, there is an optimal range of applied pulling forces for central T1 transitions, . Within this range, for the central and shoulder junctions, and during the initial contracting phase, they are contractile and extensile, respectively. Outside this optimal regime, the probability for central T1 events decreases rapidly, though T1 transitions still occur elsewhere for large values of .
The orientation of the central T1 transition in the optimal region is shown in Figure 3D. It was computed from the angle of the principal direction of corresponding to the largest eigenvalue before and after the T1. One can immediately observe that the transition is highly symmetric, with the orientation of the shrinking (disappearing) junction being very close to horizontal, and the expanding (appearing) junction very close to vertical. Outside the optimal regime, this symmetry disappears, with T1 events of non-central junctions occurring in different directions, similar to the fully random case discussed below.
Figure 4B shows the time to the first transition, , measured from the point when the activity and viscoelastic relaxation were switched on. The analysis was limited to central T1 transitions at parameter values where at least 25% of simulations yield a central T1, with and . One immediately observes that , consistent with a T1 dynamics being dominated by myosin activation and viscoelastic relaxation. We find that has a minimum in the same optimal region identified in Figure 4A, with rising both for larger and smaller values of . Furthermore, increasing beyond gradually reduces , consistent with active contractions becoming stronger. For , cells shapes become increasingly distorted, suggesting that the model is no longer applicable.
We now quantify convergence–extension generated by the model. The shear component of the total integrated strain given in Equation 24 was computed for the 14 cells comprising the central and buffer regions (Figure 2A). Figure 4C shows as a function of , evaluated at and averaged over . This point approximately corresponds to the empirically determined peak of convergence–extension in the optimal region of applied pulling forces (see sample time traces in Figure 3—figure supplement 1). From Figure 4C it is evident that the model generates pronounced convergence–extension. Without activity, that is, for , the system extends in the direction of the applied pulling force and contracts perpendicular to it with and, as expected, it increases with . For , decreases, indicating that the activity acts against the applied pulling force. As increases beyond a critical value, , active forces are strong enough to counteract the pulling forces and the system shrinks against the external load and extends in the perpendicular direction, that is, . During this process there are no significant changes of the area, that is, . This mechanism is, however, effective only for a range of values of . If is insufficiently strong, the myosin–tension feedback loop does not fully activate (Figure 2—video 4). Conversely, if is too strong, the feedback loop is active, but all junctions are activated, stiffening the tissue (Figure 2—video 5).
We conclude the analysis of a hexagonal tissue patch by investigating the influence of and (Figure 2—videos 7–13). Figure 5A shows the probability of a central T1 for , , that is, deep in the optimal region, as a function of and . Here, we only consider the biologically plausible regime with , with proportionally scaled myosin noise , and we excluded very small values of where noise dominates. We find that the mechanism for T1 events is very robust over 2–3 orders of magnitude in both and , with a guaranteed T1 transition in most of the parameter space. The only exception is the regime , that is, very slow myosin dynamics compared to viscous relaxation, where the system fails to develop anisotropy. All simulations generated at least one T1, and there is no equivalent of the red contour in Figure 4A. In Figure 5B, we show the timescale of the T1 transition. We find the T1 timescales as and . This influence of both timescales is consistent with the complex interplay between myosin activation on central and shoulder junctions. Finally, in Figure 5C, we show the convergence–extension strain as a function of for a range of values of . The effectiveness of the T1 mechanism is largely independent of and , and we have , the same value as in Figure 4C. This additional robustness of the model to varying timescales is likely due to being in a quasistatic regime where the elastic deformation timescale is much shorter than any other timescale.

Robustness of the T1 mechanism as a function of and for and .
(A) Probability of a central T1 event, averaged over simulations with different realisations of the myosin noise. The probability of any T1 event is 1 throughout. (B) Contraction time to collapse for the central T1 as a function of , for different values of , for points where a central T1 event occurred in at least 25% of simulations. (C) Peak of the total convergence extension strain , showing very weak dependence on viscoelastic and myosin timescales. Shading in panels (B) and (C) indicates the standard error of the mean.
Convergence–extension in a fully active random patch
The hexagonal tissue patch is convenient to analyse isolated active T1 events. Cells in real epithelia, however, do not have regular shapes packed in crystalline order. We, therefore, studied a patch of 520 active and 80 passive cells generated from a centroidal Voronoi tessellation starting from points placed at random in the simulation box. As in the case of the hexagonal patch, we applied pulling forces of constant magnitude on left and right boundaries in order to generate internal anisotropic stresses that mimic the mechanical conditions in the sickle region of the embryo (see next section). First, we investigated the occurrence of active T1 transitions and the emergence of convergence–extension as a function of and , at fixed and , that is, in the same region of parameter space as in Figure 4 (see also Figure 6). To be consistent with the hexagonal patch, the unit of length here is set by the length of a regular hexagon of area , where is the initial patch size (see ‘Materials and methods’).

Disordered active tissue at time as a function of the magnitude of the pulling force and activity .
The region of convergence–extension is at the centre of the diagram, around and . The remaining parameters are the same as in Figure 2.
We observed T1 transitions for all simulated systems, though their rate rapidly increased when either or were increased. Unlike in the case of the active inclusion in a hexagonal passive patch, there is no clear threshold for active T1 transitions. Instead, parts of the tissue are activated, and we observed the emergence of pronounced myosin cables and accompanying tension chains (e.g. Figure 6, middle row). The system starts to experience significant flow and convergence–extension from , significantly below the values observed in the hexagonal patch with a single active inclusion where, depending on the magnitude of the applied force, T1 events start to appear for above . This suggests cooperative rearrangements in the tissue, and we indeed see evidence of serial active T1 transitions along tension chains (Figure 6—video 2). This suggests that the individual junction feedback mechanism together with mechanical (as opposed to chemical) propagation of myosin activation is a key ingredient in the formation of the myosin cables that have been observed to accompany convergence–extension flow in chick embryo gastrulation (Rozbicki et al., 2015) and in Drosophila germ band extension (Jacinto et al., 2002).
The system starts to rearrange from , and for , one observes a highly active state with many uncorrelated rearrangements and implausible cell shapes. The region with realistic cell shapes is significantly below the hexagonal patch with a single active inclusion, where, depending on the magnitude of the applied force, T1 events start to appear for above .
In Figure 7A, we show a snapshot of a random patch long after () activity was switched on for and , in what emerges to be the optimal region for convergence–extension. Oriented myosin cables accompanied by tension chains form in the initial stage of the simulation at after activity was turned on. At longer times, the pronounced orientation of myosin and tension decreases (but does not disappear entirely), and convergence–extension stops (see ‘Materials and methods’ for details). The convergence–extension process is accompanied by active T1 events that predominately occur perpendicular to the direction of the external pulling, as shown in the orientational histogram (Figure 7B).

Characterisation of convergence-extension in a random tissue patch.
(A) Snapshot of a random tissue patch for and at ,that is during the convergence–extension flow. Red arrow indicates that a constant pulling force is applied throughout the entire simulation. (B) Angular histogram of T1 events for same values of and . Cells and junctions are coloured as in Figure 6. (C) Magnitude of the convergence–extension deformation as a function of characterised by measuring induced by the T1 transition. (D) Anisotropy along (, solid line) and perpendicular to (, dashed line) the direction of the external pulling force for myosin (green), mechanical stress (red), and shape tensor (black) as functions of for at . In (C) and (D), each point was averaged over independent samples and the error bar is smaller than the symbol size.
Figure 7C quantifies the amount of convergence–extension by measuring the difference of total integrated strain in directions along and perpendicular to the direction of the applied pulling force, that is, . It is evident that like the hexagonal case, the random patch undergoes substantial convergence–extension over a range of activities. The process in accompanied by spatially anisotropic distribution of myosin, mechanical stress, and cell shapes (Figure 7D), which we measured through the eigenvalues of the myosin (Equation 26) and tension (Equation 25) tensors, and the eigenvalues of the shape tensor (Equation 20).
Comparison with myosin-driven intercalations observed in the sickle-shaped mesoderm precursor domain of chick embryo at the onset of gastrulation
Previous experiments have shown that the embryo scale tissue flows driving the lateral to medial convergence and posterior anterior extension of the posterior sickle-shaped mesendoderm precursor domain that results in its transformation into the streak is driven by directed cell intercalations. These cells undergoing directed intercalations are characterised by an anisotropic cell shape and an anisotropic distribution of active myosin II in their cell junctions, organised in supercellular chains of variable length (Figure 8A and A’; Rozbicki et al., 2015; Chuai et al., 2023). We here make a qualitative comparison of the myosin-driven intercalations in the model with the myosin-dependent directed cell intercalations driving tissue deformations in the sickle-shaped mesendoderm precursor region of the early gastrulation stage chick embryo. Figure 8A shows an image of a typical embryo prior to primitive streak formation with a patch magnified in panel A’. Both the myosin (green) and cell shapes exhibit a clear anisotropy (indicated by green arrows) along the lateral–medial axis of the mesendoderm precursor domain. The myosin cables are believed to reflect and generate the tension chains observed in simulations.

Analysis of the tissue flows in the early-stage chick embryo.
(A) Image of a typical early-stage chick embryo prior to the gastrulation (i.e. primitive streak formation). The primitive streak will form along the yellow dashed line. The direction of myosin anisotropy is shown by the green double-headed arrows, and the direction of the tissue flow is indicated by the red arrows. x-axis is chosen to coincide with the long direction of the sickle-shaped active region in the embryo’s posterior (Rozbicki et al., 2015). (A’) Zoom-in of the rectangular region on the posterior side of the embryo; myosin II (green), actin (red), and nuclei (blue). (B) Measured distribution of the orientation of T1 events in a circular patch of diameter ≈ 190 tracked over the period of ≈ 6 hr (cf. model distribution in Figure 7B). Blue (red) denotes junctions that appear (disappear); arrows have the same meaning as in (A). (C) Angle (red) and magnitude (green) of tissue shape anisotropy (dots) and tissue flow (crosses) for rectangular patches along the sickle with corresponding anisotropy and flow patterns shown in (C’). Components of the elastic strain tensor (D) and the total integrated strain tensor (E; definition in Equation 24) as a function of time during first 4 hr of the streak formation for two central regions of the sickle (yellow stripes in C). Details of the analysis are given in ‘Materials and methods’.
We have quantified the experimental tissue convergence–extension as follows. Figure 8B shows an orientational histogram of T1 events in a rectangular region of size along the sickle tracked over a period of approximately 6 hr. The events were identified by calculating the tensor using segmented images (see ‘Materials and methods’ and Figure 8—video 1), and validated and corrected manually. In Figure 8C, we analysed tissue patches of diameter ≈ 190 μmm chosen sequentially along the sickle-shaped region. We computed tensors that measure shape anisotropy and strain rate (Graner et al., 2008) by tracking the patches along the tissue flow over approximately 4 hr (‘Materials and methods’). The anisotropy remains constant at around 20% and its direction is at the angle close to 90°, that is, along the sickle and orthogonal to the streak (Figure 8C’, green arrow). We also computed the mean flow magnitude and direction from the eigenvalues and eigenvectors of the integrated tensor. There is a pronounced spatial pattern to its direction, pointing towards the streak and parallel to the direction of anisotropy on outer parts of the sickle, and orthogonal to the sickle and its orientation in the middle of the sickle (Figure 8C’, red arrows). At the same time, the magnitude of flow peaks in the middle part of the sickle. This is consistent with the incipient flow to create the streak. In Figure 8D and E, we show the time dependence of spatial anisotropy and total strain, for two central patches. We see that anisotropy is along (i.e. perpendicular to the streak) and flow is along (i.e. along the streak), corresponding to convergence–extension flow. For comparison, we show the behaviour in an isotropic region anterior to the streak that shows non-directional intercalations and absence of tissue deformation (Figure 8—figure supplement 1, Figure 7—figure supplement 3C, and Figure 8—video 2). These observations show a close correspondence with the results of the model in the absence of imposed stress anisotropy, where the myosin-mediated active intercalations emerge spontaneously, leading to disorganised tissue patches with random intercalation directions (shown in Figure 7—figure supplement 3B).
Discussion
There are two key features of our model with active junctions. First, the anisotropy of myosin distribution is induced by anisotropic mechanical tension. Second, active contractions are triggered by tension-sensitive accumulation. The feedback loop between tension and myosin motor activity leads to contraction against and extension perpendicular to tension. This cellular mechanism has been suggested to be driving the tissue flows during primitive streak formation in avian embryos (Rozbicki et al., 2015), where there is no clear evidence of chemical prepatterning. Although it is yet to be experimentally confirmed, it is plausible that the symmetry breaking event that induces the initial myosin anisotropy occurs as a result of anisotropic tension combined with cell differentiation early in development. For example, in the chick embryo, Gdf3 is expressed in the sickle-shaped region in the posterior epiblast, and it is believed to play a role in triggering the contractions that initiate the large-scale tissue flows that subsequently lead to the formation of the primitive streak (Serrano Nájera and Weijer, 2020). The key conclusion of this study is that once the process has been initiated, chemical anisotropy emerges spontaneously and there is no need to impose it.
Although the model investigated here generates active T1 events, the process loses coherence after several T1 events, resulting in biologically implausible tissue shapes. This is very prominent in the toy case of regular hexagonal patch. While the problem is to some extent alleviated in patches of randomly shaped active cells, generating robust convergence–extension flows that would span scales of the entire embryo will be hard to achieve with the current model. One source of instability is likely that the post-T1 expansion of the junction, currently effectively modelled as a passive process, is not properly captured by the model and requires additional sources of activity to be considered. At the scale of the entire embryo, other cellular processes such as cell division, differentiation, and ingression all play non-trivial roles. These events have not been considered here.
Furthermore, recent results confirm that the mechanics of vertex models is complex (Bi et al., 2016; Yan and Bi, 2019; Tong et al., 2021), which makes the tissue-scale flows that emerge from active elements coupled to the vertex model hard to predict using continuum approaches. Active nematic flows have been observed and modelled in in vitro epithelial tissues (Saw et al., 2017), and due to the locally contractile and extensile dynamics it is plausible that large-scale flows predicted by this model belong to the same class of models. In real tissues, however, the mechanisms that regulate how cells coordinate their internal mechanical stress directions in order to produce a stable flow pattern are unknown. Active vertex models therefore provide valuable new insights into the intricate interplay between mechanical and biochemical processes that control the collective cell behaviours in epithelial tissues.
We also briefly discuss how this work relates to other recent models for active junction contractions and convergence extension. The single-junction model shares a lot of common features with the model of Dierkes et al., 2014. The key difference, however, is a different myosin–tension feedback mechanism and that the dashpot in Dierkes et al., 2014 was replaced with a Maxwell element. This was inspired by laser tweezer measurements of the response of cell–cell junctions in the Drosophila embryo to applied pulling force (Clément et al., 2017) which showed that the cellular junctions behave as a Maxwell viscoelastic material. The presence of an elastic spring attached to the dashpot introduces a viscoelastic timescale, which leads to the suppression of the oscillatory behaviour seen in Dierkes et al., 2014.
The model of Staddon et al., 2019 considered a Maxwell element subject to active contraction, with the spring constant of the cellular junctions that constantly remodels itself to match strain in the junctions. The tension, however, remodels only if the strain exceeds a threshold value. These two features combine to provide a simple mechanism by which the junction can undergo ratchet-like behaviour and contract to a T1 event, and where T1 events can be triggered by applying external forces. Staddon et al., however, do not explicitly include kinetic equations for molecular motors.
Furthermore, the model of Curran et al., 2017 introduced myosin-driven tension fluctuations in cell–cell junctions to cell ordering observed in experiments on the fly notum. They argued that myosin fluctuations combined with the isotropic distribution of myosin can either drive or inhibit T1 transitions. This suggests a robust mechanism by which epithelia can tune their properties. Similar conclusions about cell ordering have also been reached by Krajnc et al., 2021 using a model with a mechanical feedback between junction contraction and force generation. Neither of those studies, however, explored whether myosin fluctuations could lead to directional deformations as observed during convergence–extension.
The two-dimensional models of Wang et al., 2012 and Lan et al., 2015 provide detailed descriptions for coupling between chemical signalling and corresponding mechanical responses. While they were able to produce T1 events, chemical anisotropy in the cell was externally imposed by tuning concentrations of relevant molecular species based on the origination of the junctions. This was also the case for the two-dimensional version of the model of Staddon et al., 2019, which was able to produce convergence–extension flows, albeit with time-dependent activity imposed in a given direction. Meanwhile, Noll et al., 2017 have also introduced a vertex model with junctions that incorporate generic active feedback in a model for tissue contraction in Drosophila gastrulation. They did not, however, consider active T1 transitions.
Finally, we remark that pulling experiments on suspended cultured Madine Darby Canine Kidney (MDCK-II) cells (Harris et al., 2012) and on the Drosophila wing disk (Duda et al., 2019) do not show any topological transitions even for strains that exceed 50%. Instead, the deformation is accommodated by changes of cell shapes. In the case of MDCK layers, directional cell divisions also play an important role is releasing some of the mechanical stress introduced by the pull (Wyatt et al., 2015). While it is possible that T1 transitions would appear given sufficiently long time, it is unlikely that the mechanism proposed in this study would apply to these tissues. This is not unexpected since MDCK cells are transformed cells derived from adult tissues and likely not comparable to embryonic tissues that need to undergo large-scale highly coordinated shape changes involving massive cell rearrangements. The same argument holds for fly wing disk that also does not show large-scale anisotropic strain rates, driven by directed cell intercalations as are typical for embryonic tissues during gastrulation.
In summary, in this study we have introduced a mechanochemical model that describes the dynamics of active T1 transforms, that is, cell intercalation events that occur perpendicular to the externally applied mechanical stress. Such processes are believed to play a key role in the primitive streak formation in avian embryos. Crucially, this study suggests that mechanical propagation of activation of myosin is a key ingredient in formation of the myosin tension chains that have been observed to accompany convergence–extension flow in chick embryo gastrulation (Rozbicki et al., 2015) and in Drosophila germ band extension (Jacinto et al., 2002). Results of this study show a good qualitative agreement with measurements on early-stage chick embryos.
We conclude by observing that it is remarkable that models that share the common assumption of a feedback between activity and mechanical tension are able to describe a range of markedly different biological processes in different organisms. Although the details depend on the specific biological system and its molecular details, this suggests that there is a set of universal physical mechanisms that govern tissue-scale behaviours.
Materials and methods
Model setup and analysis
Single active junction
Request a detailed protocolTo understand the mechanism that couples the kinetics of myosin motors to the local mechanical tension and leads to the activation of contractility in cell–cell junctions, we first analyse a model for a single junction. The surrounding tissue is abstracted by assuming that it provides an elastic, tension-generating background against which the junction actively contracts. There are two key ingredients that make an active junction. First, the junction is viscoelastic. Pull-release optical tweezer experiments on cell–cell junctions in Drosophila and chick embryos have shown that cell–cell junctions have a viscoelastic response (Clément et al., 2017; Ferro et al., 2020). This means that the junction is able to remove imposed tension by remodelling itself. Second, the junction can generate tension. Tension is generated by myosin motors that form mini filaments slide actin filaments past each other. The single-junction model provides insight into the conditions under which the junction length can contract to zero and trigger a T1 transition.
Specifically, the single junction (Figure 1A) is modelled as an active mechanochemical system comprising three components connected in parallel (Figure 1—figure supplement 1A): (1) a viscoelastic SLS element, (2) a viscous dashpot, and (3) a tension-sensitive force-generating motor. The junction is subject to external tension produced and transmitted by the surrounding cells. The SLS element consists of an elastic spring of stiffness and rest length connected in parallel with a Maxwell element containing a spring of stiffness and rest length attached in series to a dashpot of viscosity (Larson, 1999). The spring captures the passive elastic response of the junction and models the effects of the surrounding tissue. In the two-dimensional model discussed below, this term arises naturally and accounts for changes in the cell area and perimeter. It is referred to as the elastic barrier. The Maxwell element models the viscoelastic nature of the junction and the ratio of viscosity and stiffness sets its relaxation timescale . This is the timescale over which the junction remodels and adjusts its length to that imposed by the external load. The dashpot with viscosity models dissipation with the environment. Finally, the junction is equipped with an active source of tension, which models the action of myosin motors.
Each molecular motor produces a force of magnitude . motors attached to actin filaments of the junction, therefore, generate tension . If the maximum possible number of attached motors is , , with and . The junction is, however, a part of the tissue and in the absence of perturbations the average steady-state value of motors attached to all junctions . It is, therefore, appropriate to model the active tension as . This expression can also be understood as the leading-order term in the expansion of around (Dierkes et al., 2014). For , , that is, tension acts to extend the junction. While the term may appear counterintuitive, it reflects the fact that if motors attached to the junction are depleted, contractions of the surrounding junctions produce stronger pull on it than it can resist, resulting in elongation.
Tension feeds back on the kinetics of association and dissociation of myosin to actin filaments, leading to the kinetic equation for ,
where is the association rate constant, assumed to be tension-independent with no limit of the total available myosin, and is the dissociation rate constant, assumed to be a monotonously decaying function of tension . The kinetics of actin-bound myosin, therefore, resembles that of a catch bond. This assumption is motivated by measurements of binding and unbinding rates of myosin motors on single actin filaments and has been shown to be described as a simple negative exponential dependence of the dissociation rate on applied tension (Veigel et al., 2003; Kovács et al., 2007). The dynamics of the junction is given by the following set of equations,
with a natural choice being a sigmoid curve,
where is the contribution to the myosin dissociation that does not depend on tension. The first equation describes the time evolution of the junction length due to the internal tension, , and external tension, . is the threshold tension, and controls the steepness of the curve in the vicinity of (Figure 1—figure supplement 1C). The minus sign in front of the first term on the right-hand side indicates that corresponds to a junction that is contracting, that is, for . The second equation accounts for the viscoelastic nature of the junction (Clément et al., 2017), that is, the rest length relaxes towards the actual length with a characteristic timescale . Finally, the third equation was obtained by dividing Equation 6 by , where is the timescale of myosin association. In general, binding and unbinding of molecular motors is a stochastic process and the third equation should also include stochastic terms. For simplicity, such terms were omitted here, but were included in the two-dimensional model. Furthermore, it is assumed that and are comparable in magnitude, that is, that which makes using a linear spring model appropriate despite the total length of the junction changing significantly as the junction collapses.
Finally, in all simulations of the single-junction model , , and were kept fixed and, therefore, length is measured in units of , time in units of , and force in units of . Parameters and their values used in the analysis of the single junction are listed in Table 1.
Vertex model with active junctions
Request a detailed protocolThe mechanical response of the tissue is modelled with the vertex model (Farhadifar et al., 2007; Fletcher et al., 2014). The associated mechanical energy is a function of the cell area and perimeter,
where and are the area and the perimeter of cell , respectively, and are the preferred area and perimeter, respectively (assumed to be the same for all cells), and the sum is over all cells. The first term in Equation 9 accounts for three-dimensional incompressibility of cells, and is the corresponding elastic modulus of the cell . The second term in Equation 9 contains a combination of actomyosin contractility in the cell cortex and intercellular adhesions, where is the contractility modulus of cell (Farhadifar et al., 2007).
All inertial effects were neglected, and, in line with the existing literature, it is assumed that the friction can be modelled as viscous drag on each vertex. The equation of motion for vertex is, therefore, a balance between friction and mechanical forces,
where is the friction coefficient, and accounts for all active forces. Stochastic forces are, however, omitted since those do not qualitatively affect the dynamics at timescales of interest. Inserting Equation 9 into Equation 11 leads to
with and being the magnitudes of, respectively, the area and perimeter contributions to the force due to junction . The subscript () denotes the cell to the left (right) of the junction when facing in the direction of . is the unit-length vector perpendicular to the plane of the tissue and the vector points along the junction away from vertex , with , is the hydrostatic pressure on the cell C, and . Finally, the sum is over all junctions that originate at vertex and terms in the sum appear in counterclockwise order (Figure 2—figure supplement 1A).
Activity is introduced by assuming that each junction contains two active elements supplied by the two cells sharing it. Furthermore, each cell is assumed to have a finite pool of myosin, . The finite pool of myosin acts to introduce correlations in the distribution of junctional myosin within each cell. In other words, it models the mechanism by which, for example, depletion of myosin on a given junction is correlated to myosin accumulation on its neighbouring junctions. This is of central importance to establish myosin anisotropy within the cell. There are clearly many alternative ways to model coupling of myosin on different junctions within a cell. The model used here, however, requires a minimal number of parameters while being biologically plausible.
Of this total myosin, , a fraction is assumed to be activated, that is, bound the junctions, and thus depleted from the pool. Here, is the number of junctions shared by the cell . The association rate of myosin to an individual junction is proportional to . As in the single-junction model, the dissociation rate is proportional to the amount of myosin bound to the junction , modulated by a tension-dependent function . To match the steady-state value of myosin of the single junction where , the prefactor of the unbinding term also needs to be . Then, to lowest order in the contributions of the coupled junctions, . We can recover at by choosing small but finite values of ( in simulations). This leads to the kinetic equation for the myosin motor attached to the junction by cell ,
where it is assumed that and are same for all cells. Without loss of generality, it is possible to set . is the inverse rate of myosin binding, that is, the timescale of attachment of myosin motors and is a random white noise with zero mean and variance
which accounts for the stochastic nature of myosin binding and unbinding. The noise term is important for the system to be able to break the symmetry imposed by using regular hexagonal tilings. It is, however, not strictly necessary to introduce noise to the myosin kinetics, but instead consider, for example, that mechanical properties of the cells are randomly distributed. While the quantitative results would be affected, such a model is not expected to have qualitatively different behaviour compared to what is discussed here.
Equation 9 models the passive elastic response of the tissue. In order to achieve an active T1 event, remodelling needs to be present, that is, the system must be viscoelastic. There are various ways to include viscoelastic effects into the vertex model. In order to be consistent with the single-junction model, it is further assumed that the junction has a viscoelastic contribution to the tension, , where is the spring constant analogue to the elastic part of the Maxwell element in the single-junction model and is the time-dependent rest length with dynamics,
where is the characteristic timescale for viscoelastic remodelling. The full expression for the tension of junction is
As above, the superscript () denotes the cell to the left (right) of the junction when facing in the direction of . is the activity of the junction produced by the cell , that is, it is a constant with units of force that measures the strength of the mechanochemical coupling and m0 has the same meaning as in the single-junction model. For simplicity, for all cells and for all junctions and both parameters were kept constant in all simulations. The unit of length, , was chosen to be the length of the side of a regular hexagon, which allows us to set the unit of time , and the unit of force . Equation 11, with given by Equation 15, and Equation 12 and Equation 14 describe the dynamics of the vertex model with active junctions.
The equations of motion for the two-dimensional model were integrated numerically for a rectangular patch made of hexagonal cells and patch made of randomly shaped cells using open boundary conditions. Mechanical anisotropy was created by applying a force to the left and right boundary vertices (Figure 2A), where the positive (negative) sign corresponds to the right (left) boundary. The initial pull was applied for with both activity and viscoelastic relaxation switched off, sufficient to reach mechanical equilibrium. Once the system reached an equilibrium stretched state, activity was switched on in 14 central cells in the case of the hexagonal patch and for 520 cells in the case of the random tissue patch. For the hexagonal case, activity was set to in 4 cells and to in 10 ‘buffer’ cells surrounding them, in order to suppress numerical instabilities at contacts between active and passive cells. For the random patch, activity was set to in all cells except for a single-cell thick layer of boundary cells that were kept passive to prevent artefacts due to tension chains reaching the sample boundary. An external pulling force of constant magnitude was applied throughout the entire simulation. The active system was simulated for using time step . In the hexagonal case, the dynamics of the central horizontal junction and the shoulder junctions shared by the four central active cells was monitored. The orientation of the hexagonal lattice was chosen such that the central active junction was parallel to the direction of the applied external force. In the case of the random patch, the dynamics of all junctions shared by active cells was monitored.
Parameters and their values used in the analysis of the vertex model with active junctions are listed in Table 2.
Characterisation of T1 transitions and tissue flow
Request a detailed protocolThe T1 events were implemented following the procedure proposed by Spencer et al., 2017, where a junction shorter than collapses into a fourfold vertex. The fourfold vertex either remains stable or it is resolved into two threefold vertices based on the sum of the forces acting along the four junctions connected to it. Importantly, the direction of the new junction is not imposed, and this procedure does not generally lead to a new junction orthogonal to the collapsed one. For simplicity, vertices with connectivity greater than four were not considered.
If a T1 transition occurs, it is necessary to assign myosin to the newly created junction and redistribute the myosin associated to the collapsed junction to the surrounding junctions. As illustrated in and using the notation introduced in Figure 2—figure supplement 1B, the myosin of the collapsing junction is stored right before the junction collapses. After the T1 transition, the myosin on the new junction was set to . Both inner shoulder junctions increased their myosin by , while the myosin on the outer junctions was reduced by , where m2 is its myosin right before the T1 transition and the minimum function ensures that myosin remains . The myosin redistribution procedure reduces artificial jumps of the junctional myosin as the system progresses through T1 transition and also provides a natural way to handle the finite myosin pool. After the transition, of the new junction is set to .
To identify and characterise T1 transitions, and quantify the associated deformation of the model tissue, we used the analysis method introduced by Graner et al., 2008. For cellular patterns, three tensors are defined, texture (), geometrical texture change (), and topological texture change (). Tensor describes the shape of the current cell configuration, and it is defined as
where () is the () component of the vector connecting centroids of two neighbouring cells at positions and , respectively, is an average over pairs of neighbours, that is, cell–cell contacts, and . Using the same notation, the tensor describes shape changes of the cell configuration during a time interval , and it is defined as
where the superscript T denotes the matrix transpose and
Finally, the tensor identifies T1 transitions by quantifying topological changes of the cell configuration in the time interval via tracking appearance and disappearance of contacts between cells. It is defined as
where () is the average over contacts that appeared (disappeared) during the time interval . The total number of contacts that appear and disappear is typically much smaller than , which means that data can be quite noisy. While all three tensors can be calculated for individual cells, the averaging is meant to be carried over a mesoscopic region. We, therefore, averaged over the 14 central active cells (Figure 2A, dark-shaded cells), as well as over an ensemble of noise realisations. For the disordered tissue, we averaged over all active cells.
Tensors , , and all involve averaging over cell–cell contacts and, therefore, describe a discrete system. In order to make connections to continuous deformations of the entire tissue, one introduces their continuous counterparts, the statistical strain tensor (, the velocity gradient tensor (, and the tensor of the rate plastic deformations (i.e. topological rearrangements rate) () (Graner et al., 2008). These tensors are defined as
where is the texture tensor of an arbitrary reference configuration for the statistical relative strain – we chose the initial undeformed configuration. Further,
and
In the case when variations in can be neglected, one can show that (Graner et al., 2008)
where is the corotational derivative (Larson, 1999). This equation just states that the velocity gradient is a sum of two contributions, reversible changes of the internal strain and the rate of irreversible plastic rearrangements. The significance of Equation 23 is that if one assumes that there are no plastic events other than T1 transitions, it is possible to obtain the total strain of the system by integrating the tensor over time, that is,
where t0 is the time when the activity is switched on and is the total simulation time. The difference of the and components of (i.e. total strains in the and direction, respectively) was used as the measure of the amount of convergence–extension in the tissue induced by the active T1 transition. Figure 3—figure supplement 1 shows ensemble-averaged time traces of and the time-integrated total strain and plastic strain through the active T1 transition and the subsequent propagation phase.
Convergence–extension in a patch of randomly shaped active cells
Request a detailed protocolWe used the two-dimensional model with parameters for the cell mechanics and the myosin feedback loop as defined in Table 2 on a tissue patch with cells of random shapes. The patch was generated using an iterative procedure. First, points were placed at random in a square region of size chosen such that the average cell area is equal to that of the hexagonal patch, that is, . During the initialisation process, we ensured that no two points are closer than to each other. Once all points are placed in the box, the Voronoi diagram was constructed and centroids of each cell of the Voronoi diagrams were used as the seeds for constructing Voronoi diagram in the next iteration. This was repeated until no centroid moved more than between two consecutive iterations, resulting in a centroidal Voronoi tessellation. Centroidal Voronoi tessellations have several convenient properties. For example, typically, there are no outliers (i.e. very large or very small cells are unlikely) and the distribution of cell neighbours is remarkably similar to actual epithelia. Finally, to avoid the spontaneous formation of an actomyosin cable at the outside border, we used a setup where the interior cells are active, and a one-cell-thick layer of passive cells forms the boundary.
We used the same stretching protocol as for the ordered patch. At every vertex on the left and right boundaries, the force was exerted during the entire duration of the simulation. The passive tissue was first made anisotropic by stretching it for with activity and the viscoelastic relaxation turned off. We then turned on activity and viscoelasticity and simulated the system for . This protocol was repeated for a range of activities , viscous relaxation times , and myosin times (not shown since like we observed for the single active T1, results for convergence–extension were largely independent of and and quantitatively similar to Figure 7). All results were averaged over different random initial configurations and with different random number generator seeds for the myosin noise. The results shown in Figure 7C and D are reported with a 95% confidence interval computed using bootstrapping by resampling 103 times.
Figure 7 shows tissue dynamics for and , which is the optimal region for convergence–extension. We measured convergence–extension using the integral of , where the calculation was done over the entire active region, and the starting point was the time when the activity was turned on. For the elastic strain , we used the undeformed disordered initial condition to define the reference strain . Figure 7—figure supplement 1 shows the time traces of the ensemble averaged tensors for elastic strain , integrated total strain , and integrated plastic strain .
We also quantified the amount of tension and myosin anisotropy using the tissue-averaged stress tensor components
and
To measure convergence–extension strain and anisotropy of shape, myosin, and tension, we chose , which is the time point where convergence–extension and anisotropy stabilise, giving the results shown in Figure 7C and D.
To measure orientations of T1s, we first identified cells involved in a T1 transition by diagonalising the tensor associated with active cells at a given instance in time. The signature of a T1 event is a non-zero , and the sign of its trace determines if a junction appeared or disappeared. The eigenvector corresponding to the largest eigenvalue determines the direction of the event. We tracked all of such events over the simulation runtime and generated polar histograms such as the one shown in Figure 7B. As shown in Figure 7—figure supplement 2A for and , there are temporally strongly correlated back-and-forth T1 transitions at the same angle. These correspond to four cells flipping back and forth through a T1 transition, or ‘flickering’ in the videos. As these events are artefacts of the simulation, we excluded them from the data. Figure 7—figure supplement 2B, which was averaged over independent simulations, shows that flickering events have mostly been filtered out of the dataset.
Figure 7—figure supplement 2 shows filtered T1 histograms for other mechanical conditions, outside the parameter region of where convergence–extension occurs, averaged over independent simulations. In particular, we also include a passive tissue that flows in the direction of the pulling, and an active isotropic tissue without applied forces where the T1 distribution is isotropic. This last situation strongly resembles the observed T1 distribution in the anterior region of the streak (Figure 7—figure supplement 3C). As shown in Figure 8—figure supplement 1, this is also a mostly isotropic tissue.
Experimental data analysis
Experimental data
Request a detailed protocolActive myosin (phosphorylated myosin light chain) and actin staining in fixed embryos was performed as described in Rozbicki et al., 2015. Embryos of a transgenic chick line with cell membranes of all cells in the embryonic and extra embryonic tissues labelled with a green fluorescent protein tag (myr-EGFP) were live-imaged using a dedicated light-sheet microscope as described previously (Rozbicki et al., 2015). The microscope produces cell-resolution images of the entire embryo with time resolution of 3 min between frames for periods up to 16 hr (stage EGXIII-HH4). rectangular areas of interest of size ≈ 255 × 220 were chosen to lie next to each other along the sickle-shaped mesendoderm precursor region in the embryo’s posterior, perpendicular to the direction of the forming primitive streak. The two central sections were chosen to lie in the middle of the sickle region that initiates the formation of the streak. Each area of interest was tracked for ≈ 4.6 hr covering the onset of the flows driving streak formation. The average motion of the area of interest, determined via particle image velocimetry using PIVlab (Thielicke and Sonntag, 2021), was used to track its displacement to be able to follow the same patch of cells over time. After this, these image time series were bandpass filtered to remove some noise followed by segmentation using the watershed algorithm in MATLAB. To follow the same cells over time, their centroid positions calculated form the segmentation of the first image of the time series are projected forward to the next frame using a newly calculated high-resolution velocity field between these successive images and using these as seed points for the segmentation of the next image (Rozbicki et al., 2015). This procedure allowed us to track individual cells between consecutive time frames and determine changes in neighbours over time as well as determine the directions of appearing and disappearing junctions associated with T1 transitions. The , , and tensors in a particular area of interest were averaged for all cells in centred circular domains of 190 µm diameter between successive pairs of segmented images and averaged over 30 min time intervals.
The spatially averaged texture tensor allowed us to both compute and also to directly quantify shape anisotropy from the eigenvalues and eigenvectors of the time-averaged , where and label the large and small components, respectively. We defined the dimensionless shape anisotropy shown in Figure 8C in the main text as
and the shape anisotropy direction as the angle makes with the lab frame x-axis.
Unlike in the simulation, in the real embryo, cells divide, ingress, and flow in and out of the region of interest. Therefore, it is not immediately clear what to use as the reference texture tensor when computing tensor. For simplicity, we chose the isotropic tensor constructed from the time-averaged eigenvalues as
To compute the flow anisotropy, we measured the integrated total strain tensor
analogous to the simulated tissue patches. Similar to , we can then compute the eigenvalues and eigenvectors of the final , where we again label the large and small components and , respectively. Here, we now typically have and , corresponding to the convergence and extension directions of the tissue, respectively. We compute the total flow magnitude as
and the flow anisotropy direction as the angle makes with the lab frame x-axis (Figure 8C).
For comparison, in Figure 8—figure supplement 1, we show the integrated and the tensor for a region in the anterior of the embryo which does not undergo convergence–extension.
The phosphomyosin light chain staining patterns shown in Figure 8A are representative of the patterns observed in over 100 embryos of similar developmental stages. The intercalation data shown in Figure 8C are based on analysis of a light-sheet microscopy sequence of a single embryo, representative of over 50 experiments of control embryos.
Fertile eggs (Shaver brown) were obtained from Henri Stewart & Co Ltd, UK. Fertile Myr-GFP eggs were obtained from the National Avian Research Facility at the Roslin Institute, University of Edinburgh, UK.
Data availability
The current manuscript is primarily a computational study, so no data have been generated for this manuscript. Modelling code is publically (GNU public license v 2.0) available on GitHub at: https://github.com/sknepneklab/ActiveJunctionModel (copy archived at Sknepnek, 2023). The experimental data presented in Figure 8 and Figure 8—figure supplement 1 has been generated as described in the methods section.
References
-
Physical models of collective cell migrationAnnual Review of Condensed Matter Physics 11:77–101.https://doi.org/10.1146/annurev-conmatphys-031218-013516
-
A density-independent rigidity transition in biological tissuesNature Physics 11:1074–1079.https://doi.org/10.1038/nphys3471
-
Motility-driven glass and jamming transitions in biological tissuesPhysical Review. X 6:021011.https://doi.org/10.1103/PhysRevX.6.021011
-
Local and tissue-scale forces drive oriented junction growth during tissue extensionNature Cell Biology 17:1247–1258.https://doi.org/10.1038/ncb3226
-
The reaction-limited kinetics of membrane-to-surface adhesion and detachmentProceedings of the Royal Society of London. Series B, Biological Sciences 234:55–83.https://doi.org/10.1098/rspb.1988.0038
-
Spontaneous oscillations of elastic contractile materials with turnoverPhysical Review Letters 113:14.https://doi.org/10.1103/PhysRevLett.113.148102
-
Myosin II dynamics are regulated by tension in intercalating cellsDevelopmental Cell 17:736–743.https://doi.org/10.1016/j.devcel.2009.09.003
-
Vertex models of epithelial morphogenesisBiophysical Journal 106:2291–2304.https://doi.org/10.1016/j.bpj.2013.11.4498
-
Discrete rearranging disordered patterns, part I: robust statistical tools in two or three dimensionsThe European Physical Journal E 25:349–369.https://doi.org/10.1140/epje/i2007-10298-8
-
Patterned mechanical feedback establishes a global myosin gradientNature Communications 13:7050.https://doi.org/10.1038/s41467-022-34518-9
-
Coming to consensus: a unifying model emerges for convergent extensionDevelopmental Cell 46:389–396.https://doi.org/10.1016/j.devcel.2018.08.003
-
Mechanisms of convergence and extension by cell intercalationPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 355:897–922.https://doi.org/10.1098/rstb.2000.0626
-
Forces directing germ-band extension in Drosophila embryosMechanisms of Development 144:11–22.https://doi.org/10.1016/j.mod.2016.12.001
-
Active instability and nonlinear dynamics of cell-cell junctionsPhysical Review Letters 127:19.https://doi.org/10.1103/PhysRevLett.127.198103
-
Unjamming and cell shape in the asthmatic airway epitheliumNature Materials 14:1040–1048.https://doi.org/10.1038/nmat4357
-
Theoretical aspects of the biological catch bondAccounts of Chemical Research 42:693–703.https://doi.org/10.1021/ar800202z
-
Nature and anisotropy of cortical forces orienting Drosophila tissue morphogenesisNature Cell Biology 10:1401–1410.https://doi.org/10.1038/ncb1798
-
Cellular processes driving gastrulation in the avian embryoMechanisms of Development 163:103624.https://doi.org/10.1016/j.mod.2020.103624
-
SoftwareAJM (active junction model), version swh:1:rev:02b61a131f725276e83974c1fe19401074e2ad8bSoftware Heritage.
-
Vertex stability and topological transitions in vertex models of foams and epitheliaThe European Physical Journal. E, Soft Matter 40:2.https://doi.org/10.1140/epje/i2017-11489-4
-
Mechanosensitive junction remodeling promotes robust epithelial morphogenesisBiophysical Journal 117:1739–1750.https://doi.org/10.1016/j.bpj.2019.09.027
-
Mechanics and remodelling of cell packings in epitheliaThe European Physical Journal E 33:117–127.https://doi.org/10.1140/epje/i2010-10677-0
-
Particle image velocimetry for Matlab: accuracy and enhanced algorithms in pivlabJournal of Open Research Software 9:12.https://doi.org/10.5334/jors.334
-
Biophysics of catch bondsAnnual Review of Biophysics 37:399–416.https://doi.org/10.1146/annurev.biophys.37.032807.125804
-
Linear viscoelastic properties of the vertex model for epithelial tissuesPLOS Computational Biology 18:e1010135.https://doi.org/10.1371/journal.pcbi.1010135
-
A cell-level biomechanical model of Drosophila dorsal closureBiophysical Journal 103:2265–2274.https://doi.org/10.1016/j.bpj.2012.09.036
-
Multicellular rosettes drive fluid-solid transition in epithelial tissuesPhysical Review X 9:9.011029.https://doi.org/10.1103/PhysRevX.9.011029
Decision letter
-
Karsten KruseReviewing Editor; University of Geneva, Switzerland
-
Aleksandra M WalczakSenior Editor; CNRS, France
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Decision letter after peer review:
Thank you for submitting your article "Generating active T1 transitions through mechanochemical feedback" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by Aleksandra Walczak as the Senior Editor. The reviewers have opted to remain anonymous.
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
Essential revisions:
1) Improve the comparison with experiments (see reports).
2) Argue more convincingly for biological relevance.
3) Put the positive feedback tension/myosin more into context and discuss possible effects of mechanisms that are an alternative to the one you implemented (see reports).
Reviewer #1 (Recommendations for the authors):
This paper presents a numerical study of tissue convergent extension during epithelial morphogenesis. It builds upon several earlier works, including so-called vertex models, and single junction mechanical models. The specificity of the model presented lies in its active component: a positive feedback exists between Myosin activity and mechanical tension in each junction. This leads to interesting tissue-scale dynamics. Upon tissue pulling, the authors report conditions for which intercalation occurs perpendicular to the pulling direction. They provide a rather detailed analysis of the model and its numerical outcome, and a qualitative comparison with convergent extension in avian embryos.
Strengths:
The model is presented in a rigorous manner and with a lot of details and pedagogy. It is discussed in comparison with other models of the field, and the similarities and differences are clearly and honestly presented by the authors in their discussion.
The numerical study provides interesting and possibly counter-intuitive predictions in which tissue convergent extension can occur perpendicular to an external pulling axis, due to a positive feedback between myosin motors recruitment and tension which leads to the formation of Myosin chains (or cables) along the pulling axis.
The conditions in which this regime of tissue remodeling can occur (under the hypotheses of the model) are thoroughly investigated. Finally the authors go a bit beyond the sole numerical study, and analyze experimental data of the polarized tissue flow occurring upon primitive streak formation in avians.
Weaknesses:
Although the numerical investigation is thorough and the model hypotheses are clearly stated, the paper falls a bit short of convincing readers that the emergent behaviors observed in the model are of actual biological relevance.
The hypotheses and ingredients of the model are clearly stated (1) junctions are viscoelastic, (2) junctions are active, i.e. they can generate tension through myosin, and (3) positive feedback between tension and myosin. Yet the 3rd hypothesis, which gives its singularity to the model, is not discussed enough, and its relevance to well-known convergent extension systems is not analyzed. In the paper, the positive feedback arises from the decrease of unbinding rate upon tension increase. This "catch bond" hypothesis (note: the term catch-bond only appears once, in the methods) is at the core of all the results obtained in the paper, and this should be clearly stated and discussed. As for the positive feedback, one could also imagine that it is not the unbinding rate that decreases upon tension increase, but rather the myosin recruitment rate that increases. Would that be formally equivalent? The authors should address this point and its biological implications.
The comparison between the model outcome (notably, the T1 transitions preferably oriented in the direction perpendicular to pulling) and existing experimental data is a bit elusive. The authors propose a semi-quantitative comparison with convergent-extension occurring upon primitive streak formation in avians. The data and conclusion is not fully convincing, partly because of the polarization data itself (Figure 8C), and partly because it's not clear what would be the equivalent of the external pulling force in this self-polarizing flow.
The authors also mention the germband of Drosophila, but the model is, I believe, not relevant for that system. Although there is indeed an external force pulling the tissue (pulling posterior midgut), polarization is perpendicular and extension parallel to the pulling force (the opposite of the model outcome).
In that spirit, they should also discuss pulling experiments (fly wing disc / Mao lab, suspended epithelia / Charras lab) where it's not so clear that topological transitions occur upon pulling.
Finally, the paper builds upon already existing models of junction and tissue mechanics. Beyond the vertex model itself, the hypothesis of a feedback between myosin activity and tension has also been around in the field for some time, including in modeling papers (which authors acknowledge honestly).
In conclusion:
The paper presents a thorough numerical analysis of T1s and convergent extension upon tissue pulling when a positive feedback exists between Myosin and junction tension. Yet, the fact that this type of model and hypotheses are not completely new, combined to a lack of comparison to experimental data (except for avian gastrulation) makes it a bit difficult for the reader to assess the novelty of the findings presented in the paper.
I liked the presentation of the model a lot, the details given, the pedagogy, the comparison with what already existed (Discussion section), and the thorough analysis of the model outcome.
On the other hand, I find it a bit difficult to see whether the novelty is strong enough for publication in eLife:
First, because somewhat similar models exist in the field. According to the authors, the specificity here is the inclusion of the positive feedback for the active part of the model. Yet, other models that the authors cite include a similar type of feedback, although they might not specifically focus on topological transitions and convergent extension. In addition, the existence of such a feedback (and its possible consequences on tissue morphogenesis) has been around in the field for a while, under various forms (here, the authors chose a catch-bond type dissociation rate for myosin).
Second, because the comparison (even qualitative) with existing experimental data remains elusive.
Hence the major finding of this paper is that a tissue in which a positive feedback between tension and motor recruitment exists can display active intercalations and convergent extension perpendicular to an external pulling force. In which experimental systems this prediction could be relevant remains to be determined.
– Does the positive feedback tension/myosin need to stem from the catch bond model? Could one imagine increased motor recruitment open junction/cell stretch (This would be in line with some experimental data, including stretching experiments). Would that be formally equivalent to the catch bond, or would that change the outcome of the model? If yes, how?
– In which systems do people observe tissue extension along the direction perpendicular to pulling? The authors only refer to the avian embryo in which the nature of the flow (embryo-scale) makes it a bit difficult (at least for me) to interpret what is the chicken and what is the egg. A lot of data is now available out there for analysis. It would be great to see the same kind of analyses in other experimental systems. Notably, the fly germband (that the authors rapidly mention) seems to behave completely differently, with cables of Myosin perpendicular to the pulling axis (pulling force being exerted by the posterior midgut). How about convergent extension in the zebrafish or Xenopus embryos?
– What happens at the boundaries? Here we look at a central patch of cells. And somewhat related question, how would that work in a closed geometry (such as the fly embryo). I feel that the avian discussion (paragraph 1 of discussion) could be expanded with considerations on how the feedback can produce long range flows even though only a small patch of cells is initially "activated".
– Authors should cite Duda et al. (Dev Cell 2019) from the Mao group. In this article it is shown that Myosin polarizes and form cables upon external tissue stretch in the fly wing disc. This is fully relevant to their model and it is somewhat surprising that they do not cite it. Disclaimer: I'm not on that paper
– The authors present their model as "one of the simplest descriptions" that features the three main ingredients required (viscoelastic junctions, active junctions, positive feedback myosin/tension). What would be the simplest and why didn't they choose it?
– Authors should better explain the rationale of the "elastic barrier" term B(l-a) , as compared to the k(l-l0) term.
Reviewer #2 (Recommendations for the authors):
Sknepnek et al. study a model for junction dynamics in a two-dimensional vertex model, intended to describe epithelial mechanics. In the model, a concentration of myosin on each half-junction evolves according to a chemical balance equation and induces an active tension in the junction. The chemical balance equation incorporates a mechanical feedback, whereby myosin unbinding is inhibited at high tension. The authors first study a single junction and consider thresholds for the junction to fully contract. They then investigate the behaviour of a few « active » junctions in a passive network and a fully active vertex model. They find that at intermediate external pulling force and myosin-induced tension, the simulated epithelium can undergo convergence-extension against the pulling force.
The manuscript describes a work which is seriously and rigorously performed. The simulations are analysed thoroughly with a serious quantification effort. The authors have made an important effort of clarity in their study.
My central question is about the motivation of the manuscript. The central theme of the manuscript authors study « contrarian » active T1 transitions which induce convergence extension in the direction orthogonal to an externally applied stress. I would disagree with the authors that this is the only type of active T1 transition. For instance polarised transitions which are oriented by the epithelium chemical polarity and induce converge-extension in the absence of external tension, as is thought to be the case during germ-band elongation, could also be called « active ». From the manuscript it is not clear to me why inducing convergence extension as a response to external tension, against the external tension, is important or biologically relevant. I think the authors should significantly clarify the rationale for their study.
In that respect I am not sure that the comparison with experiment is very useful. My understanding is that essentially myosin cables in the chick embryos appear oriented along the direction of tissue contraction, which seems to make sense if these cables are actively contracting. I am not sure that this observation strongly supports the much more detailed model proposed by the authors; notably the notion that myosin polarises in response to external tension.
Judging from Figure 6, the authors also observe relatively limited convergence extension happening against the external force. I wonder if this may be due to myosin molecules playing both the role of the force-exerting molecule and the sensor in the model? For instance, possibly if a secondary molecule responding to external tension, and itself inducing myosin polarisation, a more permanent convergence-extension could be induced?
– I think the authors may have overlooked the reference « active instability and nonlinear dynamics of cell-cell junctions », Krajnc, et al, PRL, 2021. The approach and questions seem very close to the current manuscript. I think the results and approach of this study should be contrasted to the current work.
– When comparing Figure 7C and Figure 6, I had the impression that Figure 7C showed significant « negative » convergence extension for large values of β (0.75), why the authors say that this is not the case in the caption of Figure 6?
– The model assumes a constant number of available myosin molecules in each cell. It was not clear to me if this assumption is important for the results of the model?
– I do not understand the reason for the factor « z » in the second term of the right-hand side of Equation 5.
Reviewer #3 (Recommendations for the authors):
The authors propose a mechanism for the dynamics of cell-cell junctions under an applied external mechanical load. In their theoretical analysis they assume that myosin recruitment and therefore active contractile stress in a junction increases with an applied external load. They find that this mechanism can induce a T1 transition through which adjacent cells exchange neighbours such that an originally aligned junction is reorientated perpendicularly to an applied uniaxial extensile stress. By using a generalised vertex model they also find that this effect on an isolated cell-cell junction can lead to convergence-extension of a tissue, where the extension occurs perpendicularly to the applied uniaxial extensile stress. In both cases, the phenomenon is present for a broad range of parameter values. Eventually, the authors compare their theoretical results to convergence extension of a part of a tissue in an early chick embryo.
The theoretical analysis is carried out mostly through simulations and is done in a convincing way. The dynamics assumed for a single junction is based mostly on reasonable assumptions and differs from earlier works. It appears though that the assumption of a monotonic decrease of the myosin unbinding rate with the applied force needs more justification. I would assume that for sufficiently high forces, this rate will increase as is the case for catch bonds. Taking such an effect into account might limit the range of parameter values for which stress-induced convergence-extension is observed. Also, it was unclear to me whether all ingredients of the model where necessary to obtain the observed effect. In particular, I wonder what happens if either of the dashpots or the springs are eliminated? Finally, I do not understand why m_0 should describe effects of the surrounding. Shouldn't this effect be included in T_ext? Otherwise the theoretical analysis is convincing and justifies the conclusions of the authors.
I am less convinced by the comparison with experiments on gastrulation in early-stage chick embryos. Although the relative orientation between myosin cables and the direction of extension agrees with that in the theory, I fail to see where the external uniaxial extensile stress should come from. It thus appears that an essential element of the theory is not checked for in the experiments. The authors write themselves in ll 395: "Although it is yet to be experimentally confirmed, it is plausible that the symmetry breaking event that induces the initial myosin polarity occurs as a result of anisotropic tension combined with cell differentiation early in development." The authors should clarify this point, which seems to be central to the authors' idea that mechanical signalling drives convergence-extension in this case. Otherwise, the experimental data do not add much to the work.
The authors frequently use the term "polarisation" or "polarised" when referring to the alignment of the myosin cables, for example. This notion implies directionality, for example, a difference between left and right in Figure 2, which is not present. I would suggest that the authors rather use 'anisotropy'.
The text exhibits some jumps between figures that are not always referred to in order, for example between Figures2 and 3 or Figure 4C and 4B. The authors might want to adapt their text to the figures or vice versa.
In Figure 1A the term \β m is hard to see. Please, improve.
In Figure 3 the colour code is confusing. In A black, red, blue correspond to different junction classes, in B and C to different mechanical quantities and in D again to different junction classes but not the ones in A. Please, improve.
Ll 49 "With cellular behaviours being coordinated over thousands of cells in the case of the chick embryo, biochemical signalling alone is unlikely to account for the observed motion patterns." Even though I have a guess, I do not really understand this argument. Can you add some words to explain, why the long-range coordination is at odds with pure biochemical signalling?
L 218 side a -> side OF a
L 256 along the direction the central junction -> along the direction OF the central junction
Ll 274 could you add a panel to show this?
In Figures6 and 7A, convergence-extension is not really visible assuming that the initial state was square. Maybe you can show a sequence similar to Figure 2?
[Editors' note: further revisions were suggested prior to acceptance, as described below.]
Thank you for resubmitting your work entitled "Generating active T1 transitions through mechanochemical feedback" for further consideration by eLife. Your revised article has been evaluated by Aleksandra Walczak (Senior Editor) and a Reviewing Editor.
The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:
Your model doesn't seem to apply to the fly germband elongation, in which the invaginating midgut (pulling external force) and the T1s occur along the same axis, with cables perpendicular to the pulling axis. In your response, you insist that the pulling force is rather the ventral invagination than the posterior midgut (hence a force perpendicular to tissue elongation and T1s, and parallel to the cables, as in their model). This overlooks that the germband does not extend without the posterior midgut pulling (Torso), while it does without the ventral mesoderm invagination. In the paper, claims about their finding applying to germband extension are found in the intro and discussion, without other justification. These comments should be removed unless very strongly justified, as it seems that they are based on a direct misinterpretation of previous observations.
When the germband is mentioned you usually cite a paper that does NOT deal with the germband (Jacinto et al. 2002, which is about dorsal closure). And you also cite Duda et al., even though it's a wing disc paper, in which cables form parallel to the pulling force, (as in their model, but in contrast with the germband). This could be confusing and even detrimental to your work when readers familiar with fly morphogenesis read it.
https://doi.org/10.7554/eLife.79862.sa1Author response
Essential revisions:
1) Improve the comparison with experiments (see reports).
In the introduction of the revised manuscript, we have highlighted the experimental data on large scale tissue flows and the underlying oriented myosin driven intercalations observed in the chick embryo. These experiential observations directly inspired the tension-sensitive myosin recruitment model studied in this work. Furthermore, we have included the discussion about the catch-bond mechanism that has here been proposed as a simple mechanism of tension dependent myosin accumulation. We, however, also explicitly discuss that the critical requirement for this model to lead to active contraction is tension sensitive myosin accumulation and that mechanisms other than the proposed catch bond mechanism could be involved. Resolving the actual mechanism in detail is, however, beyond the scope of this work, and will have to await further experiments.
2) Argue more convincingly for biological relevance.
We have extended the explanation for the biological relevance for tension-dependent myosin mediated control of the direction and patterns of cell intercalations. This explanation is primarily based on the observation that the intercalations need to be organised robustly over very large (several millimetres) domains. In chicken embryos, there is a clear spatial pattern of propagation of the onset of motion, which suggests a rapidly propagating signal as well as the dynamic patterns of cell deformations that we have observed. Those correlate closely with the formation and direction of supercellular myosin cables.
3) Put the positive feedback tension/myosin more into context and discuss possible effects of mechanisms that are an alternative to the one you implemented (see reports).
In the revised manuscript, we discuss the catch bond idea in more detail and mention explicitly that there could be additional and/or other molecular mechanisms of implementing a tension dependent myosin recruitment.
Reviewer #1 (Recommendations for the authors):
This paper presents a numerical study of tissue convergent extension during epithelial morphogenesis. It builds upon several earlier works, including so-called vertex models, and single junction mechanical models. The specificity of the model presented lies in its active component: a positive feedback exists between Myosin activity and mechanical tension in each junction. This leads to interesting tissue-scale dynamics. Upon tissue pulling, the authors report conditions for which intercalation occurs perpendicular to the pulling direction. They provide a rather detailed analysis of the model and its numerical outcome, and a qualitative comparison with convergent extension in avian embryos.
We thank the review for their insightful comments on our work. Below, we address the specific issues raised by the reviewer.
Strengths:
The model is presented in a rigorous manner and with a lot of details and pedagogy. It is discussed in comparison with other models of the field, and the similarities and differences are clearly and honestly presented by the authors in their discussion.
The numerical study provides interesting and possibly counter-intuitive predictions in which tissue convergent extension can occur perpendicular to an external pulling axis, due to a positive feedback between myosin motors recruitment and tension which leads to the formation of Myosin chains (or cables) along the pulling axis.
The conditions in which this regime of tissue remodeling can occur (under the hypotheses of the model) are thoroughly investigated. Finally the authors go a bit beyond the sole numerical study, and analyze experimental data of the polarized tissue flow occurring upon primitive streak formation in avians.
Weaknesses:
Although the numerical investigation is thorough and the model hypotheses are clearly stated, the paper falls a bit short of convincing readers that the emergent behaviors observed in the model are of actual biological relevance.
In the introduction of the revised manuscript, we have extended the discussion about experimental observations of supercellular myosin cables that are believed to drive contractions and active cell intercalation events in early-stage chick embryos. These experiments have been the direct motivation for the model. The model, however, does not aim to be a specific, detailed description of the chick embryo, but instead aims to explore whether a relatively simple and generic mechanism of tension-controlled activation of myosin motors can lead to the observed convergence-extension tissue deformations.
The hypotheses and ingredients of the model are clearly stated (1) junctions are viscoelastic, (2) junctions are active, i.e. they can generate tension through myosin, and (3) positive feedback between tension and myosin. Yet the 3rd hypothesis, which gives its singularity to the model, is not discussed enough, and its relevance to well-known convergent extension systems is not analyzed. In the paper, the positive feedback arises from the decrease of unbinding rate upon tension increase. This "catch bond" hypothesis (note: the term catch-bond only appears once, in the methods) is at the core of all the results obtained in the paper, and this should be clearly stated and discussed. As for the positive feedback, one could also imagine that it is not the unbinding rate that decreases upon tension increase, but rather the myosin recruitment rate that increases. Would that be formally equivalent? The authors should address this point and its biological implications.
The comparison between the model outcome (notably, the T1 transitions preferably oriented in the direction perpendicular to pulling) and existing experimental data is a bit elusive. The authors propose a semi-quantitative comparison with convergent-extension occurring upon primitive streak formation in avians. The data and conclusion is not fully convincing, partly because of the polarization data itself (Figure 8C), and partly because it's not clear what would be the equivalent of the external pulling force in this self-polarizing flow.
The authors also mention the germband of , but the model is, I believe, not relevant for that system. Although there is indeed an external force pulling the tissue (pulling posterior midgut), polarization is perpendicular and extension parallel to the pulling force (the opposite of the model outcome).
We hope that the revised manuscript now makes it clear that there may be external forces as exerted by the Drosophila midgut invagination during germband extension, but the idea for the chick embryo is that the intercalations start somewhere in the tissue (e.g. possibly by localised ingressions of cells in the sickle region). These intercalations exert mechanical stresses on neighbouring cells, which develop anisotropy themselves, contract, and activate their neighbours, leading to a chain-like sequence of intercalations. This is similar to the Drosophila case, where the apical contraction of the mesoderm precursor cells in the ventral midline resulting in their invagination in the ventral furrow, exert a pulling force on the lateral mesoderm. Also, in this case that cables formed are oriented towards the ventral furrow, resulting in contraction of this tissue along the direction of the myosin cables and extension perpendicular to those.
In that spirit, they should also discuss pulling experiments (fly wing disc / Mao lab, suspended epithelia / Charras lab) where it's not so clear that topological transitions occur upon pulling.
In the revised manuscript, we have pointed out that the fly wing disk of the Mao lab and the suspended epithelia of the Charras lab are epithelia likely made up of cells that are in a different more mature differentiation state and may not be able to undergo the same amount of intercalations as the mesendoderm precursor tissues in the chick and the fly. Neither pulling experiments on cultured MDCK epithelial cell monolayers (Harris, et al. 2012) nor on the fly wing disk (Duda, et al. 2019) show any topological transitions even at stains above 50%. In the case of the MDCK monolayers, it is attributed to the slow turnover of keratins and desmosomes (~1h), reflecting their more mature differentiation state. Therefore, while it is possible that topological transitions would occurs over long times, the mechanism presented in this work would not likely apply to stretching experiments of these tissues.
Finally, the paper builds upon already existing models of junction and tissue mechanics. Beyond the vertex model itself, the hypothesis of a feedback between myosin activity and tension has also been around in the field for some time, including in modeling papers (which authors acknowledge honestly).
In conclusion:
The paper presents a thorough numerical analysis of T1s and convergent extension upon tissue pulling when a positive feedback exists between Myosin and junction tension. Yet, the fact that this type of model and hypotheses are not completely new, combined to a lack of comparison to experimental data (except for avian gastrulation) makes it a bit difficult for the reader to assess the novelty of the findings presented in the paper.
I liked the presentation of the model a lot, the details given, the pedagogy, the comparison with what already existed (Discussion section), and the thorough analysis of the model outcome.
On the other hand, I find it a bit difficult to see whether the novelty is strong enough for publication in eLife:
First, because somewhat similar models exist in the field. According to the authors, the specificity here is the inclusion of the positive feedback for the active part of the model. Yet, other models that the authors cite include a similar type of feedback, although they might not specifically focus on topological transitions and convergent extension. In addition, the existence of such a feedback (and its possible consequences on tissue morphogenesis) has been around in the field for a while, under various forms (here, the authors chose a catch-bond type dissociation rate for myosin).
We agree with the reviewer that similar models already exist in the literature. There is growing evidence that many tissue-scale processes are driven by myosin activity that manifests itself as active contractions and expansions of cell-cell junctions that lead to cell rearrangements, tissue ordering, etc. Assuming some form of tension-dependent myosin dynamics (e.g. catch-bonds), it is not surprising that many authors have developed similar models.
It is, however, surprising that models that share common basic principles can capture, at least qualitatively, rather different biological processes in different organisms. This suggests a set of universal physical processes that govern tissue scale behaviours.
In this work, we have applied a version of a model that treats cell-cell junctions as active entities to the previously unexplored problem of T1 transitions that occur against externally applied tension, a process of importance in understanding early development of avian embryos. We, therefore, argue that this works makes a significant contribution to the existing literature of understanding the role of activity in epithelial tissue mechanics. In addition, our computational model is able to treat the whole sequence from the application of external forces to the induced tissue flow through T1 transitions, i.e. the full emergent active rheology, a significant advance compared to other models.
In the revised manuscript, we make the connection to other models clearer.
Second, because the comparison (even qualitative) with existing experimental data remains elusive.
Hence the major finding of this paper is that a tissue in which a positive feedback between tension and motor recruitment exists can display active intercalations and convergent extension perpendicular to an external pulling force. In which experimental systems this prediction could be relevant remains to be determined.
In the revised manuscript, we make this connection clearer by referring to experiments on primitive streak formation in chick embryos. Current state of the art experiments, however, cannot yet follow myosin dynamics in-vivo, and that is why we focus on cell intercalations that, we argue, are a good readout of myosin driven active contractions. We also reiterate that the other main qualitative evidence for the tension feedback state, tissue extension is orthogonal to the main direction of cell elongation, is seen in both model (Figure 7C) and experiment (Figure 8D-E). We emphasise, however, that the aim of this work is to study a generic biologically plausible model for tension-driven activation of cell-cell junctions, and its effects on cellular behaviours.
– Does the positive feedback tension/myosin need to stem from the catch bond model? Could one imagine increased motor recruitment open junction/cell stretch (This would be in line with some experimental data, including stretching experiments). Would that be formally equivalent to the catch bond, or would that change the outcome of the model? If yes, how?
While the proposed model is coarse and does not make any specific assumptions about the molecular mechanisms behind the positive tension/myosin feedback, the catch-bond mechanism would be a plausible molecular realisation for the proposed model. The assumption that the myosin association rate is tension independent and that the dissociation rate decays with tensions is based on experiments of Kovács et al. (PNAS 104, 9994, 2007) on non-muscle myosin II that is believed to be the key molecular motor for active processes studied in this work.
Due to the coarse-grained nature of the model, assuming that motor recruitment depends on tension would not lead to qualitatively different results compared to the current model. The reason is that the central assumption is that external tension leads to accumulation of motors and, hence, increased contractility, without focusing on the specific molecular origin of such build-up.
We have revised the manuscript to make this point clearer.
– In which systems do people observe tissue extension along the direction perpendicular to pulling? The authors only refer to the avian embryo in which the nature of the flow (embryo-scale) makes it a bit difficult (at least for me) to interpret what is the chicken and what is the egg. A lot of data is now available out there for analysis. It would be great to see the same kind of analyses in other experimental systems. Notably, the fly germband (that the authors rapidly mention) seems to behave completely differently, with cables of Myosin perpendicular to the pulling axis (pulling force being exerted by the posterior midgut). How about convergent extension in the zebrafish or Xenopus embryos?
This model has been directly inspired by experiments on avian (chicken, to be precise) embryos, where there has been no clear evidence of chemical pre-patterning. The model analysed in this study proposed a plausible cell-level mechanism for observed cell intercalation events that occur along the direction of the actomyosin cables, believed to be a proxy for the direction of mechanical tension. This mechanism is also applicable to the germband extension in the fly, where the apical contraction of the cells in the ventral midline also exert pulling forces on the lateral mesoderm, resulting in myosin cables driving junction contractions of intercalating cells towards the ventral midline.
We briefly note that direct comparison with Xenopus and the zebrafish would not be possible. In those organisms, the morphogenic flows involve mesenchymal cells (not epithelial cells) and they do exhibit a different mechanism of intercalation movement, i.e. their migrate on top of each other, which is different from junction shortening and lengthening. It still involves polarisation, extension of cellular process and myosin II driven contraction. The detailed mechanism is, however, less clear in this case (Pfister, et al., Development 143.4, 715, 2016).
– What happens at the boundaries? Here we look at a central patch of cells. And somewhat related question, how would that work in a closed geometry (such as the fly embryo). I feel that the avian discussion (paragraph 1 of discussion) could be expanded with considerations on how the feedback can produce long range flows even though only a small patch of cells is initially "activated".
The boundaries could indeed play an important role, especially in the case, as pointed out by the reviewer, of a closed geometry such as the embryo. This is clearly the case both for the fly and the chick embryo, for which in the early stages of the development, the central embryonic area is surrounded by the area opaca (a thicker region of cells attached to the vitelline membrane) that provides tension to the central region. Our attempts to simulate such geometries (not shown) led to unsatisfactory results. The posterior sickle-shaped active region indeed starts to contract and exhibits motion reminiscent of convergence-extension but loses coherence and does not lead to a streak-like structure. We speculate that presence of cell divisions and cell ingressions plays an important role in large scale flows (as suggested by recent experiments Chuai, et al. Science Advances 9.1, eabn5429, 2023), which were not included in the model. One major effect of division and ingression is that the tissue is always a fluid without a yield stress (Matoz-Fernandez, et al. Soft Matter 13, 3205, 2017), which will indeed facilitate activity-induced flows. Such effects will, however, be explored in detail elsewhere. The goal of this study is primarily to showcase that under simple assumption of tension dependent myosin kinetics, it is possible to have cell intercalations that are perpendicular to the direction of external tension.
– Authors should cite Duda et al. (Dev Cell 2019) from the Mao group. In this article it is shown that Myosin polarizes and form cables upon external tissue stretch in the fly wing disc. This is fully relevant to their model and it is somewhat surprising that they do not cite it. Disclaimer: I'm not on that paper
We thank the reviewer for pointing out the Duda, et al. paper. It is indeed very relevant for our work, and it was omitted by from the original submission by mistake. We refer to it in the revised manuscript.
– The authors present their model as "one of the simplest descriptions" that features the three main ingredients required (viscoelastic junctions, active junctions, positive feedback myosin/tension). What would be the simplest and why didn't they choose it?
This is statement has been poorly worded and it is meant to convey that we have chosen a minimal model that is able to describe a cell-cell junction that has a viscoelastic, an elastic and an active component, and acts against the elastic background of the surrounding tissue that is under tension. In the revised manuscript, we have rephrased the sentence to read: “In order to describe the three key ingredients that characterise a cell-cell junction, while taking into account the effects of the elastic background, we adopt a minimal description, i.e. we assume that the junction consists of four elements connected in parallel:”.
– Authors should better explain the rationale of the "elastic barrier" term B(l-a) , as compared to the k(l-l0) term.
The k(l – l0) term describes the viscoelastic properties of the junction itself, which is encoded in the time-dependent value of rest length l0 that relaxes towards the actual junction length. The elastic barrier, on the other hand, refers to the presence of the surrounding tissue that is assumed to be under tension and therefore, its response contains a purely elastic component, which acts against junction remodelling. In the language of soft matter physics, it represents the yield stress of the surrounding tissue. Clearly, there are multiple possible approaches to modelling the effects of such elastic background. We opted for a simplistic description that models the elastic background as the additional elastic spring of stiffness B. While this is a rather crude approximation, it is reasonable to assume that more sophisticated approaches would not qualitatively change the results. In the revised manuscript we have add the following sentences: “Furthermore, epithelial tissues, and early-stage embryos, in particular, are often under mechanical tension. This tension generates an elastic background against which the junction contracts and expands. We refer to that background elasticity as the elastic barrier since it is assumed that it acts to prevent junction remodelling. ”
Reviewer #2 (Recommendations for the authors):
Sknepnek et al. study a model for junction dynamics in a two-dimensional vertex model, intended to describe epithelial mechanics. In the model, a concentration of myosin on each half-junction evolves according to a chemical balance equation and induces an active tension in the junction. The chemical balance equation incorporates a mechanical feedback, whereby myosin unbinding is inhibited at high tension. The authors first study a single junction and consider thresholds for the junction to fully contract. They then investigate the behaviour of a few « active » junctions in a passive network and a fully active vertex model. They find that at intermediate external pulling force and myosin-induced tension, the simulated epithelium can undergo convergence-extension against the pulling force.
The manuscript describes a work which is seriously and rigorously performed. The simulations are analysed thoroughly with a serious quantification effort. The authors have made an important effort of clarity in their study.
We thank the reviewer for insightful comments and suggestions. We address each point below.
My central question is about the motivation of the manuscript. The central theme of the manuscript authors study « contrarian » active T1 transitions which induce convergence extension in the direction orthogonal to an externally applied stress. I would disagree with the authors that this is the only type of active T1 transition. For instance polarised transitions which are oriented by the epithelium chemical polarity and induce converge-extension in the absence of external tension, as is thought to be the case during germ-band elongation, could also be called « active ». From the manuscript it is not clear to me why inducing convergence extension as a response to external tension, against the external tension, is important or biologically relevant. I think the authors should significantly clarify the rationale for their study.
We agree with the reviewer that T1 transitions studied here are not the only possible type of active T1 transitions. In the revised manuscript, we have added the following sentences: “The focus of this work is on the T1 transition that occurs perpendicular to the direction of the maximum principal mechanical stress. We refer to such T1 transition as active, and remark that other types of active T1 transition in embryonic tissues are possible, e.g. polarised transitions oriented by the epithelium chemical polarity that induce converge-extension in the absence of external tension, as is thought to be the case during germ-band elongation (Bertet, et al. (2004)).”
In that respect I am not sure that the comparison with experiment is very useful. My understanding is that essentially myosin cables in the chick embryos appear oriented along the direction of tissue contraction, which seems to make sense if these cables are actively contracting. I am not sure that this observation strongly supports the much more detailed model proposed by the authors; notably the notion that myosin polarises in response to external tension.
The main idea is that some cells start to contract and then pull on neighbouring cells as discussed in the revised text, resulting in myosin accumulation and building up of supercellular chains that then, most likely, direct the intercalation chains. The chains of intercalations that we see in the full vertex simulations have been shown to exist in the experiments.
Judging from Figure 6, the authors also observe relatively limited convergence extension happening against the external force. I wonder if this may be due to myosin molecules playing both the role of the force-exerting molecule and the sensor in the model? For instance, possibly if a secondary molecule responding to external tension, and itself inducing myosin polarisation, a more permanent convergence-extension could be induced?
The reviewer raises and interesting and important question, which, unfortunately, cannot be directly answered with the current model. We speculate that a more prominent convergence extension would require spatial and temporal coordination of active T1 events. Such coordination does not seem to arise in the current model, and we speculate that it is possibly related to the effects the reviewer has alluded to. For example, in the current model, the post-T1 expansion of the newly formed junction is effectively passive. Additionally, our model omits divisions and ingressions which necessarily induce flow in model tissues (Matoz-Fernandez, et al. Soft Matter 13, 3205, 2017) and could help break “stuck” myosin cables that halt convergence-extension flows in the simulations.
To provide a detailed answer to reviewer’s question, however, one would need to extend the part of the model that describes kinetics of myosin motors to include additional species of molecules. Such efforts have been made, albeit in a different context, e.g. by Lan, et al. (2015). It would also likely be necessary to modify the coupling between myosin kinetics and mechanical tensions. For example, recent studies by Krajnc, et al. (2021) and PérezVerdugo and Banerjee (arXiv preprint, Nov 2022) suggest that tensions fluctuations have interesting effects on the cell intercalation dynamics, among other things leading to a nonintuitive observation of enhanced occurrences of T1 transitions with increasing mechanical tension. Extending the model to include such effects, while very interesting, is beyond the scope of this work, and will be addressed elsewhere.
– I think the authors may have overlooked the reference « active instability and nonlinear dynamics of cell-cell junctions », Krajnc, et al, PRL, 2021. The approach and questions seem very close to the current manuscript. I think the results and approach of this study should be contrasted to the current work.
We thank the reviewer for pointing out the paper by Krajnc, et al. 2021. In the revised manuscript, we cite it and comment on its main conclusions.
– When comparing Figure 7C and Figure 6, I had the impression that Figure 7C showed significant « negative » convergence extension for large values of β (0.75), why the authors say that this is not the case in the caption of Figure 6?
The reason convergence extension in Figure 6 appears less significant compared to the measurements shown in Figure 7C, is due visualisation in Figure 6. In simulations, the activity was first switched off and the system was first stretched until it reached the mechanical equilibrium. The activity was then switched on and the system was allowed to actively contract. The amount of convergence-extension was measured relative to the initial stretched state when the activity was first switched on. In order to make this clear, in the revised manuscript, we have included a short supplement video to Figure 6 that shows the initial passive stretching and subsequent active contractions for the same parameters used in Figure 6.
– The model assumes a constant number of available myosin molecules in each cell. It was not clear to me if this assumption is important for the results of the model?
Conservation of myosin per cell is introduced in order provide a mechanism that would allow coupling of myosin distributions on neighbouring junctions. In other words, the assumption is that junctions adjacent to a given junction act to deplete its myosin, which is an important assumption required for the proposed mechanism to drive an active T1 transition. Clearly, there are many ways to model such mechanism and we explored several possible choices. We settled with the conserved myosin model since it has only one free parameter (i.e. the total myosin per cell), while being plausible. In the revised manuscript, we have added the following sentence to make this point clear: “The finite pool of myosin acts to introduce correlations in the distribution of junctional myosin within each cell. In other words, it models the mechanism by which, e.g. depletion of myosin on a given junction is correlated to myosin accumulation on its neighbouring junctions. This is of central importance to establish myosin anisotropy within the cell. There are clearly many alternative ways to model coupling of myosin on different junctions within a cell. The model used here, however, requires a minimal number of parameters while being biologically plausible.”
– I do not understand the reason for the factor « z » in the second term of the right-hand side of Equation 5.
The factor “z” in Equation 5 is the number of junctions of a given cell. It is there because of the assumption of conserved myosin in each cell as discussed in the answer to the reviewer’s previous question. Its origin is explained in the discussion leading to Equation 12. Its presence, however, does not qualitatively affect the behaviour of the model, and even in the case of a random tissue patch, without the loss of generality, one can set z=M, where M is the total myosin pool in the cell.
Reviewer #3 (Recommendations for the authors):
The authors propose a mechanism for the dynamics of cell-cell junctions under an applied external mechanical load. In their theoretical analysis they assume that myosin recruitment and therefore active contractile stress in a junction increases with an applied external load. They find that this mechanism can induce a T1 transition through which adjacent cells exchange neighbours such that an originally aligned junction is reorientated perpendicularly to an applied uniaxial extensile stress. By using a generalised vertex model they also find that this effect on an isolated cell-cell junction can lead to convergence-extension of a tissue, where the extension occurs perpendicularly to the applied uniaxial extensile stress. In both cases, the phenomenon is present for a broad range of parameter values. Eventually, the authors compare their theoretical results to convergence extension of a part of a tissue in an early chick embryo.
We thank the reviewer for critical reading of our manuscript and helpful and insightful comments and suggestions.
The theoretical analysis is carried out mostly through simulations and is done in a convincing way. The dynamics assumed for a single junction is based mostly on reasonable assumptions and differs from earlier works. It appears though that the assumption of a monotonic decrease of the myosin unbinding rate with the applied force needs more justification. I would assume that for sufficiently high forces, this rate will increase as is the case for catch bonds. Taking such an effect into account might limit the range of parameter values for which stress-induced convergence-extension is observed. Also, it was unclear to me whether all ingredients of the model where necessary to obtain the observed effect. In particular, I wonder what happens if either of the dashpots or the springs are eliminated? Finally, I do not understand why m_0 should describe effects of the surrounding. Shouldn't this effect be included in T_ext? Otherwise the theoretical analysis is convincing and justifies the conclusions of the authors.
As shown in Figure 4A, there is indeed a finite range of parameters over which the convergence-extensions occurs. This is because at very high applied tensions, passive T1 transitions in the direction of pulling additionally appear. We agree that effects such as the increase of the detachment rate at very high tensions would indeed affect the range of convergence-extension, by contributing to another increase in passive T1 transitions at very high stresses. We argue, however, that the parameter range favourable for convergence extension is still sufficiently wide to be biologically plausible. Additionally, the tissue eventually extending in the pulling direction at extreme forces is also realistic.
Regarding the other ingredients, in the revised manuscript, we explain the motivation behind and the significance of each element in the model. We argue that this is the smallest number of elements that are necessary to observe tension-driven activation and subsequent convergence-extension flows. Therefore, eliminating any of the elements would lead to qualitative changes in its behaviour. The fact that we need a second spring, m0 (i.e. the myosin reference level) and T* (i.e. the threshold tension) is a peculiarity of the vertex model. A large number of iterations over the details of the vertex model showed that the perimeter term (responsible for the second spring) is necessary to retain polygonal cell shapes. This induces self-stresses in the tissue, requiring a T* so that junctions do not randomly activate. Equally, for the active parts if we, e.g. implement a shift by setting both T* and m0 to 0, individual junctions in the tissue activate spontaneously at no applied force. Results from a continuum model (to be published separately) indicate that in the absence of the local self-stresses induced by a disordered vertex model, T* and the second spring can indeed be omitted.
I am less convinced by the comparison with experiments on gastrulation in early-stage chick embryos. Although the relative orientation between myosin cables and the direction of extension agrees with that in the theory, I fail to see where the external uniaxial extensile stress should come from. It thus appears that an essential element of the theory is not checked for in the experiments. The authors write themselves in ll 395: "Although it is yet to be experimentally confirmed, it is plausible that the symmetry breaking event that induces the initial myosin polarity occurs as a result of anisotropic tension combined with cell differentiation early in development." The authors should clarify this point, which seems to be central to the authors' idea that mechanical signalling drives convergence-extension in this case. Otherwise, the experimental data do not add much to the work.
In the revised manuscript, we make the connection with the experiments clearer. Namely, we point out that mechanism does not need specific external forces, and the pulling forces can be generated within the tissue as described several times in the text and above. In order to avoid the complicated problem of the origin of pulling forces in real biological system (which is also beyond the scope of this work), in the model we used external forces solely for the purpose of generating the initial anisotropy that activates the system. We have also added an example of a central myosin pulse triggering a chain contraction in a one-dimensional system (new Figure 1 – supplement figure 2).
The authors frequently use the term "polarisation" or "polarised" when referring to the alignment of the myosin cables, for example. This notion implies directionality, for example, a difference between left and right in Figure 2, which is not present. I would suggest that the authors rather use 'anisotropy'.
We agree with the reviewer that the term “polarisation” might be confusing. It indeed is meant to mean “anisotropy”. In the revised manuscript, we have used the term “anisotropy” instead.
The text exhibits some jumps between figures that are not always referred to in order, for example between Figures2 and 3 or Figure 4C and 4B. The authors might want to adapt their text to the figures or vice versa.
We have updated Figure 2 to have only two panels (A and B). We have swapped panels B and C in Figure 4, so they are no referenced in order in the text.
In Figure 1A the term \β m is hard to see. Please, improve.
We have updated Figure 1A such that βm is clearly visible.
In Figure 3 the colour code is confusing. In A black, red, blue correspond to different junction classes, in B and C to different mechanical quantities and in D again to different junction classes but not the ones in A. Please, improve.
We have modified colours in Figure 3 and Figure 2B —figure supplement, and updated figure captions, such that there is no longer ambiguity about what each colour represents.
Ll 49 "With cellular behaviours being coordinated over thousands of cells in the case of the chick embryo, biochemical signalling alone is unlikely to account for the observed motion patterns." Even though I have a guess, I do not really understand this argument. Can you add some words to explain, why the long-range coordination is at odds with pure biochemical signalling?
Long range biochemical signalling requires either setting up and maintaining gradients of signalling molecules, which is a slow process over larger distances. These gradients would be deformed continuously by the tissue flows making it difficult to see how this would be maintained. Alternatively, some kind of small molecule cell to cell signal relay mechanism, such as the cAMP relay mechanism that controls the chemotactic aggregation of Dictyostelium cells, could be present. However, so far, no evidence for any kind of long-range relay mechanism exists. It seems much more logical that stresses generated locally due to specific cell behaviours, and which spread rapidly and over longer distances into surrounding tissues, influence directly biochemical signalling via mechanosensitive biochemical reactions. These have now been abundantly demonstrated to exist, especially mechanical signalling to the cytoskeleton that directly affect cellular function and behaviours such as adhesion and motility.
L 218 side a -> side OF a
Fixed.
L 256 along the direction the central junction -> along the direction OF the central junction
Fixed.
Ll 274 could you add a panel to show this?
The behaviour is very similar to what is shown in Figure 7 —figure supplement 3B for the random patch. The regime where this happens is not biologically relevant and the model tissue made of hexagonal cell only is rather artificial. Therefore, we do not believe that adding a panel would be warranted. In the revised manuscript, we, however, alert the reader to the discussion in the section devoted to a fully random active patch.
In Figures6 and 7A, convergence-extension is not really visible assuming that the initial state was square. Maybe you can show a sequence similar to Figure 2?
We agree with the reviewer that the static images do convey convincingly convey the extent of convergence-extension. To address this, in the revised manuscript, we added a video supplement to Figure 6 that shows the initial passive stretching and subsequent active contractions, thus showing the full convergence-extension process.
[Editors' note: further revisions were suggested prior to acceptance, as described below.]
The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:
Your model doesn't seem to apply to the fly germband elongation, in which the invaginating midgut (pulling external force) and the T1s occur along the same axis, with cables perpendicular to the pulling axis. In your response, you insist that the pulling force is rather the ventral invagination than the posterior midgut (hence a force perpendicular to tissue elongation and T1s, and parallel to the cables, as in their model). This overlooks that the germband does not extend without the posterior midgut pulling (Torso), while it does without the ventral mesoderm invagination. In the paper, claims about their finding applying to germband extension are found in the intro and discussion, without other justification. These comments should be removed unless very strongly justified, as it seems that they are based on a direct misinterpretation of previous observations.
We thank the reviewers and editors for this comment. We agree that extrinsic forces associated with midgut invagination are involved in germband extension and that in the torso mutant there is no germband extension. However, in the original experiments by Irvine and Wieschaus (1994) [doi.org/10.1242/dev.120.4.827] they state that although there is no posterior movement, cell intercalations do not stop. More recent experiments by Collinet, et al. (2015) [doi.org/10.1038/ncb3226] show that both in Torso mutants as well as in embryos where a DV line of tissue in the posterior of the embryo is cauterised to the eggshell, cell intercalations and junctional myosin accumulation continue, showing that intercalations are either autonomous or dependent on other external forces. Furthermore, the same work shows via junction ablation experiments that the tension in the DV junctions is higher than in AP junctions both in wildtype as well as in Torso mutants. The DV/AP aspect ratio of the cells is initially 1.4 in the wildtype decreases to smaller than 1 over time, while in Torso mutants it stays close to 1.2. This would be in line with the observation that at least initially the DV tension dominates over the AP tension. More recently experiments by Gustafson, et al. (2022) [doi.org/10.1038/s41467-022-34518-9] have shown that the strain rate shows a DV gradient with highest value that ventral side and this correlates strongly with the junctional myosin association rate. Furthermore, they also showed that contraction of cells using optogenetic techniques can induce deformation and myosin accumulation in neighbouring cells. In twist and snail mutants where no ventral furrow forms, germband extension is greatly slowed down, again hinting at the importance of that process in driving intercalations. Furthermore, Lefebvre et al., (2023) [doi.org/10.7554/eLife.78787] have shown that during germband extension the orientation of the myosin anisotropy remains primarily DV oriented, while the pair rule gene expression domains show a high vorticity, suggesting an indirect or complex coupling of gene expression to myosin localisation. Therefore, we think that the type of feedback we propose in the model, together with morphogenetic patterning, is relevant to the process of germband extension. We have rewritten the section on Drosophila germband extension to make this more explicit. The updated text is shown in magenta.
When the germband is mentioned you usually cite a paper that does NOT deal with the germband (Jacinto et al. 2002, which is about dorsal closure). And you also cite Duda et al., even though it's a wing disc paper, in which cables form parallel to the pulling force, (as in their model, but in contrast with the germband). This could be confusing and even detrimental to your work when readers familiar with fly morphogenesis read it.
We completely agree with these comments that this was not clearly explained. The Jacinto paper was cited as another case where supercellular myosin cables were observed, while the Duda paper was cited as another case where it as been shown that tension can result in myosin recruitment. These observations, although relevant, did not relate to germband extension and we have removed them from the revised version.
https://doi.org/10.7554/eLife.79862.sa2Article and author information
Author details
Funding
Biotechnology and Biological Sciences Research Council (BB/N009789/1)
- Rastko Sknepnek
- Manli Chuai
- Cornelis Weijer
Biotechnology and Biological Sciences Research Council (BB/N009150/1-2)
- Ilyas Djafer-Cherif
- Silke Henkes
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
RS, MC, and CJW acknowledge support by the UK BBSRC (award BB/N009789/1). SH and IDC acknowledge support by the UK BBSRC (grant number BB/N009150/1-2). IDC acknowledges funding under Dioscuri, a programme initiated by the Max Planck Society, jointly managed with the National Science Centre in Poland, and mutually funded by the Polish Ministry of Science and Higher Education and German Federal Ministry of Education and Research (UMO-2019/02/H/NZ6/00003). We thank Antti Karjalainen for providing the MATLAB code used to analyse experimental data. RS thanks Andrej Košmrlj, Daniel Matoz-Fernandez, and Sijie Tong for many helpful discussions about the vertex model.
Senior Editor
- Aleksandra M Walczak, CNRS, France
Reviewing Editor
- Karsten Kruse, University of Geneva, Switzerland
Version history
- Preprint posted: June 23, 2021 (view preprint)
- Received: April 29, 2022
- Accepted: April 9, 2023
- Accepted Manuscript published: April 11, 2023 (version 1)
- Accepted Manuscript updated: April 24, 2023 (version 2)
- Version of Record published: May 3, 2023 (version 3)
- Version of Record updated: May 4, 2023 (version 4)
Copyright
© 2023, Sknepnek et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 739
- Page views
-
- 122
- Downloads
-
- 3
- Citations
Article citation count generated by polling the highest count across the following sources: PubMed Central, Crossref, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
- Physics of Living Systems
Olfactory navigation is observed across species and plays a crucial role in locating resources for survival. In the laboratory, understanding the behavioral strategies and neural circuits underlying odor-taxis requires a detailed understanding of the animal’s sensory environment. For small model organisms like Caenorhabditis elegans and larval Drosophila melanogaster, controlling and measuring the odor environment experienced by the animal can be challenging, especially for airborne odors, which are subject to subtle effects from airflow, temperature variation, and from the odor’s adhesion, adsorption, or reemission. Here, we present a method to control and measure airborne odor concentration in an arena compatible with an agar substrate. Our method allows continuous controlling and monitoring of the odor profile while imaging animal behavior. We construct stationary chemical landscapes in an odor flow chamber through spatially patterned odorized air. The odor concentration is measured with a spatially distributed array of digital gas sensors. Careful placement of the sensors allows the odor concentration across the arena to be continuously inferred in space and monitored through time. We use this approach to measure the odor concentration that each animal experiences as it undergoes chemotaxis behavior and report chemotaxis strategies for C. elegans and D. melanogaster larvae populations as they navigate spatial odor landscapes.
-
- Physics of Living Systems
- Structural Biology and Molecular Biophysics
Non-membrane-bound biomolecular condensates have been proposed to represent an important mode of subcellular organization in diverse biological settings. However, the fundamental principles governing the spatial organization and dynamics of condensates at the atomistic level remain unclear. The Saccharomyces cerevisiae Lge1 protein is required for histone H2B ubiquitination and its N-terminal intrinsically disordered fragment (Lge11-80) undergoes robust phase separation. This study connects single- and multi-chain all-atom molecular dynamics simulations of Lge11-80 with the in vitro behavior of Lge11-80 condensates. Analysis of modeled protein-protein interactions elucidates the key determinants of Lge11-80 condensate formation and links configurational entropy, valency, and compactness of proteins inside the condensates. A newly derived analytical formalism, related to colloid fractal cluster formation, describes condensate architecture across length scales as a function of protein valency and compactness. In particular, the formalism provides an atomistically resolved model of Lge11-80 condensates on the scale of hundreds of nanometers starting from individual protein conformers captured in simulations. The simulation-derived fractal dimensions of condensates of Lge11-80 and its mutants agree with their in vitro morphologies. The presented framework enables a multiscale description of biomolecular condensates and embeds their study in a wider context of colloid self-organization.