# A biochemical network controlling basal myosin oscillation

Xiang Qin, Edouard Hannezo, Thomas Mangeat, Chang Liu, Pralay Majumder, Jiaying Liu, Valerie Choesmel-Cadamuro, Jocelyn Mcdonald, Yiyao Liu, Bin Yi, et al.

► To cite this version:

Xiang Qin, Edouard Hannezo, Thomas Mangeat, Chang Liu, Pralay Majumder, et al. A biochemical network controlling basal myosin oscillation. Nature Communications, 2018, 9 (1), pp.1210. 10.1038/s41467-018-03574-5 . hal-04783417

HAL Id: hal-04783417

https://hal.science/hal-04783417v1

Submitted on 14 Nov 2024

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

DOI: 10.1038/s41467-018-03574-5

OPEN

# A biochemical network controlling basal myosin oscillation

Xiang Qin<sup>1,2,3</sup>, Edouard Hannezo<sup>1,4,5</sup>, Thomas Mangeat<sup>1,2</sup>, Chang Liu<sup>1,2,6</sup>, Pralay Majumder<sup>7</sup>, Jiaying Liu<sup>1,2</sup>, Valerie Choesmel-Cadamuro<sup>1,2</sup>, Jocelyn A. McDonald<sup>8</sup>, Yiyao Liu<sup>3</sup>, Bin Yi<sup>6</sup> & Xiaobo Wang<sup>1,2</sup>

The actomyosin cytoskeleton, a key stress-producing unit in epithelial cells, oscillates spontaneously in a wide variety of systems. Although much of the signal cascade regulating myosin activity has been characterized, the origin of such oscillatory behavior is still unclear. Here, we show that basal myosin II oscillation in *Drosophila* ovarian epithelium is not controlled by actomyosin cortical tension, but instead relies on a biochemical oscillator involving ROCK and myosin phosphatase. Key to this oscillation is a diffusive ROCK flow, linking junctional Rho1 to medial actomyosin cortex, and dynamically maintained by a self-activation loop reliant on ROCK kinase activity. In response to the resulting myosin II recruitment, myosin phosphatase is locally enriched and shuts off ROCK and myosin II signals. Coupling *Drosophila* genetics, live imaging, modeling, and optogenetics, we uncover an intrinsic biochemical oscillator at the core of myosin II regulatory network, shedding light on the spatiotemporal dynamics of force generation.

<footnote n="1">¹Université de Toulouse, UPS, F-31062 Toulouse, France. ²CNRS, LBCMCP, F-31062 Toulouse, France. ³Department of Biophysics, School of Life Science and Technology, University of Electronic Science and Technology of China, Chengdu, 610054 Sichuan, P. R. China. ⁴The Wellcome Trust/Cancer Research UK Gurdon Institute, University of Cambridge, Cambridge CB2 1QN, UK. ⁵Institute of Science and Technology Austria, Am Campus 1, A-3400</footnote>

<footnote>Klosterneuburg, Austria. ⁶Department of Anesthesia, Southwest Hospital, Third Military Medical University, Chongqing 400038, P.R. China. ⁷Department of Life Sciences, Presidency University, Kolkata 700073, India. ⁸Division of Biology, Kansas State University, Manhattan, KS 66506, USA. These authors contributed equally: Xiang Qin, Edouard Hannezo. Correspondence and requests for materials should be addressed to Y.L. (email: liuyiyao@uestc.edu.cn) or to B.Y. (email: yibin1974@163.com) or to X.W. (email: xiaobo.wang@univ-tlse3.fr)</footnote>
Tissue morphogenesis, the acquisition of the three-dimensional (3D) shape of tissues and organs, is driven by biological forces that are typically generated by molecular motors, such as non-muscle myosin II (Myo-II) pulling filamentous actin (F-actin)¹,². The activity of the Myo-II molecule, composed of two heavy chains (MHC), two regulatory light chains (MRLC), and two essential light chains, is regulated by dynamic phosphorylation and dephosphorylation of MRLC³,⁴. Phosphorylation and dephosphorylation of MRLC, mainly by Rho-associated protein kinase (ROCK) and myosin light-chain phosphatase (MLCP), respectively, govern an active or inactive state of the Myo-II molecule, thus controlling both its interaction with actin filaments and the generation of contractile force⁴,⁵.

Actomyosin networks are very dynamic and often undergo cycles or pulses of assembly and disassembly¹. Pulses of actomyosin networks have been observed in several *Drosophila* tissue morphogenesis processes, including apical constrictions of invaginating mesoderm cells and delaminating neuroblasts, ectoderm cell intercalation during germ band extension, follicle cells during oogenesis, and amnioserosa cells during dorsal closure, but also during *Xenopus* gastrulation and mouse embryo compaction⁶⁻¹⁷. Several competing explanations for this phenomenon have been put forward, including mechanosensitive assembly¹⁰, balance between contraction forces and turnover¹¹, or phosphorylation-dephosphorylation cycles of MRLC¹⁸. Consistent with a control via phosphorylation-dephosphorylation cycles, spatio-temporal correlation between ROCK-MLCP and Myo-II has been observed in various morphogenetic processes¹⁸⁻²². However, how such cycles are generated and regulated is still unclear.

Basal actomyosin pulses, also essential for tissue morphogenesis, occur at the basal domain of ovarian epithelial follicle cells during *Drosophila* oogenesis⁹, a biological process to give rise to mature eggs. This basal actomyosin oscillation is thought to create an isotropic mechanical force, along the dorsal-ventral (DV) axis, in order to constrict the underlying tissue and thus eventually give the mature egg its characteristic anterior-posterior (AP) elongation⁹. Studies⁹,¹⁵,²³ have demonstrated that Rho1 and ROCK positively control basal Myo-II intensity, while Flapwing (Flw), one catalytic subunit of MLCP, negatively governs both basal Myo-II intensity and the initiation of basal Myo-II oscillation.

Here, through a combination of time-lapse live cell imaging, genetic, and optogenetic tools, and a physical based reaction diffusion modeling, we demonstrate that basal Myo-II oscillation is controlled by a biochemical oscillatory network involving spatio-temporal patterns of ROCK and MLCP, giving rise to activator/inhibitor dynamics.

## Results

### Basal shift of Myo-II regulatory network during elongation.
Basal Myo-II oscillation starts from the developing early stage 9 and it lasts until the end of stage 10B during *Drosophila* oogenesis⁹. During these periods, Myo-II displays a prominent shift from the apical to basal domain. Since ROCK and myosin-binding subunit (MBS, the regulatory subunit of MLCP²⁴) are well known activator and inhibitor of Myo-II phosphorylation and its mediated contraction, we assessed the spatio-temporal distribution of both proteins in follicle cells. Consistent with Myo-II patterns, both ROCK and MBS proteins display prominent and gradual shifts from apical to basal domain from early stage 9 to stage 10B (Fig. 1a–d). As basal Myo-II activity plays a critical role in controlling tissue elongation⁹, we tested whether tissue elongation would be affected upon changes of ROCK or MBS activity. Firstly, we confirmed the knockdown efficiency of

MBS RNAi in follicle cells (Supplementary Fig. 1a, b). Consistent with their predicted effect on Myo-II activity, the activation of ROCK (by overexpression of ROCK active form) and the inhibition of MBS (by MBS RNAi expression) both significantly enhance the elongation of egg chambers, while the inhibition of ROCK (by ROCK RNAi expression) and the activation of MBS (by overexpressing MBS active form²⁵) both result in rounder egg chambers (Fig. 1e, f). These results indicate that ROCK and MBS are involved in regulating Myo-II activity and governing tissue elongation. We therefore asked whether ROCK, MBS, and other Myo-II regulators were involved in the oscillatory properties of Myo-II.

