Collective cell migration (CCM) and epithelial–mesenchymal transition (EMT) are common to cancer and morphogenesis, and are often considered to be mutually exclusive in spite of the fact that many cancer and embryonic cells that have gone through EMT still cooperate to migrate collectively. Here we use neural crest (NC) cells to address the question of how cells that have down-regulated cell–cell adhesions can migrate collectively. NC cell dissociation relies on a qualitative and quantitative change of the cadherin repertoire. We found that the level of cell–cell adhesion is precisely regulated by internalization of N-cadherin downstream of lysophosphatidic acid (LPA) receptor 2. Rather than promoting the generation of single, fully mesenchymal cells, this reduction of membrane N-cadherin only triggers a partial mesenchymal phenotype. This intermediate phenotype is characterized by an increase in tissue fluidity akin to a solid-like–to–fluid-like transition. This change of plasticity allows cells to migrate under physical constraints without abolishing cell cooperation required for collectiveness.

## Introduction

Epithelial–mesenchymal transition (EMT) is essential during embryo development and found in common pathologies such as organ fibrosis and in the initiation of metastasis for cancer progression. EMT is a process that converts an epithelium into individual mesenchymal cells. Cells lose their apico-basal polarity and cell–cell adhesion, and gain migratory and invasive properties to become mesenchymal cells (Thiery et al., 2009; Hanahan and Weinberg, 2011; Lim and Thiery, 2012). However, not all EMTs go to completion, and cells can have various degrees of mesenchymal phenotypes. In particular, cell–cell adhesion can be partially conserved. Interestingly, the ability to retain stable cell–cell contacts does not correlate with the capability of undergoing collective cell migration (CCM), a process during which a group of cells cooperate to migrate in a coordinated manner. Indeed, collective behavior can be found in cells that have been described as epithelial, mesenchymal, or as having an intermediate phenotype (Rørth, 2009; Friedl et al., 2012; Theveneau and Mayor, 2013). It is unclear what such intermediate phenotypes represent and what advantage, if any, they would confer on cells compared with fully epithelial or mesenchymal phenotypes. In particular, this raises the question of the role of cell–cell adhesion remodeling during EMT, especially when the cell population that activates an EMT program has to subsequently undergo CCM.

Here we use the neural crest (NC) cell population to (1) explore how cell–cell adhesion is regulated in a collectively migrating cell population and to (2) assess the implication of maintaining or disrupting cell–cell adhesion during collective migration.

NC cells are a highly migratory and multipotent embryonic cell population, whose invasive behavior has been likened to malignant invasion (Mayor and Theveneau, 2013; Powell et al., 2013). It has been well characterized that the initiation of NC migration during embryo development requires activation of an EMT program, which involves a qualitative and quantitative change of cell adhesion (Sauka-Spengler and Bronner-Fraser, 2008; Duband, 2010; Theveneau and Mayor, 2012). Migratory Xenopus laevis NC cells have been described as a pseudoepithelial cell population that progressively disassemble their cell–cell junctions (Alfandari et al., 2010). In this system, cells become fully migratory before complete cell–cell dissociation, which allows us to address specifically the role of cell–cell dissociation during CCM in vivo.

Looking for candidate regulators of cell–cell adhesion, we found incipient data linking lysophosphatidic acid (LPA) signaling with changes in cadherin function during EMT in both cancer and NC cells (Smicun et al., 2007; Groysman et al., 2008; Kam and Quaranta, 2009; Huang et al., 2012; Liu et al., 2012). The cellular activities controlled by LPA signaling are diverse, including proliferation, cell motility, chemotaxis, tumor invasion, gap-junction closure, tight junction opening, etc. (Mills and Moolenaar, 2003). This diversity of biological functions, as well as some apparent different cellular responses triggered by LPA, is likely related to the fact that LPA can bind any of six distinct receptors (Lin et al., 2010). In addition, some level of redundancy has been described in mammalian embryos (Contos et al., 2000a,b, 2002), making impossible to characterize the biological activity of each LPA receptor using knockout mice.

Here we studied the specific role of LPA via LPAR2 during Xenopus NC cell development. Our results show that LPAR2 is essential for NC migration in vivo by controlling the internalization of N-cadherin. Importantly, LPAR2-induced cell–cell dissociation is critical only when cells are challenged into migrating through interstices in vivo or in vitro. Conversely, forcing cells to fully abolish cell–cell adhesion blocks migration in vivo but enhances dispersion in vitro. Thus down-regulation of cell–cell contacts has to be mild to endow the NC population with enough plasticity to invade its surroundings without abolishing cell cooperation required for collectiveness. LPA has been widely studied as a regulator of guanine nucleotide exchange factors (GEFs) for RhoA (Moolenaar et al., 2004; Xiang et al., 2013) but LPA function as a negative regulator of cell–cell adhesion by controlling N-cadherin endocytosis has never been shown before. Similarly, the role of cadherin endocytosis during EMT has been well characterized (Ogata et al., 2007), but its ability to control tissue plasticity necessary for in vivo CCM has not been previously revealed.

Our work offers novel insight on how cells can balance the need to migrate as individual cells, in an environment where physical constraints are present, with the need to maintain a certain degree of cooperation via cell–cell adhesion. We anticipate that similar fine-tuning of cell–cell contact during collective migration in vivo occurs in other systems such as carcinoma cells or gastrulating mesoderm.

## Results

We first looked at the expression of several components of the LPA signaling pathway and found that LPA receptor 2 (LPAR2) is expressed in the whole NC region at premigratory stages as compared with the NC marker FoxD3 (Fig. 1 A), and seems to be progressively lost at late migratory stages, as it is detected only in subpopulations of NC cells after these cells have reach the ventral part of the head (Fig. 1 B and Fig. S1 A). LPAR1 is found in the otic vesicle and in the ectoderm, but not in the NC (Fig. S1 B). Autotaxin (ATX), one of the extracellular enzymes required for LPA synthesis (Umezu-Goto et al., 2002; Houben and Moolenaar, 2011), is weakly expressed at premigratory stages and found in a subpopulation of migratory NC cells (Fig. S1 C). We designed an antisense Morpholino against LPAR2 and tested its ability to silence a reporter containing its target sequence (Fig. 1, C and D). Blocking LPAR2 (Fig. 1, E and F) or that of ATX (Fig. S1 D) impaired NC migration in vivo without affecting NC induction (Fig. S1, E and F). We then checked that LPAR2 was required only in NC cells and not in their surrounding tissues by grafting morphant cells into control host embryos, or the converse experiment of grafting control cells into an LPAR2MO-injected host (Fig. 1 G). Cells injected with Fluorescein-dextran (FDx) as a control were able to migrate after being grafted onto a host (Fig. 1 H, left). LPAR2MO cells did not migrate into the host embryo (Fig. 1 H, middle), whereas LPAR1MO cells migrated as well as control cells (Fig. 1 H, right). Finally, control NC cells grafted into LPAR2MO embryos migrated normally (Fig. 1 I). These results show that LPA signaling within the NC via LPAR2 (but not LPAR1) is required for NC migration.

