Unleashing shear: Role of intercellular traction and cellular moments in collective cell migration.

In the field of endothelial biology, the term "shear forces" is tied to the forces exerted by the flowing blood on the quiescent cells. But endothelial cells themselves also exert physical forces on their immediate and distant neighbors. Specific factors of such intrinsic mechanical signals most relevant to immediate neighbors include normal (Fn) and shear (Fs) components of intercellular tractions, and those factors most relevant to distant neighbors include contractile or dilatational (Mc) and shear (Ms) components of the moments of cytoskeletal forces. However, for cells within a monolayer, Fn, Fs, Mc, and Ms remain inaccessible to experimental evaluation. Here, we present an approach that enables quantitative assessment of these properties. Remarkably, across a collectively migrating sheet of pulmonary microvascular endothelial cells, Fs was of the same order of magnitude as Fn. Moreover, compared to the normal components (Fn, Mc) of the mechanical signals, the shear components (Fs, Ms) were more distinctive in the cells closer to the migration front. Individual cells had an innately collective tendency to migrate along the axis of maximum contractile moment - a collective migratory process we referred to as cellular plithotaxis. Notably, larger Fs and Ms were associated with stronger plithotaxis, but dilatational moment appeared to disengage plithotactic guidance. Overall, cellular plithotaxis was more strongly associated with the "shear forces" (Fs, Ms) than with the "normal forces" (Fn, Mc). Finally, the mechanical state of the cells with fast migration speed and those with highly circular shape were reminiscent of fluid-like and solid-like matter, respectively. The results repeatedly pointed to neighbors imposing shear forces on a cell as a highly significant event, and hence, the term "shear forces" must include not just the forces from flowing fluid but also the forces from the substrate and neighbors. Collectively, these advances set the stage for deeper understanding of mechanical signaling in cellular monolayers.

Introduction

Endothelial cells possess the cytoskeletal apparatus that is capable of generating and transmitting contractile forces [1-3]. Across an endothelial monolayer, such forces propagate from each cell to its substrate and neighbors [4]. Such force propagation induces tractions across the cell-substrate and cell-cell interfaces. From a mechanical viewpoint, the distribution of tractions across the cell-substrate and cell-cell interfaces will have a strong influence on immediate neighbors, whereas resultant magnitude and moment of the tractions will have a strong influence on distant neighbors. A comprehensive assessment of tractions involves measurement of their magnitude and the direction of their action. Currently, however, such comprehensive assessment is only possible for the substrate tractions and not for the intercellular tractions and the moments of these tractions [5]. The inability to quantify such basic mechanical properties has significantly narrowed the scope of the studies focused on mechanical signaling processes in cellular monolayers.

The substrate and intercellular tractions induce mechanical stresses across the monolayer. Monolayer Stress Microscopy (MSM) enabled a comprehensive assessment of thickness averaged monolayer stresses across such a monolayer [6,7]. At any subcellular region within the monolayer, MSM can identify maximum principal orientation along which normal stress (tensile or compressive stress) is highest and shear stress is zero [7]. Such orientation is central to the mechanical regulation of cellular motion, reorientation, and intercellular gap formation [6,8,9]. For example, in an advancing cellular monolayer, local cellular motion tends to be aligned with the local maximum principal orientation [6]. However, this phenomenon, described as plithotaxis, was identified by comparing motion and stresses defined over subcellular regions. The validity of plithotaxis has not been tested for maximum principal orientation and motion defined over an entire cell. Moreover, the physical significance of intercellular tractions in plithotactic guidance remains unexplored.

To address these knowledge gaps, we focused on an advancing sheet of pulmonary microvascular endothelial cells. We then described the approach that enables quantitative assessment of intercellular tractions and cellular moments [10]. To identify the physical significance of the intercellular tractions and cellular moments, we examined their distribution across the cellular monolayer and their relationship with cellular shape, speed, and plithotaxis. Surprisingly, compared to the normal components of intercellular forces and cellular moments, their shear counterparts were more strongly related to cellular shape, speed, and plithotaxis. Collectively, these findings open promising avenues for future research focused on the significance of fluid-like and solid-like behavior in the mechanics of cell monolayer.

Materials and methods

Cell culture and image acquisition