### Oscillatory phase shift within Myo-II regulatory network.
We started by assessing the spatio-temporal dynamics of the Rho1-ROCK-MLCP signaling cascade. Firstly, we noticed that Rho1 intensity and activity are strongly distributed at basal junctional membrane, with almost undetectable signal at basal medial cortex, where oscillating Myo-II is present (Fig. 2a, c and ref. ²³). In addition, we did not detect a prominent pulsatile pattern of basal membrane Rho1 intensity and activity, during periods of Myo-II pulses (Fig. 2b, d, k and Supplementary Movie 1). We next assessed the dynamics of ROCK and MLCP signals. Basal ROCK levels oscillate with the same period as basal Myo-II levels, but preceding Myo-II by nearly 1 min (Fig. 2e, f, j, l and Supplementary Movie 2). Next, we assessed the dynamic correlation between Myo-II and MLCP. Here, we tracked the spatio-temporal dynamics of MBS and Flw, the respective regulatory and catalytic subunits of MLCP²⁶. To validate whether our MBS fluorescence reporters reflect with precision the endogenous MBS signals during basal Myo-II oscillation, we compared MBS fluorescence reporters and MBS antibody staining at different stages of Myo-II oscillations, and found that MBS reporters correlated strongly with MBS antibody stainings (Supplementary Fig. 1c-i). Importantly, although the total MBS protein levels per basal area remain constant throughout Myo-II oscillations, MBS displays periodic cycles of redistribution towards the medial pool where Myo-II is located, with the same period, but prominently trailing behind by almost 1 min (Supplementary Fig. 1j, k, 2a-c, g, i and Supplementary Movie 3). Similar to the spatio-temporal dynamics of MBS signals, endogenous Flw signals, detected by the protein-trap line *flw-YFP-159*²⁷, are spatially colocalized with, and temporally synchronized, with MBS signals (Supplementary Fig. 2d-f, h, i and Supplementary Movie 4). Thus, Flw also displays periodic cycles of redistribution toward the medial Myo-II enriched pool, with the similar period and 1 min-delay (Fig. 2g-l and Supplementary Movie 5). Taken together, our results demonstrate that ROCK, Myo-II, and local MLCP display similar, but temporally shifted oscillatory patterns at the basal medial cortex (Fig. 2j, l), while membrane Rho1 level and activity remain constant (Fig. 2k). This indicates that, similar to observations in the *Drosophila* embryo²⁰, Myo-II basal oscillations are not dictated by an upstream Rho1 oscillation, thus leading us to ask whether oscillations could be intrinsically encoded in the signaling network downstream of Rho1.

### Oscillation is not dependent on actomyosin contraction force.
A key question is whether such oscillations derive from the mechanical or chemical properties of actomyosin networks. Using Y27632 to inhibit the ROCK kinase activity and thus Myo-II contraction²⁰ cannot answer this question precisely, as it is expected to affect both mechanics, via changing the contractility of Myo-II, and the chemical network by blocking Rock kinase
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

activity. Instead, to test whether oscillations derive from the
mechanical properties of the cytoskeleton, as hypothesized in
previous studies¹⁰,¹¹, we mosaically expressed a dominant nega-
tive form of MHC²⁸ in follicle cells. In control cells, Myo-II

oscillation is followed by periodic change of basal membrane area
reduction (by around a minute, Fig. 3a, b, e), as expected from the
contractile properties of actomyosin to reduce basal membrane
area⁹. Considering its role in blocking the Myo-II loading to F-
actin<sup>28</sup>, overexpression of dominant negative MHC inhibits the actomyosin contractility, and accordingly, no prominent change of basal area reduction can be observed in this condition (Fig. 3c, d). Strikingly however, this inhibition of actomyosin contractility did not affect basal Myo-II oscillation; strong Myo-II oscillations are still present, with periods and intensities undistinguishable from those of control cells (Fig. 3c-g). These results confirmed that actomyosin contractility does not control basal Myo-II oscillation, and thus rather suggest that, in this system, membrane area changes are a passive consequence of Myo-II oscillation, rather than an active feedback, in contrast to a previous model<sup>10</sup>. This led us to investigate whether basal Myo-II oscillation could instead arise from a biochemical oscillator.

**Oscillations initiate from a ROCK flow beginning at junction.** First question that arises is the nature of the Myo-II activator initiating a cycle of oscillation. As Rho1, the better characterized upstream activator of Myo-II, is restricted to the membrane, we asked whether and how Rho1 could provide cues to the medial-restricted Myo-II pool. Interestingly, we observed that ROCK oscillation starts at the basal junctional membrane, where Rho1 resides, and then propagates inward (Fig. 4a, b and Supplementary Movie 6). To quantify this, we applied a “mean field spatial analysis” method named as “optical flow” analysis after denoising of raw images (see Methods section), and detected strong flow movements for both ROCK and Myo-II signals in basal domain of oscillating follicle cells (Fig. 4c). ROCK flow movements occur earlier than those of Myo-II around 1 min (Fig. 4d, e), which is similar to the intensity temporal delay between these two proteins at the medial basal pool (Fig. 2j, l). Interestingly, this is in stark contrast to intercalating ectoderm cells, characterized by the inverse patterns between flow movement and signal intensity<sup>20</sup>.

This suggests that the initiation of intensity increase and movement could be dependent on the ROCK interaction with membrane active Rho1, which unfolds and activates ROCK<sup>29</sup>. To test this hypothesis, we used the chemical inhibitor Rhosin to block Rho1 activity in order to determine the effect on both intensity change and flow movement of ROCK signals. Compared with the control, treatment of Rhosin strongly results in the loss of ROCK distribution at basal junction (Fig. 4f-h). Consistently, follicle cells with the Rho1 inhibition have much weaker intensity increase and flow movement of ROCK signals (Fig. 4g, i). Thus, the interaction of ROCK with membrane active Rho1 appears to initiate both signal increase and flow movement. However, as Rho1 activity and concentration do not oscillate, such an interaction alone would simply produce a uniform concentration of ROCK. From the theoretical perspective, activator/inhibitor type dynamics have the potential to create such self-organized oscillatory patterns, even in the absence of possible upstream regulators (other than Rho1) that could be oscillatory<sup>30</sup>.

**Positive and negative regulators of basal ROCK signals.** Therefore, we checked how activating or inhibiting the members of the Rho1 to MLCP signaling cascade affected the levels and

dynamics of wild type (WT) ROCK. As expected, active and dominant negative forms of Rho1 respectively increase and decrease ROCK concentration in the basal medial pool (Fig. 5a and Supplementary Fig. 3a). It is believed that ROCK acts as an inhibitor of MBS, thus relieving the inhibitory effect of MLCP on Myo-II<sup>31</sup>. Surprisingly, inhibiting (resp. constitutively activating) MBS increased (resp. decreased) basal medial ROCK concentration 2-fold, arguing for the converse negative feedback, either direct or indirect, of MBS on ROCK (Fig. 5a and Supplementary Fig. 3a). Moreover, overexpression of an active form of ROCK strongly enhanced the accumulation of WT ROCK in the basal medial pool, indicating that open and active ROCK might be able to amplify the accumulation and flow movement of WT ROCK molecules. To validate these findings, we repeated the experiments while measuring ROCK levels via antibody staining instead, and found the same results (Supplementary Fig. 4).

Next, we hypothesized that this ROCK amplification could be either via direct control of ROCK accumulation, i.e., via an oligomerization mechanism with auto or trans-phosphorylation between different ROCK molecules<sup>29</sup>, or via the indirect control of Myo-II phosphorylation, in which the formation of Myo-II mini-filaments would serve as a scaffold for additional ROCK accumulation. To test the latter hypothesis, we investigated the effect of Myo-II phosphorylation-mimic mutant (Myo-II EE) on ROCK levels, and found on the contrary an inhibitory effect compared to control and unphosphorylatable mutant (Myo-II WT form and Myo-II AA form, in Fig. 5a and Supplementary Figs. 3a-f, 4e, f). A recent study demonstrated that Myo-II EE had much weaker motor activity and thus less contractile property<sup>32</sup>. However, in the presence of WT genetic background, over-expression of Myo-II EE in follicle cells strongly enhanced the formation of stress fibers and also the elongation of egg chambers, while overexpression of Myo-II AA significantly reduced both (Supplementary Fig. 3g-j). The different effects of Myo-II EE vs. AA forms on ROCK levels, together with on stress fibers and actomyosin contractility, contradicted the hypothesis of the Myo-II mini-filament-mediated ROCK signal amplification. Interestingly, although the aforementioned overexpression of dominant negative MHC inhibited the actomyosin contractility (Fig. 3c, d), we found that this did not affect the basal medial ROCK levels (Supplementary Fig. 4e, f). This further indicates that actomyosin contractility is not necessary for ROCK signal self-amplification. These data thus hint to the fact that ROCK kinase activity could self-amplify ROCK signals via its oligomerization mechanism.

