Main

The single adherent cell moves by the action of two synchronized cycles, one involving extension and contraction of its cytoskeleton and the other involving formation and detachment of its adhesions13,14. Although this complex process remains a matter of intense research14,15,16, it is now well established that a fundamental aspect of the motility mechanism is the transmission of contractile forces to the surrounding matrix at the cell’s leading and trailing edges17,18. In contrast with the case of migration of the single cell studied in isolation14,15,16, the case of collective migration of cells within a contiguous cell sheet has more physiological relevance but is substantially less well understood19. Within an advancing epithelial cell sheet, for example, each individual cell is physically constrained by its neighbours, and cell–cell signalling through biochemical and biophysical pathways may influence the collective motion of the group20,21. Do leader cells at the advancing front edge of the sheet exert physical forces locally that are transmitted rearward, from cell-to-cell, and thus act to pull along those cells in the ranks behind5,6,8,9? Or instead is each individual cell in the sheet mechanically self-propelled21? Or does cell proliferation expand the cell colony and thereby push the advancing front forward? Or is the correct answer none of the above? For more than a century these fundamental questions have been debated intensively5,22 and, using a variety of methods in vivo23, in vitro4,9,24 and in silico21, much conflicting evidence has accumulated. This conflicting evidence has been in most cases indirect or inferential, however, because within the cell sheet the physical forces themselves have remained largely inaccessible to direct experimental observation.

Here, we report by direct measurement the first explicit maps of those physical forces and their distribution. To do this within an advancing cell sheet, we used Fourier-transform traction microscopy together with a balance of forces that is demanded by straightforward application of Newton’s laws of motion (see Supplementary Information S2 and S3). To address the case of an advancing cell sheet, however, traction microscopy as described originally18,25 or as modified subsequently26,27,28 is inadequate and therefore required fundamental reformulation (see Supplementary Information S2 and S3). We seeded a small number of Madin–Darby canine kidney epithelial cells (5,000) at the centre of a soft collagen-coated polyacrylamide gel (Young’s modulus of 1.3 kPa). The cells adhered readily to this substrate and within 24 h formed a confluent colony. With time, the colony expanded radially outward, thus providing a simple model of collective migration without need of damaging the monolayer as in classical scratch-wounding experiments (see Supplementary Information S1). Growth of the colony was largely insensitive to the stiffness of the underlying substrate (see Supplementary Information S4). After allowing the colony to expand for at least 72 h, we mapped the traction forces that marginal and submarginal cells exerted on their underlying matrix.

We first assessed the locus of traction forces in the proximity of the leading edge (Fig. 1, Supplementary Movie S1). Maps of tractions normal () and parallel (T) to the leading edge show that tractions are not restricted to cells at the leading edge or even to cells located 2–3 rows behind it, as is commonly emphasized5,6,29. Instead, large tractions are applied by cells many cell rows behind the edge. Independent of the distance from the edge, both and T exhibited broad non-Gaussian distributions characterized by exponential tails (Fig. 2a, b). The distribution of was skewed towards positive tractions at the leading edge, whereas the distribution of T was symmetric. Both traction distributions narrowed as the distance from the leading edge increased. Taken together, these data are inconsistent with the existence of two populations of cells, each with a distinct mechanical phenotype, one corresponding to active mechanical leaders at the leading edge and the other to passive mechanical followers. Instead, our data show a single distribution, the tails of which were roughly exponential rather than Gaussian, revealing probabilities of high tractions much larger than would be predicted according to the central limit theorem for independent and identically distributed random variables. Exponential distributions have previously been reported at the level of the single focal adhesion30,31,32, thereby indicating that this particular kind of distribution might underlie tissue behaviour over multiple length scales.

Figure 1: Traction forces generated by a collectively migrating cell sheet.
figure 1

a, Phase contrast image. b, Tractions normal to the edge. c, Tractions parallel to the edge. The field of view is 750 μm×750 μm. T and were calculated from Tx and Ty and from the local normal vector to the cell edge (see Supplementary Methods).

Figure 2: Traction force distributions at different distances from the leading edge.
figure 2