The cell culture procedure is the same as described in an earlier publication [11]. Briefly, rat pulmonary microvascular endothelial cells (rat 1, passage 11) were acquired from the cell culture core of the Center for Lung Biology at the University of South Alabama and cultured in 10% fetal bovine serum in a standard tissue culture environment. The data was acquired from cellular passages 12 through 16. Images were acquired using an inverted wide-field fluorescence microscope (Leica, DMI 600B) and confocal microscope (Nikon A1R). The data sets include three long (more than 940 min) and five short (more than 30 min) time-lapse sequences, each acquired at a 5 min interval.

Polyacrylamide hydrogel

The hydrogel preparation procedure was also the same as described in an earlier publication [11]. Briefly, all experiments were conducted in 35 mm glass-bottom dishes containing collagen-coated polyacrylamide hydrogels (1250 Pa shear modulus, approximately 100 μm thickness, 0.5 μm diameter embedded fluorescent beads).

Assessment of subcellular forces and motion

The subcellular forces and motion were also assessed as described in an earlier publication [11]. Briefly, components of substrate traction (T) were quantified using Fourier Transform Traction Microscopy (FTTM), monolayer stresses (σ) were quantified using Monolayer Stress Microscopy (MSM) and the cellular motion was quantified using Particle Image Velocimetry (PIV) [5-7]. Since the substrate tractions and monolayer stresses are two-dimensional, the suffix i and j can be either x or y. The cross-correlation window size and spacing between adjacent windows were respectively 15.8 μm and 7.9 μm for PIV data analysis, and respectively 7.9 μm and 2.6 μm for both FTTM and MSM data analysis.

Assessment of intercellular tractions and cellular moments

The location of intercellular boundaries was detected from a phase-contrast image of the cellular monolayer using a custom image segmentation software (Fig. 1a) [11]. Through the pixels that defined intercellular boundary, we fitted a Fourier series with eight harmonics. This curve fitting enabled straightforward computation of unit normal (n) and tangent (t) vectors along the intercellular boundary (Fig. 1a). Normal (Fn) and shear (Fs) components of intercellular traction were computed using equations, Fn = σnn and Fs = σnt (Fig. 1b). Although these traction components were defined on a plane that is perpendicular to the monolayer, this approach can be applied to any angle between the plane that separates adjoining cells.

Fig. 1.

Approach to quantify intercellular tractions and cellular moments, and distribution of these properties across the monolayer.

a. Image of the cellular monolayer with intercellular boundaries (blue lines) identified using a custom segmentation software. The middle sub-panel shows pattern of sampling points where MSM provides T, and σ. The bottom sub-panel shows unit vectors along normal (red arrows) and tangential (blue arrows) orientation. b. Fn (red arrows), and Fs (blue arrows), for a cell within the monolayer. c. Image of the cellular monolayer that illustrates including Mmax (green arrow), Mmin (yellow arrow), v (magenta arrow), cellular orientation (white dashed line), and angles θMv, θCv, and θCM (Table 1). The cumulative probability distribution of d. cellular spread area, e. cellular circularity, f. TRMS, g. , h. , i. , j. Mc, k. Ms, l. vRMS. All cumulative distributions are plotted for quintiles of the distance of the centroid of the cell from the migration front. The lowest quintile is shown in blue and progressively higher quintiles are shown in progressively warmer colors. Arrow in panel d indicates progression from a lower quintile to a higher quintile.

Components of the cellular moment were computed using , where r was the x or y coordinate where tractions were quantified and is the area and is the boundary of the cell as observed in the phase-contrast image. This definition of moment provides a positive number when moments are contractile representing tension in the cell, and a negative number when moments are dilatational representing compression of the cell. To compute principal moments Mmax and Mmin, we diagonalized the moment tensor using , where p and q are the eigenvectors.

The orientation of maximum principal moment, Mmax, was computed using θM = arctan(p/p). The contractile or dilatational component of the cellular moment was computed using Mc = (Mmax + Mmin)/2. Maximum shear or anisotropy of the cellular moment was computed using Ms = (Mmax − Mmin)/2. It is important to note that for the lowest quintile of Mc, the median value was dilatational; but for the top four quintiles, the median values were contractile.

Physical properties considered in the data analysis

The circularity of the cell was computed using the formula, Circularity = 4πArea/Perimeter2. To examine the significance of intercellular traction, we defined a ratio between its two components, R = |Fs|/Fn. This ratio assesses how other physical properties of individual cells are associated with (1) the relative magnitude of Fn and Fs, and (2) the orientation of intercellular traction. For each of the examined physical properties, we generated cumulative probability distribution from 300 randomly chosen frames each containing close to 300 cells.

