Pull-Push Action Helps Neural Crest Cells Migrate Collectively
Neural crest cells (NCCs) are produced at the neural plate during early stages of vertebrate
embryogenesis. Then they migrate in different streams to various parts of the body that are taking shape,
the cephalic stream to the head to form various facial features, and the enteric and cardiac streams toward
the caudal end, and then from the dorsal to the ventral side to form the various organs (Fig. 1). This
migration is rapid and precisely controlled by the various chemical signals at different locations, some
attracting the NCCs while others repelling them. Thus, NCC migration is a prime example of chemotaxis [1, 2].
Figure 1: A schematic showing the migration of neural crest cells (NCCs) during embryonic
development. Adapted from Ref.  with permission, ©2013, The Company of Biologists.
What is of interest to us in this article, however, is not the chemotaxis but the peculiar behavior of
spontaneous collective migration, i.e. collective migration of a NCC cluster in the absence of external
chemical cues. This has been demonstrated mostly by in vitro experiments. Mayor and coworkers [4, 5]
have shown that an NCC cluster spontaneously migrates down a corridor coated with the adhesive
molecule fibronectin, over the course of several hours. In vivo evidence also points to spontaneous
collective migration at work [6, 7], although the evidence tends to be more circumstantial as it is difficult
to ensure the absolute absence of an external chemoattractant. Spontaneous collective migration has
been documented for other cell types in conﬁned geometries, e.g. bovine capillary endothelial
cells  and Madin-Darby canine kidney cells . The question is: how do cells interact
within a cluster to maintain persistent directional migration over long distances and long time
Two critical factors in spontaneous collective migration have been identified. The ﬁrst is contact
inhibition of locomotion (CIL), the tendency of NCCs to move away from each other upon contact, a
process mediated by N-cadherins and the non-canonical Wnt signalling pathway [10, 11]. Second,
co-attraction (COA) describes how NCCs attract each other through the autocrine production of the
short-ranged chemoattractant C3a and its receptor C3aR . Computer models have successfully
reproduced spontaneous migration of a group of NCCs under the simultaneous action of CIL and COA [12, 13]. However, we ask how CIL and COA, neither having any inherent directional preference, could
produce persistent directional migration. What keeps the cell cluster moving down the corridor, instead of
reversing course? We started out by building our own computer model that encodes CIL and COA, and it
failed to recapitulate persistent collective migration. Once away from the end of the corridor,
the group loses direction; its centroid wanders back and forth in essentially a 1-D Brownian
From these observations we came up with the hypothesis that a symmetry-breaking mechanism has to
be at work, in addition to CIL and COA. This we have called the persistent of polarity (POP). In recent
work published in Developmental Biology , we have tested this hypothesis on the basis of the
biochemical signaling pathways, as opposed to postulating rules governing cell interaction as commonly
done in agent-based or rule-based modeling [12, 13]. This online article summarizes the main findings of
. First, CIL and COA arise naturally from the underlying reaction and diffusion of Rho GTPases.
Second and more importantly, we demonstrate that CIL and COA, acting on a randomization scheme for
cell polarity, produce persistence of cell polarity, and consequently spontaneous collective
migration of NCCs. As it turns out, POP is not an additional rule to be posed alongside
CIL and COA. It is in fact a natural outcome of CIL and COA, as well as a conduit through
which these two fundamental mechanisms give rise to the observed spontaneous collective
Our general conceptualization of the NCC collective migration is as follows. Polarization and protrusion of
individual cells are governed by Rho GTPases on the membrane , subject to turnover
between the membrane-bound and cytoplasmic forms of the signaling proteins. In particular,
Rac1 and RhoA have been identiﬁed as the key proteins modulating CIL and COA . Rac1
promotes F-actin assembly in lamellipodia at the protrusive front of the cell, while RhoA
enhances myosin-induced cell contraction at the rear . For simplicity, our model only
accounts for Rac1 and RhoA and omits other GTPases such as Cdc42 and various downstream
regulators. The level of active Rac1 and RhoA on the membrane determines the protrusive and
contractile forces on the membrane and in turn the deformation and movement of the cell. Cell-cell
interaction occurs through modulating each other’s Rac-Rho biochemistry. For example, cell-cell
contact inhibits Rac1 and elevates RhoA at the site of contact in both cells. Thus, the cell
protrusions retract and the cells move apart. Finally, the Rac-Rho dynamics is subject to a
random noise so as to produce the tortuous trajectory commonly seen for single migrating
Figure 2: A model cell represented by a polygon of elastic edges.
The vertices bear active and inactive levels of Rho GTPases and move according to the various forces acting on them. The edge between vertices
and carries an elastic tension of .
The normal forces include pressure
and a protrusion or contraction force .
Each cell is represented by a polygon of
vertices connected by elastic edges (Fig. 2), similar to the vertex models widely used
for epithelial morphogenesis [17, 18]. Rac1 and RhoA levels are deﬁned on the
vertices as well as in the cytoplasm. Through “mesh refinement”, we have chosen
simulations in the rest of the paper . Each vertex on the cell membrane is subject to 2D forces acting
in the plane: a pressure force from the cytoplasm enclosed by the cell membrane, cortical tension in the
membrane segments, and a protrusion or contraction force determined by the Rac1 and RhoA levels on
the vertex. As a result, the vertices move, without inertia, at a speed determined by the resultant force
and a friction factor. We do not explicitly account for focal adhesions between the cells and a
The biochemical and mechanical components of the model are intimately coupled. For the biochemical
component, the model tracks three forms of the Rac1 and RhoA proteins: the active membrane-bound
form, the inactive membrane-bound form, and the inactive cytosolic form. Note that the membrane-bound
forms are deﬁned on the cell boundary vertices, and may exhibit spatial distributions. In fact, cell
polarization will be marked by spatially segregated distributions of active Rac1 and RhoA.
The cytosolic levels are functions of time but not space; we assume the cytosol to be well
mixed due to fast bulk diffusion [19, 20]. The total amounts of Rac1 and RhoA are each
On the membrane, the active and inactive forms interconvert with activation and deactivation rates,
for RhoA. Only the inactive form of protein may dissociate from the membrane to diffuse
within the cytosol. The membrane association and dissociation rates are denoted by
and by and
for RhoA. Using Fick’s law to
compute the 1D diffusion flux
of the active and inactive Rac on the membrane, we discretize the reaction-diffusion equations for Rac1 as
where the subscript
indicates the vertex on
the cell membrane, is the
edge length between vertex
and , and
is the diffusivity on the
the diffusive ﬂux from vertex
to vertex , its
for the active and inactive forms of Rac. The amounts of Rac1
normalized by the total amount of Rac1 in the cell. Similar equations can be written for the normalized amounts
of RhoA ,
, but we
omit them for brevity. These equations imply conservation of the total amount of Rac1 and
The biochemical interactions are encoded in the activation and inactivation rates
according to biological observations and prior modeling in the literature . They give rise to cell
polarization, stochastic repolarization (i.e. random changes in migration direction), CIL and
- Polarization. To capture cell polarity, the activation rates of Rac1 and RhoA each reﬂect
the species’ autocatalytic capacity, while their de-activation rates reﬂect mutual inhibition on
the cell membrane . Following [19, 21], we represent the Rac and Rho auto-activation
and mutual inhibition through Hill functions, which allow spontaneous polarization of cells
with Rac1 and RhoA peaking on opposite sides of a cell. This polarity is the precursor of cell
- Stochastic repolarization. NCCs produce random Rac1-mediated protrusions  that compete
with existing protrusive fronts and potentially change the cell’s polarization and direction
of motion. To model this, we randomly select a subset of vertices and up-regulate the Rac1
activation rate on them every
minutes. The model also accounts for competition between protrusions by a negative feedback
of the cortical tension on Rac1 activation . Thus, hotspots of Rac1 activity compete
with each other on the membrane, allowing an upstart to replace an existing protrusion on
occasion. This concept has been confirmed by recent biological observations on chick embryos .
- Contact inhibition of locomotion. Contact between two NCCs is known to activate the non-canonical
Wnt signalling pathway, which results in the down-regulation of Rac1 and the up-regulation of
RhoA , leading to CIL. This is effected in the kinetic equations by CIL factors that elevate
the Rac1 deactivation rate
and the RhoA activation rate
on any vertex that has come sufficiently close to a neighboring cell.
- Co-attraction. We realize COA by up-regulating the Rac1 activation rate on any vertex of
a cell that is sufficiently close to a neighboring cell to “sense” the C3a distribution of the
latter . This enhances protrusion toward each other between neighbors and produces COA.
Since C3a diffuses fast through the extracellular matrix, we assume a steady state exponential
distribution of C3a surrounding each NCC .
As in previous vertex models, the mechanical model represents the cell’s shape and movement via the position
of each vertex
and its speed
is the linearly elastic tension along the edge between vertices
is a cytoplasmic pressure that
resists cell area reduction,
is the protrusion or retraction force due to actin ﬁlaments or myosin motors on the membrane, and
is the unit outward
normal vector at vertex
(Fig. 2). Details of the mechanical model, along with the evaluation of parameters and a study of the
model’s sensitivity to the parameters, can be found in the full article .
The biochemo-mechanical model described above is able to reproduce single-cell behavior such
as polarization, motility and stochastic repolarization. Movie A gives an overall impression.
Besides, contact inhibition and co-attraction have been captured, as illustrated in Movies B and C, respectively. In this context, POP emerges naturally from the
simultaneous action of CIL and COA, as illustrated by the two-cell migration simulation discussed
3.1 The origin of POP
Movie D depicts the persistent directional migration of two cells. At the start of the simulation,
the cells are assigned random initial Rac1 and RhoA distributions, and placed next to the
left end of the corridor, whose boundaries confine the cells via CIL in a similar way as
cell-cell contact. Because of this confinement and the initial asymmetry in the geometry, the
two cells develop protrusion fronts toward the right, and the cells start to move to the right
( min). What
is somewhat surprising is that this motion is sustained for the entire duration of the simulation (10 hours) across
cell diameters, the initial portion of the migration being shown in Fig. 3. This contrasts the behavior of a
single cell that, after the initial movement away from the left end, essentially executes a 1D random walk
in the corridor.
Figure 3: Persistence of polarity exhibited by two cells migrating within a 1-cell-diameter wide
corridor. The lengths of line segments pointing in or out are proportional to the local levels of
active GTPases (outwards blue lines for Rac1, inwards red lines for RhoA). Highly persistent group
motion develops due to the trailing cell constantly neutralizing any new protrusions on the rear end
of the leading cell through CIL, and the leading cell reinforcing forward protrusions on the trailing
cell through the effect of COA.
We note first that COA keeps the two cells within close proximity for the entire duration. This ensures
the continual interaction between the two through CIL and COA. The effect of CIL on the leading cell is
such that it never develops a viable new protrusion in the rear. Any such hotspot, as appears at 82 min, is
quickly extinguished by CIL as the budding protrusion comes into contact with the trailing cell (116 min).
If such contact weakens the protrusion front of the trailing cell, this effect is short-lived as a slow-down of
the trailer will end the contact. The nascent Rac1 hotspot on the rear of the trailing cell, visible at 133
min, cannot compete with the forward protrusion that is reinforced by COA. As a result,
the pair continues its directional migration, and similar cycles of interaction repeat in time
(e.g. 200 min). Thus, COA and CIL act together to suppress the bursts of Rac1 up-regulation
on the cell membrane, which would have produced random repolarization on a single cell.
POP arises naturally from the pull-push action due to COA and
CIL, and perpetuates the initial asymmetric motion of the cells endowed by the geometric
This insight answers the main questions that had motivated our model. If CIL and COA are
postulated as ad hoc rules on the supracellular scale, they are insufficient for spontaneous collective
migration. From this we have hypothesized that POP is a necessary third ingredient. Having built our
model from the underlying GTPase dynamics on the intracellular scale, however, we find that POP
emerges from the collaboration between COA and CIL, and need not be added separately on this scale.
What has been added to our model, and missing from prior rule-based models [12, 13], is the Rac-Rho
biochemistry that allows a description and rational explanation of random repolarization. The fact that
Rac1 suppression promotes persistence in cell motility is well established experimentally for several cell
types [25, 26, 27, 23], and indeed this was what motivated our modeling on the level of intracellular
It is important to note that POP is stochastic in nature. If, by chance, new Rac1 hotspots appear
simultaneously on the rear end of both cells, they may overcome the forward polarity and cause the pair
to reverse course. Then POP would fail. Whether this occurs depends on the initial Rac1 and RhoA
distributions, and on the number of cells in the cluster. We will discuss the robustness of POP
3.2 Spontaneous collection migration of larger clusters
The key features of the persistent migration of two cells carry over to the collective migration of clusters
comprising more cells. Figure 4 and Movie E demonstrate the spontaneous collective migration of 49 cells
down a corridor. The cluster migrates persistently as a whole for about 20 cell diameters during a 10-hour
period. Although the cluster area fluctuates, COA keeps the coherence of the group so that the mean area
does not increase over time.
Figure 4: Spontaneous collective migration of a 49-cell cluster down a corridor of width
denotes the cluster centroid location and the snapshots show that the cluster retains its integrity
while migrating collectively. See Movie E for a more detailed view of the migration.
To quantify the cluster size effect on the spontaneous collective migration, we have to
distinguish it from the confinement effect of the corridor. With a larger number of cells inside the
same corridor, the cluster will be more crowded and this amounts to an effectively stronger
confinement. To separate the cluster size and confinement effects, we vary the width of the corridor
, scaled by the cell
diameter , with
the cell number
Thus, we always have a square initial array of cells enclosed on three sides by walls, as in Fig. 4,
This gives us, in a sense, “equal confinement” among the different
Figure 5: Effect of the cluster size
on key features of spontaneous collective migration. For each
value, the Tukey box-whisker represents 20 independent runs from different random
initial GTPase distributions. The top and bottom of the box mark the first and third
quartiles of the data, with the median indicated by the bar inside. The top and bottom
of the whiskers mark 1.5 times the interquartile range from the edges of the box. (a) As
increases from 4 to 49, we have modulated the strength of
COA to maintain roughly the same nearest neighbor separation
throughout the simulation. (b) The group persistence ratio
over 10 hours, computed from the distance between the start and end point of the centroid
of the cluster divided by the arc length of the centroid trajectory. (c) The persistence time
of the centroid of the cluster, deﬁned as the characteristic time for the
decay of the directional autocorrelation . (d) The average migration speed
of the centroid of the cluster. (e) The migration intensity
which gives us a characteristic frequency at which the cells would pass an observer fixed in the
Figure 5 depicts the variation of several indicators of collective migration with the cluster size
. In panel (a), the averaged
cell-cell separation is maintained
roughly constant for different
values by adjusting the strength of co-attraction. As co-attraction from different neighbors is additive, having a larger
and thus more
neighbors will elevate the COA effect. This is normalized by reducing the pairwise COA strength so as to produce the same
. Panels (b–e) shows the
persistence ratio , persistence
time , average centroid
migration speed and the
migration intensity .
Their definitions are given in the caption to Fig. 5. Taken together, these data
show that the intensity of persistent collective migration increases with increasing
, rapidly for
small clusters (4 to 9) and mildly for the larger ones.
Figure 6: The trajectories of the cluster centroid
for various group sizes. The cell diameter .
Each plot includes 20 runs starting from different random initial Rac1 and
RhoA distributions under otherwise identical conditions. Some of the trajectories
suffer reversals, but the collective migration becomes more robust with increasing
The explanation lies in the fact that for small clusters, POP has a greater probability of
failure, when the cluster reverses course as a whole. For larger clusters, the chance of POP
failure is much reduced. Figure 6 illustrates the stochastic nature of POP fallibility. Take
example. Among the 20 trajectories, the 10 that have traveled the farthest have never suffered a reversal,
for the entire duration. For the rest, reversal may first set in as early as
60 min or
as late as 520 min. The failure of POP requires a large portion of the group to develop rear protrusions
simultaneously. If only one or a few cells acquire rear hotspots among many, they will be quick extinguished by
their neighbors through CIL. This explains why the probability of POP failure declines with increasing
all, as POP arises from cell-cell interactions, it should be more robust in larger clusters with more frequent
and reliable interactions.
Finally, we can make some quantitative comparisons between our model prediction and experimental data.
For NCCs migrating in vivo in Xenopus laevis embryos, Szabo et al. 13] reported a collective migration
/min. In vitro, the
NCCs migrate at
down a corridor coated with fibronectin and conﬁned by borders that are rich in versican, an extracellular
protein that repels NCCs. With the parameters used here , our model predicts an average speed
/min for the larger
clusters () (Fig. 5d),
comparable with the experimental data. Note that in our model, we have adjusted the contractile force factor in
(Eq. 5) to produce the
experimental speed of 3 /min
for a single cell in the rectilinear phase of motility. The agreement in cluster migration speed
is not the result of parameter fitting. For the persistence ratio
, Szabo et al. 
documented over a time
period of 4 hours in vivo, and
in vitro. In Fig. 5, we have predicted an average persistence of
for the larger
over 10 hours, comparable to the experimental values. A caveat about comparison with experiment is that
the model has a large number of parameters, not all of which can be ascertained directly from
In the original article , we have also explored how the confinement in the corridor affects the
migration of the cell cluster, and found that optimal migration efficiency is achieved when the corridor width
Interested readers can consult the full paper.
Our goal is to provide an explanation of spontaneous collective migration of neural crest cells in terms of
key GTPases that regulate cell polarization and protrusion. From a broad perspective, we can summarize
this work by the following three points.
First, we have adopted a modeling approach based on the biochemical signalling pathways as an
alternative to the existing paradigms of continuum and agent-based modeling. Continuum models seek a
“mean-ﬁeld” description of the average orientation and movement of many cells without resolving the
scales of individual cells . However, it is challenging to represent intercellular interaction within the
continuum formulation. Agent-based models address this concern by assigning rules on individual agents
that recapitulate known cell-level behavior [8, 12, 13]. While helpful when little is known about the
underlying biology, these models tend simply to give back the behavior designed for. The
biochemically-based model uses a finer level of resolution. Moreover, it seeks to connect cell-level behavior
to intracellular signaling pathways that have been established experimentally. Thus, it can
be more general and even simpler than models based on postulated rules. For example, our
model predicts cell behaviors ranging from polarization to spontaneous collective migration
from a few well-establish principles governing GTPases  and cell mechanics . We
should note that in recent years, biochemically-based modeling of cell mechanics has gained
increasing interest and currency [32, 33]. This is evidently motivated and enabled by growing
experimental knowledge of the role of chemical signaling in complex behaviours such as collective
Second, we have identified the persistence of polarity, or POP, as an essential factor in enabling
spontaneous collective migration of cell clusters. Prior to this work, the extensive experimental studies of
Mayor and coworkers had identified CIL and COA as key to spontaneous collective migration. In their
recent models, Woods et al.  treated the cells as particles that move with inertia, according to
Newton’s second law of motion. Szabo et al.  used lattice-based rules that preferred the direction of
motion of the preceding Monte Carlo step in a cellular Potts model. Although POP appears merely as an
incidental feature of these models, we have come to suspect that it may in fact be essential for the
appearance of spontaneous collective migration. To test this hypothesis, we built a simpler model that
recapitulates CIL and COA but with POP intentionally suppressed by erasing the cell polarity
periodically and replacing it by a random “initial distribution” of the Rho GTPases. This
model does not exhibit spontaneous collective migration in a corridor; the centroid of the cell
cluster executes a 1D random walk once the cluster has moved away from the end of the
Third, we have demonstrated that POP emerges naturally from the simultaneous action of COA and
CIL. Thus, we have not only validated the existing proposal that COA and CIL cause spontaneous
collective migration [10, 4], but also placed it on the concrete biological basis of Rac-Rho
signaling. COA maintains the integrity of a cell cluster and ensures continual proximity and CIL
interaction among neighbors. The two cooperate to suppress Rac1 and perpetuate an initial
polarity and direction of migration that may have arisen from asymmetric boundary conditions
and confinement. This is consistent with experimental observations that Rac1 suppression by
drugs or Syndecan-4 enhances persistence in the polarity and motility in fibroblast and NCCs [25, 26, 27]. Thus, the model offers an explanation for the origin of POP; it is an emergent
behavior based on COA and CIL rather than a separate mechanism to be postulated in addition
to these two. Conceivably, this hypothesis can be tested by experiments that use drugs to
partially suppress various proteins (or their activation levels) in the Rho GTPase signaling
It is interesting to put the current model in the context of directional motion and symmetry-breaking
in vastly different systems, e.g. the active swimming of bacteria and flocking of insects or birds .
While alignment and directionality in macroscopic flocks rely on visual cues or pressure waves, we have
shown that in the NCC context, it emerges from the pull-push actions of COA and CIL. Generally, one
could view our POP as a device for alignment similar to the previously proposed directional persistence of
random walk , persistence due to particle inertia , “persistence decay” for cell polarity  and
collision-based alignment .
 R. McLennan, L. J. Schumacher, J. A. Morrison, J. M. Teddy, D. A. Ridenour, A. C.
