Pulsed actomyosin contractility underlies diverse modes of tissue morphogenesis, but the underlying mechanisms remain poorly understood. Here, we combined quantitative imaging with genetic perturbations to identify a core mechanism for pulsed contractility in early Caenorhabditis elegans embryos. We show that pulsed accumulation of actomyosin is governed by local control of assembly and disassembly downstream of RhoA. Pulsed activation and inactivation of RhoA precede, respectively, the accumulation and disappearance of actomyosin and persist in the absence of Myosin II. We find that fast (likely indirect) autoactivation of RhoA drives pulse initiation, while delayed, F-actin–dependent accumulation of the RhoA GTPase-activating proteins RGA-3/4 provides negative feedback to terminate each pulse. A mathematical model, constrained by our data, suggests that this combination of feedbacks is tuned to generate locally excitable RhoA dynamics. We propose that excitable RhoA dynamics are a common driver for pulsed contractility that can be tuned or coupled differently to actomyosin dynamics to produce a diversity of morphogenetic outcomes.
Pulsed contractility is a widespread mode of actomyosin contractility expressed by many non-muscle cells in which transient accumulations of F-actin and Myosin II accompany local contractions of the cell surface. Pulsed contractions were first identified in the polarizing Caenorhabditis elegans zygote (Munro et al., 2004) and have now been documented in a wide variety of embryonic and extra-embryonic epithelia (Martin et al., 2009; Solon et al., 2009; Blanchard et al., 2010; David et al., 2010; He et al., 2010; Rauzi et al., 2010) and mesenchymal cells (Kim and Davidson, 2011). A similar phenomenon, known as cell shape oscillations, has been observed in many cultured cells (Salbreux et al., 2007; Kapustina et al., 2008; Sedzinski et al., 2011). Pulsed contractions produce transient shape changes that can be biased or rectified in different ways to produce distinct morphogenetic outcomes, such as tissue invagination (Martin et al., 2009), tissue elongation (He et al., 2010; Rauzi et al., 2010; Levayer and Lecuit, 2013), epithelial tissue closure (Solon et al., 2009), and wound healing (Razzell et al., 2014). During embryonic development, pulsed contractions may represent an adaptation to accommodate rapid cell and tissue deformations while maintaining overall tissue integrity (Vasquez et al., 2014). In other contexts, such as in many cultured cells, shape oscillations may represent an aberrant behavior that manifests when cells lose normal adhesion to their substrates (Paluch et al., 2005; Salbreux et al., 2007), when microtubules are depolymerized (Bornens et al., 1989; Werner et al., 2007; Kapustina et al., 2008, 2013; Piekny and Glotzer, 2008; Rankin and Wordeman, 2010), or when contractile tension is very high during cytokinesis (Sedzinski et al., 2011).
Despite the widespread occurrence of pulsed contractions and increasing evidence for their functional relevance, the mechanisms that initiate and terminate pulsed contractions remain poorly understood. From a dynamic perspective, pulsed contractions represent a form of excitable behavior, exemplified by action potentials in neuronal cells (Izhikevich, 2007) or pulses of intracellular calcium release observed in many cell types (Goldbeter, 1996). Theoretical studies highlight two key ingredients for excitability: positive feedback to drive a rapid upswing in activity and delayed negative feedback to bring it back down again. A key challenge is to identify the specific modes of positive and negative feedback that drive pulsed contractions.
Multiple forms of positive feedback could contribute to initiating pulsed contractions. For example, local actomyosin-based contraction could promote further accumulation of actomyosin through mechanosensitive motor-filament binding (Effler et al., 2006; Fernandez-Gonzalez et al., 2009; Ren et al., 2009; He et al., 2010; Schiffhauer et al., 2016), by enhancing actin filament assembly and/or stability (Hayakawa et al., 2011; De La Cruz and Gardel, 2015) or by transporting and concentrating actomyosin and/or its upstream activators (White and Borisy, 1983; Dierkes et al., 2014; Munjal et al., 2015). Alternatively, dynamic clustering of F-actin and/or Myosin II by scaffolding proteins such as Anillin could promote Myosin II recruitment and focal contraction (Maddox et al., 2005). Finally, autocatalytic activation of upstream regulators such as RhoA could drive local excitation independent of, or in addition to, myosin-based tension or network contraction (Bement et al., 2015; Munjal et al., 2015; Zhang and Glotzer, 2015; Graessl et al., 2017). Similarly, multiple forms of delayed negative feedback could contribute to terminating pulses, including progressive buildup of steric or elastic resistance to further contraction (Dierkes et al., 2014), contraction-mediated disassembly, or delayed recruitment of disassembly factors or inhibitors of Myosin II or RhoA (Kasza and Zallen, 2011; Bement et al., 2015; Munjal et al., 2015; Graessl et al., 2017).
Here, we combined quantitative imaging with experimental manipulations and mathematical modeling to identify the dynamic basis for pulsed contractility in the early C. elegans embryo. Using single-molecule imaging and particle-tracking analysis, we provide definite evidence that the initiation of pulsed contractions does not involve or require local redistribution of actomyosin or its upstream activators. Instead, pulsed contractions are driven by local pulses of RhoA activity, which feed forward to control local accumulation of downstream targets F-actin, Myosin II, and Anillin. We present evidence that pulsed accumulation of RhoA is governed by locally excitable RhoA dynamics: Local autocatalytic activation of RhoA drives the rapid upswing of RhoA activity during pulse initiation, while F-actin–dependent recruitment of the redundantly acting RhoA GTPase-activating proteins (GAPs) RGA-3/4 provide delayed negative feedback to terminate the pulse. A minimal model, sharply constrained by our experimental data, suggests that this combination of feedbacks is sufficient to generate locally excitable or oscillatory RhoA dynamics and to explain quantitatively the temporal dynamics of RhoA activation and RGA-3/4 accumulation during each pulse. We propose that excitable RhoA dynamics define a core mechanism for pulsed contractility and suggest that this mechanism may be tuned or filtered through downstream effectors to control the size or spacing or lifetime of pulsed contractions.
Pulsed contractions were originally described in C. elegans during interphase in the polarizing one-cell zygote (P0; Fig. 1, A–C [individual pulses indicated by white arrowheads in Fig. 1 B]; and Video 1, top; Munro et al., 2004). In these cells, pulsed contractions are associated with transient deep invaginations of the cell surface (Fig. 1, A and B, magenta arrows); this makes it more difficult to quantify local changes in the density of cortical factors during individual pulses because these measurements could be confounded by movements of the cortex in and/or out of the plane of focus. Therefore, we focused on pulsed contractions that occur at the two-cell stage in the anterior blastomere (AB; Fig. 1, D–F [individual pulses indicated by white arrowheads in Fig. 1 E]). As in P0, pulsed contractions in the AB involve transient accumulations of F-actin and Myosin II; they are associated with transient local contractions of the actomyosin cortex (Fig. 1 F), but they are not associated with pronounced invaginations of the cell surface.
Single-molecule analysis of actomyosin dynamics during pulsed contractions
As a key step toward distinguishing different mechanisms for pulsed contractions, we used single-molecule imaging and single-particle–tracking analysis to quantify the relative contributions of local turnover and redistribution to changes in F-actin and Myosin II density during individual pulses. As described previously (Robin et al., 2014), we used RNAi against GFP to obtain embryos expressing single-molecule levels of Actin::GFP or of the non-muscle myosin heavy chain fused to GFP (NMY-2::GFP) over the endogenous proteins (Fig. 2 A and Videos 2 and 4). We combined near-total internal reflection fluorescence (TIRF) imaging with single-molecule detection and tracking to measure the appearance, motion, and disappearance of single-molecule speckles (Fig. 2, A–C; Robin et al., 2014). We assumed, as did others (Watanabe and Mitchison, 2002; Vallotton et al., 2004), that single-molecule appearance and disappearance events report directly on local rates of filament assembly and disassembly. We have shown previously that rates of turnover measured by single-molecule tracking agree well with those measured from single-molecule data by fitting kinetic models to photobleaching curves (Materials and methods; Robin et al., 2014).
We then devised new methods to measure simultaneously (a) single-molecule appearance rates, disappearance rates, and densities, and (b) local rates of cortical deformation on a moving and deforming patch of cortex during individual pulsed contractions (Fig. S1; see Materials and methods for details). In brief, we identified a reference frame for each pulse near the onset of contraction; within that frame, we identified a polygonal region of interest (ROI) containing the contracting patch (Fig. 2 A, dashed blue polygon, hereafter referred to as “the patch”). We propagated the patch forward and backward in time by extrapolating the displacements of tracked particles on or near its boundary (Fig. 2 D, Fig. S1, and Video 3). We then measured local deformation of the patch as frame-to-frame changes in patch area or by estimating a local strain rate from frame-to-frame displacements of the individual particles, with similar results in each case (Fig. 2 E and Fig. S2, A–C; Materials and methods). At the same time, we measured the number of molecules and single-molecule appearance and disappearance rates within the patch over time (Fig. 2, F–H). Finally, we aligned single-molecule measurements with respect to the onset or termination of individual contractions to produce a dynamic signature of actin assembly, disassembly, and deformation over the lifetime of a pulse (Fig. 3, A–D; and Fig. S2, D–F). These measurements allowed us to cleanly distinguish changes in single-molecule densities due to local assembly and disassembly from those due to local contraction (or expansion) of the cortical patch.
If pulses are initiated by positive feedback in which local contraction concentrates actomyosin and/or its upstream regulators, then the onset of actomyosin accumulation should coincide with the onset of contraction. Contradicting this expectation, we found that, on average, Actin:GFP began to accumulate ∼5 s before the onset of contraction (Fig. 3, B and F; and Fig. S2, A–C), during a period of time in which the cortex was locally expanding (Fig. 3 A). Approximately 30% of the total increase in Actin::GFP single-molecule density measured during a pulse occurred before the onset of contraction (Fig. S2, A–C). This initial accumulation was entirely due to a net imbalance of assembly and disassembly (Fig. 3 E): Before the onset of contraction, assembly rates increased (Fig. 3 C) and disassembly rates decreased (Fig. 3 D), leading to a sharp increase in the net rate of single-molecule accumulation that peaked at the onset of contraction (Fig. 3 E). During the contraction phase itself, the rate of change in single-molecule densities was determined almost entirely by a net imbalance of assembly/disassembly, with a very minor (6%) contribution from contraction itself (Fig. 3 E). Assembly rates decreased steadily, and disassembly rates increased steadily, such that a transition from net assembly to net disassembly (and from increasing density to decreasing density) occurred ∼7 s after the onset of contraction (Fig. 3 E). We obtained very similar results in embryos depleted of ARX-2, an essential subunit of the ARP2/3 complex (Fig. S3), suggesting that our results are not strongly biased by selective incorporation of Actin::GFP into branched versus unbranched F-actin (Chen et al., 2012).
Single-molecule analysis of GFP-tagged Myosin II (NMY-2::GFP) revealed local assembly/disassembly dynamics that were strikingly similar to those measured for GFP::Actin. (Fig. 3, G–K). On average, the density of single molecules of NMY-2::GFP began to increase ∼6 s before the onset of contraction during a period of local cortical expansion (Fig. 3 H), and ∼50% of this increase occurred before the onset of contraction. As observed for GFP::Actin, the upswing in Myosin II before the onset of a contraction was associated with both a sharp increase in appearance rates and a sharp decrease in disappearance rates (Fig. 3, I and J); the net rate of increase peaked at the onset of contraction, and during the contraction phase, the rate of change in single-molecule densities was determined largely by a net imbalance of assembly/disassembly, with a minor contribution from contraction (Fig. 3 K).
In summary, we found that changes in the actomyosin density during pulsed contractions are governed primarily by a dynamic local imbalance of F-actin and Myosin II appearance and disappearance rates. A large fraction of the increase in F-actin and Myosin II density during each pulse occurs before the onset of contraction, and local contraction accounts for only a minor fraction of the subsequent density increase during the contraction phase itself. We conclude that changes in actomyosin density during pulsed contractions are governed primarily by dynamic modulation of assembly and disassembly, not by local clustering of these factors or by dynamic coupling of contraction and advection.
Pulsed activation of RhoA drives the pulsed accumulation of F-actin and Myosin II
The observation that F-actin and Myosin II accumulate with very similar kinetics during pulsed contractions suggests that their accumulation is driven by a common upstream regulator. An obvious candidate is the small GTPase RhoA (encoded by rho-1 in C. elegans), which recruits and/or activates downstream effectors including formins, Rho Kinase, and Anillin to control F-actin assembly and Myosin II activation in a variety of cell types (Jaffe and Hall, 2005; Piekny and Glotzer, 2008). RhoA activity is required for pulsed contractions in P0 (Motegi and Sugimoto, 2006; Schonegg and Hyman, 2006; Tse et al., 2012), and a biosensor for active RhoA containing the RhoA-binding domain of Anillin (GFP::AHPH) accumulates locally during pulsed contractions in the zygote (Tse et al., 2012; Nishikawa et al., 2017).
To determine the relative timings of RhoA activation and actomyosin accumulation during pulsed contractions, we used a strain coexpressing GFP::AHPH (Tse et al., 2012) with an RFP-tagged version of the myosin heavy chain (NMY-2::RFP) to comonitor RhoA activity and Myosin II accumulation during individual pulses in AB. We observed a striking correlation between pulsed accumulation of GFP::AHPH and NMY-2::RFP during individual pulsed contractions (Fig. 4, A–C; and Video 5). GFP::AHPH accumulated rapidly within a broad domain that prefigured the initial accumulation of NMY-2::RFP, reached a peak near the onset of visible contraction, and then began to disappear before NMY-2::RFP (Fig. 4, B and C). NMY-2::RFP accumulated as discrete particles that increased in number and size, while contracting together into a smaller and tighter central domain. The initial accumulation of GFP::AHPH was diffuse, but over time, a fraction of the total GFP::AHPH became enriched within brighter punctae that colocalized with NMY-2::RFP (Fig. 4 B, yellow arrows). During the falling phase of each pulse, the diffuse pool of GFP::AHPH disappeared rapidly, while a smaller and more persistent fraction of GFP::AHPH remained colocalized with NMY-2::RFP particles (Fig. 4 B, yellow arrows).
To quantify these observations, we aligned data for multiple pulses from embryos coexpressing NMY-2::RFP and GFP::AHPH as follows (Fig. S3; Materials and methods). For each pulse, we smoothed and thresholded the NMY-2::RFP signal to identify a ROI containing high levels of NMY-2::RFP just before the onset of contraction (Fig. S3 A; Materials and methods). We propagated this ROI forward and backward in time (Fig. S3 A). We measured the mean intensities of the NMY-2::RFP and GFP::AHPH signals within the ROI before, during, and after the pulse and normalized these data with respect to the minimum (precontraction) and maximum intensities measured during this interval (Fig. S3 B). Then we aligned data for multiple pulses with respect to the time point at which NMY-2::RFP reached 25% of its maximum intensity (Fig. 4 D and Fig. S3 C). To distinguish different GFP::AHPH fractions, we thresholded the raw NMY-2::RFP signal to create a binary mask that distinguishes regions of high and low Myosin II intensity within the ROI (Fig. S3 D). Then we used this mask to decompose the normalized AHPH::GFP signal into “diffuse” and “myosin-colocalized” fractions (Fig. S3, E and F).
These aligned data confirm that sharp increases and decreases in total GFP::AHPH intensity precede, respectively, the appearance and disappearance of NMY-2::RFP (Fig. 4 D). On average, total GFP::AHPH reached 25% of its maximum intensity 8.6 ± 3.9 s before NMY-2::RFP (Fig. 4 F) and fell below 75% of its maximum intensity 11.1 ± 3.5 s before NMY-2::RFP (Fig. 4 G). Consistent with direct observations, the mean intensity of the diffuse fraction rose in phase with, but peaked earlier than, the total GFP::AHPH signal, while the mean intensity of the myosin-colocalized fraction rose in phase with NMY-2::RFP and peaked later than the total GFP::AHPH signal (Fig. 4, E–G).
We used the same approach to align data for embryos coexpressing NMY-2::RFP and the F-actin–binding domain of Utrophin (UTR) fused to GFP (GFP::UTR), a marker for F-actin (Burkel et al., 2007; Tse et al., 2012). Using NMY-2::RFP as the common reference to coalign data for NMY-2::RFP, GFP::AHPH, and GFP::UTR, we found that, like Myosin II, F-actin accumulates and dissipates during pulsed contractions with a significant delay relative to GFP::AHPH (Fig. 4, D, F, and G). Thus, local activation and inactivation of RhoA precedes the accumulation and disappearance of its downstream targets.
Finally, we used the time points at which NMY-2::RFP intensities (Fig. 4 D) and single-molecule densities of Myosin::GFP (Fig. 3 H) reached 25% of their peak values to align the time course of GFP::AHPH accumulation with respect to the onset of contraction as measured by single-molecule imaging. This analysis shows that the GFP::AHPH signal peaks just before (diffuse signal) or just after (total signal) the onset of contraction (Fig. 4 E, gray box). Thus, the local concentration of active RhoA by advection-contraction (Munjal et al., 2015) cannot explain the rising phase of RhoA activation in C. elegans embryos. Moreover, the diffuse GFP::AHPH signal peaks and begins to fall at a point where F-actin and Myosin II disappearance rates are at a minimum (Fig. 3, D and J); thus, factors other than cortical actomyosin disassembly drive the disappearance of active RhoA at the end of a pulse.
We used the same two-color image analysis to characterize coaccumulation of GFP::AHPH, NMY-2::RFP, and GFP::UTR during pulsed contractions in the P0. To remove the potentially confounding effects of large-scale cortical flows that occur during zygote polarization, we performed these measurements in embryos depleted of the centrosomal factor SPD-5 (Hamill et al., 2002; Munro et al., 2004), which exhibit pulsed contractions but lack cortical flows. In spd-5(RNAi) P0s, as in wild-type AB cells, increases in local GFP::AHPH intensity preceded the rise of NMY-2::RFP (Fig. S4, B–D) and to a lesser extent GFP::UTR (Fig. S4 D). As in AB, the initial accumulation of GFP::AHPH was broad and diffuse, while a more persistent pool of punctate GFP::AHPH colocalized with NMY-2::RFP within a more central region of the initial domain (Fig. S4 B, yellow arrows). The duration of GFP::AHPH and NMY-2::RFP accumulation was significantly longer in spd-5(RNAi) P0s than in wild-type AB cells. However, decomposition of the GFP::AHPH signal into diffuse and myosin colocalized fractions (Figs. S3 F and S4 E) revealed that the diffuse GFP::AHPH signal rose and fell with very similar temporal dynamics in AB and P0, while the myosin colocalized signal was far more persistent in P0 than in AB (Fig. S3, compare E and F).
Pulsed activation of RhoA does not require Myosin II
The observation that active RhoA peaks before the onset of contraction rules out models in which local contraction concentrates RhoA, or its upstream activators, to initiate pulses (Munjal et al., 2015). However, Myosin II could act in other ways to promote or shape pulsed activation of RhoA. To test this possibility, we used RNAi to deplete NMY-2 in a strain coexpressing transgenic GFP::AHPH and NMY-2 tagged with mKate2 at the endogenous locus by CRISPR-mediated homologous recombination (a gift of Dan Dickinson, University of Texas, Austin, TX). We found that local pulses of GFP::AHPH accumulation persisted in zygotes almost completely depleted of NMY-2, such that fewer than two NMY-2::mKate2-containing particles could be detected during a given pulse (Fig. 5, A and B; and Video 6). Focal accumulations of GFP::AHPH in nmy-2(RNAi) zygotes were comparable in size and spacing to those observed in spd-5 (RNAi) zygotes. However, the GFP::AHPH signal was entirely diffuse; the punctate signal that colocalizes with Myosin II in wild-type embryos was absent, and the mean time course for GFP::AHPH accumulation and disappearance was very similar to that measured for the diffuse pool of GFP::AHPH both in wild-type AB cells and in spd-5(RNAi) zygotes (Fig. 5 C). Assuming the diffuse AHPH signal reflects active RhoA, these data reveal the existence of a mechanism for locally pulsatile activation of RhoA, with similar kinetics in P0 and AB cells, that does not require the presence or activity of Myosin II. Whether the more persistent pool of myosin-colocalized AHPH reflects myosin-dependent recruitment and/or stabilization of active RhoA or binding of the Anillin AHPH domain to other cortical factors remains unclear (see Discussion).
Another key regulator of pulsed contractions is the scaffold anillin, which binds F-actin, Myosin II, active RhoA, and other cortical factors. The C. elegans anillin ANI-1 promotes focal accumulations of Myosin II during polarization and cytokinesis (Maddox et al., 2007; Tse et al., 2011), but whether it does so by promoting RhoA activation, or by recruiting/stabilizing Myosin II downstream of RhoA or both remains unclear. To distinguish these possibilities, we first used two-color image analysis to examine the relative timing of GFP::AHPH and GFP::ANI-1 accumulation during individual pulses using NMY-2::RFP as a common reference. In both AB and P0 cells, GFP::ANI-1 accumulated slightly before NMY-2::GFP, and thus with a pronounced delay relative to GFP::AHPH, although the delay was smaller in P0 cells (Fig. 6, A–H). In AB cells, the GFP::AHPH signal reached its maximum rate of accumulation before a detectable increase in GFP::ANI-1 (Fig. 6 C), suggesting that a local increase in ANI-1 is not required for rapid activation of RhoA during pulse initiation.
To test this further, we depleted zygotes of ANI-1 by RNAi in two different ways: RNAi against the full-length ani-1 sequence [ani-1FL(RNAi)] produced strong depletion of both ANI-1 and GFH::AHPH, while RNAi against a shorter N-terminal portion of the ani-1 sequence (ani-1N(RNAi)) produced weaker depletion but without affecting GFP::AHPH expression. In ani-1FL(RNAi) zygotes, we still observed pulsed accumulations of NMY-2::RFP, but the pulse duration was sharply reduced relative to wild-type or spd-5(RNAi) zygotes (Fig. S5, I, K, L, and N), consistent with previous reports (Maddox et al., 2007; Tse et al., 2011). In ani-1N(RNAi) zygotes, the duration of NMY-2::RFP pulses was also sharply reduced, although to a lesser extent than in ani-1FL(RNAi) (Fig. 6, J and M). The duration of GFP::AHPH pulses was also strongly reduced, mainly due to the reduction of the myosin-colocalized signal (Fig. 6, J and M). In contrast, the diffuse GFP::AHPH signal rose and fell with kinetics that were very similar to those observed in spd-5(RNAi) zygotes (Fig. 6, compare L and M). These data argue against an important role for ANI-1 in promoting activation of RhoA during pulse initiation and instead support a key role for ANI-1 in stabilizing Myosin II downstream of RhoA.
RhoA feeds back locally to promote its own activity, and this is required for pulse initiation
What drives the rapid rise in RhoA activity during pulse initiation? Plotting the rate of change in GFP::AHPH intensity versus intensity during the rising phase of individual pulses in either P0 or AB cells revealed a sharp increase in the rate of RhoA activation with increasing RhoA (Fig. 7 A). This is consistent with a scenario in which active RhoA feeds back positively to promote further activation of RhoA. However, it could also reflect pulsed activation of RhoA (without feedback) by an upstream activator. To distinguish these possibilities, we used RNAi to progressively deplete embryos of RhoA. If the time course of RhoA activation is dictated by an upstream activator, we should observe pulsed accumulation of GFP::AHPH as long as it remains detectable at the cortex. In contrast, if positive feedback of RhoA onto itself drives pulse initiation, then there should be an abrupt loss of pulsing below a threshold level of RhoA. Consistent with the latter expectation, we observed an abrupt transition from pulsed to nonpulsed RhoA accumulation after ∼12 h of feeding (Fig. 7, B and C; and Video 7). In zygotes that lacked pulsed RhoA accumulation, we could still readily detect robust RhoA-dependent cortical flows (Motegi and Sugimoto, 2006; Schonegg and Hyman, 2006) during polarity establishment (Fig. 7 B, dashed yellow lines) and localized accumulation of active RhoA before cytokinesis (Fig. 7 B, cyan arrowheads). Together, these observations support the idea that RhoA feeds back positively to amplify its own activation and that sufficiently strong feedback is required to generate local pulses of high RhoA activity.
Delayed accumulation of the Rho GAPs RGA-3/4 underlies pulse termination
What terminates RhoA activity during each pulse? Our results imply that local termination of RhoA activity at the end of a pulse does not require Myosin II activity or (by extension) its local inhibition by Myosin phosphatase (Piekny and Mains, 2002; Piekny et al., 2003; Diogon et al., 2007), nor is it timed by cortical disassembly. Previous studies identified the redundant RhoA GAPs RGA-3 and RGA-4 as inhibitors of RhoA activity during polarization and cytokinesis (Schmutz et al., 2007; Schonegg et al., 2007; Tse et al., 2012; Zanin et al., 2013). A YFP-tagged RGA-3 transgene accumulates at the cortex in early embryos (Schonegg et al., 2007), and simultaneous depletion of RGA-3 and RGA-4 leads to hyper activation of RhoA and hypercontractility during zygotic polarization (Schmutz et al., 2007; Schonegg et al., 2007; Tse et al., 2012). We wondered, therefore, if RGA-3/4 could provide delayed negative feedback to terminate RhoA activity during individual pulses.
To test this possibility, we first imaged embryos coexpressing GFP::RGA-3 (Zanin et al., 2013) and NMY-2::mKATE. Focusing on AB, and using two-color analysis as above, we confirmed that GFP::RGA-3 is present throughout the cortex but accumulates locally during individual pulsed contractions (Fig. 8, A–C; and Video 8). Significantly, GFP::RGA-3 and NMY-2::mKATE accumulated with very similar timing (Fig. 8 B). Using NMY-2::mKATE and NMY-2::RFP as common signals to co-align data for GFP::AHPH and GFP::RGA-3, we inferred that GFP::RGA-3 accumulates with an ∼6-s delay relative to GFP::AHPH, the rate of RhoA activation peaks before the onset of GFP::RGA-3 accumulation, and rapid accumulation of GFP::RGA-3 coincides with deceleration and then reversal of RhoA activation (Fig. 8 C, bottom). Together, these observations suggest that delayed accumulation of RGA-3/4 plays a key role in terminating each pulse of RhoA activity. To test this further, we created a strain in which GFP::AHPH and NMY-2::mKATE were coexpressed in rga-3; rga-4 (rga-3/4) double-mutant embryos (Zanin et al., 2013). Consistent with previous reports (Schmutz et al., 2007; Schonegg et al., 2007; Tse et al., 2012; Zanin et al., 2013), during polarity establishment in rga-3/4 double-mutant zygotes, we observed hyperaccumulation of GFP::AHPH and hypercontractility, characterized by convulsive contractions of the anterior cortex and rapid anterior-directed cortical flows (Fig. 8 D, second column; and Video 9). However, we could no longer detect local pulses of GFP::AHPH in these embryos. To exclude the possibility that rapid flows sequester factors required for pulsed contractility to the extreme anterior pole, we used partial depletion of the myosin regulatory light chain (MLC-4) to attenuate contractility and cortical flows in rga-3/4 double-mutant zygotes or in control zygotes that were doubly heterozygous for rga-3 and rga-4 (Fig. 8 D and Video 9). In rga-3/4 heterozygotes partially depleted of MLC-4, cortical flows were sharply reduced, but pulsed accumulation of GFP::AHPH could be readily detected (Fig. 8 D, third column; and Video 9). In contrast, in rga-3/4 double-mutant zygotes partially depleted of MLC-4, cortical flows during the polarity establishment phase were slower than observed in wild-type embryos; GFP::AHPH was uniformly enriched, but we did not observe local pulses of GFP::AHPH accumulation (Fig. 8 D, fourth column; and Video 9). Together, these data suggest that negative feedback through delayed accumulation of RGA-3/4 plays a key role in terminating local pulses of RhoA activity.
RGA-3/4 colocalizes with F-actin and requires F-actin for cortical association
Recent work suggests that F-actin accumulation inhibits RhoA activity to promote excitable dynamics in echinoderm and frog oocytes and embryos (Bement et al., 2015). We wondered if F-actin might play a similar role in C. elegans embryos by promoting recruitment of RGA-3/4. Consistent with this possibility, two-color live imaging of embryos coexpressing GFP::RGA-3 and mCherry::Lifeact (a marker for F-actin; Pohl et al., 2012) revealed extensive overlap of GFP::RGA-3 and mCherry::Lifeact signals throughout the cortex (Fig. 9, A and D) and during individual pulses (Fig. 9, B and E) in both AB and P0 cells. Scatterplots showing a pixelwise comparison of the two signals (Fig. 9, C and F) and measurements of Pearson’s r and Mander’s (MRGA) correlation coefficients over the entire cortex (P0: r = 0.9, MRGA = 0.98; AB: r = 0.8, MRGA = 0.96) or during individual pulses (P0: r = 0.84 ± 0.03, MRGA = 0.97 ± 0.02; AB: r = 0.65 ± 0.05, MRGA = 0.89 ± 0.03; mean ± SD; n = 10 pulses) confirmed strong colocalization.
To test whether F-actin is required for cortical localization of GFP::RGA-3, we treated permeabilized zygotes (Carvalho et al., 2011; Olson et al., 2012) coexpressing mCherry::Lifeact and GFP::RGA-3 with Latrunculin A (LatA) to rapidly depolymerize F-actin. Acute treatment with 10 µM LatA led to a rapid collapse of the cortical F-actin network and complete disappearance of F-actin from the cortex within ∼60 s. GFP::RGA-3 disappeared from the cortex with nearly identical timing and remained closely associated with the F-actin network as it collapsed (Fig. 9 G and Video 10; n = 6 embryos). These data suggest that cortical localization of RGA-3 depends on continuous physical association with F-actin. In contrast, upon acute depolymerization of F-actin, RFP::NMY-2, GFP-AHPH, GFP::Anillin, and GFP::Septin became rapidly enriched within dense punctate structures (Fig. 9 H; Fig. S5, A and B; and data not shown), similar to what has previously been described in other contexts (Hickson and O’Farrell, 2008), suggesting that none of these other downstream targets of RhoA are sufficient for cortical recruitment of RGA-3 (at least in the absence of an intact cortical actin network).
We could no longer detect pulsed accumulation of AHPH after acute treatment of permeabilized zygotes with LatA. However, this could be due to the sequestration of RhoA itself, or the biosensor, into dense clusters, as described above. We found that GFP::AHPH was no longer enriched in dense clusters after LatA treatment in embryos mutant for the Septin unc-59 (Fig. 9 I). Pulsed accumulation of GFP::AHPH and NMY-2::RFP persisted in unc-59-mutant zygotes (Fig. S5, D and E). However, acute treatment of permeabilized unc-59 zygotes with LatA lead to rapid accumulation of a diffuse pool of GFP::AHPH throughout the cortex and a loss of pulsing (Fig. 9 I), analogous to what we observed in rga-3/4 double-mutant zygotes.
Fast positive and delayed negative feedback involving RhoA and RGA-3/4 can account quantitatively for locally pulsatile RhoA dynamics
Our data suggest that in the early C. elegans embryo, locally excitable RhoA dynamics arise independently of myosin activity through a combination of fast positive feedback on RhoA activity and delayed negative feedback via local F-actin–dependent recruitment of RGA-3/4 (Fig. 10 A). To test this idea, we asked if a simple model based on these feedbacks with parameters constrained by experimental measurements predicts excitable dynamics. We formulated the model as a pair of ordinary differential equations, describing local rates of change in active RhoA and RGA-3/4, based on the following assumptions: (a) RhoA is activated at a basal rate and feeds back positively to promote further RhoA activation, (b) RhoA feeds forward through F-actin assembly to promote local, reversible association of RGA-3/4 with the cortex, and (c) RGA-3/4 acts as a GAP to promote local inactivation of RhoA (see Materials and methods for modeling details). Consistent with our experimental observations (Fig. 6 A and Fig. 10, B and C), we assumed that autoactivation of RhoA is a saturating function of RhoA activity, represented by a Hill function with Hill coefficient n = 1. We assumed that inactivation of RhoA by RGA-3/4 obeys Michaelis-Menten kinetics. To account for the observed delay between an increase in RhoA and the sharp onset of RGA-3/4 accumulation (Fig. 7 C), we assumed an ultrasensitive dependence of RGA-3/4 accumulation rate on RhoA, with the steepness of the response governed by an exponent m (see Materials and methods for mathematical details).
We set values for basal RhoA activation and RGA-3/4 recruitment rates to very low values, based on the slow rates of increase in GFP::AHPH and GFP::RGA-3 observed before the sharp upswing of each pulse (Fig. 7 C). Then we estimated values for the model’s remaining parameters by fitting the right-hand sides of the model equations to the relationships between dRhoA/dt and dRGA-3/4/dt and RhoA and RGA-3/4, inferred from coaligned, averaged, and normalized GFP::AHPH and GFP::RGA-3 intensities (Fig. 10, B–G; see Materials and methods for details).
Given these parameter estimates, and depending on the exact values chosen for the basal RhoA activation and RGA-3/4 association rates, we observed one of two closely related behaviors: excitable dynamics, characterized by a stable rest state (low RhoA and low RGA), from which the system can be induced to generate a transient pulse by a small decrease in RGA (Fig. 9 B) or a small increase in RhoA (not shown), and oscillatory dynamics, in which the stable rest state is lost and the system undergoes transient pulses at regular intervals (Fig. 9 C). The transition from excitable to oscillatory dynamics could be induced by either a small increase in basal RhoA activation or a small decrease in basal RGA activation. This is consistent with our observations in nmy-2(RNAi) embryos that some patches of cortex are quiescent while others exhibit repeated pulses of activity at regular intervals (Fig. 5 A, yellow arrows) and suggests that, in the absence of contractility, the C. elegans cortex may be poised close to a transition between excitable (with a low threshold for excitation) and oscillatory dynamics. We conclude that a simple combination of positive and negative feedback loops, coupling local RhoA activity and RGA-3/4 accumulation, is in principle sufficient to explain pulsatile RhoA dynamics in early C. elegans embryos, independent of actomyosin contractility.
Pulsed contractility is a widespread mode of actomyosin contractility, but its mechanistic basis has remained poorly understood (Levayer and Lecuit, 2012; Gorfinkiel, 2016). Current models for pulsed contractility invoke mechanochemical feedback in which contractile forces produced by Myosin II couple in different ways with actomyosin assembly/disassembly to drive excitable or oscillatory dynamics. Proposed feedback mechanisms include tension-dependent motor-binding kinetics (Effler et al., 2006; Ren et al., 2009; He et al., 2010; Luo et al., 2012), tension-dependent filament assembly/stabilization (Hayakawa et al., 2011; De La Cruz and Gardel, 2015) or disassembly (Machado et al., 2014), tension-dependent activation of Myosin II via, e.g., Ca2+ (Kapustina et al., 2008) or RhoA (Koride et al., 2014), or modes of feedback in which local contraction advects and concentrates actomyosin and/or its upstream activators (Bois et al., 2011; Kumar et al., 2014; Munjal et al., 2015; Nishikawa et al., 2017). Here, we have identified a mechanism for pulse generation in early C. elegans embryos that does not require force production or redistribution of cortical factors by Myosin II. Using single molecule imaging and particle tracking analysis, we have shown that the rapid initial accumulation of F-actin and Myosin II begins well before the onset of contraction, at a time when the cortex is locally expanding. Redistribution of actomyosin by local contraction makes a minor contribution to the overall accumulation of actomyosin during each pulse. Instead, our data show that pulsed accumulation and disappearance of F-actin and Myosin II are determined primarily by local modulation of their assembly/recruitment and disassembly. Using two-color imaging, we have shown that, during each pulse, a biosensor for active RhoA begins to accumulate well before RhoA’s downstream targets F-actin, Myosin II, and Anillin. Active RhoA nearly reaches its peak level before the onset of contraction (Fig. 4 D), and then it begins to disappear well before its downstream targets. Significantly, locally pulsed activation of RhoA continues to occur on patches of cortex that contain only a few (one to two) particles of Myosin II, too few to produce significant local contractile stress. Thus a Myosin-independent RhoA pulse generator underlies pulsed contractility in early C. elegans embryos.
One important caveat is that the ANI-1 AHPH domain may bind septin, lipids, or other cortical targets in addition to active RhoA (Piekny and Maddox, 2010). Indeed, we observe what appears to be two phases of AHPH accumulation—a diffuse phase that accumulates early and a second phase that accumulates later and colocalizes with NMY-2. The diffuse pool is likely to reflect RhoA-GTP that is available to bind RhoA affectors. Whether the myosin-colocalized pool bids active RhoA or other cortical factors remains unclear. Importantly, the diffuse pool has very similar kinetics in wild-type AB and in polarization-deficient or myosin-depleted P0 cells, supporting the conclusion that a common myosin-independent mechanism for pulsed activation of RhoA operates in both P0 and AB cells.
The pulses of RhoA activity described here and in other contexts (Bement et al., 2015; Munjal et al., 2015; Mason et al., 2016) are strikingly reminiscent of excitable behaviors found in other systems, such as action potentials in neuronal (Izhikevich, 2007) or cardiac cells (Luo and Rudy, 1991), transient pulses of intracellular calcium release (Goldbeter, 1996), or pulses of actin assembly observed in motile cells (Weiner et al., 2007). Theoretical studies highlight two key ingredients for excitable dynamics: fast positive feedback and delayed negative feedback (Strogatz, 1994). The sharp acceleration of active RhoA accumulation that we observed during the early rising phase of individual pulses is a dynamic signature of fast positive feedback in which RhoA promotes its own activity. Stronger evidence that RhoA participates in a positive feedback loop that is essential for pulsing comes from our observation that depletion of RhoA below a certain threshold leads to an abrupt loss of pulsed contractions while having minimal effects on other RhoA-dependent functions, such as cortical flow during polarization (Motegi and Sugimoto, 2006; Schonegg and Hyman, 2006) or cytokinesis (Loria et al., 2012).
The mechanism for this feedback remains unclear. In AB cells, the rapid acceleration of RhoA activation during pulse initiation occurs before any visible accumulation of Myosin II, Anillin, or F-actin, and in P0, depletion of NMY-2 or ANI-1 has little effect on the rising phase of RhoA activation. Thus, with the possible exception of F-actin in P0, it seems very unlikely that accumulation of these downstream targets makes a significant contribution to positive feedback. A more likely possibility is that RhoA feeds back through one or more of its upstream activators, such as ECT-2, CYK-4, and NOP-1 (Tse et al., 2012). In the zygote, pulsing is sensitive to depletion of ECT-2 and NOP-1 (Tse et al., 2012; Naganathan et al., 2014). During cytokinesis, active RhoA can act as a cofactor to promote transactivation of the RhoGEF ECT-2 by the RhoGAP CYK-4 (Zhang and Glotzer, 2015). While CYK-4 is not required for pulsed activation of RhoA during polarization, it is possible that RhoA could feed back through NOP-1, a protein of unknown activity that is required for RhoA activation during interphase in P0 and AB (Tse et al., 2012). Identifying the molecular mechanism(s) for this feedback is an important goal for future studies.
Regarding delayed negative feedback, direct comparison of our single-molecule and two-color imaging data (compare Figs. 3 and 4) argues strongly against a mechanism in which pulses are terminated by actomyosin disassembly. Both total and diffuse GFP::AHPH signals peak and begin to fall sharply at a time when F-actin and Myosin II disassembly rates are lowest (compare Figs. 3 and 4). Indeed, our data suggest that local accumulation of ANI-1 downstream of RhoA stabilizes Myosin II (and possibly also F-actin) to prolong the contraction phase of each pulse. Instead, our data suggest that the redundantly acting RhoGAPs RGA-3/4 play a key role in providing the delayed negative feedback that terminates RhoA pulses. RGA-3/4 act as GAPs toward RHO-1 in vitro (Schonegg et al., 2007), and loss of RGA-3/4 leads to hyperactivation of RhoA in vivo (Tse et al., 2012). We find that during each pulse, RGA-3 accumulates with a delay of ∼6 s relative to active RhoA. Significantly, the rate of active RhoA accumulation peaks and begins to fall just as GFP::RGA-3 begins to accumulate, suggesting that rapid accumulation of RGA-3/4 plays a key role in timing the end of each RhoA pulse. Consistent with this possibility, RhoA activity is uniformly high, and pulsatile RhoA activation is completely abolished, in rga-3/4 double-mutant zygotes, even when contractility is attenuated to prevent sequestration of RhoA activators by hyperactive cortical flow. Whether factors other than RGA-3/4 contribute to terminating RhoA pulses remains unclear.
Together, these data suggest a model for locally excitable RhoA dynamics in which RhoA feeds back positively to promote its own activation and feeds back negatively with a delay through RGA-3/4 to promote its own inactivation. Indeed, when we formulate this model mathematically and constrain model parameter values to match the local dependencies of RhoA and RGA-3/4 accumulation rates on levels of RhoA and RGA-3/4 inferred from two-color imaging data, the model predicts locally pulsatile RhoA dynamics, and small tunings of the model’s parameters mediate interconversion between excitable dynamics and spontaneous oscillations. This simple modeling exercise establishes an internally consistent hypothesis for pulsatile contractility that can be confirmed and extended by future experiments.
What governs the recruitment of RGA-3/4 during each pulse? We found that GFP::RGA-3 colocalizes broadly and extensively with cortical F-actin in both P0 and AB. RGA-3/4 accumulates with the same timing as F-actin during each pulse, and depolymerizating F-actin abolishes this accumulation. This suggests a specific mechanism for delayed recruitment of RGA-3/4 in which RhoA promotes increased local F-actin assembly (potentially through the formin Cyk-1 [Severson et al., 2002]), and F-actin in turn recruits RGA-3/4. Interestingly, a recent study (Bement et al., 2015) suggests that RhoA and cortical F-actin form an excitable circuit, with RhoA as activator and F-actin as inhibitor, that propagates cortical waves of RhoA activity and F-actin assembly in oocytes and embryonic cells of frogs and echinoderms. The mechanism by which F-actin feeds back to inactivate RhoA in these cells remains unknown, but our observations in C. elegans support the idea, proposed by Bement et al. (2015), that a RhoGAP homologous (or analogous) to RGA-3/4 may be recruited by F-actin to mediate negative feedback in frog and echinoderm cells. A similar mechanism has recently been reported in the context of secreting in Drosophila melanogaster larval salivary glands (Segal et al., 2018), and a similar circuit design may underlie the propagation of actin waves observed in many motile cells (Allard and Mogilner, 2013).
It is also interesting to compare our observations to those made recently in the Drosophila germband (Munjal et al., 2015). In germband cells, pulsed accumulation of a biosensor (GFP fused to the AHPH domain of Drosophila anillin) appears to coincide with the local accumulation of F-actin and Myosin II and with the onset of contraction, and it is abolished by inhibition of Rho Kinase, an upstream activator of Myosin II. In C. elegans, by contrast, pulsed accumulation of GFP::AHPH precedes actomyosin accumulation and the onset of contraction by many seconds and persists in the almost complete absence of Myosin II.
To some extent, these differences could reflect the imaging methods used to detect the RhoA biosensor. Using near-TIRF imaging in C. elegans embryos, we detected two pools of the biosensor: a diffuse pool that begins to accumulate well before Myosin II and a second more punctate pool whose distribution strongly overlaps with Myosin II (Fig. 4 B; Fig. S3, F and G; and Fig. S5). Based on studies in other cells (Weiner et al., 2007), the diffuse pool may be more difficult to detect with the use of confocal microscopy. Thus, it remains possible that a diffuse pool of active RhoA accumulates before actomyosin in the Drosophila germband but escapes detection by confocal microscopy.
An alternative idea is that pulsed contractility is governed by locally excitable RhoA dynamics in both systems, but that different forms of positive feedback may contribute differently to driving the rapid upswing of RhoA activity, and that different mechanisms may operate to trigger pulses (by driving RhoA activity above a threshold for excitation). For example, in the Drosophila germband, a mode of feedback in which local contraction advects and concentrates active RhoA (or upstream activators) may be required to initiate pulses, whereas in C. elegans, local fluctuations in RhoA or RGA-3/4 levels may be sufficient to do so in the absence of contractility. Importantly, in the Drosophila germband, as in C. elegans, advection/contraction coupling accounts for only a fraction of the total accumulation of active RhoA during each pulse; thus, other modes of positive feedback must also make a significant contribution.
More generally, we hypothesize that the different modes of RhoA excitability that have been described in frog, echinoderm, C. elegans and Drosophila embryos share a deeper underlying mechanistic origin. We suggest that a comparative analysis of mechanisms for pulsing in these and other systems will be a very fruitful way to uncover core conserved circuitry for pulse generation and to understand the ways in which this core circuitry is tuned or accessorized in different contexts to achieve different functional outcomes.
Materials and methods
C. elegans culture and strains
We cultured C. elegans strains at 22°C under standard conditions (Brenner, 1974). Table 1 lists the mutations and transgenes used in this study. Unless otherwise specified, strains were provided by the Caenorhabditis Genetics Center, which is funded by the National Institutes of Health National Center for Research Resources.
RNAi was performed by the feeding method as previously described (Timmons et al., 2001). Bacteria targeting nmy-2, spd-5, rho-1, perm-1, arx-2, ani-1, and mlc-4 were obtained from the Kamath feeding library (Kamath et al., 2003). The L4417 plasmid targeting the entire GFP sequence (generated by the Fire laboratory and available at Rauzi et al. , Levayer and Lecuit , and He et al. ) was transformed into HT115(DE3) bacteria. For the Myosin depletion experiments, L4 larvae coexpressing GFP::AHPH and NMY-2::mKate2 were transferred to nmy-2 RNAi feeding plates 24–30 h before imaging. Strong depletion of myosin was verified by strong loss of cortical NMY-2::mKate2. For the ANI-1 depletion experiments, L4 larvae coexpressing GFP::AHPH and NMY-2::RFP, were transferred to ani-1 RNAi feeding plates 30–36 h before imaging. For experiments involving spd-5 RNAi, L4 larvae were transferred to feeding plates for 24–30 h before imaging. For the RhoA depletion experiments, synchronized young adults were transferred to rho-1 RNAi plates 8–16 h before imaging. For experiments involving mlc-4 RNAi, synchronized young adults were transferred to feeding plates for 12–16 h before imaging. For the LatA experiments, late L4 larvae were transferred to perm-1 RNAi plates 16–24 h before imaging. For experiments involving arx-2 RNAi, L4 larvae were transferred to feeding plates for 30–36 h before imaging.
We mounted embryos as described previously (Robin et al., 2014) on glass slides under #1.5 coverslips in 3–5 µl of standard Egg Salts containing ∼100 uniformly sized polystyrene beads (18.7 ± 0.03 µm diameter; no. NT29N; Bangs Laboratories). The beads acted as spacers and allowed us to achieve uniform compression of the embryo surface across experiments (Robin et al., 2014).
We performed all imaging at 21–23°C on a Nikon ECLIPSE-Ti inverted microscope equipped with a Ti-ND6-PFS Perfect Focus Unit. A laser merge module (Spectral Applied Research) controlled fast, tunable delivery of 481-nm and 561-nm laser excitation from 50 mW solid-state lasers (Coherent Technology) to a motorized TIRF illuminator. We adjusted the laser illumination angle to achieve near-TIRF illumination (Tokunaga et al., 2008). We collected images using a Nikon CFI Apo 1.45 NA oil immersion TIRF objective combined with 1.5× intermediate magnification onto an Andor iXon3 897 EMCCD camera. All image acquisition was controlled by using Metamorph software. We used ImageJ to set minimum and maximum pixel values and perform gamma adjustments on the original 16-bit image data before converting images to 8-bit red, green, blue format or grayscale format for display. We performed these operations identically for all images that are compared directly.
We performed single-molecule imaging as described previously (Robin et al., 2014). For NMY-2::GFP, we used a combination of RNAi against GFP and mild photobleaching in wide-field-illumination mode to reduce surface densities of GFP-tagged transgenic proteins to single-molecule levels. For GFP::Actin, which is expressed at very low levels in the strain that we used, we used mild photobleaching alone. For GFP::Actin and NMY-2::GFP, we imaged single molecules using 10% laser power (∼0.16 µW/µm2), with 100-ms exposures in continuous streaming mode (GFP::Actin and NMY-2::GFP), yielding an approximate photobleaching rate of 0.05/s (Robin et al., 2014).
Analysis of F-actin and Myosin II turnover
In previous work, we compared two methods for estimating local F-actin disassembly rates from single-molecule data (Robin et al., 2014). The first method (smPReSS) estimates average disassembly rates in a local region by fitting kinetic models to the approximately exponential decay in particle densities measured during photobleaching, assuming steady-state conditions. The second method relies on single-molecule detection and tracking and infers appearance and disappearance events directly from single-molecule trajectories. We showed that, under steady-state conditions and when Myosin II is inhibited to remove the effects of local contraction and cortical flow, these two methods yield estimates of local disassembly that agree to within 20%. During pulsed contractions, the steady-state assumption is not valid, and the effects of cortical flow cannot be ignored. Therefore, in this work, we relied exclusively on the particle-tracking method to measure local appearance, disappearance, and motion of single molecules.
In preliminary analyses, we found that single molecules of GFP::Actin and NMY-2::GFP move sufficiently slowly during pulsed contractions that we could obtain marginally better results by pre-averaging 10 consecutive frames of raw data to produce sequences of images at 1-s intervals. We performed single-molecule detection and tracking on this pre-averaged data using a MATLAB implementation of the Crocker-Grier method (Crocker, 1996; Pelletier et al., 2009). We then inferred single-molecule appearance and disappearance events and frame-to-frame single-molecule displacements directly from the single-molecule trajectories.
Measuring local deformations from single-molecule data
The key goal of our single-molecule analysis was to distinguish the relative contributions of local assembly/disassembly and local deformation/flow to changes in local density during pulsed contractions. To do so, it was essential to follow dynamic changes in assembly/disassembly on a moving and contracting patch of cortex, i.e., in a material (Lagrangian) frame of reference. We used single-molecule displacements to identify and track regions of cortex undergoing pulsed contractions as follows.
Tracking a polygonal ROI during individual pulses
For each pulse, we identified a reference time point at/near the onset of contraction by visual inspection of the time-lapse sequence. At this time point, we identified manually an elliptical region containing the patch of cortex undergoing contraction (Fig. S1 A). We computed the smallest polygon (the “convex hull”) containing all the particles detected within the elliptical region (Fig. S1 B). Each vertex of the reference polygon was thus associated with a single molecule detected on the cell surface. To propagate the polygonal ROI forward and backward in time, we computed the frame-to-frame displacement of each of its vertices, either from the displacement of a vertex-associated molecule or (once the molecule disappears) from a weighted average of the frame-to-frame displacements of nearby particles (Fig. S1 C). We then measured local deformation and turnover within this polygonal ROI as follows.
Measuring local deformation and turnover during individual pulses
We compared three different measures of local compression (or dilation) within the polygonal ROI from frame to frame (Fig. S2 A): the change in normalized surface area, a particle-based strain rate, and a material strain rate. We computed the change in surface area sA as the time-derivative of the normalized area of the polygonal ROI:
where is the mean surface area taken over all frames in the pulse sequence. To compute a particle-based strain rate, for each particle in the polygonal ROI, we computed the average normalized change in distance between that particle and its near-neighbors:
where is the distance between particle i and a near-neighbor particle j at time t, and the sum is taken over all neighbor particles within a disk of radius 20 pixels centered on particle i. We then averaged over all particles in the ROI to obtain a particle-based strain rate for the entire ROI:
Finally, to compute a material strain rate sM, we used a linear least squares regression method to estimate the local gradient of particle velocities (Landau and Lifshitz, 1987). We then decomposed the resulting velocity gradient tensor into antisymmetric (rotation) and symmetric (strain rate) components. We then took one-half the trace of the symmetric strain-rate tensor as a measure of local compressive strain. In practice, we found that all three methods yielded very similar results regarding the magnitude and timing of contractions (Fig. S2 A). We report results based on the particle-based strain rate in Figs. 2 and 3.
To quantify turnover rates for F-actin or Myosin II within the polygonal ROI, for each time point t, we measured the area of the ROI (At), the number of particles Nt, their density Dt = Nt/At,and the number of appearance and disappearance events that occurred within the ROI between time t and t + Δt(ΔNt+, ΔNt–). We quantified the mean appearance rate and the mean disappearance rates within the ROI as and . We computed the change in actin density within a polygonal ROI at time t as We estimated the contribution to the change in density from deformation of the ROI to be , where sp is the particle-based strain rate measured as described above, and the contribution from turnover (i.e., net imbalance of appearance and disappearance) to be , such that ΔDt = ΔDdeformation + ΔDturnover.
Two-color imaging, pulse tracking, and analysis
We performed two-color imaging using the imaging system described above with near-TIRF illumination. We performed the initial steps of image processing, pulse identification and extraction, using the software package ImageJ (National Institutes of Health). For all subsequent steps, including pulse tracking, intensity measurements, data normalization, and alignment across multiple pulses, we used custom functions written in MATLAB (Mathworks).
We imaged embryos coexpressing GFP- and RFP-tagged transgenes by alternating 100-ms exposures with 488-nm and 561-nm excitation, thus giving five two-color frames per second. We used 25% maximum laser power ≈0.4 µW/µm2 for each channel. For subsequent analysis, we averaged more than five consecutive frames to obtain a single image for each channel at 1-s intervals. We limited our analysis to individual pulses that moved very little during the period leading up to the onset of contraction. We used ImageJ to extract a subregion containing each pulse of interest. With the exception of Myosin-depleted embryos, we used NMY-2::RFP as a reference signal to track the location of the pulse through time. For the analysis of Myosin-depleted embryos, we used GFP::AHPH as the reference signal.
We used the reference signal to identify and track a moving ROI associated with each pulse as follows (Fig. S3 A). First, we smoothed each frame of the image sequence using a gaussian filter (sigma = 2–3 µm). We then thresholded the smoothed image to identify ROIs associated with the pulse in consecutive frames. We used the same value of sigma and the same threshold for all frames, and we chose these values such that each ROI in the sequence was simply connected and such that the largest ROI in the sequence was approximately the same size as the region of strong signal accumulation near the peak of the pulse in the unprocessed data, as viewed by eye.
To measure signal intensities versus time during a pulse, we first used MATLAB to determine the centroid of each of the ROIs identified above to obtain a sequence of centroid positions Ct = (xt, yt). We extended this sequence backward in time using the first measured centroid position Cfirst and extended it forward in time using the last measured centroid position Clast. We then used a single reference ROI, centered on positions to measure a sequence of GFP and RFP intensities over time. We compared two different reference ROIs: (a) the largest ROI measured in the sequence, which corresponds roughly to peak accumulation of Myosin II and (b) a “bounding box,” defined to be the smallest square region containing this largest ROI. Because the results did not differ significantly, we used the bounding box ROI for all measurements reported here.
For each probe, we measured the mean signal intensity in the entire ROI over time. We then normalized these data for individual pulses using the equation
where Imean(t) is the mean intensity in the ROI at time t, Imin is the minimum mean intensity in the ROI before the onset of contraction, Imax is the maximum mean intensity in the ROI measured over the entire sequence . Finally, we aligned data across multiple pulses with respect to the time point at which NMY-2::RFP crossed 25% of its normalized maximum intensity (Fig. S3 C). The mean was calculated with a 95% confidence interval.
To partition the AHPH::GFP signal into diffuse and myosin colocalized fractions, we defined a myosin colocalized mask to be the collection of pixels within the ROI in which the raw NMY-2::RFP signal was greater than a threshold percentage (75–95%) of its peak intensity, and the diffuse mask to be its complement (Fig. S3 D, top). We used these masks to decompose the total GFP::AHPH signal into diffuse and myosin-colocalized pools (Fig. S3 D, bottom). We then partitioned the total normalized intensity within the ROI into diffuse and myosin-colocalized signals given by
where FD(t) and FC(t) are the fraction of pixels, and ID(t) and IC(t) are the mean intensities, in the diffuse and myosin-colocalized regions respectively (Fig. S3, E and F). The timing with which the diffuse and myosin-colocalized signals rose and fell was insensitive to the exact threshold. We used 95% for the results reported here.
We performed several additional controls to assess the sensitivity of our results to variation across strains, or the details of pulse identification, tracking and intensity measurements. First, we compared the kinetics of Myosin II accumulation during pulses in embryos coexpressing NMY-2::GFP and NMY-2::RFP and confirmed that the dynamics of accumulation were essentially identical after normalizing for differences in expression level and probe brightness (Fig. 4, E and F; and Fig. S3, G and H). Second, we confirmed that the dynamics of Myosin::RFP accumulation was essentially identical across the different two-color strains that we used (Fig. S3, G and H). Finally, we verified that our measurements of the rate and relative timing of accumulation of different signals were largely insensitive to differences in the size of the box/blobs used (Fig. S3 and data not shown).
To produce the kymographs shown in Figs. 5,6,7,8,9, and S4, we aligned images so that the anterior–posterior axis of the embryo coincided with the horizontal (x) image axis. We selected rectangular regions aligned with the x image axis, whose width (in x) coincided with the embryonic ROI and whose height (in y) was 10–20 pixels. From the original image stack, we extracted an xyt substack corresponding to this rectangular region; we used ImageJ’s reslice tool to reslice this stack with respect to the xt plane, and then we used a maximum-intensity projection to collapse the individual slices in y to obtain a kymograph in x versus t. For the kymographs in Fig. 4 C and S4 C, we selected the indicated regions aligned with the vertical (y) axis and resliced with respect to the yt plane.
We analyzed colocalization of mCherry::Lifeact and GFP::RGA-3 signals using the ImageJ plugin Coloc_2 (https://imagej.net/Coloc_2). Given a pair of image stacks representing the two signals for a sequence of frames within a ROI, we first scaled the data for each stack to lie between 0 and 255, then used Coloc_2 to compute Pearson’s correlation coefficient and an associated significance value and to estimate Mander’s correlation coefficients with automated threshold determination (Costes et al., 2004).
Statistical analysis methods
We computed 95% confidence intervals for pulse intensity data assuming a normal distribution. We tested differences between samples using a two-tailed t test, assuming normal distributions with unequal variances.
Modeling RhoA pulse generation
We built a simple ordinary differential equation model for RhoA pulse generation at a single point in space based on autocatalytic activation of RhoA and delayed negative feedback via RhoA-dependent recruitment of RGA-3/4. We started with the following assumptions:
(1) RhoA is activated at a constant basal rate.
(2) Active RhoA feeds back to promote further RhoA activation at a rate that can be described as a Hill function of RhoA density.
(3) RGA-3 and RGA-4 can be treated as a single species (RGA-3/4) that acts as a GAP to promote local inactivation of RhoA.
(4) Active RhoA promotes local F-actin assembly; RGA-3/4 binds F-actin from an abundant cytoplasmic pool, and dissociates from F-actin at a constant rate. Because RGA-3 and F-actin accumulate with very similar kinetics, we did not model F-actin directly. Instead, we assumed that RGA-3/4 binds the cortex at a constant basal rate plus a rate that depends on the local density of active RhoA and that RGA-3/4 dissociates from the cortex at a constant rate.
(5) To capture the observed delay between the sharp upswing in RhoA activity and the onset of F-actin and RGA-3/4 accumulation (Fig. 4 D and Fig. 7 B), we assumed ultrasensitive dependence of RGA-3/4 recruitment rate on RhoA activity.
With these assumptions, letting ρ represent the density of RhoA and r represent the density of RGA-3/4, we write a pair of ordinary differential equations for ρ and r:
Estimating model parameter values
We estimated the values for model parameters; we first constructed empirical relationships between intensities of GFP::AHPH and GFP::RGA-3 (hereafter RhoA and RGA) and their time derivatives during pulses of RhoA activity. Intensity data for RhoA and RGA were normalized, aligned, and averaged over many individual pulses, then coaligned using Myosin::RFP as a common reference as described above (Fig. 7 C). We smoothed these data using a five-frame moving average and then estimated the time derivatives of signal intensity as the difference in intensities between consecutive frames. We then fit model equations to these data as follows: First, we fixed reference values for basal RhoA activation and RGA recruitment , consistent with the slow rates of increase in RhoA and RGA observed before the sharp upswing of each pulse (Fig. 7 A). To estimate parameters that govern RhoA dynamics, we first considered the rising phase of RhoA activation (colored segments of the curves in Fig. 10 B), when RGA is ∼0, and used nonlinear least squares regression to fit
to data for dRhoA/dt versus RhoA (Fig. 10 C). We set n = 1, based on the observed form of dependence of dRhoA/dt on RhoA (Fig. 10 C), fixed the value of , and used model fits to estimate values for . Next, holding values for fixed, we fit the full equation for dp/dt to data for dRhoA/dt versus RhoA and RGA for the termination phase (colored segments of the curves in Fig. 10 D) to estimate values for kGAP and KGAP (Fig. 10 E). These fits yielded a value for KGAP that was consistently negative and close to zero, corresponding to a scenario in which RGA-3/4 operates near saturation on active RhoA. Therefore, in subsequent analyses, we set KGAP to a small constant positive value (KGAP = 0.001) and used the model fit to estimate a value for kGAP. To estimate parameters that govern RhoA dynamics, we fit the full equation for dr/dt to data for dRGA/dt versus RhoA and RGA over the entire pulse (Fig. 10 G).
For each set of parameters determined as above, we set the initial values for r and p to zero, simulating a scenario in which RhoA is minimally active and a small perturbation reduces RGA-3/4 to a minimally observed level. We then solved the equations numerically using MATLAB to determine if this initial perturbation would result in either a single pulse of RhoA activity followed by a return to a stable inactive state (excitability) or a train of pulses (oscillatory dynamics).
Online supplemental material
Fig. S1 shows the schematic overview of methods for tracking a moving and deforming patch of cortex from single molecule data. Fig. S2 shows the comparison of different methods to quantify local deformation (strain rate) and to align data across multiple pulses. Fig. S3 shows the measurement, alignment, and decomposition of fluorescence intensity data during pulsed contractions. Fig. S4 shows the two-color analysis of pulsed contractions in P0. Fig. S5 shows Septin localization and contributions to pulsed contractions in one-cell embryos. Video 1 shows the actomyosin pulse dynamics in wild-type cells expressing GFP::UTR and NMY-2::RFP. Video 2 shows a two-cell embryo expressing GFP::ACT-1 at single-molecule levels and imaged in near-TIRF mode. Video 3 shows dynamic tracking of a cortical patch from single-molecule data. Video 4 shows a two-cell embryo expressing NMY-2::GFP at single-molecule levels and imaged in near-TIRF mode. Video 5 shows the pulse dynamics in a two-cell embryo expressing GFP::AHPH and NMY-2::RFP. Video 6 shows the pulse dynamics in a P0 expressing GFP::AHPH and NMY-2::mKate and subjected to strong nmy-2 RNAi. Video 7 shows the pulse dynamics in zygotes expressing GFP::AHPH and subjected to rho-1 RNAi for different periods of time. Video 8 shows the dynamics of GFP::RGA-3 and NMY-2::RFP accumulation during pulsing in AB cells. Video 9 shows the absence of RhoA pulsing in rga-3; rga-4 double mutants. Video 10 shows the coordinated response of cortical F-actin and RGA-3 to acute treatment with LatA.
We thank Amy Maddox, Katrina Longhini, Dan Dickinson, Esther Zanin, and Karen Oegema for worm strains; Michael Glotzer, Margaret Gardel, Chip Ferguson, and members of the Munro laboratory for valuable discussions; and Katrina Longini for technical assistance.
This work was supported by National Institute of General Medical Sciences grants GM09844 (to E.M. Munro) and T32GM007197 (to J.B. Michaux). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
The authors declare no competing financial interests.
Author contributions: E.M. Munro conceived the project. J.B. Michaux and F.B. Robin shared equally in its intellectual development. F.B. Robin developed and applied the single molecule imaging and particle tracking analysis and performed all the work reported in Fig. 2 and 3. J.B. Michaux developed the multicolor imaging approach to analysis of pulsed contractions and performed the experiments reported in Fig. 4, 5, 6, 7, 8, and 9. W.M. McFadden and E.M. Munro performed the modeling shown in Fig. 10. J.B. Michaux, F.B. Robin, and E.M. Munro prepared figures and wrote the paper.
J.B. Michaux and F.B. Robin contributed equally to this paper.
F.B. Robin’s present address is Sorbonne Universités, University Pierre and Marie Curie Paris 06, Centre National de la Recherche Scientifique, Institut National de la Santé et de la Recherche Médicale, Biologie du Développement Paris Seine—Institut de Biologie Paris Seine, Paris, France.