Results

Distribution of tractions, moments, and other properties across the cellular monolayer

Compared to cells far from the migration front, cells near the migration front had several distinctive properties. Their morphological distinction included larger spread area and slightly larger circularity (Fig. 1d,e). Their distinctive mechanical state included larger TRMS, , Ms, and vRMS (Fig. 1f,h,k,l and Table 1). However, these cells did not have a distinctive , , RSD, or Mc (Figs. 1g,i,j and S1, and Table 1). Finally, the influence of cellular position with respect to the migration front was least striking for cellular circularity but most striking for vRMS (Fig. 1d,e).

Table 1

List of commonly used cellular properties and mathematical operators.

Property	Description (units)
P	Cumulative probability distribution (%)
vi	Components of the average velocity of a cell (μm/min)
Ti	Components of the tractions exerted by the cell on the substrate (Pa)
σij	Components of monolayer stresses (Pa)
Fs	Intercellular traction component oriented parallel to the intercellular interface (Pa)
Fn	Intercellular traction component oriented perpendicular to the intercellular interface (Pa)
R	Ratio of intercellular traction components, |Fs|/Fn
Ms	Shear moment or moment anisotropy of the cell (pN-m)
Mc	Contractile or dilatational moment of the cell (pN-m)
Mmax	Maximum principal moment of the cell (pN-m)
θMv	Angle between the orientation of Mmax of the cell and the axis of net cellular velocity (°)
θCv	Angle between the orientation of the cell and the axis of net cellular velocity (°)
θCM	Angle between the orientation of the cell and the orientation of Mmax of the cell (°)
(…)RMS	Root mean square value across the cell
|…|	Absolute value of the local property
(…)SD	Standard deviation across the cell
(…˜)	Median value across the cell

Relationships of intercellular tractions and cellular moments with cellular shape and speed

and were weakly related to cellular circularity and speed (Fig. 2a–d). Both and were larger in the cells that had lower circularity and faster speed. and were independent of cellular spread area, but Mc and Ms were larger for the cells with bigger spread area (Fig. S2).

Fig. 2.

Relationship of intercellular tractions and cellular moments with cellular shape and speed.

The cumulative distribution of a, b. , and c, d. , e, f. Mc/Area, and g, h. Ms/Area. Cumulative distributions in the top row of panels are for quintiles of cellular circularity and those in the bottom row of panels are for quintiles of vRMS. The lowest quintile is shown in blue and progressively higher quintiles are shown in progressively warmer colors.

Cellular circularity and speed were not as strongly related to Mc as they were with the density of moment defined as Mc/Area (Figs. 2e,f and S3a,b, and Table 1). Similar enhancement of relationship was observed in Ms/Area versus cellular circularity but was not observed in Ms/Area versus cellular speed (Figs. 2g,h and S3c,d). Mc/Area and Ms/Area were larger in the cells with lower circularity (Fig. 2e,g). Relationship between Mc/Area and vRMS was opposite to that between Ms/Area and vRMS (Fig. 2f,h). In that, faster cells had smaller Mc/Area but larger Ms/Area.

Relationships of intercellular traction orientation with cellular speed and moments

The distribution of was almost independent of vRMS and top four quintiles of Mc, but systematically larger for higher quintiles of Ms (Fig. 3a–c). The lowest quintile of Mc which had a dilatational median value and large contributions from compressive Fn showed a distinct relationship (Fig. 3b, blue curve). Unlike the relationship with , vRMS and Mc had a stronger relationship with RSD (Fig. 3d,e). By contrast, Ms was more strongly related with than with RSD (Fig. 3c,f). Interestingly, Mc and Ms had an opposite relationship with RSD (Fig. 3e,f).

Fig. 3.

Significance of median value and variability in intercellular traction orientation.

The cumulative probability distribution of , for quintiles of a. vRMS, b. Mc, and c. Ms. The cumulative probability distribution of RSD for quintiles of d. vRMS, e. Mc and f. Ms. The lowest quintile is shown in blue and progressively higher quintiles are shown in progressively warmer colors.

Role of tractions, cellular moments, cellular speed, and cellular orientation in cellular plithotaxis

Smaller θMv was associated with larger , , , Mc, and Ms (Fig. 4a–c,e,f, and Table 1). The relationship was strongest for (Fig. 4b). In contrast, the relationship was absent in the cells with compressive or dilatational Mc (Fig. 4a,e, blue curve).