Box, C. L. Semerad, H. Li, W. McDowell, D. Kay, P. K. Maini, R. E. Baker, P. M. Kulesa,
Vegf signals induce trailblazer cell identity that drives neural crest migration, Developmental
biology 407 (1) (2015) 12–25.
 A. Shellard, R. Mayor, Chemotaxis during neural crest migration, Semin. Cell Dev. Biol.
55 (2016) 111–118.
 R. Mayor, E. Theveneau, The neural crest, Development 140 (11) (2013) 2247–2251.
 C. Carmona-Fontaine, E. Theveneau, A. Tzekou, M. Tada, M. Woods, K. M. Page,
M. Parsons, J. D. Lambris, R. Mayor, Complement fragment c3a controls mutual cell
attraction during collective cell migration, Dev. Cell 21 (6) (2011) 1026–1037.
 A. Szabó, R. Mayor, Modelling collective cell migration of neural crest, Curr. Opin. Cell
Biol. 42 (2016) 22–28.
 E. F. Boer, E. D. Howell, T. F. Schilling, C. A. Jette, R. A. Stewart, Fascin1-dependent filopodia are required for directional migration of a subset of neural crest cells, PLoS Genet.
11 (1) (2015) e1004946.
 A. J. Burns, J.-M. M. Delalande, N. M. Le Douarin, In ovo transplantation of enteric