a, Normal tractions. b, Parallel tractions. Red circles: first cell row from the leading edge; green circles: row 2; blue circles: rows 3–5; purple circles: rows 6–11; orange circles: rows 12–17. For computational simplicity, each row was assumed to be 19.2 μm in radial dimension. Data were pooled from n=4 different cell sheets at four different time points for each well. The tails of each distribution appear straight in a semilog plot, showing the exponential nature of the distributions. A Gaussian fit to the distribution of parallel tractions for rows 12–17 is plotted as a reference (dashed grey line). c, The average normal traction decayed slowly with distance from the edge (filled symbols), whereas the average parallel traction was negligible and independent of the distance from the edge (open symbols). Error bars indicate standard errors. d, Stress within the cell sheet increased as a function of the distance from the leading edge. Error bars indicate standard errors. e, Schematic diagram illustrating the computation of stress within the cell sheet. The average stress 〈σx x〉 normal to a plane perpendicular to the substrate and parallel to the leading edge can be calculated by integration of tractions Tx between the edge and the plane.

Submarginal cells have previously been shown to extend cryptic lamellipodia beneath cells in front of them24. Regardless of the extent to which cryptic lamellipodia are mechanically active and represent a locus of force generation, traction forces generated by these submarginal cells are seen to be comparable to those at the leading edge. A more important question, however, is whether these forces are balanced locally, as if each cell were self-propelled, or whether they are transmitted from cell to cell in a cooperative manner. We note that no amount of kinematic data or structural data, no matter how detailed, and no molecular manipulation, no matter how sophisticated, can ever suffice to resolve this question. Instead, we provide here a conclusive answer based solely on the application of Newton’s laws. We begin by computing the spatial averages (denoted by 〈 〉) of and T (Fig. 2c). As expected by symmetry, 〈T〉 was near zero. was maximum at the leading edge and progressively decayed with distance from the edge before changing sign in the opposite half of the sheet. In contrast to 〈T〉, however, the spatial average far away from the edge remained weakly but systematically greater than zero, with typical tugging tractions of the order of 5 Pa. In the context of cell mechanics, a regional stress of magnitude in this range is often regarded as being small, and is certainly small compared with the fluctuations that we observed. One might therefore conclude that its effects are essentially negligible; we come to a quite different conclusion, however.

Within the field of measurement—here spanning less than half the sheet diameter—the average stress at the interface between the cell base and the cell substrate, , was regionally positive everywhere. This stress acts systematically in a direction that pulls the sheet towards the leading edge. The question then arises, how are these stresses balanced, as is required by Newton’s laws? The simple and inescapable answer is this: at any arbitrary given distance L remote from the leading edge, the sum of the traction stresses perpendicular to the edge from x=0 up to x=L must be balanced by forces carried within the cell sheet at position L (Fig. 2e). At every position within the sheet, therefore, the accumulated traction must be balanced by local cell-borne stresses that are transmitted along the cell sheet by the cytoskeleton within cells and by cell–cell junctions between cells. Using σx x to denote stresses within the cell sheet on a plane perpendicular to the substrate and parallel to the edge, as distinct from cell tractions T at the cell–substrate interface, force balance demands that these be related by

where we take cell height, hz, to be roughly 5 μm (ref. 33), and hy is the length of the field of view.

We note, first, that if each cell were entirely self-propelled, then 〈σx x〉 would be identically zero everywhere. This possibility can now be ruled out (Fig. 2d). Instead, 〈σx x〉 increased steadily with x, reaching a stress 〈σx x300 Pa within the first 10 cell rows (Fig. 2d). Second, although they contribute little to the overall balance of forces, the pull that leading cell rows exert may be sufficient to direct the ranks immediately behind; in that sense they may still be ‘leaders’. To study further the spatial extent of this tugging mechanism of force transmission, we obtained traction maps spanning the whole diameter of the cell colony (Fig. 3). Approaching the centre of the cell colony, 〈σx x〉 increased steadily and reached values that dominated traction fluctuations. Therefore, as the cell sheet grows, it exists in a global state of tensile stress. Such tensile stress establishes definitively that a build-up of pressure caused by proliferation is not the driving mechanism that underlies expansion of the cell sheet, for if this were the case, stresses within the monolayer would be compressive and tractions would point outwards. In both regards, we found systematic evidence to the contrary. Furthermore, our data show that cell density increased proportionally to the tensile stress σx x; for reasons that remain unclear, proliferation during tissue growth seems to be regulated in such a way that the ratio between density and tensile stress is invariant (Fig. 3e). Shraiman34 proposed that lateral stresses may act as a local feedback signal to regulate the rate of tissue growth, although his model deals with monolayer compression and buckling as opposed to the tensile stresses applicable here.

