{'fraction_non_alphanumeric': 0.028595531766348617, 'fraction_numerical': 0.015621363742145683, 'mean_word_length': 4.717112922002328, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
An Oscillatory Contractile Pole-Force Component Dominates the Traction Forces Exerted by Migrating Amoeboid Cells

We used principal component analysis to dissect the mechanics of chemotaxis of amoeboid cells into a reduced set of dominant components of cellular traction forces and shape changes. The dominant traction force component in wild-type cells accounted for ~40% of the mechanical work performed by these cells, and consisted of the cell attaching at front and back contracting the substrate towards its centroid (pole-force). The time evolution of this pole-force component was responsible for the periodic variations of cell length and strain energy that the cells underwent during migration. We identified four additional canonical components, reproducible from cell to cell, overall accounting for an additional ~20% of mechanical work, and associated with events such as lateral protrusion of pseudopodia. We analyzed mutant strains with contractility defects to quantify the role that non-muscle Myosin II (MyoII) plays in amoeboid motility. In MyoII essential light chain null cells the polar-force component remained dominant. On the other hand, MyoII heavy chain null cells exhibited a different dominant traction force component, with a marked increase in lateral contractile forces, suggesting that cortical contractility and/or enhanced lateral adhesions are important for motility in this cell line. By compressing the mechanics of chemotaxing cells into a reduced set of temporally-resolved degrees of freedom, the present study may contribute to refined models of cell migration that incorporate cell-substrate interactions. Electronic supplementary material The online version of this article (doi:10.1007/s12195-011-0184-9) contains supplementary material, which is available to authorized users.


INTRODUCTION
Amoeboid cell migration is central to many physiological processes, including wound healing, immune response, or embryonic development. 20,22,29,30 Its deregulation is a very important factor for serious pathological processes, as varied as tumor invasion, vascular diseases, or mental retardation. 9,24,37 Cell migration requires the coordinated regulation of many molecules, but occurs through a small repertoire of shape changes consisting mainly of a sequence of anterior protrusions and posterior retractions. 2 There is ample evidence indicating that these morphological shape changes are mechanically coupled to the biochemical modulation of the adhesions via the traction stresses exerted by the cell on the extracellular matrix (ECM). 6,13,26,34,39 However, there is a need to better understand the coupling of the global cell shape changes, the distribution of traction stresses and the associated biochemical processes.
Cell locomotion has been extensively studied using cell shape as readout. Wessels and Soll 36,40 developed a framework for the 2D and 3D analysis of shape changes during cell locomotion, which enabled them to track the motion of amoeboid cells 20 through the observation of the time evolution of descriptors such as the cell length or the velocity of the cell centroid. Machacek and Danuser 25 identified different protrusion phenotypes using level set methods. Killich et al. 19 and Keren et al. 18 used Principal Component Analysis (PCA) (also known as Karhunen-Loe`ve decomposition or Proper Orthogonal Decomposition) to quantify cell shape changes in migrating Dictyostelium cells and epithelial fish keratocytes. These studies employed a version of PCA that optimized the variance of the measured contour of the cell, to obtain a reduced set of dominant components of cell shape changes during migration. However, there is still the need to understand the dominant traction force patterns, and the associated biochemical processes, that the cell employs to produce mechanical work (energy) during migration.
Measuring the mechanical energy spent by the cell during migration requires having experimental access to the traction forces imposed by the cells on their ECM. 4 Harris et al. 14 identified qualitatively the nature of the stresses exerted by cells on the ECM through the observation of the wrinkling patterns cells induced over a thin layer of silicone rubber substratum. Dembo et al. 11 and Peterson 28 first described the use of elastic deformable substrates to measure the traction stresses exerted by cells attached to the substrate surface. This technique has since been applied to study the traction stresses exerted by different types of cells, including locomoting keratocytes, 11 fibroblasts 28 and Dictyostelium. 10,23,26 The objective of this work is to refine the spatiotemporal description of motility by obtaining a statistically significant characterization of the relationship between traction forces and cell shape in migrating amoeboid cells, using PCA as primary tool. For this purpose, we carried out simultaneous measurements of the time evolution of cell shape and traction forces exerted by isolated Dictyostelium cells undergoing chemotactic migration over a flat elastic substrate. In order to provide a physically relevant description, we applied a modified version of PCA to the measurements of the traction forces that optimized the strain energy accounted by each component. This analysis allowed us to determine the principal traction force patterns and shape changes implemented by migrating cells according to their importance for strain energy production. These patterns have a clear physical interpretation because the strain energy represents the amount of mechanical work performed by the cell against its surroundings while it is moving.
PCA effectively ''compressed'' the mechanics of migrating Dictyostelium cells by identifying a minimal set of traction force components that accounted the work performed by the cells. We found that the most important of these components accounted for~40% of the traction work, and consists of a simultaneous contraction of the substrate from front and back towards the center of the cell. Its temporal evolution was oscillatory and correlated well with that of the strain energy and the cell length. In addition, we found that a set of five canonical principal components, exhibiting a high degree of repeatability from cell to cell, accounted for virtually all of the traction work exerted by WT cells. We also identified the most relevant principal components in two mutant strains with contractility defects: MyoII Essential Light Chain Null (mlcE À ) and MyoII Heavy Chain Null (mhcA À ) cells and compared it to the dominant components in WT.
The identification of a single most relevant component that accounts for a large portion of the strain energy opens new venues to develop semiempirical modeling approaches 3,27,32,33 to describe the motion of the cell.

Dictyostelium Culture and Microscopy
Axenically grown Dictyostelium wild-type (WT) and mutant cells were prepared for chemotaxis and seeded onto a flat elastic gelatin gel as described previously. 10,26 Time-lapse sequences of chemotaxing cells were acquired on an inverted microscope controlled by Metamorph software (Molecular Devices, Downingtown, PA).

Cell Shape Segmentation
Cell outlines where determined from differential interference contrast microscopy images using standard segmentation techniques. 10

Traction Cytometry Experiments
The traction stresses exerted by cells were determined from measurements of the substrate deformation after solving the equation of static equilibrium for a linearly-elastic substrate in Fourier space, as described by del Á lamo et al. 10 (see Fig. 1a). The substrate deformation was obtained from the lateral displacements u ¼ ðu; vÞ of 0.1-lm fluorescent latex beads embedded in the gel. These displacements were measured by cross-correlating each instantaneous image with a reference image of the relaxed substrate. 41 The Young's modulus of the gel was determined from the indentation of a tungsten carbide sphere. 17 In the calculations, the tensile stresses at the free surface of the substrate were assumed to be negligible (s zz ðz ¼ hÞ $ 0), so that the stresses at the z ¼ h could be characterized by the 2D vector s ¼ ðs zx ; s zy Þj z¼h :

Cholesky Decomposition of the Strain Energy U S Exerted on a Flat Substrate
The strain energy exerted by the cell on the surface of the substrate at a given instant of time was 4 where ðÁÞ represents scalar product, R A dx represents a surface integral in a box enclosing the cell, and fg represents an inner product. Note that U S is equivalent to the mechanical work the cell is exerts on the substrate.
We applied PCA to obtain a reduced set of traction force patterns which the cell used to exert mechanical work on the substrate. PCA works by optimizing the modal decomposition of a set of observations according to the norm associated to a predefined inner product. We then rewrote the strain energy as the norm U S ¼ fh; hg; where hðxÞ is the strain energy square root or U S -square root. This decomposition is called the Cholesky decomposition and is easier to perform in the Fourier domain, where b u ab are the Fourier coefficients of the lateral displacement vector, A a;b is the Fourier norm matrix and ðÞ Ã denotes Hermitian transposition. In Fourier space, the problem of finding the functional Cholesky decomposition of the norm operator is reduced to finding the algebraic Cholesky decomposition of the norm matrix, A ab ¼ Q Ã ab Q ab ; where Q ab is a lowertriangular matrix. Introducing this decomposition in Eq. (2), one obtains that b h ab ¼ Q Ã ab b u ab ; and inversion of the Fourier transform yields Cell-based Coordinate System Prior to applying PCA, we mapped the U S -square root hðxÞ in a non-dimensional cell-based coordinate system that took into account the changes in cell shape and orientation (see Fig. 1b and del Alamo and co-workers 10,26 ). The origin of this coordinate system was the cell centroid, ðx c ðtÞ; y c ðtÞÞ; the coordinates were rescaled with the half-length of the cell, LðtÞ=2; and the axes were parallel to the principal axes of the cell. The coordinates in this cell-based reference system were where ðx; yÞ were the coordinates in the laboratory reference frame, and /ðtÞ was the instantaneous angle between the longitudinal axis of the cell and the x-axis of the laboratory reference frame (see Fig. 1b). In order to ensure that the magnitude of both the traction forces and the strain energy were conserved upon transformation of reference frame, the traction stresses and the U S -square root were redefined respectively as s ? ðx ? Þ ¼ L 2 sðx ? Þ=4 and h ? ðx ? Þ ¼ Lhðx ? Þ=2: For simplicity, hereafter we will refer to x ? as x and to h ? ðx ? Þ as hðxÞ: The size of the computational box where traction forces were calculated was 6 non-dimensional units, and was discretized using 129 Â 129 nodes.

Scalar Field Mapping of Cell Shape Allows to Relate it to Traction Forces
The instantaneous cell shape is mapped using the non-dimensional cell-based coordinate system defined in Eq. (4) as a 2D function dðxÞ such that dðxÞ ¼ 1 if x laid within the contour of the cell, and dðxÞ ¼ 0 otherwise (Fig. 2a). This function was discontinuous at the contour of the cell, which could lead to undesired spurious oscillations in its modal decomposition (the Gibbs phenomenom, see Canuto et al. 5 ), and deteriorates the convergence of PCA. This was avoided by using a continuous scalar field sðxÞ to represent cell shape where cðxÞ was the distance from each point inside the cell x to the cell contour (see Fig. 2b). Because our aim was to determine how the cell shape associates with strain energy and traction forces at each instant of time, we applied PCA on the composite three dimensional function which grouped both the U S square root and the shape of the cell. The value of parameter e was set very low (e ¼ 10 À9 À 10 À12 ) to ensure that the results from PCA only depended on the strain energy and were independent on shape variance. We verified that the use of the small parameter ; as well as its selected magnitude, did not affect the calculation of the components of hðxÞ: Principal Component Analysis (PCA) on the Strain Energy Square Root and Associated Cell Shape Given a set of N experimental recordings of cell shape and traction forces (i.e. at times t ¼ t k for k ¼ 1 . . . N), we applied PCA 1,15 to the composite function wðx; tÞ (Eq. (6)), which allowed us to express it as the weighted sum of principal components where the basis functions w k ðxÞ ¼ /ðxÞ; evðxÞ ð Þ k ; contain the spatial structure of cell shape and tracion forces of each principal component and are denoted principal functions, and a k is the weight coefficient of each component. To facilitate the interpretation of the principal functions, we transformed them into their traction force equivalents, where hi denotes temporal average and hLi is the average cell length. The principal functions are mutually orthonormal with respect to the inner product (i.e. fw i ; w j g ¼ 1 if i = j and zero otherwise), so that the total strain energy is given by This property allows us to evaluate the instantaneous contribution of each term of the sum to the strain energy, U k S ðtÞ ¼ ja k ðtÞj 2 : The time-averaged contribution of each component is usually referred to as the principal value associated to the component, k k : The principal components are arranged in decreasing order of k k ; (i.e. k kþ1 <k k ). The relative contribution of each principal component to the strain energy is obtained from the ratio m k ¼ k k =hU S i: The defining property of PCA is that, when applied to the U S square-root (Eq. (3)), it distributes the maximum amount of strain energy in the fewest possible number of principal components. 15 As a closing remark, we note that in many applications it is customary to subtract the average of the observations prior to performing PCA. This procedure is well suited for systems whose dynamics can be well represented as a fluctuating process superimposed on a steady state. However, this is not the case of chemotaxing amoeboid cells as previous studies have shown cell shape and traction forces show a marked periodic behavior in this type of cell migration. 10,26 Thus, we applied PCA without subtracting the average of the observations, a procedure already reported in the literature, 12 and which in our case led to the maximization of the strain energy accounted for by each mode. We hypothesize that the dominant mode obtained by this procedure will capture the temporal periodicity of cell shape and traction forces. This hypothesis is tested in ''One Principal Component Captures the Temporal Periodicity of Cell Length and Strain Energy During Amoeboid Cell Migration'' section.

Individual and Ensemble PCA
We applied PCA to the time history of hðx; tÞ for 1) each single cell (individual PCA), and 2) for a pooled set of observations coming from all the cells recorded in our experiments (ensemble PCA). The high computational cost of performing ensemble PCA to all the cells (10468 observations from N = 24 cells) led us to consider 10 different sets, each of them containing 600 observations evenly distributed among each cell. Comparing the principal components obtained from individual and ensemble PCA allowed us to determine the cell-to-cell reproducibility of the principal components of shape and strain energy.

Traction Force Components in Wild-type Cells and Reproducibility
We applied PCA (see ''Materials and Methods'') to determine the strain energy components exerted by WT Dictyostelium cells migrating chemotactically. The five most energetic components were highly reproducible from cell to cell, allowing for their visual identification. Their spatial distributions are shown in Fig. 3. We will refer to them as the canonical principal components of the strain energy (CPCs). The reproducibility of these five components was proven by the fact that they were ranked as the dominant ones by ensemble PCA in all data pools analyzed (see Fig. 1   Evolution of Traction Forces During Pseudopod Protrusion and Retraction'' section shows that five principal components provided just enough time resolution to capture the evolution of the traction forces exerted by the cells during migration. Altogether, the five CPCs accounted for m 1À5 $ 60% of the strain energy exerted by migrating, WT cells (Fig. 4), resulting in an average of 12% energy captured per mode. The narrow distribution of m k for ensemble PCA suggested that each of the pools considered were representative of the whole cell line. Because our analysis employed time-lapse recordings containing 600 traction force maps (see ''Individual and Ensemble PCA'' section), PCA identified 600 strain energy components. The data in Fig. 4 indicate that the 595 components other than the CPCs account for remaining 40% of the strain energy exerted WT cells, leading to a an average percent energy of 0.07%, which is 150 times lower than that of the CPCs.
The most important canonical principal component, CPC1, accounts for~40% of the total strain energy exerted by WT cells (Fig. 4). Figure 3a displays CPC1 as a 1 j 50 n 1 ðxÞ; where n 1 ðxÞ is the traction force corresponding to CPC1 (see Eq. (8)) and a 1 j 50 is the median of its weight coefficient. CPC1 depicts the cell attaching at both its anterior and posterior parts and contracting the substrate inwards towards the centroid of the cell. The coefficient a 1 ; which modulated the magnitude of the traction forces in CPC1, was skewed towards positive values in 98.5% of the observations (see Fig. 5a), indicating that CPC1 accounted for a persistent contraction of the substrate from front and back. Unlike the first principal component, the weight coefficients of CPC2-5 were centered around zero (see Figs. 5b-5e). This result motivated the distinction between (+) and (À) configurations of these components (see Figs. 3b1-3e1 and 3b3-3e3). Because of both the high proportion of the strain energy U S accounted for by CPC1 and the marked skewness of its weight coefficient towards positive values, it was useful to determine the perturbation effect that each CPCk (k = 2 À 5) had on the first principal component. and Figs. 3b4-3e4 display the addition of CPCk (À) to CPC1 as a 1 j 50 n 1 ðxÞ þ a k j 10 n k ðxÞ % a 1 j 50 n 1 ðxÞ À a k j 90 n k ðxÞ: ð11Þ The percentiles 10 and 90 of a k were chosen to facilitate the visualization of the effect that CPC2-5 had in the overall traction force pattern, which was dominated by CPC1. CPC2 and CPC5 were respectively responsible for 6% and~3% of the strain energy exerted by WT cells (Fig. 4) and accounted for lateral asymmetries in their traction forces. These components were associated with the bending of cell shape due to projection of pseudopods not aligned with the longitudinal axis of the cell (Figs. 3b and 3e, respectively). The overall effect of CPC2 (+)/(À) (Figs. 3b1/3b3) was to displace the region of application of traction forces towards the right/left at the back of the cell, and towards the left/ right at the front. At the front, the cell shape displaced towards the right/left, not exerting significant forces, thus suggesting that the front was not attached. CPC5 (+)/(À) was similar to CPC2 (+)/(À), but presented an net enhancement of the contractile traction forces at the front (Figs. 3e2 and 4), which suggests that the front was attached to the substrate, in contrast to CPC2.
CPC3 was responsible for~6% of the strain energy exerted by cells (Fig. 4). The (+)/(À) configuration of CPC3 depicted the contraction/stretch of the substrate from center to back and stretch/contraction of it from front to center (Figs. 3c1/3c3). The effect of CPC3 on the traction force pattern was easier to understand when added up to CPC1 (Figs. 3c2/3c4). CPC3 (+)/ (À) displaced the areas of contraction of the substrate  towards the back/front of the cell. The cell shape associated with CPC3 (+)/(À) was a slight elongation/ dilation (see dashed black line in Fig. 3c). CPC4 (Fig. 3d) accounted for~5% of the total strain energy (Fig. 4). The (+)/(À) configuration of CPC4 accounted for an increase/decrease in the strength of the contractile forces at front and back, and a decrease/increase of the contractile lateral forces at the sides of the cell. CPC4 (+)/(À) was associated with the elongation/dilation of the cell shape.

One Principal Component Captures the Temporal Periodicity of Cell Length and Strain Energy During Amoeboid Cell Migration
Previous studies showed that U S and L oscillate periodically and in phase during amoeboid cell migration, 10,26 consistent with the cyclic implementation of anterior protrusions and posterior retractions by the cell. 20 In this section, we tested the hypothesis that these periodic oscillations had one degree of freedom and that they were controlled by the first CPC only. For this purpose, we analyzed the Spearman's rank correlation between the weight coefficient of CPCk, a k ; and either U S or L (r U S Àa k or r LÀa k ; k ¼ 1 À 5), computed from individual PCA (see Fig. 6).
The PDF of r U S Àa 1 was clearly displaced towards high positive values (hr U S Àa 1 i $ 0:9; p<10 À5 ; Fig. 6b), indicating that, as we hypothesized, the temporal evolution of U S was accurately described by the temporal dynamics of CPC1. The weight coefficients a k for components k = 2 À 5 showed no correlation with U S ; collapsing to zero for the asymmetric components k = 2 and 5. The PDF of r LÀa 1 showed a statistically significant positive correlation (p = 0.007, Fig. 6c), which together with the high correlation between U S and CPC1, indicates that the modulation of the inward contraction of the substrate from front and back of the cell was responsible for the positive correlation between U S and L. 10,26 Weight coefficient a 4 and L were also positively correlated (p ¼ 4 Â 10 À5 ), due to the fact that CPC4 was an aspect ratio component that accounted for changes in traction forces related to the dilation and elongation of the cell (Eq. (4), Fig. 3d). The 2D-PDFs of a k (k = 2 À 5) vs. a 1 (Fig. 7) showed that the weight of CPC2-5 was bounded by the weight of CPC1. This result confirmed that the first principal component governed the mechanics of amoeboid motility.

Five Principal Components Resolve the Spatio-temporal Evolution of Traction Forces During Pseudopod Protrusion and Retraction
In the previous section, we showed that the first CPC alone describes the periodic behavior of cell length and strain energy during migration. However, it is known that key events of cell migration are completed in shorter times than the full oscillations of U S and L. Meili et al. 26 reported that both pseudopod protrusion and retraction take on average~22 s in chemotaxing WT cells, whereas the period of the oscillations in cell length was found to be~94 s. This raised two questions: 1) does the first CPC provide enough temporal resolution to represent the mechanical stresses that cells exert on their substrate during migration? and 2) if not, how many CPCs are needed to obtain the required time resolution? We estimated the temporal resolution of each principal component of the strain energy from the temporal memory of the weight coefficients obtained by individual PCA. That is, weight coefficients with long memory have a poor temporal resolution and viceversa. The temporal memory was quantified using the autocorrelation of the weight coefficients of each component, R a k ðsÞ ¼ ha k ðtÞa k ðt þ sÞi: We obtained a parabolic fit to the first 10 s of positive lag of R a k ; and calculated the time t R0 at which that parabola crossed the horizontal axis R a k ¼ 0; providing with an estimated memory time for component k (Fig. 8a). Figure 8b shows that the memory time decreased with the order of the principal component. In particular, the first CPC had a memory time t R0 ðk ¼ 1Þ ¼ 45 s. The corresponding time resolution is too coarse to reproduce the dynamics of events such as protrusion or retraction, which have a duration of 22 s as reported by Meili et al. In order to determine how many principal components are required to resolve these events, we calculated the index k 22 above which t R0 ðkÞ becomes <22 s. Analysis of k 22 revealed that five principal components are required to achieve a memory time of 22 s (Fig. 8c). Thus, we concluded that five principal components are needed to resolve the spatiotemporal evolution of traction forces during pseudopod protrusion and retraction.

Temporal Coordination Between Traction Forces and Lateral Protrusions
The two asymmetric canonical principal components, CPC2 and CPC5, were found to be associated with lateral protrusion of pseudopods. CPC5 accounts for the cell attaching the lateral protrusion to the substrate (Fig. 3b), whereas in CPC2 the lateral protrusion is not attached (Fig. 3b). The aim of this section was to analyze the possible temporal coordination between these two components, and how they are coordinated with the resulting bending of cell shape. The results from this analysis showed that the protrusion of lateral pseudopods takes about 12 s in average, starting with the lateral protrusion leading to the bending of cell shape and finishing with the attachment of the pseudopod to the substrate.
Application of PCA only to cell shape, similar to Killich et al. 19 and Keren et al., 18 provided with a set of cell shape components. One of them accounts directly for cell bending and is temporally described by its weight coefficient a s b ðtÞ: The PDFs of the Spearman's rank correlation coefficient between a s b and either a 2 or a 5 (Fig. 9b) proved that the association between cell shape bending and both CPC2 and CPC5 was positive and statistically significant. We studied the temporal coordination of cell shape bending, CPC2, and CPC5 through analysis of the cross-correlation of their weight coefficients and their temporal delays (Figs. 9a,  9c), showing that CPC5 is delayed~12 s with respect to CPC2 (p = 0.004) (Figs. 9a, 9c, first boxplot) and that the physical bending of the cell precedes CPC5 by~7 s (p = 0.0015) (Figs. 9a and 9c, second boxplot). Although Fig. 9c (third box-plot) suggests that CPC2 precedes the bending of the cell by~4 s, the statistical significance was marginal (p = 0.09), probably because this delay was close to the temporal resolution of our experiments, 2 s.

Traction Force Components in Contractility Mutants are Similar to those in Wild-type Cells
To investigate the role that MyoII contractility plays in the strain energy components we applied PCA to recordings of 1) MyoII essential light chain null cells with altered MyoII motor function 7 (mlcE À , N = 13), and 2) MyoII null cells lacking MyoII cross-linking and motor function 31 Applying both ensemble and individual PCA (see ''Materials and Methods'') to recordings of traction forces and cell shape for both mhcA À and mlcE À cells, we identified five CPCs for each mutant cell line. These principal components, specially in mlcE À ; were remarkably similar to those identified for WT cells (Fig. 14 for mlcE À ; and 15 for mhcA À in the Supporting Material), and accounted for~55% of U S exerted by mlcE À and mhcA À (Fig. 13A in the Supporting Material). As in WT cells, CPC1 accounted for a substantial proportion of the strain energy (~40%) in both mlcE À and mhcA À cells, indicating that temporal evolution of the mechanical interaction between contractility mutants and their substrate was also dynamically simple in these mutants. However, in mhcA À cells the principal components with k > 1 were less reproducible from cell to cell (Fig. 13C in the Supporting Material) than in WT, indicating the difficulty to define a clear ''mechanical phenotype'' for this cell line.
The first principal component in mlcE À was very similar to that in WT cells (Figs. 3a1, 10a), with traction forces focused at two areas at the posterior and anterior parts of the cell. However, this was not the case in mhcA À cells, where CPC1 presented a large amount of lateral forces (Fig. 10b). The prominent exertion of lateral forces by mhcA À is present throughout all the components (Fig. 15 in the Supporting Material).

DISCUSSION
Dr. Chien has been a major force in laying the foundation of mechanobiology, an emerging branch of cell biology that studies the interplay between mechanical forces and molecular biology. An essential element of mechanobiology is the variety of processes by which the cell responds to mechanical stresses, remodels its cytoskeleton and regulates its function in physiology and disease. Among his many contributions, Dr. Chien has studied how external forces applied to the cell affect cell shape or cell migration. 8,16,21,35,38,42 The present paper has been greatly inspired by Dr. Chien's pioneering work in this area.
In motile cells, shape changes are strongly related with the modulation of cell adhesions to the substrate via the forces they exert on it. 6,13,26,34,39 In order to study this modulation, we measured the traction forces and cell shape of Dictyostelium cells migrating chemotactically over a flat elastic substrate. We then applied a tailored formulation of Principal Component Analysis (PCA) that optimized the mechanical work exerted by these cells. In addition, we used a 2D mapping of cell shape to obtain the different shape changes associated with each strain energy component.
Principal component analysis determined a set of five canonical principal components of the strain energy (CPCs) that were reproducible from cell to cell, (a) (b) (c) FIGURE 9. Analysis of the temporal coordination between cell shape bending and the two asymmetric canonical principal components of the strain energy, CPC2 and CPC5. (a) For a WT cell, cross-correlations between the bending of the cell (weight coefficient a s b ), CPC2 (weight coefficient a 2 ) and CPC5 (weight coefficient a 5 ), used to quantify the delay between the signals (indicated with solid circles). (b) Box-plot of the Spearman's rank correlation coefficient between: a 2 and a s b ; and a 5 and a s b : (c) Boxplot of the time delay, in s, between different weight coefficients, as extracted from their cross-correlations (see panel (a)) (a x ! a y = time delay of a y with respect to a x ). Stars indicate significance of testing the null hypothesis that the statistical distributions in panels (b-c) come from a continuous distribution with zero median: *p < 0.05 and **p < 0.01. and which captured~60% of the total mechanical work exerted by cells. Furthermore, the five canonical components were shown to provide enough temporal detail to resolve the mechanics of pseudopod protrusions and retractions. This result opens a promising venue to the development of reduced models because it reveals that the complex dynamics of amoeboid cell migration can be encapsulated within the relatively simple dynamics of just five degrees of freedom. Figure 11 and Movie 1 in the Supporting Material illustrates this idea by comparing measurements of the time evolution of the shape and traction forces in a chemotaxing WT Dictyostelium cell with their projections onto the 5 CPCs. The agreement between the full measurements and the data obtained from the 5 CPCs is remarkable. Furthermore, the restriction to 5 CPCs reduces the noise of the traction force measurements. The most important of these five components accounted for~40% of U S ; presenting an oscillatory temporal evolution which correlated with both the strain energy and global cell shape. This dominant component depicted the cell attaching at front and back, contracting the substrate inward towards its center, and it is consistent with previous traction force measurements. 10,23,26 The rest of the components were substantially less important in terms of mechanical energy, accounting for less than 7% each.
Keren et al. 18 showed that a low number of components described the shape variability of migrating keratocytes within a population. In their cells, the most important component behaved steadily in time. Our study illustrates that PCA is also able to compress the more complex dynamics of amoeboid migrating cells using a reduced number of degrees of freedom with coherent temporal evolution. For instance, our PCA identified two modes related with the initiation and termination of lateral protrusions. These components were separated by a period of~12 s, during which an increase of traction forces revealed the formation of new adhesions.
The canonical principal components were found to be similar in the case of mutants with contractile deficiencies (mlcE À and mhcA À ), although these components were less reproducible in mhcA À cells, indicating they have a less clear or less coherent ''mechanical phenotype''. The first component in mhcA À cells presented high lateral forces, suggesting the proliferation of lateral adhesions and cortical contractility as the origin of the traction forces in the absence of MyoII.
We have presented here the basics of an improved PCA analysis that allows for the dissection of cell motility into a reduced set of principal components. The application of this PCA method to mutant cell lines with adhesion, contractility, and actin polymerization defects represents a very attractive venue to improve our understanding on the role that the different cytoskeleton and adhesion proteins play in the generation of traction forces and the work needed for locomotion.

ELECTRONIC SUPPLEMENTARY MATERIAL
The online version of this article (doi:10.1007/ s12195-011-0184-9) contains supplementary material, which is available to authorized users.
(a) (b) FIGURE 11. Comparison between instantaneous cell shape and traction force fields exerted by two chemotaxing cells (a) and reconstruction obtained with the 5 CPCs obtained from individual PCA (b). The colormap indicates the magnitude (colorbars in Pa) and the arrows the direction of the traction forces. The solid black contours represent measured (a) and reconstructed cell contours (b). The blue line indicates the cell trajectories, and the black bar is 10 lm long. Reconstructed cell contours were obtained as the contour of the reconstructed cell shape field function whose inscribed area equaled the instantaneous cell area.