Correlating cell shape and cellular stress in motile confluent tissues

Xingbo Yang, Dapeng Bi, Michael Czajkowski, Matthias Merkel, M. Lisa Manning, and M. Cristina Marchetti
PNAS November 28, 2017 114 (48) 12663-12668; published ahead of print November 14, 2017 https://doi.org/10.1073/pnas.1705921114

  1. Edited by William Bialek, Princeton University, Princeton, NJ, and approved October 11, 2017 (received for review April 14, 2017)

Significance

Using a self-propelled Voronoi model of epithelia known to predict a liquid–solid transition, we examine the interplay between cell motility and cell shape, tuned by cortex contractility and cell–cell adhesion, in controlling the mechanical properties of tissue. Our work provides a unifying framework for existing, seemingly distinct notions of stress in tissues and relates stresses to material properties. In particular, we show that the temporal correlation function of shear stresses can be used to define an effective tissue viscosity that diverges at the liquid–solid transition. This finding suggests a unique way of analyzing traction force microscopy data that may provide information on tissue rheology.

Abstract

Collective cell migration is a highly regulated process involved in wound healing, cancer metastasis, and morphogenesis. Mechanical interactions among cells provide an important regulatory mechanism to coordinate such collective motion. Using a self-propelled Voronoi (SPV) model that links cell mechanics to cell shape and cell motility, we formulate a generalized mechanical inference method to obtain the spatiotemporal distribution of cellular stresses from measured traction forces in motile tissues and show that such traction-based stresses match those calculated from instantaneous cell shapes. We additionally use stress information to characterize the rheological properties of the tissue. We identify a motility-induced swim stress that adds to the interaction stress to determine the global contractility or extensibility of epithelia. We further show that the temporal correlation of the interaction shear stress determines an effective viscosity of the tissue that diverges at the liquid–solid transition, suggesting the possibility of extracting rheological information directly from traction data.

It is now broadly recognized that the transmission of mechanical forces can be as important as genetics and biochemistry in regulating tissue organization in many developmental processes, including embryogenesis, morphogenesis, wound healing, and cancer metastasis (1⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓–12). To make quantitative predictions for large-scale cell remodeling in tissues, we must understand their material properties, such as stiffness and viscosity, as well as the forces that build up inside them, characterized by local pressures and stresses. Motivated by experiments highlighting the slow glassy dynamics of dense epithelia, work by us and others has suggested that monolayers of motile cells may form glassy or jammed states and that a relatively small change of parameters may trigger a change from an elastic response to a state with fluid-like behavior (13⇓–15). This liquid–solid transition is tuned by the interplay of cell–cell adhesion and cortex contractility, manifested in cellular shape, and by cell motility. This suggestion has been verified experimentally in specific cell types (16), indicating that the paradigm of tissues as active materials may be a useful way of organizing experimental data and classifying large-scale tissue behavior in terms of a few effective parameters.

In contrast, quantifying stresses and pressures in active, motile tissues is largely an open problem. One possible reason is that there are different definitions of stresses and pressures that arise naturally in different experiments and simulations, and it is not immediately clear how they are related to one another or under which conditions each definition applies. For example, traction force microscopy (TFM) is a powerful tool that probes the dynamic forces cells exert onto soft substrates by measuring the substrate deformations (10, 11, 17⇓–19). In some experiments, intercellular stresses are extracted from the traction forces, using finite-element analysis under the condition that traction forces are balanced by cellular interactions. The resulting stress maps reveal a highly dynamical and heterogeneous mechanical landscape, characterized by large spatial and temporal fluctuations in both normal and shear stresses (6, 10, 11, 18, 20). A key assumption used in the TFM approach to infer stresses from tractions is that the cell layer can be described as a continuum linear elastic material (18).

A second important set of mechanical inference methods also predicts cellular stresses, but in static tissues. These methods rely on advances in imaging techniques that provide a spatially resolved view of tissue development during morphogenesis, with visualization of the cell boundaries of 2D cell sheets (21⇓⇓⇓–25). Assuming mechanical equilibrium, one can then infer the tensions along cell edges and pressures within each cell from the cell configurations. This method provides a spatial distribution of intercellular stresses and directly couples mechanics to morphology, but it does not capture temporal stress fluctuations arising from the nonequilibrium nature of the tissue. It has been used successfully to characterize cell morphology in the development of the Drosophila wing disk (23, 26), where cellular rearrangements are slow on the timescales of interest, demonstrating that the analysis of cell shapes can provide fundamental insight on the mechanical state of tissues in developmental processes.

A third line of research has focused on the homeostatic pressure that tissues exert on their containers or surroundings. The homeostatic pressure has been proposed as a quantitative measure of the metastatic potential of a tumor (27, 28). It is defined as the force per unit area that a confined tissue would exert on a moving piston permeable to fluid, and hence it represents an active osmotic cellular pressure. The existing literature has focused on the contribution to the homeostatic pressure from tissue growth due to cell division and death, but in general other active processes, such as cell motility and contractility, will also contribute to the forces exerted by living tissues on confining walls. A related body of work has investigated the pressures generated by motile particles, such as active colloids. In these highly nonequilibrium systems that break time-reversal symmetry, there is a contribution to the total pressure called the swim pressure generated entirely by persistent motility (29⇓–31).