LPAR2 phosphorylation increases at the onset of NC migration (Fig. 2, A and B). Although this increase in LPAR2 phosphorylation is quite reproducible, we were not able to show a statistically significant difference. Therefore we decided to test for a role of this phosphorylation on NC migration by mutating the phosphorylation sites. Two main putative phosphorylation sites are present in the cytoplasmic tail of LPA receptors: Tyrosine-317 (Y317) and Serine-324 (S324). However, S324 is the only site conserved across both LPAR1 and LPAR2 (Fig. S1, G and H), making it likely to be the most relevant site for signaling. Importantly, a single mutation targeting S324 in the cytoplasmic tail of LPAR2 (LPAR2SA) was sufficient to abolish phosphorylation of LPAR2 (Fig. 2 C). The inhibition of NC migration by LPAR2MO was rescued by overexpression of an mRNA encoding full-length LPAR2, which does not hybridize with LPAR2MO, but not by the nonphosphorylatable form LPAR2SA (Fig. 2, D and E). Together these results indicate that LPA signaling and phosphorylation of LPAR2 are essential for in vivo NC migration.

To better characterize which aspect of cell migration was affected in the absence of LPAR2, we performed in vitro cell migration assays and found that control and LPAR2MO cells had similar cell motility properties (Fig. 3, A and B). We then tested whether LPAR2MO cells might have an impaired response to signals known to guide NC migration in vivo. We performed chemotaxis assays where cells were exposed to an external source of stromal cell derived factor-1 (SDF1), a well-known NC chemoattractant (Belmadani et al., 2005, 2009; Olesnicky Killian et al., 2009; Kasemeier-Kulesa et al., 2010; Theveneau et al., 2010, 2013; Fig. 3, C and D). Both control and LPAR2MO-injected cells were able to respond to SDF1 and migrate toward the source with similar efficiency (Fig. 3 E and Video 1). These results show that blocking LPA signaling does not block NC cell motility or their ability to respond to chemotactic cues. Our data also show that in 2D culture assays, LPAR2-depleted NC cells are able to migrate normally. This is in striking contrast to the situation in the 3D environment of the living embryo, in which NC migration requires LPAR2. One possible explanation for this discrepancy may lie in the physical nature of the embryo: in vivo, NC cells encounter physical barriers (i.e., the eye, otic vesicle, epibranchial placodes) and chemical inhibitors (i.e., class 3 semaphorins, Eph/ephrins) that restrict migration into specific territories. Such constraints do not exist in an open 2D culture.

To test the idea that LPAR2 is required for cells to migrate into physical constraints, we developed a chemotaxis assay in which cells are forced to navigate through tunnels. The dimensions of these tunnels corresponded to the size of the NC streams observed in vivo (Fig. S2). NC cells were cultured on microfluidic chips containing tunnels of various widths corresponding to one nucleus diameter (10 µm), one cell diameter (25 µm), two cell diameters (50 µm), and six cell diameters (150 µm). All tunnels were one nuclei diameter in height (10 µm), thereby dictating that cells within a multilayered explant must rearrange themselves to squeeze into the tunnel (Fig. S2). We found that cells in control explants were capable of invading all channel widths (Fig. 4, A, C, and E), which indicates that substantial cell rearrangements had taken place. In contrast, LPAR2MO-injected cells consistently stalled at the entrance or in the first half of the tunnels (Fig. 4, B, D, and E; and Video 2), which suggests that LPAR2MO cells had a low plasticity and were deficient in their ability to exchange positions with their neighbors.

To test the idea that a lack of LPAR2 impairs the ability of cells to exchange neighbors, we quantified the relative motion of cells within a group in control and LPAR2MO conditions. We cultured cells for 1–3 h on a 2D substrate and characterized the degree of alignment of velocity vectors by calculating the correlation function of the cell velocity field (see Materials and methods). Alignment of velocity vectors is expressed as a function of the relative position of the cells averaged over time (Fig. 4, F and G) or as a function of radial distance between cells and time (Fig. 4, H and I). Velocity vectors are perfectly aligned where this function is 1, and randomly aligned where the function is 0. In other words, cells move in the same direction where this function is 1, whereas their movements are not correlated where this function is 0. In vitro and in silico studies (Kabla, 2012; Ladoux and Nicolas, 2012; Doxzen et al., 2013) have shown that a population appears solid-like at scales smaller than the correlation length, but free to reshape at larger scales. Therefore, the characteristic length at which this function decays (the distance between two cells at which coordination is lost) qualitatively defines the response of the cell population to confinement. The ability to flow in order to take on the shape of a container is a characteristic of fluids, whereas solids have a structural rigidity making them resistant to such changes of shape. Hence, cell populations behave as solid-like when exposed to structures that are smaller than their correlation length whereas they behave as fluids in spaces that are bigger than their correlation length.

We found that, in control groups, movement of cells located anywhere within the group showed little correlation with neighbors, with a correlation length on the order of 25 µm, roughly corresponding to one and a half cell diameters (Fig. 4 G, left) in accordance with our observation that control NC cells were able to reorganize to pass through tunnels of 25 µm. In contrast, LPAR2MO cells displayed a solid-like behavior even over distances between 150 and 200 µm, roughly six to eight cell diameters (Fig. 4 G, right). This result means that LPAR2MO cells may need an infinite amount of time to reorganize in order to take on a shape smaller than 150 µm. This can explain why LPAR2MO cells were only able to partially enter tunnels of 150 µm but not smaller ones. These dynamics measurements were also seen to be constant over time (Fig. 4, H and I). This result indicated that LPAR2MO cells, despite being fully motile and responding accurately to chemotactic cues, tend to behave as a solid at the scale of hundreds of microns, making them unable to rearrange to invade narrow spaces. By comparison, control cells emerge as a fluid-like material at the scale of only a few tens of microns and as such are able to reorganize to pass through interstices. Thus, the acquisition of a fluid-like behavior allows CCM into interstices, whereas inhibition of this phase transition by blocking LPAR2 makes cell clusters too solid to pass.

Such an increase in the solid-like behavior observed in LPAR2MO cells could be caused by a decrease in cell motility or an increase in the cohesion of the tissue. Because cell motility of single or clustered cells was not affected by LPAR2 depletion (Fig. 3, A and B), we hypothesized that the observed phenotype might be caused by increased cell–cell adhesion slowing down cell rearrangements within the group. To determine whether cells lacking LPAR2 displayed stronger cell–cell adhesions, we performed a 2D culture dispersal experiment using control and LPAR2MO-injected cells (Fig. S3). Under these conditions, control NC cells are known to progressively spread in a radial fashion with cells at the border of the group, eventually breaking away as single cells within 5 h (DeSimone et al., 2005). This dispersion is a direct consequence of the loss of cell–cell adhesion and can be quantified by measuring the area between the nearest neighbors using Delaunay triangulation (Carmona-Fontaine et al., 2011; Barriga et al., 2013). Control cells were able to efficiently disperse, whereas LPAR2MO cells remained as a group (Fig. S3 and Video 3). Altogether, our results show that LPAR2MO cells are as migratory and as responsive to external guidance cues as control cells but fail to dissociate and thus represent a partial EMT phenotype where motility is intact but cell–cell dissociation does not occur.

To further assess whether the altered N-cadherin traffic might explain the change of NC migratory abilities observed when LPA signaling is blocked, we tried to rescue LPAR2MO cells by manipulating N-cadherin binding, expression, or endocytosis. We first performed a cell dispersion assay (Fig. 7, A–G; and Video 5) in which NC cells injected with a control Morpholino normally break away as single cells (Fig. 7 A), whereas LPAR2MO cells fail to do so (Fig. 7 B). Dispersion defects could be rescued in LPAR2MO cells by lowering calcium (Fig. 7 C), blocking N-cadherin expression (Fig. 7 D), reducing N-cadherin binding in trans with a blocking antibody (Fig. 7 E), or increasing endocytosis with an activated Rab5 (Fig. 7 F and Video 5). Importantly, spatial correlation kymography analyses demonstrated that these treatments were also sufficient to dramatically lower the correlation length. This indicated that blocking N-cadherin in LPAR2MO cells was sufficient to restore a fluid-like behavior at scales similar to control cells (Fig. 7, H and I).