nervous system precursors from vagal to sacral neural crest results in extensive hindgut
colonisation, Development 129 (12) (2002) 2785–2796.
 S. Huang, C. Brangwynne, K. Parker, D. Ingber, Symmetry-breaking in mammalian
cell cohort migration during tissue pattern formation: Role of random-walk persistence,
Cytoskeleton 61 (4) (2005) 201–213.
 S. R. K. Vedula, M. C. Leong, T. L. Lai, P. Hersen, A. J. Kabla, C. T. Lim, B. Ladoux,
Emerging modes of collective cell migration induced by geometrical constraints, Proc. Natl.
Acad. Sci. U.S.A. 109 (32) (2012) 12974–12979.
 C. Carmona-Fontaine, H. K. Matthews, S. Kuriyama, M. Moreno, G. A. Dunn,
M. Parsons, C. D. Stern, R. Mayor, Contact inhibition of locomotion in vivo controls neural
crest directional migration, Nature 456 (7224) (2008) 957–961.
 R. Mayor, E. Theveneau, The role of the non-canonical wnt–planar cell polarity pathway
in neural crest migration, Biochem. J 457 (1) (2014) 19–26.
 M. L. Woods, C. Carmona-Fontaine, C. P. Barnes, I. D. Couzin, R. Mayor, K. M. Page,
