Abstract
Morphogenetic processes often involve the rapid rearrangement of cells held together by mutual adhesion. The dynamic nature of this adhesion endows tissues with liquid-like properties, such that large-scale shape changes appear as tissue flows. Generally, the resistance to flow (tissue viscosity) is expected to depend on the cohesion of a tissue (how strongly its cells adhere to each other), but the exact relationship between these parameters is not known. Here, we analyse the link between cohesion and viscosity to uncover basic mechanical principles of cell rearrangement. We show that for vertebrate and invertebrate tissues, viscosity varies in proportion to cohesion over a 200-fold range of values. We demonstrate that this proportionality is predicted by a cell-based model of tissue viscosity. To do so, we analyse cell adhesion in Xenopus embryonic tissues and determine a number of parameters, including tissue surface tension (as a measure of cohesion), cell contact fluctuation and cortical tension. In the tissues studied, the ratio of surface tension to viscosity, which has the dimension of a velocity, is 1.8 µm/min. This characteristic velocity reflects the rate of cell-cell boundary contraction during rearrangement, and sets a limit to rearrangement rates. Moreover, we propose that, in these tissues, cell movement is maximally efficient. Our approach to cell rearrangement mechanics links adhesion to the resistance of a tissue to plastic deformation, identifies the characteristic velocity of the process, and provides a basis for the comparison of tissues with mechanical properties that may vary by orders of magnitude.
INTRODUCTION
Among multicellular organisms, animals are unique in their ability to perform rapid cell rearrangements based on a combination of dynamic, flexible cell adhesion and the modulation of the cytoskeleton at cell contacts. This feature is essential for many morphogenetic processes, and is extensively studied at the subcellular level during epithelial cell rearrangement. It also underlies the liquid-like behaviour of tissues (Gordon et al., 1972; Forgacs et al., 1998; Marmottant et al., 2009), which invites the analysis of cell rearrangement at a more generalised, mechanical level.
A core mechanical parameter of a liquid-like tissue is its cohesion: the work required to separate the tissue into two parts by breaking cell adhesions that hold it together. In analogy to liquids, cohesion can be measured by its effect at the tissue surface, as surface tension (Foty and Steinberg, 2005; Manning et al., 2010) (see Box 1). While cohesion ensures tissue integrity, it also confers resistance to large-scale movement. This resistance is measured as tissue viscosity (Gordon et al., 1972; Marmottant et al., 2009; Rieu and Sawada, 2002; Mombach et al., 2005). To understand the basic mechanics of cell rearrangements in liquid-like tissues, the relationship between tissue cohesion and viscosity needs to be clarified.
Box 1. Cell adhesion and cortical tensions
In an isolated cell, cortical tension minimises the cell surface by rounding it up into a sphere. For cells to attach to each other, cortical tension β at the exposed surfaces of two adjacent cells has to be balanced by a reduced tension β* in each cell at the area of contact, as shown in the drawing. For this equilibrium, β*/β=cosθ, where θ is the contact angle between the cells. Thus, without any reduction of tension at contacts, i.e. β*=β, the contact angle is zero and cells touch at one point only. By contrast, complete reduction to β*=0 gives θ=90°, and two cells attach maximally. These cellular tensions are linked to a tissue-level parameter – the tissue surface tension σ – by the approximation σ≈β−β* (Brodland et al., 2009; Manning et al., 2010) (see supplementary mathematical analyses, Section 1). This allows us to calculate β=σ/(1−cosθ) and β*=β cosθ from measurements of σ and θ.
Cell aggregation is mediated by adhesion molecules such as cadherins. Cadherin engagement has two effects. First, release of the binding energy Γ favours adhesion by reducing tension at contacts. Second, tension in the cortical cytoskeleton is also reduced, to a degree that depends on cadherin density (see main text). The total reduced tension amounts to β*=βc−Γ/2 in each cell, in which βc is a cell cortical tension at a contact (Manning et al. 2010). Importantly, for cadherins the binding energy Γ is small compared with cortical tensions (Maître et al., 2012; Youssef et al., 2011) (supplementary mathematical analyses, Section 2), and β* corresponds practically to the reduced cortical tension βc.
With binding energy Γ being negligible, the strength of adhesion, or adhesion energy, is determined by the exposed surface cortical tension of a cell, and by the degree of its reduction at contacts. Tensions β and β*, as well as surface tension σ, all have the dimension of [energy/area], and separating two cells corresponds to generating new surface with surface energy β from previous surface with energy β*, which requires an energy per unit area of new surface of σ≈β−β*. Thus, tissue surface tension is a measure of the cell-cell adhesion energy per area or, equivalently, of tissue cohesion. The value for cohesion or adhesion energy is twice the respective surface tension value, as two surfaces, each with a surface energy corresponding to surface tension, are generated when a tissue is split into two. However, we may use these terms interchangeably when the difference is not relevant, and will colloquially refer to adhesion energy as the strength of cell adhesion.
Given the importance of tension reduction at cell contacts, we define a relative adhesiveness α as the reduction in tension at contacting interfaces, relative to exposed surface cortical tension, α=(β−β*)/β. For β=β* there is essentially no adhesion, and to the degree that Γ is negligible, α≈0; if β*=0 then maximal adhesion ensues and α=1. With α=(β−β*)/β=σ/β, the strength of adhesion can be expressed as σ=αβ. Moreover, as α=(β−β*)/β=1−(β*/β)=1−cosθ, the relative adhesiveness α depends only on the contact angle, and can be simply determined from its measurement.
As a first step in this analysis, we show that for different tissues of the Xenopus gastrula, and for other tissues from various species, these two parameters are directly proportional to each other. We then determine a number of additional cellular parameters, in part derived from an analysis of cell adhesion, and use these to show that a recent model of tissue viscosity (Marmottant et al., 2009) predicts and explains the linear relationship between tissue surface tension and viscosity. These findings have far-reaching implications for the mechanics of cell rearrangement. We argue that rearrangements are characterised by a constant rate, and that tissues are tuned for maximally efficient cell movement over a wide range of cohesiveness.
RESULTS
Relationship between tissue surface tension and tissue viscosity
We quantified cohesion and viscosity for various regions of the Xenopus gastrula. Tissue explants, after having rounded up under the influence of surface tension, deviate from a spherical equilibrium shape due to gravity. We used this effect to determine tissue surface tension σ (Del Rio and Neumann, 1997; David et al., 2009; Luu et al., 2011) as a measure of cohesion (Box 1, Fig. 1A, Table 1; supplementary material Fig. S2, Table S1 and mathematical analyses, Section 1).
Tissue surface tension and viscosity. (A) Scanning electron microscopy image of gastrula, indicating the regions explanted for measurements. LEM, leading edge mesendoderm; PCM, prechordal mesoderm; CM, chordamesoderm; endo, endoderm; ecto, ectoderm. Surface tension values are given (mJ/m2) with s.d. (B-E) Tissue viscosity measurement. (B) Ectodermal explant at beginning of rounding (left) and 4 h later. Major and minor axes of fitted ellipses are indicated (red lines). (C) Cell rearrangement during explant rounding. Intercalating cells are outlined in the high-magnification frames to the right. Arrows indicate the long axis of the explant. (D) Timecourse of the axis ratio (blue) of LEM explant fits to a dying exponential (red). (E) Proportionality of surface tension and viscosity (tissues as in A). Bars indicate 95% confidence intervals. Slope of the regression line (red), 1.60 µm/min; correlation coefficient, 0.9988. (F) Log-log plot of surface tensions and viscosities of tissues from E and from the literature (supplementary material Table S2). Slope of the regression line (red), 1.78 µm/min; correlation coefficient, 0.9923. (G) Experimentally treated tissues can be off the regression line (red) taken from E.
Tissue parameters
To measure tissue viscosity η, we analysed the timecourse of surface tension-driven, off-equilibrium shape changes of explants. (a list of the symbols used for the main parameters can be found in supplementary material Table S4). A commonly employed method, the fusion of two explants (Gordon et al., 1972), underestimated viscosity. Explants not only merged by flowing into each other as expected, but also extended protrusions across the space between them and zippered up (supplementary material Fig. S1 and mathematical analyses, Section 2). We therefore used the rounding of explants instead (Gordon et al., 1972) (Fig. 1B; supplementary mathematical analyses, Section 2). Observation started about 20 min after excision, when explants had taken on ellipsoidal shapes, and elastic (von Dassow and Davidson, 2009; Luu et al., 2011) and wounding (Li et al., 2013) reactions had subsided. Rounding was associated with cell rearrangement (Fig. 1C; supplementary material Movie 1) and was followed for over 3 h (supplementary material Movie 2). Surface tension was stable over this period (supplementary material Fig. S2, Movie 3), and evolving explant shapes fit well to the rigorous solution for a viscous liquid (Fig. 1D). Tissue viscosities η were calculated from the ratio σ/η obtained from the rounding rates for individual explants and the average surface tension σ of the examined tissue, since a parallel determination of σ for each explant was technically not feasible (Fig. 1E, Table 1; supplementary material Table S1). This slightly overestimates viscosities.
For Xenopus gastrula tissues, surface tension and viscosity are directly proportional over a 10-fold range of values (Fig. 1E). Even more strikingly, when published data from various other tissues are included in the analysis, surface tension and viscosity each vary by three orders of magnitude, and proportionality between the two parameters holds for viscosity values ranging between 2 and 200 kPa·s (Fig. 1F; supplementary material Table S2). The ratio σ/η has the dimension of a velocity and, in the proportional range, for tissues from cnidarians, fish, frog, chicken and mouse, this characteristic velocity ν is νc=1.8±0.9 µm/min. We will refer to the subclass of liquid-like tissues thus defined as νc tissues. The value of νc is reminiscent of common cell migration velocities; we will return to this observation below.
Intuitively, an increase of viscosity in proportion with cohesion might be expected. However, this is not a mechanical necessity: first, above 200 kPa·s viscosity increases without an increase in surface tension (Fig. 1F); second, viscosity and surface tension can easily be decoupled (Fig. 1G). C-cadherin is essential for gastrula cell adhesion (Heasman et al., 1994), and C-cadherin morpholino antisense oligonucleotides (C-cad-MOs) (Ninomiya et al., 2012) reduced surface tension but, counterintuitively, increased viscosity. Expression of a mutant protocadherin, M-PAPC (Kim et al., 1998), or of a cytoskeletal regulator, the kinase Pak1 (Bokoch, 2003), did not change surface tension, but altered viscosity. Thus, proportionality between cohesion and viscosity is not a trivial feature, but appears to be maintained in νc tissues. To elucidate the link between these two mechanical properties, we will derive viscosity from cell-level parameters. For this purpose, we will first analyse cell adhesion, in particular its dependence on cell cortical tension.
Cadherin-dependent cell adhesion
In recent studies, an intimate link between cell adhesion and cell cortical tension became apparent (Box 1) (Brodland et al., 2009; Manning et al., 2010; Youssef et al., 2011; Maître et al., 2012). The cytoskeleton of a cell cortex is usually under tension, which causes isolated cells to take on a spherical shape. Adhesion molecules such as cadherins link cells together, but, as it turned out, the binding energy released upon cadherin engagement is insufficient to overcome the rounding effect of cortical tension. Instead, the cortical tension β itself has to be reduced at contact sites in order for cells to spread on each other (Box 1; supplementary mathematical analyses, Section 3). This implies that the strength of adhesion is mostly determined not by the cadherin binding energy, but by the magnitude of basic cortical tension of a cell, and by how much this tension is reduced at cell-cell contacts (Box 1).
Cell cortex tensions are also related to tissue surface tension (Brodland et al., 2009; Manning et al., 2010), and we calculated cortical tension β and the reduced tension at contacts β* from the surface tension σ and the contact angle between adjacent cells (Box 1). Moreover, we defined a relative adhesiveness α to write surface tension as σ=αβ: surface tension, which represents tissue cohesion, is the product of a dimensionless factor α and the cortical tension at free cell surfaces. α varies between 0 and 1 and expresses the relative reduction of cortical tension at cell contacts (Box 1). Thus, two components contribute in principle to the variation of surface tension between tissues: surface tension is increased or decreased by an increase or decrease in cortical tension β or in relative adhesiveness α.
Relative adhesiveness α only depends on the contact angle between cells and is easily measured (Box 1). In normal gastrula tissues, α is 0.6-0.8 (Fig. 2A, Table 1), indicating a fairly constant reduction of tensions at contacts to about a quarter of exposed-surface cortical tension. As expected from α being constant between tissues, cortical tension β increased with increasing surface tension (Table 1; supplementary material Table S1). Differences between the surface tensions of normal tissues were thus primarily due to differences in cortical tension. However, variation of relative adhesiveness α could be induced by altering cadherin expression. C-cad-MO injection diminished surface tension (supplementary material Fig. S3) by decreasing relative adhesiveness while cortical tension remained constant (Table 1, Fig. 2B). In endoderm, lower cadherin density, as compared with ectoderm, was also accompanied by lower relative adhesiveness (supplementary material Fig. S3; Table 1), but overall the modulation of cadherin-dependent relative adhesiveness was less important in normal tissues than changes in cortical tension.
Cell adhesion and fluctuations. (A) Relative adhesiveness α as a function of surface tension of gastrula tissues. (B) Relative adhesiveness α as function of cadherin density in normal or C-cad-MO-injected (3 and 10 ng/embryo) ectoderm from early gastrula. Averages and s.e.m. from three experiments. (C) Ectoderm cells expressing LifeAct (green) during contact formation. Yellow arrowheads indicate actin cortex at contact; red arrowheads indicate exposed surface cortex. (D) Ectoderm cell pairs more than 20 min after contact formation, with randomly varying small (left) or large (right) contact angle, for comparison of phalloidin-stained F-actin cortex density at exposed (red arrowheads) and contact (yellow arrowheads) surfaces. (E) Ratio of cortex density at contacts (half the measured density, as both cells contribute) versus free surface as a function of relative adhesiveness α=1−cosθ. Cell pairs as in D from uninjected, C-cad-MO-injected (10 ng/embryo) and C-cadherin mRNA-injected ectoderm from early gastrula are compared. The average ratio of cytoplasmic to exposed cortex F-actin intensity is 0.42; a decrease of contact-to-exposed ratios below this level is difficult to detect. (F) Cell-cell contact length fluctuation (red in explanatory diagram) in an ectoderm cell. The inverse of the time between peaks (asterisks) corresponds to the fluctuation frequency. (G) Spontaneous T1 cell rearrangement in an explant of prechordal mesoderm, viewed at incident light. Two cells labelled by red dots separate, allowing cells with yellow dots to establish new contact. (H) Model showing cortical density changes during cell adhesion, cell detachment and in aggregates, and cell detachment due to myosin-dependent boundary contraction.
Cortical tension is actomyosin dependent (Pasternak et al., 1989; Tinevez et al., 2009; Stewart et al., 2011) and is regulated by myosin activity (Maître et al., 2012) and the amount of F-actin at cell contacts (Krieg et al., 2008; Yamada and Nelson, 2007; Hidalgo-Carcedo et al., 2011). In ectoderm cells in the process of attachment, the actin cortex was diminished within minutes of contact (Fig. 2C). In already established cell pairs, the ratio of F-actin at contacting versus exposed surfaces decreased with increasing, cadherin-dependent adhesiveness α (Fig. 2D,E). Our results suggest that cadherin reduces tension at cell contacts through a downregulation of cortex thickness or density (Clark et al., 2013), with the residual cortex being sufficient for cadherin anchoring (Maître et al., 2012). Such a catalytic role of cadherins in adhesion can be mediated by α-catenin (Stirbat et al., 2013). It is the reduced tension β* at contacts that becomes relevant in the following section for the calculation of tissue viscosity (see Fig. 2H).
Viscosity as a function of tissue surface tension
Marmottant et al. (2009) modelled tissue flow by starting from the plastic flow of foams. In a foam, the liquid walls surrounding the air-filled cells minimise their surface due to liquid-air interfacial tension. A small stress leads to a small elastic deformation of the foam; it returns to its initial state when the stress is relieved. However, with increasing stress more and more air cells begin to rearrange by the concordant shrinkage and disappearance of existing walls and the expansion of new walls, a process termed T1 transition (Box 2). At these transitions, energy is dissipated and the non-reversible rearrangements amount to a plastic flow. The resistance to this flow is determined by the energy barrier for T1 transitions. In analogy, tissue flow was treated as occurring through T1 cell rearrangements, with the energy barrier being proportional to the tension at cell-cell contacts and the contact area (Box 2).
Box 2. Tissue viscosity
Marmottant et al. (2009) describe the viscosity η of a tissue due to cell rearrangement as:
We substituted into this equation parameters that we could determine experimentally. ΔεT1 is the energy barrier for a T1 event (see drawing), and is given by the tension at cell-cell contacts multiplied by the contact area (Marmottant et al. 2009). As shown in Box 1, the tension at contacts is 2β* (tensions of both cells added) and with σ=β−β* to replace β* we obtain 2β*=2(β−σ). Regarding contact area, the surface of a cell increases about 1.1-fold when the initially spherical cell is squeezed into tight contact with neighbours to assume an approximately dodecahedral shape of the same volume. Moreover, the number of contacts per cell in a random foam of equally sized cells is 13.4 (Weaire and Rivier, 2009), and we approximated contact area per cell pair amax by dividing total cell surface area by 13.4. From the radii of dissociated cells, we calculated the cell surface area of tightly packed cells, and from this we estimated the average contact area under maximal contact (no gaps) as amax≈0.328πr2. To correct for the reduction of contact area by interstitial gaps between cells, we determined gap sizes (supplementary material Fig. S4B-D and mathematical analyses, Section 4). The ratio of the diameter of the corrected contact area to the diameter corresponding to amax is d=dc/dmax and hence a/amax=d2.
In tissues, cell boundary fluctuations partly overcome T1 energy barriers and are thus essential determinants of tissue viscosity (Marmottant et al. 2009). ξ is the amplitude of energy fluctuations at cell contacts. In analogy to Marmottant et al. (2009) we propose that ξ is proportional to contact area fluctuations by ξ=(qa)(2β*), with q being an average relative area change derived from the standard deviation s of contact lengths dc (i.e. diameter of contact area). For small s<dc/4, Δa=π(dc/2)2−π((dc/2)−s)2≈πdcs=4as/dc, and therefore q=4s/dc. For larger s, as in chordamesoderm (supplementary material Table S1), q is overestimated by about 20%, which barely affects viscosity estimates (see supplementary material Fig. S5). f is the attempt frequency of fluctuation-driven rearrangement (Fig. 2F). Cell volume V (in the equation of Marmottant et al.) is given by the radius r of dissociated cells; the relative deformation during a rearrangement step for a hexagonal array is δε≈0.423.
Inserting these parameters into the above equation of Marmottant et al (2009) yields Eqn (1), which connects tissue viscosity η to surface tension σ. Viscosity and surface tension are also linked through the ratio ν=σ/η, with ν having the dimension of a velocity. This characteristic velocity takes on the constant value νc for some tissues.
An important difference between tissues and foams is that energy barriers in tissues are partly overcome by metabolically driven cell boundary fluctuations. In explants, cell contact lengths fluctuated substantially (Fig. 2F; supplementary material Movie 4), with occasional spontaneous cell rearrangements (Fig. 2G). This effect leads to Newtonian tissue flow at sufficiently low stress and, together with the energy barrier, it determines viscosity (Marmottant et al., 2009) (Box 2). We calculated tensions at cell contacts β* as described above, and contact areas from cell sizes and the sizes of interstitial gaps (Box 2, Table 1; supplementary material Fig. S4, Table S1 and mathematical analyses, Section 4). We also measured the amplitudes and frequencies of contact area fluctuations from time-lapse recordings (Fig. 2F, Table 1; supplementary material Table S1). Substituting these parameters into the equation of Marmottant et al. (2009) yielded (Box 2):
1The first term in Eqn (1) amalgamates all dimensionless factors and contains the cell radius r and the correction for interstitial gaps d, linking viscosity to tissue geometry. The second term uses the relative amplitude q (Box 2) and the frequency f to describe contact area fluctuations. The third term expresses the mechanical tension β* as the difference between cell cortical and tissue surface tension.
We calculated the average viscosity η from Eqn (1) for several gastrula tissues, and found good agreement with measured values (Table 1). We also used a calculation of error propagation to estimate how the within-tissue variation of measured parameters translates into that of viscosity, and found that calculated and measured standard deviations were similar (Table 1). The overall agreement between calculation and measurement suggested that viscosity is indeed largely determined by the cell rearrangement process. Cytoplasmic flow during cell deformations at and between T1 events could contribute to viscosity. However, Xenopus egg cytoplasmic viscosity is 10-30 mPa·s (Valentine et al., 2005), which is orders of magnitudes lower than measured tissue viscosities.
Eqn (1) can be used to analyse the proportionality between tissue surface tension and viscosity. For gastrula tissues, the fluctuation term is nearly constant, as the fluctuation amplitude has only a small effect on viscosity (supplementary material Fig. S5) and the fluctuation frequency is similar in all tissues (Table 1; supplementary material Table S1). Cell size differs significantly between endoderm and other tissues, but in general we expect that cell radius will vary much less than viscosity, i.e. less than 1000-fold. We therefore propose that it is the third term, the tension term of Eqn (1), that is mainly responsible for the proportionality of cohesion and viscosity in νc tissues.
With the tension term β*=β−σ, viscosity is obviously not proportional to surface tension directly, but to the difference between cortical and surface tension. Therefore, a requirement for proportionality is that cortical and surface tension themselves vary proportionally, i.e. that relative adhesiveness α=σ/β remains constant. This is indeed the case for normal gastrula tissues (Table 1). When, instead, cortical tension β is kept constant, the tension term increases with decreasing surface tension σ. This counterintuitive relationship explains why C-cadherin knockdown reduced tissue cohesion but increased viscosity (Table 1): cortical tension at cell boundaries is less strongly reduced during adhesion (i.e. β* is higher), hence cells can be separated more easily (i.e. cells are less adhesive), whereas energy barriers for cell rearrangement are higher, in proportion to the increased β* (see Fig. 2H).
Interpretation of the constant νc=1.8 µm/min: flow in response to external stress
Eqn (1) explains how proportionality between tissue surface tension and viscosity can be achieved, but not why proportionality is actually implemented in νc tissues. It is also unclear whether the constant νc takes on the value of 1.8 µm/min ‘accidentally', or whether this reflects some basic property of cell movement in all νc tissues. To derive the cell biological meaning of νc, we will relate tissue flow to cell-level events. We will argue that 1.8 µm/min is the rearrangement velocity of a cell in all tissues, and that νc tissues are tuned for maximal efficiency of rearrangement, which in turn implies proportionality between surface tension and viscosity. We will first consider tissue flow in response to external stress; for example, the compression of a cell aggregate between plates or explant rounding driven by tissue surface tension.
According to Marmottant et al. (2009), the occurrence of T1 events, i.e. the fraction of cells rearranging at any time point, increases in proportion to the applied stress. We therefore model tissue flow assuming that, for normal, moderate stress levels, only a fraction of cells is moving directionally at any time, with individual cells translocating at a fixed velocity relative to their neighbours for the duration of a T1 event, and pausing between events. This rearrangement velocity Vr depends on the contraction and expansion rates of cell boundaries during T1 rearrangements, and on a geometric factor that accounts for the orientation of cell boundaries relative to the applied stress (Box 3).
Box 3. Cell rearrangement velocity
In the schematic drawing, it is seen that a cell moves about one cell diameter (2r) during a T1 step of duration T. During the process, the boundary shown in red contracts and the blue boundary expands. The rearrangement velocity Vr of a cell depends on the contraction rate c1 of a respective cell boundary during the first phase of a T1 transition and the expansion rate c2 of the new boundary during the second phase. Both rates are very similar, and we used the harmonic mean c of c1 and c2 to describe boundary length changes (Table 2).
To derive Vr from c, we considered that the movement of a cell relative to its neighbour during a T1 event occurs maximally at velocity Vrmax=c cos30°, assuming an angle of 120° at cell edges, as shown in the drawing. With random packing, planes of cell boundaries are variably oriented relative to the direction of stress, and contribute to different degrees to velocity. To correct for this effect, we considered a tilted cell group, as shown on the righthand side of the drawing, to estimate the minimal velocity Vrmin during T1. At sixfold rotational symmetry, the orientation of the cell pair shown relative to the direction of stress maximally deviates by 30° in a plane perpendicular to the cell edge, in which case the component of Vrmax in the direction of the overall resultant movement is Vrmax cos30°. Applying this also to the remaining, perpendicular plane, the minimal velocity contributed by unfavourably oriented contacts will be Vrmin=Vrmax cos230° and with Vrmax=c cos30°, Vrmin=c cos330°. We approximated Vr by the average of Vrmax and Vrmin to obtain Vr≈1.8 µm/min (Tables 1 and 2). This coarse estimate can be refined; for example, by considering that cells often elongate to varying degrees in the direction of movement. This effect should tend to increase Vr. With k the fraction of cells moving via T1 transitions at a given instant, or equivalently with k the likelihood of a cell moving at that instant, the long-term average relative velocity is Vr=kVr.
Contraction and expansion rates were similar in all gastrula tissues and not significantly affected by treatments such as cadherin inhibition (Table 2), although a slight tendency to increase with cell size was noted (correlation coefficient 0.186). The average rate of contraction or expansion was 2.3 µm/min. During Drosophila germ band extension, myosin II-controlled contraction rates of 2.5 and 1.6 µm/min have been measured (Rauzi et al., 2010), and, for a similarly rearranging Xenopus chordamesoderm cell, rates varied around 2.6 µm/min (Shindo and Wallingford, 2014). Rates are of the same order in cell contractions during micropipette aspiration-induced cell rearrangements in aggregates (Guevorkian et al., 2011). We propose that the similarity of these rates reflects a constraint of the molecular mechanisms of cell boundary dynamics, which sets a time-scale for rearrangement. From the rate for Xenopus tissues, we estimated a rearrangement velocity Vr of 1.8 µm/min (Box 3, Table 2), and expect it to be similar in other tissues.
Measured parameter values for boundary contraction and extension
The two characteristic velocities that we have identified, Vr and νc, are of the same magnitude, and we asked how they are related. We had defined ν=σ/η. To express Vr similarly, we used the equation for the shear movement of a Newtonian fluid, which relates the applied stress to viscosity and to the velocity gradient generated in the fluid (Box 4). The average speed of a cell layer relative to an adjacent layer, which depends on Vr, was used to derive the velocity gradient. Moreover, we expressed stress as a fraction or multiple p of a standard stress conveniently defined in terms of tissue surface tension (Box 4). This gave Vr=(σ/η)(1/kσ), and thus ν=kσVr, with kσ denoting the fraction of cells moving at standard stress. In other words, the ratio ν corresponds to the rearrangement velocity of 1.8 µm/min multiplied by the fraction of cells moving under standard conditions. When ν=νc=1.8 µm/min, all cells should be moving, i.e. kσ=1; at lower ν values, the fraction of moving cells at a comparable stress is correspondingly lower, i.e. kσ<1. Apparently, kσ measures the efficiency of rearrangement, and νc tissues are tuned for maximal efficiency. The constancy of νc implies also a constant ratio σ/η for these tissues, i.e. proportionality.
Box 4. Tissue flow
According to Marmottant et al. (2009), tissues behave like Newtonian fluids below a critical stress, i.e. their shear movement is described by τ=η(dV/dx), where τ is an applied stress and dV/dx a velocity gradient (see drawing, left part). To generate in different tissues the same shear rate, or velocity gradient, the applied stress must in each case be proportional to the viscosity of the tissue, τ∼η. For νc tissues in which viscosity is proportional to surface tension (i.e. η∼σ), this stress must be proportional to tissue surface tension, that is, τ∼σ.
This consideration motivated us to express stress in terms of surface tension. Surface tension has the dimension of [energy/area] or equivalently [force/length]; to obtain the dimension of stress, i.e. [force/area], it has to be divided by a length. We defined σ/L as a standard stress τσ with L a characteristic cell length, such as the average centre-to-centre distance between cells (see drawing). The actual applied stress can then be expressed as a fraction or multiple p (p being a dimensionless factor) of the standard stress, τ=pτσ=pσ/L. The long-term, averaged velocity of cells relative to neighbours, Vr, divided by length L, represents the shear rate of the flow, i.e. the average speed of a cell layer moving relative to an adjacent layer (see drawing). With Vr=kVr (Box 3), the shear rate becomes Vr/L=kVr/L. Taken together, we relate stress to shear rate as pσ/L=ηkVr/L. At standard stress, i.e. at p=1, and with kσ denoting k under the condition of standard stress, the equation becomes σ=ηkσVr and thus σ/η=kσVr=ν. We conclude that the ratio ν=σ/η corresponds to the fraction of cells moving multiplied by their rearrangement velocity, and for ν=νc it follows that kσ=νc/Vr. If νc=Vr, then kσ=1. We can interpret kσ as indicating the efficiency of cell rearrangement.
If movement is driven by active, controlled cell-cell boundary constrictions instead of an external force, the dimensionless factor p becomes the fraction of cells actively engaged in rearrangement at a given time point. To separate two cells, a controlled increase in tension in each cell at the contact zone equivalent to surface tension σ, i.e. from β* to β, is sufficient (Fig. 2H). The stress generated in this way at the cell level amounts to σ/L, which corresponds to the standard stress τσ as defined above. The overall stress at the tissue level is approximately proportional to this standard stress, and to the fraction p of cells activated to constrict, that is, τ≈pσ/L, assuming that the contracting boundaries are sufficiently well oriented to the direction of movement. Altogether, this leads to the same equation that describes rearrangement due to external stress, only with p denoting now the fraction of activated cells instead of the fraction of cells prompted to move passively.
For normal gastrula tissues ν is indeed around 1.8 µm/min, and with similar values calculated for Vr we obtained an average kσ of 1.0 (range 0.8-1.3; Table 1), suggesting maximal efficiency of movement. However, C-cadherin knockdown led to a ν of only 0.3 µm/min and hence to a kσ of 0.2 (Table 1). Similarly, the adhesion-diminishing C-cadherin construct EPΔC reduced ν to 0.8 µm/min and kσ to 0.6 (Table 1). In such tissues, a stress considerably higher than standard stress would be required to generate the same shear rate. Notably, the stress that drives tissue flow is normally much lower than the standard stress. In the gastrula, the actual average cell velocity in tissues is 0.2-0.5 µm/min instead of 1.8 µm/min (supplementary mathematical analyses, Section 5), i.e. the fraction of cells rearranging at any time is only 0.1-0.3. Since for these tissues kσ≈1, the effective driving stress should be one-tenth to one-third of the standard stress.
To see how the efficiency factor kσ is determined, we derived Eqn (S1) (supplementary mathematical analyses, Section 6), which relates kσ to relative adhesiveness α. It shows that for gastrula tissues to have a kσ≈1 requires that 0.6<α<0.8 (supplementary mathematical analyses, Section 6), in agreement with the measured values for relative adhesiveness. Apparently, only a narrow range of α values is compatible with maximal efficiency of rearrangement: at low α, kσ is expected to be smaller than 1. In fact, C-cadherin knockdown reduced α to 0.38, for which Eqn (S1) predicts a kσ of 0.3, close to the value of 0.2 that was determined independently from Vr and νc (Table 1).
Interpretation of the constant νc: flow driven by internally generated stress
An important difference between foams and tissues is that tissues can actively generate internal stress by mechanochemical processes to drive controlled cell rearrangements. To see whether the concepts derived above for external stress apply here as well, we now consider active cell intercalation, as seen for example during Xenopus convergent extension (Keller et al., 2008; Shindo and Wallingford, 2014) or Drosophila germ band extension (Bertet et al., 2004). We replace the externally applied stress by an active, regulated constriction of cell boundaries that is superimposed on the random fluctuations. The boundaries are preselected according to their proper spatial orientation by some patterning mechanism. For example, in convergent extension and in germ band extension, anterior and posterior contacts of cells are determined to contract by an activin signalling gradient (Ninomiya et al., 2004) or striped pair-rule gene expression (Zallen and Wieschaus, 2004) to ensure anteroposterior tissue elongation.
To separate two cells in such a controlled T1 event, an increase in cortical tension in the contact area equivalent to surface tension is sufficient: this would correspond to the cell attachment process running in reverse (Box 4, Fig. 2H). Cell boundary constriction can be generated by myosin activation (Lecuit and Lenne, 2007; Skoglund et al., 2008; Shindo and Wallingford, 2014). Alternatively, a transient suspension of the process that normally downregulates cortical density at contacts could ensure cell separation regardless of adhesion strength. Regardless, the stress generated at the cell level will correspond to the standard stress, identifying it as the average stress required per cell to drive a T1 event (Box 4).
In a simple model, the overall stress at the tissue level under these conditions is proportional to the fraction p of cells actively engaged in constriction at a given instant, and the same equation as for the flow due to external stress applies (Box 4). Consequently, when kσ≈1, essentially all cells attempting a controlled boundary contraction at a given time would indeed be moving through a T1 process, whereas in tissues with low kσ only a small fraction of attempts would be productive. Note that progression through T1 can be stepwise when contractions are pulsed (Skoglund et al., 2008; Shindo and Wallingford, 2014), such that short episodes of cell movement at velocity Vr alternate with pauses, and the overall movement becomes smoother. Irrespective, here too kσ is a measure of the efficiency of cell rearrangement. Overall, the available data suggest that νc tissues are tuned for maximal efficiency: with a Vr of 1.8 µm/min for these tissues, kσ≈1 regardless of whether rearrangement is intrinsically controlled or a response to external stress. This propensity for maximal efficiency can explain the constancy of the ratio νc in many tissues.
Whether the same tissue is deformed by an external stress or by an internal stress of the same magnitude, the same viscosity should be measured in both cases. We confirmed this prediction for chordamesoderm, which lengthens by active cell rearrangement. From its rate of elongation (Ninomiya and Winklbauer, 2008) and from the driving force measured by Moore et al. (1995), we calculated a viscosity of 5.7 kPa·s (supplementary mathematical analyses, Section 7), which is close to the 6.9 kPa·s determined from surface tension-driven explant rounding. Thus, internally generated stress produces the same shear rate as if stress of the same magnitude were applied externally. As a principle, this allows the calculation from simple, external stress-driven movements, such as explant rounding, of mechanical parameters that are relevant for active cell rearrangement, for example through the determination of the characteristic velocity ν.
DISCUSSION
We found that the ratio of tissue surface tension to viscosity is remarkably constant over a wide range of values, and identifies the νc subclass of liquid-like tissues. We further showed that tissue viscosity can be predicted quantitatively from a model of tissue flow based on foam rheology (Marmottant et al., 2009). Building on this model, we argued that cell rearrangement is characterised by a constant derived from the contraction and expansion rates of cell boundaries during T1 events: these rates translate into the velocity by which individual cells move past each other during a rearrangement step. It is this velocity that corresponds to the ratio of surface tension to viscosity in νc tissues. As a consequence, cell rearrangement is predicted to be maximally efficient, i.e. the maximal number of cells is moving simultaneously at a given applied stress.
This analysis implies that movement is due to cell rearrangements via T1 processes, as is indeed common in epithelia. This also seems likely for rearrangements driven by external forces. However, models of active cell intercalation in three-dimensional tissues usually employ the concept of intercellular migration (Gumbiner, 2005), in which cells crawl across each other's surfaces with the use of locomotory protrusions. For example, this has been the predominant model for convergent extension of the Xenopus chordamesoderm (Keller et al., 2008). However, an epithelial-type T1 mechanism has been proposed recently for this process (Shindo and Wallingford, 2014), and this kind of rearrangement might be more common than previously thought. Moreover, the limiting steps of intercellular migration are probably also cell contact-changing processes, such as membrane detachment and reattachment, the efficiency of which could depend on the ratio ν similarly to cell boundary contraction and expansion.
The present analysis draws attention to the striking variability in mechanical parameters both within and between tissues. Strong within-tissue variation often represents true variability and not measurement errors (von Dassow and Davidson, 2007, 2009). For some parameters, the origin of the variation is obvious. For example, contact angles between cells vary heavily due to cell boundary fluctuations, whereas the measurement error is small. Fluctuation frequency also shows strong within-tissue variability as it does in other systems, e.g. in pulsating actomyosin networks (Martin et al., 2009; David et al., 2010). The measurement of tissue surface tension by different methods (parallel plate compression versus drop shape analysis) and in tissues from different species revealed a strikingly similar variability in each case (Davis et al., 1997; Schötz et al., 2008). The stiffness of a tissue can also vary substantially, without affecting morphogenesis (von Dassow and Davidson, 2007, 2009).
While most mechanical parameters show considerable within-tissue variation, some also differ enormously between different types of tissue whereas others do not. For νc tissues, we predict that relative adhesiveness, magnitude of cell boundary fluctuations and velocity of cells in T1 events will be nearly constant. This would keep the ratio ν constant as tissue surface tension and viscosity vary widely. Moreover, constant relative adhesiveness implies a similar effect of cadherins on cortical tension reduction in different νc tissues. Furthermore, we predict that the vast range of surface tension values is matched by a similar range of cell cortical tensions, as cortical tension must always be larger than surface tension (Brodland et al., 2009; Manning et al., 2010). A range of cortical tension values consistent with this assumption is indeed observed for νc tissues (supplementary material Table S3). As cortical tension is related to the strength of cell-cell adhesion, its variation could be explained by requirements for high or low tissue cohesion. Notably, however, cortical tension also varies strongly among solitary cells (supplementary material Table S3).
At its highest level, viscosity increases without a corresponding increase in tissue surface tension. This could be due to high cortical tensions, which are no longer proportionally reduced at contacts but by a fixed amount that corresponds to a limiting, high value of tissue surface tension. Another possibility is that the amplitude and/or frequency of cell boundary fluctuations are decreased to raise tissue viscosity. For example, a slight reduction in the relative fluctuation amplitude q below 0.5 (supplementary material Fig. S5) would cause a sharp increase in viscosity, generating a tissue with more solid-like behaviour. Lastly, when a stiff extracellular matrix is present, the time required for its enzymatic remodelling could determine flow rates. This condition would mark the transition to tissues that behave essentially as solids and show predominantly elastic responses to mechanical stress even at long time-scales (Levental et al., 2007).
By applying concepts of fluid flow to cell rearrangement, we gained basic insights into the mechanics of cell movement in tissues. In particular, by building on the work of Manning et al. (2010) and Marmottant et al. (2009) we identified cell-level parameters that determine the mechanical properties of tissues relevant for cell rearrangement. We did not characterise the relationship between fluctuation frequency and amplitude, cell boundary contraction rate and cell size in full detail, but these parameters are clearly not independent. Furthermore, cortical tensions at exposed and contact surfaces and relative adhesiveness are connected. We propose that when choosing boundary contraction/expansion rate c, frequency f and cell radius r as easily measurable variables to represent the kinetic aspect of cell rearrangement, and cortical tension β and relative adhesiveness α to represent adhesion strength, these five variables are sufficient to express tissue-level properties, such as viscosity or surface tension, and important tissue characteristics, such as ν and kσ. In tissues with significant interstitial gaps a respective correction factor would have to be included.
Our analysis draws attention to numerous cell biological questions. For example, it would be interesting to see how the orders-of-magnitude differences in cortical tension between cell types are implemented at the cytoskeletal level, or, on the other hand, which cellular processes limit rates for cell boundary contraction and expansion to a narrow range. The relationship between cadherin density and the reduction of cortical tension (Maître et al., 2012), and the molecular basis of the random cell contact fluctuations, which are so important for the facilitation of cell rearrangement, are currently also not understood. Finally, the concepts developed here can guide our understanding of the large-scale design of tissues of vastly varying mechanical properties, and provide a framework for their analysis and comparison.
MATERIALS AND METHODS
Embryos and explants
Xenopus laevis embryos were obtained, explants were made and cultured in Modified Barth's Solution (MBS) or dissociated into single cells in Ca2+/Mg2+-free MBS, as described (Winklbauer, 1990).
Measurements of cell and tissue parameters
Axisymmetric drop shape analysis (ADSA) (Del Rio and Neumann, 1997) was used to determine tissue surface tensions for regions of the early gastrula, as described (Luu et al., 2011). For unknown reasons, surface tension of the ectoderm, but not of the other tissues, alternated between high and low states over the course of many months (although apparently not in a seasonal rhythm). Tissue viscosity was determined by measuring the rate of rounding up of ellipsoidal gastrula explants (Gordon et al., 1972) (for details see supplementary mathematical analyses, Section 1). To view cell rearrangement, rounding explants were slightly compressed under a coverslip to keep sufficient surface area in focus at the high magnification required. As a measure of cell size, the diameter of dissociated single cells was determined. To quantitate cell contact fluctuation, contact lengths were measured at 1 min intervals in time-lapse recordings of explants filmed under indirect illumination. Contact angles between cells were measured in cell pairs or at the periphery of tissue explants (Stirbat et al., 2013) fixed and fractured after rounding up for 0.5-1 h. All measurements used the Carl Zeiss Axiovision program; accuracy was ±0.2° (s.d.; n=54) for angles, ±0.28 µm (s.d.; n=21) for lengths. A list of the symbols used for the main parameters can be found in supplementary material Table S4.
Microscopy
For immunostaining, embryos were fixed with 2% paraformaldehyde and post-fixed with Dent's fixative; cryostat sections were immunostained with β-catenin antibody (Santa Cruz Biotechnology, H-102; dilution 1:200). For F-actin staining, dissociated cells were seeded on BSA-coated dishes, fixed with 4% paraformaldehyde/0.1% Triton X-100 and incubated with Alexa Fluor 488 phalloidin (Invitrogen). Pixel intensity was determined using the Carl Zeiss Axiovision program (pixel size 625 nm2). Actin dynamics were viewed in LifeAct mRNA-injected cells using a Carl Zeiss LSM 510 confocal microscope. For transmission electron microscopy, gastrula explants were fixed in 2.5% glutaraldehyde/2% formaldehyde, post-fixed with 1% osmium tetraoxide and infiltrated with Spurr's resin as described (Kurth et al., 2010). Ultrathin sections were stained with uranyl acetate and post-stained with Reynold's lead citrate.
Microinjections
Embryos were microinjected in 4% Ficoll solution with mRNA or MOs (GeneTools) at the four-cell stage, and kept in 1/10 MBS at 15°C until the required stage. Embryos were injected with previously characterised C-cad-MO (Ninomiya et al., 2012), standard control MO, C-cadherin mRNA or EPΔC mRNA (Ninomiya et al., 2012), and with M-PAPC mRNA (Kim et al., 1998), Pak1 mRNA or LifeAct mRNA. The experiments involving X. laevis embryos conform to the regulatory standards of the University of Toronto, Canada.
Acknowledgements
We thank B. Gumbiner and E. M. DeRobertis for plasmids; and A. McGuigan, A. Bruce, T. Harris, U. Tepass, D. Barua and S. Parent for suggestions to improve the manuscript.
Footnotes
Competing interests
The authors declare no competing financial interests.
Author contributions
R.W. conceived the study, designed experiments and wrote the manuscript with contributions from R.D. and the other authors. R.D. developed the tools for data analysis. R.W. and R.D. analysed and interpreted the data. R.D., O.L., E.W.D., J.W.H.W., M.N. and R.W. performed the experiments.
Funding
Funding was provided by the Canadian Institutes of Health Research [MOP-53075].
Supplementary material
Supplementary material available online at http://dev.biologists.org/lookup/suppl/doi:10.1242/dev.104315/-/DC1
- Received September 27, 2013.
- Accepted July 17, 2014.
- © 2014. Published by The Company of Biologists Ltd
References
Article navigation
Related articles
Cited by...
More in this TOC section
Similar articles
Other journals from The Company of Biologists
An interview with Swathi Arur
Swathi Arur joined the team at Development as an Academic Editor in 2020. Her lab uses multidisciplinary approaches to understand female germline development and fertility. We met with her over Zoom to hear more about her life, her career and her love for C. elegans.
Jim Wells and Hanna Mikkola join our team of Editors
We are pleased to welcome James (Jim) Wells and Hanna Mikkola to our team of Editors. Jim joins us a new Academic Editor, taking over from Gordan Keller, and Hanna joins our team of Associate Editors. Find out more about their research interests and areas of expertise.
New funding scheme supports sustainable events
As part of our Sustainable Conferencing Initiative, we are pleased to announce funding for organisers that seek to reduce the environmental footprint of their event. The next deadline to apply for a Scientific Meeting grant is 26 March 2021.
Read & Publish participation continues to grow
“I’d heard of Read & Publish deals and knew that many universities, including mine, had signed up to them but I had not previously understood the benefits that these deals bring to authors who work at those universities.”
Professor Sally Lowell (University of Edinburgh) shares her experience of publishing Open Access as part of our growing Read & Publish initiative. We now have over 150 institutions in 15 countries and four library consortia taking part – find out more and view our full list of participating institutions.
Upcoming special issues
Imaging Development, Stem Cells and Regeneration
Submission deadline: 30 March 2021
Publication: mid-2021
The Immune System in Development and Regeneration
Guest editors: Florent Ginhoux and Paul Martin
Submission deadline: 1 September 2021
Publication: Spring 2022
Both special issues welcome Review articles as well as Research articles, and will be widely promoted online and at key global conferences.
Development presents...
Our successful webinar series continues into 2021, with early-career researchers presenting their papers and a chance to virtually network with the developmental biology community afterwards. Here, Michèle Romanos talks about her new preprint, which mixes experimentation in quail embryos and computational modelling to understand how heterogeneity in a tissue influences cell rate.
Save your spot at our next session:
10 March
Time: 9:00 (GMT)
Chaired by: Thomas Lecuit
Join our mailing list to receive news and updates on the series.