So far, there is no unifying theory for these seemingly distinct notions of pressure and stress or their relationship to material properties. In this paper, we show that a recently proposed self-propelled Voronoi (SPV) model of epithelia (Fig. 1) (15) provides a natural framework for unifying these ideas. One of the benefits of the SPV model is that it explicitly accounts for the forces that motile cells exert on the substrate. This allows us to develop a generalized mechanical inference method to infer cellular stresses from traction forces and to show that these match the stresses calculated from instantaneous cell shape, relating TFM data and mechanical inference techniques in motile tissues. Additionally, our method provides absolute values for junctional tensions and pressure differences. This is in contrast to equilibrium mechanical inference, which yields only relative forces (21, 23).

Fig. 1.
Fig. 1.

(A) Illustration of the SPV model, where cells are represented by polygons obtained via a Voronoi tessellation of initially random cell positions, with a self-propulsion force applied at each cell position. (B) Phase diagram in the (P0,v0) plane based on the value of the shape parameter q (color scale), with the phase boundary (red crosses) determined by q=3.813.

There are two additive contributions to the mechanical stress that describe the forces transmitted in a material across a bulk plane. The first one represents the flux of propulsive forces through a bulk plane carried by particles that move across it. The second one describes the flux of interaction forces across a bulk plane. We demonstrate that the generalized mechanical inference measurements probe the latter, which we denote interaction stresses. The former, which we denote the tissue swim stress, approximates the contribution from cell motility to the osmotic pressure generated by cells immersed in a momentum-conserving solvent on a semipermeable piston and hence to the tissue homeostatic pressure. The tensorial sum of the swim stress and the interaction stress is the total stress. The normal component of the total stress determines whether a tissue will tend to exert extensile or contractile forces on its environment, which is an important consideration in wound healing and cancer tumorigenesis.

An obvious open question, then, is how these stresses vary as a function of material properties. We find that the normal component of the interaction stress is contractile in both the solid and the liquid due to the contractility of the actomyosin cortex, although much more weakly so in the liquid state. In contrast, the normal component of the motility-induced swim stress is always extensile, corresponding to a positive swim pressure, although its magnitude depends on the phase: In a solid the swim pressure is negligible, while in the fluid it can be significant. This can result in a change in sign of the total mean stress: Indeed, we find it is always contractile in the solid state but becomes extensile deep in the liquid state when cell motility exceeds actomyosin contractility.

Because the transition from contractile to extensile does not coincide with the fluid to solid transition, it is natural to ask whether the stress displays any signatures of the fluid–solid transition. We develop a definition for the effective viscosity of the tissue that can be extracted from the temporal correlation of the interaction shear stress and find that it diverges as the tissue transits from the liquid state to the solid state. Importantly, this theoretical prediction suggests that TFM combined with mechanical inference can provide rheological information about the tissue and could be tested by a new analysis of experimental data.

Results and Discussion

SPV Model.

The SPV model describes an epithelium as a network of polygons. Each cell i is endowed with a position vector 𝒓i, and cell shape is defined by the Voronoi tessellation of all cell positions (Fig. 1), which has been shown to provide a good representation of some real epithelia, such as the blastoderm of the red flour beetle Tribolium castaneum and the fruit fly Drosophila melanogaster (32). Like for vertex models, tissue forces are obtained from an effective energy functional E({𝒓i}) for N cells, given by (7, 15, 26, 33, 34)E=∑i=1NEi,Ei=KA(Ai−A0)2+KP(Pi−P0)2,[1]with Ai and Pi the cross-sectional area and perimeter of the ith cell, respectively. The first term in Eq. 1 arises from incompressibility of the layer in three dimensions and its resistance to height fluctuations, with A0 a preferred cross-sectional area. The second term represents the competition between cortical tensions from the actomyosin network at the apical surface and cell–cell adhesions from adhesive complexes at intercellular junctions (26), with P0 a preferred perimeter resulting from this competition. We simulate N cells in a square box of area AT, with A¯=AT/N the average cell area and with periodic boundary conditions. The system is initialized with a set of N random cell positions, independently drawn from a uniform distribution. Throughout the simulations, we set A¯=A0=1 unless otherwise noted, and KA=KP=1.

Each Voronoi cell is additionally endowed with a constant self-propulsion speed v0 along the direction of polarization 𝒏^i=(cos⁡θi,sin⁡θi) describing cell motility. The dynamics of each Voronoi cell are governed by∂t𝒓i=μ𝑭i+v0𝒏^i,∂tθi=2Drηi(t),[2]where 𝑭i=−∇iE is the force on cell i and μ is the mobility. The direction of cell polarization is randomized by orientational noise of rate Dr, with ⟨ηi(t)⟩=0 and <ηi(t)ηj(t′)>=δijδ(t−t′). The timescale τr=1/Dr controls the persistence of single-cell dynamics. As in self-propelled particle (SPP) models, an isolated cell performs a persistent random walk with a long-time translational diffusivity D0=v02/(2Dr) (29, 35, 36). After each time step, a new Voronoi tessellation is generated based on the updated cell positions. The cell shapes are determined in the process and the exchange of cell neighbors occurs naturally through topological transitions (13).

