Search tips
Search criteria 


Logo of bmcbiophysBioMed CentralBiomed Central Web Sitesearchsubmit a manuscriptregisterthis articleBMC BiophysicsJournal Front Page
BMC Biophys. 2012; 5: 8.
Published online May 11, 2012. doi:  10.1186/2046-1682-5-8
PMCID: PMC3414788
Differences in adhesion and protrusion properties correlate with differences in migration speed under EGF stimulation
Yue Hou,1 Sarah Hedberg,3 and Ian C Schneidercorresponding author1,2
1Department of Chemical and Biological Engineering, Iowa State University, Iowa, USA
2Department of Genetics, Development and Cell Biology, Iowa State University, Iowa, USA
3Department of Biochemical Engineering, University College London, London, UK
corresponding authorCorresponding author.
Yue Hou: arielyh/at/; Sarah Hedberg: s.hedberg/at/; Ian C Schneider: ians/at/
Received December 23, 2011; Accepted May 11, 2012.
Cell migration plays an essential role in many biological processes, such as cancer metastasis, wound healing and immune response. Cell migration is mediated through protrusion and focal adhesion (FA) assembly, maturation and disassembly. Epidermal growth factor (EGF) is known to enhance migration rate in many cell types; however it is not known how FA maturation, FA dynamics and protrusion dynamics are regulated during EGF-induced migration. Here we use total internal reflection fluorescence (TIRF) microscopy and image analysis to quantify FA properties and protrusion dynamics under different doses of EGF stimulation.
EGF was found to broaden the distribution of cell migration rates, generating more fast and slow cells. Furthermore, groups based on EGF stimulation condition or cell migration speed were marked by characteristic signatures. When data was binned based on EGF stimulation conditions, FA intensity and FA number per cell showed the largest difference among stimulation groups. FA intensity decreased with increasing EGF concentration and FA number per cell was highest under intermediate stimulation conditions. No difference in protrusion behavior was observed. However, when data was binned based on cell migration speed, FA intensity and not FA number per cell showed the largest difference among groups. FA intensity was lower for fast migrating cells. Additionally, waves of protrusion tended to correlate with fast migrating cells.
Only a portion of the FA properties and protrusion dynamics that correlate with migration speed, correlate with EGF stimulation condition. Those that do not correlate with EGF stimulation condition constitute the most sensitive output for identifying why cells respond differently to EGF. The idea that EGF can both increase and decrease the migration speed of individual cells in a population has particular relevance to cancer metastasis where the microenvironment can select subpopulations based on some adhesion and protrusion characteristics, leading to a more invasive phenotype as would be seen if all cells responded like an “average” cell.
Cell migration plays an important role in tumor progression [1]. During invasion and metastasis, migration is driven by soluble extracellular cues like epidermal growth factor (EGF). EGF is a well-known chemoattractant [2-4]; however, uniform doses also stimulate chemokinetic responses. EGF’s control of cell motility originates from its regulation of adhesion and protrusion [5-7]. This occurs through altering adhesive attachments called focal adhesions (FAs) [8-10] as well as actin cytoskeleton organization.[8,11,12] The response to EGF at the level of cell migration is dose dependent, but there exists a range of maximal stimulation concentrations. Often migration saturates at 2–10 nM EGF [2,13,14], but some of the studies showed an inhibition of migration at EGF concentrations >2-10 nM [7,15]. This is in agreement with other work demonstrating that in certain contexts, EGF can inhibit migration [16,17]. Within each study there is wide diversity in migration behavior, even among cells observed during the same experiment [14]. Interestingly, the distribution in migration speed and persistence time appears to be dependent on EGF stimulation [18], suggesting that EGF controls not only the mean response, but also the amount of cell-to-cell variability. Cell-to-cell variability has been widely observed, and has drawn much attention due to its influence on physiology [19], pathology[20] and pharmacology [21]. Consequently, mathematical models [19,20] have been used to show that even small changes in the distribution of protein concentrations yield enhanced wound healing or metastasis due to the selection of an optimal subpopulation. When the subpopulation is defined based on migration speed, it will not only be beneficial to examine the distribution of protein concentrations, but also higher level characteristics like focal adhesion (FA) properties and cytoskeletal dynamics.
FAs are dynamic, macromolecular structures that serve as both mechanical linkages and centers of intracellular signal transduction [22-24]. They assemble as nascent adhesions, mature into focal complexes, focal adhesions and fibrillar adhesions and disassemble [23]. Consequently, FAs exhibit different morphological maturation states throughout their lifetime and this is thought to regulate their behavior. For example, small, nascent FAs, transmit strong forces and serve as traction points for propulsive forces to move the cell body forward [25,26]. They also generate signals for protrusion by activating actin accessory proteins [27-31]. Under tension, these small FAs can mature into larger focal complexes, focal adhesions and fibrillar complexes with different force transmission characteristics and propensities for protrusion signaling [32-34]. Several morphological characteristics have been used to predict traction force and cell migration speed including FA protein density, number per cell, sliding speed, lifetime, size and elongation [22,23,25,33-35]. These morphological characteristics have begun to be quantitatively measured [36,37] and the distributions properly quantified [38]. However, their direct correlation to migratory states as well as their response to extracellular cues like EGF is unknown.
Protrusion is mediated by actin polymerization, whereas retraction is driven through myosin II activity and actin depolymerization [39,40]. Protrusion and retraction can either occur continuously in spatially confined regions as in keratocyte migration or it can occur in cycles or waves of protrusion that move laterally along the edge [41-43]. This has been characterized in several cell types when cells are either spreading [44] or migrating [41-43]. In fact a recent paper has shown that slower migrating keratocytes employ lateral protrusion waves [43]. While the timing of the cycles and the propagation of the waves is dependent on intracellular pathways, very little work has been done to examine how protrusion is quantitatively altered in response to extracellular stimuli like EGF.
To understand the relationship between EGF-stimulated cell migration, FA properties and protrusion dynamics, we imaged metastatic (MTLn3) and non-metastatic (MTC) cell lines. We analyzed the cell migration speed and persistence under various EGF stimulation conditions and found that EGF moderately increased the median migration rate and persistence of MTLn3 cells, whereas it had no significant effect on the speed and persistence of MTC cells. Interestingly, higher concentrations of EGF broadened the distributions and increased the coefficient of variation of both the migration rate and persistence of MTLn3 cells, but not MTC cells. When data was binned based on EGF stimulation conditions, FA intensity and FA number per cell showed the largest difference among stimulation groups. FA intensity decreased with increasing EGF concentration and FA number per cell was highest under intermediate stimulation conditions. No difference in protrusion behavior was observed. However, when data was binned based on cell migration speed, FA intensity and not FA number per cell showed the largest difference among groups. FA intensity was lower for fast migrating cells. Additionally, waves of protrusion tended to correlate with fast migrating cells. Consequently, low FA intensity and waves of protrusion are markers for fast migrating cells, but these characteristics are only partially predictive of EGF stimulation conditions because of the large cell-to-cell variability in response to EGF.
EGF stimulation broadens the distributions of migration speed and persistence of MTLn3 cells
In many cell types EGF has been reported to enhance the mean migration speed. However, the long term migration response of individual cells after challenge with EGF in this model system is not known. Consequently, we examined cell migration under various doses of EGF in both adenocarcinoma cells (MTLn3) and non-metastatic cells taken from the same tumor (MTC) (Figure 1). In vivo, EGF concentrations in serum can be between 0.1-2 nM, with local tissue concentrations as high as ~20 nM [13]. Consequently, MTLn3 and MTC cell lines were stimulated with a wide range of EGF concentrations (0–100 nM). Cell tracks were generated (Additional file 1: Figure S1) and cell speed and persistence were calculated based on fitting the mean squared displacement with a random walk model [45,46]. Contrary to previous reports in other cell lines, EGF stimulation only increased the median speed and decreased the median persistence slightly in MTLn3 cells and acted more like an on-off switch between no EGF stimulation and EGF stimulation (Figure 1A and B). In addition, there was no dose response in either median speed or persistence of MTC cells (Figure 1A and B). Cell persistence decreased somewhat with increasing cell speed in MTLn3 cells, but MTC showed no such correlation, populating a much lower range of migration speeds (Figure 1C) with roughly the same range of persistence times. Interestingly, at higher EGF concentrations, MTLn3 cells both migrated faster and slower than those at lower EGF concentrations (Figure 1A). On the other hand, MTC cells did not show this same behavior. The coefficient of variation (standard deviation/mean) increased in a dose dependent manner for both cell speed and persistence in MTLn3 cells (Figure 1 D and E).
Figure 1
Figure 1
Increasing EGF broadens the distribution of both cell speed and persistence of MTLn3 cells while not changing MTC cell migration. A. cell speed and (um/hr) in the axis description overlaps. B. Cell persistence of MTLn3 (left) and MTC (right) under different (more ...)
Given that EGF did not dramatically affect the median migration, we grouped cells according to no (0 nM EGF), low (0.01 and 0.1 nM) and high (1, 10 and 100 nM) EGF stimulation. Additionally, since the variability in cell migration speed seems to be an additional feature of the data, we also grouped cells based on cell migration speed. A k-means clustering algorithm for a cluster number equal to two was applied to the migration speeds of all cells under different EGF concentrations. The cutoff speed between slow and fast migration cells was found to be 42 μm/hr. Thus we assigned cells with speeds of greater than 42 μm/hr to the fast migrating group and cells with speeds of less than 42 μm/hr to the slow migrating group. Having grouped cells in two different ways, we wanted to examine FA characteristics and protrusion dynamics to see if certain signatures were exhibited by EGF stimulated cells or fast moving cells.
The distribution of FA characteristics differ under different EGF stimulation conditions and between fast and slow migrating cells
In order to measure the FA characteristics of FA intensity, number per cell, size, speed, lifetime and elongation, we used TIRF microscopy to observe FAs within a single cell transfected with paxillin-EGFP, a component that marks FAs throughout their entire lifetime. Expression of paxillin-EGFP did not seem to alter the migration speed, nor was the expression dependent on EGF stimulation (Additional file 2: Figure S2). In MTLn3 cells, many FAs assembled, matured and disassembled over several minutes, so images were taken every 10s for 40–60 minutes (Figure 2A), with little influence of photobleaching (Additional file 3: Figure S3). A segmentation and tracking algorithm was used to quantify FA characteristics (Figure 2B-D) and time-resolved data of focal adhesion characteristics ( Additional file 4: Figure S4) [37]. We categorized by eye the FA tracking results of different EGF concentrations and scored them as poor, good and excellent. The example shown in Figure 2 was scored as good. Most tracks at each EGF condition resulted in more than 70 % that were either good or excellent (Table 1).
Figure 2
Figure 2
Time-lapse series of FA dynamics in MTLn3 cells. MTLn3 cell expressing paxillin-EGFP and stimulated with 0.01 nM EGF is shown. Images were taken at 0, 10, 20, 30 and 40 min. A. Original total internal reflection fluorescence (TIRF) images of cells (more ...)
Table 1
Table 1
Qualitative assessment of tracking results
FA characteristics can be ordered based on the magnitude in the difference between either the no, low and high EGF stimulation conditions or between slow and fast migrating cells. This magnitude was quantified by the Kolmogorov-Smirnov statistic (Figure 3). When this statistic is large, it is more likely that there is a difference in distributions between groups. FA intensity and number per cell showed the largest values, so we decided to focus on these characteristics. Distributions of all other FA characteristics are shown as Additional file 5: Figure S5 and Additional file 6: Figure S6 and a summary of relevant ranges of these FA characteristics under different conditions is shown in Table 2. Most FA characteristics fit best to either lognormal or Weibull probability distribution functions. As EGF concentration increased from no to low to high, FA intensity decreased (Figure 4). FA number per cell on the other hand showed highest numbers at low concentrations of EGF. Both FA intensity and number per cell showed strong differences between EGF concentration groups (Figure 3). When cells were grouped based on migration speed, FA intensity was lower for fast migrating cells (Figure 5). However, distributions for FA number per cell were now much less, indicating that this characteristic shows poorer correlation with migration speed (Figure 3). Having identified some FA characteristics that correlate with either EGF stimulation conditions or cell speed, we were interested if any protrusion characteristics showed difference among groups.
Figure 3
Figure 3
Quantification of the difference between experimental distributions of FA properties. Cells were binned based on EGF concentration (no (0 nM), low (0.1 and 0.01 nM) and high (100, 10 and 1 nM)) or cell speed (low (<42 μm/hr) or (more ...)
Table 2
Table 2
Summary of FA characteristics in fast migrating cells and those stimulated with low and high EGF concentrations
Figure 4
Figure 4
FA intensity decreases with increasing EGF concentration and number per cell is maximal at intermediate EGF concentrations. Histograms of FA intensity A.-C. and number per cell E.-G. were generated by dividing cells into three EGF stimulation groups (A., (more ...)
Figure 5
Figure 5
FA intensity is lower in fast migrating cells and FA number per cell is slightly higher in fast migrating cells. Histograms of FA intensity A.-B. and FA number per cell D.-E. were generated by dividing cells into two cell migration speed groups (A., D. (more ...)
Unique spatial organization of protrusion and retraction is exhibited in fast migrating cells
We assessed differences in protrusion and retraction behavior under different EGF stimulation conditions and between the slow and fast migrating cells. The cell-average protrusion and retraction velocities would be faster in the fast migrating cells resulting in the increased migration rate; however different patterns of protrusion could lead to the same average value, so we analyzed local protrusion behavior. One prominent feature that we observed was traveling waves of protrusion along the edge of the cell. This traveling wave behavior had the effect of broadening of the protrusion velocity distribution. Upon qualitative examination, traveling waves did not seem to be linked to EGF stimulation conditions. Additionally, slow migrating cells usually showed large quiescent areas (green) and random, disorganized protrusion and retraction behavior (Figure 6A). Only fast migrating cells showed traveling waves of protrusion (Figure 6B). We examined the difference in protrusion velocity distributions in the same way that we examined distributions of FA properties (Figure 3). The Kolmogorov-Smirnov statistic was smaller when comparing EGF stimulation conditions than it was when comparing cell migration speeds (Figure 6C). This suggested that waves as described by a wide protrusion velocity distribution correlate with differences in migration speed and not EGF concentration. We computed the fraction of cells with waves and measured the standard deviation of the protrusion velocity distribution in slow and fast migrating cells. Both waves and high standard deviations were features of fast migrating cells (Figure 6D, E and F). Consequently, fast cells tend to organize their protrusion in a qualitatively different way than slow migrating cells and this does not necessarily correlate with EGF stimulation.
Figure 6
Figure 6
The spatial control of protrusion differs between slow and fast migrating cells. A. Protrusion velocity map for slow migrating cells at 0.01, 1, and 100 nM EGF from left to right. B. Protrusion velocity map for fast migrating cells at 0.01, 1, and 100 (more ...)
Variability in cell response to environmental cues is becoming a more appreciated phenomenon that can drive how populations of cells respond to their environment. Cell-to-cell variability can arise from heterogeneity in protein level [47,48] or organization of cellular structures such as the membrane [49] or the cytoskeleton [50]. Interestingly, this variability can be enhanced by extracellular stimuli [51]. The idea that variability can be enhanced under certain conditions sets up the interesting possibility that the mean response is a relatively poor statistical metric. Rather, the distribution itself or the standard deviation or another parameter that characterizes the distribution may be more appropriate. The obvious result of this dependence on the distribution is a sensitizing of a subpopulation of cells to particular environments. This is acutely evident in pathologies such as cancer metastasis, where subpopulations of cells are selected based on differing responses to the tumor microenvironment. Therefore, the fastest cells most likely drive metastasis, whereas the average cell migration rate might be less important. We showed that the distribution of cell migration speed and persistence is very much regulated under EGF stimulation, even though the average response differs marginally. Indeed, this has been demonstrated previously. [18] Ware et al. generated distributions of migration rate in response to no EGF or high EGF concentration. However, the focus of that paper was primarily on the changes in the average migration response and the widening of the distribution in response to EGF was evident, but not discussed. What causes this widening? Heterogeneity in the local ECM concentration might play a role. We have examined collagen coverage and it tends to be fairly homogeneous at the resolution of the light microscope (~100 nm) and we observed cells in close proximity that varied greatly with respect to their migration speed. However, ECM inhomogeneity cannot be fully dismissed as a possible cause for the cell-to-cell variability. Another cause of the cell-to-cell variability might be autocrine or paracrine signaling. MTLn3 cells are known to secrete other EGF receptor ligands, namely TGF-α [52]. However, we did not observe clustering of migration speeds around sources. Often cells in the same clusters showed distinct behavior. A third possibility is that concentrations of signaling, adhesion or cytoskeletal regulatory proteins might contribute to the heterogeneity. This might be the most probable cause of the cell-to-cell variability; however determining which specific components might contribute to this is the subject of further investigation.
EGF does seem to regulate some FA characteristics, namely FA intensity and number per cell. FA intensity decreases as EGF stimulation increases. FA number per cell is highest at low EGF concentrations, suggesting that either the assembly is maximized or disassembly is minimized at this point. EGF is known to alter actin cytoskeleton dynamics, perhaps resulting in enhanced assembly dynamics. Alternatively, EGF is also known to upregulate calpain, a protease involved in disassembly, which might be largest at high EGF concentrations. Interestingly, FA number per cell does not seem to affect cell migration rate. Rather, low FA intensity seems to correlate with fast migrating cells. This might be a direct link between EGF stimulation and cell migration speed regulation (Figure 7). Zaidel-Bar et al. observed that the localization of paxillin in large FAs did not affect the rate of protrusion of the nearby lamella. However, paxillin association with focal complexes was inversely correlated with the rate of local protrusion. Thus, focal complexes containing relatively low levels of paxillin were found in fast protrusions [29]. While seemingly less important as measured by the Kolmogorov-Smirnov statistic, we also found that fast migrating cells contain FAs with intermediate sizes and FAs with intermediate sliding speeds and FAs with short lifetimes (Table 2 and Additional file 6: Figure S6). The fact that fast migrating requires intermediate sized FAs (0.3-3 μm2) is not surprising. Focal complexes, small FAs are traditionally thought to occupy this range of areas [24]. These smaller FAs are usually located near the leading edge and transmit strong propulsive traction forces needed during fast migration. Larger, mature FAs exert weaker forces [25] and “supermature” fibrillar adhesions [53] are involved in ECM remodeling, both processes that are typically seen in slower migrating cells. Fast migrating cells also contained FAs with intermediate speed. FA speed affects cell speed in complicated ways due to its spatial regulation. For example, Smilenov et al. found that fibroblasts with stationary FAs tend to transmit large forces and result in migratory cells [35]. However, Diener et al. found that FAs moved with a sliding speed of 4 μm/hr in migrating human osteosarcoma cells on collagen-coated coverslips [54] and FAs at the trailing edge are pulled forward at rates of >5-10 μm/hr [55]. Lifetime was minimal in fast migrating cells. Others have shown that FAs with short lifetimes correlate with fast migrating cells, in line with what we observe [28].
Figure 7
Figure 7
Different adhesion and protrusion characteristics correlate with EGF stimulation and cell speed. Cells can either be grouped based on EGF concentration or cell speed. EGF concentration is considered an input that acts to regulate adhesion and protrusion (more ...)
Given that local protrusion is linked to FA intensity and that FA intensity was lowest in fast migrating cells we were encouraged to examine the protrusion dynamics under different stimulation conditions. These cells are known to respond acutely to EGF stimulation with two peaks of barbed end formation resulting in a robust protrusion response [56]. However, cells are often not exposed to these acute signals in vivo and so we asked how protrusion changes under chronic EGF stimulation. We found that while EGF stimulation condition correlated poorly with lateral waves generated in cells, fast migrating cells usually generated lateral waves of protrusion as has been seen elsewhere [41-43,57].The existence of lateral protrusion waves suggests locally activated feedback loops that travel laterally along the edge of the cell [43,58]. This positive feedback loop operates through adhesion signaling for protrusion and protrusion resulting in more adhesions [59]. How does this behavior relate to migration rate? Barnhart et al. noticed that keratocytes migrating on more adhesive substrates generated these lateral waves and migrated with a slower speed [43]. We see an opposite relationship, where high speeds result in lateral waves of protrusion. This difference may be related to the differences in cytoskeletal organization and in fact morphology between these cells. Keratocytes adopt highly regular persistent cytoskeletal structure and cellular morphology resulting in extremely fast migration speeds (500–600 μm/hr). MTLn3 cells on the other hand have a varied cytoskeletal structure and cellular morphology and are much slower (<100 μm/hr). Consequently, highly organized, persistent protrusion that is seen in keratocytes results in the fastest migrating cells. Less efficient, but somewhat organized lateral protrusion seen in both keratocytes and MTLn3 cells results in intermediate speeds. Poorly organized protrusion seen in MTLn3 cells results in slow speeds. Local differences in ECM in our system might explain why EGF is not a primary driver for fast migrating cells or lateral protrusion waves, leading high cell-to-cell variability (Figure 7).
EGF was found to broaden the distribution of cell migration rates, generating both faster and slower cells, but not dramatically affecting the average response. Several different adhesion and protrusion characteristics correlated with EGF stimulation and cell migration speed, however there is a hierarchy of these correlations. FA intensity and number per cell correlate with EGF stimulation conditions. FA intensity decreases with increasing EGF stimulation and FA number per cell is highest at low EGF stimulation conditions. In contrast, FA intensity and not number per cell as well as protrusion waves correlate with cell speed. Fast cells are marked by low FA intensity and protrusion waves. Consequently, while EGF stimulation could regulate FA intensity to modulate cell speed directly or by partially activating protrusion waves, other environmental factors most likely lead to protrusion waves. Adhesion and protrusion characteristics that do not correlate with EGF stimulation condition but do correlate with cell migration speed constitute the most sensitive outputs for identifying why cells respond differently to EGF. The idea that EGF can both increase and decrease the migration speed of individual cells in a population has particular relevance to cancer metastasis where the microenvironment can select subpopulations based on some adhesion and protrusion characteristics, leading to a more invasive phenotype as would be seen if all cells responded like an “average” cell.
Cell culture media was α-MEM medium with L-glutamine (Invitrogen) containing 5 % fetal bovine serum (Invitrogen) and 1 % penicillin-streptomycin (Invitrogen). Collagen and poly-L-lysine (PLL) solution contained 3 μg/ml of rat tail collagen I (Invitrogen) and 2 μg/ml of poly-L-lysine hydrochloride (Sigma), dissolved in 0.5 M acetic acid (Fisher) and sterilized under ultraviolet light for 30 minutes. Serum free imaging media was α-MEM medium without phenol red (Invitrogen) containing 1 mg/ml bovine serum albumin (Sigma), 12 mM HEPES (Fisher), and 1 % penicillin-streptomycin (Invitrogen), adjusted to pH 7.4 and filtered through 0.22 μm pore size filter (Millipore, Fisher).
Cell culture
Rat mammary adenocarcinoma cell lines (metastatic MTLn3 and non-metastatic MTC) were obtained from Dr. Jeffrey E. Segall (Albert Einstein college of Medicine). Cell lines were derived from the 13762NF rat mammary adenocarcinoma tumor. Cells were maintained in cell culture media at 37 °C in 5 % CO2 and were passed every 2 or 3 days. Collagen and PLL solution was incubated on 22 × 22 mm squeaky cleaned coverslips (Corning, Fisher) at room temperature for 1 hour. Cells were seeded on coverslips with collagen and PLL and incubated for 24 ~ 48 hours at 37 °C in 5 % CO2 (50,000 ~ 100,000 cells/coverslip).
Cell migration assay
MTLn3 and MTC cells were incubated on coverslips with collagen and PLL for 48 hours and were switched to serum free imaging media for 2 hours. Coverslips were mounted onto glass slide chambers in serum free imaging media with different concentrations of EGF (0, 0.01, 0.1, 1, 10 and 100nM EGF). Chambers were maintained at 37 °C for 2 hours and then imaged on a heated stage every 2 minutes for 8 hours. Phase contrast time-lapse images were captured at 20× (NA 0.50, Nikon) with a charge-coupled device (CoolSNAP HQ2, Photometrics) attached to an inverted microscope (Eclipse Ti, Nikon). Cell centroids were identified and tracked manually by MTrackJ plugins of ImageJ. Single cell speed, S, and directional persistence time, P, were obtained by fitting these to the persistent random walk equation [45,46]:
equation M1
using a non-linear least squares regression analysis. The sampling time is every two minutes for 6–8 hours. The mean-squared displacement was constructed using non-overlapping time intervals. Consequently, the model was fitted to data up to a 30 min time lag due to the small number of displacements (<12-16) at time lags greater than 30 min. To quantify protrusion rate we used a constrained optimization program to measure the protrusion and retraction rates from masked images as done previously [42]. The cell edge was segmented into 100 sectors. The average protrusion rate in these sectors was calculated over time.
Fluorecence imaging
MTLn3 cells were incubated on coverslips with collagen and PLL for 24 hours and transfected with paxillin-EGFP and Fugene 6 (Roche) according to the manufacturer’s protocol (6 μl of Fugene 6 and 3 μg of EGFP-paxillin). After one hour transfection, the media was changed to cell culture media and the transfected cells were maintained at 37 °C in 5 % CO2 for 23 hours. Then the cells were switched to serum free imaging media for 2 hours. Coverslips were mounted onto glass slide chambers in serum free imaging media with different concentrations of EGF (0, 0.01, 0.1, 1, 10 and 100 nM EGF). Chambers were maintained at 37 °C for 2 hours and then imaged on a heated stage every 10 seconds for 40 ~ 60 minutes. TIRF images were captured at 60× oil objective (NA 1.49, Nikon) equipped with a TIRF illuminator and fiber optic-coupled laser illumination. The 488 nm laser line of an air-cooled tunable Argon laser (Omnichrome Model 543-AP-A01, Melles Griot) was reflected off a dichroic mirror (89000 ET-QUAD, Chroma). Camera and shutter were controlled by μManager 1.3. An automated segmentation and tracking algorithm was utilized for large-scale analysis of FA dynamics [37]. FAs smaller than 0.05 μm2 and larger than 10 μm2 were excluded from our analysis because they represent either FAs consisting of less than three pixels or several FAs clustered together. FA fluorescence intensities were calibrated to the standard condition of 1 mW laser power with a 300 ms exposure time, so FA intensity should be directly proportional to protein level across all samples.
Statistical analysis
All graphs and statistical analyses were done using JMP and Matlab software. Distributions of FA properties were constructed in the following ways. FA number as described in the results section is more precisely a FA number per cell and consequently the distribution was generated by using the calculated FA number per cell at each time point during the experiment for each cell. Consequently, the number of measurements of FA number per cell is the product of the average number of frames and the cell number. All the other distributions of FA properties were generated by using the time-averaged FA property for each FA in each cell. Consequently, the number of measurements of FA properties is the product of the average FA number and the cell number. Differences between conditions under various EGF concentrations and cell migration speeds were quantified by calculating the Kolmogorov-Smirnov statistic using the Matlab function . Model distributions were fitted by minimizing the Kolmogorov-Smirnov statistic using the Matlab function between the experimental distribution and the model distribution. To determine the statistical differences of the mean values between conditions under various EGF concentrations and cell migration speeds, a student’s t-test was utilized and a p < 0.01 was deemed significant.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
YH carried out the bulk of the experiments. SH conducted the long term migration experiments. YH and ICS analyzed the results and wrote the paper. All authors read and approved the final manuscript.
Additional file 1
Figure S1.Migration trajectories of typical cells under different EGF concentrations. Three cell tracks were chosen randomly under each EGF concentrations and labeled with different colors. All trajectories were aligned to the starting point (0, 0).
Additional file 2
Figure S2.Paxillin-EGFP expression levels for cells with different cell speeds and under different EGF stimulation conditions. Mean intensity of individual cells as a function of A. cell speed, and B. EGF concentrations. Ncell = 56.
Additional file 3
Figure S3.Average paxillin-EGFP intensity in cells as a function of time. A. Mean FA intensity and B. mean intensity of the whole cell for all cells as a function of time. Mean FA Intensity: Ncells = 55. Mean Intensity Ncells = 53.
Additional file 4
Time-resolved data of mean FA intensity and FA number per cell. A. Mean FA intensity and C. FA number per cell for the cell shown in Figure 2 under 0.01 nM EGF stimulation as a function of time. B. FA mean intensity and D. number per cell of five typical cells under different EGF concentrations as a function of time.
Additional file 5
Figure S5.Distributions of FA size, speed, lifetime and elongation for different EGF stimulation conditions. Histograms of FA size A.-C., speed E.-G., lifetime I.-K. and size M.-O. were generated by dividing cells into three EGF stimulation groups (A., E., I. and M., no EGF (0 nM), B., F., J. and N. low EGF (0.01 and 0.1 nM) and C., G.,K. and O. high EGF (1, 10 and 100 nM)). Histograms were fitted with lognormal probability distributions except for speed, which was fitted with a Weibull probability distribution. The mean values of FA D. size, H. speed, L. lifetime and P. elongation are also shown. The number of measurements of FA properties is the product of the average FA number and the cell number. Size, speed, lifetime and elongation: Ncell,no = 5, NFA,no = 1963, Ncell,low = 20, NFA,low = 13,024, Ncell,high = 30, NFA,high = 14,648. Error bars are 95 % confidence intervals and asterisks denote p < 0.01.
Additional file 6
Figure S6.Distributions of FA size, speed, lifetime and elongation for slow and fast migrating cells. Histograms of FA size A.-B., speed D.-E., lifetime G.-H. and size J.-K. were generated by dividing cells into slow and fast migrating groups (A., D., G. and J. slow migrating cells and B., E., H. and K. fast migrating cells). Histograms were fitted with lognormal probability distributions except for speed, which was fitted with a Weibull probability distribution. The mean values of FA C. size, F. speed, I. lifetime and L. elongation are also shown. The number of measurements of FA properties is the product of the average FA number and the cell number. Size, speed, lifetime and elongation: Ncell,slow = 21, NFA,slow = 12,773, Ncell,fast = 34, NFA,fast = 16,842. Error bars are 95 % confidence intervals and asterisks denote p < 0.01.
We thank the Roy J. Carver Charitable Trust for funding ICS and YH. We also thank Dr. Jeffrey Segall for cell lines, Dr. Antonio Sechi and Thomas Wurflinger for help with the FA analysis; Dr. Gaudenz Danuser and Shann-Ching Chen for help with the protrusion analysis and Nick Romsey for help with the cell migration analysis.
  • Friedl P, Wolf K. Plasticity of cell migration: a multiscale tuning model. J Cell Biol. 2010;188(1):11–19. doi: 10.1083/jcb.200909003. [PMC free article] [PubMed] [Cross Ref]
  • Segall JE, Tyerech S, Boselli L, Masseling S, Helft J, Chan A, Jones J, Condeelis J. EGF stimulates lamellipod extension in metastatic mammary adenocarcinoma cells by an actin-dependent mechanism. Clin Exp Metastasis. 1996;14(1):61–72. doi: 10.1007/BF00157687. [PubMed] [Cross Ref]
  • Sawyer C, Sturge J, Bennett DC, O'Hare MJ, Allen WE, Bain J, Jones GE, Vanhaesebroeck B. Regulation of breast cancer cell chemotaxis by the phosphoinositide 3-kinase p110 delta. Cancer Res. 2003;63(7):1667–1675. [PubMed]
  • Wang SJ, Saadi W, Lin F, Nguyen CMC, Jeon NL. Differential effects of EGF gradient profiles on MDA-MB-231 breast cancer cell chemotaxis. Exp Cell Res. 2004;300(1):180–189. doi: 10.1016/j.yexcr.2004.06.030. [PubMed] [Cross Ref]
  • Harms BD, Bassi GM, Horwitz AR, Lauffenburger DA. Directional persistence of EGF-Induced cell migration is associated with stabilization of lamellipodial protrusions. Biophys J. 2005;88(2):1479–1488. doi: 10.1529/biophysj.104.047365. [PubMed] [Cross Ref]
  • Bailly M, Condeelis JS, Segall JE. Chemoattractant-induced lamellipod extension. Microsc Res Tech. 1998;43(5):433–443. doi: 10.1002/(SICI)1097-0029(19981201)43:5<433::AID-JEMT9>3.0.CO;2-2. [PubMed] [Cross Ref]
  • Hinz B, Alt W, Johnen C, Herzog V, Kaiser H-W. Quantifying Lamella Dynamics of Cultured Cells by SACED, a new computer-assisted motion analysis. Exp Cell Res. 1999;251(1):234–243. doi: 10.1006/excr.1999.4541. [PubMed] [Cross Ref]
  • Schneider IC, Hays CK, Waterman CM. Epidermal growth factor-induced contraction regulates paxillin phosphorylation to temporally separate traction generation from de-adhesion. Mole Biol Cell. 2009;20(13):3155–3167. doi: 10.1091/mbc.E09-03-0219. [PMC free article] [PubMed] [Cross Ref]
  • Katz M, Amit I, Citri A, Shay T, Carvalho S, Lavi S, Milanezi F, Lyass L, Amariglio N, Jacob-Hirsch J. et al. A reciprocal tensin-3-cten switch mediates EGF-driven mammary cell migration. Nat Cell Biol. 2007;9(8):961–969. doi: 10.1038/ncb1622. [PubMed] [Cross Ref]
  • Xie H, Pallero MA, Gupta K, Chang P, Ware MF, Witke W, Kwiatkowski DJ, Lauffenburger DA, Murphy-Ullrich JE, Wells A. EGF receptor regulation of cell motility: EGF induces disassembly of focal adhesions independently of the motility-associated PLC gamma signaling pathway. J Cell Sci. 1998;111:615–624. [PubMed]
  • Rijken PJ, Hage WJ, Henegouwen PMPVE, Verkleij AJ, Boonstra J. Epidermal growth-factor induces rapid reorganization of the actin microfilament system in human A431 Cells. J Cell Sci. 1991;100:491–499. [PubMed]
  • Chan AY, Raft S, Bailly M, Wyckoff JB, Segall JE, Condeelis JS. EGF stimulates an increase in actin nucleation and filament number at the leading edge of the lamellipod in mammary adenocarcinoma cells. J Cell Sci. 1998;111:199–211. [PubMed]
  • Li Y, Fan J, Chen M, Li W, Woodley DT. Transforming growth factor-alpha: a major human serum factor that promotes human keratinocyte migration. J Invest Dermatol. 2006;126(9):2096–2105. doi: 10.1038/sj.jid.5700350. [PubMed] [Cross Ref]
  • Ando Y, Jensen PJ. Epidermal growth-factor and insulin-like growth Factor-I enhance keratinocyte migration. J Investig Dermatol. 1993;100(5):633–639. doi: 10.1111/1523-1747.ep12472297. [PubMed] [Cross Ref]
  • Price JT, Tiganis T, Agarwal A, Djakiew D, Thompson EW. Epidermal growth factor promotes MDA-MB-231 breast cancer cell migration through a phosphatidylinositol 3-kinase and phospholipase C-dependent mechanism. Cancer Res. 1999;59(21):5475–5478. [PubMed]
  • Maheshwari G, Wells A, Griffith LG, Lauffenburger DA. Biophysical integration of effects of epidermal growth factor and fibronectin on fibroblast migration. Biophys J. 1999;76(5):2814–2823. doi: 10.1016/S0006-3495(99)77435-7. [PubMed] [Cross Ref]
  • Maheshwari G, Wiley HS, Lauffenburger DA. Autocrine epidermal growth factor signaling stimulates directionally persistent mammary epithelial cell migration. J Cell Biol. 2001;155(7):1123–1128. doi: 10.1083/jcb.200109060. [PMC free article] [PubMed] [Cross Ref]
  • Ware MF, Wells A, Lauffenburger DA. Epidermal growth factor alters fibroblast migration speed and directional persistence reciprocally and in a matrix-dependent manner. J Cell Sci. 1998;111:2423–2432. [PubMed]
  • Monine MI, Haugh JM. Cell population-based model of dermal wound invasion with heterogeneous intracellular signaling properties. Cell Adh Migr. 2008;2(2):137–146. doi: 10.4161/cam.2.2.6511. [PMC free article] [PubMed] [Cross Ref]
  • Anderson AR, Weaver AM, Cummings PT, Quaranta V. Tumor morphology and phenotypic evolution driven by selective pressure from the microenvironment. Cell. 2006;127(5):905–915. doi: 10.1016/j.cell.2006.09.042. [PubMed] [Cross Ref]
  • Sorger PK, Niepel M, Spencer SL. Non-genetic cell-to-cell variability and the consequences for pharmacology. Curr Opin Chem Biol. 2009;13(5–6):556–561. [PMC free article] [PubMed]
  • Zaidel-Bar R, Geiger B. The switchable integrin adhesome. J Cell Sci. 2010;123(9):1385–1388. doi: 10.1242/jcs.066183. [PubMed] [Cross Ref]
  • Parsons JT, Horwitz AR, Schwartz MA. Cell adhesion: integrating cytoskeletal dynamics and cellular tension. Nat Rev Mol Cell Biol. 2010;11(9):633–643. doi: 10.1038/nrm2957. [PMC free article] [PubMed] [Cross Ref]
  • Gardel ML, Schneider IC, Aratyn-Schaus Y, Waterman CM. Mechanical integration of actin and adhesion dynamics in cell migration. Annu Rev Cell Dev Biol Vol 26. 2010;26:315–333. [PubMed]
  • Beningo KA, Dembo M, Kaverina I, Small JV, Wang Y-. Nascent focal adhesions are responsible for the generation of strong propulsive forces in migrating fibroblasts. J Cell Biol. 2001;153(4):881–888. doi: 10.1083/jcb.153.4.881. [PMC free article] [PubMed] [Cross Ref]
  • Gardel ML, Sabass B, Ji L, Danuser G, Schwarz US, Waterman CM. Traction stress in focal adhesions correlates biphasically with actin retrograde flow speed. J Cell Biol. 2008;183(6):999–1005. doi: 10.1083/jcb.200810060. [PMC free article] [PubMed] [Cross Ref]
  • Serrels B, Serrels A, Brunton VG, Holt M, McLean GW, Gray CH, Jones GE, Frame MC. Focal adhesion kinase controls actin assembly via a FERM-mediated interaction with the Arp2/3 complex. Nature Cell Biol. 2007;9(9):1046–1056. doi: 10.1038/ncb1626. [PubMed] [Cross Ref]
  • Nayal A, Webb DJ, Brown CM, Schaefer EM, Vicente-Manzanares M, Horwitz AR. Paxillin phosphorylation at Ser273 localizes a GIT1-PIX-PAK complex and regulates adhesion and protrusion dynamics. J Cell Biol. 2006;173(4):587–589. doi: 10.1083/jcb.200509075. [PMC free article] [PubMed] [Cross Ref]
  • Zaidel-Bar R, Ballestrem C, Kam Z, Geiger B. Early molecular events in the assembly of matrix adhesions at the leading edge of migrating cells. J Cell Sci. 2003;116(22):4605–4613. doi: 10.1242/jcs.00792. [PubMed] [Cross Ref]
  • DeMali KA, Barlow CA, Burridge K. Recruitment of the Arp2/3 complex to vinculin: coupling membrane protrusion to matrix adhesion. J Cell Biol. 2002;159(5):881–891. doi: 10.1083/jcb.200206043. [PMC free article] [PubMed] [Cross Ref]
  • Valles AM, Beuvin M, Boyer B. Activation of Rac1 by paxillin-Crk-DOCK180 signaling complex is antagonized by Rap1 in migrating NBT-II cells. J Biol Chem. 2004;279(43):44490–44496. doi: 10.1074/jbc.M405144200. [PubMed] [Cross Ref]
  • Geiger B, Zamir E, Katz BZ, Aota S, Yamada KM, Kam Z. Molecular diversity of cell-matrix adhesions. J Cell Sci. 1999;112(11):1655–1669. [PubMed]
  • Stricker J, Aratyn-Schaus Y, Oakes PW, Gardel ML. Spatiotemporal constraints on the force-dependent growth of focal adhesions. Biophys J. 2011;100(12):2883–2893. doi: 10.1016/j.bpj.2011.05.023. [PubMed] [Cross Ref]
  • Balaban NQ, Schwarz US, Riveline D, Goichberg P, Tzur G, Sabanay I, Mahalu D, Safran S, Bershadsky A, Addadi L. et al. Force and focal adhesion assembly: a close relationship studied using elastic micropatterned substrates. Nature Cell Biol. 2001;3(5):466–472. doi: 10.1038/35074532. [PubMed] [Cross Ref]
  • Smilenov LB, Mikhailov A, Pelham RJ, Marcantonio EE, Gundersen GG. Focal adhesion motility revealed in stationary fibroblasts. Science. 1999;286(5442):1172–1174. doi: 10.1126/science.286.5442.1172. [PubMed] [Cross Ref]
  • Berginski ME, Vitriol EA, Hahn KM, Gomez SM. High-resolution quantification of focal adhesion spatiotemporal dynamics in living cells. PLoS One. 2011;6(7):e22025. doi: 10.1371/journal.pone.0022025. [PMC free article] [PubMed] [Cross Ref]
  • Wurflinger T, Gamper I, Aach T, Sechi AS. Automated segmentation and tracking for large-scale analysis of focal adhesion dynamics. J Microsc. 2011;241(1):37–53. doi: 10.1111/j.1365-2818.2010.03404.x. [PubMed] [Cross Ref]
  • Welf ES, Ogunnaike BA, Naik UP. Quantitative statistical description of integrin clusters in adherent cells. IET Syst Biol. 2009;3(5):307–316. doi: 10.1049/iet-syb.2009.0009. [PubMed] [Cross Ref]
  • Small JV. Dicing with dogma: de-branching the lamellipodium. Trends Cell Biol. 2011;20(11):628–633. [PMC free article] [PubMed]
  • Pollard TD, Borisy GG. Cellular motility driven by assembly and disassembly of actin filaments. Cell. 2003;112(4):453–465. doi: 10.1016/S0092-8674(03)00120-X. [PubMed] [Cross Ref]
  • Dobereiner HG, Dubin-Thaler BJ, Hofman JM, Xenias HS, Sims TN, Giannone G, Dustin ML, Wiggins CH, Sheetz MP. Lateral membrane waves constitute a universal dynamic pattern of motile cells. Phys Rev Lett. 2006;97(3):4. [PubMed]
  • Machacek M, Danuser G. Morphodynamic profiling of protrusion phenotypes. Biophys J. 2006;90(4):1439–1452. doi: 10.1529/biophysj.105.070383. [PubMed] [Cross Ref]
  • Barnhart EL, Lee KC, Keren K, Mogilner A, Theriot JA. An adhesion-dependent switch between mechanisms that determine motile cell shape. Plos Biol. 2011;9(5):e1001059. doi: 10.1371/journal.pbio.1001059. [PMC free article] [PubMed] [Cross Ref]
  • Giannone G, Dubin-Thaler BJ, Dobereiner HG, Kieffer N, Bresnick AR, Sheetz MP. Periodic lamellipodial contractions correlate with rearward actin waves. Cell. 2004;116(3):431–443. doi: 10.1016/S0092-8674(04)00058-3. [PubMed] [Cross Ref]
  • Othmer HG, Dunbar SR, Alt W. Models of dispersal in biological systems. J Math Biol. 1988;26(3):263–298. doi: 10.1007/BF00277392. [PubMed] [Cross Ref]
  • Dunn GA. Characterising a kinesis response: time averaged measures of cell speed and directional persistence. Agents Actions Suppl. 1983;12:14–33. [PubMed]
  • Yuan TL, Wulf G, Burga L, Cantley LC. Cell-to-Cell variability in PI3K protein level regulates PI3K-AKT pathway activity in cell populations. Curr Biol. 2011;21(3):173–183. doi: 10.1016/j.cub.2010.12.047. [PMC free article] [PubMed] [Cross Ref]
  • Rinott R, Jaimovich A, Friedman N. Exploring transcription regulation through cell-to-cell variability. Proc Natl Acad Sci. 2011;108(15):6329–6334. doi: 10.1073/pnas.1013148108. [PubMed] [Cross Ref]
  • Wieser S, Weghuber J, Sams M, Stockinger H, Schutz GJ. Cell-to-cell variability in the diffusion constants of the plasma membrane proteins CD59 and CD147. Soft Matter. 2009;5(17):3287–3294. doi: 10.1039/b902266j. [Cross Ref]
  • Lacayo CI, Pincus Z, VanDuijn MM, Wilson CA, Fletcher DA, Gertler FB, Mogilner A, Theriot JA. Emergence of large-scale cell morphology and movement from local actin filament growth dynamics. PLoS Biol. 2007;5(9):e233. doi: 10.1371/journal.pbio.0050233. [PMC free article] [PubMed] [Cross Ref]
  • Colman-Lerner A, Gordon A, Serra E, Chin T, Resnekov O, Endy D, Gustavo Pesce C, Brent R. Regulated cell-to-cell variation in a cell-fate decision system. Nature. 2005;437(7059):699–706. doi: 10.1038/nature03998. [PubMed] [Cross Ref]
  • Goswami S, Sahai E, Wyckoff JB, Cammer N, Cox D, Pixley FJ, Stanley ER, Segall JE, Condeelis JS. Macrophages promote the invasion of breast carcinoma cells via a colony-stimulating factor-1/epidermal growth factor paracrine loop. Cancer Res. 2005;65(12):5278–5283. doi: 10.1158/0008-5472.CAN-04-1853. [PubMed] [Cross Ref]
  • JrmM G, Pittet P, Csucs G, Lussi JW, Meister J-J, Hinz B. Focal adhesion size controls tension-dependent recruitment of alpha-smooth muscle actin to stress fibers. J Cell Biol. 2006;172(2):259–268. doi: 10.1083/jcb.200506179. [PMC free article] [PubMed] [Cross Ref]
  • Diener A, Nebe B, Luthen F, Becker P, Beck U, Neumann HG, Rychly J. Control of focal adhesion dynamics by material surface characteristics. Biomaterials. 2005;26(4):383–392. doi: 10.1016/j.biomaterials.2004.02.038. [PubMed] [Cross Ref]
  • Palecek SP, Schmidt CE, Lauffenburger DA, Horwitz AF. Integrin dynamics on the tail region of migrating fibroblasts. J Cell Sci. 1996;109:941–952. [PubMed]
  • Mouneimne G, Soon L, DesMarais V, Sidani M, Song XY, Yip SC, Ghosh M, Eddy R, Backer JM, Condeelis J. Phospholipase C and cofilin are required for carcinoma cell directionality in response to EGF stimulation. J Cell Biol. 2004;166(5):697–708. doi: 10.1083/jcb.200405156. [PMC free article] [PubMed] [Cross Ref]
  • Giannone G, Dubin-Thaler BJ, Rossier O, Cai Y, Chaga O, Jiang G, Beaver W, Döbereiner H-G, Freund Y, Borisy G. et al. Lamellipodial actin mechanically links myosin activity with adhesion-site formation. Cell. 2007;128(3):561–575. doi: 10.1016/j.cell.2006.12.039. [PubMed] [Cross Ref]
  • Enculescu M, Sabouri-Ghomi M, Danuser G, Falcke M. Modeling of protrusion phenotypes driven by the actin-membrane interaction. Biophys J. 2010;98(8):1571–1581. doi: 10.1016/j.bpj.2009.12.4311. [PubMed] [Cross Ref]
  • Cirit M, Krajcovic M, Choi CK, Welf ES, Horwitz AF, Haugh JM. Stochastic model of integrin-mediated signaling and adhesion dynamics at the leading edges of migrating cells. Plos Comput Biol. 2010;6(2):e1000688. doi: 10.1371/journal.pcbi.1000688. [PMC free article] [PubMed] [Cross Ref]
Articles from BMC Biophysics are provided here courtesy of
BioMed Central