Early in the development of the central nervous system, progenitor cells undergo a shape change, called apical constriction, that triggers the neural plate to form a tubular structure. How apical constriction in the neural plate is controlled and how it contributes to tissue morphogenesis are not fully understood. In this study, we show that intracellular calcium ions (Ca2+) are required for Xenopus neural tube formation and that there are two types of Ca2+-concentration changes, a single-cell and a multicellular wave-like fluctuation, in the developing neural plate. Quantitative imaging analyses revealed that transient increases in Ca2+ concentration induced cortical F-actin remodeling, apical constriction and accelerations of the closing movement of the neural plate. We also show that extracellular ATP and N-cadherin (cdh2) participate in the Ca2+-induced apical constriction. Furthermore, our mathematical model suggests that the effect of Ca2+ fluctuations on tissue morphogenesis is independent of fluctuation frequency and that fluctuations affecting individual cells are more efficient than those at the multicellular level. We propose that distinct Ca2+ signaling patterns differentially modulate apical constriction for efficient epithelial folding and that this mechanism has a broad range of physiological outcomes.
The construction of multicellular biological systems requires substantial tissue movement based on the collective motility of single cells. Problems such as congenital malformations, breakdown of homeostasis and cancer can arise when the mechanisms that guide the concerted movement of single cells are dysfunctional (Guillot and Lecuit, 2013). In vertebrates, the central nervous system (CNS) is established during early embryogenesis by the collective closing movement of an epithelial sheet of neural progenitor cells in a process called neural tube closure (NTC), and NTC failure causes a common human birth defect (Copp et al., 2003; Suzuki et al., 2012). During NTC, neuroepithelial cells on the dorsal side contract and stabilize their apical sides to change from a columnar to a wedge shape. As a common strategy for folding an epithelial sheet, this cellular morphogenesis, called apical constriction (AC), generates physical force to bend the neural plate (NP) inward (Martin and Goldstein, 2014; Sawyer et al., 2010; Suzuki et al., 2012). The lateral edges of the bending neural plate lift above the dorsal surface, eventually come together along the dorsal midline and fuse to form the neural tube.
AC involves the accumulation of actin filaments (F-actin) at the apical side, where non-muscle myosin II generates contractile force. Thus, AC is thought to be controlled by a coordinated action of cytoskeletal networks, cell-cell adhesions and their regulators (Martin and Goldstein, 2014; Sawyer et al., 2010; Suzuki et al., 2012). However, how AC in the individual cell is controlled collectively within the tissue and how AC contributes to tissue remodeling are not fully understood. Recent studies have revealed that the AC in Drosophila and Caenorhabditis elegans epithelial invaginations is highly dynamic and is governed by oscillations of the apical actomyosin network, in which cell contraction is followed by stabilization (Martin and Goldstein, 2014; Mason et al., 2013; Sawyer et al., 2010). Hence, we hypothesized that the signaling pathway that induces rapid cellular behaviors controls AC during NTC.
In this study, we examined the roles of intracellular calcium ion (Ca2+) signaling in NTC. Generally, the intracellular Ca2+ concentration is maintained at a low level under basal conditions and is transiently increased by Ca2+ influx through membrane-localized Ca2+ channels or release from the endoplasmic reticulum via channels such as the inositol triphosphate receptor (IP3R) (Clapham, 2007). Ca2+ binds to and activates its target proteins that control diverse biological phenomena, such as cell proliferation, muscle contraction, gene expression and cellular movement (Clapham, 2007; Markova and Lenne, 2012). Inhibiting Ca2+ influx causes NTC to fail (Lee and Nagele, 1986; Smedley and Stanisstreet, 1986), while inducing an increase in intracellular Ca2+ promotes folding morphogenesis of the NP (Ferreira and Hilfer, 1993; Moran and Rice, 1976). Since these effects at the tissue level are accompanied by changes in cell shape (Ferreira and Hilfer, 1993), intracellular Ca2+ signaling has been proposed to regulate AC, although the detailed molecular mechanism is unclear.
Physiological significance of intracellular Ca2+ in NTC
We first examined the role of intracellular Ca2+ signaling during NTC in the African clawed frog, Xenopus laevis (Fig. S1) by using pharmacological treatments. Perturbing the Ca2+ influx with 2-aminoethoxydiphenyl borate (2APB, which blocks IP3R and other membrane-localized Ca2+ channels) or with nifedipine (which blocks voltage-dependent Ca2+ channels) inhibited NTC (Fig. 1A). We confirmed these results by measuring the width of NP using expression of the pan-neural marker Sox2 (Fig. 1B). Phalloidin staining revealed that AC in NP cells was inhibited, thus delaying folding morphogenesis of the NP (Fig. 1C,D). Since these inhibitors did not change the expression patterns of marker genes (Fig. 1E), these results suggest that Ca2+ signaling regulates AC during NTC independently of neural development.
Live-imaging analyses of Ca2+ fluctuation
In a number of organisms, including Xenopus, active changes of intracellular Ca2+ concentration have been observed during development (Hunter et al., 2014; Leclerc et al., 2000; Markova and Lenne, 2012; Shindo et al., 2010; Wallingford et al., 2001). Therefore, we hypothesized that Ca2+ signaling is spatially and temporally correlated with cellular morphogenesis. We investigated the patterns of intracellular Ca2+ concentration changes during NTC using live cell-imaging analyses with transiently expressed genetically encoded Ca2+ indicators for optical imaging (GECOs) (Fig. S2A) (Zhao et al., 2011). The intensity of R-GECO1.0 fluorescence fluctuated actively throughout the NP at both single-cell and whole-tissue levels (Fig. 2; Movie 1). In this study, we defined a single increase and the following decrease in GECO fluorescence intensity as a Ca2+ ‘transient’. We observed two spatial patterns of Ca2+ transients, one that occurred at the single-cell level and lasted less than 40 s (Fig. 2A,C) and another that originated from one or a few cells and propagated radially, wave-like, across neighboring cells until anywhere from several cells to hundreds of cells were involved (Fig. 2B; Movie 2). These multicellular Ca2+ fluctuations lasted up to 100 s (Fig. 2C). The fold changes in fluorescence intensity were larger at the multicellular level than those at the single-cell level (Fig. 2C). The multicellular waves spread at a speed of 3.3±0.94 µm/s (mean±s.d.), and larger waves spread at higher speeds (Pearson's r=0.499, P=0.0213, Fig. 2D), suggesting that a relationship between the wave propagation and the spatial strength of the Ca2+ signaling contributes to diverse patterns of the Ca2+ fluctuation.
Quantitative analyses of the Ca2+ fluctuation
To characterize the spatial and temporal patterns of these Ca2+ fluctuations more precisely, we conducted quantitative image analyses based on ratiometric imaging (Fig. S2B; Materials and Methods). After extracting and visualizing the Ca2+ fluctuation profiles, we confirmed that most of the Ca2+ transients occurred in the NP rather than in the non-neural ectoderm (Christodoulou and Skourides, 2015), which is lateral to the NP and did not undergo folding morphogenesis (Fig. 3A). We quantified the number of Ca2+ transients as a function of time and found that the number of single-cell transients particularly increased in the last 100 min before NTC was completed (Fig. 3B,C; Fig. S3A). Therefore, we divided the NTC process into an early and a late phase (Fig. 3B). The number of multicellular transients was one order of magnitude smaller than that of single-cell transients throughout the NTC process, and the number of transients remained constant through both the early and late phases (Fig. 3D). Taken together, these results suggest that single-cell Ca2+ transients were more frequent in the late phase than in the early phase, whereas the spatial scale (the propagative property) of the Ca2+ transients was larger in the early phase than in the late phase. The technical limitations of our method did not allow us to segment individual cell outlines, so we could not compare the total number of cells contributing to the Ca2+ transients. Therefore, we determined the total area of each type of transient and confirmed that the ratio value of the multicellular transients to the single-cell transients decreased from 0.55 to 0.18 in the late phase of NTC (Fig. 3E). These results suggest that the effect of Ca2+ fluctuations on NTC was derived from both single-cell and multicellular transients in the early phase, but mostly from single-cell transients in the late phase.
Ca2+ fluctuation precedes NTC
In our observations, the multicellular wave-like Ca2+ propagations were associated with a slight deformation of the neural plate tissue (Movies 1 and 2). Together with the observations that perturbation of Ca2+ influx inhibits NTC, it is suggested that intracellular Ca2+ signaling is essential for the movements that close the neural tube. We compared the temporal pattern of NTC movements with the total intensity values of R-GECO1.0 at tissue level. As expected, peaks in the R-GECO1.0 temporal profile overlapped with peaks in the speed of the closing movements, both in the early and late phase (Fig. 3F). Cross-correlation analyses revealed a maximum correlation when the R-GECO1.0 value was shifted 20 s later in time (Fig. 3G). Thus, the Ca2+ fluctuation slightly and consistently preceded the tissue movements to close the neural tube.
Regulatory mechanisms upstream of the Ca2+ fluctuation
We examined the effect of the pharmacological treatment on Ca2+ fluctuations. 2APB suppressed both the single-cell and multicellular transients in the early phase, whereas nifedipine did not clearly exhibit suppression effects on the multicellular transients (Fig. 3A-D; Fig. S3A). These results suggest that the IP3R pathway regulates both types of Ca2+ transient, whereas the voltage-dependent Ca2+ channel preferentially regulates Ca2+ transients at the single-cell level. Cross-correlation analyses revealed that 2APB completely abolished the temporal correlation between Ca2+ fluctuations and movements of the NP toward closing the neural tube compared with nifedipine (Fig. 3G). This difference might be due to residual multicellular transients in the nifedipine-treated embryo that induced a transient acceleration in the closing movement of the NP (Fig. 3D). Treatment with 2APB or nifedipine did not affect the basal intracellular Ca2+ level, as estimated using the ultrasensitive Ca2+ indicator yellow Cameleon-Nano (YC-Nano) (Horikawa et al., 2010) (Fig. S3B). This result suggests that, at the concentrations used, these inhibitors selectively affected the active Ca2+ fluctuations, and therefore, that the active Ca2+ signaling is required for NTC.
During Xenopus gastrulation, a purinergic receptor (Corriden and Insel, 2010; Schwiebert and Zsembery, 2003) for extracellular ATP (eATP) is required for frequent Ca2+ elevations (Shindo et al., 2010; Wallingford et al., 2001), and early ectodermal cells respond to eATP causing transient contractions of their apical side (Kim et al., 2014). To examine the possibility that eATP is involved in Ca2+ fluctuations in the NP, we exhausted the eATP by overexpressing Xenopus ectonucleoside triphosphate diphosphohydrolase 1 (E-NTPDase1) (Massé et al., 2007). E-NTPDase1 decreased the frequency of both single-cell and multicellular Ca2+ transients compared to an enzymatically inactive variant of E-NTPDase1 (ΔACR) (Fig. 4A-C; Fig. S3C, Movie 3). Moreover, overexpressing E-NTPDase1 abolished the temporal correlation between Ca2+ fluctuations and the closing movements of the NP, as seen with the 2APB treatment (Fig. 4D). We further examined the involvement of eATP in Ca2+ fluctuations using the fluorescence resonance energy transfer (FRET)-based ATP indicator ATeam3.10 (see Supplementary Materials and Methods) (Imamura et al., 2009). We found that extracellularly expressed ATeam3.10 increased the ratio of yellow-to-cyan fluorescence on the cell membrane compared with a mutant form of ATeam3.10 that is unable to bind to ATP (RK mutant) (Fig. 4E,F). Co-expressing E-NTPDase1 suppressed this increase in the ratio of yellow-to-cyan fluorescence but co-expressing the ΔACR variant did not (Fig. 4E,F). Taken together, these data suggest that eATP exists in the NP and regulates the Ca2+ fluctuations.
Ca2+ transients induce cortical F-actin remodeling and AC
Ca2+ fluctuations in the embryo regulate various cellular behaviors, such as cell flattening (Markova and Lenne, 2012), cell intercalation and polarization (Shindo et al., 2010; Wallingford et al., 2001), and cell contraction (Hunter et al., 2014). In addition, during wound healing in the developing epidermis, multicellular wave-like Ca2+ transients accelerate the closure process through cell shape changes (Antunes et al., 2013; Herrgen et al., 2014), suggesting that Ca2+-mediated cellular mechanisms underlie tissue development and homeostasis through epithelial remodeling. We investigated the relationship between Ca2+ fluctuations and NTC at the cellular level by simultaneously imaging R-GECO1.0 and the F-actin marker Lifeact-EGFP, which marks the cell periphery (Riedl et al., 2008). We found that immediately after a Ca2+ transient, mesh-like F-actin structures developed in the center of the cell and were maintained for several minutes, even after the Ca2+ had dropped back to basal levels within a few tens of seconds (Fig. 5A, +60 s to +160 s, and B; Movie 4). Furthermore, the area of the apical side of the NP cells decreased rapidly and then stabilized within a minute after a Ca2+ transient (Fig. 5A,B; Movie 4). Cross-correlation analyses among the R-GECO1.0 values, medial F-actin and cell area dynamics revealed that the maximum correlation with the F-actin signals was obtained at time lag 0, whereas with the cell area dynamics, it occurred when the R-GECO1.0 value was shifted later in time (Fig. 5C), which corresponded with tissue-level analysis (Fig. 3G). These data suggest that the rapid activation and prolonged F-actin remodeling induced by Ca2+ may be the mechanism of AC and NTC by which relaxation to the original state (i.e. before Ca2+ transients occurred) is suppressed. We could not determine whether the Ca2+-induced F-actin remodeling was due to overall polymerization or to the relocalization of F-actin to the apical side; both processes are implicated in the rapid cell contractions during Drosophila mesodermal invagination (Mason et al., 2013). In addition, it is also possible that non-muscle myosin II is activated by Ca2+ transients through known pathways such as calmodulin-dependent kinase (CaMKII)-mediated RhoA (Murakoshi et al., 2011) and myosin light chain kinase (MLCK) activation (Vicente-Manzanares et al., 2009).
We addressed the causal relationship between the Ca2+ transient and AC by using 1-(2-nitrophenyl)ethyl (NPE)-caged inositol triphosphate (IP3) and 1-(4,5-dimethoxy-2-nitrophenyl)ethyl (DMNPE)-caged ATP to experimentally induce Ca2+ transients. Brief ultraviolet (UV) illumination applied to an embryo that had been injected with NPE-caged IP3 induced a localized Ca2+ transient (Movie 5), followed by a decrease in the apical-side area of the UV-illuminated NP cells compared with that of unilluminated cells (Fig. 5D). Furthermore, by uncaging DMNPE-caged ATP applied to the ringer solution, we induced an intense multicellular Ca2+ propagation wave followed by the shrinkage of NP cells at the tissue level (Fig. 5E; Movies 6 and 7). Since F-actin in the medial region is linked to the junctional cadherin complexes during AC (Martin and Goldstein, 2014; Mason et al., 2013; Sawyer et al., 2010), we focused on N-cadherin (cdh2), which is required for F-actin assembly in the NP (Morita et al., 2010; Nandadasa et al., 2009). We found that blocking the function of N-cadherin in DMNPE-caged ATP-stimulated NP cells inhibited the decrease in apical area (Fig. 5E; Movies 8 and 9). These data suggest that the Ca2+ transients induce AC in NP cells, and that N-cadherin plays a role in the Ca2+-induced AC.
Mathematical analysis of NTC based on pulsed apical contraction
Our experimental results raise the possibility that the Ca2+ fluctuations function to decrease the total apical area of the NP. It is possible that Ca2+ fluctuations modulate cellular properties or mechanical processes. Indeed, we have demonstrated that Ca2+ fluctuations appear to regulate the localization of F-actin, which is involved in regulating both cell properties and mechanical processes. There are two types of transient Ca2+ fluctuations: those that occur at the single-cell level, and those that occur at the multicellular level. In order to evaluate differences in the effect and function of each type of Ca2+ fluctuation at the cellular level on tissue deformation, we constructed a multi-cell-based mechanical model, in which the Ca2+ fluctuations modulate a cellular mechanical parameter of a simplified epithelial tissue (Fig. 6A; Supplementary Materials and Methods). Our model is based on the vertex model, which is a well-known mathematical framework used to describe multicellular tissue dynamics (Farhadifar et al., 2007; Nagai and Honda, 2001; Okuda et al., 2013; Rauzi et al., 2008). To examine how the transient and patterned modulation of mechanical processes or properties in the cell results in the persistent shrinkage of an epithelial sheet, we introduced two components to our model. The first is the mechanical effect of the Ca2+ transient as a modification of the line tension of the cellular edges over a short time period (Fig. 6B), which is expected to decrease the cell surface area (see Supplementary Materials and Methods for definition). The second is the constrictive nature of the apical cell surface with a ratchet effect (Fig. 6C,D): the cells' natural surface area and the cells' natural perimeter are permitted to decrease but not to increase (see Supplementary Materials and Methods for definition). In our model, by regulating the spatial and temporal pattern of the modification of the line tensions, the effect of the Ca2+ transient at the single-cell and the multicellular levels can be simulated. Furthermore, the regulations of the cells' natural surface area and the cells' natural perimeter enabled us to introduce a ratchet-like mechanism (Martin and Goldstein, 2014; Mason et al., 2013; Sawyer et al., 2010) that has not been previously introduced (Spahn and Reuter, 2013).
Using this model, we first examined the combined effect of the Ca2+ transient (pulse) and the constrictive nature of apical cell surface (Acn) on the modeled epithelial tissue. Without either of these components and with the pulse component alone, the tissue size rarely changed at the end of the simulation time (Fig. 7A,B). Interestingly, adding the Acn component caused the modeled tissue to gradually decrease in size, and the combination of pulse and Acn components accelerated such decreasing behavior (Fig. 7A,B; Movie 10). At the single-cell level, transiently increased potential energies of apical surface area and perimeter, which were induced by modification of line tension during a single pulse, led to contraction of apical surface area and perimeter, respectively (Fig. 7C; Fig. S4). When the Acn component existed, those contractions further resulted in the irreversible decrease of the natural surface area and the natural perimeter (Fig. 7C; Fig. S4). Interestingly, after the contractions of apical cell surface area and perimeter that were induced by the pulse, a considerable relaxation event also occurred (Fig. 7C; Fig. S4), which were indeed observed during an AC event after a single Ca2+ transient in vivo (Fig. 5B). These results suggest that our model reproduces the in vivo situation, in which Ca2+ transients accelerate AC and, hence, decreases the tissue size when the stabilization mechanisms are present.
Next, in order to examine the role of Ca2+ transients in decreasing tissue size qualitatively, we gradually changed the pulse number and found that the effect became more pronounced as the total number of pulses was increased independently of the proportion of activated edges per cell (Fig. 7D). Under simulation with sparse pulses, the relationship between the number of pulses and the constriction rate showed a nearly linear dependence (Fig. S6). Next, we concentrated on the spatial and temporal frequencies of the pulse to examine the collective effect. In comparing dense versus sparse pulses, we found that dense pulses induced tissue deformation more rapidly within the active period (Fig. 7E,G). However, over the entire course of the simulation, dense pulses were less effective than sparse pulses when the same total number of pulses was used (Fig. 7E,F; Movie 11). The negative effect of pulse collectivity on the decrease in tissue size was attenuated as the proportion of activated edges per cell was decreased (Fig. 7F). However, because almost all of the segments contract in the in vivo situation (Fig. 5A), these data suggest that random Ca2+ fluctuations at the single-cell level decrease the tissue size more effectively than multicellular collective Ca2+ fluctuations.
In order to examine the robustness of the above results, we analyzed the modified models with varied assumptions (Fig. S5). When the coefficients of elasticity of a cell's apical surface and apical perimeter were assumed to decrease during the pulse, the effectiveness of sparse pulses on decreasing tissue size was observed in the presence of the Acn component (Fig. S5A), similar to when the pulse was linked to activation of line tension (Fig. 7E). On the other hand, when the tension fluctuation was induced during the pulse, the decrease in tissue size was not accelerated (Fig. S5B). As another attempt, we analyzed the two modified models in which the natural/preferred area was allowed to increase. The results showed that the effectiveness of sparse pulses on decreasing tissue size was also observed in the presence of the Acn component, although the differences in the constriction rate became smaller (Fig. S5C,D). These data suggest that under the conditions that recapitulate the combinational effects in vivo, the prediction regarding the effectiveness of sparse pulses on decreasing tissue size was robust.
In order to validate the prediction from the model simulations in vivo, we fitted a linear mixed model with single-cell and multicellular Ca2+ transients as explanatory variables and the amount of NTC as a response variable to data of untreated embryos, and obtained estimated contributions of single-cell Ca2+ transients and multicellular Ca2+ transients to the amount of NTC (Fig. S7 and Supplementary Materials and Methods). As is shown in Table 1, the coefficient of single-cell Ca2+ transients (0.072±0.011; estimated mean±s.e.m.) was significantly larger than that of multicellular Ca2+ transients (0.028±0.011) (P-value=7.231×10−3, Student's t-test). These results support the prediction from the model simulations that the Ca2+ fluctuations at the single-cell level more effectively accelerate AC and NTC than multicellular Ca2+ fluctuations do in vivo.
In this study, we show that two types of active Ca2+ signaling patterns modulate AC via cortical F-actin remodeling, and thereby contribute differently to efficient NTC. Our present study and a recent paper (Christodoulou and Skourides, 2015) propose a new role of active Ca2+ fluctuation in accelerating AC. The trigger for Ca2+ transients and the mechanism that allow it to propagate through a tissue remain to be investigated. We did not detect any bias in the direction of wave propagation relative to the body axes, whereas the wave propagation in the wound healing response has previously been shown to originate from wounded cells (Antunes et al., 2013; Herrgen et al., 2014). Multicellular Ca2+ transients in the NP produced a higher intensity of R-GECO1.0 fluorescence per cell than did single-cell level transients (Fig. 2); therefore, the wave-like propagation of Ca2+ transients at the tissue level might be induced above a particular threshold Ca2+ signal strength. Our study also suggests that eATP regulates both the single-cell and the multicellular Ca2+ transients (Fig. 4). eATP activates downstream signaling in the cell that secreted the ATP (autocrine signaling) and in neighboring cells (paracrine signaling) (Corriden and Insel, 2010; Schwiebert and Zsembery, 2003). Several pathways, including ATP-permeable release channels and ATP-filled vesicles have been proposed to secrete cytosolic ATP, which is triggered by mechanical stimulation or extracellular biochemical cues (Corriden and Insel, 2010; Schwiebert and Zsembery, 2003). Hence, ATP release induced by spontaneous mechanical stretch or compression might trigger Ca2+ transients and possibly the wave propagation in the developing NP, in which dynamic autonomous and non-autonomous tissue movements take place (Morita et al., 2012; Wallingford and Harland, 2001). In order to dissect more precisely the regulatory mechanisms upstream of Ca2+ transients, functional studies using genome editing, for example, would be useful since pharmacological approaches can be misleading due to off-target effects (Bootman et al., 2002; Xu et al., 2005).
Although multicellular Ca2+ transients have been shown to accelerate closure process during wound healing (Antunes et al., 2013; Herrgen et al., 2014), the contribution of random single-cell Ca2+ fluctuations to epithelial remodeling has been elusive, probably because the effect of these transients is too small at the tissue level to be examined experimentally. Moreover, since the temporal patterns of single-cell Ca2+ fluctuations varied among embryos (Fig. 3B-E), it might be regulated by mechanisms that are inherently noisy or heterogenous, further complicating efforts to draw general conclusions about its functional significance. In this study, we used a cell-vertex mathematical model in surface view to investigate how distinct Ca2+ fluctuation patterns affect AC and persistent epithelial remodeling (Figs 6 and 7). Computational simulations using this model showed that pulsed contractions accelerate the decrease in tissue size, even when pulses occurred with a sparse frequency, suggesting that Ca2+-induced AC always contributes positively to NTC. Our mathematical and statistical analyses also suggest that spatially and temporally random Ca2+ transients at the single-cell level function to reduce the overall tissue size more effectively than multicellular collective Ca2+ fluctuations, making them cost-effective for normal NTC. On the other hand, the multicellular wave-propagated Ca2+ transients might be useful when a tissue must be closed quickly, as in wound healing. Since the number of single-cell Ca2+ transients increases in the late phase of NTC, a genetic program that represses collectivity and promotes the randomization of Ca2+ fluctuations might exist to ensure that embryos are able to undertake their primitive CNS formation under physiological conditions. In summary, the Ca2+-dependent mechanisms highlighted in this study might be acting in parallel with known pathways that control AC, such as the Shroom3 or Wnt/PCP pathway, as a compensatory system to avoid incomplete NTC, which is the second most frequent congenital malformation in humans (Copp et al., 2003).
MATERIALS AND METHODS
All protocols for animal care and experiments were approved by the Institutional Animal Care and Use Committee of the National Institutes of Natural Sciences, Japan, and were performed in accordance with institutional guidelines for the care and use of laboratory animals.
Injections, morpholinos and RNA expression constructs
Experiments with Xenopus laevis embryos were performed as previously described (Suzuki et al., 2010). After injection, the embryos were cultured in 3% Ficoll with 0.1× Steinberg's solution until stage 9, then washed and cultured in 0.3× Marc's Modified Ringer's medium (MMR) until the appropriate stage for analysis (Nieuwkoop and Faber, 1967).
Antisense morpholino (Gene Tools) was injected into embryos at the four-cell stage (17 ng per embryo). The morpholino sequences were as follows: N-cadherin morpholino 5′-GAAGGGCTCTTTCCGGCACATGGTG-3′ (Nandadasa et al., 2009); 5mis-N-cadherin morpholino, 5′-GAACGGGTCTTTGCGCCACATCGTG-3′.
Capped sense RNAs were synthesized with the mMESSAGE mMACHINE kit (Ambion), purified on a NICK column (GE Healthcare) and injected into embryos at the four-cell stage. All open reading frames were cloned into the pCS2p+ vector. Embryos received injections of RNA constructs in the following amounts: R-GECO1.0 (Zhao et al., 2011) and G-GECO1.2 (Zhao et al., 2011), 1 ng per embryo; YC-Nano50 (Horikawa et al., 2010), YC-Nano30 (Horikawa et al., 2010) and YC-Nano15 (Horikawa et al., 2010), 500 pg per embryo. Cell surface-targeted ATeam3.10 (AT3.10-GPI) and its RK mutant were generated by fusing ATeam3.10 (Imamura et al., 2009) or its RK mutant (Imamura et al., 2009) with the signal peptide sequence and the GPI-anchor domain of Xenopus glypican4 (Ohkawara et al., 2003), and 1 ng was injected per embryo. EGFP was injected at 100 pg per embryo. Membrane-targeted EGFP and mRFP were generated by fusing the farnesylation signal of c-Ha-Ras to the C terminus of EGFP and mRFP, respectively, and 100-200 pg was injected per embryo. Lifeact-EGFP (Riedl et al., 2008) was injected at 180 pg per embryo. Xenopus E-NTPDase1 and its mutant version, in which the enzymatic active site had been deleted (ΔACR), were injected at 100-200 pg per embryo.
For inhibitor treatments, 100× stocks of 2APB (D9754; Sigma) and nifedipine (145-05781; Wako) in DMSO were added to the medium at stage 13, to concentrations of 25 µM and 200 µM, respectively. For uncaging experiments, NPE-caged IP3 (I23580; Molecular Probes) was injected into embryos at the four-cell stage (5 pmol per embryo). DMNPE-caged ATP (A1049; Molecular Probes) was added to the medium to a concentration of 100 µM.
In situ hybridization and F-actin staining
In situ hybridization was performed as previously described (Suzuki et al., 2010). The following plasmids were used for probe synthesis: Sox2 (XL039o24; XDB3); N-cadherin (XL289n05ex; XDB3); and Epidermal keratin (XL056e18; XDB3). For F-actin visualization, embryos were fixed in MOPS, EGTA, magnesium sulfate, formaldehyde buffer (MEMFA) and stained with Alexa 546-phalloidin (A22283; Molecular Probes).
Devitellinized embryos were mounted on an agarose-coated glass-based dish (3910-035; AGC, Japan) filled with 0.3× MMR. Confocal fluorescence images were acquired as a z-series of 512×512 pixels using a spinning-disc confocal unit (CSU-X1; Yokogawa, Japan) with an EMCCD camera (iXon3; Andor) and 455-, 488- and 561-nm lasers (Andor) on an inverted microscope (IX81; Olympus) with a 10× (UPlanSAPO 10×/0.4; Olympus) or 20× (UPlanAPO 20×/0.7; Olympus) objective lens at a room temperature adjusted to 18-20°C. Unless otherwise noted, we observed the R-GECO1.0 fluorescence at 20 or 40 s intervals for up to 6 h.
For UV illumination in uncaging experiments, mercury lamp-based epifluorescent optics with a band-pass filter (AT350/50×; Chroma) were inserted into the path of the excitation light, and the duration of illumination was controlled by iQ2 software (Andor). For the photolysis of NPE-caged IP3, the samples were illuminated for 3 s by a UV light reduced with a 6% neutral-density filter (32-ND6; Olympus). For the photolysis of DMNPE-caged ATP, the samples were illuminated for 20 s by unreduced UV light.
Image processing and analyses were performed using ImageJ (NIH) and R (www.r-project.org) software. All quantifications used the maximum intensity projection of the z-series. Details are provided in Supplementary Materials and Methods.
The mathematical model constructed in this study is based on the vertex model (Farhadifar et al., 2007; Nagai and Honda, 2001; Okuda et al., 2013; Rauzi et al., 2008). The simulations were programmed in C. Details are provided in Supplementary Materials and Methods.
Illustrations and statistical analyses
Figures were assembled with Photoshop (Adobe), Illustrator (Adobe) and PowerPoint (Microsoft). Data plots and statistical analyses were performed in Excel (Microsoft) and R software. Data sets were subjected to Shapiro–Wilk's test to determine distributions. No statistical method was used to predetermine sample size. The experiments were not randomized. The investigators were not blinded to allocation during experiments and outcome assessment.
We thank the Data Integration and Analysis Facility, and the Spectrography and Bioimaging Facility in the National Institute for Basic Biology Core Research Facilities for technical support; members of the N.U. laboratory, S. Nonaka laboratory and N. Shiina laboratory for valuable discussions and comments.
M. Suzuki, M. Sato, H.K., Y.H., T.F. and N.U. designed research; M. Suzuki, H.K., Y.H. and N.Y. performed research; M. Suzuki, M. Sato., H.K., K.H., H.I., T.N. and R.E.C. contributed new reagents and analytical tools; M. Suzuki., M. Sato, H.K., Y.H. and N.Y. analyzed the data; and M. Suzuki, M. Sato, H.K., R.E.C. and N.U. wrote the paper.
This work was supported by KAKENHI [23770261 and 15KT0154 to M.S., 22127007 and 15H05865 to N.U.] from the Japan Society for the Promotion of Science (JSPS).
The authors declare no competing or financial interests.