We showed in ref. 15 that the SPV model exhibits a transition from a solid-like state to a fluid-like state upon increasing the single-cell motility v0, the persistence time τr, or the cell shape parameter P0/A0 that characterizes the competition between cell–cell adhesion and cortical tension. The phase diagram in the (P0,v0) plane is reproduced in Fig. 1B. The transition is identified by setting the shape index q=⟨Pi/Ai⟩ to the value q= 3.813, where ⟨…⟩ denotes the average over all cells. It was shown in ref. 15 that the transition line located by q= 3.813 coincides with the one based on the vanishing of the effective diffusivity obtained from the cellular mean-square displacement. Note that for fixed system size AT and cell number N, the preferred cell area A0 does not affect the interaction forces or cellular shapes. Hence the solid–fluid transition is insensitive to A0, as shown analytically in SI Text. The preferred area A0 only shifts the total pressure of the tissue by a constant.

Developing and Validating Traction-Based Mechanical Inference.

It is well established that in a model tissue described by the energy Eq. 1 , the mechanical state of cell i is characterized by a local stress tensor σαβ(i)int given by (21, 23, 37)σαβ(i)int=−Πiδαβ+12Ai∑ab∈iTabαlabβ,[3a]Πi=−∂E∂Ai, 𝑻ab=∂E∂𝒍ab,[3b]where Πi is the hydrostatic cellular pressure and 𝑻ab=Tab𝒍^ab is the cell-edge tension, with 𝒍^ab=𝒍ab/|𝒍ab|. Here we use Roman indexes i,j,k,… to label cells and a,b,c,… to label vertices, while Greek indexes denote Cartesian components. The summation in Eq. 3a runs over all edges of cell i, and 𝒍ab is the edge vector joining vertices a and b when the perimeter of cell i is traversed clockwise, as shown in Fig. S1A. The factor of 1/2 in the second term on the right-hand side of Eq. 3a is due to the fact that each edge is shared by two cells. We have used the convention that the cellular stress is positive when the cell is contractile and negative when the cell is extensile. Contractile stress means that if the cell is cut off from its neighbors, it tends to contract, consistent with the contractility of the actomyosin network within the cell.

Note that in a vertex model the interaction stress as defined in Eq. 3a is indeed the stress acting on the tissue boundary. This is, however, not the case for the Voronoi model because the Voronoi construction introduces constraints not present in the vertex model. We have verified that the differences are small (SI Text) and in the following use Eq. 3a as a good approximation for the Voronoi model.

Our goal is to obtain the distribution of cellular stresses in a layer of motile cells, where cellular configurations do not minimize the tissue energy, but are governed by the dynamics described by Eq. 2. In this case, as discussed in the Introduction, the local cellular stress can be written as the sum of contributions from interactions and propulsive forces asσαβi=σαβ(i)int+σαβ(i)swim,[4]where, following recent work on active colloids (29, 30),σαβ(i)swim=−v0μAinαirβi[5]describes the flux of propulsive force v0μ𝒏i across a boundary, calculated from a virial expression. The negative sign in Eq. 5 ensures that the swim stress follows the same convention as the interaction stress, i.e., is positive for contractile stress and negative for extensile stress. The swim stress is proportional to the cell motility v0 and vanishes for nonmotile cells. The contribution σαβ(i)int is still given by Eq. 3a, but depends implicitly on cell motility through the instantaneous values of Πi and Tab that are determined not by energy minimization, but by the system dynamics governed by Eq. 2. From the local stresses one can then obtain the total mean stress in the tissue as (SI Text) σαβ=σαβint+σαβswim, withσαβint=1AT∑iAiσαβ(i)int,σαβswim=1AT∑iAiσαβ(i)swim.[6]In simulations where the energy functional is known, it is simple to directly extract the instantaneous pressures and tensions from cell shapes to calculate the interaction contribution σαβ(i)int. The definitions, Eq. 3b, giveΠi=−2KA(Ai−A0),Tab=2KP[(Pj−P0)+(Pk−P0)],[7]where ab is the cell–cell interface that separates cells j and k. Both Ai and Pi are obtained at every time step of the simulation and implicitly depend on cell motility, parameterized by speed v0 and persistence τr. This method directly infers the interaction stress from instantaneous cell shape fluctuations through Eq. 3a, yielding what we call shape-based stresses. While easily implemented numerically, this method is of limited use in experiments where the energy functional is not known and likely more complicated.

For this reason we develop a new mechanical inference method for motile monolayers that attempts to approximate the interaction stresses using only information that is accessible in experiments. Specifically, the proposed traction-based mechanical inference infers tensions and pressures from segmented images of cell boundaries and traction forces obtained by TFM. In the SPV model, we define the traction force at each vertex as the gradient of the tissue energy with respect to the vertex position 𝒕a=−∇aE, which balances the interaction force 𝑭a. Equilibrium mechanical inference methods express the interaction force 𝑭a=−∇aE at each vertex in terms of cellular pressures and edge tensions and the measured geometry of the cellular network (Eqs. S6–S8). Pressures and edge tensions are then obtained by inverting the equations 𝑭a({Πi},{Tab})=0. For a nonequilibrium epithelial layer of motile cells we invert the force balance equations𝑭a({Πi},{Ti})=𝒕a,[8]where the edge tensions Tab have been written as the sum of cortical tensions Ti of adjacent cells (Eq. S7), reducing the number of independent unknowns. The interaction contribution to the local cellular stresses is then calculated using Eq. 3a. A constraint counting yields 4N force balance equations for 2N variables, rendering the system overdetermined, which requires the implementation of a least-squares minimization for the mechanical inference (SI Text).