Fig. 4.

Role of intercellular tractions, cellular moments, orientation, speed and traction in cellular plithotaxis.

The cumulative probability distribution of θMv for quintiles of a. , b. , c. , d. TRMS, e. Mc, f. Ms, and g. vRMS. The cumulative probability distribution of θCv for quintiles of h. θMv. The lowest quintile is shown in blue and progressively higher quintiles are shown in progressively warmer colors.

Unlike the intercellular tractions, the substrate tractions did not have a systematic relationship with plithotaxis (Fig. 4a–d). Faster cells were associated with stronger plithotaxis, but stronger plithotaxis was not associated with smaller θCv (Fig. 4g,h). However, cellular alignment with the axis of the principal moment and the axis of cellular motion appeared to be strongly proportional (Fig. S4).

Discussion

Visualizing fundamental intrinsic mechanical signals within individual cells of a monolayers

Both intrinsic and extrinsic mechanical signals induce structural and biochemical changes in individual cells of an adherent endothelial monolayer [12]. Early quantitative assessments of intrinsic mechanical signals involved measurement of isometric tension across the cellular monolayer [13]. Application of FTTM to cellular monolayer enabled quantitative assessment of the local adhesive forces transmitted to the extracellular matrix [5]. MSM enabled quantitative assessment of local mechanical stresses across the monolayer [6]. Scalar changes in tension across specific adhesion molecules, and net moments within an isolated individual cell, a doublet, and a small island have been measured [14-19]. However, for individual cells within a monolayer, two fundamental intrinsic signals remained hidden. First, individual components of the adhesive forces across the intercellular junction (Fn and Fs), and second, individual components of the moments across the cell (Mc and Ms). The approach described here quantifies these fundamental mechanical signals and enables the application of the principles of elasticity to better understand the mechanical interactions between immediate and distant neighbors.

The intercellular shear traction is a dominant mechanical signal

Substrate tractions were higher near the migration front than those far from it (Fig. 1f). Such traction distribution is observed in a variety of cell types [6]. A widely appreciated interpretation comes from a comparison of these tractions with force-balance in a tug-of-war [5]. But such comparison is made by adding tractions along a particular coordinate axis. The interpretation suggests that the monolayer is in a state of a “global” tug-of-war where the tension increases continuously from the advancing front to the center of the monolayer [5]. Importantly, in this interpretation, monolayer shear stresses along the same coordinate axis are assumed to be negligible.

While such interpretation is useful to appreciate the existence of intercellular force transmission, the inference that monolayer tension increases with the distance from the migration front can be misleading. In contrast with such inference, a rigorous force balance across the cellular monolayer revealed rather occasional regions – each spanning tens of cells – engaged in a “local” tug-of-war [6]. A striking contrast with the “global” tug-of-war model can be observed by comparing distributions of TRMS, , and Mc as a function of distance from the migration front (Fig. 1f,g,j). While TRMS was systematically larger closer to the migration front, both and Mc had no systematic gradients. Moreover, compared to Fn, Fs was often smaller but far from negligible (Fig. 1i). Remarkably, was larger for the cells located closer to the migration front (Fig. 1h). Such distribution revels the key source of the misleading inference form the “global” tug-of-war model – large portion of the interpreted tension is balanced closer to the migration front.

The cellular shear moment is also a dominant mechanical signal

Ms was of a similar order of magnitude as Mc (Fig. 1j,k). Although both Mc and Ms have contributions from all traction components, the relationship between and Ms was strongest (Fig. S5). In that the cells within the lowest quintile of had negligible Ms and the cells within the highest quintile of had a much larger Ms (Fig. S5e).

The mechanical state of highly circular cells includes signatures of a solid-like behavior

The mechanical state of the cells with higher circularity included lower , , Mc/Area, and Ms/Area (Fig. 2a,c,e,g). Interestingly, the cellular circularity was independent of RSD, θCM, and vRMS, and higher circularity was associated with weaker plithotaxis (Fig. S6). Coincidentally, cellular circularity increasing and approaching the value of 0.866 is a structural signature that the monolayer approaching a jamming transition, where individual cells are able to fluctuate but unable to change places with immediate neighbors [20]. Increased cellular density makes cellular motion structurally correlated over a wider area, an indicator of an approach to jamming transition [6]. Collectively, these properties are reminiscent of the particles of a jammed matter that withstands small shear stress [21]. Overall, the mechanical state of the cells with higher circularity indicates the need to use the approach described here combined with assessment of the cooperativity of cellular motion, and shear modulus of the cells to critically examine their solid-like behavior [6,22].