Further, we performed an invasion assay using microfluidic tunnels and found that LPAR2MO cells that were preincubated with the blocking antibody against N-cadherin (NCD2) were able to invade narrow spaces (Fig. 8, B, D, and E) in a manner similar to that seen in control cells (Fig. 8, A, B, and E; and Video 6). Finally, we tested the ability of a mild inhibition of N-cadherin or increased endocytosis to rescue migration of LPAR2MO cells in vivo. We first grafted LPAR2MO NC cells directly into control hosts or first incubated them with the blocking antibody against N-cadherin (NCD2) before transplantation (Fig. 8 F). LPAR2MO cells were unable to migrate after being grafted into hosts (Fig. 8, H and K), whereas LPAR2MO cells treated with NCD2 (Fig. 8, I and K) migrated as efficiently as control cells (Fig. 8, G and K). Similarly, coinjection of LPAR2MO with an mRNA encoding constitutively active Rab5 was sufficient to restore in vivo NC migration (Fig. 8, J and K). These results show that inhibition of N-cadherin or activation of N-cadherin endocytosis in LPAR2-deficient cells is sufficient to restore the plasticity required for NC migration under physical constraints in vitro and in vivo.

We have previously shown that cell–cell interactions are required for in vivo CCM of NC cells (Carmona-Fontaine et al., 2008; Theveneau et al., 2010; Carmona-Fontaine et al., 2011). Cell contacts polarize cells, rendering them able to sense chemotactic cues in a process coined collective chemotaxis (Theveneau et al., 2010). Therefore, lowering cell–cell adhesion even further than in control cells should favor cell dispersion over CCM by abolishing cell cooperation and thus inhibit NC migration in vivo but not in vitro. To test this idea, we analyzed the effect of reducing cell–cell adhesion in a 2D in vitro dispersion assay (Fig. 9 A) and found that NC cells cultured under low calcium/magnesium conditions (orange curve) or injected with a dominant-negative form of N-cadherin (blue curve) dissociate faster than control cells (green curve). Then we performed a similar experiment in vivo by reducing cell–cell adhesion by specifically blocking N-cadherin expression or using deletion constructs that are known to act as dominant-negatives (Fig. 9, B–E). Indeed, blocking N-cadherin expression (Fig. 9 B) or overexpressing its truncated forms lacking the cytoplasmic or β-catenin binding domains (Fig. 9, C and D) impaired in vivo NC migration (Fig. 9 E). These observations, together with previously published results (Theveneau et al., 2010), show that cells with low cell–cell adhesion dissociate faster but are unable to cooperate properly in order to undergo the directional collective migration characteristic of the cephalic NC population.

## Discussion

Thus, such a change of plasticity appears to be an essential part of the EMT program, which allows invasion into regions imposing physical constraints without impairing the cell–cell interactions needed for CCM. Similar fine-tuning of the cell–cell adhesion properties might happen in other systems where EMT and CCM coexist, such as cancer metastasis or gastrulation. In addition, the need to maintain efficient cell interactions to allow for collective migration might in part explain the counterintuitive observation that many tumors have de novo expression of cell–cell adhesion molecules while undergoing EMT (Wheelock et al., 2008; Thiery et al., 2009). Down-regulation of cadherins favors emigration of the primary tumor cells from their host organ (Nguyen et al., 2009; Polyak and Weinberg, 2009) but might also impair cooperation among cancer cells. Alternatively, altering cell–cell adhesion, without cancelling it, would allow cancer cells to both successfully detach from a tissue and preserve cooperative behaviors required for long-distance migration (i.e., contact-dependent polarity and survival).

The role of mechanical forces is highly debated in cancer metastasis (Wirtz et al., 2011). This led some to theoretically explore analogies between tumor metastasis and phase transition from solid to liquid to gas (Davies et al., 2011). The results described here provide experimental support and a molecular mechanism for such ideas, and indicate that in vivo CCM is the result of a balance between physical interactions among migratory cells and the mechanical properties of the local environments. Overall, our findings contribute to a better understanding of the relative importance of cell–cell versus cell–environment interactions during collective movements.

Previous studies have shown that inhibition of Rho/ROCK signaling in zebrafish impairs NC migration (Clay and Halloran, 2013). An opposite role for Rho/ROCK signaling controlled by LPA has been proposed in chick NC. Groysman et al. (2008) showed that adding LPA to trunk NC cultures activated RhoA/B and blocked cell emigration. It is not known which LPA receptors are present in trunk NC cells in chicks, and the phenotype induced by adding LPA in the culture medium is likely to be the result of the activation of multiple signaling cascades. Here we specifically blocked LPAR2 and showed an impairment in cell–cell dissociation. We do not know the effect that activation or inhibition of other LPA receptors could have on NC migration in our system.

Genetic studies with LPA receptor–null mice have suggested some degree of redundancy, as LPAR2-null mice have no obvious phenotype, but the LPAR1/LPAR2 double-null mice have a higher lethality than LPAR1-null mice (Contos et al., 2000a, 2002). It is common to find this degree of redundancy between signaling pathways in mice and not in other animal models like Xenopus, which allowed using this animal model to discover this novel function for LPAR2. Previous studies in Xenopus embryos have shown that LPAR2 is involved in wound healing by controlling the actin cytoskeleton (Tao et al., 2007). Whether the actin cytoskeleton is modulated by LPAR2 signaling in our cells remains to be investigated. Interestingly, our findings on the role of ATX and subsequent LPA signaling via LPAR2 are consistent with the biological activity reported in cancer models (Houben and Moolenaar, 2011), and suggest that Xenopus NC cells could be an adequate model to study cell behaviors important for cancer metastasis.

## Materials and methods

### PCR primers

#### For LPAR2 mutations.

Forward primer: 5′-TGAATTCATGAAcCAaTGtTTcTAcAATGAGAGTATTAAGTTC-3′. Mutations are represented by lower-case italicized letters. For the Ser324 point mutation of LPAR2, pCS2+/5 mm LPAR2 whole plasmid was amplified using the following primers. S324>A forward (Fw), 5′-ACTAATACAGAGAGGAGAGTCTTGGCTGAC-3′; S324>A reverse (Rv), 5′-TGCCTTAAGGGAAGCCGGTGAGTATTTGG-3′. For making 5′ UTR-LPAR2-GFP shown in Fig. S4, 5′ UTR-LPAR2 ORF cDNA was amplified by 5′ UTR-LPAR2 forward, 5′-TCGAATTCGACAGGCGGAGGC-TGAGCTGTA-3′; and AgeI-LPAR2 reverse primer, 5′-CGACCGGTGCAAGACTGCATTCTCTCATTG-3′.

N-cad Fw, 5′-CAGCATCGATGTGCCGGAAAGAGCCCTTCC-3′; Rv, 5′-AGTCTAGATTCAGTCGTCGCTCCCTCCGTAC-3′. First, the restriction sites for HaloTag fusion protein were introduced into XN-cadherin by PCR with the following sets of primers. N-cadSPFw, 5′-CGTTTATGCAAGACAGGATTTCTTGAGGATG-3′; N-cadSPRv, 5′-AGTTCCCCCGAGTGCTTCAACTGGTCCCT-3′; HaloTagFw, 5′-AGTATGGCAGAAATCGGTACTGGCTT-3′; HaloTagRv, 5′-CGTGGTTGGCTCGCCGGAAATCTC-3′. Then the pro-domain of N-cadherin was removed by another PCR, N-cadEC1Fw: 5′-AATCACGCGTGATTGGGTTATTCCACCAAT-3′.

