| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
* Laboratoire de Physique de la Matière Condensée et Nanostructures, Université Claude Bernard Lyon 1 and CNRS, 69622 Villeurbanne Cedex, France; Department of Developmental Biology and Neurosciences, Graduate School of Life Sciences, Tohoku University, Aoba, Sendai 980-8578, Japan; and Tohoku Institute of Technology, Taihaku, 983, Sendai, Japan
Correspondence: Address reprint requests to Jean-Paul Rieu, E-mail: rieu{at}lpmcn.univ-lyon1.fr.
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAPPENDIX: AN ITERATIVE METHOD...SUPPLEMENTARY MATERIALACKNOWLEDGEMENTSREFERENCES |
|---|
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAPPENDIX: AN ITERATIVE METHOD...SUPPLEMENTARY MATERIALACKNOWLEDGEMENTSREFERENCES |
|---|
). The social amoeba Dictyostelium discoideum provides an experimentally accessible and simple model system to investigate chemotaxis and cell signaling, and to compare individual and coordinated cell movements (2). When the environment becomes depleted of nutrients, amebas stop dividing and aggregate to form a mound of cells. The mound elongates to become a finger-like structure called slug that migrates in response to light or heat. Both the aggregation of free living single amebas and the movements in multicellular aggregates including slugs are controlled by the concentration gradient of the chemoattractant cAMP produced and relayed by the cells (3–4).
Over the last decade, a large amount of information has been obtained on the distribution of cell movements along slug axis both in normal three-dimensional (3D) slugs (5–7) and in two-dimensional (2D) slugs formed at the oil water/interface (8,9). There is a characteristic pattern of movements in slugs: cells in the anterior prestalk zone show vigorous lateral movement by constantly changing their relative position, whereas cells in the central prespore zone move straight forward in the direction of slug migration in a generally periodic fashion (5–7,9). In 3D slugs, motion of prestalk cells is often helical around the central core of the tip (5,7). Siegert and Weijer (5) proposed then that the prestalk cell movement is organized by a rotating scroll wave of cAMP, which serves as pacemaker for the formation of planar cAMP waves, which in turn direct periodic forward movement of prespore cells. Periodic motions and optical waves with the same period have been observed in 3D slugs of many Dictyostelium strains (10) and in 2D slugs (9).
The distribution of mechanical forces exerted by the migrating slug was never directly measured. Hence the mechanisms by which motive force is transmitted to the substrate and their location have been subject to numerous speculations and hypotheses. According to Dormann et al. (7), anterior prestalk cells do not contribute to slug migration because the tip is often raised, only prespore cells propel the slug forward due to their close contact to the substrate. On the other hand, other studies suggested that the anterior cells exert larger forces than the posterior ones (11–13) and this is proposed to be the primary cause of the anteroposterior pattern of prestalk and prespore cells in the slug (14,15). Some models assumed that migration involves coordinated force of a special group of peripheral cells (16,17). The only indirect measurement of motive force was made by Inouye et al. (11,18,19). Assuming a viscous drag in their analysis of the experimental data, the total motive force was found proportional to slug volume suggesting that the sum of the crawling movements of the whole cells is propelling the slug forward. Many models of slug migration have since then simply postulated that the motive force is volumetric (14,15,20,21). However, the mechanisms by which interior cells can transfer forces to the substrate remain unclear. Recently, Dallon and Othmer (22) predicted that only cells in contact with the substrate can reasonably gain traction to produce a motive force for the slug. All these models take friction forces (cell/cell and cell/substrate) proportional to the relative velocity between object considered. This assumption has never been verified experimentally. The role of extracellular matrix surrounding the slug (slime sheath) made of cellulose and glycoproteins (23) was also never investigated in models.
It has been shown experimentally that the slug velocity is correlated to the slug length (14). One of the motivations of this work came from our own measurements of this relation over extended ranges of lengths, using both 2D and 3D slugs (9). To understand the origin of this relation and more generally the mechanisms of the slug migration, one must have an idea of the force distribution along slug axis which is clearly not the case presently. Therefore, we apply here for migrating slugs the flexible substrate method which was so far only used to measure mechanical forces exerted by single cells like fibroblasts or keratocytes (24). The method lies on the observation of the deformations of fluorescent beads embedded inside an elastomer substrate. Complex calculations are necessary to invert the system of coupled integral equations given by linear elasticity theory that relate the substrate deformations to the forces (25). It was demonstrated that this inverse problem is ill-posed (i.e., the solution is highly sensitive to small changes in the deformation data) for usual levels of noise and that regularization in general cannot be neglected (26). Regularization is the process of solving these problems numerically by introducing some additional information about the solution, such as an assumption on its localization, on its smoothness, or a bound on the norm.
In the Appendix of this article, we show using simulations that for slug force fields and for our range of experimental parameters (low noise level and close plane of recorded deformations), regularization is not necessary using our iterative method. Taking the advantage that slugs often have a linear steady trajectory, the deformation field is averaged in the moving slug frame. This average reduces the displacement error drastically especially for long slugs (L > 1 mm) for which the lower magnification increases bead position error. We discuss the relevance of the averaging procedure in the results section. We measure the forces exerted by migrating 3D slugs of various lengths ranging between 400 and 1100 µm. The force patterns confirm our preliminary measurements on 2D slugs (9). We find resistive forces in the slug tip and tail, traction in the central prespore area and large perpendicular forces on the sides. In addition, we find that traction and friction stresses are decreasing function of slug velocity.
|
MATERIALS AND METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAPPENDIX: AN ITERATIVE METHOD...SUPPLEMENTARY MATERIALACKNOWLEDGEMENTSREFERENCES |
|---|
). For development of three-dimensional slugs, 10 µl of a suspension at 107 cells/ml was dropped on polyacrylamide gel and incubated for 16 h before observations at 21°C.
Polyacrylamide gels
Flexible polyacrylamide gel sheets were prepared with 10% acrylamide, 0.03% bis, ammonium persulfate (10% w/v solution, 1:138 v/v), TEMED (1:1380 v/v; all products of Bio-Rad, Hercules, CA), and either 1-µm fluorescent beads (1:54 v/v) or 4-µm fluorescent beads (1:9 v/v; all beads are 2% solid; products of Molecular Probes, Eugene, OR). The mixture (450 µl) was poured on treated glass slides (25) and the droplet was flattened using a nontreated slide glass and 400-µm spacers. After polymerization (15–30 min), the nontreated slide glass was removed, the elastomer was covalently coated with type I collagen (Wako Chemicals, Osaka, Japan) using Sulfo-SANPAH (Pierce Chemical, Rockford, IL) (25). The Young's modulus of the elastomer was characterized using the method described in Wang et al. (27). We found values in the range 5–8.5 kPa. The Poisson ratio v of polyacrylamide gels was taken as 0.5 (25).
Detection of substrate deformations
We visualized simultaneously with a confocal microscope (Olympus IX70-KrAr-SPI, Tokyo, Japan) the ventral portion of the slug (transmission channel) and the beads (fluorescence channel, 488-nm line). We selected a field of view with the slug approaching on the side and we took the initial image as the undisturbed position to calculate the displacement vector of each bead. The focus plane ZM was fixed just underneath the elastomer/cell surface and was carefully measured at the end of the experiments. After thresholding and binarizing images, all bead centroids and sizes (fluorescent area) were first recorded at every time using Scion-Image (http://www.scioncorp.com, Scion, Frederick, MD). The threshold was chosen to eliminate beads <2 pixels (or sometimes larger depending on experimental conditions). We then reconstructed bead trajectories using our own C codes. Briefly, for each bead, we examine all possible corresponding beads next image (typically 1 min later) with the following criteria: i), the fluorescent area change of the bead should be <75% to eliminate beads coming into contact; ii), relative bead displacement between two successive images should be less than a value determined for each experiment (typically 2 µm, i.e., less than the larger absolute bead displacements between a given time and the initial undisturbed position because it takes several minutes to reach maximal substrate deformation at a given location); iii), if several bead pairs fulfill the previous criteria, we stop at that time the bead trajectory to avoid mixing trajectories.
The efficiency of this code depends greatly on the image resolution that is fixed by the microscope objective, the bead size, the confocal scan, and photomultiplier settings. In case of the slug A displayed in Figs. 1 and 2, which is recorded with a 20x objective and a large density of 1-µm beads (>3200 beads are in the field of view), 580 beads (i.e., 18%) could be followed until t = 45 min (Fig. 1 B). However, many beads were lost in the inner slug area where large displacements induce blinking or collisions between beads. It is possible to mark manually several hundreds of beads in the slug area using the brush tool of Scion image for interesting experiments. In case of the slug B recorded with a 40x objective, the efficiency of the tracking code is very good because the bead density is lower and the focus plane is closer to the surface (Fig. 3, A and B). More than 75% of the 440 beads could be followed during the whole experiment. Bead centroids are measured with an accuracy of about 
XY = 0.15 
0.5 µm depending on the experimental conditions. The error 
B on the displacement vector (i.e., nearly twice) was measured experimentally from nondeflected beads far from the slug.
|
|
|
= 20 µm side whose origin is fixed at slug tip. The grid is moving and crossing the recorded field of view. The mean displacement in each cell unit corresponds then to an average displacement over the set of beads NB belonging at a given time to the moving unit cell (typically NB 
30). This average over a large number of beads reduces greatly the error XY on the average displacement vector in each unit cell as 
Experimentally, we measured typically XY 0.15 µm far from the slug. There is also an uncertainty in the relative bead depth that affects the Green functions (see below). This error, 
arises from the finite depth of field of the objective (2ZF) and from the eventual geometrical roughness of the elastomer surface (ZG). With ZF 2 µm for a 20x objective and a probably overestimated roughness of the elastomer surface (at the millimeter scale) ZG 5 µm, we obtain Z 1.7 µm.
Calculations of forces
In the framework of linear elasticity theory (25), the deformation field u(r) inside a semiinfinite elastic medium caused by a distribution of forces F(r) on the surface is described by a Fredholm integral equation of the first kind:
 | (1) |
,28). If one makes the reasonable assumption that there are no vertical forces, there are four Green functions Gxx, Gxy, Gyx, and Gyy. Our method is an unconstrained method (29): we calculate the forces at the sites where deformations are measured but at finite depth ZM and we do not impose the calculated force to be zero outside the slug area. The integral equation from Eq. 1 becomes a set of linear equations, u = GF, in which u = (ux(r1), uy(r1), ux(r2), uy(r2), . . .), and F = (Fx(r1), Fy(r1), Fx(r2), Fy(r2), . . .) are 2N-vectors and G is a 4N2-matrix (typically N 2 x 103 and the matrix has >107 terms). To solve these equations, we run a numerical program based on iterative biconjugate gradient method generally used to invert large sparse matrix (30). Of course, here the matrix is dense but this method is known to give good results especially for 2D problems because, in this case, direct solution of Eq. 1 is hopeless because the matrix G is too large (30). It takes typically 5 min to invert a 107-terms matrix using a personal computer with a 1.6-GHz microprocessor.
We start from a zero initial force solution at iteration zero (#0) and the solution progressively builds up with increasing iterations. At every iteration (#), we record the force field, the mean error per site on the calculated displacements Er = |u – GF| (in micrometers or nanometers) and the mean force per unit area |Pout| and |Pin| outside and inside the slug area. Experimentally, the optimal # is chosen when Er reaches the experimental error for /ZM is large (typically 
4), or when Er is minimal within the range #1–10 for /ZM low. In the Appendix, we show using simulations that, for slug force pattern (extended force area with rather smooth variations), the minimum of øext = |Pout|/|Pin| provides a good estimation of the optimal # and of the accuracy of the calculated force pattern. Indeed, the study performed on simulated force pattern (see Appendix) shows that øext F where F = |Freal – F(#)|slug/|Freal|slug is the mean force deviation from the original real simulated pattern in the slug area. Typically it takes <10 iterations to reach the optimal #. At higher #, chaotic solutions may exist for large recorded bead depth (see Appendix), but they are easy to distinguish from regular solutions.
Deformation and stress profiles
When forces are calculated from time-averaged deformations in the moving slug frame, forces are divided by the lattice site area to obtain stress components. Parallel and perpendicular stress profiles are obtained by averaging stress over slug width W and a central band half-slug length, respectively. Parallel direction is the direction of the slug migration. When forces are calculated from raw data, parallel and perpendicular stress profiles in slug frame are calculated by first summing up forces over slug width W and a central band half-slug length, respectively, and then by dividing these sums by xW and yL/2, respectively, where x = 25 µm and y = 15 µm are the bin sizes. Error on stress profiles is determined in the Appendix using the bootstrap analysis (25,26). We have computed the parallel and perpendicular deformation components in the slug frame using the same definition. The statistical error (number of beads in each bin) and the geometrical error due to the slug angle choice are used to estimate the error on the deformation profiles. We recorded as a function of slug velocity the following characteristic features of stress profiles: the mean net traction per unit area T (often simply referred to as traction) is the average traction over the length LT of the traction area where stress is negative, including arches when present; the tip and tail frictions per unit area fP and fR are also averaged over the length LP and LR, respectively, of the friction peaks (peaks of positive stress in the tip and tail parts).
|
RESULTS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAPPENDIX: AN ITERATIVE METHOD...SUPPLEMENTARY MATERIALACKNOWLEDGEMENTSREFERENCES |
|---|
). We have recorded the deformations created by 3D slugs ranging between 400 µm and 1100 µm. These slugs induce large displacements of fluorescent beads. Fig. 1, A and B, shows the displacement fields for intermediate length slug A at two different times (see also Supplementary Material (Movie 1)). Substrate deformations present for all slugs a characteristic pattern in the slug frame. The bead displacements are oriented in the direction of migration in the tip region and opposite to this direction in the central region with a perpendicular component directed centrifugally (Fig. 1 A). In Fig. 2 E, we present the deformation field averaged over time in the slug frame using a 20-µm grid spacing. The linear portion of the trajectory between the dotted lines (see in Fig. 1 B) was used for the average. The noise level is estimated to B = 0.5 µm for nonaveraged deformations (Fig. 1, A and B) and to XY 0.15 µm for time-averaged deformations (Fig. 2 E). The maximal deformations are umax 9 and 6.6 µm for nonaveraged and time-averaged deformations, respectively. In Fig. 1, C and D, we computed the parallel and perpendicular component of these deformations. The parallel deformation presents some important fluctuations larger than the error bar, the tip friction peak is sometimes sharper (i.e., at 22 min) but the characteristic shape of the profile described by the time average is well conserved (Fig. 1 C). As a function of the perpendicular direction (Fig. 1 D), the parallel deformation is maximal in the inner area corresponding to the lighter area of the slug clearly visible in the transmission images (Fig. 1, A and B). This inner lighter area is presumably the region of close contact between cells and the substrate. For the perpendicular deformation, instantaneous and averaged profiles are all identical (Fig. 1 D). It is maximal just after the inner lighter outline. After that maximum, both the parallel and perpendicular components of the deformation field decrease continuously with perpendicular distance with no sign of break-off at slug external black outline. Far behind the slug, the beads returned slowly to their initial positions in the parallel direction but not in the perpendicular direction (see arrowheads in Fig. 1 B). The transmission image shows that the elastomer presents wrinkles in this area (see Supplementary Material (Movie 1, late times)). This surprising result indicates that the slime trail has enough mechanical stiffness to keep the substrate in a stretched state.
Typical migrating slug force pattern
The force field of the slug A is first calculated from raw data at different times when the slug is mostly entirely within the field of view (Fig. 2, A and B). The force pattern has the same characteristic properties as the force field calculated from the averaged deformation field (Fig. 2 F): large perpendicular forces, friction in the tip and in the tail, and traction in the central prespore area. It is easy to demonstrate using simulations that all these ingredients are absolutely necessary to retrieve the recorded deformation patterns (not shown). Unlike the deformation field, force (Fig. 2, A and B) or stress vectors (Fig. 2 F) are almost absent outside the slug. Fig. 2 G shows the force field recalculated after adding supplementary in-plane and vertical Gaussian noises (bootstrap method) with zero mean ± SD XY 0.15 and Z 1.7 µm (values discussed in the Materials and Methods section). The similarity between Fig. 2, F and G, demonstrates the robustness of the solution.
Largest forces are exerted at the outline of the lighter close contact area perpendicular to slug axis (Fig. 2 D) where the bead deformations are larger (Fig. 1 D). These perpendicular forces decrease almost to zero at the external black outline of the slug. The parallel stress profiles (Fig. 2 C) present large fluctuations with time but at every time we found two positive peaks in the anterior (tip) and very posterior region of the slug indicating resistive forces there. Although smaller than in the tip, parallel and perpendicular stresses in the posterior region extend beyond the tail of the slug in the trail left by the collapsed sheath (see arrowheads in Fig. 2 F) and are responsible for the deformations already commented in this area (Fig. 1 B). In average, the remaining parallel component in the trail is always resistive (see Fig. 2 F and also 3D).
In the central prespore zone of the slug, parallel stress is negative in average indicating dominance of traction here (Fig. 2 C). Note also the presence of two distinctive traction zones at 75 and 600 µm from the tip separated by a central zone with almost vanishing parallel forces. This traction free zone is often but not always present in other investigated slugs. It corresponds probably to the presence of a steady arch as found previously (9). The overall sum of parallel or perpendicular stresses is not significantly different from zero as expected for a slowly moving body (25). However, to satisfy exactly the force balance it would be necessary to sum up the forces along the full slug trajectory due to the remaining friction forces in the trail.
Fluctuations, correlations, and waves of forces
When we averaged deformations over time in the moving slug frame, we selected portions of their trajectory where slugs are moving straight at constant velocity with a constant shape. However, before presenting the slug velocity dependence of time-averaged forces, we must question the statistical validity of our averaging procedure.
Only two times are represented in Fig. 2, A–D, but we calculated at 12 different times the force field from raw displacement data with a slug almost entirely in the field of view. We found always a perfect correspondence between instantaneous and time-averaged perpendicular forces (Fig. 2 D). Instantaneous parallel force profiles although presenting fluctuations follow also the characteristic averaged profile described above (Fig. 2 C). In particular, we always find the two distinctive traction zones separated by a central traction free zone that is also present in the force calculated from the time-averaged displacements. When averaging the force profiles obtained at different times, central fluctuations are smoothed and we obtain almost exactly the force calculated from the time-averaged displacements. For our purpose of measuring the mean traction and friction forces (averaged over traction and friction areas, respectively) as a function of slug velocity, we can state that the time-averaged deformations and subsequent force calculations are well representative of slug behavior. They are also robust as we used different time intervals and different portions of the slug trajectory, and we never observed significant difference between the resulting force or deformation patterns.
We may now ask whether the large fluctuations present in the parallel deformations and stresses are correlated (Figs. 1 C and 2 C). Do they correspond to waves of forces? We know that optical waves and periodic motions with the same period have been observed in 3D slugs of many Dictyostelium strains (10) and in 2D slugs (9). Although a direct chemical demonstration is still required, it was proposed that these waves of motion correspond to periodic cAMP waves originating in the prestalk tip and progressing backward toward the tail (5). Whatever their exact origin, waves may just polarize the cell deformations in the migrating direction or increase the motive force of the cells or both.
Obviously, to detect waves of forces, it is necessary first that some periodicity could be detected in the deformation raw data. We performed for that purpose an autocorrelation-function calculation for the parallel component of deformation. The results indicate that a weak oscillatory component of a 6-min period may exist (Fig. 1 E). However, by identifying the amplitude of the oscillation with the calculated amplitude of the autocorrelation function of u = udc + uac cos(
t), we found that the experimental uac is as small as 10% of udc. This estimation is comparable with the estimated error bars: in the central traction area of Fig. 1 C, mean deformation is 
u
= 2.7 µm and error that originates from experimental resolution is u = 0.24 µm. Fluctuations larger than 10% from the time-averaged deformation exist in the parallel direction (Fig. 1 C) but they are not correlated temporally.
Assuming, however, that the 10% uac oscillation is meaningful, can we calculate the oscillatory part fac of the force from which it originates? For that, we used a realistic force pattern keeping the same features with the experimental pattern of Fig. 2 (same tip and tail friction; same perpendicular forces) but symmetric with respect to slug to avoid any effect from the perpendicular direction. We introduced a parallel force modulation f = fdc + fac cos(2
x/
+ 
) in the central traction area, keeping the total traction force constant, equal to our experimental result from time-averaged data (see supplemental Fig. 9, Supplementary Material). We tested several ratio p = fac / fdc and found that udc/uac = 10% corresponds roughly to p = 25%. This estimation is comparable to the estimated error on force F from time-averaged data (case /ZM = 2; see Appendix) but lower than the one for nonaveraged data. From the series of the 12 different times we calculated forces for the slug A, we never observed any indication of correlated or propagating forces. Therefore, at that moment, there are no clues for the existence of oscillatory forces due to chemical waves. However, if they exist they are certainly smaller than the uncorrelated forces. The irregular tip lift-landing events clearly visible in the supplemental Movie 1 (see Supplementary Material) seem responsible for these uncorrelated fluctuations.
Measurements of the stress field as a function of the slug velocity
The stress fields for a long fast slug (slug B, L = 1085 µm) and a small slow slug (slug C, L = 592 µm) are displayed in Fig. 3. The deformation fields (at a given time for slug B in Fig. 3 A; averaged over time in Fig. 3, C and E) have been recorded at a similar mean bead depth (ZM = 3 and 5 µm, respectively) with the same grid spacing = 20 µm. The force profiles are qualitatively similar to those of the slug of intermediate length (Fig. 2, C and D). Friction in the tip is generally larger than in the tail. Perpendicular stresses are larger than parallel ones. But, surprisingly, both the parallel and the perpendicular forces are larger for the smaller slug (Fig. 3, G and H).
We have reported in Fig. 4 the characteristic features of the parallel stress profiles as a function of the slug velocity for nine investigated slugs. Both the absolute value of the negative traction T and frictions are decreasing functions of slug velocity (Fig. 4, A and B). The reported friction is the mean overall friction on both tip and tail defined as F = (LP x fP + LR x fR)/(LP + LR), where LP and LR are the lengths of the tip and the tail + trail friction areas. Traction is best fitted by a power law relation V–1.36 (solid line in Fig. 4 A) but, due to the large experimental error bars, an inverse velocity form V–1 (dotted line) that is used in the discussion to model the length dependence of the slug velocity also fits the data. Tip and tail frictions may be fitted by the power law decreasing functions V–1.31 and V–0.96 (not shown), respectively. The overall friction scales as V–1.26 (see solid line in Fig. 4 B).
|
20 and 12%, respectively, on the relative error on mean parallel stress components: mean traction (Fig. 4 A) and tip and tail frictions. The overall friction F (Fig. 4 B) has a supplementary error due to the determination of LP and LR. The amplitude of perpendicular stress, calculated from the maximum of the vertical stress profile, is also decreasing function of slug velocity (not shown). But due to the difficulty of finding sharp forces when noise increases (as discussed in the second section of the Appendix), perpendicular forces suffer large errors. We also measured the lengths of the traction area (LT), of the tip (LP) and of the tail + trail (LR) friction areas (Fig. 5). LT is clearly linearly increasing with slug length L (not with slug velocity V). LP and LR values are close to the slug width W (thick gray line). But the large scatter in the data for LP, LR, W, and the overall friction length LF = LP + LR prevents conclusion about the dependence of these quantities with L (or with V).
|
|
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAPPENDIX: AN ITERATIVE METHOD...SUPPLEMENTARY MATERIALACKNOWLEDGEMENTSREFERENCES |
|---|
) is confirmed by this work for all 3D slugs investigated (Fig. 3 A). Symmetric lateral movements from the slug axis toward periphery have never been reported in 3D or 2D slugs. In 2D slugs, cells presented sometimes alternatively elongated shapes but mostly in the slug migration direction (9). Therefore, the observed perpendicular forces do not result from the lateral movements of the cells; they are likely the result from the mechanical forces transmitted by the slime sheath to the substrate. New surface sheath made of cellulose and glycoproteins is continuously formed in the tip (31). It is smooth and extremely thin, only 10–50 nm as measured by electron microscopy, but acts as a protective shell for the slug. Once it makes contact with the substrate, it attaches strongly to it and is left behind as a collapsed tube (slug trail). It is known that newly formed sheath in the tip part of the slug can be deformed by anterior cells but in the posterior part cells move in is a rigid and immobile tube (32). Our current view of sheath/substrate force transmission mechanism is described in Fig. 6. As the slug advances, sheath production and sheath/substrate anchorage in the tip create a resistive force fP pushing the elastomer forward in the direction of slug migration and possibly also with a component in a direction normal to the hemispherical tip (Fig. 6, A and B). Although sheath is not moving, it is probably stretched and pulled forward by the moving cell mass. This forward stretching force may cause the parallel friction fR in the tail (Fig. 6, A and B). We can postulate that the motive force of the slug is supplied by cell traction T mostly in the central prespore region in the inner lighter contact area between slug and substrate, and that this traction is essentially independent from perpendicular forces P and from the sheath. This hypothesis is consistent with the calculated forces (Figs. 2 and 3) and with the measured localization of maximum parallel and perpendicular bead deformations (Fig. 1 D). We will examine below in more details the force balance equation between parallel traction and friction forces. Before that, we are now presenting possible mechanisms responsible for the large perpendicular forces.
|
). We assumed that due to its tension, the sheath can pull out the elastic substrate in the direction tangential to the slug profile where it is anchored in the vicinity of the contact line (Fig. 6 C). There are several possible origins of the sheath tension. One possibility is that frustrated cell movements in the lateral direction cause an internal pressure on the sheath. Another is that slug/substrate adhesion induces sheath tension in a mechanism analogous to the lipid membrane tension induced by the adhesion of vesicles to surfaces (33). As the elastic substrate assays do not visualize vertical bead displacements, we visualize only the in-plane component of the stretching forces. In case of slugs, there is evidence that the slug has a flattened cylindrical shape from the visualization of the slug ventral profile using the confocal microscope (see also the relative size of the lighter inner and dark areas of the slug in Fig. 1). If slug/substrate contact angle 
(Fig. 6 C) is much larger than 90°, an in-plane stretching force directed outward exists. The localization of deformations and perpendicular forces in the vicinity of the outline of the inner lighter contact area (Figs. 1 D and 2 D) supports this hypothesis. Note that a broad distribution of attachment points between sheath and substrate in this area may spread the perpendicular forces. However, this mechanism does not explain the remaining perpendicular deformations and forces often observed far from the slug in the collapsed sheath trail (arrowheads in Figs. 1 B and 2 F). Here, it is possible the substrate is kept stretched by the frozen sheath (32).
Forces are decreasing with slug velocity
The most important result of this article is that traction and friction forces per unit area are decreasing with the slug length or the slug velocity.
We can draw several important consequences from these measurements. First, traction is not constant with slug length or slug velocity (Fig. 4 A) contrary to what is usually assumed in models (14,15,20–22). As a result, the total motive force TLTW is roughly constant at –5 µN (not shown). Second, friction is not a viscous fluid friction increasing with slug velocity as always assumed in these models, but it is probably close to solid friction, stick-slip, or other mechanisms (see below). Moreover, this finding has a very important consequence on the interpretation of the experimental measurements of motive force of Inouye et al. (11,18,19). They indeed assumed a fluid-type resistive force to extrapolate the force. As we show that friction is decreasing with the slug velocity (Fig. 4 B), we believe that the measured motive force is erroneous and the conclusion that motive force is volumetric is not founded. We cannot conclude, however, from our measurements whether traction is a volume or a surface force; it will be necessary for that to measure carefully forces as a function of height for a given slug length.
A velocity dependence of cell traction forces was never directly reported to our knowledge, however, several independent observations suggest that for single cells the same relation probably holds. During ameboid migration, traction forces are generated by the actin-myosin cytoskeleton (24) and coordinated with other events: extension of protrusions at the leading edge, attachment/detachment to the substrate, translocation of the cell body (34). In slowly moving fibroblasts (V 0.5 µm/min), the force pattern has been resolved to a resolution of 2 µm (25). Traction at the leading edge is balanced by friction in the tail caused by the passive stretching of stationary adhesive linkages as the cell moves forward. Mean traction stress is 2 kPa. In fast-moving keratocytes (V 30 µm/min), significantly smaller traction forces are detected (i.e., 0.2 kPa) (35). Moreover, it was found that fibroblast cell speed and traction force depend on the rigidity of the substrate. Fibroblasts move slower but exert a significantly higher traction force on stiffer polyacrylamide substrates (36).
Dictyostelium single cells are also fastly moving but forces are more difficult to resolve due to the small cell size and the fast cell shape changes. It was recently observed by Uchida and Yumura that cell velocity is inversely proportional to the number of active cell/substrate adhesion sites referred to as actin foci (37). They suggested that actin foci constitute the active feet of Dictyostelium cells. The parallel between these observations and ours is striking. The fact that friction and traction forces are decreasing with slug velocity (Fig. 4, A and B) suggest that the dynamics of the cell/substrate adhesion sites may regulate the slug velocity. Too many sites create large traction but also large friction; as a result the cell is slowly moving like fibroblasts. On the other hand, once the cell is moving fast, the probability for an actin bundle to make an adhesion site with the surface is inversely proportional to the velocity.
In dry solid frictional sliding between a rough surface and a smooth surface, one observes often a stick-slip regime in the range of low velocities. In this regime, the friction coefficient is described by a logarithmic decreasing function of V due to a logarithmically strengthening of adhesive joints when aging at rest (stick) and their more frequent rejuvenation when slipping fast (38). At higher velocity (V > 60 µm/min), the friction coefficient increases again logarithmically. The fact that we are measuring a friction decreasing with slug velocity in the range 10 30 µm/min suggests that cell/substrate and/or sheath/substrate friction is compatible with a stick-slip mechanism.
Both in the study of actin foci or of solid friction cited here, the velocity of the moving or sliding object is playing a central role, but not the macroscopic size of the object. Substrate/moving object interactions indeed always occur at the scale of the micrometric size asperities or adhesion sites. This justifies fully our assumption for the scaling relations derived below that forces depend on slug velocity but not on slug length.
Scaling relations of slug migration
One of the important questions related to slug migration is why the slug velocity V correlates with the slug length L (14). We have reported in Fig. 7 our measurements of 2D and 3D slug velocities on acrylamyde (9) as well as 3D results on agar (14). Several numerical models addressed this point. They found that slug speed increases with slug length L but saturates when L is larger than the wavelength of chemoattractant waves (21,39). However experimentally, slug speed still increases by a factor of 8 10 between the wavelength, i.e., 250 µm (9), and the larger 3D slugs observed (Fig. 7).
|
 | (2) |
). In dry solid friction, friction decreases and then increases logarithmically with V (38). Here, we found that the simplest form that captures the saturation of the slug velocity is a decreasing function of V, which saturates at high velocities. The relation F = f0 + f1/V = 60 + 1036/V is therefore used (f1 is a free parameter of the fit and f0 was chosen to retrieve the experimental saturation of the velocity).
We postulate now that the lengths of the traction and friction area are proportional to L (not to V). Experimentally, LT shows an almost perfect linear increase with L (Fig. 5, LT = BTL = 0.75L, in micrometers) not with V (not shown). LF is probably also governed by slug geometry. Tip friction occurs where the newly produced sheath attaches to the substrate, in the contact area between the flattened hemispherical tip and substrate. Tail and trail frictions are caused by collapsed sheath resistance that should occur on a length roughly proportional to the slug lateral dimension. Experimentally, LP, LR, and slug width W are clearly of the same order but it is hard to find a correlation with L due to the large scatter of the data (Fig. 5). Extensive studies of slug morphology for a large population of NC-4 slugs showed that W is increasing with L (40). The fact that W is not clearly increasing with L for the selected slugs of this study is not surprising with respect to the small number of experiments analyzed and the large scatter inherent to such measurements (40). Finally, we use the following relation for the friction length LF = AF + BFL = 120 + 0.25L (in micrometers) where AF was fitted (see Fig. 5) and BF was chosen to satisfy a sum LT + LF –L = 120 µm independent of L. This value represents the characteristic friction length remaining outside the slug in the trail of the slime sheath.
By injecting the experimental fits of LF(L), LT(L), F(V), and T(V) into Eq. 2, we obtain the length dependence of the velocity (solid line in Fig. 7). The agreement with the experimental velocity data for the full range of slug lengths is reasonably good except for very small slugs. The reason why the velocity saturates for large slugs (L > 3 mm) is mainly due to the dissymmetry at large V between slug friction that is constant (F
= f0) and traction that is decreasing to zero. The maximum velocity (V = [BT x t1 – BF x f1]/BFf0 50 µm/min) is positive because the traction force decreases faster than the friction forces when the slug size increases (i.e., BT x t1 > BF x f1). Of course, the choice of the scaling relations is not unique. In particular, logarithmic or power-law decreasing functions of the stresses may satisfy as well the force balance equation (Eq. 2). The important point to be kept in mind is that forces decreasing with slug velocity satisfy both the force balance equation and the relation between slug length and slug velocity.
In conclusion, we directly measured for the first time the distribution of mechanical forces of migrating slugs of different sizes. We found that traction forces as well as friction forces decrease with the slug velocity. In addition, we found a large force perpendicular to slug migration direction. Our measurements and scaling relations propose a new explanation of the well-known relation between slug velocity and slug length. We suggest that it will be interesting to integrate these new findings in numerical models of slug migration taking explicitly forces into account (15,20–22).
|
APPENDIX: AN ITERATIVE METHOD TO CALCULATE FORCES FROM ACCURATELY COLLECTED DEFORMATIONS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAPPENDIX: AN ITERATIVE METHOD...SUPPLEMENTARY MATERIALACKNOWLEDGEMENTSREFERENCES |
|---|
) to solve the coupled equations of elasticity (Eq. 1). In this Appendix, we address the details of our computational method. We show first how the iteration and the solution are chosen and we then study the influence on the solutions of the experimental parameters such as the noise on bead deformations XY (simply referred to as "noise") or the mean recorded bead depth ZM. To find the validity range of our computational method, i.e., the range of the experimental parameters for which our method gives robust and accurate solution, we simulated two very different force patterns.
The first artificial force pattern mimics the experimental force pattern of migrating Dictyostelium slug (supplemental Fig. 9 A, Supplementary Material). In the traction area, we introduced a parallel force modulation f = fdc + fac cos(2x/ + ) (fac/fdc = 0.5), keeping the total traction force constant and balanced by friction forces. One test of our force reconstruction will be whether these "waves of forces" can be retrieved when changing experimental parameters. We took a Young modulus of 5 kPa giving a maximal deformation umax = 8.6 µm at ZM = 5 µm.
The second simulated force pattern mimics forces exerted by a fish keratocyte with sharp high traction forces directed inward at the edge of the cell and low forces in the middle separated by 3 µm (supplemental Fig. 10, Supplementary Material). This second pattern allows an exploration of the applicability of our method to measure forces exerted by single cells. We also took here a Young modulus of 5 kPa giving a maximal deformation umax = 3 µm at ZM = 0.75 µm with the typical forces measured in Doyle et al. (35).
Starting from the simulated force pattern, we use Eq. 1 to calculate the displacement field at different ZM. We add Gaussian noise of zero mean ± SD XY to the displacement field and then we run our numerical program to reconstruct the forces. This process is often termed as a "Monte Carlo" simulation. In simulations, contrary to experiments, the force deviation F is exactly known.
Choice of the optimal iteration
For all simulated force patterns and for all slug experiments, we always found two classes of behavior of the solutions as a function of the iteration number (#). These two classes depend mostly on the mean plane of recorded bead deformations ZM.
The first class is represented in Fig. 8 A and corresponds to the case of low ZM. Actually the important parameter is the ratio between the spatial resolution (the distance between two neighboring sites of deformation) and ZM, which should be larger than or equal to /ZM 4 (the precise estimate depends on the type of the force pattern). In that case, the displacement error Er decreases monotonically to zero; it crosses the experimental level of noise XY approximately when the force deviation F presents its minimum, typically between #5 and #10 depending on the level of noise. The fact that Er is continuously decreasing toward zero with # indicates that the code is able to find a solution that accommodates any added noise. Of course, such a solution is not physical and we should not exceed the iteration at which Er = XY. Er = XY defines then the criterion to choose the optimal # in experiments when /ZM is large. Here, it gives the #8, which is slightly different than #6 where F is minimal. However, solutions at iterations 5–8 are really difficult to distinguish from each other (supplemental Fig. 11, A and B, Supplementary Material).
|
F is. With XY = 0.15 µm (umax/XY = 60), /ZM = 4 and in the case of the slug pattern, the optimal solution at #6 (as given by the minimum of the force deviation F) has a F = 15%. This solution (Fig. 8 C) is very close to the original one (supplemental Fig. 9 A, Supplementary Material). It is robust when changing or increasing the noise. The stress noise outside the slug area is low and the slug boundary is sharp. At low # (i.e., #1–3), the stress pattern is still not fully reconstructed and resembles the deformation pattern in the sense that there is no sharp vanishing of the stress vectors after the slug boundary (see the stress profiles in supplemental Fig. 11, A and B, Supplementary Material). The mean stress inside |Pin| the slug area increases up to #6–7 and then remains constant. On the other hand, F increases until #13–14 because directional noise on stress vector increases although |Pin| is nearly constant (Fig. 8 A). The remaining stress outside the slug area |Pout| also increases until #13–14. As a result the proportion of external stress øext shows a minimum at 5 iterations, increases, and finally stabilizes at large iterations (Fig. 8 A). Although the minima are not exactly at the same #, the similarity between the curves of øext and F is striking and is observed for all slug-type force patterns and all ZM. It indicates that the minimum of øext provides a good estimate of the force deviation at the optimal #. Note finally that after #15 and until the last iteration we could investigate (#190 where Er = 10–27 µm !), the code converged to a stable solution as |Pin|, |Pout| and the calculated pattern no more change.
The error landscape is very different at same noise for a larger ZM = 10 µm (/ZM = 2; Fig. 8 B). Er is a nonmonotonic function of the # and presents a series of minima separated by large peaks (note the log scale in Fig. 9 B). The first minimum of Er at #6 corresponds to the lower F in the range #0–30. At this optimal # given by F, the reconstructed stress pattern (supplemental Fig. 8 D, Supplementary Material) is reasonably close to the original one (supplemental Fig. 9 A, Supplementary Material). The force deviation F = 25.4% is larger than the one found for ZM = 5 µm and Er does not reach the added noise XY = 0.15 µm. It is possible to find large # for which the absolute minimum of Er is smaller. But after #12, the solutions are in fact completely erratic (Fig. 8 E). The criterion to find the optimal # in experiments when /ZM is low is therefore to take the minimum of Er within the range 1–10 iterations. We indeed never found optimal solutions at higher iterations for a level of noise consistent with experiments. Again, here for large ZM, øext F and its minimum provides good estimate of the experimental accuracy without any prior knowledge of F.
To show that these results do not depend on a given type of pattern, we reconstruct now a very different keratocyte-type force pattern (supplemental Fig. 10 A, Supplementary Material). The error landscape (supplemental Fig. 10 B, Supplementary Material) is qualitatively similar than the one of the slug at similar geometrical and signal/noise ratios (/ZM = 4 and umax/XY = 60). Although less marked, a minimum F = 13.4% at #11 exists, and Er decreases monotonically well below the noise level. The criterion Er = XY gives #7, but F is only 15.2% at this # and the reconstructed force patterns are not distinguishable within this # range #6–. Both |Pin| and |Pout| are indeed perfectly constant after #12. The real force vectors of smaller amplitude within the cell area can be clearly distinguished from external force noise at a lower noise XY. At /ZM = 2, even for much lower noise, the sharp traction forces at the edges of the keratocyte are dramatically spread over the two or three next neighboring sites. The required range of parameters of our method for single-cell experiments seems presently difficult to fulfill experimentally, contrary to the slug case. Details and a more complete discussion will be given in a forthcoming article.
To summarize this first Appendix section, our iterative algorithm is numerically stable in a large range of parameters. For signal/noise ratio umax/XY > 20, and a geometrical ratio /ZM 2 for slugs and 4 for keratocytes, it gives always an optimal solution characterized by a minimum of the force deviation F from the original simulated force in the range 10–30%. This solution is robust when adding or changing noise contrary to that we could expect for erratic and nonregularized solutions. Therefore, we confirm here what was already implicitly demonstrated that regularization is less relevant for low noise level and lower spatial resolution (26).
Accuracy of slug force measurements
We already briefly discussed that there is little quantitative difference in terms of F between the optimal and the neighboring iterations. However, the parallel profile shows little difference between all # investigated in supplemental Fig. 11 A (Supplementary Material) (except #3 in the tip friction peak). A small difference on the perpendicular stress profiles between #5 and optimal #6 is on the other hand detectable and the difference is huge with #3 (supplemental Fig. 11 B, Supplementary Material). Note also that the perpendicular profile of #6 itself is slightly more spread than the real simulated pattern at the slug external boundary. This means that in the regions of the field where
).
), tail friction area (
), and sum of tip + tail friction area (