Altogether, our findings demonstrate that ROCK can both display self-amplification, and be repressed by MBS.

**Modeling of oscillatory characteristics and phase shift.** Given the uncovered regulatory interactions, we next asked whether these would be sufficient, from a theoretical perspective, to give rise to self-organized Myo-II oscillations. We therefore mathematically modeled the evolution in time of ROCK, Myo-II, and MBS/flw concentrations (resp. denoted *R*, *m*, and *M*). Before asking about spatial flow patterns, we first made the simplifying

**Fig. 1** ROCK, Myo-II, and MBS signals have an apical to basal distribution switch during the periods of basal Myo-II oscillation and control morphogenetic tissue elongation. **a, b** Apical and basal views of follicle cells from egg chambers at early stage 9, middle stage 9, late stage 9, stage 10A, and stage 10B, marked by the presence of MyoII-mCherry (red) with ROCK-GFP (green) (**a**), and MyoII-mCherry (red) with MBS-antibody staining (green) (**b**), respectively. Both scale bars are 30 µm. **c, d** Quantifications of relative apical and basal ROCK/Myo-II intensities (**c**) and MBS/Myo-II intensities (**d**) at different stages. **e** Quantification of the A-P to D-V length ratio in early stage 9, middle stage 9, late stage 9, and stage 10 egg chambers expressing the indicated transgenes. **f** Morphology of early stage 9, middle stage 9, late stage 9, and stage 10 egg chambers expressing the indicated transgenes, staining by Armadillo and DAPI, 4',6-diamidino-2-phenylindole. All scale bars are 50 µm. *n* is the number of samples (cells in **c** and **d** from four independent egg chambers) analyzed. Error bars indicate ±s.d. NS means not significant, *p* < 0.001 means significant difference by Student's *t*-test
<figure><img src="image_7.png" /><figcaption>**Fig. 2** Myo-II pulses with its regulators ROCK and Flw at the medial basal cortex but not with basal junction Rho1. **a, c, e, g** Time-lapse series of the representative oscillating follicle cells, labeled with Rho1-GFP and Myoll-mCherry (**a**), Rho FRET activity and Myoll-mCherry (**c**), ROCK-GFP and Myoll-mCherry (**e**), Flw-YFP and Myoll-mCherry (**g**), respectively. In **c**, Top panel is the processed Rho FRET signal; middle panel is YFP channel only; bottom panel is Myo-II signal. All scale bars are 5 µm. **b, d, f, h, i** Quantifications of the dynamic change of Rho1-GFP and Myoll-mCherry intensities (**b**), Rho FRET activity and Myoll-mCherry intensity (**d**), ROCK-GFP and Myoll-mCherry intensities (**f**), Flw-YFP and Myoll-mCherry intensities (**i**) in one oscillating cell, and Flw-YFP and Myoll-mCherry (**h**) in a local basal medial region of one oscillating cell (marked by dotted square in **g**), respectively. Intensity of each channel is normalized to its mean in **b, d, f, h, i**, except that Rho FRET activity is original value in **d, j, k** Average temporal cross-correlation of Myoll-mCherry with Myoll-GFP (cyan), ROCK-GFP (magenta) and local Flw-YFP (yellow) (**j**), Myoll-mCherry with Myoll-GFP (cyan), total Flw-YFP (magenta), Rho FRET activity (yellow) and Rho1-GFP (blue) (**k**). **l** Time offset calculated from the cross-correlation analysis. **n** is the number of samples analyzed. All error bars indicate ±s.d. *p* < 0.001 means significant difference by Student's t-test</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Fig. 3** Myo-II undergoes periodic accumulation and disassembly no matter whether Myo-II is able to be loaded on F-actin or not. **a, c** Time-lapse series of the representative wild type (**a**) and dominant negative MHC-expressing (**c**) follicle cells, labeled with mCD8-GFP and MyoII-mCherry. Both scale bars are 10 µm. The digitized cell contour is color-coded on the basis of the percentage decrease in surface area relative to the maximum area captured during imaging, as indicated in the heat map. **b, d** Quantifications of the dynamic change of relative Myo-II intensity and rate of basal area reduction from the same wild type (**b**) and dominant negative MHC-expressing (**d**) follicle cells. **e** Myo-II intensity-area reduction cross-correlation coefficient in the wild type vs. dominant negative MHC-expressing follicle cells. **f, g** Quantifications of relative Myo-II intensities (**f**) and oscillating cycle time period (**g**) in the indicated transgene-expressing GFP-positive cells compared with wild-type cells in the same tissue sample. **n** is the number of samples analyzed. All error bars indicate ±s.d. **NS** means no significant difference by Student's t-test</figcaption></figure>

assumption of uniform basal concentrations which in non-dimensional form reads

$$ \frac{dR}{dt} = \frac{R}{\tau_R} (1 - M) $$

$$ \frac{dm}{dt} = - \frac{m}{\tau_m} (1 - AR + BM) + \frac{J}{\tau_m} $$

$$ \frac{dM}{dt} = \frac{(m - M)}{\tau_M} $$

where the model has been reduced to three equations with five
<figure><img src="image_9.png" /><figcaption>**Fig. 4** Myo-II pulsatility follows a flow movement of ROCK initiated at basal junction. **a** Time-lapse series of one representative follicle cell labeled with ROCK-GFP and Rho1-mCherry, followed by the time-lapse series analysis of flow velocity of ROCK signals. **b** Quantification of the dynamic change of mean ROCK intensity, which is normalized to its minimal value. The rectangle region is the period when ROCK is in close proximity to Rho1 at and near the basal junction. **c** Dynamics of ROCK and Myo-II signals by the analysis of “Optical flow” method. Upper panels show the time-lapse series of ROCK-GFP and Myo1-mCherry signals in one representative oscillatory follicle cell basal domain. Lower panels show the time-lapse series analysis of flow velocity of ROCK and Myo-II signals. **d** Quantification of the dynamic change of mean Myo-II and ROCK flow velocities in one oscillating follicle cell. Velocity of each channel is normalized to its mean. **e** Average temporal cross-correlation of Myo-II velocity with ROCK velocity. **f, g** Time-lapse series of the representative oscillating follicle cell labeled with ROCK-GFP and Rho1-mCherry, in either control (f) or Rhosin treatment condition (g). All scale bars are 5 µm. **h** Quantification of the percentage of cells with ROCK membrane distribution in control vs. Rhosin treatment conditions. **i** Quantification of the percentage of cells with ROCK flow movement in control vs. Rhosin treatment conditions. **n** is the number of samples analyzed. All error bars indicate ±s.d. *p* < 0.001 means significant difference by Student's t-test</figcaption></figure>
free parameters: the three characteristic time scales of ROCK, Myo-II, and MBS/flw recruitment ($τ_R, τ_m, τ_M$), and two interactions coefficients A and B, respectively describing the effect of ROCK and MBS/flw on Myo-II. As detailed in Supplementary Note 2, we found that the aforementioned combination of ROCK self-activation and ROCK inhibition by MBS/flw was sufficient to generically generate biochemical oscillations. Moreover, the model qualitatively reproduced well the time-delay between ROCK, Myo-II, and MBS/Flw successive oscillations (Fig. 5b, d and Supplementary Note 2). To be able to make more quantitative predictions, we wished to contain the parameter set, and quantified, in particular, the recovery times of all three components. We thus performed Fluorescence Recovery after Photobleaching (FRAP) experiments, with confocal variable volumes (see Methods section and Supplementary Note 1) to calculate diffusion coefficients and characteristic recovery times from each protein of interest. A classical diffusion fitting law showed different recoveries times $τ_i$ for Myo-II (half-life time of recovery: 60 ± 30 s), ROCK (15 ± 7 s), and Flw (38 ± 15 s) (Fig. 5c and Supplementary Fig. 5). We thus used these as inputs for our model. Importantly, a theoretical prediction is that the delays between ROCK and Myo-II (resp. Myo-II and MBS/Flw) should be closely related to the characteristic recovery time of Myo-II (resp. MBS/Flw), see Supplementary Note 2 and Supplementary Fig. 6a-e, l for a phase diagram. Generically, a theoretical sensitivity analysis showed that the period of the oscillation $T_p$, as well as the delays between ROCK, Myo-II ($t_1$), and MBS/Flw ($t_2$) were robust with respect to the remaining two parameters, and fitted well the data (Fig. 5d and Supplementary Note 2 for details, $T_p \approx 7$ min, $t_1 \approx 50$ s, $t_2 \approx 40$ s), so that the model reproduces accurately the observed oscillations and time sequence between the three components.