#### For semiquantitative RT-PCR.

RNAs were extracted from whole embryos (2–3 embryos) or 20 pieces of dissected cranial NC cells using RNeasy Protect Mini kit (QIAGEN). The cDNA was reverse-transcribed by SuperScript II Reverse transcription (Life Technologies). For semiquantitative PCR, PCR cycles were decided by pilot experiment using the whole embryos’ cDNA. The program for amplification was 96°C for 1 min, 96°C for 15 s, 55°C for 30 s, and 72°C for 45 s × cycles as follows: Sox9 (29 cycles), Xβ-Cat (27 cycles), N-cadh (27 cycles), ornithine decarboxylase 1 (ODC; 23 cycles). Primers for RT-PCR for Ep-Ker, ODC, and muscle actin were described in Xenbase (Bowes et al., 2008), and Sox9 (Monsoro-Burq et al., 2003), Snail2 (Aybar et al., 2003), and FoxD3 (Sasai et al., 2001) have been described previously.

### Reagents and solutions

The following reagents and solutions were used. Agarose (BP1356-100; Thermo Fischer Scientific), BSA (A4503; Sigma-Aldrich), fibronectin (F1141; Sigma-Aldrich), FDx (D1821; Sigma-Aldrich), gelatin (G1890; Sigma-Aldrich), methylcellulose (M0387-100; Sigma-Aldrich), Mowiol 40-88 (81386; Fluka), penicillin-streptomycin (P4458; Sigma-Aldrich), Rhodamine-dextran (D1824; Sigma-Aldrich), sucrose and (27480.294; VWR). Danilchick’s medium 1×: 53 mM NaCl, 5 mM NA2CO3, 4.5 mM potassium gluconate, 32 mM sodium gluconate, 1 mM MgSO4-7H20, 1 mM CaCl2, and 0.1% BSA. MEMFA: 1 mM MOPS, 2 mM EGTA, 1 mM MgSO4, and 3.7% formaldehyde. Normal amphibian medium (NAM): 110 mM NaCl, 2 mM KCl, 1 mM Ca(CO3)2, 1 mM MgSO4, 0.1 mM EDTA, 1 mM NaHCO3, and 2 mM sodium phosphate. NTMT: 0.1 M NaCl, 0.1 M TrisHCl, pH 9.5, 50 mM MgCl2, and Tween 0.1%. PBS 1×: 137 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, 1.8 mM KH2PO4, 1 mM CaCl2-2H2O, and 0.5 mM MgCl2-6H2O.

### Microfabrication of microfluidic tunnels

Two-layer chips were designed using the AutoCAD 2009 software (Autodesk). High-resolution soda-lime masks for the microchannel layer, and emulsion masks for the chamber layers were obtained from JD Photo-Tools. 3-inch silicon wafers were ordered from Virginia Semiconductor Inc. Molds were processed in two steps in the cleanroom facility of the London Centre for Nanotechnology. AZ5214E resist (MicroChemicals GmbH) was spun at 5,000 rpm and processed as a negative resist using standard procedures to create the mask for the microchannel layer. The patterned wafer was then etched for 10 or 20 µm using a deep reactive ion etcher (Surface Technology Systems ICP DRIE; SPTS Technologies), plasma cleaned in a TePla plasma cleaner (PVA TePla) and further cleaned with fluoric acid. A second 150-µm-thick SU-8-2100 photoresist (MicroChem Corp.) layer was then deposited and processed using standard microfabrication techniques (García and Henry, 2006) to generate the mold for the explant culture chamber and the chemoattractant channel. The mold was finally cleaned once more in the TePla plasma cleaner. It was then plasma activated (using a Plasma Clean 4; Ilmvac), plunged in 1% tridecafluoro-1,1,2,2-tetrahydrooctyl-1-dimethylchlorosilane (13F; MCC) in toluene for 30 min, washed, dried, and baked at 55°C for 15 min. Polydimethyl methylhydrogen siloxane (PDMS) casts were prepared by mixing 1:10 parts of curing agent to base Sylgard 184 elastomer (Dow Corning Corp.). The PDMS mix was degassed and poured onto the 13F-coated patterned wafer and cured at 55°C overnight. Then, input/output were created using a biopsy punch and the culture chamber was opened by cutting out a square of PDMS with a surgical scalpel. PDMS casts were further exposed to plasma for 30 s and irreversibly bonded to 35-mm-diameter glass-bottom dishes (WillCo; Intracell). Finally, the assembled chips were baked for 15 min at 75°C.

### Spatial correlation analysis

This method allows us to assess the degree of alignment between cells within a group and to express this alignment as a function of distance. Nuclei are tracked over time and their displacement between two consecutive time points give a series of velocity vectors that can then be systematically compared with one another. For statistical comparisons, one has to choose a threshold of alignment. Then, the mean distance between cells at which this mean alignment is reached is used to compare experimental conditions. We chose the distance at which the mean alignment of velocity vectors across the population is 0.5 to compare experimental conditions. Further detail of the strategy to obtain velocity vectors is provided hereafter.

To quantify the spatial correlations of cell displacements, the velocity fields are first calculated using an optical flow method on images of stained nuclei (cells were injected with mRNA for nuclear-mCherry). Displacements are quantified from pairs of images 5 min apart, at a resolution of 2.05 µm/pixel. The spatial correlation function of the velocity field is then calculated at each time step using the methodology developed by Kabla (2012). The 2D correlation function g(δr) is introduced to characterize the mean degree of coordination between a cell moving at the origin moving along the x axis, and another cell located at a position δr with respect to the first. This enables one to distinguish correlations along the direction of motion and correlations in the transverse direction. Quantitatively, the correlation function at any instant is calculated from:

$g(δr)=〈v(r)·v(r+R(r)δr)〉r〈v2〉r,$

where R(r) is the transformation rotating the velocity vector v(r) along the horizontal x axis. The normalization of the function ensures that g = 1 where cells’ velocities are fully aligned, and 0 where orientations are random.

Spatial maps were obtained by averaging over all time steps between 0 and 75 min. Time-resolved maps were obtained by calculating, at each time point, the radial mean of g(δr) over 200 min (3 h, 20 min). The distance at which such function decays can be interpreted as the correlation length of the cell population.

### Statistical analysis

Comparison of percentages was performed using contingency tables as described previously (Taillard et al., 2008). Two datasets were considered significantly different (null hypothesis rejected) if T > 3.841 (*, α = 0.05), T > 6.635 (**, α = 0.01), or T > 10.83 (***, α = 0.001). Normality of datasets was tested using the Kolmogorov-Smirnov’s test, d’Agostino and Pearson’s test, and Shapiro-Wilk’s test using Prism 4 (GraphPad Software). A dataset was considered normal if found as normal by all three tests. Datasets following a normal distribution were compared with a Student’s t test (two-tailed, unequal variances) in Excel (Microsoft) or a one-way analysis of variance (ANOVA) with a Dunnett’s multiple comparisons post-test in Prism 4. Datasets that did not follow a normal distribution were compared using a Mann Whitney’s test or a nonparametric ANOVA (Kruskal-Wallis with Dunn’s multiple comparisons post-test) using Prism 4. Cross-comparisons were performed only if the overall P-value of the ANOVA was <0.05. Scatter plots show individual values for each individual experiment in each condition, and error bars show the mean and the standard deviation. Box and whiskers plot: the box extends from the 25th to the 75th percentile; the whiskers show the extent of the whole dataset. The median is plotted as a line inside the box.