Figure 3: The state of stress of the whole expanding colony is set by a global tug-of-war.
figure 3

a, Radial component of tractions along a diameter of the expanding cell colony (bar=200 μm). b, Phase contrast image. Note the presence of multicellular protrusions on the left margin of the colony. This particular shape is reminiscent of that caused by fingering instabilities in fluids9. c, Tug-of-war model. A small portion of the traction that each cell generates is transmitted to the cell behind. As such, tension in the cytoskeleton and cell–cell junctions increases towards the centre of the cell colony (red arrows). d, Once integrated over many cell rows, this small portion becomes dominant over traction fluctuations (see Fig. 2). Dotted lines indicate mean± standard error. Inset: Representative measurement of σx x over the whole diameter of a colony. σx x reaches a maximum roughly at the centre of the colony. The fact that σx x does not decrease exactly to zero at the end of the monolayer points to the existence of weak shear stresses that also contribute to the force balance. e, The ratio between radial stress and cell density across the monolayer was roughly constant. Error bars indicate standard errors.

To assess further the nature of force transmission, we considered the time evolution of tractions at specific points of the traction map. An illustrative example of such time evolution is shown in Fig. 4 (see also Supplementary Movie S2). As the first leader cell of the sheet advanced over the traction sensing point (that is, a particular pixel of the traction map), a strong pull was first observed followed by a sharp decay. This decay, however, did not drop significantly below zero as would be the case for a self-propelled cell in isolation. Instead, a second sharp pulse was observed with a maximum corresponding to the boundary between the third and fourth cells. In this particular case, the first significantly negative force was observed after the fourth row, implying that the force generated by the first three rows was almost entirely transmitted to the following rows. In every case, traction fluctuations occurred at timescales longer than the time for a whole cell to move one cell length, which further demonstrates the existence of mechanical cooperativity and force transmission at length scales larger than the single cell.

Figure 4: Time fluctuations of tractions show cell–cell force transmission.
figure 4

a, Traction Tx at a fixed point of the gel (red circle in bg) as cells crawl over that point. bg, The right-to-left advance of the cell sheet at the time points labelled in a. Note that Tx remains mostly positive throughout the observation period, which implies stress transmission. (See Supplementary Movie S2.)

In this connection, an exponential distribution of forces as reported here has been considered to be the signature of the force distributions that arise in jammed granular materials such as a pile of sand, grain in a silo or coffee beans stuck in a chute10,11,12, although the precise asymptotic form remains an open question35. The physics of jammed materials remains poorly understood, but the exponential nature of the force distribution is thought to arise from the combination of three generic features: close packing, structural disorder and long-range force transmission36. Both close packing and structural disorder are obvious properties of the cell sheet (Supplementary Movies S1 and S2), and here we provide evidence that transmission of force is long-ranged with exponentially distributed tails. Mechanics of jammed granular matter is governed by compressive stresses, however, whereas that of the cell sheet is governed by tensile stresses as demonstrated here. Although the connection between them remains a matter of speculation, the analogy between jammed inert materials and collective migration of cells is striking and suggests that these systems may share common mechanisms of long-range force transmission that are yet to be fully understood.

In summary, we present here definitive evidence establishing that collective motion in an advancing epithelial cell sheet results neither from leader cells dragging those behind, nor from cells that are individually self-propelled. Instead, each individual cell, both at the leading edge and well inside the sheet, engages in a global tug-of-war that integrates local force generation into a global state of tensile stress. Such a mechanism is innately integrative and would be inapparent in any cell studied in isolation. Whether this integrative mechanism is specific to certain tissues or instead is applicable generally during development, tissue healing and disease remains an open question, but one that is now accessible to direct experimental attack.

Methods

Cell culture.