**Quantitative prediction of perturbation on the oscillator.** Next, we further tested our biochemical oscillator model. In particular, we sought to predict quantitatively the mutant experiments via in silico considerations. We first tested the effect of modifying the activity of Rho1, and found that above a certain value, it desensitizes the system, destroying oscillations, reminiscent of the effect of the external potential in the FitzHugh-Nagumo model of spiking neurons³³ (Supplementary Fig. 6h and Supplementary Note 2). Importantly, experimental activation and inhibition of Rho1 strongly increase and decrease basal Myo-II intensity respectively, both with Myo-II oscillation becoming undetectable (Fig. 6a-h). We then investigated theoretically the effect of modifying the activity of ROCK and MBS/Flw, keeping all other parameters constant (Supplementary Note 2 for details), and found for both, a robust positive correlation between the resulting period and amplitude of Myo-II oscillations (Fig. 5e-g and Supplementary Fig. 6a-e). Importantly, we tested experimentally this unexpected prediction, and found that activation and inhibition of ROCK respectively increase and decrease both Myo-II oscillatory amplitude and cycle time period (Fig. 6i-p). Strikingly, activation and inhibition of MBS respectively decrease and increase both Myo-II oscillatory amplitude and cycle time (Fig. 6q-x). This confirmed experimentally the model prediction of a robust positive correlation between the period of the oscillation and its amplitude (Fig. 5g). One should note that in the model, the amplitude of oscillation is more sensitive to change in ROCK/MBS activity than the experimental counter-part, which could be due to saturation effects in the kinetic coefficients entering in our simplified model, and/or additional secondary feedbacks.

**Modeling of spatio-temporal oscillations.** In a second model step, we wondered whether we could reproduce the observed diffusive flow of ROCK signals in the model scheme. For this, we introduced the spatial counter-part to the above model, assuming the exact same reaction kinetics, but additionally allowing ROCK, Myo-II, and Flw to diffuse spatially on the basal surface, as observed from the aforementioned FRAP experiments (Supplementary Fig. 5; see Methods section). Interestingly, we found that Myo-II recovery displays minimal 2D diffusion in the basal plane, and is thus dictated by turnover from the cytoplasmic pool (Supplementary Fig. 5). For ROCK and Flw, the recovery is with turnover from the cytoplasmic pool, but also a significant 2D diffusion along the basal area (Supplementary Fig. 5), which provides a potential explanation for the centripetal flows of ROCK and Flw observed during oscillations. Importantly, as Rho1 is only present at the edge of cell, numerical simulations of the model predict that ROCK is first activated at the boundary, and then propagates as a kinematic wave towards the cell basal center at a constant velocity, mirroring the data, before being “caught up” by MBS/Flw at the edge (Fig. 5h-j for a 1D numerical integration and Supplementary Fig. 6i-k for a 2D numerical integration).

**Confirmation of the biochemical oscillator by optogenetics.** To further and directly probe the spatio-temporal dynamics of the biochemical oscillator, we introduced an optogenetic method named as GFP-LARIAT (Supplementary Fig. 7a), which we recently established in *Drosophila* in vivo to create a light-inducible GFP trap and thus disturb the localization and function of GFP-tagged proteins²³,³⁴. Here, we used GFP-LARIAT to induce the clustering trap of ROCK-GFP, Flw-YFP, or MBS-GFP in order to assess the effect on Myo-II. Firstly, we confirmed that the clustering trap of ROCK-GFP or Flw-YFP/MBS-GFP by light-induced GFP-LARIAT reduces or enhances the phosphorylation of MRLC, respectively, in follicle cells, compared with control cells (Fig. 7a-d and Supplementary Fig. 8a-c). Clustering trap of ROCK-GFP gradually and strongly reduces pulsatile Myo-II signals, while clustering trap of Flw-YFP or MBS-GFP gradually and significantly enhances Myo-II oscillations, compared with control conditions (Fig. 7e-k, Supplementary Figs. 7b-g, 8d-i, and Supplementary Movies 7-12). Before light illumination, both ROCK and Flw/MBS are highly co-localized with Myo-II. With the clustering progression, both ROCK and Flw/MBS tend to separate from Myo-II, indicating the functional loss of GFP-tagged proteins in control of downstream signals, compared with control (Supplementary Figs. 7h, i, 8j). Therefore, our optogenetic results fully confirmed in a more direct manner the main conclusion acquired by genetics and modeling.

## Discussion

Here, we have demonstrated that a biochemical oscillator involving ROCK and MLCP, two key regulators of Myo-II, controls basal Myo-II oscillation, which is important in normal *Drosophila* egg chamber elongation. ROCK functionally links its upstream control (active Rho1 at and near basal junctional membrane) with its downstream substrates (Myo-II at basal medial region). This link is dependent on a directional flow movement of ROCK from junctional membrane to medial cortical region. As ROCK shuttles medially, the signals of ROCK get self-amplified by its kinase activity-dependent accumulation, which leads to phosphorylation and mini-filament accumulation of Myo-II. This in turn governs the redistribution of MLCP signals toward the regions of Myo-II
<figure><img src="image_10.png" /><figcaption>**Fig. 5** In silico predictions of basal Myo-II oscillations and flow movements. **a** Quantification of relative ROCK-GFP intensity in the basal medial pool from the indicated transgene-expressing mCD8RFP-positive cells compared with wild-type cells in the same sample. The representative images have been shown in Supplementary Fig. 3a, from which ROCK-GFP signals at the basal medial region have been analyzed. **n** is the number of samples analyzed. All error bars indicate ±s.d. NS means no significant difference, while *p* < 0.001 means significant difference by Student's t-test. **b** Summary schematics of the inferred interaction diagram between ROCK, Myo-II, and MBS/Flw, showing both positive and negative interactions. **c** The half-life recovery time of Flw-YFP, ROCK-GFP, and MyoII-GFP after photobleaching. **d** In silico predictions on the dynamics of the ROCK/Myo-II/Flw system, based on the interactions of **b** and the measured FRAP recovery kinetics in **c**, matching the experimentally observed oscillations. **e, f** Increasing (resp. decreasing) the activity of ROCK theoretically increases (resp. decreases) both the period and the amplitude of the oscillation (the control parameters were A = 0.8 and B = 0.5). **g** Predicted correlation between oscillation amplitude and period upon manipulation of either ROCK (red) or MBS/Flw (blue) activity. **h-j** Schematics (h) of adding a spatial component to the model; we assume that ROCK is preferentially activated on the junctional membrane due to Rho1 localization, and corresponding simulated kymographs for ROCK (i) and Myo-II (j) showing that oscillations still occur, but start at the junctional membrane, resulting in a diffusive flow from membrane to center</figcaption></figure>
accumulation, in order to inactivate ROCK and Myo-II, and restart a cycle (Supplementary Fig. 9).

This biochemical oscillator controlling basal Myo-II con-
tractility shares some similarities with the recently reported