Directional collective cell migration emerges as a property of cell interactions, PloS One 9 (9)
 A. Szabó, M. Melchionda, G. Nastasi, M. L. Woods, S. Campo, R. Perris, R. Mayor,
In vivo confinement promotes collective migration of neural crest cells, J. Cell Biol. 213 (5)
 B. Merchant, L. Edelstein-Keshet, J. J. Feng, A Rho-GTPase based model explains spontaneous collective migration of neural crest cell clusters, Dev. Biol. (in press 2018). URL https://doi.org/10.1016/j.ydbio.2018.01.013
 M. M. Zegers, P. Friedl, Rho GTPases in collective cell migration, Small GTPases 5 (3)
 A. J. Ridley, Rho GTPase signalling in cell migration, Curr. Opin. Cell Biol. 36 (2015)
 A. G. Fletcher, M. Osterﬁeld, R. E. Baker, S. Y. Shvartsman, Vertex models of epithelial
morphogenesis, Biophys. J. 106 (2014) 2291–2304.
 H. Lan, Q. Wang, R. Fernandez-Gonzalez, J. J. Feng, A biomechanical model for cell
polarization and intercalation during drosophila germband extension, Phys. Biol. 12 (2015)
 Y. Mori, A. Jilkine, L. Edelstein-Keshet, Wave-pinning and cell polarity from a bistable