The equations developed thus far require knowledge of the tractions at each vertex, which is again not realistic in experiments. Therefore, we have developed and implemented a coarse-grained version of this approach that uses experimentally accessible traction forces averaged over a square grid, with a grid spacing of the order of a cell diameter. Pressures and tensions (Πi,Ti) are then calculated by inverting the force balance equation at the center of each grid element,𝑭grid({Πi},{Ti})=𝒕grid.[9]An outline of the coarse-graining procedure is given in SI Text. We refer to stresses inferred from Eq. 9 as traction-based stresses. We emphasize that the traction-based stress from mechanical inference is generally different from the intercellular stress obtained with monolayer stress microscopy (MSM), which rests upon the assumption that the tissue is an isotropic, homogeneous, and linearly elastic material (18). The mechanical inference does not make such assumptions and is compatible with any epithelia whose cell–cell interactions can be decomposed into tensions at cell junctions and pressures within cell bodies. Examples include Drosophila ectoderm (9) and mesoderm (25), Drosophila wing disk (26), and Madin–Darby canine kidney (MDCK) epithelium (12).

Within the framework of the SPV model, we have validated the coarse-grained method by showing that the resulting traction-based stresses agree with the shape-based stresses computed exactly from the simulations (Fig. 2 C and F).

Fig. 2.
Fig. 2.

Comparison of shape-based and coarse-grained traction-based stress. (A–C) Solid state at v0=0.5, P0=3.3. (D–F) Liquid state at v0=0.5, P0=3.8. (A and D) Interaction normal stress σn(i)int calculated from the instantaneous cell shapes obtained from Eqs. 3a and 7. Red denotes positive (contractile) stress and blue negative (extensile) stress. (B and E) Interaction normal stress σn(i)int calculated using the coarse-grained traction-based mechanical inference by inverting Eq. 9 and using Eq. 3a. The arrows denote the traction forces. (C and F) The coarse-grained traction-based mechanical inference is validated by plotting the traction-based stress against the shape-based stress in the solid (C) and in the liquid (F) state. The data are for 400 cells in a square box of side L=20 with Dr=0.1 and with periodic boundary condition.

Stress Characterizes Rheological Properties of the Tissue.

To study the mechanical properties of motile confluent tissues, we simulate a confluent cell layer with periodic boundary conditions, using the SPV model. By examining the temporal correlations of the mean stress in the tissue, as defined in Eq. 6, we show that the tissue displays distinct mechanical properties in the liquid and the solid states. Thus, mechanical measurements such as those provided by TFM can be used to characterize the rheological properties of the tissue.

The stress tensor σαβ is symmetric and has three independent components in two dimensions. Both the mean and local stress components are most usefully expressed in terms of normal stress σn, shear stress σs, and normal stress difference σd, withσn,d=12(σxx±σyy),σs=12(σxy+σyx).[10]Each of the interaction and swim contributions can similarly be split in normal, shear, and normal difference components. Below we focus on normal and shear stresses. TFM probes the forces exchanged between tissue and substrate, which by force balance are determined entirely by intercellular forces and hence by interaction stresses. In contrast, the swim components of stress and pressure cannot be probed in TFM, but contribute to the pressure Π=−σn=Πint+Πswim that the tissue would exert laterally on a confining piston. As we show below, the swim contribution dominates the pressure in the liquid state.

Using the expression for the local stress obtained from cell shapes, the mean interaction normal stress of the tissue can be expressed entirely in terms of area and perimeter fluctuations in a virial-like form (SI Text)σnint=1AT∑i[2KAAi(Ai−A0)+KPPi(Pi−P0)].[11]The first term represents the interaction contribution from the pressures within the cells. The second term is the contribution from the competition between actomyosin contractility and cell–cell adhesion that controls the cortical tensions. In our simulation, the cellular pressure is suppressed by setting A¯=A0 and the normal stress comes mainly from the cortical tensions. Eq. 11 then provides a way for extracting mechanical information directly from cell shape based on snapshots of segmented cell images.

Normal stresses are contractile in the solid phase and may become extensile deep in the liquid phase.

We show in Fig. 2 snapshots of the local interaction normal stress in the solid state (A–C) and in the liquid state (D–F). In both the solid and the liquid the interaction normal stress is on average contractile (red), with relatively weaker spatial fluctuations, but much larger mean value in the solid, where contractile cortical tension exceeds cell–cell adhesion.