oscillatory signaling cascade that regulates apical Myo-II con-
traction in the intercalating ectoderm cells²⁰; both oscillatory
signaling cascades contain ROCK, Myo-II, and MBS, and there is
a significant temporal delay between Myo-II and MBS signals,

<figure><img src="image_11.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>i</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>j</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>m</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>n</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>q</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>r</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>u</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>v</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>g</figcaption></figure>

<figure><img src="image_24.png" /><figcaption>k</figcaption></figure>

<figure><img src="image_25.png" /><figcaption>o</figcaption></figure>

<figure><img src="image_26.png" /><figcaption>s</figcaption></figure>

<figure><img src="image_27.png" /><figcaption>W</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>h</figcaption></figure>

<figure><img src="image_30.png" /><figcaption>l</figcaption></figure>

<figure><img src="image_31.png" /><figcaption>p</figcaption></figure>

<figure><img src="image_32.png" /><figcaption>t</figcaption></figure>

<figure><img src="image_33.png" /><figcaption>x</figcaption></figure>

<figure><img src="image_34.png" /></figure>
which not only allows the formation of biochemical oscillation but also shut off all biochemical signals. However, these two oscillatory systems are quite different in several aspects. Firstly, Rho1 seems to be only transiently in contact with ROCK but not to Myo-II or MBS in basal Myo-II oscillation, while Rho1 moves together with ROCK and actomyosin networks in apical Myo-II oscillation. This difference might stem from the subcellular localization variation of the Rho1–Myo-II signal cascade. Secondly, Myo-II signals have prominent spatio-temporal delay from ROCK signals during basal Myo-II oscillation, while ROCK and Myo-II signals are synchronized in ectoderm cells. Thirdly, we propose that ROCK/Myo-II flows in our system are diffusive and biochemical in nature, as the inhibition of actomyosin contractility (via the expression of the dominant negative form of MHC) did not affect ROCK levels or dynamics, providing a complementary mechanism to the force-driven physical advertisement of the actomyosin contraction networks. A possible explanation for the large difference in time cycle between these two Myo-II oscillations would be that the basal condensed stress fibers might block the physical advection of biochemical molecules and thus slow down the speed of flow movement and oscillating cycle period. It was striking in this regard to observe from FRAP experiments that Myo-II does not significantly diffuse laterally, with most of its dynamics associated with exchange kinetics with the cytoplasmic pool.

In the oscillatory ectoderm cells, apical E-cadherin difference along DV and AP axis membrane plays a major role in controlling the specific direction of Myo-II flow movement from the junctional membrane to medial apical cortex⁸,³⁵. However, basal E-cadherin levels are quite equal along DV vs. AP membrane²³, indicating that basal Myo-II flow movement is not due to the E-cadherin tension variance. We also excluded the possibility that the contraction of medial Myo-II motors triggers the flow of ROCK signals. Our studies indicate that signaling spatial difference is critically important in the formation of flow movements and signal oscillations (Figs. 2 and 5h–j). How exactly ROCK signals get self-amplified by its kinase activity, and then disassembled by MBS during flow movements, would be an important next step to further dissect molecularly this oscillation. In particular, the amplification and disassembly of ROCK signals might be via the direct interaction of ROCK/ROCK or ROCK/MBS, or via another intermediate factor which is linked with and possibly controlled by both ROCK and MBS.

Interestingly, although biomechanical effects, including actomyosin contractility, are dispensable for basal Myo-II oscillation, they still have key roles for morphogenesis. It is noteworthy that the occurrence of this biochemical oscillator is concomitant to strong tissue shape change. Indeed, any misregulation of the factors in this oscillator strongly changes the underlying tissue shape. Thus, it remains unclear if the oscillator has a defined

function for morphogenesis, or if it is a necessary consequence of increased contractility, due to the structure of the ROCK/Myo-II/ MBS regulatory network. Interestingly, however, different from the defective phenotype induced by exacerbated apical Myo-II contractility, enhanced basal Myo-II contractility results in the stronger phenotype of egg chamber elongation⁹,²³. It indicates that basal pulsatile actomyosin contractility might function to create a gradual, but not acute, tissue elongation during morphogenesis so that a functional organ is precisely developed.

## Methods

**Drosophila stocks and genetics.** The following fly stocks were used: Sqh::RLCmyosinII-mCherry and Sqh::RLCmyosinII-GFP⁶ (from Eric E. Wieschaus), DN Zipper (from Daniel Kiehart), Ubi::ROCK-GFP³⁶ (from Yohanns Bellaiche), UAS-ROCKRNAi, UAS-MBSRNAi (from Vienna Drosophila RNAi center), UAS-Sqh WT, UAS-Sqh EE, UAS-Rho CA, UAS-ROCK CA, UAS-MBS N300, Rho1-GFP, Rho-mCherry, UAS-Rho DN, ROK[1]/FM7 (from Bloomington Drosophila stock center), hsGal4/CyO, MKRS/TM6B was used to express UAS lines in follicle cells in the experiments of Rho FRET dynamics. For ROCK/myo-II and ROCK/Rho1 dynamics, Ubi::ROCK-GFP on the third chromosome were combined with ROK[1]/FM7 mutant flies in order to prevent the side effect of ectopic ROCK overexpression. Clones were generated using FLP-OUT technique by crossing UAS transgenic flies with: (1) P[hspr70-flp]; Sqh::Sqh-mCherry; UAS-mCD8GFP, AyGal4; (2) P[hspr70-flp];+/-; UAS-mCD8GFP, AyGal4; (3) P[hspr70-flp];+/-; Ay(CD2)Gal4. To detect ROCK signal in mosaic clones, Ubi::ROCK-GFP on the X chromosome were needed to be combined with UAS transgenic flies first, and then clones were generated using FLP-OUT technique by doing the second cross with P [hspr70-flp];+/-; UAS-mCD8RFP, AyGal4. All stocks and crosses were maintained at room temperature. For signal analysis in mosaic clones, hsFLPase was induced for 1 h at 37 °C twice with a 5 h interval, then flies were kept at 18 °C for 2 days and then fattened at 25 °C for overnight before dissection. For the analysis of tissue elongation, P[hspr70-flp];+/-; Ay(CD2)Gal4 was used for cross and later heat shock treatment, which can induce more than 90% clone cells in egg chamber. hsFLPase of this mosaic system was induced for 1 h at 37 °C once, then the flies were kept at 18 °C for 1–2 days and then fattened at 25 °C overnight before dissection. For the experiments of Rho FRET dynamics, hsGal4 flies were incubated at 37 °C for 1 h and the flies were kept at 25 °C for 5–6 h before dissection. For LARIAT experiments, hsGal4 was used to induce LARIAT expression in ovary follicle cells, and all steps were carried on in dark condition, including cross maintenance, and heat shock. Flies were treated with heat shock at 37 °C for 1 h, and then incubated 2–3 h at 25 °C before dissection. *Drosophila* ovaries were dissected and egg chambers were mounted under red light before blue light illumination.

**DNA constructs and transgenic fly generation.** Full length MBS cDNA was obtained from Genebank and amplified by primers listed below. The plasmids of Ubi-MBS-GFP and Ubi-MBS-RFP were generated by in-fusion cloning (Clontech). The flies were generated by Bestgene Inc. using w1118 fly.

**PCR to amplify MBS.** Forward 5'-
GGTGGAGGTGGTGGTATGTCCTCGCTGGACGC-3'
Reverse 5'- AATTGGGGTACGTCTAGACCAACTGGTAATGGTAGCGAC-3'

**PCR to amplify GFP.** Forward 5'-
CCGGGCTGCAGGAATTCATGGTGAGCAAGGGCGAG-3'
Reverse 5'- ACCACCACCTCCACCCTTGTACAGCTCGTCCATGC-3'