### Delaunay triangulation

In brief, individual nuclei were first automatically detected using the Particle Analysis tool in ImageJ. Then, xy positions of all nuclei were used to build triangles between the nearest neighbors with the Delaunay Triangulation plug-in. Areas of individual triangles are extracted. An average triangle area is calculated for each explant and the values were then normalized according to the mean of the control batch of each experiment. It is important to note that triangle area distribution is exponential. Thus, to perform statistical analysis, the datasets are converted into logarithms of areas that follow a normal distribution. Paired and multiple statistical comparisons were made as described in the Statistical analysis section.

### Cell tracking

Cell tracking was performed with the Manual Tracking plug-in in ImageJ or with spot tracking function in Imaris (Bitplane). In all cases nuclear staining was used to follow individual cells over time. Directionality, velocity, and chemotaxis index (also known as the forward motion index) were extracted using the Chemotaxis Tool plugin in ImageJ or calculated from data generated by the Imaris tracking tool.

### NC grafts

The appropriate morpholino oligo and FDx or rhodamine-dextran (RDx; 1 µg/µl, mol wt 1,000; Molecular Probes) were injected into two-cell-stage embryos. Embryos at stage 18 were immobilized in a Petri dish filled with plasticine. The pigmented ectoderm layer located above the NC region was carefully removed using an eyebrow knife. NC cells can be mechanically detached from their surrounding tissues by applying gentle pressure on the side of the neural fold, passing the eyebrow knife from lateral to medial a few times. To perform the graft, a given NC explant has to first be taken from a donor embryo, then a similar procedure is immediately done on the host embryo. The NC explant coming from the donor is then placed into the wound of the host and kept in place by placing a square glass coverslip on top for 15 min. The coverslip is then removed and embryos are left to heal. Further details on grafting procedures can be found elsewhere (DeSimone et al., 2005; Theveneau et al., 2013).

### Gene knockdown

Gene knockdown of Autotaxin, LPAR1, LPAR2, and N-cadherin (Nandadasa et al., 2009) was done by injection of morpholino oligonucleotides (MOs; Gene Tools, LLC). Morpholinos sequences were as follows XlAutotaxin, 5′-CTGAATCCGTTTTTCATTGCCATAG-3′; XlLPAR1, 5′-ACTCTGAAAGAGAAGCCATTTTAGC-3′; and XlLPAR2, 5′-ATAAAAGCACTGATTCATCTTGAGT-3′. ATX morpholino was provided by L. Dale (University College London, London, England, UK).

### In situ hybridization

Embryos were fixed overnight at 4°C in MEMFA before being dehydrated by several baths in pure methanol. Embryos were then rehydrated by solutions of decreasing methanol concentration, washed in PBS, and bleached in a peroxide solution to remove skin pigments. After bleaching a short postfixation in formaldehyde 3.7% was performed. Embryos were then incubated overnight at 65°C in formamide-based hybridization buffer containing a digoxigenin-labeled antisense probe against the gene of interest. Probes are washed in Formamide-based washing solutions, then washed in PBS 0.1% Tween, incubated in a serum-based blocking solution, and left overnight at 4°C in blocking solution containing the anti-digoxigenin antibody coupled with alkaline phosphatase. Staining was achieved by incubating embryos in a pH 9.5 solution containing NBT and BCIP in equal proportions.

### Immunostaining

Cells were cultured on fibronectin-coated dishes, left to migrate for a few hours, fixed in 4% PFA for 30 min, blocked in a serum-based solution, incubated for 2 h with a primary antibody, quickly washed in PBS, and incubated 1 h with a secondary antibody mixed with DAPI or phalloidin if necessary. Primary antibodies: rat anti–N-cadherin (MNCD2; Developmental Studies Hybridoma Bank), rabbit anti–β-catenin (ab6302; Abcam), and mouse anti-vinculin (V9131; Sigma-Aldrich). Secondary antibodies were all Alexa Fluor–coupled immunoglobulins from Invitrogen (Alexa Fluor 488 or 555).

### NC cell culture

Xenopus cranial NC cells were dissected as described in DeSimone et al. (2005). In brief, at around stage 18 the pigmented epidermal layer is removed, then NC cells are gently taken out by microdissection. When needed, cell dissociation was performed by putting the NC explants in Ca2+/Mg2+-free medium for a few minutes before transferring them to normal culture medium. All cell cultures were performed on Fibronectin-coated Petri dishes.

### cDNAs and expression constructs

cDNA clones were obtained from the I.M.A.G.E. consortium. Xenopus autotaxin, I.M.A.G.E. 5570505; Xl-LPAR2 (EDG4), I.M.A.G.E. 6864256. Xenopus LPAR1 cDNA was a gift from J. Chun (The Scripps Research Institute, San Diego, CA). Specificity of LPAR2 MO was tested with a mutant LPAR2 including five mismatches in MO target sequences, and different LPAR2 mutants were also created by PCR (sequence given in the Gene knockdown section). For HaloTag–N-cadherin, Xenopus N-cadherin ORF was amplified by primers shown in the PCR primers section, and subcloned into pCSf107mT vector. HaloTag cDNA was inserted between N-cadherin prodomain and EC1.

### Pulse-chase assay (HaloTag in situ labeling)

At the eight-cell stage, HaloTag N-cadherin mRNA (500 pg/embryo) and membrane-mCherry mRNA (300 pg/embryo) were injected, and the embryos were cultured until stage 19. HaloTag-N-cadherin–expressing NC explants were dissected and preincubated for 30 min at room temperature with Alexa Fluor 488–nonpermeable HaloTag ligand solution in 1× Marc’s modified Ringer (MMR; 1:1,000, Promega; Los and Wood, 2007), and the excess ligands were washed. Then labeled NC cells were cultured on fibronectin-coated coverslips (Tissue Tek II; Nalgene), and time-lapse imagining was performed (LSM 510 Meta; Carl Zeiss) immediately after cell attachment. The lens used for this imaging was a Plan-Apochromat 63×/1.4 oil objective. The temperature was ∼23°C. The imaging medium was Immersol 518F (Carl Zeiss). Images were obtained and analyzed with LSM 510 version 3.2 microscope operation software and Zeiss Image Browser (both from Carl Zeiss).

### Phostag PAGE and Western blotting

Phos-tag PAGE and Western blotting were performed as described previously (Kinoshita et al., 2006; Kuriyama and Mayor, 2009; Tanaka et al., 2012). The protein samples for phosphorylation analyses were prepared from Xenopus tissues lysed in PLC buffer with 0.5 mM sodium orthovanadate as a phosphorylation inhibitor. Phos-tag gel contains 50 µM phos-tag chemical compound in 12.5% acrylamide gel (SuperSep Phos-tag; Wako Chemicals USA). After electrophoresis, Phos-tag PAGE gels were rinsed in Tris-glycine blotting buffer containing 10 mM EDTA for 10 min three times, and twice without EDTA. For the control of protein dephosphorylation of LPAR2 cytoplasmic domain, we performed alkaline phosphatase treatment of protein lysates for 2 h at 50 C°. Antibodies for WB were as follows: monoclonal anti-Myc antibody (1:2,000 dilution, 9E10; Santa Cruz Biotechnology, Inc.), polyclonal rabbit anti-LPAR2 antibody (1: 375, HPA019616; ATLAS Antibodies), monoclonal anti–N-cadherin antiserum (1:200, MNCD2; Developmental Studies Hybridoma Bank), and monoclonal anti–α-tubulin antibody (1:5,000, B-5-1-2; Sigma-Aldrich).