The mechanical state of fast-moving cells includes signatures of a fluid-like behavior

The mechanical state of the fast-moving cells included higher , and Ms/Area (Fig. 2d,h). These properties are reminiscent of an external shear stress-driven motion of a viscous fluid. In addition, fast-moving cells are located closer to the migration front where cellular density is low (Fig. 1d,l). Collectively, these properties are reminiscent of the particles of a jammed matter gaining the ability to change places with immediate neighbors when exposed to high shear stress [6,21]. In addition, the association of faster motion with lower Mc/Area is reminiscent of the cells that generate lower contractile forces due to cytoskeletal fluidization (Fig. 2f) [23]. Overall, the mechanical state of fast-moving cells indicates the need to use the approach described here combined with visualization of subcellular motion and cellular contractile apparatus including cytoskeletal architecture to critically examine their fluid-like behavior [24,25].

The potential use of intercellular traction orientation to develop predictive models for monolayers

RSD distinguished vRMS and Mc more strongly than (Fig. 3a,b,d,e). By contrast, distinguished Ms more strongly than RSD (Fig. 3c,f). As such, a mechanical or biochemical intervention that increases RSD will tend to increase vRMS, reduce Mc, and increase Ms. But an intervention that does not increase RSD but instead increase will tend to have little effect on vRMS and Mc, but it will tend to increase Ms. These findings and expectations can guide the development of theoretical models that can predict the mechanical state of the cell in response to mechanical and biochemical interventions or predict the behavior of a cell that is in a homogeneous versus heterogeneous microenvironment [11].

Cellular plithotaxis

The key observation of plithotaxis is that in a crowded environment of a monolayer, each constituent cell tends to migrate along the local axis of maximum principal stress. Plithotactic tendency is strongest at the locations where the anisotropy of the monolayer stress, or equivalently, the local maximum shear stress is highest [6]. However, along the axis of maximum principal stress, shear stress is necessarily zero. Hence, in plithotaxis, the collective significance of the role of zero shear stress along the preferred axis of motion and large maximum shear stress at the mechanically-steered locations remain unclear. One key question is whether the intercellular interfaces coincide with local axis of zero shear stress. The source of such knowledge gaps is that the picture of plithotaxis is identified using local subcellular analysis [6]. But access to cellular moments enabled us to reframe plithotaxis as a property of individual cells. In this process, we replaced the axis of maximum principal stress with the axis of Mmax, and also replaced the orientation of local velocity vector with the orientation of v. This cellular plithotaxis is an innately collective tendency of individual cells in a monolayer to migrate along the axis of maximum contractile moment. Mapping Fs showed that the intercellular interfaces rarely coincided with the local axis of zero shear stress. On the contrary, larger was associated with stronger plithotaxis (Fig. 4c).

Cellular plithotaxis revealed additional useful insights. While cellular plithotaxis was associated with intercellular tractions, it was not associated with substrate tractions (Fig. 4a–d). Cells exhibiting stronger plithotaxis moved faster and had lower circularity (Figs.4g and S6d). However, the faster motion was not related to cellular alignment with either the axis of motion or the axis of the maximum principal moment (Fig. S7). Finally, plithotaxis was not associated with cells moving along their long axis, but cellular alignment along the maximum principal moment and cellular alignment with velocity went hand-in-hand (Figs. 4h and S4).

Unleashing shear

Cells in an advancing monolayer tended to generate high intrinsic “shear forces” (Fs and Ms). These shear forces were strongly associated with cellular plithotaxis and tended to increase cellular speed. In a way, individual cells within a monolayer are unleashed, i.e., they move faster, when neighbors unleash, i.e., impose, shear forces on them. In addition to such unleashing of and by shear forces, the data suggest the need to unleash the term “shear forces”, i.e., expand its scope. In the field of endothelial biology, the term “shear forces” has long been used exclusively for the forces exerted by viscous fluid flowing over the cells. However, the current study suggests making that term inclusive of intercellular shear and apply it to guide our understanding of intercellular interactions.

The framework presented here can be used to examine how the relationships described here hold in other cellular systems or in quiescent monolayers, and how the three shear stresses – fluid shear stress, cell-substrate tractions, and intercellular shear traction – compete in disrupting or restoring endothelial homeostasis.