reaction-diffusion system, Biophys. J. 94 (9) (2008) 3684–3697.
 B. Vanderlei, J. J. Feng, L. Edelstein-Keshet, A computational model of cell polarization
and motility coupling mechanics and biochemistry, Multiscale Model. Simul. 9 (4) (2011)
 B. Huang, M. Lu, M. K.
Jolly, I. Tsarfaty, J. Onuchic, E. Ben-Jacob, The three-way switch operation of Rac1/RhoA
GTPase-based circuit controlling amoeboid-hybrid-mesenchymal transition, Sci. Rep. 4 (2014)
 A. R. Houk, A. Jilkine, C. O. Mejean, R. Boltyanskiy, E. R. Dufresne, S. B. Angenent,
S. J. Altschuler, L. F. Wu, O. D. Weiner, Membrane tension maintains cell polarity by
confining signals to the leading edge during neutrophil migration, Cell 148 (1) (2012) 175–188.
 M. A. Genuth, C. D. C. Allen, T. Mikawa, O. D.
Weiner, Chick cranial neural crest cells migrate by progressively refining the polarity of their
protrusions, https://doi.org/10.1101/180299, accessed: 2017-10-10.
 A. Roycroft, R. Mayor, Molecular basis of contact inhibition of locomotion, Cell. Mol. Life