### VinTS FRET

Samples for analysis of FRET by acceptor photobleaching were imaged using a laser scanning confocal microscope (LSM 510 Meta) and a 25× NA 1.4 Ph2 or 63× Plan Apochromat NA 1.4 Ph3 oil objective lens as specified. The CFP and YFP channels were excited using the 440-nm diode laser and the 514-nm argon lines, respectively. The two emission channels were split using a 545-nm dichroic mirror, which was followed by a 475–525-nm bandpass filter for CFP and a 530-nm longpass filter for YFP (Chroma Technology Corp). Pinholes were opened to give a depth of focus of 2 mm for each channel. Scanning was performed on a sequential line-by-line basis for each channel. The gain for each channel was set to ∼75% of dynamic range (12-bit, 4,096 gray levels) and offset such that backgrounds were zero. Time-lapse mode was used to collect one prebleach image for each channel followed by bleaching with a minimum of 50 iterations of the 514-nm argon laser line at maximum power (to bleach YFP). A second postbleach image was then collected for each channel. Control nonbleached areas were acquired for all samples in the same field of view as bleached cells to confirm specificity of FRET detection. Pre- and post-bleach CFP and YFP images were then imported into ImageJ for processing (AccpbFRET plugin; Roszik et al., 2008). Within this, all pixels were reregistered (using the Fast Hartley Transform algorithm) to correct for any shifts in the xy plane; bleed-through between donor and acceptor channels and efficiency of acceptor photobleaching were all accounted for in the analysis.

### Invasion assay in microfluidic tunnels

The closed compartment of the microchip was loaded with a solution of SDF-l (0.25 µg/ml), and NC explants were labeled with nuclear-mCherry and membrane-GFP, or nuclear and membrane-mCherry were seeded in front of each tunnel. The chip was sealed and imaged on an upright microscope (i80; Nikon) with a 20× lens.

### Online supplemental material

Fig. S1 shows that LPAR2 and Autotaxin are required for NC migration. Fig. S2 details the use of a microfluidic chamber to test 2D migration under physical constraints. Fig. S3 provides the results of Delaunay triangulation analyses of a cell dispersion assay. Video 1 shows how control and LPAR2MO-injected NC cells respond to Sdf-1 chemoattraction. Video 2 shows how NC cells migrate through various sizes of microfluidic tunnels. Video 3 shows a dispersion assay with controls and LPAR2MO cells. Video 4 monitors N-cadherin internalization using HaloTag pulse labeling. Video 5 shows that the impaired LPAR2MO cell dispersion is rescued by modification of cell–cell adhesions. Video 6 shows that functional blocking antibody of N-cadherin rescues LPAR2MO phenotype in microfluidic tunnel.

## Acknowledgments

We thank J. Leslie, A. Yap, and J.P. Thiery for comments on the manuscript. We thank Elena Scarpa for providing data that contributed to Fig. 5 J.

This study was supported by grants from the Medical Research Council (MR/J000655/1) and the Biotechnology and Biological Sciences Research Council to R. Mayor; by The Ministry of Education, Culture, Sports, Science and Technology (MEXT) Grants-in-Aid for Scientific Research (KAKENHI; 25111702), Japan Society for the Promotion of Science (JSPS) KAKENHI (24700962), the Uehara memorial foundation, and the Asahi Glass foundation to S. Kuriyama; and JSPS KAKENHI (25290042) to M. Tanaka. E. Theveneau was supported by the Wellcome Trust Value in People Award. M. Parsons was supported by a Royal Society University Research fellowship. G. Charras was supported by a Royal Society University Research Fellowship. A. Benedetto was supported by a Biotechnology and Biological Sciences Research Council grant (BB/F019769/1) to G. Charras.

The authors declare no competing financial interests.

Author contributions: S. Kuriyama and E. Theveneau performed the experiments. S. Kuriyama, E. Theveneau, and R. Mayor designed the experiments. S. Kuriyama, E. Theveneau, and R. Mayor analyzed and interpreted most of the data. A. Benedetto, E. Theveneau, and G. Charras designed the microfluidic chips. A. Benedetto performed the microfluidic chip fabrication. M. Parsons and E. Theveneau performed the FRET experiments. M. Tanaka and S. Kuriyama designed the fusion proteins. A. Kabla performed the correlation length analysis. E. Theveneau and R. Mayor organized the figures and supplementary information. E. Theveneau, R. Mayor, S. Kuriyama, and A. Kabla wrote the manuscript. All authors commented on the manuscript.

## References

