Steering cell migration by alternating blebs and actin-rich protrusions
BMC Biology volume 14, Article number: 74 (2016)
High directional persistence is often assumed to enhance the efficiency of chemotactic migration. Yet, cells in vivo usually display meandering trajectories with relatively low directional persistence, and the control and function of directional persistence during cell migration in three-dimensional environments are poorly understood.
Here, we use mesendoderm progenitors migrating during zebrafish gastrulation as a model system to investigate the control of directional persistence during migration in vivo. We show that progenitor cells alternate persistent run phases with tumble phases that result in cell reorientation. Runs are characterized by the formation of directed actin-rich protrusions and tumbles by enhanced blebbing. Increasing the proportion of actin-rich protrusions or blebs leads to longer or shorter run phases, respectively. Importantly, both reducing and increasing run phases result in larger spatial dispersion of the cells, indicative of reduced migration precision. A physical model quantitatively recapitulating the migratory behavior of mesendoderm progenitors indicates that the ratio of tumbling to run times, and thus the specific degree of directional persistence of migration, are critical for optimizing migration precision.
Together, our experiments and model provide mechanistic insight into the control of migration directionality for cells moving in three-dimensional environments that combine different protrusion types, whereby the proportion of blebs to actin-rich protrusions determines the directional persistence and precision of movement by regulating the ratio of tumbling to run times.
Efficient directed migration is assumed to rely on high directional persistence [1–3]. Indeed, in a stable chemotactic gradient, straight trajectories allow to reach the target in a minimal time. In contrast, lower directional persistence has been associated with poorly directed migration such as in the absence of chemotactic cues or in shallow chemotactic gradients [2, 3]. For instance, the persistence of fibroblasts and dendritic cells has been shown to decrease in presence of a uniform concentration of chemoattractant when compared to migration of the same cells in a chemotactic gradient . Yet, cells undergoing directed migration in vivo often display trajectories with frequent direction changes and low persistence compared to directed migration in culture [5–7]. Such trajectories have been described as biased random walks or as series of runs and tumbles, i.e., alternating phases with high and low directional persistence [8–11]. In zebrafish primordial germ cells, whose chemotactic migration during development can be described as a succession of run and tumbles, low persistence and frequent direction changes associated with tumbling have been proposed to fine-tune the migration of these cells, as they progress to intermediate targets during development [9, 12]. However, the cellular mechanisms controlling directional persistence during animal cell migration in vivo are poorly understood, and the functional importance of a proper control of this parameter remains elusive.
Here, we investigate the cellular control and function of directional persistence during cell migration in vivo. We use zebrafish early mesendoderm progenitor cells, which, during early gastrulation, predominantly migrate as single cells and display frequent direction changes . We have previously shown that mesendoderm progenitors can form different protrusion types, including blebs and actin-polymerization driven ones, and that enhancing the formation of blebs decreases migration directional persistence . Therefore, we reasoned that mesendoderm progenitors represented a good model for investigating migration directionality in vivo.
We first show, using an unbiased trajectory analysis algorithm, that lateral progenitors migrating towards the forming body axis alternate run and tumbling phases. We then employ a transplantation assay to investigate how protrusion formation relates to migration directionality during single cell migration of progenitor cells. Using custom-made cell segmentation and protrusion detection software, we show that run phases correlate with the formation of directed actin-rich protrusions, while enhanced blebbing is observed during tumbles. Changing the proportion of blebs to actin-rich protrusions changes the ratio of tumbling to run times. Strikingly, we observe that both decreasing and increasing the ratio of tumbling to run times increase cell dispersion during migration, indicative of reduced migration precision. A theoretical model quantitatively recapitulating the characteristics of progenitor cell migration indicates that an optimal tumbling-to-run ratio enhances migration precision in a changing environment. Together, our experiments and model suggest that the precision of mesendoderm progenitor cell migration depends on the ratio of tumbling to run times, and that this ratio is controlled by adjusting the proportion of blebs to actin-rich protrusions formed by these cells.
Zebrafish lateral mesendoderm progenitors display run-and-tumbling during directed migration
In order to investigate how migration directionality is determined in zebrafish mesendoderm progenitors, we transplanted mesendodermal cells (cells expressing the Nodal-ligand Cyclops (Cyc), to induce mesendoderm cell fate ) injected with a fluorescent histone in a wild type (wt) host (Fig. 1a). The transplanted cells displayed mostly single cell migration, with only sporadic interaction with neighboring mesendoderm progenitors, for at least 3 hours following transplantation (from 30 min before shield to 70 % epiboly), as previously reported . Cell nuclei were tracked for over 2 hours during mid gastrulation stages (~6–8 hours post-fertilization (hpf), starting 30 min to 1 hour post-transplantation) (Fig. 1b). We found that the trajectories of transplanted mesendoderm progenitors displayed a mean persistence, i.e., ratio of the net displacement to cell trajectory length, of 0.68 ± 0.13 (mean ± standard deviation (SD), n = 18 cells), lower than the typical persistence values observed during chemotaxis in vitro [15, 16]. An unbiased analysis of the trajectories’ cell scaled speed (S) and alignment index (a measure of the local persistence, A) revealed that the cells displayed a multi-modal behavior that can be described as alternating phases of relatively straight migration (run phases) and phases of slowed and poorly directed movement (tumble phases). Accordingly, the cell trajectories could be divided into run and tumble phases, where the cut-off between phases was determined automatically, based on a quantitative analysis of the local persistence and speed of the cells (Fig. 1c, d and Additional file 1: Supplementary Methods for details). This automated analysis yielded an average ratio of tumbling to run times in mesendodermal progenitors of 0.58 ± 0.34 (mean ± SD, n = 18 trajectories). The relatively large SD reflects the fact that both run and tumble times displayed exponential distributions, which are characterized by SDs of the order of the mean (Additional file 2: Figure S1). Instantaneous cell speed, measured with a 1.5 min time interval, was approximately 1.8 times higher during run phases compared to tumble phases (Fig. 1e). Finally, tumbles usually resulted in a significant direction change, with an average angle between successive runs of 56 ± 34 degrees (mean ± SD, n = 18 trajectories).
Even though lateral progenitors display mostly single cell migration in early gastrulation , they still transiently interact with neighboring mesendoderm progenitors, which could influence their trajectories. To investigate the migration of these cells in an in vivo environment while avoiding any influence of transient contacts with neighboring cells, we transplanted single mesendoderm cells, into the lateral side of maternal zygotic oep (MZoep) mutant embryos, which lack mesendoderm progenitors . Transplanted cells display directed migration between the yolk and the overlying ectoderm towards the dorsal side of the embryo, as their wt counterparts, but do not have neighboring cells to interact with . Thus, they represent a good model system for the study of single cell migration in a complex in vivo environment. We acquired trajectories of mesendoderm progenitors injected with a fluorescent histone transplanted into MZoep hosts and applied the same automated analysis as described above to their trajectories. We found that, similarly to progenitors transplanted into wt hosts, the cells displayed multi-modal trajectories that can be described as successions of run and tumble phases (Fig. 1f–h). Similar to progenitors migrating in wt hosts, the average ratio of tumbling to run times was 0.68 ± 0.38 (mean ± SD, n = 23 trajectories), instantaneous cell speed was approximately 1.8 times higher during run phases compared to tumble phases (Fig. 1e), and tumbles resulted in a significant direction change, with an average angle between successive runs of 68 ± 37 degrees (mean ± SD, n = 23 trajectories).
Taken together, our analysis indicates that zebrafish mesendoderm progenitors alternate phases of directed migration (runs) and reorientation events (tumbles) during directed migration in vivo.
Protrusion formation during run and tumbling phases
We have previously observed that enhancing bleb formation while reducing actin-rich protrusions in mesendoderm progenitors decreases the directional persistence of their migration . We thus asked how the formation of different protrusion types relates to the run-and-tumbling behavior of mesendoderm progenitor cells. We acquired 10–30 min high-resolution two-photon microscopy movies of transplanted mesendoderm cells injected with Alexa594-Dextran to mark the cytoplasm and expressing Lifeact-GFP  to follow filamentous actin (Fig. 2a, b and Additional file 3: Movie 1). We observed that, similarly to collectively migrating prechordal plate cells , single mesendoderm progenitors formed blebs (spherical protrusions initially devoid of actin) and actin-rich protrusions (protrusions containing actin throughout their expansion) (Fig. 2b and Additional file 3: Movie 1).
To analyze the orientation of each protrusion type with respect to the direction of cell migration, we developed a new software package for three-dimensional (3D) cell and protrusion segmentation and automated detection and identification of individual protrusions (Automated Protrusion Analyzer (APA), Fig. 2a–c and Additional file 4: Figure S2). Protrusion identification and classification is based on detection of changes in cell surface curvature and morphological differences between protrusion types. APA identifies two types of protrusions: blebs and actin-rich protrusions (Fig. 2b). Actin-rich protrusions are distinguished from blebs by the presence of actin (labeled with Lifeact) in all phases of their expansion (Additional file 3: Movie 1), and by a higher curvature than blebs (Additional file 1: Supplementary Methods). Using APA, we could monitor the center of mass of the cells and each protrusion formed, as well as the intensity of actin in actin-rich protrusions during 3D migration (Fig. 2b, c). As lamellipodia size and actin content have been shown to correlate with migration speed , we analyzed the angle distribution of actin-rich protrusions weighted with the total intensity of the Lifeact signal in the protrusion. Thus, this weighted distribution mostly reflects the orientation of larger actin-rich protrusions. The overall orientation of a specific protrusion type was quantified using the polar order parameter (POP). The POP magnitude indicates how sharply focused the protrusion angle distribution is (Additional file 1: Supplementary Methods).
We then used these automated analysis tools to relate protrusion formation to mesendoderm progenitors’ run-and-tumbling behavior. Run-and-tumbling was evident in 11 out of 17 two-photon high-resolution timelapses (Fig. 2d); in the remaining timelapses, cells displayed directed motion only, likely because the shorter (10–30 min long) high-resolution movies necessary for protrusion analysis are sometimes too short to capture the tumbling behavior. Analysis of the timelapses where run-and-tumbling could be quantified showed that, during run phases, mesendoderm cells formed actin-rich protrusions in the direction of migration (Additional file 5: Movie 2, Fig. 2d–f) and poorly oriented blebs, as evidenced by the clear difference in POP between the two protrusion types (POP = 0.444 ± 0.151 for actin-rich protrusions vs. 0.187 ± 0.197 for blebs in run phases, mean ± standard error of the mean (SEM), Fig. 2f). In contrast, tumble phases were associated with the formation of an increased number of randomly oriented blebs (Fig. 2e) and a decrease in the focus of actin-rich protrusion formation (POP = 0.158 ± 0.132 for actin-rich protrusions formed during tumble phases, mean ± SEM, Additional file 5: Movie 2, Fig. 2f). In about 15 % of the tumble events, less blebbing was observed and a change in direction was achieved by the formation of a new leading edge actin-rich protrusion (corresponding to the two cells labeled as blue data points in Fig. 2e, Additional file 6: Movie 3). Taken together, our observations suggest that actin-rich protrusions may drive directed migration of mesendoderm progenitors whereas blebs mainly contribute to cell re-orientation.
Modulating the proportion of blebs to actin-rich protrusions changes the ratio of tumbling to run times without affecting protrusion orientation
To test whether the proportion of blebs to actin-rich protrusions formed by mesendoderm progenitors determines their run-and-tumbling behavior, we aimed to change the frequency of bleb formation. We increased bleb formation by reducing membrane-to-cortex attachment using a morpholino (MO) against ezrin , a protein that binds the actin cortex to the plasma membrane. Consistent with our previous observations in the prechordal plate , we found that single transplanted mesendoderm cells with reduced Ezrin activity showed a strong increase in the frequency and size of blebs and a reduction in actin-rich protrusions (Fig. 3a–c, Additional file 7: Figure S3A and Additional file 8: Movie 4). We previously showed that enhancing bleb formation by reducing Ezrin activity (either by expressing a dominant negative version of Ezrin or using a MO against ezrin) significantly reduces migration directional persistence, leading to less straight cell migration tracks in transplanted mesendoderm cells . We thus asked whether the decrease in directional persistence in ezrin-MO cells was due to increased tumbling. Alternatively, reduced directional persistence could result from a change in the focus of protrusion expansion, as Ezrin depletion affects the entire cell and could affect overall cell polarity. To distinguish between these two possibilities, we analyzed protrusion orientation in ezrin morphant cells. We observed that the angle distributions of blebs and actin-rich protrusions were not affected by Ezrin depletion (Fig. 3d and Additional file 7: Figure S3B–D). We then analyzed the trajectories of transplanted progenitor cells during mid gastrulation stages (6–8 hpf) for control cells and ezrin morphant cells. We found that enhanced bleb formation in ezrin morphant mesendoderm progenitors significantly increased the ratio of the time spent tumbling to the time spent in run phases (Fig. 3e). This increase was due to a decrease in the duration of run phases (on average 5 min in control runs, n = 209, vs. 3.8 min in ezrin-MO runs, n = 231), while the duration of individual tumble phases was not significantly changed (on average 3.1 min in control tumbles, n = 216, vs. 3 min in ezrin-MO tumbles, n = 234).
We next sought to investigate how increasing the formation of actin-rich protrusions at the expense of blebs affects the run-and-tumbling behavior of mesendoderm progenitors. To this end, we increased membrane-to-cortex attachment by expressing a constitutively active version of Ezrin (CAEzrin, T564D ). CAEzrin-expressing transplanted single mesendoderm cells showed a strong decrease in blebbing activity and an increase in formation of actin-rich protrusions (Fig. 4a–d and Additional file 9: Movie 5). We then investigated how expression of CAEzrin affected the migratory trajectories of single mesendoderm progenitors transplanted into MZoep hosts from mid-to-late gastrulation stages (6–8 hpf). We observed that single CAEzrin expressing mesendoderm progenitors showed an increase in migration directional persistence and net speed, while their instantaneous speed remained unchanged compared to co-transplanted control cells (Fig. 4e, f). We first checked whether this increase in directional persistence could result from an overall increase in the focus of protrusion formation upon expression of CAEzrin. We found that the angle distribution of actin-rich protrusion formation was less focused in CAEzrin-expressing cells than in control cells, indicating that the observed increase in cell directional persistence does not result from more focused actin-rich protrusions (Fig. 4g, Additional file 7: Figure S3C, D and Additional file 9: Movie 5). Bleb formation was rarely observed and only a few events could be analyzed (Fig. 4c and Additional file 7: Figure S3B). We then investigated whether expression of CAEzrin affected the run-and-tumbling behavior of mesendoderm progenitors, and found that the ratio of tumbling to run times was decreased in progenitors expressing CAEzrin (Fig. 4h). This decrease was due to an increase in the duration of run phases (on average 5 min in control runs, n = 209, vs. 6.4 min in CAEzrin runs, n = 102), while the duration of individual tumble phases was not significantly affected (on average 3.1 min in control tumbles, n = 216, vs. 3 min in CAEzrin tumbles, n = 104). Together, these observations suggest that the proportion of blebs to actin-rich protrusions controls the directional persistence of cell migration in mesendoderm progenitors by modulating the ratio of tumbling to run times.
Modulating the ratio of tumbling to run times affects migration precision
Frequent direction changes have been proposed to enhance the precision of cell migration in complex environments, particularly during directed migration where the chemotactic target is moving or changing over time as might be the case during zebrafish gastrulation [9, 21]. Indeed, considering that mesendoderm cells migrate dorsally and vegetally towards the forming body axis, it is commonly believed that they follow a chemotactic signal from the epiboly front. We thus asked whether changing directional persistence affects the overall precision of mesendoderm progenitor migration. We assessed the precision of cell migration by quantifying the spatial dispersion after approximately 2 hours of migration of cells that were co-transplanted at the same location at 50 % epiboly, for cells with different levels of Ezrin activity. Interestingly, we found that both the cells displaying enhanced blebbing and tumbling, and the cells displaying enhanced formation of actin-rich protrusions and running, had a significantly higher spatial dispersion than control cells (Fig. 5a). These observations suggest that both decreasing and increasing the ratio of tumbling to run times in mesendoderm progenitors decreases the precision of cell migration.
To test whether the ratio of tumbling to run times observed in mesendoderm progenitors might indeed optimize migration precision, we developed a stochastic model of cells migrating towards a target moving at constant speed. We represented the moving cells by active Brownian particles randomly switching between run and tumble phases (Fig. 5b, Additional file 1: Supplementary Methods, Additional file 10: Figure S4 and Additional file 11: Figure S5). During run phases cells perform directed active Brownian motion with stochastic speed and a direction fluctuating around a mean value oriented towards the target with a detection error. During tumble phases cells are randomly moving without any preferred direction. We constrained the model parameters by comparing characteristic observables of motion obtained from simulated tracks (analyzed with the same procedure as applied to the experimental data) to experimental measurements. Specifically, several parameters describing cell velocity, as well as run and tumble durations were compared between simulations and experiments. A parameter search yielded a set of parameters very accurately accounting for measured experimental values in control mesendodermal cells (Additional file 1: Table S2 and Additional file 1: Supplementary Methods for details). We found that, with this selected set of parameters, the combined 2D distribution of alignment and cell speed, and the probability distribution of cell speeds during run and tumble phases were well captured by the simulations without further fitting (Fig. 5c, compare to Fig. 1g, and Fig. 5d, e). These observations indicate that the numerical model accurately captures the aspects of cell migration relevant to the observed progenitor trajectories.
Using the estimated parameters, we then systematically varied the run time of the model cells and assessed the precision of cell migration by computing the distance to target and the dispersion of the cell population at the end of the experiment (te = 1.5 h). We found that the distance to target decreased as a function of the run time, indicating that longer runs are more favorable for overall cell velocity. Strikingly, cell dispersion showed a clear minimum around the mean run time measured for control mesendoderm progenitors. This prediction is consistent with the increased cell dispersion measured for CAEzrin and ezrin-MO cells (Fig. 5a), which display run times longer and shorter than control cells, respectively. Taken together, our experiments and model thus indicate that the ratio of tumbling to run times is a critical factor controlling the precision of cell migration in vivo.
Low directional persistence is often thought to be a consequence of a shallow chemotactic gradient resulting in the formation of unfocused protrusions [1, 3]. Here, we show that the directional persistence of zebrafish mesendoderm progenitors migrating in vivo does not depend on the directional focus of protrusion formation, but rather is determined by the ratio of persistent run phases to tumble phases associated with cell reorientation. Interestingly, progenitor cells appear to control the ratio of tumbling to run times by adjusting the proportion of blebs to actin-rich protrusions formed during migration. Blebs have previously been implicated in mediating directed migration of primordial germ cells during zebrafish embryogenesis , and of a number of cancer lines in culture and in vivo [23, 24]. In zebrafish primordial germ cells, bleb growth appears to expand the cell body forward, and subsequent anchoring of the bleb neck to the substrate by adhesive contacts to surrounding cells is thought to drive cell migration . Our finding that blebs in mesendoderm progenitor cells are predominantly associated with tumbling reorientation events suggests that, in these cells, blebs are primarily used for exploring the environment, whereas actin-rich protrusions drive directed migration during run phases. Specifically, undirected bleb formation, as observed during tumble phases, induces displacement of the cell towards random directions and might thus provide a stochastic way of exploring the environment. This difference in bleb function between primordial germ cells and mesendoderm cells may be due to the fact that mesendoderm progenitors form directed actin-rich protrusions, whereas primordial germ cell migration appears to rely exclusively on blebs .
The run and tumbling behavior of control mesendoderm progenitors appears highly similar for cells in wt and in MZoep hosts. Furthermore, our experiments indicate that the ratio of run and tumbling can be modulated in single transplanted cells by tuning the amount of Ezrin activity. To account for experimental variability between different embryos, cells with increased or decreased Ezrin activity were always co-transplanted with control cells in the same MZoep embryo (internal controls) (see also ). These observations indicate that run and tumbling is largely a cell autonomous behavior. Nonetheless, it remains to be investigated whether extracellular factors, such as the distribution, organization and nature of extracellular matrix or the proximity to the chemotactic signal followed by the cells, influence run and/or tumbling in zebrafish mesendoderm progenitors.
Run-and-tumbling is a common feature of bacterial chemotaxis, where it is a strategy for efficient gradient sensing , but has also been observed in a variety of eukaryotic motile cells, including primordial germ cells , chlamydomonas , and mammary epithelial cells . Bacteria are too small to accurately measure a chemoattractant gradient without moving, and use temporal comparisons instead, leading to a biased random walk with longer run phases in the direction of the chemotactic gradient. Animal cells are large enough to polarize in a gradient without motion  and thus alternating run and tumbling phases during migration is likely to serve a different function than in bacterial chemotaxis. It has been speculated that tumble-associated direction changes might increase the precision of chemotactic cell migration in animal cells [12, 21]. Our observation that changing the ratio of tumbling to run times impairs the focus of cell migration provides direct experimental evidence supporting this hypothesis. Indeed, both increasing and decreasing the tumbling to run ratio by modulating the bleb-to-actin-rich protrusion ratio led to impaired cell migration precision (Fig. 5a). Distinct molecular pathways regulate the formation of blebs and actin-rich protrusions [23, 30], suggesting that the ratio between the two protrusion types could be readily tuned. Such a sub-specialization of protrusion function would allow cells to easily modulate the frequency of re-orientation events during migration in complex and changing environments. Our theoretical model, which recapitulated key features of mesendoderm progenitor migration, predicts that an optimal tumble to run ratio enhances migration precision. Indeed, too long runs increase cell dispersion by overly amplifying initial errors in migration direction, whereas too short runs increases cell dispersion because frequent direction changes enhance heterogeneity in direction between cells. Furthermore, it is possible that alternating run and tumbles enhances the robustness of migration to noise in, for example, lamellipodium orientation .
Our experiments and model indicate that mesendoderm progenitors may be operating close to an optimum tumbling to run ratio for precise migration in the in vivo context of the developing zebrafish embryo. Taken together, our data suggest that, by adjusting the proportion of blebs to actin-rich protrusions, mesendoderm cells modulate the ratio of run to tumbling times, and thereby control the precision of their migration. A number of cell types have been reported to combine blebs and actin-rich protrusions during migration [32–35]. Future studies will have to investigate whether blebs and actin-rich protrusions also have distinct functions in these cells types.
Embryo staging and maintenance
MRNA, morpholino, and dye injection
mRNA was synthesized as previously described . For single cell transplantation, wt TL embryos were injected with 50 pg of Lifeact-GFP , 3.25 ng of Dextran Alexa Fluor-595 (D22913, Invitrogen), and 100 pg of cyc alone (control) or together with 4 ng of ezrin-UTR-MO , to generate ezrin-MO cells or 150 pg of CAEzrin mRNA (T564D of Danio rerio’s gene as in ) to generate CAEzrin cells.
For tracking of cell nuclei in low magnification transplantation experiments, wt donor embryos were injected with 100 pg of cyc together with Alexa Fluor-488 conjugated histone H1 (H13188, Invitrogen) (control), or 100 pg of histoneH2Azf::mcherry plus 150 pg of CAEzrin mRNA (CAEzrin cells). MZoep host embryos were injected with Dextran Alexa Fluor-647 (D22914, Invitrogen) (see also ).
Transplantation experiments, cell imaging, and bleb size measurements
For transplantation experiments, wt and experimental TL donors and MZoep dharma::GFP host embryos were dechorionated with Pronase (2 mg/mL in E2) and transferred onto an agarose plate with E3 medium. Two to three cells were taken from control and experimental donor embryos at dome stage (4.5 hpf) and co-transplanted into the emerging lateral mesendoderm of a host embryo labeled with Dextran Alexa Fluor-647 at 50 % epiboly (5.5 hpf).
For low magnification experiments, time-lapse images were obtained with an upright Leica SP5 confocal microscope equipped with a 20× water immersion lens, using 488-nm Argon, DPSS 561 nm, and 633-nm HeNe laser lines. Frames were captured at 90 s intervals for 3 h (~5.5–8.5 hpf). The temperature was kept constant in all videos (28 °C).
For big magnification transplantation experiments, images were obtained with a Zeiss 710 two-photon microscope equipped with a 63×/1.2 objective, using 910 nm wavelength of the Chamaleon laser. Frames were captured at 10–25 s intervals for 10–30 min, between 6 and 8 hpf.
For bleb size measurements, the projected area of each bleb at its maximal extension was measured using ImageJ and normalized to the projected area of the whole cell.
For cell dispersion measurements, pictures were taken with a dissecting microscope (Olympus SZX 12) equipped with a QImaging Micropublisher 5.0 camera approximately 2 hours post-transplantation.
For single transplanted cells in low magnification movies nuclei tracking in three dimensions (x, y, and z) was performed with Imaris 7.3.0 software. The instantaneous and net speeds, as well as directional persistence (ratio of the net displacement to the distance actually traveled by the cells), were extracted from the tracks.
Analysis of the directions of protrusion formation in combination with cell tracking in higher magnification movies was performed using the APA software, described in Additional file 1: Supplementary Methods.
t tests were performed after the data were confirmed to have normal distribution and equal variance; otherwise, Mann–Whitney U tests were applied. P values were computed in R. For low magnification cell transplantation experiments and variance of cell position (used to assess cell dispersion), one-sided t test was used, which compared experimental data points to an equal sized group of 1. We also computed the P values with ttest2 from Matlab, which compared experimental data points with a random distribution of numbers around one with the same standard deviation as our data. ttest2 yielded similar results and conclusions.
To numerically describe the angular distribution of protrusions, we used the polar order parameter (POP), as explained in detail in Additional file 1: Supplementary Methods. We consider two POP values to be significantly different when their SEMs do not overlap.
Definition of run-and-tumbling phases
For longer trajectories (Figs. 1 and 4e, f), a timeframe of 1.5 min was used as it maximized the amount of embryos we could image simultaneously without a change in the run-to-tumble behavior or in the instantaneous speed. Run-and-tumbling phases were automatically extracted using and unbiased procedure described in Additional file 1: Supplementary Methods [39, 40]. For the analysis of short cell trajectories (timeframe ~10 s, Figs. 2 d–f, 3e, 4 h), “runs” were defined as phases where the trajectory does not deviate by more than 45 degrees from the direction at the beginning of the run or if a change in direction larger than 45 degrees persists for less than 5 timeframes. “Tumbles” were defined as phases where a change of direction higher than 45 degrees occurs and persists for longer than 5 timeframes.
Measurements of cell dispersion
Cell dispersion was assessed using cell position variance, as measured by adding the variances in x and y of the positions of control and experimental cells approximately 2 hours after they were co-transplanted at the same location in a host embryo at 50 % epiboly. Only embryos with at least three control cells and three experimental cells were considered. The ratio has been normalized to transplanted control cells in the same embryo (internal controls) to account for experimental variability between individual transplantation experiments.
Alignment index (a measure of the local persistence)
- MZoep :
Maternal zygotic oep
Green fluorescent protein
Automated protrusion analyzer
Polar order parameter
Standard error of the mean
- CAEzrin :
Constitutively active version of Ezrin
King JS, Insall RH. Chemotaxis: finding the way forward with Dictyostelium. Trends Cell Biol. 2009;19:523–30.
Theisen U, Straube E, Straube A. Directional persistence of migrating cells requires Kif1C-mediated stabilization of trailing adhesions. Dev Cell. 2012;23:1153–66.
Petrie RJ, Doyle AD, Yamada KM. Random versus directionally persistent cell migration. Nat Rev Mol Cell Biol. 2009;10:538–49.
Arrieumerlou C, Meyer T. A local coupling model and compass parameter for eukaryotic chemotaxis. Dev Cell. 2005;8:215–27.
Arboleda-Estudillo Y, Krieg M, Stühmer J, Licata NA, Muller DJ, Heisenberg C-P. Movement directionality in collective migration of germ layer progenitors. Curr Biol. 2010;20:161–9.
Woo S, Housley MP, Weiner OD, Stainier DYR. Nodal signaling regulates endodermal cell motility and actin dynamics via Rac1 and Prex1. J Cell Biol. 2012;198:941–52.
Liu L, Das S, Losert W, Parent CA. mTORC2 regulates neutrophil chemotaxis in a cAMP- and RhoA-dependent fashion. Dev Cell. 2010;19:845–57.
Miller MJ, Wei SH, Parker I, Cahalan MD. Two-photon imaging of lymphocyte motility and antigen response in intact lymph node. Science. 2002;296:1869–73.
Reichman-Fried M, Minina S, Raz E. Autonomous modes of behavior in primordial germ cell migration. Dev Cell. 2004;6:589–96.
Textor J, Peixoto A, Henrickson SE, Sinn M, von Andrian UH, Westermann J. Defining the quantitative limits of intravital two-photon lymphocyte tracking. Proc Natl Acad Sci U S A. 2011;108:12401–6.
Sarris M, Masson J-B, Maurin D, Van der Aa LM, Boudinot P, Lortat-Jacob H, et al. Inflammatory chemokines direct and restrict leukocyte migration within live tissues as glycan-bound gradients. Curr Biol. 2012;22:2375–82.
Minina S, Reichman-Fried M, Raz E. Control of receptor internalization, signaling level, and precise arrival at the target in guided cell migration. Curr Biol. 2007;17:1164–72.
Diz-Muñoz A, Krieg M, Bergert M, Ibarlucea-Benitez I, Muller DJ, Paluch E, et al. Control of directed cell migration in vivo by membrane-to-cortex attachment. Plos Biol. 2010;8:e1000544.
Link V, Carvalho L, Castanon I, Stockinger P, Shevchenko A, Heisenberg C-P. Identification of regulators of germ layer morphogenesis using proteomics in zebrafish. J Cell Sci. 2006;119:2073–83.
White DP, Caswell PT, Norman JC. alpha v beta3 and alpha5beta1 integrin recycling pathways dictate downstream Rho kinase signaling to regulate persistent cell migration. J Cell Biol. 2007;177:515–25.
Xu J, Wang F, Van Keymeulen A, Rentel M, Bourne H. Neutrophil microtubules suppress polarity and enhance directional migration. Proc Natl Acad Sci U S A. 2005;102:6884.
Gritsman K, Zhang J, Cheng S, Heckscher E, Talbot WS, Schier AF. The EGF-CFC protein one-eyed pinhead is essential for nodal signaling. Cell. 1999;97:121–32.
Riedl J, Crevenna A, Kessenbrock K, Yu J, Neukirchen D, Bista M, et al. Lifeact: a versatile marker to visualize F-actin. Nat Methods. 2008;5:605–7.
Keren K, Pincus Z, Allen G, Barnhart E, Marriott G, Mogilner A, et al. Mechanism of shape determination in motile cells. Nature. 2008;453:475–80.
Gautreau A, Louvard D, Arpin M. Morphogenic effects of ezrin require a phosphorylation-induced transition from oligomers to monomers at the plasma membrane. J Cell Biol. 2000;150:193–203.
Erban R, Othmer HG. Taxis equations for amoeboid cells. J Math Biol. 2007;54:847.
Blaser H, Eisenbeiss S, Neumann M, Reichman-Fried M, Thisse B, Thisse C, et al. Transition from non-motile behaviour to directed migration during early PGC development in zebrafish. J Cell Sci. 2005;118:4027–38.
Bergert M, Chandradoss SD, Desai RA, Paluch E. Cell mechanics control rapid transitions between blebs and lamellipodia during migration. Proc Natl Acad Sci U S A. 2012;109(36):14434–9.
Tozluoğlu M, Tournier AL, Jenkins RP, Hooper S, Bates PA, Sahai E. Matrix geometry determines optimal cancer cell migration strategy and modulates response to interventions. Nat Cell Biol. 2013;15:751–62.
Kardash E, Reichman-Fried M, Maître J-L, Boldajipour B, Papusheva E, Messerschmidt E-M, et al. A role for Rho GTPases and cell-cell adhesion in single-cell motility in vivo. Nat Cell Biol. 2010;12:47–53.
Sourjik V, Wingreen NS. Responding to chemical gradients: bacterial chemotaxis. Curr Opin Cell Biol. 2012;24:262–8.
Polin M, Tuval I, Drescher K, Gollub JP, Goldstein RE. Chlamydomonas swims with two “gears” in a eukaryotic version of run-and-tumble locomotion. Science. 2009;325:487–90.
Potdar AA, Lu J, Jeon J, Weaver AM, Cummings PT. Bimodal analysis of mammary epithelial cell migration in two dimensions. Ann Biomed Eng. 2009;37:230–45.
Servant G, Weiner OD, Herzmark P, Balla T, Sedat JW, Bourne HR. Polarization of chemoattractant receptor signaling during neutrophil chemotaxis. Science. 2000;287:1037–40.
Ridley AJ. Life at the leading edge. Cell. 2011;145:1012–22.
Romanczuk P, Salbreux G. Optimal chemotaxis in intermittent migration of animal cells. Phys Rev E. 2015;91:042720–10.
Trinkaus J. Surface activity and locomotion of Fundulus deep cells during blastula and gastrula stages. Dev Biol. 1973;30:69–103.
Tyson RA, Zatulovskiy E, Kay RR, Bretschneider T. How blebs and pseudopods cooperate during chemotaxis. Proc Natl Acad Sci U S A. 2014;111:11703–8.
Imamura F, Mukai M, Ayaki M, Takemura K, Horai T, Shinkai K, et al. Involvement of small GTPases Rho and Rac in the invasion of rat ascites hepatoma cells. Clin Exp Metastasis. 1999;17:141–8.
Grinnell F. Migration of human neutrophils in hydrated collagen lattices. J Cell Sci. 1982;58:95–108.
Westerfield M. The zebrafish book: a guide for the laboratory use of zebrafish (Danio rerio). 4th ed. Eugene: University of Oregon Press; 2000.
Kimmel C, Ballard W, Kimmel S, Ullmann B, Schilling T. Stages of embryonic development of the zebrafish. Dev Dyn. 1995;203:253–310.
Montero J, Carvalho L, Wilsch-Brauninger M, Kilian B, Mustafa C, Heisenberg C-P. Shield formation at the onset of zebrafish gastrulation. Development. 2005;132:1187–98.
Blackman RB, Tukey JW. The measurement of power spectra from the point of view of communications engineering—Part II. Bell System Technical J. 1958;37:185–282.
Scott DW. Multivariate Density Estimation. New York: Wiley; 2015.
Gardiner C. Stochastic Methods. Berlin: Springer; 2009.
Romanczuk P, Bär M, Ebeling W, Lindner B, Schimansky-Geier L. Active Brownian particles. Eur Phys J Special Topics. 2012;202(1):1–162.
Romanczuk P, Schimansky-Geier L. Brownian motion with active fluctuations. Phys Rev Lett. 2011;106:230601.
van den Broeck C. On the relation between white shot noise, Gaussian white noise, and the dichotomic Markov process. J Stat Phys. 1983;31:467–83.
Hansen N. The CMA Evolution Strategy: A Comparing Review. In: Lozano JA, Larraaga P, Inza I, Bengoetxea E, editors. Towards a New Evolutionary Computation. Number 192 in Studies in Fuzziness and Soft Computing. Heidelberg: Springer; 2006. p. 75–102.
We thank K. Lee, C. Norden, A. Webb, and the members of the Paluch lab for comments on the manuscript. We are grateful to P. Rørth and Peter Dieterich for discussions, S. Ares, Y. Arboleda-Estudillo and S. Schneider for technical help, M. Biro for help with programming, and the BIOTEC/MPI-CBG and IST zebrafish and imaging facilities for help and advice at various stages of this project.
This work was supported by the Max Planck Society, the Medical Research Council UK (core funding to the MRC LMCB), and by grants from the Polish Ministry of Science and Higher Education (454/N-MPG/2009/0) to EKP, the Deutsche Forschungsgemeinschaft (HE 3231/6-1 and PA 1590/1-1) to CPH and EKP, a A*Star JCO career development award (12302FG010) to WY and a Damon Runyon fellowship award to ADM (DRG 2157-12). This work was also supported by the Francis Crick Institute which receives its core funding from Cancer Research UK (FC001317), the UK Medical Research Council (FC001317), and the Wellcome Trust (FC001317) to GS.
ADM, PR, CPH, GS, and EP designed the research and wrote the manuscript. ADM performed most of the experiments and the analysis of protrusions; MB helped with analysis and KY helped with experiments. PR performed the run and tumble data analysis and simulations. PR and GS developed the theoretical model. WY developed APA. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Supplementary Methods. Supplementary methods file with detailed description of: 1. Data analysis: (A) Automatic Protrusion Analyzer (APA) software. (B) Polar Order Parameter (POP) used to characterize the distribution of orientations of protrusions of the model. (C) Automatic detection of run and tumbles. 2. Model of cell migration: (A) Cell migration model description [41–44]. (B) Parameterization of the computational model using experimental measurements . (C) Model Predictions – Distance to target and position variance . (PDF 203 kb)
Durations of run and tumble phases. Distribution of run and tumble durations as detected by our algorithm from low-magnification movies of wt cells transplanted in a MZoep host. Frames were captured at 90 s intervals for 3 h (~5.5–8.5 hpf). The inset shows the same data in a semi-logarithmic plot, where a line indicates an exponential relationship. Exponentially distributed times are expected in a run and tumble trajectory if we assume a simple two-state Markov processes, namely a process that undergoes transitions from one state to another where the probability distribution of the next state depends only on the current state and not on the sequence of events that preceded it . (PDF 103 kb)
Accuracy of protrusion detection using the Automatic Protrusion Analyzer (APA). APA accurately detects 83 % of the blebs formed by the cells as validated by two independent experimentalists. The majority of false positive detected blebs consisted of the same bleb being detected at more than one time point. Number of analyzed cells = 28. Manual control was also performed for actin-rich protrusions and detection was found to be accurate for more than 99 % of the protrusions and for this reason we do not provide the results of the manual segmentation here. (PDF 148 kb)
Frequency of actin-rich protrusions, orientation of cell protrusions, and summary of POP values. (A) Frequency of actin-rich protrusions in control and ezrin-MO-injected mesendoderm cells. Arbitrary units (AU) are used as actin-rich protrusions are weighted with the intensity of the Lifeact signal in the protrusion. (B) Orientation of blebs with respect to the local migration axis. POP: Mean ± SEM of the magnitude of the polar order parameter. Number of analyzed cells = 17 for wt, 6 for ezrin-MO, and 6 for CAEzrin. Number of blebs n = 349 for wt, 163 for ezrin-MO, and 6 for CAEzrin. Statistical significance was determined comparing the mean ± SEM of the magnitude of the POP of the angular distributions. (C) Orientation of bleb and actin-rich protrusion formation with respect to the Yolk-ectoderm axis. For all experimental conditions, protrusions are almost exclusively oriented perpendicular to the Yolk-ectoderm axis, indicating that protrusions are formed into the extracellular space between the Yolk cell and the overlaying ectoderm layer. (D) Mean values ± SEM of the magnitude of the POP of the angular distributions of the analyzed protrusions. Green shaded area covers the wt mean ± SEM. Number of analyzed cells = 17 for wt, 6 for ezrin-MO, and 6 for CAEzrin. Number of actin-rich protrusions = 10853 for wt, 1501 for ezrin-MO, 1160 and 2549 for CAEzrin. Number of blebs n = 349 for wt and 163 for ezrin-MO. For CAEzrin only 6 blebs were observed so the POP was not calculated. (PDF 600 kb)
Computational model description and supplementary results. (A) Schematic visualization of the simulation set-up: three-dimensional model where the cell is confined in the z-direction to the mesoendodermal layer of height ΔME. The red arrow indicates the direction of polarization within the layer. (B) Schematic visualization of example model trajectories in the mesoendodermal plane (xy) in the tumble and run phases. The arrow in the run panel indicates the direction of the run determined by the polar angle 'φD = 0, which sets the average direction of migration during the run. The dashed lines indicate the instantaneous movement angle 'φ(t0) at time t0. Impact of orientation error ε on the model results (C–E). Distance to target and spatial variance of the cells at t = 90 min as a function of the scaled avg. run time τr for different values of the reorientation error towards the target ε (value estimated for wt cells and used in the main text ε = 0:2, Fig. 5); (C) ε = 0:1, (D) ε = 0:3, and (E) ε = 0:4. Average relative error versus run time (F). Average relative error E rel defined as the mean observable deviation (Additional file 1: Table S2) between model simulations and experimental results (wt cells in MZoep host). All model parameters as fitted to experimental data with the average run time being varied around the fitted value as in Fig. 5b. The vertical line represents the threshold value 0.1, which defines the shaded ‘consistency region’ in Fig. 5b. (PDF 1475 kb)
Comparison of experimental and simulated cell trajectories (A). Examples of cell trajectories from wild-type cells in an MZoep host from low-magnification experiments (top) and from model simulations with the fitted parameters (bottom). All experimental (three-dimensional) trajectories were projected to the primary two-dimensional plane of their motion and rotated so that their average direction of migration aligns with the y-axis as in the simulations. Frames were captured at 90 s intervals for 3 h (~5.5–8.5 hpf). The red symbols indicate the points along the trajectory identified as tumbling events using the algorithm introduced in Sec. I C. Simulated cell trajectories corresponding to different experimental conditions (B). Examples of simulated cell trajectories for different values of model parameter τr with other parameters as fitted for wt cells in an MZoep embryo. (left) τr as fitted to the wt cells (τr = 8 min). (center) For long runs (τr = 25 min), which yield T r = T t - values consistent with results obtained from CAEzrin cells, and (right) short runs (τr = 4 min), which yield T r = T t , values consistent with the ones measured for ezrin-MO cells. (PDF 1651 kb)
About this article
Cite this article
Diz-Muñoz, A., Romanczuk, P., Yu, W. et al. Steering cell migration by alternating blebs and actin-rich protrusions. BMC Biol 14, 74 (2016). https://doi.org/10.1186/s12915-016-0294-x
- Cell migration
- Protrusion orientation
- Actin-rich protrusion
- Run and tumble