Sci. 73 (6) (2016) 1119–1130.
 R. Pankov, Y. Endo, S. Even-Ram, M. Araki, K. Clark, E. Cukierman, K. Matsumoto,
K. M. Yamada, A Rac switch regulates random versus directionally persistent cell migration,
J. Cell Biol. 170 (5) (2005) 793–802.
 M. D. Bass, K. A. Roach, M. R. Morgan, Z. Mostafavi-Pour, T. Schoen, T. Muramatsu,
U. Mayer, C. Ballestrem, J. P. Spatz, M. J. Humphries, Syndecan-4–dependent Rac1
regulation determines directional migration in response to the extracellular matrix, J. Cell
Biol. 177 (3) (2007) 527–538.
 H. K. Matthews, L. Marchant, C. Carmona-Fontaine, S. Kuriyama, J. Larraín, M. R.
Holt, M. Parsons, R. Mayor, Directional migration of neural crest cells in vivo is regulated
by Syndecan-4/Rac1 and non-canonical Wnt signaling/RhoA, Development 135 (10) (2008)
 R. Gorelik, A. Gautreau, Quantitative and unbiased analysis of directional persistence in
cell migration, Nat. Protoc. 9 (8) (2014) 1931–1943.
 J. Löber, F. Ziebert, I. S. Aranson, Collisions of deformable cells lead to collective