Fig. 3 displays the total mean normal stress (the separate contributions from interaction and swim stress are shown in Fig. S2) across the solid–liquid transition. The color map shows that the total normal stress is contractile in the solid and across the transition line (Fig. 3A, black crosses), but changes sign and becomes extensile deep in the liquid. While the interaction stress is always positive due to cell contractility and consistent with experimental observations (1⇓–3, 38, the change in sign of the total stress is due to the swim stress that is zero in the solid and always negative in the liquid (Fig. S2), indicating that motility induces extensile stresses, tending to stretch the tissue. The total normal stress is analogous to the stress on a wall confining an active Brownian colloidal fluid (29, 30). We speculate that its change in sign could lead to an expansion of the tissue if released from confinement due to substrate patterning or to surrounding tissue and may contribute to epithelia expansion in wound-healing assays. In our model confinement is provided by the periodic boundary conditions.

Fig. 3.
Fig. 3.

Mean total normal stress. (A) Heat map of the mean total normal stress of the tissue in the (v0,P0) plane. The black crosses outline the solid–liquid phase boundary determined by q=3.813. Red indicates contractile stress and blue extensile stress. (B) Mean total normal stress as a function of P0 at various v0, showing a change in sign deep in the liquid state (400 cells for T=1,000 and Dr=0.1 with periodic boundary condition).

The tissue effective shear viscosity diverges at the liquid–solid transition.

While the local shear stress averages to zero in both the liquid and the solid states, its temporal correlations provide a distinctive rheological metric for distinguishing the liquid from the solid and identifying the transition. The time autocorrelation function of the interaction shear stress,Css(τ)=⟨σsint(t0)σsint(t0+τ)⟩t0,[12]where ⟨…⟩t0 denotes the average over the length t0 of the simulation, is shown in Fig. 4A for various v0 across the liquid–solid transition. Shear stress correlations decay in the liquid state and slow down as the transition is approached from the liquid side, becoming frozen in the solid. To quantify this we have defined the correlation time τm as the time when the correlation has decreased below 1% of its initial value. This correlation time shown in Fig. 4B diverges at the liquid–solid phase transition, suggesting that shear stress autocorrelations can provide a robust metric for the transition. Our work therefore suggests that TFM combined with mechanical inference can provide a tool for the measurement of tissue rheology. Correlating such measurements with cell shape data will provide a stringent test for our theory. Finally, in the liquid state we define an effective viscosity η𝑒𝑓𝑓, using a Green–Kubo-type relation by integrating the correlation function over the duration of the correlation time (41, 42),η𝑒𝑓𝑓=ATkbT𝑒𝑓𝑓∫0τmCss(τ)dτ,[13]where we have used the ideal gas effective temperature kBT𝑒𝑓𝑓=v02/(2μDr) so that η𝑒𝑓𝑓 has dimensions of a shear viscosity in two dimensions. Of course the Green–Kubo relation is based on the existence of a fluctuation–dissipation theorem, which does not hold in active systems. For small values of the persistence time τr, however, the orientational noise in the SPV becomes uncorrelated in time and can be mapped onto thermal noise at an effective temperature T𝑒𝑓𝑓. In this limit we expect η𝑒𝑓𝑓 to indeed play the role of a shear viscosity. Remarkably, the effective viscosity shown in Fig. 4 C and D grows as the tissue approaches the solid state from the liquid side and diverges at the transition. The effective viscosity quickly approaches zero deep in the liquid phase, suggesting that the system behaves as a gas of uncorrelated cells in this region.

Fig. 4.
Fig. 4.

Time autocorrelation of interaction shear stress and effective tissue viscosity. (A) Time autocorrelation function of the mean interaction shear stress for P0=3.45 at various v0. (B and D) Heat maps of the correlation time τm (B) and of the effective viscosity (D). The red crosses denote the liquid–solid phase transition boundary determined by q=3.813. The white region corresponds to regions where correlation time τm and effective viscosity η𝑒𝑓𝑓 diverge. Note that these regions largely coincide with the solid regime. (C) The effective viscosity as a function of P0 at various v0. The dashed lines correspond to the critical P0 where the liquid–solid phase transition occurs. All results are for 100 cells, Dr=1, and T=40,000 with periodic boundary condition.

Discussion and Conclusions

Using the SPV model we have formulated a unifying framework for quantifying the contributions from cell shape fluctuations and cell motility to mechanical stresses in an epithelial tissue. Cell shape fluctuations from actomyosin contractility and cell–cell adhesion control the interaction stress, while cell motility determines the swim stress that is generically present in all self-propelled systems (29, 30).

Unlike monolayer stress microscopy that computes interaction stress from traction forces by assuming the tissue to be a continuum linear elastic material (18), the traction-based mechanical inference method developed here incorporates spatial and temporal deformations of the tissue due to actomyosin contractility and cell–cell adhesion and can be generalized to account for cell division, apoptosis, and nematic/polar order of the tissue. In contrast to equilibrium mechanical inference techniques (21, 23), our approach does not require cells to be in or close to static mechanical balance, and it also provides the absolute scale of the junctional tensions and pressure differences. This can for instance be important for testing hypotheses involving mechanosensitive biomolecules. Experimentally, our method provides unique ways to extract intercellular interaction stresses from existing traction force data and segmented cell images.

The swim stress, on the other hand, cannot be measured using TFM as it represents the flux of propulsive forces across a bulk plane in the tissue. It contributes to the homeostatic pressure at the lateral boundary of the tissue. The sum of the swim stress and the interaction stress approximates the total stress at the tissue boundary, which is generally contractile but can become extensile when the tissue is deep in the liquid state and cell motility exceeds actomyosin contractility. The exact location of the transition from contractile to extensile depends on the average cellular pressure, which we fix by setting the average cell area to the preferred area, as discussed in SI Text. This change in sign may be observable in wound-healing assays where the transition from contractile to extensile behavior can result in tissue expansion upon removal of confinement by neighboring tissue.

We have extracted an effective tissue viscosity from the temporal correlation of the interaction shear stress. The correlation time and effective viscosity display a slowing down and arrest at the transition to the solid, thus serving as a direct probe of tissue rheology. Moreover, we observed a similar behavior for the temporal correlations of traction forces as demonstrated in SI Text. Therefore, our work suggests that TFM measurement (11) combined with mechanical inference could provide information on tissue rheology. To our knowledge, this has not been attempted yet on experimental data.

Our work sets the stage for examining the feedback between cell activity and tissue mechanics that is apparent in many tissue-level phenomena. Recent work has shown that mechanical stresses influence cell proliferation in tumor spheroids (28) and regulate cell growth in the developing Drosophila wing (43). Regulation of cell motility, as in contact inhibition of locomotion, has been proposed to explain stress patterns during collective cell migration (44). TFM has revealed the tendency of cells to move along the direction of minimal shear stress, a phenomenon termed “plithotaxis” (11). Our model provides a unifying framework for quantifying the relative roles of various cell properties, such as shape, motility, and growth, on the mechanics of the tissue.

Acknowledgments

We thank Daniel Sussman for valuable discussions. This work was supported by the Simons Foundation through Targeted Grant Award 342354 (to M.C.M. and M.C.) and Investigator Award 446222 (to M.L.M., M.C., and M.M.) in the Mathematical Modeling of Living Systems; by the National Science Foundation (NSF) through Awards DMR-1305184 (to M.C.M. and X.Y.), DMR-1609208 (to M.C.M.), DMR-1352184 (to M.L.M. and D.B.), and the Integrative Graduate Education and Research Traineeship (IGERT) Grant DGE-1068780 (to M.C.M. and M.C.); by the NIH through Grant R01GM117598-02 (to M.L.M.); and by the computational resources provided by Syracuse University and through NSF Grant ACI-1541396. All authors acknowledge support from the Syracuse University Soft Matter Program.

Footnotes

  • ↵1To whom correspondence should be addressed. Email: xingbo_yang{at}fas.harvard.edu.

  • Author contributions: X.Y., D.B., M.L.M., and M.C.M. designed research; X.Y., D.B., M.C., M.M., M.L.M., and M.C.M. performed research; X.Y., D.B., and M.C. analyzed data; and X.Y., D.B., M.L.M., and M.C.M. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1705921114/-/DCSupplemental.

Published under the PNAS license.

References

  1. ↵
    1. Lecuit T,
    2. Lenne PF
    (2007) Cell surface mechanics and the control of cell shape, tissue patterns and morphogenesis. Nat Mol Cell Biol 8:633–644.
    .
    OpenUrl
  2. ↵
    1. Heisenberg CP,
    2. Bellaiche Y
    (2013) Forces in tissue morphogenesis and patterning. Cell 153:948–962.
    .
    OpenUrlCrossRefPubMed
  3. ↵
    1. Murrell M,
    2. Oakes PW,
    3. Lenz M,
    4. Gardel ML
    (2015) Forcing cells into shape: The mechanics of actomyosin contractility. Nat Rev Mol Cell Biol 16:486–498.
    .
    OpenUrlCrossRefPubMed
  4. ↵
    1. Etournay R, et al.
    (2015) Interplay of cell dynamics and epithelial tension during morphogenesis of the drosophila pupal wing. eLife 4:e07090.
    .
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Shraiman BI
    (2005) Mechanical feedback as a possible regulator of tissue growth. Proc Natl Acad Sci USA 102:3318–3323.
    .
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Serra-Picamal X, et al.
    (2012) Mechanical waves during tissue expansion. Nat Phys 8:628–634.
    .
    OpenUrlCrossRef
  7. ↵
    1. Staple DB, et al.
    (2010) Mechanics and remodelling of cell packings in epithelia. Eur Phys J E Soft Matter 33:117–127.
    .
    OpenUrlCrossRefPubMed
  8. ↵
    1. Guillot C,
    2. Lecuit T
    (2013) Mechanics of epithelial tissue homeostasis and morphogenesis. Science 340:1185–1189.
    .
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Rauzi M,
    2. Verant P,
    3. Lecuit T,
    4. Lenne PF
    (2008) Nature and anisotropy of cortical forces orienting drosophila tissue morphogenesis. Nat Cell Biol 10:1401–1410.
    .
    OpenUrlCrossRefPubMed
  10. ↵
    1. Trepat X, et al.
    (2009) Physical forces during collective cell migration. Nat Phys 5:426–430.
    .
    OpenUrlCrossRef
  11. ↵
    1. Tambe DT, et al.
    (2011) Collective cell guidance by cooperative intercellular forces. Nat Mater 10:469–475.
    .
    OpenUrlCrossRefPubMed
  12. ↵
    1. Puliafito A, et al.
    (2012) Collective and single cell behavior in epithelial contact inhibition. Proc Natl Acad Sci USA 109:739–744.
    .
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Bi D,
    2. Lopez JH,
    3. Schwarz JM,
    4. Manning ML
    (2014) Energy barriers and cell migration in densely packed tissues. Soft Matter 10:1885–1890.
    .
    OpenUrlCrossRefPubMed
  14. ↵
    1. Bi D,
    2. Lopez JH,
    3. Schwarz JM,
    4. Manning ML
    (2015) A density-independent rigidity transition in biological tissues. Nat Phys 11:1074–1079.
    .
    OpenUrlCrossRef
  15. ↵
    1. Bi D,
    2. Yang X,
    3. Marchetti MC,
    4. Manning ML
    (2016) Motility-driven glass and jamming transitions in biological tissues. Phys Rev X 6:021011.
    .
    OpenUrlCrossRef
  16. ↵
    1. Park JA, et al.
    (2015) Unjamming and cell shape in the asthmatic airway epithelium. Nat Mater 14:1040–1048.
    .
    OpenUrlCrossRefPubMed
  17. ↵
    1. Soiné JRD, et al.
    (2014) Model-based traction force microscopy reveals differential tension in cellular actin bundles. PLoS Comput Biol 11:e1004076.
    .
    OpenUrl
  18. ↵
    1. Tambe DT, et al.
    (2013) Monolayer stress microscopy: Limitations, artifacts, and accuracy of recovered intercellular stresses. PLoS One 8:e55172.
    .
    OpenUrlCrossRefPubMed
  19. ↵
    1. Zimmermann J, et al.
    (2014) Intercellular stress reconstitution from traction force data. Biophys J 107:548–554.
    .
    OpenUrlCrossRefPubMed
  20. ↵
    1. Kim JH, et al.
    (2013) Propulsion and navigation within the advancing monolayer sheet. Nat Mater 12:856–863.
    .
    OpenUrlCrossRefPubMed
  21. ↵
    1. Ishihara S,
    2. Sugimura K
    (2012) Bayesian inference of force dynamics during morphogenesis. J Theor Biol 313:201–211.
    .
    OpenUrlCrossRefPubMed
  22. ↵
    1. Nier V, et al.
    (2016) Inference of internal stress in a cell monolayer. Biophys J 110:1625–1635.
    .
    OpenUrlCrossRefPubMed
  23. ↵
    1. Chiou KK,
    2. Hufnagel L,
    3. Shraiman BI
    (2012) Mechanical stress inference for two dimensional cell arrays. PLoS Comput Biol 8:e1002512.
    .
    OpenUrlCrossRefPubMed
  24. ↵
    1. Noll N,
    2. Mani M,
    3. Heemskerk I,
    4. Streichan S,
    5. Shraiman BI
    (2015) Active tension network model of epithelial mechanics. arXiv:1508.00623.
    .
  25. ↵
    1. Brodland GW, et al.
    (2010) Video force microscopy reveals the mechanics of ventral furrow invagination in Drosophila. Proc Natl Acad Sci USA 107:22111–22116.
    .
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Farhadifar R,
    2. Roper JC,
    3. Aigouy B,
    4. Eaton S,
    5. Julicher F
    (2007) The influence of cell mechanics, cell-cell interactions, and proliferation on epithelial packing. Curr Biol 17:2095–2104.
    .
    OpenUrlCrossRefPubMed
  27. ↵
    1. Basan M,
    2. Risler T,
    3. Joanny JF,
    4. Sastre-Garau X,
    5. Prost J
    (2009) Homeostatic competition drives tumor growth and metastasis nucleation. HFSP J 3:265–272.
    .
    OpenUrlCrossRefPubMed
  28. ↵
    1. Montel F, et al.
    (2011) Stress clamp experiments on multicellular tumor spheroids. Phys Rev Lett 107:188102.
    .
    OpenUrlCrossRefPubMed
  29. ↵
    1. Yang X,
    2. Manning ML,
    3. Marchetti MC
    (2014) Aggregation and segregation of confined active particles. Soft Matter 10:6477–6484.
    .
    OpenUrlCrossRefPubMed
  30. ↵
    1. Takatori SC,
    2. Yan W,
    3. Brady JF
    (2014) Swim pressure: Stress generation in active matter. Phys Rev Lett 113:028103.
    .
    OpenUrl
  31. ↵
    1. Solon AP, et al.
    (2015) Pressure is not a state function for generic active fluids. Nat Phys 11:673–678.
    .
    OpenUrlCrossRef
  32. ↵
    1. Vazquez-Faci T,
    2. van Drongelen R,
    3. van der Zee M,
    4. Idema T
    (2017) Mechanics of epithelial tissue formation in early insect embryos. arXiv:1705.06205.
    .
  33. ↵
    1. Barton DL,
    2. Henkes S,
    3. Weijer CJ,
    4. Sknepnek R
    (2016) Active vertex model for cell-resolution description of epithelial tissue mechanics. PLoS Comput Biol 13:e1005569.
    .
    OpenUrl
  34. ↵
    1. Su T,
    2. Lan G
    (2016) Overcrowding drives the unjamming transition of gap-free monolayers. arXiv:1610.04254.
    .
  35. ↵
    1. Henkes S,
    2. Fily Y,
    3. Marchetti MC
    (2011) Active jamming: Self-propelled soft particles at high density. Phys Rev E 84:040301.
    .
    OpenUrl
  36. ↵
    1. Fily Y,
    2. Marchetti MC
    (2012) Athermal phase separation of self-propelled particles with no alignment. Phys Rev Lett 108:235702.
    .
    OpenUrlCrossRefPubMed
  37. ↵
    1. Nestor-Bergmann A,
    2. Goddard G,
    3. Woolner S,
    4. Jensen O
    (2016) A vertex-based model relating cell shape and mechanical stress in an epithelium. arXiv:1611.04744.
    .
  38. ↵
    1. Bambardekar K,
    2. Clement R,
    3. Blanc O,
    4. Chardes C,
    5. Lenne PF
    (2015) Direct laser manipulation reveals the mechanics of cell contacts in vivo. Proc Natl Acad Sci USA 112:1416–1421.
    .
    OpenUrlAbstract/FREE Full Text
    1. Rauzzi M,
    2. Lenne PF,
    3. Lecuit T
    (2010) Planar polarized actomyosin contractile flows control epithelial junction remodelling. Nature 23:1110–1114.
    .
    OpenUrl
    1. Kasza KE,
    2. Farrell DL,
    3. Zallen JA
    (2014) Spatiotemporal control of epithelial remodeling by regulated myosin phosphorylation. Proc Natl Acad Sci USA 111:11732–11737.
    .
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Hansen JP,
    2. McDonald I
    (2006) Theory of Simple Liquids (Academic, Cambridge, MA), 3rd Ed.
    .
  42. ↵
    1. Levashov VA,
    2. Morris JR,
    3. Egami T
    (2011) Viscosity, shearwaves, and atomic-level stress-stress correlations. Phys Rev Lett 106:115703.
    .
    OpenUrlPubMed
  43. ↵
    1. Pan Y,
    2. Heemskerk I,
    3. Ibar C,
    4. Shraiman BI,
    5. Irvine KD
    (2016) Differential growth triggers mechanical feedback that elevates hippo signaling. Proc Natl Acad Sci USA 113:E6974–E6983.
    .
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Zimmermann J,
    2. Camley BA,
    3. Rappel WJ,
    4. Levine H
    (2016) Contact inhibition of locomotion determines cell-cell and cell-substrate forces in tissues. Proc Natl Acad Sci USA 113:2660–2665.
    .
    OpenUrlAbstract/FREE Full Text
  45. MathWorks (2017) Matlab R2017b documentation: Moore-Penrose pseudoinverse. https://www.mathworks.com/help/matlab/ref/pinv.html. Accessed November 7, 2017.
    .
View Abstract
Back to top
Article Alerts
Email Article
Citation Tools
Request Permissions
Correlating cell shape and cellular stress
Xingbo Yang, Dapeng Bi, Michael Czajkowski, Matthias Merkel, M. Lisa Manning, M. Cristina Marchetti
Proceedings of the National Academy of Sciences Nov 2017, 114 (48) 12663-12668; DOI: 10.1073/pnas.1705921114
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo

You May Also be Interested in

Researchers report biparental inheritance of mitochondrial DNA in 17 members of three unrelated multigeneration families, paving the way for insights into alternative mechanisms for the treatment of inherited mitochondrial diseases. Image courtesy of Pixabay/PublicDomainPictures.
Human mitochondrial DNA can be paternally inherited
Researchers report biparental inheritance of mitochondrial DNA in 17 members of three unrelated multigeneration families, paving the way for insights into alternative mechanisms for the treatment of inherited mitochondrial diseases.
Image courtesy of Pixabay/PublicDomainPictures.
Researchers report a machine-learning approach to identify land plants at risk of extinction, suggesting that the approach can be used to guide policies aimed at allocating resources for biodiversity conservation. Image courtesy of Abbie Zimmer (Ohio State University, Columbus, OH).
Machine learning identifies plants at risk of extinction
Researchers report a machine-learning approach to identify land plants at risk of extinction, suggesting that the approach can be used to guide policies aimed at allocating resources for biodiversity conservation.
Image courtesy of Abbie Zimmer (Ohio State University, Columbus, OH).
A study explores how cats groom fur using fine structures called papillae on the surface of the tongue and presents a biologically inspired hairbrush to remove allergens from cat fur and apply medications on cat skin. Image courtesy of Candler Hobbs (Georgia Institute of Technology, Atlanta).
Fine analysis of fur grooming in cats
A study explores how cats groom fur using fine structures called papillae on the surface of the tongue and presents a biologically inspired hairbrush to remove allergens from cat fur and apply medications on cat skin.
Image courtesy of Candler Hobbs (Georgia Institute of Technology, Atlanta).
Proceedings of the National Academy of Sciences: 115 (50)
Current Issue

Submit

Jump to section