References
Alfandari
,
D.
,
H.
Cousin
, and
M.
Marsden
.
2010
.
Mechanism of Xenopus cranial neural crest cell migration
.
4
:
553
560
.
Aybar
,
M.J.
,
M.A.
Nieto
, and
R.
Mayor
.
2003
.
Snail precedes slug in the genetic cascade required for the specification and migration of the Xenopus neural crest
.
Development.
130
:
483
494
.
Barriga
,
E.H.
,
P.H.
Maxwell
,
A.E.
Reyes
, and
R.
Mayor
.
2013
.
The hypoxia factor Hif-1α controls neural crest chemotaxis and epithelial to mesenchymal transition
.
J. Cell Biol.
201
:
759
776
.
,
A.
,
P.B.
Tran
,
D.
Ren
,
S.
Assimacopoulos
,
E.A.
Grove
, and
R.J.
Miller
.
2005
.
The chemokine stromal cell-derived factor-1 regulates the migration of sensory neuron progenitors
.
J. Neurosci.
25
:
3995
4003
.
,
A.
,
H.
Jung
,
D.
Ren
, and
R.J.
Miller
.
2009
.
The chemokine SDF-1/CXCL12 regulates the migration of melanocyte progenitors in mouse hair follicles
.
Differentiation.
77
:
395
411
.
Bowes
,
J.B.
,
K.A.
Snyder
,
E.
Segerdell
,
R.
Gibb
,
C.
Jarabek
,
E.
Noumen
,
N.
Pollet
, and
P.D.
Vize
.
2008
.
Xenbase: a Xenopus biology and genomics resource
.
Nucleic Acids Res.
36
:
D761
D767
.
Carmona-Fontaine
,
C.
,
H.K.
Matthews
,
S.
Kuriyama
,
M.
Moreno
,
G.A.
Dunn
,
M.
Parsons
,
C.D.
Stern
, and
R.
Mayor
.
2008
.
Contact inhibition of locomotion in vivo controls neural crest directional migration
.
Nature.
456
:
957
961
.
Carmona-Fontaine
,
C.
,
E.
Theveneau
,
A.
Tzekou
,
M.
,
M.
Woods
,
K.M.
Page
,
M.
Parsons
,
J.D.
Lambris
, and
R.
Mayor
.
2011
.
Complement fragment C3a controls mutual cell attraction during collective cell migration
.
Dev. Cell.
21
:
1026
1037
.
Clay
,
M.R.
, and
M.C.
Halloran
.
2013
.
Rho activation is apically restricted by Arhgap1 in neural crest cells and drives epithelial-to-mesenchymal transition
.
Development.
140
:
3198
3209
.
Contos
,
J.J.
,
N.
Fukushima
,
J.A.
Weiner
,
D.
Kaushal
, and
J.
Chun
.
2000a
.
Requirement for the lpA1 lysophosphatidic acid receptor gene in normal suckling behavior
.
97
:
13384
13389
.
Contos
,
J.J.
,
I.
Ishii
, and
J.
Chun
.
2000b
.
Lysophosphatidic acid receptors
.
Mol. Pharmacol.
58
:
1188
1196
.
Contos
,
J.J.
,
I.
Ishii
,
N.
Fukushima
,
M.A.
Kingsbury
,
X.
Ye
,
S.
Kawamura
,
J.H.
Brown
, and
J.
Chun
.
2002
.
Characterization of lpa(2) (Edg4) and lpa(1)/lpa(2) (Edg2/Edg4) lysophosphatidic acid receptor knockout mice: signaling deficits without obvious phenotypic abnormality attributable to lpa(2)
.
Mol. Cell. Biol.
22
:
6921
6929
.
Davies
,
P.C.
,
L.
Demetrius
, and
J.A.
Tuszynski
.
2011
.
Cancer as a dynamical phase transition
.
Theor. Biol. Med. Model.
8
:
30
.
DeSimone
,
D.W.
,
L.
Davidson
,
M.
Marsden
, and
D.
Alfandari
.
2005
.
The Xenopus embryo as a model system for studies of cell migration
.
Methods Mol. Biol.
294
:
235
245
.
Doxzen
,
K.
,
S.R.
Vedula
,
M.C.
Leong
,
H.
Hirata
,
N.S.
Gov
,
A.J.
Kabla
,
B.
, and
C.T.
Lim
.
2013
.
Guidance of collective cell migration by substrate geometry
.
Integr Biol (Camb).
5
:
1026
1035
.
Duband
,
J.L.
2010
.
Diversity in the molecular and cellular strategies of epithelium-to-mesenchyme transitions: Insights from the neural crest
.
4
:
458
482
.
Foty
,
R.A.
, and
M.S.
Steinberg
.
2005
.
The differential adhesion hypothesis: a direct evaluation
.
Dev. Biol.
278
:
255
263
.
Friedl
,
P.
,
J.
Locker
,
E.
Sahai
, and
J.E.
Segall
.
2012
.
Classifying collective cancer cell invasion
.
Nat. Cell Biol.
14
:
777
783
.
García
,
C.D.
, and
C.S.
Henry
.
2006
.
Micro-molding for poly(dimethylsiloxane) microchips
.
Methods Mol. Biol.
339
:
27
36
.
Grashoff
,
C.
,
B.D.
Hoffman
,
M.D.
Brenner
,
R.
Zhou
,
M.
Parsons
,
M.T.
Yang
,
M.A.
McLean
,
S.G.
Sligar
,
C.S.
Chen
,
T.
Ha
, and
M.A.
Schwartz
.
2010
.
Measuring mechanical tension across vinculin reveals regulation of focal adhesion dynamics
.
Nature.
466
:
263
266
.
Groysman
,
M.
,
I.
Shoval
, and
C.
Kalcheim
.
2008
.
A negative modulatory role for rho and rho-associated kinase signaling in delamination of neural crest cells
.
Neural Dev.
3
:
27
.
Hanahan
,
D.
, and
R.A.
Weinberg
.
2011
.
Hallmarks of cancer: the next generation
.
Cell.
144
:
646
674
.
Houben
,
A.J.
, and
W.H.
Moolenaar
.
2011
.
Autotaxin and LPA receptor signaling in cancer
.
Cancer Metastasis Rev.
30
:
557
565
.
Huang
,
R.Y.
,
C.C.
Wen
,
C.K.
Liao
,
S.H.
Wang
,
L.Y.
Chou
, and
J.C.
Wu
.
2012
.
Lysophosphatidic acid modulates the association of PTP1B with N-cadherin/catenin complex in SKOV3 ovarian cancer cells
.
Cell Biol. Int.
36
:
833
841
.
Ishiyama
,
N.
,
N.
Tanaka
,
K.
Abe
,
Y.J.
Yang
,
Y.M.
Abbas
,
M.
Umitsu
,
B.
Nagar
,
S.A.
Bueler
,
J.L.
Rubinstein
,
M.
Takeichi
, and
M.
Ikura
.
2013
.
An autoinhibited structure of α-catenin and its implications for vinculin recruitment to adherens junctions
.
J. Biol. Chem.
288
:
15913
15925
.
Kabla
,
A.J.
2012
.
Collective cell migration: leadership, invasion and segregation
.
J. R. Soc. Interface.
9
:
3268
3278
.
Kam
,
Y.
, and
V.
Quaranta
.
2009
.
Cadherin-bound β-catenin feeds into the Wnt pathway upon adherens junctions dissociation: evidence for an intersection between β-catenin pools
.
PLoS ONE.
4
:
e4580
.
Kasemeier-Kulesa
,
J.C.
,
R.
McLennan
,
M.H.
Romine
,
P.M.
Kulesa
, and
F.
Lefcort
.
2010
.
CXCR4 controls ventral migration of sympathetic precursor cells
.
J. Neurosci.
30
:
13078
13088
.
Kinoshita
,
E.
,
E.
Kinoshita-Kikuta
,
K.
Takiyama
, and
T.
Koike
.
2006
.
Phosphate-binding tag, a new tool to visualize phosphorylated proteins
.
Mol. Cell. Proteomics.
5
:
749
757
.
Kuriyama
,
S.
, and
R.
Mayor
.
2009
.
A role for Syndecan-4 in neural induction involving ERK- and PKC-dependent pathways
.
Development.
136
:
575
584
.
,
B.
, and
A.
Nicolas
.
2012
.
Physically based principles of cell adhesion mechanosensitivity in tissues
.
Rep. Prog. Phys.
75
:
116601
.
Lim
,
J.
, and
J.P.
Thiery
.
2012
.
Epithelial-mesenchymal transitions: insights from development
.
Development.
139
:
3471
3486
.
Lin
,
M.E.
,
D.R.
Herr
, and
J.
Chun
.
2010
.
Lysophosphatidic acid (LPA) receptors: signaling properties and disease relevance
.
Prostaglandins Other Lipid Mediat.
91
:
130
138
.
Liu
,
Y.
,
R.
Burkhalter
,
J.
Symowicz
,
K.
Chaffin
,
S.
Ellerbroek
, and
M.S.
Stack
.
2012
.
Lysophosphatidic Acid disrupts junctional integrity and epithelial cohesion in ovarian cancer cells
.
J. Oncol.
2012
:
501492
.
Los
,
G.V.
, and
K.
Wood
.
2007
.
The HaloTag: a novel technology for cell imaging and protein analysis
.
Methods Mol. Biol.
356
:
195
208
.
Mayor
,
R.
, and
E.
Theveneau
.
2013
.
The neural crest
.
Development.
140
:
2247
2251
.
Mills
,
G.B.
, and
W.H.
Moolenaar
.
2003
.
The emerging role of lysophosphatidic acid in cancer
.
Nat. Rev. Cancer.
3
:
582
591
.
Miyake
,
Y.
,
N.
Inoue
,
K.
Nishimura
,
N.
Kinoshita
,
H.
Hosoya
, and
S.
Yonemura
.
2006
.
Actomyosin tension is required for correct recruitment of adherens junction components and zonula occludens formation
.
Exp. Cell Res.
312
:
1637
1650
.
Monsoro-Burq
,
A.H.
,
R.B.
Fletcher
, and
R.M.
Harland
.
2003
.
Neural crest induction by paraxial mesoderm in Xenopus embryos requires FGF signals
.
Development.
130
:
3111
3124
.
Moolenaar
,
W.H.
,
L.A.
van Meeteren
, and
B.N.
Giepmans
.
2004
.
The ins and outs of lysophosphatidic acid signaling
.
Bioessays.
26
:
870
881
.
,
S.
,
Q.
Tao
,
N.R.
Menon
,
J.
Heasman
, and
C.
Wylie
.
2009
.
N- and E-cadherins in Xenopus are specifically required in the neural and non-neural ectoderm, respectively, for F-actin assembly and morphogenetic movements
.
Development.
136
:
1327
1338
.
Nguyen
,
D.X.
,
P.D.
Bos
, and
J.
Massagué
.
2009
.
Metastasis: from dissemination to organ-specific colonization
.
Nat. Rev. Cancer.
9
:
274
284
.
Ogata
,
S.
,
J.
Morokuma
,
T.
Hayata
,
G.
Kolle
,
C.
Niehrs
,
N.
Ueno
, and
K.W.
Cho
.
2007
.
.
Genes Dev.
21
:
1817
1831
.
Olesnicky Killian
,
E.C.
,
D.A.
Birkholz
, and
K.B.
Artinger
.
2009
.
A role for chemokine signaling in neural crest cell migration and craniofacial development
.
Dev. Biol.
333
:
161
172
.
Polyak
,
K.
, and
R.A.
Weinberg
.
2009
.
Transitions between epithelial and mesenchymal states: acquisition of malignant and stem cell traits
.
Nat. Rev. Cancer.
9
:
265
273
.
Powell
,
D.R.
,
A.J.
,
S.G.
Britt
, and
K.B.
Artinger
.
2013
.
Riding the crest of the wave: parallels between the neural crest and cancer in epithelial-to-mesenchymal transition and migration
.
Wiley Interdiscip Rev Syst Biol Med.
5
:
511
522
.
Rørth
,
P.
2009
.
Collective cell migration
.
Annu. Rev. Cell Dev. Biol.
25
:
407
429
.
Roszik
,
J.
,
J.
Szöllosi
, and
G.
Vereb
.
2008
.
AccPbFRET: an ImageJ plugin for semi-automatic, fully corrected analysis of acceptor photobleaching FRET images
.
BMC Bioinformatics.
9
:
346
.
Sasai
,
N.
,
K.
Mizuseki
, and
Y.
Sasai
.
2001
.
Requirement of FoxD3-class signaling for neural crest determination in Xenopus
.
Development.
128
:
2525
2536
.
Sauka-Spengler
,
T.
, and
M.
Bronner-Fraser
.
2008
.
A gene regulatory network orchestrates neural crest formation
.
Nat. Rev. Mol. Cell Biol.
9
:
557
568
.
Smicun
,
Y.
,
O.
Gil
,
K.
Devine
, and
D.A.
Fishman
.
2007
.
S1P and LPA have an attachment-dependent regulatory effect on invasion of epithelial ovarian cancer cells
.
Gynecol. Oncol.
107
:
298
309
.
Sumida
,
G.M.
,
T.M.
Tomita
,
W.
Shih
, and
S.
.
2011
.
Myosin II activity dependent and independent vinculin recruitment to the sites of E-cadherin-mediated cell-cell adhesion
.
BMC Cell Biol.
12
:
48
.
Taguchi
,
K.
,
T.
Ishiuchi
, and
M.
Takeichi
.
2011
.
Mechanosensitive EPLIN-dependent remodeling of adherens junctions regulates epithelial reshaping
.
J. Cell Biol.
194
:
643
656
.
Taillard
,
E.D.
,
P.
Waelti
, and
J.
Zuber
.
2008
.
Few statistical tests for proportions comparison
.
Eur. J. Oper. Res.
185
:
1336
1350
.
Tanaka
,
M.
,
S.
Kuriyama
, and
N.
Aiba
.
2012
.
Nm23-H1 regulates contact inhibition of locomotion, which is affected by ephrin-B1
.
J. Cell Sci.
125
:
4343
4353
.
Tao
,
Q.
,
S.
,
P.D.
McCrea
,
J.
Heasman
, and
C.
Wylie
.
2007
.
G-protein-coupled signals control cortical actin assembly by controlling cadherin expression in the early Xenopus embryo
.
Development.
134
:
2651
2661
.
Theveneau
,
E.
, and
R.
Mayor
.
2012
.
Neural crest delamination and migration: from epithelium-to-mesenchyme transition to collective cell migration
.
Dev. Biol.
366
:
34
54
.
Theveneau
,
E.
, and
R.
Mayor
.
2013
.
Collective cell migration of epithelial and mesenchymal cells
.
Cell. Mol. Life Sci.
70
:
3481
3492
.
Theveneau
,
E.
,
L.
Marchant
,
S.
Kuriyama
,
M.
Gull
,
B.
Moepps
,
M.
Parsons
, and
R.
Mayor
.
2010
.
Collective chemotaxis requires contact-dependent cell polarity
.
Dev. Cell.
19
:
39
53
.
Theveneau
,
E.
,
B.
Steventon
,
E.
Scarpa
,
S.
Garcia
,
X.
Trepat
,
A.
Streit
, and
R.
Mayor
.
2013
.
Chase-and-run between adjacent cell populations promotes directional collective migration
.
Nat. Cell Biol.
15
:
763
772
.
Thiery
,
J.P.
,
H.
Acloque
,
R.Y.
Huang
, and
M.A.
Nieto
.
2009
.
Epithelial-mesenchymal transitions in development and disease
.
Cell.
139
:
871
890
.
Umezu-Goto
,
M.
,
Y.
Kishi
,
A.
Taira
,
K.
Hama
,
N.
Dohmae
,
K.
Takio
,
T.
Yamori
,
G.B.
Mills
,
K.
Inoue
,
J.
Aoki
, and
H.
Arai
.
2002
.
Autotaxin has lysophospholipase D activity leading to tumor cell growth and motility by lysophosphatidic acid production
.
J. Cell Biol.
158
:
227
233
.
Urh
,
M.
, and
M.
Rosenberg
.
2012
.
HaloTag, a Platform Technology for Protein Analysis
.
Curr. Chem. Genomics.
6
:
72
78
.
Wheelock
,
M.J.
,
Y.
Shintani
,
M.
Maeda
,
Y.
Fukumoto
, and
K.R.
Johnson
.
2008
.
.
J. Cell Sci.
121
:
727
735
.
Wirtz
,
D.
,
K.
Konstantopoulos
, and
P.C.
Searson
.
2011
.
The physics of cancer: the role of physical interactions and mechanical forces in metastasis
.
Nat. Rev. Cancer.
11
:
512
522
.
Xiang
,
S.Y.
,
S.S.
Dusaban
, and
J.H.
Brown
.
2013
.
Lysophospholipid receptor activation of RhoA and lipid signaling pathways
.
Biochim. Biophys. Acta.
1831
:
213
222
.
Yonemura
,
S.
,
Y.
,
T.
Watanabe
,
A.
Nagafuchi
, and
M.
Shibata
.
2010
.
α-Catenin as a tension transducer that induces adherens junction development
.
Nat. Cell Biol.
12
:
533
542
.

Abbreviations used in this paper:

• ANOVA

analysis of variance

•
• CCM

collective cell migration

•
• EMT

epithelial–mesenchymal transition

•
• FDx

Fluorescein-dextran

•
• FRET

Förster resonance energy transfer

•
• LPA

lysophosphatidic acid

•
• MO

morpholino oligonucleotide

•
• NC

neural crest

•
• PDMS

polydimethyl methylhydrogen siloxane

•
• VinTS

vinculin-tension sensor

## Author notes

S. Kuriyama and E. Theveneau contributed equally to this paper.

E. Theveneau’s present address is Centre de Biologie du Developpement UMR5547, Centre National de la Recherche Scientifique, Université Paul Sabatier, Toulouse, France.