migration, Sci. Rep. 5 (2015) 9172.
 A. B. Jaffe, A. Hall, Rho GTPases: biochemistry and biology, Annu. Rev. Cell Dev. Biol.
21 (2005) 247–269.
 T. Wu, J. J. Feng, Modeling the mechanosensitivity of neutrophils passing through a
narrow channel, Biophys. J. 109 (11) (2015) 2235–2245.
 B. A. Camley, J. Zimmermann, H. Levine, W.-J. Rappel, Collective signal processing in
cluster chemotaxis: Roles of adaptation, ampliﬁcation, and co-attraction in collective guidance,
PLoS Comput. Biol. 12 (7) (2016) e1005008.
 J. Delile, M. Herrmann, N. Peyriéras, R. Doursat, A cell-based computational model
of early embryogenesis coupling mechanical behaviour and gene regulation, Nat. Commun. 8
 M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, R. A.
Simha, Hydrodynamics of soft active matter, Rev. Mod. Phys. 85 (2013) 1143–1189.
 L. Coburn, L. Cerone, C. Torney, I. D. Couzin, Z. Neufeld, Tactile interactions lead to
coherent motion and enhanced chemotaxis of migrating cells, Phys. Biol. 10 (4) (2013) 046002.
Supplemental Movies (mp4 format):
Movie A (283 KB): A single model cell showing random migration due to stochastic fluctuations in Rac1 activation rates. The lengths of line segments pointing in or out are proportional to the local levels of active GTPases (outwards blue lines for Rac1, inwards red lines for RhoA).
Movie B (1,142 KB): Two model cells undergoing contact inhibition of locomotion (CIL) upon contact in a 1-cell-diameter wide corridor. Co-attraction is turned off in this simulation to highlight the CIL repulsion.
Movie C (10,272 KB): Co-attraction (COA) maintains the integrity of a cell cluster, with nearest neighbors keeping a roughly constant average separation over time. Turning off COA between the cells would lead to dispersion of the group over time.
Movie D (552 KB): Two model cells exhibiting spontaneous collective migration in a 1-cell wide corridor.
Movie E (24,638 KB): A group of 49 cells exhibiting spontaneous collective migration in a 7-cell wide corridor.