Controlling contractile instabilities in the actomyosin cortex
Abstract
The actomyosin cell cortex is an active contractile material for driving cell- and tissue morphogenesis. The cortex has a tendency to form a pattern of myosin foci, which is a signature of potentially unstable behavior. How a system that is prone to such instabilities can rveliably drive morphogenesis remains an outstanding question. Here, we report that in the Caenorhabditis elegans zygote, feedback between active RhoA and myosin induces a contractile instability in the cortex. We discover that an independent RhoA pacemaking oscillator controls this instability, generating a pulsatory pattern of myosin foci and preventing the collapse of cortical material into a few dynamic contracting regions. Our work reveals how contractile instabilities that are natural to occur in mechanically active media can be biochemically controlled to robustly drive morphogenetic events.
https://doi.org/10.7554/eLife.19595.001Introduction
Alan Turing described in his seminal 1952 paper the ability of an initially homogeneous spatial system that contains diffusing and chemically interacting species to form a self-organized pattern (Turing, 1952). Turing’s original conjecture was that such processes contribute to the patterning of developing organisms. While many examples have been found that are compatible with this idea (Kondo and Asal, 1995; Müller et al., 2012; Sheth et al., 2012; Raspopovic et al., 2014), self-organized patterning in morphogenesis, however, is known to not only rely on biochemical regulation but also depend on cell-and tissue scale active mechanical processes (Turing, 1952; Howard et al., 2011). General physical mechanisms by which the interplay between regulatory and mechanical processes endows active biological materials to form self-organized spatiotemporal patterns have remained largely unexplored.
Actomyosin contractility (Bray and White, 1988; Salbreux et al., 2012) is an essential cellular mechanical process, responsible for driving many cell- and tissue scale morphogenetic events (Murrell et al., 2015; Levayer and Lecuit, 2012). The cortex consists to a large extent of actin filaments and myosin motor proteins, forming a thin layer underneath the cell membrane that can be thought of as a thin film of an active gel (Salbreux et al., 2012; Jülicher et al., 2007). Contractility by myosin motor proteins generates active tension in the gel, and gradients in active tension are known to generate cortical flows of this layer (Mayer et al., 2010). Cortical flow participates in forming the cytokinetic furrow (Bray and White, 1988; Benink et al., 2000; Yumura, 2001; Eggert et al., 2006), and drives polarization of the one-cell stage C. elegans embryo (Hird and White, 1993; Guo and Kemphues, 1996; Cheeks et al., 2004; Munro et al., 2004; Goehring et al., 2011). Highly contractile cortices, like the one driving polarization in C. elegans, tend to exhibit transient accumulations of myosin that form a pulsatile pattern. Pulsatile actomyosin patterns are ubiquitous in development (Munro et al., 2004; Martin et al., 2009; Solon et al., 2009; Rauzi et al., 2010; Roh-Johnson et al., 2012; Maître et al., 2015), and it has been suggested that they result from positive feedback and contractile instabilities (Kruse and Jülicher, 2000; Bois et al., 2011; Gowrishankar et al., 2012; Kumar et al., 2014; Munjal et al., 2015; Hannezo et al., 2015). A contractile instability causes the cortex to become inhomogeneous over space, with cortical material collapsing into contracting regions (Bois et al., 2011; Alvarado et al., 2013). Theoretical work has shown that contractile instabilities are inevitable when contractility is high enough (Bois et al., 2011), raising the question of how a system that is prone to such instabilities can reliably drive morphogenesis. Here we show that there indeed exists a contractile instability in the actomyosin cortical layer of the C. elegans zygote, and we discover that this instability is controlled by a RhoA oscillator.
Results and discussion
In order to investigate spatiotemporal patterns in the C. elegans cortex, we first sought to see whether non-muscle myosin II (NMY-2) in the C. elegans zygote displays pulsatile dynamics (Munro et al., 2004). For this, we determined the temporal derivative of the average NMY-2 intensity (Figure 1A, averaging over a region in the posterior indicated by a white box), as a proxy of myosin foci assembly and disassembly behavior (Figure 1A–C). We also quantified the time-dependence of the average speed of cortical flow in this region as determined by Particle Image Velocimetry (PIV, see Appendix for detail) (Figure 1B,C). Notably, both quantities exhibited signs of oscillatory behavior (Figure 1C) and an auto-correlation analysis revealed periodic changes in both quantities with a time constant of approximately 30 s (Figure 1D,E). To conclude, the myosin foci pattern in the C. elegans zygote exhibits pulsatile, oscillatory dynamics.
We next sought to understand where this oscillatory behavior comes from. One possibility is that positive feedback mediated by RhoA (RHO-1 in C. elegans) (Bement et al., 2015), a key activator of myosin (Jenkins et al., 2006; Motegi and Sugimoto, 2006; Schonegg and Hyman, 2006), plays a role in generating this pulsatile pattern (Munjal et al., 2015). We investigated the dynamics of active RhoA by use of a GFP fused anillin homology domain (AHPH) probe, to image the GTP-bound active form of RhoA (Tse et al., 2012). We find that active RhoA forms a dynamic, pulsatile pattern that is similar to that of myosin, with both active RhoA and myosin co-localizing in pulsatile foci (Figure 1F, Figure 1—figure supplement 1, and Video 1). We speculated that flow-based transport of an activator of myosin could give rise to positive feedback and a contractile instability (Bois et al., 2011), favouring the spontaneous formation of self-organized pulsatory patterns (Munjal et al., 2015). However, testing for this possibility requires knowledge of the kinetics of active RhoA mediated myosin recruitment coupled with a hydrodynamic description of active cortical mechanics, for evaluating if the full mechanochemically coupled system indeed is unstable.
We set out to test if coupling RhoA mediated myosin recruitment to gel flow and advection results in an intrinsically unstable cortex (Figure 2A). To this end, we sought to determine the effective reaction kinetics of the regulatory interaction between active RhoA and myosin in vivo. We developed a method of measuring the kinetic diagram of active RhoA mediated myosin recruitment (CO-moving Mass Balance Imaging; COMBI): We investigated the mass balance of both species in the comoving frame of reference of the flowing cortex, under consideration of the effects of dilution/enrichment by divergent/convergent gel flow (Figure 1G) (Vallotton et al., 2004). In the frame of reference of the embryo, concentrations of myosin and active RhoA can change due to transport by flow (advective fluxes) or due to association/dissociation (chemical fluxes). The chemical fluxes and , where r denote active RhoA and m denotes myosin, correspond to reaction terms that capture turnover and biochemically regulated recruitment effects. They can depend on the concentrations of both species.
COMBI determines the average changes of per area myosin concentration () and of active RhoA concentration ( due to turnover/regulation and as a function of the concentrations of both species. This provides us with information of the reaction kinetics of RhoA and myosin in the myosin and active RhoA concentration phase space (Figure 1G). We determined and every 5 s by spinning-disk confocal microscopy (Materials and methods). Advective fluxes account for the effects of dilution/enrichment by divergent/convergent gel flow, and these were determined by measuring the velocity field of cortical flow by particle image velocimetry (PIV), an image-based crosscorrelation analysis (Mayer et al., 2010; Raffel et al., 2007) that quantifies the movement of interrogation areas between two sequential timelapse images (Materials and methods). The spatial resolution of the velocity field is determined by the spacing of the interrogation areas which we choose as . This is sufficiently smaller than the correlation length of cortical flow (the hydrodynamic length is ∼14 μm [Mayer et al., 2010; Saha et al., 2016]), hence, COMBI can provide information on actomyosin homeostasis by determining the average reaction kinetics on a timescale of seconds and a length scale of microns.
We visualize the reaction terms determined by COMBI in a vector field that illustrates the average evolution of concentrations of both species (Figure 1H). This reveals interesting features, for example a trajectory that starts with high active RhoA and low myosin levels will overshoot in its level of myosin prior to approaching the single stable fixed point (thick black line in Figure 1—figure supplement 2A). Key aspects of these dynamics can be captured by a nullcline analysis (Izhikevich, 2007). Reaction terms are zero on a nullcline ( for active RhoA, Figure 1H, solid green line; and for myosin, Figure 1H, solid magenta line), which describes the concentration that a species would achieve when the concentration of the other species is fixed. This reveals that active RhoA recruits myosin (Figure 1H, solid magenta line) (Motegi and Sugimoto, 2006; Schonegg and Hyman, 2006) while RhoA activation kinetics is essentially independent of myosin levels (Figure 1H, solid green line). The full kinetic landscape can be linearized over its entire range (Figure 1H, dashed lines; Figure 1—figure supplement 2B–D), capturing global aspects of RhoA-based myosin recruitment. To conclude, COMBI can provide insight into the cortical kinetics over a broad range of myosin and active RhoA concentrations.
Given our kinetic analysis, we next sought to test if the full mechanochemically coupled system is indeed unstable. We describe the actomyosin cortex as a thin film of an active gel (Simha and Ramaswamy, 2002; Kruse et al., 2004; Ahmadi et al., 2006; Salbreux et al., 2009), with active tension generation by myosin under control of RhoA (Figure 2A; see Appendix for details). We measured the relevant material parameters of the gel in vivo directly from laser ablation experiments (hydrodynamic length: , and a conversion factor from NMY-2 intensity to active tension, )(Saha et al., 2016). This allowed us to perform a linear stability analysis of the homogeneous state for the full model of the mechanochemical patterning system, with the above determined and linearized reaction kinetics between active RhoA and NMY-2 (Figure 1H). Figure 2B shows the corresponding stability diagram as a function of both the hydrodynamic length of the cortex and the active tension measure (see Appendix for detail). Notably, the homogeneous state, in which all concentrations are constant in space, always becomes unstable above a critical value of the active tension. Furthermore, we find that the parameter values of the C. elegans cortex are such that the system is close to the transition line between stable and unstable, but placed within the unstable regime (Figure 2B). Hence, our analysis is consistent with the actomyosin cortex in C. elegans being unstable and poised to form a spatial pattern.
Our theory predicts that the contractile instability depends on the strength of positive feedback, and thus the amount of recruitment of myosin by active RhoA (Figure 2C,D). Hence, we asked if suppression of RhoA mediated recruitment of myosin in the C. elegans zygote prevents the instability and results in a homogeneous NMY-2 distribution. LET-502 is the Rho-associated protein kinase that phosphorylates the regulatory myosin light chain, MLC-4, to activate NMY-2 (Piekny and Mains, 2002). Hence, reducing the concentration of LET-502 by RNAi should suppress RhoA mediated recruitment of myosin to the cortex. Indeed, COMBI analysis of let-502 RNAi embryos (30 hr) revealed that RhoA mediated recruitment of NMY-2 to the cortex is reduced, since the myosin nullcline displays a significantly decreased slope as compared to the non-RNAi condition (Figure 2E,F; see Table 1). Using the non-RNAi values of and and the linearized reaction kinetics between active RhoA and NMY-2 measured by COMBI in let-502 RNAi (dark dashed lines in Figure 2F), we find that the cortex is predicted to be stable because all eigenvalues are negative (Figure 2G, compare to D; see Appendix for detail). Consistent with this prediction, we observed that let-502 RNAi embryos display a homogeneous NMY-2 distribution without pulsatory myosin foci (Figure 2H, compare to Figure 1A; Video 2). We conclude that, consistent with COMBI and theory, the actomyosin cortex can be brought into a stable regime by reducing positive feedback via suppressing myosin recruitment by active RhoA.
We next asked if the patterns that are formed in the unstable regime in our theory correspond to the pattern of myosin foci observed in the embryo. Earlier work that considers a cortical gel with a diffusible activator of myosin (Bois et al., 2011; Kumar et al., 2014) suggests that a contractile instability results in a myosin foci pattern with a spacing that is determined by the hydrodynamic length . Indeed, a numerical solution of the full mechanochemical patterning system (Figure 2A) reveals the formation of a few dynamic contracting regions which travel and are spaced about apart (Figure 3—figure supplement 1, Video 3). These traveling peaks have rapid flows converging upon them (peak flow speed: ), and they persist and are not pulsatile. This pattern is different from the myosin foci pattern observed in the C. elegans zygote, which is pulsatile and exhibits a shorter spacing between foci (, compare to ; see Figure 3—figure supplement 1) (Munro et al., 2004; Mayer et al., 2010). This suggests that our model is missing an essential feature, which is responsible for determining the myosin pattern beyond the contractile instability.
To identify the element missing in our model, we note that our theory predicts that reducing myosin recruitment by let-502 RNAi should cause both myosin and RhoA to be homogeneous and non-pulsatile (all eigenvalues are negative, see Figure 2G). However, imaging active RhoA under let-502 RNAi revealed that while the pulsatile myosin pattern is lost, the pulsatile active RhoA pattern is still present (Figure 3A,C; Video 4). Similarly, we find that 16 hr of RNAi of nmy-2 led to an almost complete loss of cortical myosin with, however, active RhoA still forming pulsatile foci (Figure 3B,D; Video 5). We conclude that, in contrast to the scenario in Drosophila germband extension (Munjal et al., 2015), active RhoA in C. elegans exhibits pulsatile foci dynamics independently of NMY-2 function. Importantly, both the characteristic spacing of the myosin-independent active RhoA pattern and its characteristic timescale were similar between nmy-2 RNAi, let-502 RNAi, and the non-RNAi condition (Figure 3E,F). Given that active RhoA in the wild-type acts to recruit myosin (Figure 1H), this raises the possibility that the myosin-independent dynamic active RhoA pattern is responsible for setting the myosin spatiotemporal pattern beyond the contractile instability. We conclude that the dynamic active RhoA pattern is generated in a manner that is independent of the myosin foci pattern, possibly through an independent RhoA spatiotemporal oscillator.
Oscillatory activities of Rho GTPases have previously been observed (Hwang et al., 2005; Miller and Bement, 2009; Das et al., 2012; Antoine-Bertrand et al., 2016). We next asked if this spatiotemporal oscillator requires ect-2, a RhoGEF, responsible in the early morphogenesis (Motegi and Sugimoto, 2006; Schonegg and Hyman, 2006). Indeed, RNAi of ect-2 leads to a complete absence of RhoA pulsation (Video 6). Furthermore, it is interesting to speculate if the myosin-independent active RhoA oscillator that we identify here is related to the RhoA/actin-based excitable oscillatory system reported previously (Bement et al., 2015; Westendorf et al., 2013). To test if the underlying mechanism to generate myosin-independent active RhoA oscillator is shared between C. elegans single-cell embryo and Xenopus embryo, we used COMBI to investigate the effective kinetic regulation between active RhoA and actin. We used LifeAct tagRFP-T as a probe for filamentous actin in the cortex (Riedl et al., 2008; Reymann et al., 2016). We determined the kinetic diagram in the active RhoA and actin concentration phase plane, to quantitatively evaluate the rate constants in the effective kinetic equations (Figure 1—figure supplement 3A). We find that the active RhoA nullcline is nearly vertical and inconsistent with actin behaving as a negative regulator of RhoA. Note that this does not exclude the general possibility of negative feedback between actin and RhoA (Robin et al., 2016), but suggests that the C. elegans cortex is normally operating in a regime where no such negative feedback is accessed. While the detailed mechanism as well as the kinetic interactions that underlie RhoA pulsation in C. elegans remain to be determined, the RhoGEF ect-2 is involved and the system appears to undergo spatiotemporal oscillations in the absence of negative feedback between actin and RhoA.
We next sought to test in our theory if it is possible that an active RhoA spatiotemporal oscillator sets the myosin pattern beyond the contractile instability (Figure 4A, left). To this end, we described the dynamical behavior of an active RhoA pacemaker by use of a generic model of spatiotemporal oscillating patterns, the complex Swift-Hohenberg Equation (Figure 4—figure supplement 1A) (Sakaguchi, 1997). Importantly, coupling in our model this generic spatiotemporal oscillator ( characteristic timescale, characteristic length scale, Figure 3E,F; see Appendix for detail) to the full mechanochemical patterning system does not destroy the active RhoA spatiotemporal oscillator pattern. Instead our model predicts that the myosin pattern (which in the absence of the generic oscillator formed a single traveling peak, see Figure 3—figure supplement 1) now follows that of the active RhoA spatiotemporal oscillator (Figure 4B left). Hence, the active RhoA oscillator can determine the myosin pattern in the unstable regime (Video 7). As a consequence, controlling the myosin pattern also results in reduced cortical flow speeds (peak flow speed: ) as compared to the case where the RhoA oscillator is absent (, see above). However, we find that the ability of the RhoA oscillator to control the myosin pattern critically depends on the level of mechanochemical feedback. We demonstrate this by reducing the hydrodynamic length in our model, which increases overall flow speeds and advection, and thereby increases the mechanochemical feedback strength. We find in our model that this change destroys the pattern of the active RhoA spatiotemporal oscillator. Both the myosin and active RhoA pattern no longer form a regular spatiotemporal oscillation (Figure 4—figure supplement 1, is reduced by to ). Instead, the system displays a dynamical state that is characterized by an irregular spatiotemporal pattern of dynamic contracting regions that move rapidly (Figure 4B right; Figure 4—figure supplement 1; Video 8). In this state, the pattern of active RhoA now depends on myosin and flows and is essentially under control of the contractile instability. Finally, flow speeds are again increased and comparable to the case when the RhoA oscillator is absent (peak flow speed: ). In conclusion, theory indicates that the active RhoA oscillator can act as a pacemaker for the system, to control the contractile instability and to prevent the formation of large and irregularly moving contracting regions of myosin.
We next sought to seek experimental evidence that the myosin pattern in the C. elegans zygote is under control of the RhoA pacemaker. To this end, we tested if increasing the level of mechanical feedback in C. elegans destroys the pattern of the active RhoA spatiotemporal oscillator as predicted from theory. For this we recorded space-time patterns of myosin in a midplane section under conditions of spd-5 RNAi. SPD-5 is a centriole constituent that is essential for centriole maturation (Hamill et al., 2002), and its RNAi leads to a delay of polarizing flows which gives us more time for an analysis of the pulsatory dynamics. To increase mechanochemical feedback we recorded space-time patterns of myosin under RNAi of the actin nucleator pfn-1 (Severson et al., 2002) for which the cortex is weakened and flow speeds are increased by a factor of three to five (Figure 4F) and for which the hydrodynamic length is decreased by to (Severson et al., 2002; Naganathan, unpublished). For spd-5 RNAi, we observed the ’normal’ pulsating pattern of myosin foci (Figure 4C left; Video 9 and Video 10). In contrast, for pfn-1 RNAi we observed large contracting regions of myosin that rapidly move in an irregular fashion (Figure 4C right, D right; Videos 9 and 11). Importantly, in pfn-1 RNAi active RhoA assembles in large and irregularly moving foci structures (Figure 4E right, compare to Figure 4E left and Figure 3A,C; Video 11). This suggests that the normal pattern of the active RhoA spatiotemporal oscillator is destroyed, and the distribution of active RhoA is now governed by the irregular dynamics of myosin. Note that pfn-1 RNAi does not destroy the general ability of RhoA to generate a pulsating pacemaker pattern, as revealed by double RNAi of pfn-1 and nmy-2 (Video 12). We conclude that reducing the hydrodynamic length increases mechanochemical feedback and advection. This causes the RhoA pacemaker to lose the ability to control the contractile instability. Consistent with the predictions of our theory, this destroys the RhoA pacemaker pattern and causes the system to undergo a transition to irregular behavior with large and rapidly moving contracting regions (Figure 4C–E). Interestingly, the uncontrolled gel is no longer capable to drive coherent flows of the cortex over large distances (Figure 4G), and the embryo fails to polarize (Severson et al., 2002). Taken together, our quantitative analysis is consistent with the interpretation that the spatiotemporal RhoA oscillator acts as a pacemaker in C. elegans, controlling the contractile instability of the actomyosin cortex.
We have here investigated the mechanisms of pattern formation in an active system that combines the contractile force generation and flow with regulation and advection. For this, we introduced the COMBI method to directly infer reaction kinetics without relying for example on photobleaching (Sprague et al., 2004). We determined the effective reaction kinetics of myosin and active RhoA in the actomyosin cortex with COMBI. This allowed us to build a quantitative model of mechanochemical patterning in the actomyosin layer. By use of linear stability analysis, we found that the actomyosin cortex is unstable and spontaneously forms a self-organized pattern. We speculate that during embryogenesis cells need high cortical contractility to drive morphological changes. This can lead them near or beyond contractile instabilities, leading to dynamics characterized by strong fluctuations and irregular behavior, possibly exhibiting active turbulence (Giomi, 2015). We suggest that such instabilities are inevitable in dynamic systems that are highly contractile. We discovered a spatiotemporal RhoA oscillator that determines the myosin pattern even beyond the contractile instability, thereby controlling the contractile instability. The independent biochemical RhoA oscillator endows the cell with the ability to use an intrinsically unstable active contractile medium for driving morphogenetic processes such as polarization. To conclude, our work paves the way for understanding pattern formation in active biological materials that utilize potentially unstable contractile processes.
Materials and methods
Worm strains, maintenance, and sample preparation
Request a detailed protocolThe following transgenic lines were used in this study. SWG003: nmy-2(cp8[nmy-2::GFP + unc-119(+)]) I; unc-119(ed3) III; gesIs002[Ppie-1::Lifeact::tagRFP-T::pie-1UTR + unc-119(+)], for imaging of GFP labelled NMY-2 (the images shown in Figure 1A and B and in Figure 4C). SWG012: nmy-2(ges6[nmy-2::tagRFP-T + unc-119(+)]) I; xsSi5[cb-UNC-119 (+) GFP:: ANI-1 (AH+PH)] II; unc-119(ed3) III, for imaging of tagRFP-T labelled NMY-2 and GFP labelled AHPH for a probe of active RhoA in the cortex (the images shown in Figure 1F, in Figure 2H, and in Figure 3A and B).
Worm strains were maintained at 20°C, and shifted to 24°C for 24 hr before imaging. Embryos were dissected in M9 buffer and mounted onto agar pads (2% agarose in water) to squish the embryos gently. All experiments were performed at 23–24°C. RNA interference experiments were performed by feeding as described in Naganathan et al. (2014). Feeding times for RNAi experiments were 16–18 hr for nmy-2, 23–25 hr for spd-5, 19–21 hr for pfn-1, 19–21 hr for pfn-1;nmy-2 double knockdown and 29–31 hr for let-502. Feeding clones were obtained from the Hyman lab (MPI-CBG, Dresden, Germany).
Imaging
Request a detailed protocolOne-cell stage embryos were observed under the inverted fluorescence microscope (Axio Observer Z1, Zeiss) using a Zeiss C-Apochromat 63 water immersion lens, equipped with a spinning disc confocal unit (Yokogawa, CSU-X1) and AOTF laser combiner (Andor, ALC). Fluorescence images were acquired by a sCMOS camera (Hamamatsu, ORCA flash 4.0) at 5 s time intervals for non-RNAi, let-502 RNAi, nmy-2 RNAi, and pfn-1;nmy-2 RNAi embryos. For pfn-1 RNAi embryos, images were taken every 3 s. Pixel size was 0.105 0.105 , all devices were controlled through manager(Edelstein et al., 2014). Fluorescence images of GFP and tagRFT-T labeled proteins in the embryos were excited by 488 and 561 nm lasers, respectively.
Image analysis
Request a detailed protocolPrior to COMBI analysis, images were filtered using the nonlocal means method (Buades et al., 2005), reducing spatially uncorrelated noise while preserving finer structures. Filtering was performed by averaging fluorescence intensities on the basis of the similarity between the fluorescence intensity profile in the interrogation area and the intensity profile in the neighboring region, i.e., a searching window. We set the size of the interrogation area and the searching window, to be pixels and pixels, corresponding to 0.525 0.525 and 2.625 2.625 , respectively. A filtering parameter, , was set to be for NMY-2 images, for AHPH images, where denotes the standard deviation of the fluorescence intensity in each image. We performed the filtering using a freely available code from MATLAB central (Fast Non-Local Means 1D, 2D Color and 3D by Kroon). Note that it is important to remove spatially uncorrelated noise prior to the computation of the spatial derivatives, since differential value is affected by spatially uncorrelated noises.
To perform COMBI, we first determined the cortical flow velocity, , by Particle Image Velocimetry (PIV) using a freely available PIV algorithm, PIVlab 1.32 (available from http://pivlab.blogspot.de/). PIV was performed on NMY-2 images by setting the interrogation area as 24 pixels with a step of 12 pixels. Velocity vectors were then interpolated to single pixel resolution for determining and and for both the background subtracted, active RhoA and myosin intensities (denoted by and , respectively). Intensity background levels were obtained by averaging the intensity in the region outside the embryo. and were then determined for each pixel throughout the cortical plane, by the use of the mass balance equations given in Figure 1G. We obtained a kinetic diagram (Figure 1H) by averaging and in boxes located in the anterior region. We determined average values for each embryo by averaging over the first 36 frames after the start of polarizing flow. For the non-RNAi case (Figure 1H), we report the average kinetic diagram from embryos, for let-502 RNAi (Figure 2F) we averaged over embryos. Note that the active RhoA and NMY-2 concentrations were normalized by the respective mean intensities of active RhoA and NMY-2 under non-RNAi conditions.
Correlation analysis
Request a detailed protocolTo characterize the myosin intensity change, and the cortical flow speed in a box of size in the posterior region, we determined the spatial average over the box. We then computed the autocorrelation function
where or , denote the standard deviation of , and denotes the mean of , averaged over time , where represents an average over time. is mean-subtracted and normalized by the variance of . The period of oscillation was determined by the peak position in the autocorrelation function. For a precise detection of oscillatory behavior, we removed from our analysis embryos in which the second peak in the autocorrelation function of the time course was undetectable (12 out of 25 cases for non-RNAi embryos, 5 out of 12 cases in let-502 embryos, 4 out of 12 embryos for nmy-2 RNAi embryos, respectively).
We obtained the characteristic length of the spatial pattern of myosin and RhoA by detecting the location of the first peak in the radial spatial intensity correlation function. The spatial intensity autocorrelation function of intensity, , was obtained by
where denotes the spatial average of , see Figure 3—figure supplement 2C. This function was radially averaged, and the first peak was detected (Figure 3—figure supplement 2D) in each time point. The radii of first peak were then averaged over time in each embryo. For the determination of the spacing in AHPH foci, the contrast of the fluorescence images was enhanced using the Contrast Limited Adaptive Histogram Equalization method using Matlab (Mathworks)(see Figure 3—figure supplement 2).
To characterize the spatial coherence of the velocity field, we evaluated the spatial correlation of the normalized velocity vectors, . The spatial correlation function was computed by,
where represents scalar product, and denotes the spatial average. Note that the coordinate transformation from Cartesian to polar coordinates of the orientation vectors, , provides a simpler representation of as,
where denotes the anti-clockwise angle from x-axis of . The above expression shows that the provides spatial correlation of the cosine similarity between and . Therefore, the characteristic length of the decay of represents the loss of correlation between the directions of velocity vectors, away. Larger characteristic length of the decay demonstrates the large-scale flow of the cortex.
For visualizing purpose, we transformed the coordinate system from to polar and then determined the average over the angle, , to plot as a function of the radius, , e.g. Figure 4G.
Appendix
Active gel description of the cortex
We described the dynamics of the cortical layer on the surface of the zygote by use of a full mechanochemically coupled system containing an active gel description for the cortical layer(Bois et al., 2011; Kumar et al., 2014) together with a reaction-diffusion-advection equation system for active RhoA and NMY-2. For the latter, the concentrations of active RhoA and NMY-2 in one dimension are given by and , respectively, where denotes position and denotes time. The dynamics of the 1-D surface concentration vector is given by
where , with denoting NMY-2 and denoting active RhoA. Here, represents the gel velocity field determined by active cortical mechanics, see below. and denote the respective fluxes of NMY-2 and active RhoAand from the cytosol onto the surface via turnover and recruitment. and denote the respective diffusion coefficients of active RhoA and NMY-2. We estimated the surface diffusion coefficients of both active RhoA and myosin, by analyzing FRAP recovery in ANI-1::GFP and NMY-2::GFP by use of a method considers of both turnover and lateral diffusion (Goehring et al., 2010). This revealed in both cases that lateral (surface) diffusion was below the detection limit and undistinguishable from zero. Hence, we expect surface diffusion of myosin and active RhoA to not significantly impact the dynamics of the cortical layer. However, for the purpose of preventing sharp peaks in concentration fields in our numerical simulations, we here set . Note that the general results of our numerical analysis do not change when increasing or decreasing either diffusion constant by a factor of 10. Importantly, the fastest growing mode occurs for wavelengths of the order of in all cases (Figure 2—figure supplement 1).
We describe forces and flows within the cortical layer in the framework of active gel theory (Bois et al., 2011; Kumar et al., 2014; Kruse et al., 2005; Mayer et al., 2010). We consider the cortex to be a thin film active viscous fluid in 1D, where NMY-2 generates active tension (Bois et al., 2011; Kumar et al., 2014; Mayer et al., 2010). Hence, total tension
in the layer is given by a sum of viscous tension arising characterized by a bulk viscosity , and active tension (or contractility) that depends on myosin concentration acording to , with and increasing function of the concentration and a coefficient that determines the magnitude of active tension.
The constitutive equation for an active viscous fluid is then given by Mayer et al. (2010)
We choose for a saturating Hill-function, (Bois et al., 2011; Kumar et al., 2014), to limit the active stress in our numerical simulations. Here, denotes the stationary concentration of NMY-2. Note that choosing a linear dependence without saturation does not significantly change the stability diagram, see Figure 2—figure supplement 2.
The force balance equation for an active viscous fluid in the presence of friction with a coefficient between the cortex and its surrounding cytoplasm and cell membrane (Mayer et al., 2010; Bois et al., 2011; Kumar et al., 2014) is given by
Equations 6 and 7 provide an equation of motion
where represents the hydrodynamic length of the cortex, which sets a correlation length for the velocity field. We determined the parameters and that characterize the active viscous fluid with a method that compares the relaxation dynamics of the cortex in response to cortical laser ablation (COLA) between experiment and theoretical predictions (, ) (Saha et al., 2016). Note that we ignored the polar or nematic order of actin filament in the cortex, which can introduce the tension anisotropy. Based on the recent study by Reymann et al. (Reymann et al., 2016), the alignment is mainly due to large-scale cortical flow in one-cell stage C. elegans embryo, which has the highest compressible component and nematic order parameter in the mid-zone of posterior.
Rescaling
We rescale the spatial coordinate with respect to , and rewrite Equations 5–7 according to
with the following rescaled quantities
We performed numerical integrations of Equation 9 with periodic boundary conditions by use of the pseudospectral method (Boyd, 2001).
Linear stability analysis
The homogeneous state with the concentrations and a vanishing flow field can be a stationary state of the system. To test the stability of this stationary state, we apply a small perturbation of , where and denote spatial frequency and eigenvalue, respectively. Inserting into Equation 9 and retaining linear terms only, the linear stability matrix with respect to the perturbation becomes
where . To determine the linearized matrix that characterizes biochemical regulation of NMY-2 and active RhoA, we performed a linear regression over the whole range of and values determined by COMBI (Figure 1H in main text) using
Here, and denote the normalized intensities of AHPH and NMY-2, respectively, with . We get
We determined the matrix directly from our experimental data (Figure 1H in main text) by robust regression and using standard deviations for weights, and find
Please see Figure 1—figure supplement 2 for a comparison between the linearized and full landscapes. The corresponding eigenvalues of are . Both of these are negative and real, indicating that the fixpoint is a sink.
The stability diagram shown in Figure 2B in the main text was obtained in the plane of hydrodynamic length and the active tension (or contractility) measure , where denotes the stationary concentration of in the homogeneous state.
Pulsatile dynamics of active RhoA
Phenomenological description of the pulsatory dynamics
To investigate if an active RhoA pacemaking oscillator can control the contractile instability and determine the myosin pattern, we describe active RhoA pulsation dynamics by an oscillating standing wave. For this, we utilize a generic spatiotemporal oscillator, the complex Swift-Hohenberg equation (Cross and Hohenberg, 1993; Sakaguchi, 1997). We introduce a complex variable , the real part of which represents the concentration of active RhoA, . We write a conservation law for according to
with
Here, is the pacemaking oscillator term using the complex Swift-Hohenberg equation, while and are phenomenological coefficients that determine the spatiotemporal dynamics of .
We chose a set of parameter values that is based on (Sakaguchi, 1997) and that gives rise to an oscillating pattern that is similar to the active RhoA oscillation shown in Figure 3 of the manuscript, with a characteristic length-scale of and a characteristic time-scale of 60 s ( i,e. the half-period corresponds to the active RhoA oscillation period). We find that the occurrence of a transition from a pacemaker- entrained state (i.e. the myosin pattern is determined by the pattern of the active RhoA pacemaker) to a detrained state (i.e. the active RhoA pattern is destroyed and now determined by cortical mechanics and myosin) is not sensitive to the detailed characteristics of the pacemaking oscillator, e.g. its amplitude, period, and spacing. For large hydrodynamic lengths , the myosin pattern is entrained to the active RhoA pulsatility. As we decrease , the period of myosin pattern oscillation increases, due to increasing amounts of accumulation of myosin and active RhoA at the myosin foci by advection. At intermediate values of around , we observed a transition from the entrained state to the detrained state. Here, the dynamical pattern is characterized by irregular movements and ‘mechanochemical turbulence’ (Figure 4 in main text, Figure 4—figure supplement 1), and depends on a both flow and advection as well as the active RhoA pacemaking activity.
Associated time-scales
To provide an intuitive understanding of entrainment-detrainment transition, we compare characteristic time-scales of advection and biochemical regulations. First we obtain an advection time-scale, , by solving Equation 8 for , where is Dirac’s delta function. For the infinite interval, , and a no-flux boundary condition, , the velocity profile, , is given by
hence, the velocity decay over space is characterized by the hydrodynamic length , while the magnitude of velocity is set by . We obtain a time-scale for advection as
which specifies the time required for the thin film active fluid to advect biochemical regulators over the characteristic length .
Next, we introduce a characteristic time-scale of biochemical regulation as . If , the RhoA pacemaking oscillator controls the dynamics of the full mechanochemical patterning system, and advection plays a minor role. In this case, NMY-2 dynamics are entrained to the pulsatory dynamics of the active RhoApacemaker. In contrast to this, flow and advection dominate the dynamics of the full mechanochemical patterning system when . Here, the regularly oscillating pacemaker pattern of active RhoAis destroyed due to increased advection, leading to irregular spatiotemporal dynamics and ‘mechanochemical turbulence’.
In the zygote, we find that for the non-RNAi case that is larger than . Thus the myosin pattern is entrained to the active RhoA pulsatory dynamics. On the other hand, we find that for the pfn-1 RNAi case that , which is smaller than . Hence, we observe irregular behavior (Figure 4B, right, in the main text).
References
-
Analysis of cortical flow models in vivoMolecular Biology of the Cell 11:2553–2563.https://doi.org/10.1091/mbc.11.8.2553
-
Pattern formation in active fluidsPhysical Review Letters 106:028103.https://doi.org/10.1103/PhysRevLett.106.028103
-
ConferenceA non-local algorithm for image denoisingIEEE Computer Society Conference on Computer Vision and Pattern Recognition. pp. 60–65.https://doi.org/10.1109/CVPR.2005.38
-
Pattern formation outside of equilibriumReviews of Modern Physics 65:851–1112.https://doi.org/10.1103/RevModPhys.65.851
-
Advanced methods of microscope control using μmanager softwareJournal of Biological Methods 1:e10.https://doi.org/10.14440/jbm.2014.36
-
Animal cytokinesis: from parts list to mechanismsAnnual Review of Biochemistry 75:543–566.https://doi.org/10.1146/annurev.biochem.74.082803.133425
-
Geometry and topology of turbulence in active nematicsPhysical Review X 5:031003.https://doi.org/10.1103/PhysRevX.5.031003
-
Cortical and cytoplasmic flow polarity in early embryonic cells of Caenorhabditis elegansThe Journal of Cell Biology 121:1343–1355.https://doi.org/10.1083/jcb.121.6.1343
-
Turing's next steps: the mechanochemical basis of morphogenesisNature Reviews Molecular Cell Biology 12:400–406.https://doi.org/10.1038/nrm3120
-
Oscillatory ROP GTPase activation leads the oscillatory polarized growth of pollen tubesMolecular Biology of the Cell 16:5385–5399.https://doi.org/10.1091/mbc.E05-05-0409
-
Active behavior of the cytoskeletonPhysics Reports 449:3–28.https://doi.org/10.1016/j.physrep.2007.02.018
-
Asters, vortices, and rotating spirals in active gels of polar filamentsPhysical Review Letters 92:078101.https://doi.org/10.1103/PhysRevLett.92.078101
-
Generic theory of active polar gels: a paradigm for cytoskeletal dynamicsThe European Physical Journal E 16:5–16.https://doi.org/10.1140/epje/e2005-00002-5
-
Actively contracting bundles of polar filamentsPhysical Review Letters 85:1778–1781.https://doi.org/10.1103/PhysRevLett.85.1778
-
Pulsatory patterns in active fluidsPhysical Review Letters 112:208101.https://doi.org/10.1103/PhysRevLett.112.208101
-
Biomechanical regulation of contractility: spatial control and dynamicsTrends in Cell Biology 22:61–81.https://doi.org/10.1016/j.tcb.2011.10.001
-
Pulsatile cell-autonomous contractility drives compaction in the mouse embryoNature Cell Biology 17:849–855.https://doi.org/10.1038/ncb3185
-
Regulation of cytokinesis by Rho GTPase fluxNature Cell Biology 11:71–77.https://doi.org/10.1038/ncb1814
-
Forcing cells into shape: the mechanics of actomyosin contractilityNature Reviews Molecular Cell Biology 16:486–498.https://doi.org/10.1038/nrm4012
-
Rho-binding kinase (LET-502) and myosin phosphatase (MEL-11) regulate cytokinesis in the early Caenorhabditis elegans embryoJournal of Cell Science 115:2271–2282.
-
Lifeact: a versatile marker to visualize F-actinNature Methods 5:605–607.https://doi.org/10.1038/nmeth.1220
-
Determining physical properties of the cell cortexBiophysical Journal 110:1421–1429.https://doi.org/10.1016/j.bpj.2016.02.013
-
Standing wave patterns for the complex Swift-Hohenberg equationProgress of Theoretical Physics 98:577–585.https://doi.org/10.1143/PTP.98.577
-
Actin cortex mechanics and cellular morphogenesisTrends in Cell Biology 22:536–545.https://doi.org/10.1016/j.tcb.2012.07.001
-
Hydrodynamics of cellular cortical flows and the formation of contractile ringsPhysical Review Letters 103:058102.https://doi.org/10.1103/PhysRevLett.103.058102
-
Analysis of binding reactions by fluorescence recovery after photobleachingBiophysical Journal 86:3473–3495.https://doi.org/10.1529/biophysj.103.026765
-
The chemical basis of morphogenesisPhilosophical Transactions of the Royal Society B: Biological Sciences 237:37–72.https://doi.org/10.1098/rstb.1952.0012
-
Myosin II dynamics and cortical flow during contractile ring formation in Dictyostelium cellsThe Journal of Cell Biology 154:137–146.https://doi.org/10.1083/jcb.200011013
Article and author information
Author details
Funding
Deutsche Forschungsgemeinschaft (SPP 1782)
- Stephan W Grill
European Research Council (281903)
- Stephan W Grill
Human Frontier Science Program (RGP0023/2014)
- Stephan W Grill
European Commission (ITN grant - 281903)
- Stephan W Grill
Max-Planck-Gesellschaft
- Stephan W Grill
European Commission (ITN grant - 641639)
- Stephan W Grill
Deutsche Forschungsgemeinschaft (GSC 97)
- Stephan W Grill
Deutsche Forschungsgemeinschaft (GR 3271/2)
- Stephan W Grill
Deutsche Forschungsgemeinschaft (GR 3271/3)
- Stephan W Grill
Deutsche Forschungsgemeinschaft (GR 3271/4)
- Stephan W Grill
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We are grateful to Michael Glotzer and Daniel Dickinson for providing C. elegans strains. We acknowledge Stefan Eimer for the gift of tagRFP-T gene. We thank G Salbreux, JS Bois, PGross and P Gönczy for fruitful discussions. SWG was supported by the DFG (SPP 1782, GSC 97, GR 3271/2, GR 3271/3, GR 3271/4), the European Research Council (grant No 281903), ITN grants 281903 and 641639 from the EU, the Max-Planck-Society as a Max-Planck-Fellow, and the Human Frontier Science Program (RGP0023/2014).
Copyright
© 2017, Nishikawa et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,476
- views
-
- 1,047
- downloads
-
- 86
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Physics of Living Systems
Many proteins have been recently shown to undergo a process of phase separation that leads to the formation of biomolecular condensates. Intriguingly, it has been observed that some of these proteins form dense droplets of sizeable dimensions already below the critical concentration, which is the concentration at which phase separation occurs. To understand this phenomenon, which is not readily compatible with classical nucleation theory, we investigated the properties of the droplet size distributions as a function of protein concentration. We found that these distributions can be described by a scale-invariant log-normal function with an average that increases progressively as the concentration approaches the critical concentration from below. The results of this scaling analysis suggest the existence of a universal behaviour independent of the sequences and structures of the proteins undergoing phase separation. While we refrain from proposing a theoretical model here, we suggest that any model of protein phase separation should predict the scaling exponents that we reported here from the fitting of experimental measurements of droplet size distributions. Furthermore, based on these observations, we show that it is possible to use the scale invariance to estimate the critical concentration for protein phase separation.