Madin–Darby canine kidney cells (strain II) were cultured on plastic flasks in MEM with Earle’s salts supplemented with 5% FBS, 2 mM L-glutamine, 100 U ml−1 penicillin and 100 μg ml−1 streptomycin. To seed a cell colony for experiments, a 5 μl drop of supplemented media containing 5,000 cells was added to the centre of the gels. Surface tension limited the drop to approximately 1 mm2 of the central region of the gel. The cells were allowed to adhere to the gel for 30 min at 37 C and 5% CO2 before 2 ml of supplemented media was added to cover the whole surface of the dish.

Experiments.

All experiments were conducted 3–4 days after seeding the cells. Three different images were collected every 60 s, one imaging cells in phase contrast, one imaging the layer of beads located immediately below the cells and one imaging the diffraction rings of the layer of beads attached to the glass underlying the gel (Supplementary Information S2). A typical experiment lasted 6–24 h. At the end of the time course experiment, cells were trypsinized using isotonic 10× trypsin for 1 h and a stack of 30 reference images of both layers of beads was acquired. To obtain traction maps of the whole cell diameter, we merged a series of overlapping images using a correlation-based algorithm. In this case, registration was achieved by equalling displacement fields over the overlapping regions and imposing zero traction outside the colony. All experiments were conducted in the presence of serum, at 37 C and 5% CO2.

Preparation of polyacrylamide gel substrates.

Polyacrylamide substrate preparation was adapted from previously published protocols18,37 to enable image registration and improve resolution of the displacement fields (see Supplementary Information S2). The concentrations of crosslinker and polymer were adjusted for a Young’s modulus of 1.3 kPa (ref. 38). The details of the protocol for polyacrylamide gel preparation including all modifications from previous publications are provided below. Step 1: A few drops of 0.1 M NaOH were added to the centre of each 35 mm dish (glass bottom, uncoated, no.0; MatTek). The dishes were air-dried overnight. Step 2: 2–3 drops of 97% 3-aminopropyltrimethoxylsilane were added over the NaOH-stained circular regions from the previous step. The dishes were then washed and the glass surface was scrubbed with a foam swab to remove debris. The dishes were washed again and 400 μl of a solution containing 0.0001% yellow fluorescent carboxylate-modified beads (2 μm diameter, Fluospheres, Invitrogen) was added to each well. These beads were used for image registration. After the dishes dried, 0.5% glutaraldehyde in PBS was added to the central region in each dish for 30 min. The dishes were subsequently washed and air-dried overnight. Step 3: 20 μl of an acrylamide / bis-acrylamide mixture dissolved in ultrapure water containing 0.01% of 0.5-μm-diameter red fluorescent carboxylate-modified beads (Fluospheres, Invitrogen), 0.5% of ammonium persulphate and 0.05% TEMED (Bio-Rad) was added to the centre of each dish. This gel mixture was covered with glass cover slips (18 mm diameter; VWR). To ensure that all red beads laid in the top plane of the gel, the dishes were centrifuged at 500 r.p.m. for 15 min during gelation. Once polymerization was completed, the cover slips were removed. The surface was activated by adding 225 μl of a solution containing 4 μM sulphosuccinimidyl-6-(4-azido-2-nitrophenylamino) hexanoate (Sulfo-SANPAH; Pierce) dissolved in 0.1 M HEPES buffer. The dishes were then exposed to ultraviolet light for 10 min, washed twice with 0.1 M HEPES solution, washed once with PBS, coated with 1 ml of type-I collagen solution (0.1 mg ml−1; Inamed Biomaterials) and stored at 4 C. On the day before seeding the cells, the gels were washed, and incubated overnight with 3 ml of MEM with Earle’s salts supplemented with 5% FBS. Then the gel surface was allowed to dry at room temperature for 2 h immediately before seeding the cells.

Fourier-transform traction microscopy.

A new algorithm of traction microscopy was developed to account for finite substrate thickness and force imbalance within the microscope field of view (see Supplementary Information S3). Gel displacements were computed using correlation-based particle image velocimetry. To reduce systematic biases in subpixel resolution and peak-locking effects, we implemented an iterative process (n=4 iterations) based on a continuous window shift technique39. The interrogation windows were 25.6 μm on a side and window overlap ranged from 3/4 to 7/8.