**Fig. 6** Experimental perturbation confirmation of in silico prediction on basal Myo-II oscillations. **a, e, i, m, q, u** Basal views of follicle cell clones expressing the indicated transgenes, marked by coexpression of mCD8GFP. Myo-II signals are monitored by MyoII-mCherry. All scale bars are 10 µm. **b, f, j, n, r, v** Quantifications of relative Myo-II oscillation amplitude in the indicated transgene-expressing GFP-positive cells compared with wild-type cells in the same sample. **c, g, k, o, s, w** Quantifications of relative Myo-II oscillating cycle time period in the indicated transgene-expressing GFP-positive cells compared with wild-type cells in the same sample. For both Rho CA-expressing and DN Rho-expressing cells, basal Myo-II oscillations are almost undetectable. Thus, relative Myo-II oscillation amplitude is zero in these two conditions. Symbol ∞ means that oscillating cycle time period is not prominent, also in these two conditions. The relative Myo-II average intensity from these analyzed basal Myo-II oscillations are 1 ± 0.23 vs. 2.98 ± 0.79 in control vs. Rho CA cells, 1 ± 0.20 vs. 0.19 ± 0.09 in control vs. DN Rho cells, 1 ± 0.27 vs. 2.17 ± 0.45 in control vs. ROCK CA cells, 1 ± 0.16 vs. 0.27 ± 0.11 in control vs. ROCK RNAi cells, 1 ± 0.23 vs. 1.95 ± 0.52 in control vs. MBS RNAi cells, and 1 ± 0.21 vs. 0.36 ± 0.12 in control vs. MBS N300 cells, respectively. **d, h, l, p, t, x** Quantifications of the dynamic change of medial basal mean Myo-II intensity in the indicated transgene-expressing cells compared with one wild-type cell in the same sample. Each pink curve is the dynamic change from the average of Myo-II intensities in different cells with a range marked by pink color, and blue curve is the dynamic change of Myo-II intensity in one representative wild-type cell. **n** is the number of samples analyzed. All error bars indicate ±s.d. *p* < 0.001 means significant difference by Student's *t*-test
**PCR to amplify mCherry.** Forward 5'-

CCGGGCTGCAGGAATTCATGGTGAGCAAGGGCGAGG-3'

Reverse 5'- ACCACCACCTCCACCCTTGTACAGCTCGTCCATGC-3'

**Production of ROCK antisera.** To produce anti-ROCK antibody, a peptide corresponding to the C-terminus of *Drosophila* Rok-PA, amino acids 1350–1370 (CYNNNSTDGSKISPSQSTRSSY), was synthesized by YenZym Antibodies, LLC (South San Francisco, CA). The N-terminal C was added to couple the peptide to the KLH carrier for immunization. Two rabbits were immunized by YenZym,

<figure><img src="image_35.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_36.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_37.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_38.png" /><figcaption>e</figcaption></figure>

With LARIAT

<figure><img src="image_39.png" /><figcaption>g</figcaption></figure>

<figure><img src="image_40.png" /><figcaption>f</figcaption></figure>

With LARIAT

<figure><img src="image_41.png" /><figcaption>h</figcaption></figure>

<figure><img src="image_42.png" /><figcaption>i</figcaption></figure>

<figure><img src="image_43.png" /><figcaption>j</figcaption></figure>

<figure><img src="image_44.png" /><figcaption>k</figcaption></figure>

<figure><img src="image_45.png" /></figure>
YZ3971, and YZ3972. YZ3972 had the best response by ELISA and showed specific signal on fixed tissue at 1:2000 and recognized a predicted 160 kDa band.

**Dissection and mounting of Drosophila egg chamber.** One to three-day-old females were fattened on yeast with males for 1–2 days before dissection. *Drosophila* egg chambers were dissected and mounted in live imaging medium (Invitrogen Schneider's insect medium with 20% FBS and with a final pH adjusted to 6.9), by using a similar version of the protocol described in ref. ³⁷. Different from normal mounting condition, the egg chambers were slightly compressed to overcome the endogenous curvature. In this condition, basal oscillation pattern, intensity, and period were similar to those in the condition without compression.

**Imaging and FRET.** Time-lapse imaging was performed with a Zeiss LSM710 or Leica SP8 confocal microscope with a ×40, numerical aperture 1.3 inverted oil lens, a 488-nm argon laser, and a 568-nm green HeNe laser. The basal focal plane, which is about 1 µm beneath the basal surface, was selected during live imaging to maximize the basal Myo-II intensity. For the dynamics of ROCK-GFP, MBS-GFP, MBS-RFP, or Flw-YFP, the similar basal focal plane was selected to maximize the basal signal intensity, as basal MyoII dynamic imaging did. The same microscope setup was used when comparing intensity between different samples. To view the signals at different focal planes, images were taken at different Z-stack layers from the basal surface to the apical side.

FRET images of live cultured egg chambers were acquired with a Zeiss LSM710 microscope, by using a similar version of the protocol described in ref. ³⁸. A 458 nm laser was used to excite the sample. CFP and YFP emission signals were collected through channel I (470–510 nm) and channel II (525–600 nm), respectively. To capture single, high-resolution, and stationary images, a ×40/1.3 oil inverted objective was used. CFP and YFP images were acquired simultaneously for most of the experiments. Sequential acquisition of CFP and YFP channels in alternative orders were tested and gave the same result as simultaneous acquisition.

**FRAP and photomanipulation.** For FRAP experiment, ROCK-GFP, Flw-YFP, and MyoII-GFP at the medial basal cortical regions were photobleached at around 4.2 micron size (488 nm laser at 100% power, 40–50 iterations) using a Zeiss LSM710 confocal microscope with ×40, numerical aperture 1.3 inverted oil lens. Following photobleaching, confocal images were acquired at the plane of medial basal cortex every 2 or 15 s.

For photoexcitation experiment, live-cell imaging was performed using a Zeiss LSM710 confocal microscope with ×40, numerical aperture 1.3 inverted oil lens, with a 488-nm argon laser and a 568-nm green HeNe laser. LARIAT clustering system was effectively induced by the blue light wavelength (400–510 nm), and thus here 488-nm argon laser was used to do the photoexcitation. To avoid the strong photobleaching effect on both GFP and RFP signals during photoexcitation, the 488 nm laser was set at 6% power level to do both GFP signal scanning and pulsed photoexcitation of LARIAT optogenetic system in a time-lapse imaging acquisition taken every 30 s.

**Drug treatments.** Egg chambers were dissected in live imaging medium, and then incubated with Rho inhibitor Rhosin³⁹ (Merk) at 250 µM for 20 min before being mounted for imaging.

**Image processing and data analysis.** Images were processed with MATLAB and Image J. For all images, the background (intensity of area without sample) was subtracted.

Image J were used to calculate the intensity of an individual cell as the average value of all pixels within the cell area. In time-lapse experiments, images were processed by MATLAB to correct photobleaching automatically. For dual-color imaging, the intensities were calculated from manually outlined cell areas if membrane-fluorescent protein was not present to mark cell boundaries.

For Rho FRET image, CFP and YFP images were first processed by ImageJ software. A background region of interest was subtracted from the original image. The YFP images were registered to CFP images using the TurboReg plugin. Gaussian smooth filter was then applied to both channels. The YFP image was

thresholded and converted to a binary mask with background set to zero. The final ratio image was generated with the MATLAB program, during which only the unmasked pixel was calculated³⁸. To determine the FRET signal ratio between basal junction and medio-basal region, each follicle cell was separated into junctional region and medial region (based on the Myo-II distribution region in control cells), and then the FRET signals in either region were analyzed using MATLAB.

The distribution of oscillation periods was generated by measuring the intervals between each pair of adjacent peaks. We applied autocorrelation to calculate the period of time series with different time offsets. This method averages out irregularities in the sequence and gives a similar average period. We found that autocorrelation was more robust and provided better results in analyzing irregular signals with a small amplitude, in some genetic backgrounds with strong reduction of Myo-II intensity compared with control⁹. To produce the average curves which show the dynamic changes of medial basal mean Myo-II intensity in Fig. 6d, h, l, p, t, x, we first synchronized all oscillating cycles, from different follicle cells with the same genetic backgrounds in several egg chambers, before this calculation. Each pink curve is the dynamic change from the average of Myo-II intensities in different cells with a range marked by pink color.

To analyze the dynamics of ROCK-GFP, Flw-YFP, and MyoII-GFP, we used the methods proposed from EMBL *FrapCalc* software for quantitative analysis of frap data (the software is available at https://github.com/framasoft/framacalc). The detailed analysis method has been outlined at Supplementary Information.

About the analysis of the ROCK-GFP and MyoII-RFP flow movements, there are two class of methods to estimate motion without analysis of single molecule fluorescent trajectories. The first class of methods is based on direct 2D cross correlation (in space domain or frequency domain) end are used on Particle Intensity Velocimetry (PIV) or in spatio-temporal correlation analysis like STICCS⁴⁰. These methods are insensitive to noise from the image, because noise is an uncorrelated information frame by frame. But the methods need a good window variable size analysis and we had some difficulty to find a good estimation of sub pixel motion (due to interpolation method of 2D correlation peak in the presence of noise). We thus found that using the other family of classical motion estimation, based on optical flow analysis, provided a better read-out. We used the classical Horn-Schunck algorithm⁴⁰,⁴¹ to estimate the flow of myosin from open source image analysis platform (ICY software). The algorithm assumes smoothness in motion over the area of analysis and converges to solutions which show more smoothness. But it is more sensitive to noise than 2D cross correlation methods and the restauration of the image need to be done before. Thus, a 3D median filter with the kernel size close to PSF size has been performed before the estimation of the optical flow. This step reduced a lot the sensitivity of the Horn-Schunk methods to salt and pepper noise from confocal microscopy. Furthermore, an equalization of the histogram was performed to reduce fluctuation of intensity frame by frame. The largely known contrast limited adaptive histogram equalization (CLAHE) method has been used from the open source image analysis platform (ICY software)⁴².

Tissue elongation was measured by the AP to DV length ratio of S9 and S10 egg chambers.

For quantification of cluster formation during photoexcitation of LARIAT system in *Drosophila* follicle cells, all time series of GFP images were automatically adjusted by MATLAB to the initial image intensity in order to reduce the effect of photobleaching. Then background noise was subtracted from the adjusted images before the measurement of clusters. Clusters were defined as discrete puncta of fluorescence with criteria of fluorescence intensity (1200–4095 arbitrary units), size (≥0.2 µm²), and circularity (0.35–1.0 arbitrary units). The area of clusters per cells was measured with MATLAB.

**Immunohistochemistry.** *Drosophila* ovaries were dissected in Schneider's medium and fixed with 4% formaldehyde for 20 min for most of antibodies, except of p-MLC, in which samples were fixed with 8% formaldehyde. After fixation, the egg chambers were rinsed with PBST (PBS with 0.3% Triton X-100) three times. The egg chambers were incubated with various first antibodies normally overnight in cold room. Antibody staining was performed as described previously. Anti-ROCK antibody (rabbit, 1:1000 dilution) was produced in this study. Anti-MBS antibody (rabbit, 1:400 dilution) was from Change Tan²⁴. Anti-Armadillo antibody (mouse

**Fig. 7** Basal Myo-II oscillations are inhibited by the light-induced ROCK-GFP clustering while enhanced by Flw-YFP clustering. **a, b** Basal views of ROCK-GFP/p-MLC signal (**a**) and Flw-YFP/p-MLC signal (**b**) marked by p-MRLC antibody staining in the LARIAT-expressing follicle cells, under the dark condition or after 1 h photoexcitation by visible light. Both scale bars are 10 µm. **c, d** Quantifications of average GFP/YFP clustering area (**c**) and relative p-MRLC intensity (**d**) before and after photoexcitation in the LARIAT-expressing follicle cells. **e, f** Time-lapse series of the representative LARIAT-expressing follicle cells, labeled with ROCK-GFP and MyoII-mCherry (**e**) and Flw-YFP and MyoII-mCherry (**f**), under blue light excitation. Both scale bars are 5 µm. **g, h** Quantifications of the dynamic change of relative Myo-II intensity and relative area of ROCK clusters (**g**) and Flw clusters (**h**) in the LARIAT-expressing follicle cells, under photoexcitation. **i, j** Quantifications of relative Myo-II intensity (**i**) and average GFP/YFP clustering area (**j**) before and after photoexcitation in the LARIAT-expressing follicle cells. **k** Quantification of oscillating cycle time period of Myo-II signals in the LARIAT-expressing follicle cells, labeled with MyoII-mCherry together with either ROCK-GFP or Flw-YFP, before and after photoexcitation. **n** is the number of samples analyzed. Error bars indicate ±s.d. *p* < 0.001 means significant difference by Student's t-test
N27A1, 1:50 dilution) was from the Developmental Studies Hybridoma Bank. Anti-p-MLC (rabbit, 1:10 dilution) was from Cell Signaling (catalogue number 3671). Secondary antibodies conjugated with Alex-488 and Alexa-561 (Molecular Probes) were used in 1:400 dilutions. Samples were imaged on a Zeiss 710 or SP8 confocal microscope.

**Statistics.** All data are presented as mean ± standard deviation. Statistical analysis to compare results among groups was carried out by Student's *t*-test with two distribution tails. A value of *P* < 0.05 was considered to be statistically significant, while a value of *P* < 0.001 was considered to be highly statistically significant.

**Code availability.** The codes used for analyses of different images (including FRET and LARIAT) are available from the corresponding author on reasonable request.

**Data availability.** The data sets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Received: 31 August 2017 Accepted: 20 February 2018

Published online: 23 March 2018

## References

1. Heisenberg, C. P. & Bellaiche, Y. Forces in tissue morphogenesis and patterning. *Cell* **153**, 948–962 (2013).

2. Kasza, K. E. & Zallen, J. A. Dynamics and regulation of contractile actin-myosin networks in morphogenesis. *Curr. Opin. Cell Biol.* **23**, 30–38 (2011).

3. Warrick, H. M. & Spudich, J. A. Myosin structure and function in cell motility. *Annu. Rev. Cell Biol.* **3**, 379–421 (1987).

4. Vicente-Manzanares, M., Ma, X., Adelstein, R. S. & Horwitz, A. R. Non-muscle myosin II takes centre stage in cell adhesion and migration. *Nat. Rev. Mol. Cell Biol.* **10**, 778–790 (2009).

5. Dawes-Hoang, R. E. et al. Folded gastrulation, cell shape change and the control of myosin localization. *Development* **132**, 4165–4178 (2005).

6. Martin, A. C., Kaschube, M. & Wieschaus, E. F. Pulsed contractions of an actin-myosin network drive apical constriction. *Nature* **457**, 495–499 (2009).

7. Solon, J., Kaya-Copur, A., Colombelli, J. & Brunner, D. Pulsed forces timed by a ratchet-like mechanism drive directed tissue movement during dorsal closure. *Cell* **137**, 1331–1342 (2009).

8. Rauzi, M., Lenne, P. F. & Lecuit, T. Planar polarized actomyosin contractile flows control epithelial junction remodelling. *Nature* **468**, 1110–1114 (2010).

9. He, L., Wang, X., Tang, H. L. & Montell, D. J. Tissue elongation requires oscillating contractions of a basal actomyosin network. *Nat. Cell Biol.* **12**, 1133–1142 (2010).

10. Koride, S. et al. Mechanochemical regulation of oscillatory follicle cell dynamics in the developing *Drosophila* egg chamber. *Mol. Biol. Cell* **25**, 3709–3716 (2014).

11. Dierkes, K., Sumi, A., Solon, J. & Salbreux, G. Spontaneous oscillations of elastic contractile materials with turnover. *Phys. Rev. Lett.* **113**, 148102 (2014).

12. An, Y. et al. Apical constriction is driven by a pulsatile apical myosin network in delaminating *Drosophila* neuroblasts. *Development* **144**, 2153–2164 (2017).

13. Simoes, S., Oh, Y., Wang, M. F. Z., Fernandez-Gonzalez, R. & Tepass, U. Myosin II promotes the anisotropic loss of the apical domain during *Drosophilaneuroblast* ingression. *J. Cell Biol.* **216**, 1387–1404 (2017).

14. Maitre, J. L., Niwayama, R., Turlier, H., Nedelec, F. & Hiiragi, T. Pulsatile cell-autonomous contractility drives compaction in the mouse embryo. *Nat. Cell Biol.* **17**, 849–855 (2015).

15. Gorfinkiel, N. From actomyosin oscillations to tissue-level deformations. *Dev. Dyn.* **245**, 268–275 (2016).

16. Skoglund, P., Rolo, A., Chen, X., Gumbiner, B. M. & Keller, R. Convergence and extension at gastrulation require a myosin IIB-dependent cortical actin network. *Development* **135**, 2435–2444 (2008).

17. Kim, H. Y. & Davidson, L. A. Punctuated actin contractions during convergent extension and their permissive regulation by the non-canonical Wnt-signaling pathway. *J. Cell Sci.* **124**, 635–646 (2011).

18. Vasquez, C. G., Tworoger, M. & Martin, A. C. Dynamic myosin phosphorylation regulates contractile pulses and tissue integrity during epithelial morphogenesis. *J. Cell Biol.* **206**, 435–450 (2014).

19. Valencia-Exposito, A., Grosheva, I., Miguez, D. G., Gonzalez-Reyes, A. & Martin-Bermudo, M. D. Myosin light-chain phosphatase regulates basal actomyosin oscillations during morphogenesis. *Nat. Commun.* **7**, 10746 (2016).

20. Munjal, A., Philippe, J. M., Munro, E. & Lecuit, T. A self-organized biomechanical network drives shape changes during tissue morphogenesis. *Nature* **524**, 351–355 (2015).

21. Kasza, K. E., Farrell, D. L. & Zallen, J. A. Spatiotemporal control of epithelial remodeling by regulated myosin phosphorylation. *Proc. Natl Acad. Sci. USA* **111**, 11732–11737 (2014).

22. Duque, J. & Gorfinkiel, N. Integration of actomyosin contractility with cell-cell adhesion during dorsal closure. *Development* **143**, 4676–4686 (2016).

23. Qin, X. et al. Cell-matrix adhesion and cell-cell adhesion differentially control basal myosin oscillation and *Drosophila* egg chamber elongation. *Nat. Commun.* **8**, 14708 (2017).

24. Tan, C., Stronach, B. & Perrimon, N. Roles of myosin phosphatase during *Drosophila* development. *Development* **130**, 671–681 (2003).

25. Lee, A. & Treisman, J. E. Excessive Myosin activity in mbs mutants causes photoreceptor movement out of the *Drosophila* eye disc epithelium. *Mol. Biol. Cell* **15**, 3285–3295 (2004).

26. Hartshorne, D. J., Ito, M. & Erdodi, F. Myosin light chain phosphatase: subunit composition, interactions and regulation. *J. Muscle Res. Cell Motil.* **19**, 325–341 (1998).

27. Yamamoto, S., Bayat, V., Bellen, H. J. & Tan, C. Protein phosphatase 1ss limits ring canal constriction during *Drosophila* germline cyst formation. *PLoS ONE* **8**, e70502 (2013).

28. Franke, J. D., Montague, R. A. & Kiehart, D. P. Nonmuscle myosin II generates forces that transmit tension and drive contraction in multiple tissues during dorsal closure. *Curr. Biol.* **15**, 2208–2221 (2005).

29. Riento, K. & Ridley, A. J. Rocks: multifunctional kinases in cell behaviour. *Nat. Rev. Mol. Cell Biol.* **4**, 446–456 (2003).

30. Epstein I. R., Pojman J. A. *An introduction to nonlinear chemical dynamics: oscillations, waves, patterns, and chaos* (Oxford University Press, New York, NY, 1998).

31. Totsukawa, G. et al. Distinct roles of ROCK (Rho-kinase) and MLCK in spatial regulation of MLC phosphorylation for assembly of stress fibers and focal adhesions in 3T3 fibroblasts. *J. Cell Biol.* **150**, 797–806 (2000).

32. Vasquez, C. G., Heissler, S. M., Billington, N., Sellers, J. R. & Martin, A. C. *Drosophila* non-muscle myosin II motor activity determines the rate of tissue folding. *Elife* **5**, e20828 (2016).

33. Fitzhugh, R. Impulses and physiological states in theoretical models of nerve membrane. *Biophys. J.* **1**, 445–466 (1961).

34. Lee, S. et al. Reversible protein inactivation by optogenetic trapping in cells. *Nat. Methods* **11**, 633–636 (2014).

35. Levayer, R. & Lecuit, T. Oscillation and polarity of E-cadherin asymmetries control actomyosin flow patterns during morphogenesis. *Dev. Cell* **26**, 162–175 (2013).

36. Bardet, P. L. et al. PTEN controls junction lengthening and stability during cell rearrangement in epithelial tissue. *Dev. Cell* **25**, 534–546 (2013).

37. Prasad, M., Jang, A. C., Starz-Gaiano, M., Melani, M. & Montell, D. J. A protocol for culturing *Drosophila melanogaster* stage 9 egg chambers for live imaging. *Nat. Protoc.* **2**, 2467–2473 (2007).

38. Wang, X., He, L., Wu, Y. I., Hahn, K. M. & Montell, D. J. Light-mediated activation reveals a key role for Rac in collective guidance of cell movement in vivo. *Nat. Cell Biol.* **12**, 591–597 (2010).

39. Shang, X. et al. Small-molecule inhibitors targeting G-protein-coupled Rho guanine nucleotide exchange factors. *Proc. Natl Acad. Sci. USA* **110**, 3155–3160 (2013).

40. Hebert, B., Costantino, S. & Wiseman, P. W. Spatiotemporal image correlation spectroscopy (STICS) theory, verification, and application to protein velocity mapping in living CHO cells. *Biophys. J.* **88**, 3601–3614 (2005).

41. Lecomte T., Thibeaux R., Guillen N., Dufour A., Olivo-Marin J. C. Fluid optical flow for forces and pressure field estimation in cellular biology. In: *2012 19th IEEE International Conference on Image Processing*, 69–72 (2012).

42. Pujadas A. B., Manich M., Guillén N., Olivo-Marín J. C., Dufour A. C. Biophysical measurements in 2D and 3D live cell imaging using fluid dynamics and optical flow. In: *2016 IEEE 13th International Symposium on Biomedical Imaging (ISBI)* 980–983 (2016).

## Acknowledgements

We thank Adam Martin, Yohanns Bellaiche, Daniel Kiehart, Bloomington *Drosophila* stock center and Vienna *Drosophila* RNAi center for flies. We thank Change Tan for MBS antibody. We thank Tri Toulouse platform for imaging acquisition advice. We thank Francois Schweiguth and Karine Belguise for discussion of manuscript preparation. This work was supported by the Institut National de la Santé et de la Recherche Médicale [the ATIP-Avenir program (2012–2016)]; Région Midi-Pyrénées Excellence program (2013–2016); the National Natural Science Foundation of China (grant number 81670552 to B.Y., 11772088, 31470906 to Y.L.); the National Science Foundation of USA (grant number NSF 1456053 to J.A.M.); and the Indian Foundation (DST [SB/YS/LS-88/2014] and UGC [MRP-MAJOR-ZOOL-2013-16812] to P.M.
**Author contributions**

X.Q. and C.L. performed image acquisition and transgene analysis. X.Q. and J.L. pro-
cessed and analyzed images. E.H. did the mathematical modeling. T.M. conducted the
data analysis of FRAP experiments. V.C.-C. made the constructs for transgenic flies.
P.M. and J.A.M. made the antibody for ROCK. X.Q., E.H., T.M., Y. L., B.Y., and
X.W. prepared the manuscript. X.Q., E.H. and X.W. designed the experiments. All
authors participated in the interpretation of the data and the production of the final
manuscript.

**Additional information**

**Supplementary Information** accompanies this paper at https://doi.org/10.1038/s41467-018-03574-5.

**Competing interests:** The authors declare no competing interests.

**Reprints and permission** information is available online at http://npg.nature.com/
reprintsandpermissions/

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_46.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018