Research Advance

Core PCP mutations affect short-time mechanical properties but not tissue morphogenesis in the Drosophila pupal wing

Max Planck Institute of Molecular Cell Biology and Genetics, Germany
DFG Excellence Cluster Physics of Life, Technische Universität Dresden, Germany
National Phenotypic Screening Centre, University of Dundee, United Kingdom
Center for Systems Biology Dresden, Germany
Laboratoire Physico-Chimie Curie, CNRS UMR 168, Institut Curie, Université PSL, Sorbonne Université, France
Max Planck Institute for Physics of Complex Systems, Germany
Université Paris Cité, Laboratoire Matière et Systèmes Complexes, France
Mildred Scheel Nachwuchszentrum P2, Medical Faculty, Technische Universität Dresden, Germany
Biotechnologisches Zentrum, Technische Universität Dresden, Germany

Dec 20, 2023

https://doi.org/10.7554/eLife.85581

Open access
Copyright information

Version of Record

Accepted for publication after peer review and revision.

Download
Cite
Share
CommentOpen annotations (there are currently 0 annotations on this page).

Version of Record updated: February 5, 2024 (This version)
Version of Record published: February 2, 2024 (Go to version)
Accepted Manuscript published: December 20, 2023 (Go to version)
Accepted: December 18, 2023
Received: December 15, 2022
Preprint posted: December 10, 2022 (Go to version)

1. Builds upon
Interplay of cell dynamics and epithelial tension during morphogenesis of the Drosophila pupal wing

Raphaël Etournay, Marko Popović ... Suzanne Eaton

Research Article Sep 22, 2015
Further reading

Abstract
Editor's evaluation
Introduction
Results
Discussion
Materials and methods
Appendix 1
Data availability
References
Article and author information
Metrics

Abstract

How morphogenetic movements are robustly coordinated in space and time is a fundamental open question in biology. We study this question using the wing of Drosophila melanogaster, an epithelial tissue that undergoes large-scale tissue flows during pupal stages. Previously, we showed that pupal wing morphogenesis involves both cellular behaviors that allow relaxation of mechanical tissue stress, as well as cellular behaviors that appear to be actively patterned (Etournay et al., 2015). Here, we show that these active cellular behaviors are not guided by the core planar cell polarity (PCP) pathway, a conserved signaling system that guides tissue development in many other contexts. We find no significant phenotype on the cellular dynamics underlying pupal morphogenesis in mutants of core PCP. Furthermore, using laser ablation experiments, coupled with a rheological model to describe the dynamics of the response to laser ablation, we conclude that while core PCP mutations affect the fast timescale response to laser ablation they do not significantly affect overall tissue mechanics. In conclusion, our work shows that cellular dynamics and tissue shape changes during Drosophila pupal wing morphogenesis do not require core PCP as an orientational guiding cue.

Editor's evaluation

This valuable study shows that the core PCP pathway, which commonly orients morphogenetic processes in development, does not establish global cues for cellular movements in the Drosophila pupal wing. Using a combination of laser-ablation experiments and mathematical modelling, it provides compelling evidence that the core PCP pathway only affects the fast timescale cellular response, but does not appear to drive overall tissue dynamics. While the signalling pathway guiding this morphogenetic process remains to be elucidated, these relevant findings challenge the role of the core PCP pathway in morphogenesis.

https://doi.org/10.7554/eLife.85581.sa0

Introduction

The spatial–temporal pattern of mechanical deformation during tissue morphogenesis is often guided by patterns of chemical signaling. Precisely how chemical signaling couples with the mechanics of morphogenesis, however, remains an active area of research. One conserved chemical signaling pathway that is known to be patterned across tissues is the core planar cell polarity (PCP) pathway, composed of a dynamic set of interacting membrane proteins that polarizes intracellularly within the plane of a tissue. Tissue-scale alignment of this pathway is known to orient cellular structures, such as hairs and cilia, and influence dynamic cellular behaviors during morphogenesis, such as cellular movements and cell divisions, through interactions with the cytoskeleton (reviewed in Devenport, 2014; Butler and Wallingford, 2017; Deans, 2021).

Here, we examine a potential role for the core PCP pathway in the dynamics and mechanics of morphogenesis using the Drosophila pupal wing. The Drosophila wing is a flat epithelium that can be imaged at high spatial–temporal resolution in vivo during large-scale tissue flows that elongate the wing blade (Aigouy et al., 2010; Etournay et al., 2015; Guirao et al., 2015). During the pupal stage, the proximal hinge region of the wing contracts and pulls on the blade region, generating mechanical stress that is counteracted by marginal connections mediated by the extracellular matrix protein Dumpy (Etournay et al., 2015; Ray et al., 2015). As a consequence, the tissue elongates along the proximal–distal (PD) axis and narrows along the anterior–posterior (AP) axis to resemble the adult wing. Both cell elongation changes and cell rearrangements are important for tissue deformation. To some extent, mechanical stress induces these cell behaviors. However, the reduction of mechanical stress in a dumpy mutant does not completely eliminate cell rearrangements, suggesting that there could be other patterning cues that drive oriented cell rearrangements (Etournay et al., 2015). We therefore wondered whether chemical PCP systems could orient cell behaviors, such as cell rearrangements, during pupal blade elongation flows.

In the Drosophila wing, there are two PCP systems termed Fat and core PCP (Matis and Axelrod, 2013; Adler, 2012; Devenport, 2014; Butler and Wallingford, 2017). The Fat PCP system consists of two cadherins Fat and Dachsous, a cytoplasmic kinase Four-jointed, and an atypical myosin Dachs. The core PCP system is composed of two transmembrane proteins Frizzled (Fz) and Flamingo or Starry night (Fmi, Stan), the transmembrane protein Strabismus or Van Gogh (Stbm, Vang), and the cytosolic components Dishevelled (Dsh), Prickle (Pk), and Diego (Dgo).

Our group has shown that tissue-scale patterns of PCP emerge during larval stages and then are dynamically reoriented during pupal tissue flows (Sagner et al., 2012; Aigouy et al., 2010; Merkel et al., 2014). At the onset of blade elongation flows, both systems are margin oriented, however as morphogenesis proceeds, core PCP reorients to point along the PD axis, whereas Fat PCP remains margin oriented until very late, when it reorients toward veins (Figure 1—figure supplement 1A, B; Merkel et al., 2014). Whether these PCP systems and their reorientation influence tissue dynamics and mechanics during blade elongation flows is unknown.

The core PCP pathway has been shown to influence numerous processes in Drosophila tissue development. These include hexagonal cell packing in the late pupal wing (Classen et al., 2005; Sugimura et al., 2016), as well as patterning of ommatidial clusters in the developing eye (Zheng et al., 1995; Jenny, 2010), orientation of cell division in sensory organ precursors (Gho and Schweisguth, 1998), formation of joints in the legs (Capilla et al., 2012), and regulation of tracheal tube length (Chung et al., 2009). In many cases, the mechanism connecting the core PCP pathway to cell dynamics and tissue mechanics is unclear. Recent studies suggest, however, that core PCP may act in concert with Nemo kinase to regulate cell rearrangements in the eye (Mirkovic et al., 2011; Founounou et al., 2021) and with the Drosophila NuMA ortholog Mud to orient cell division orientation in the sensory organ precursors (Ségalen et al., 2010).

Here, we examine cellular dynamics in tissues mutant for core PCP and we find that they are largely unperturbed, indicating that core PCP does not have an essential role in organizing global patterns of cell rearrangements in the pupal wing. We also performed an extensive analysis of the mechanics using laser ablation, developing a rheological model to interpret the results. We find that mutants in core PCP differ from wild type in the initial recoil velocity upon laser ablation. We find, however, that this difference is produced from the very fast timescale response, which does not appear to affect morphogenesis and overall tissue stresses, consistent with the lack of phenotype in cellular dynamics.

Results

Core PCP does not guide cellular dynamics during pupal blade elongation flows

To investigate the role of core PCP in orienting cell behaviors during pupal blade elongation flows, we analyzed cell dynamics in wild type (wt) and three different core PCP mutant tissues: prickle (pk³⁰, abbreviated as pk), strabismus (stbm⁶, abbreviated as stbm), and flamingo (fmi^frz3, aka stan^frz3, abbreviated as fmi). In pk, the core and Fat PCP systems remain aligned together toward the margin and the magnitude of Stbm polarity is reduced (Merkel et al., 2014). The mutants stbm and fmi are strong hypomorphs, where the core PCP network is strongly reduced (Figure 1—figure supplement 1B; Merkel et al., 2014). We analyzed shape changes of the wing blade during blade elongation flows and decomposed these changes into contributions from cell elongation changes and cell rearrangements, which include cell neighbor exchanges, cell divisions, cell extrusions, and correlation effects (Figure 1; Etournay et al., 2015; Merkel et al., 2017).

Figure 1 with 5 supplements see all

Download asset Open asset

Core planar cell polarity (PCP) does not orient cellular behaviors and tissue reshaping during pupal blade elongation flows: (A) Cartoon of pupal wing dissection at 16 hAPF and imaging using a spinning disk microscope.

(**B–B″**) Images of a wt wing at −7, 0, and 9 hRPCE (for this movie these times correspond to 16, 23, and 32 hAPF). The green and orange regions correspond to the hinge and blade, respectively. Anterior is up; proximal to the left. Scale bar, 100 μm. (C) Schematic of the cellular contributions underlying anisotropic tissue deformation. The tissue shear rate component $v_{x x}$ , which quantifies the rate of anisotropic tissue deformation along the proximal–distal wing axis, is decomposed into deformations arising from the rate of change of cell shapes $D Q_{x x} / D t$ and the deformations arising from the cellular rearrangements $R_{x x}$ (Etournay et al., 2015; Merkel et al., 2017). Total shear is the sum of cell elongation changes (green) and cell rearrangements (magenta). (**D–D′″**) Accumulated proximal–distal (Cum PD) tissue shear during blade elongation flows in the blade region averaged for (D) wt (n = 4), (D′) pk (n = 3), (D″) *stbm* (n = 3), and (**D′″**) *fmi* (n = 2) movies. Solid line indicates the mean, and the shaded regions enclose ± standard error of the mean (SEM). Differences in total accumulated shear are not statistically significant (Figure 1—figure supplement 2C). Time is relative to peak cell elongation (hRPCE).

Figure 1—source data 1 Numerical data for Figure 1D–D′′′, accumulated proximal–distal tissue shear during blade elongation flows in the blade region for wt and core PCP mutants.: https://cdn.elifesciences.org/articles/85581/elife-85581-fig1-data1-v3.zip
Download elife-85581-fig1-data1-v3.zip

In wt, the wing blade elongates along the PD axis (blue line in Figure 1D). Cells first elongate along the PD axis and then relax to more isotropic shapes (green line in Figure 1D). Cell rearrangements, however, go the opposite direction, initially contributing to AP deformation, before turning around to contribute to PD deformation (magenta line in Figure 1D). We introduce here a relative timescale, where we measure time in hours relative to the peak in cell elongation (hRPCE). This new scale allows us to handle variation in the timing of the onset of the blade elongation flows, which we have observed recently (see Appendix 1).

In core PCP mutants, we find that the dynamics of tissue shear, cell elongation changes, and cell rearrangements, when averaged across the entire blade, occur normally (Figure 1D–D′′′). We observe that by the end of the process, only slightly lower total shear appears to occur in the core PCP mutants, caused by slightly less cell rearrangements, but these subtle changes are not statistically significant (Figure 1—figure supplement 2C). The cellular dynamics contributing to isotropic tissue deformation are also broadly the same between wt and core PCP mutant tissues (Figure 1—figure supplement 2D, E). We also looked for differences in the behavior of regions of the wing blade subdivided along the PD axis (Figure 1—figure supplement 3), as previous work has shown that distal regions of the wing blade shear more at early times, whereas proximal regions start deforming later (Merkel et al., 2017). Again, we do not find strong differences between core PCP mutants and wt when we subdivide the wing into regions along the PD axis (Figure 1—figure supplement 3, Figure 1—figure supplement 4).

From this analysis, we conclude that core PCP is not required to determine the global patterns of cell dynamics during blade elongation flows. Interestingly, core PCP mutants do have a subtle but significant phenotype in the adult wing shape: pk and stbm (but not fmi) mutant wings are slightly rounder and wider than wt (Figure 1—figure supplement 5I). In principle, these small differences could arise after the blade elongation flows studied here. However, it is also possible that the we could not reliably detect these subtle differences in pupal wings due to the small number of wings per genotype that we were able to analyze (n = 2−4). To illustrate this point, we used the pool of adult wings ( $n = 53$ for wt, $n = 47$ for pk, $n = 74$ for stbm, and $n = 56$ for fmi), where the phenotype is significant, to understand the probability that a sample of smaller size $m$ would provide a significant signal, see Figure 1—figure supplement 5J. For $m = 3$ , corresponding to the number of pupal wings we analyzed, we find that only about 20% of samples show a significant phenotype. In other words, if the same magnitude of difference occurred during the blade elongation flows as in the adult, we would have only about a 20% chance to observe it. Therefore, core PCP could subtly influence the cell dynamics occurring at this stage. To investigate this possibility, we next looked for a possible difference in mechanical stresses in core PCP mutants.

A rheological model for the response to laser ablation

We investigated cell and tissue mechanics in core PCP mutants using laser ablation in a small region of the wing blade. We used a region located between the second and third sensory organs in the intervein region between the L3 and L4 longitudinal veins, which is a region that is easy to identify throughout blade elongation flows (Figure 2A). We cut three to four cells in a line along the AP axis and measured the displacement of the tissue (Figure 2A, Video 1). We calculate the initial recoil velocity $v$ by measuring the average displacement of ablated cell membranes at the first observed timepoint after the ablation, $δ t = 0.65 s$ (see Materials and methods, Linear laser ablations to calculate the initial recoil velocity). Previously, we reported that initial recoil velocity measured along the PD axis in wt peaks around −8 hRPCE (20 hAPF in Iyer et al., 2019), and therefore we first focus on this timepoint. We find that core PCP mutants have significantly lower initial recoil velocity (Figure 2B, Figure 2—figure supplement 1A), suggesting that there is a mechanical defect in these mutants.

Figure 2 with 1 supplement see all

Download asset Open asset

Rheological model for the response to laser ablation: (A) Schematic of a wt wing at −8 hRPCE.

Linear laser ablation experiments were performed in the blade region enclosed by the red square. Dots on the wing cartoon indicate sensory organs. The red line corresponds to the ablation, and the kymograph was drawn perpendicularly to the cut (yellow). Scale bar, 5 μm. (B) Initial recoil velocity upon ablation (simplified as recoil velocity in the y-axis title) along the proximal–distal (PD) axis at −8 hRPCE for wt (gray) and pk (blue) tissues (n ≥ 9). Significance is estimated using the Mann–Whitney U-test. ***p-val ≤0.001. (C) Example of the measured displacement after laser ablation (black dots) and corresponding exponential fit of the mechanical model (red curve). The blue and green regions highlight the displacement in the fast and slow timescale, respectively. (D) Description of the mechanical model that was devised to analyze the tissue response upon laser ablation. After the cut, the spring with elastic constant $k$ is ablated (red scissor), and the tissue response is given by the combination of the two Kelvin–Voigt models arranged in series. These two correspond to the fast response given by $k_{f}$ and $η_{f}$ and the slow response given by $k_{s}$ and $η_{s}$ . The mechanical stress $σ$ is constant. The membrane displacement $Δ x (t)$ is calculated as a sum of the displacement ( $X_{f}$ ) associated with the fast timescale ( $τ_{f}$ ) and the displacement ( $X_{s}$ ) associated with the slow timescale ( $τ_{s}$ ). (**E–E′′′**) Values obtained for each of the four fitting parameters when fit to the data. (E) Displacement associated with the fast and (E′) slow timescale for wt (gray) and pk (blue). (**E′′**) Fast and (**E′′′**) slow timescale for wt (gray) and pk (blue) (n ≥ 5). Significance is estimated using the Student’s t-test. **p-val ≤0.01; ns, p-val >0.05. (F) Example of a circular laser ablation used for analysis with elliptical shape after circular ablation (ESCA). The left image shows the final shape of the ablation around 2 min after cut, and the right image shows the corresponding segmented image, where the inner and outer pieces were fit with ellipses. After the fitting, the model outputs the anisotropic and isotropic stress (equations shown on the right side). Scale bar, 20 μm. A = anterior, P* = posterior, D = distal, P = proximal. (F′) Anisotropic stress $\tilde{σ} / 2 K$ for wt (gray) and pk (blue) tissues at −8 hRPCE (n ≥ 4). Significance is estimated using the Mann–Whitney U-test. ns, p-val >0.05. (F″) Isotropic stress $σ_{0} / \bar{K}$ for wt (gray) and pk (blue) tissues at −8 hRPCE (n ≥ 4). Significance is estimated using the Mann–Whitney U-test. ns, p-val >0.05. Time is relative to peak cell elongation (hRPCE). In all plots, each empty circle indicates one cut, and the box plots summarize the data: thick black line indicates the median; the boxes enclose the first and third quartiles; lines extend to the minimum and maximum without outliers, and filled circles mark outliers.

Figure 2—source data 1 Numerical data for Figure 2B, initial recoil velocity upon ablation along the proximal–distal (PD) axis for wt and pk tissues.: https://cdn.elifesciences.org/articles/85581/elife-85581-fig2-data1-v3.zip
Download elife-85581-fig2-data1-v3.zip
Figure 2—source data 2 Numerical data for Figure 2E, values for fitted parameters of the rheological model.: https://cdn.elifesciences.org/articles/85581/elife-85581-fig2-data2-v3.zip
Download elife-85581-fig2-data2-v3.zip
Figure 2—source data 3 Numerical data for Figure 2F–F″, values for anisotropic and isotropic stress deteremined with elliptical shape after circular ablation (ESCA).: https://cdn.elifesciences.org/articles/85581/elife-85581-fig2-data3-v3.zip
Download elife-85581-fig2-data3-v3.zip

Video 1

Download asset

posterframe for video — Shown here is an example of a linear laser ablation, cutting three to four cells, in wt (left) or pk pupal wings.

The movie goes dark during the ablation itself. Thereafter, the tissue displaces. Anterior is up; proximal is left.

As initial recoil velocity is often used as a proxy for mechanical stress (e.g. Mayer et al., 2010; Etournay et al., 2015; Iyer et al., 2019; Farhadifar et al., 2007), this result seems to suggest that the PCP mutant wings generate less mechanical stress during blade elongation flows, even though the cellular dynamics are at best only subtly perturbed. To explore this phenotype in more detail, we considered that the response to laser ablation is not exactly a direct measure of mechanical stress, as it is also affected by cellular material properties. We thus further analyzed the full kinetics of the linear laser ablations, focusing on the pk mutant, and developed a rheological model to interpret the results. When plotting displacement of the nearest bond to ablation over time, we realized that a single exponential relaxation cannot account for the observed behavior (Figure 2—figure supplement 1B, right). We obtained a good fit of the data by introducing a second relaxation timescale (Figure 2C). The slow timescale (∼20 s) accounts for most of the timecourse of displacement changes, but the fast timescale (<1 s) is required to account for first 5–10 datapoints, see Figure 2—figure supplement 1B, left. We therefore developed a model consisting of two Kelvin–Voigt (KV) elements in series (Figure 2D) to represent the tissue after ablation. The two KV elements have different elastic constants ( $k_{f}$ and $k_{s}$ ) and viscosities ( $η_{f}$ and $η_{s}$ ). Before ablation, the system is subjected to a constant stress ( $σ$ ) and contains a spring with elastic constant $k$ , which represents the cell patch that will be ablated. Upon ablation, the third spring is removed which leads to change in strain of our rheological model. We represent this strain by a displacement $Δ x$ as a function of time given by

Δ x (t) = X_{f} (1 - e^{- t / τ_{f}}) + X_{s} (1 - e^{- t / τ_{s}}),

where $X_{f} = σ κ / k_{f}$ is the displacement associated with the fast timescale, $τ_{f} = η_{f} / k_{f}$ , and $X_{s} = σ κ / k_{s}$ is the displacement associated with the slow timescale, $τ_{s} = η_{s} / k_{s}$ . Here, $κ = k / (k + \bar{k})$ is the fraction of the overall system elasticity lost due to ablation (see Materials and methods, Kymograph analysis and fit to model) and $\bar{k} = k_{s} k_{f} / (k_{s} + k_{f})$ is the elasticity of the two KV elements connected in series. With this model, we presume the properties of the ablated cell itself, including its membrane, adhesion proteins, and acto-myosin cortex likely dominate the fast timescale response. The slow timescale response is a collective effect emerging from the ablated cell together with its surrounding cellular network.

We analyzed the experimentally measured displacement over time for each ablation and then fit the data to our model with four parameters ( $X_{f}$ , $X_{s}$ , $τ_{f}$ , and $τ_{s}$ ) (Equation 1, Figure 2E–E′′′). Surprisingly, we find that the only parameter that changes between pk and wt is $X_{f}$ , the displacement associated with the fast timescale (Figure 2E–E′′′). To interpret this result, we consider that these four fitted parameters constrain the five mechanical model parameters (Figure 2D) but do not provide a unique solution. Since only one measured parameter changes, we asked what is the simplest set of model parameter changes that could have such an effect. To this end, we first note that the measured values of $X_{f}$ and $X_{s}$ ( $1.8 - 2.6 μ m$ vs $6 - 8 μ m$ , respectively) indicate $k_{f} ≫ k_{s}$ and therefore the overall elasticity of our rheological model is largely determined by the elasticity of the slow relaxation $\bar{k} \approx k_{s}$ . If we also consider that the contribution to the elasticity of the cellular patch from the ablated cells, represented by $k$ in the model, is small, then we can approximate $κ \approx k / k_{s}$ and therefore $X_{f} \approx σ k / (k_{s} k_{f})$ (see Materials and methods, Kymograph analysis and fit to model). Is the observed phenotype in the fast timescale displacement $X_{f}$ due to a change in tissue stress $σ$ or a change in the elastic constants?

To address this question, we sought to probe mechanical stress in the wt and pk mutant, independent of the ablation recoil velocity. To do so, we used a method called ESCA (elliptical shape after circular ablation) (Dye et al., 2021), which uses circular laser ablation and quantifies the resulting elliptical tissue outline once the mechanical equilibrium is established (Figure 2F and Materials and methods, Elliptical shape after circular ablation). Analysis of the elliptical tissue outline provides information about two-dimensional stresses present in the tissue before the ablation. In particular, we measure the magnitude of the anisotropic shear stress tensor, normalized by the shear elastic modulus $\tilde{σ} / (2 K)$ and the isotropic stress normalized by the area elastic modulus $σ_{0} / \bar{K}$ . The stress $σ$ in the simple rheological model presented above would correspond to tissue stress normal to the linear laser ablation axis and therefore it is a linear combination of both $\tilde{σ} / (2 K)$ and $σ_{0} / \bar{K}$ . ESCA also provides an estimate of the ratio of shear and area elastic constants $2 K / \bar{K}$ .

Using ESCA, we find no significant difference between wt and pk mutants in anisotropic and isotropic stress magnitudes, nor in the ratio of elastic constants (Figure 2F–F″ and Figure 2—figure supplement 1C). Since the ratio $σ / k_{s}$ defined in the rheological model is related to the normalized tissue stresses and elastic moduli, which do not change as shown by ESCA, we conclude that that $σ / k_{s}$ is not different between wt and pk. Therefore, we account for the observed changes of fast timescale displacement $X_{f}$ in the pk mutant with a change of the single elastic constant $k_{f}$ . In this scenario, $η_{f}$ changes together with the $k_{f}$ , such that $τ_{f} = η_{f} / k_{f}$ is conserved. This suggests that fast elasticity and viscosity are not independent but stem from a microscopic mechanism that controls the relaxation timescale. An example of such mechanism is turnover of the acto-myosin network, although this mechanism would be too slow to account for the fast relaxation timescale we observe. The conclusion that only the short-time response to the ablation, and not the tissue stress, is affected in the pk mutant is consistent with the lack of a clear phenotype in the large-scale tissue flows (Figure 1).

Dynamics of stress and cell elongation throughout blade elongation flows in wild type and core PCP mutants

To examine the effect of PCP mutation throughout blade elongation flows, we aimed to simplify the time intensive segmentation of the full ablation dynamics. To this end, we measured only the initial recoil velocity at different developmental timepoints. In terms of our model, the initial recoil velocity measured during first $δ t = 0.65 s$ can be expressed as $v = (X_{f} (1 - e^{- δ t / τ_{f}}) + X_{s} (1 - e^{- δ t / τ_{s}})) / δ t$ . Since the value of $δ t$ is comparable to the fast timescale $τ_{f}$ , about 63% of the $X_{f}$ value relaxes over that time, while at the same time only about 5% of the $X_{s}$ value is relaxed. Using the measured values of $X_{f}$ and $X_{s}$ , we estimate that the fast timescale dynamics contributes about 80% of the $v$ value. Therefore, the initial recoil velocity is a good proxy for the fast displacement $X_{f}$ .

We find that the initial recoil velocity along the PD axis peaks at −8 hRPCE before declining again by 4 hRPCE (Figure 3A), consistent with previous work (Iyer et al., 2019). The behavior of the initial recoil velocity in the pk mutant is qualitatively similar throughout blade elongation flows, however, with significantly lower magnitude than wt (Figure 3A). We also observed this behavior in stbm and fmi mutant tissues (Figure 3—figure supplement 1A). This result indicates that $X_{f}$ is lower in core PCP mutants than in wt throughout blade elongation flows.

Figure 3 with 2 supplements see all

Download asset Open asset

Dynamics of stress and cell elongation throughout blade elongation flows in wt and pk mutant.

(A) Initial recoil velocity upon ablation (simplified as recoil velocity in the y-axis title) along the proximal–distal (PD) axis throughout blade elongation flows for wt (gray) and pk (blue) tissues (n ≥ 3). Significance is estimated using the Mann–Whitney U-test. ****p-val ≤0.0001; ***p-val ≤0.001; **p-val ≤0.01; ns, p-val >0.05. (B) Elliptical shape after circular ablation (ESCA) results for anisotropic stress $\tilde{σ} / 2 K$ for wt (gray) and pk (blue) tissues throughout blade elongation flows (n ≥ 3). Significance is estimated using the Mann–Whitney U-test. *p-val <0.05; ns, p-val >0.05. (B′) ESCA results for isotropic stress $σ / \bar{K}$ for wt (gray) and pk (blue) throughout blade elongation flows (n ≥ 3). Significance is estimated using the Mann–Whitney U-test. *p-val <0.05; ns, p-val >0.05. (C) Color-coded PD component of cell elongation $Q$ in the blade region between the second and third sensory organs found in the intervein region between L2 and L3. The images correspond to wt (top row) and pk (bottom row) wings throughout blade elongation flows. Scale bar, 5 μm. (D) Quantification of the PD component of cell elongation $Q$ in this region throughout blade elongation flows for wt (gray) and pk (blue) (n ≥ 3). Significance is estimated using the Mann–Whitney U-test. ns, p-val >0.05. Time is relative to peak cell elongation (hRPCE). In all plots, each empty circle indicates one experiment, and the box plots summarize the data: thick black line indicates the median; the boxes enclose the first and third quartiles; lines extend to the minimum and maximum without outliers, and filled circles mark outliers.

Figure 3—source data 1 Numerical data for Figure 3A, initial recoil velocity upon ablation along the proximal–distal (PD) axis throughout blade elongation flows for wt and pk mutant tissues.: https://cdn.elifesciences.org/articles/85581/elife-85581-fig3-data1-v3.zip
Download elife-85581-fig3-data1-v3.zip
Figure 3—source data 2 Numerical data for Figure 3B–B′, elliptical shape after circular ablation (ESCA) results for anisotropic and isotropic stress in wt and pk mutant tissues throughout blade elongation flows.: https://cdn.elifesciences.org/articles/85581/elife-85581-fig3-data2-v3.zip
Download elife-85581-fig3-data2-v3.zip
Figure 3—source data 3 Numerical data for Figure 3D, proximal–distal (PD) component of cell elongation Q throughout blade elongation flows for wt and pk mutant tissues.: https://cdn.elifesciences.org/articles/85581/elife-85581-fig3-data3-v3.zip
Download elife-85581-fig3-data3-v3.zip

We also performed ESCA at different timepoints in pk mutants and observe that anisotropic stress ( $\tilde{σ} / 2 K$ ) rises early during blade elongation flows before eventually declining (Figure 3B), whereas isotropic stress ( $σ_{0} / \bar{K}$ ) remains fairly constant (Figure 3B′). Strikingly, ESCA does not report any difference in measured stresses between pk and wt, nor in the ratio of elastic constants ( $2 K / \bar{K}$ , Figure 3—figure supplement 1B) throughout blade elongation flows. To further compare the information contained in the initial recoil velocity with the anisotropic stress measured by ESCA, we performed linear ablations also in the perpendicular orientation. With such data, we could quantify the difference in initial recoil velocity between the two orientations $δ v = v_{P D} - v_{A P}$ , which is expected to be proportional to the shear stress along the PD wing axis. We then quantified how $δ v$ evolves throughout the blade elongation flows (Figure 3—figure supplement 1C, D). Whereas ESCA clearly shows that stresses in the tissue remain the same in wt and pk throughout blade elongation flows, the difference in initial recoil velocity $δ v$ is significantly lower in pk compared to wt. This result indicates that our conclusions based on the −8 hRPCE timepoint are true throughout blade elongation flows, namely that the differences in $X_{f}$ between wt and pk stem from the fast elastic constant $k_{f}$ and not from the differences in mechanical stresses in the tissue.

To further probe the possible role of core PCP in epithelial mechanics, we also measured the dynamics of the PD component of cell elongation ( $Q$ ) in wt and core PCP mutants (Figure 3C, Figure 3—figure supplement 2E). Interestingly, in both wt and pk, anisotropic stress peaks around −6 hRPCE (Figure 3B), whereas $Q$ peaks significantly later, between −4 and 0 hRPCE. We have previously related the tissue stress and cell elongation through a constitutive relation $\tilde{σ} = 2 K Q + ζ$ , where $ζ$ represented an active anisotropic stress component (Etournay et al., 2015). The difference in timing of the peaks in stress and cell elongation indicate that the active stresses change over time. However, we observe no differences between wt and core PCP in the peak of cell elongation (Figure 3D, Figure 3—figure supplement 2E), showing that core PCP also does not affect active anisotropic stresses underlying the dynamics of cell elongation during blade elongation flows.

Discussion

Here, we used the Drosophila pupal wing as a model for studying the interplay between planar polarized chemical signaling components, specifically the core PCP pathway, and the mechanical forces underlying tissue morphogenesis.

An extensive analysis of core PCP mutants shows no significant phenotype in pupal wing morphogenesis during the blade elongation flows. We find no significant differences in overall tissue shape change, nor in the pattern or dynamics of underlying cellular contributions. Even if a larger sample size of pupal wings would reveal a statistically significant phenotype, as indicated by our analysis of adult wings, the differences to the wild type would be subtle. Furthermore, we found no significant differences in tissue mechanical stress or in cell elongation over time. Generally these results are consistent in mutants that greatly reduce core PCP polarity (stbm and fmi) or prevent its decoupling from Fat (pk).

Interestingly, we do observe a phenotype in the initial recoil velocity upon laser ablation between core PCP mutants and wt, but this is not reflected in tissue stresses or large-scale morphogenetic flows that shape the wing. A detailed analysis of wt and pk suggests that the phenotype arises from a difference in the elastic constant $k_{f}$ underlying the fast timescale response ( $τ_{f} = 0.65 s$ ) to the ablation. In our simple model, the fast and stiff spring $k_{f}$ has a small contribution to the effective tissue elasticity $\bar{k}$ , which is dominated by the slow and soft spring $k_{s}$ (See Results, A rheological model for the response to laser ablation). The observation that core PCP only affects $k_{f}$ is therefore consistent with the lack of phenotype at larger scales. What is the the biophysical nature of the fast response to laser ablation? We hypothesize that processes that react on timescaless <1 s to a laser ablation could be related to cortical mechanics of cell bonds or possibly changes in cell hydraulics, and it is unclear how core PCP would affect these processes. Whether this core PCP phenotype in $k_{f}$ leads to very subtle changes in tissue development not detected in our analyses here, or is only visible in response to a laser ablation, also remains unknown and would require a much larger sample size to address. For the adult wing, we have a sufficient sample size to reveal a weak but significant shape phenotype in core PCP mutant wings. This result suggests that a weak phenotype arises during pupal development that we could not reliably detect in our analysis of cell dynamics. We also cannot rule out the possibility that there is a compensating mechanism that prevents the phenotype from appearing at larger scales.

Initial recoil velocity after a laser ablation is often used a proxy for tissue mechanical stresses. However, our results highlight a limitation of this approach for looking at how stress changes in different genotypes, as here we show how initial recoil velocity is influenced by differences in mechanics on small scales that are not necessarily related to differences in overall tissue stress.

While we have shown that core PCP is not required to organize the dynamic patterns of cellular events underlying blade elongation flows, it might still affect later stages of wing development. Furthermore, there may still be other patterning systems acting redundantly or independently with core PCP. For example, the Fat PCP system and Toll-like receptors have been shown to influence the orientation of cellular rearrangements and cell divisions in other contexts (Bosveld et al., 2012; Mao et al., 2006; Paré et al., 2014; Lavalou et al., 2021, reviewed in Umetsu, 2022). Whether and how other polarity systems influence pupal wing morphogenesis remains unknown. Alternatively, anisotropic mechanical stress induced by hinge contraction could itself provide a polarity cue through mechano-sensitive activity of the cytoskeleton. Our recent work in the larval wing disc shows that the cell polarity that drives the patterning of cell shape and mechanical stress contains a mechano-sensitive component (Dye et al., 2021). Here, we show a detailed analysis of tissues stress dynamics and cell elongation in the pupal wing revealing that the active cellular stresses that are relevant for pupal wing morphogenesis (Etournay et al., 2015) change in time (Figure 3). Whether the same mechano-sensitive mechanism established in the larval wing can also account for the dynamics of active stresses during the pupal blade elongation flows will be an important question to answer in the future.

Materials and methods

Key resources table

Reagent type (species) or resource	Designation	Source or reference	Identifiers	Additional information
gene (Drosophila melanogaster)	w-	NA	FLYB:FBal0018186
gene (Drosophila melanogaster)	shg (shotgun; E-cadherin)	NA	FLYB:FBgn0003391
gene (Drosophila melanogaster)	pk30	NA	FLYB:FBal0101223
gene (Drosophila melanogaster)	stbm6	NA	FLYB:FBal0062423
gene (Drosophila melanogaster)	fmifrz3	NA	FLYB:FBal0143193
strain, strain background (Drosophila melanogaster, male)	wt	Other	PMID:19429710	Huang et al., 2009. Genotype: wt: w-; EcadGFP;
strain, strain background (Drosophila melanogaster, male)	pk	Bloomington Drosophila Stock Center	RRID:BDSC_44229	Gubb et al., 1999. Genotype: w-; EcadGFP, pk30;
strain, strain background (Drosophila melanogaster, male)	stbm	Bloomington Drosophila Stock Center	RRID:BDSC_6918	Wolff and Rubin, 1998. Genotype: w-; EcadGFP, stbm6;
strain, strain background (Drosophila melanogaster, male)	fmi	Bloomington Drosophila Stock Center	RRID:BDSC_6967	Wolff and Rubin, 1998. Genotype: w-; EcadGFP, fmifrz3;
chemical compound, drug	Euparal	Carl Roth	7356.1
chemical compound, drug	Holocarbon oil 700	Sigma-Aldrich	H8898
chemical compound, drug	Isopropanol (2-propanol)	Sigma-Aldrich	1.0104
software, algorithm	Fiji	Other	v. 2.0.0-rc-68/1.52e	Schindelin et al., 2012
software, algorithm	Ilastik	Other	v. 1.2.2	Berg et al., 2019
software, algorithm	MATLAB	Other	v. 9.2.0.1226206 (R2017a)	MATLAB, 2017
software, algorithm	PreMosa	Other		Blasse et al., 2017
software, algorithm	R	Other	v. 3.4.1	R Development Core Team, 2020
software, algorithm	Rstudio	Other	v. 3.6.1	RStudio Team, 2020
software, algorithm	TissueMiner	Other	v. TM_1.0.2	Etournay et al., 2016
other	Coverslip	Paul Marienfeld GmbH	107052
other	Microscope slides	Paul Marienfeld GmbH	1000200
other	Dumont #55 Forceps	Fine Science Tools	11295–51
other	Vannas Spring Scissors	Fine Science Tools	15000–08

Fly husbandry

Request a detailed protocol

Flies were maintained at 25°C and fed with standard fly food containing cornmeal, yeast extract, soy flour, malt, agar, methyl 4-hydroxybenzoate, sugar beet syrup, and propionic acid. Flies were kept at 25°C in a 12-hr light/dark cycle. Vials were flipped every 2–3 days to maintain a continuous production of pupae and adult flies. All experiments were performed with male flies, since they are slightly smaller and therefore the wings require less tiling on the microscope to be imaged than females.

Long-term timelapse imaging of pupal wing morphogenesis

Acquisition

Request a detailed protocol

White male pupae were collected, slightly washed with a wet brush, and transferred to a vial containing standard food. At 16 hAPF, the pupal case was carefully dissected so that the wing would be accessible. The pupae was then mounted onto a 0.017-mm coverslip on a self-built metal dish with a drop of Holocarbon oil 700 (Classen et al., 2008). Pupal wing morphogenesis was imaged every 5 min for approximately 24 hr, as in Etournay et al., 2015. Wings that did not develop after 4–5 hr of imaging were discarded and not analyzed.

Two different microscopes were used for acquisition of long-term timelapses. All wt, pk, and stbm movies were acquired using a Zeiss spinning disk microscope driven by ZEN 2.6 (blue edition). This microscope consists of a motorized XYZ stage, an inverted stand, a Yokogawa CSU-X1 scan head, and a temperature-controlled chamber set to 25°C. The sample was illuminated with a 488-nm laser, and the emission was collected using a 470/40 bandpass filter, through a Zeiss 63 × 1.3 W/Gly LCI Plan-Neofluar objective and a Zeiss AxioCam Monochrome CCD camera with 2 × 2 binning. The whole wing was imaged in 24 tiles with an 8 % overlap. Each tile consisted of 50–60 stacks with a Z-spacing of 1 μm. The laser power was set to 0.1 mW.

The two fmi movies were acquired with an Olympus IX 83 inverted stand driven by the Andor iQ 3.6 software. The microscope is equipped with a motorized xyz stage, a Yokogawa CSU-W1 scan head, and an Olympus 60 × 1.3 Sil Plan SApo objective. The setup was located inside a temperature-controlled chamber set to 25°C. The sample was illuminated with a 488-nm laser, and the emission was collected using a 525/50 bandpass filter. The whole wing was imaged by tiling with eight tiles with a 10 % overlap. Each tile consisted of 50–60 stacks with a distance of 1 μm between them. The laser power was set to 0.75 mW.

Table 1 summarizes the date when the long-term timelapses were acquired and the age of the pupae during the imaging.

Table 1

Date of acquisition of all long-term timelapses.

Genotype	Date of acquisition	Start [hAPF]	End [hAPF]
wt	March 30, 2016	16	39.83
	April 2, 2016	16	36.58
	April 3, 2016	16	36.50
	April 13, 2016	16	32.58
	June 20, 2018	16	41.17
pk	April 9, 2016	16	39.00
	June 28, 2016	16	36.00
	June 29, 2016	16	39.92
stbm	November 25, 2015	16	40.67
	November 28, 2015	16	35.33
	December 11, 2014	16	37.58
fmi	October 20, 2018	16	39.17
fmi	July 20, 2019	16	38.17

Processing, segmentation, tracking, and database generation

Request a detailed protocol

Raw stacks were projected, corrected for illumination artifacts, and stitched using PreMosa (Blasse et al., 2017). The stitched images of individual timepoints were cropped to fit the wing size, registered using the Fiji plugin ‘Descriptor-based series registration (2D/3D + t)’, and converted to 8 bit with Fiji (Schindelin et al., 2012). The segmentation was performed with the Fiji plugin TissueAnalyzer (Schindelin et al., 2012; Aigouy et al., 2010; Aigouy et al., 2016). Segmentation errors were identified and manually corrected by looking at the cell divisions and deaths masks.

Subsequent processing and quantifications were performed using TissueMiner (Etournay et al., 2016). Before generating the relational database, we rotated the movies so that the angle formed by a manually drawn line connecting the sensory organs would be 0. We manually defined the regions of interest, such as the blade, hinge, and the anterior and posterior regions, using the last frame of the movie. Next, we generated the relational database containing information about the cellular dynamics during blade elongation flows using TissueMiner (Etournay et al., 2016).

We queried and worked with the data using the Dockerized version of RStudio (Nickoloff, 2016), which loads all packages and functions required to work with TissueMiner. Movies were aligned by the peak of cell elongation by fitting a quadratic function around the cell elongation values 40 frames before and after the absolute maximum of cell elongation in the blade region for each movie. The maximum of this curve was identified and set as the timepoint 0 hRPCE.

Adult wing preparation and analysis of wing shape

Adult male flies were fixed in isopropanol for at least 12 hr. One wing per fly was dissected in isopropanol, transferred to a microscope slide and covered with 50% euparal in isopropanol. Wings were mounted with 50–70 μl 75 % euparal/isopropanol.

wt, pk, and stbm wings were imaged using a Zeiss widefield Axioscan Z1 microscope equipped with a Zeiss 10 × 0.45 air objective. fmi wings were imaged using a Zeiss widefield Axiovert 200 M microscope equipped with a Zeiss 5 × 0.15 Plan-Neofluar air objective.

Wing blade parameters were quantified using a custom-written Fiji macro (provided as Source code 1) (Schindelin et al., 2012). The shape or major-to-minor ratio was calculated using a custom RStudio script (R Development Core Team, 2020; RStudio Team, 2020).

Subsampling and statistical analysis

Request a detailed protocol

Random sampling was done using a custom written RStudio pipeline (R Development Core Team, 2020; RStudio Team, 2020). A group of a given sample size was randomly selected with replacement for each group (wt, pk, stbm, and fmi), and a Kruskal–Wallis test was ran to compare them. This analysis was repeated 10,000 times. The sample sizes analyzed were 3, 4, 5, 6, 7, 8, 9, 10, 20, and 40.

Quantification of the PD component of cell elongation $Q$

Request a detailed protocol

Prior to all laser ablation experiments, we acquired a stack of 50 μm thick that was projected using PreMosa (Blasse et al., 2017). We cropped a region that enclosed the region that was ablated, segmented cells using TissueAnalyzer (Aigouy et al., 2010; Aigouy et al., 2016), and generated a relational database with TissueMiner (Etournay et al., 2016).

The definition of cell elongation was first presented in Aigouy et al., 2010 and it describes the angle and magnitude of the tensor. The cell elongation tensor is given by

(\begin{matrix} ϵ_{x x} & ϵ_{x y} \\ ϵ_{x y} & - ϵ_{x x} \end{matrix}),

where

ϵ_{x x} = \frac{1}{A_{c}} \int c o s (2 ϕ) d A

and

ϵ_{x y} = \frac{1}{A_{c}} \int s i n (2 ϕ) d A .

Cell elongation is normalized by the cell area $(A_{c})$ of each cell. The magnitude of cell elongation is:

ϵ = (ϵ_{x x}^{2} + ϵ_{x y}^{2})^{\frac{1}{2}}

Here, we plot $ϵ_{x x}$ as $Q$ , which we describe as the PD component of cell elongation.

Laser ablation experiments

Pupae were dissected and mounted as described for the long-term timelapses. Ablations were always performed in the same region of the wing blade, found in the intervein region between the longitudinal veins L3 and L4 and between the second and third sensory organs. This region was chosen because these landmarks are easily visible in all timepoints. Laser ablations were performed using a Zeiss spinning disk microscope equipped with a CSU-X1 Yokogawa scan head, an EMCCD Andor camera, a Zeiss 63 × 1.2 water immersion Korr UV-VIS-IR objective, and a custom-built laser ablation system using a 355-nm, 1000-Hz pulsed ultraviolet (UV) laser (Grill et al., 2001; Mayer et al., 2010). The imaging and cutting parameters for line and circular laser ablations are shown in Table 2. All laser ablation experiments were performed between January 2018 and July 2020, after the delay in pupal wing morphogenesis was identified.

Table 2

Parameters used to perform laser ablations.

	Linear ablations	Circular ablations
Exposure time [s]	0.05	0.05
488-nm laser intensity [%]	50	50
Time interval [s]	0.09	2.55
Pulses per shot	25	25
Shots per µm	2	2
Shooting time [s]	0.67	147.28
Thickness of stack ablated [µm]	1	20

Linear laser ablations to calculate the initial recoil velocity

Request a detailed protocol

We performed both types of linear ablations in only one plane of the tissue, in order to minimize the time required for ablation and therefore be able to acquire the initial recoil velocity upon ablation (no imaging is possible during ablation). The length of the linear laser ablations was 10 μm, ablating three to four cells. We drew kymographs perpendicularly to the cut to follow the two edges of one ablated cell using Fiji (Schindelin et al., 2012). The initial recoil velocity was calculated as the average displacement of two membranes of the same cell that occurred during the black frames of the ablation itself. This calculation was made using a self-written MATLAB script (MATLAB, 2017). Scripts used to make kymographs and analyze the laser ablations are provided in Source code 2; Source code 3; Source code 4. The image acquired prior to the laser ablation was used to compute $Q$ in that region, as described in Quantification of the PD component of cell elongation Q, and the time corresponding to the maximum of cell elongation was defined as 0 hRPCE.

Elliptical shape after circular ablation

Request a detailed protocol

Circular laser ablations used for ESCA were 20 μm in radius (approximately 10 cells). This radius was selected such that it would fit into the same blade region throughout the blade elongation flows. Due to the bigger size of these cuts and the curvature of the tissue, we cut the tissue along a stack of 20 μm thick. Approximately 2 min after the ablation, we acquired a stack of 50 μm. This image was projected using PreMosa (Blasse et al., 2017) and preprocessed by applying Gausian blur (σ = 1) and background subtraction filters (rolling ball radius = 30) in Fiji (Schindelin et al., 2012). The next steps were performed as in Dye et al., 2021: the image of the final shape of the cut was segmented using Ilastik (Berg et al., 2019) by defining three regions: membrane, cell, and dark regions. The segmented image was thresholded to obtain a binary image of the final shape of the cut. We fitted two ellipses to this image: one to the inner piece and another one to the outer outline of the cut. Based on the shape of these ellipses, the method outputs the anisotropic $\frac{\tilde{σ}}{2 K}$ and isotropic stress $\frac{σ}{\bar{K}}$ as a function of their respective elastic constants, and the ratio of elastic constants $\frac{2 K}{\bar{K}}$ . A small number of experiments were fitted poorly (defined as an error per point greater than 0.3) and were therefore excluded from analysis. Prior to the circular ablation, a stack of 50 μm was acquired and used to calculate cell elongation before ablation (Quantification of the PD component of cell elongation Q). The time corresponding to the maximum of cell elongation was set to be 0 hRPCE.

Kymograph analysis and fit to model

Request a detailed protocol

The ablations used to calculate the mechanical stress along the PD axis for wt and pk were further analyzed with the rheological model. To do so, we processed the kymographs by applying a Gaussian blur (σ = 1) (Schindelin et al., 2012), and then we segmented these kymographs with Ilastik (Berg et al., 2019). Using a self-written Fiji macro (Schindelin et al., 2012), we extracted the intensity profile for each timepoint. Next, we wrote an R script (R Development Core Team, 2020; RStudio Team, 2020) to identify the membrane displacement over time and obtained a unique curve per kymograph, which could be fitted with our model. We modeled a local patch of tissue as a combination of a spring with spring constant $k$ , representing the ablated cells, and two KV elements with spring constants $k_{f}$ and $k_{s}$ and viscosity coefficients $η_{f}$ and $η_{s}$ , representing the unablated cells, as shown in Figure 2C, D. Because the local tissue strain in the experimental measurement is expressed by the displacement of the bond nearest to the ablation, in the rheological model we represent tissue strain by displacements of the two KV elements. In principle, the strain can be recovered by normalizing the displacements by the width of ablated cells. Displacements of the two KV elements are defined as a change in the distance between the end points of the KV elements $x_{i} (t)$ , relative to their initial values $x_{i} (0)$ , where $i \in {f, s}$ for fast ( $f$ ) and slow ( $s$ ) element.

Mechanical stress in the tissue is represented by the $σ$ acting on our model, and we assume that $σ$ is not changed by the ablation. Before the ablation, the model is in mechanical equilibrium and we can write

σ = (k + \bar{k}) x (0),

where $x (0)$ is the initial distance between the two end points of the model, and $\bar{k} = k_{f} k_{s} / (k_{f} + k_{s})$ is the elastic constant of the two KV elements connected in series. Upon ablation, the spring $k$ is removed and stresses in the model are imbalanced. The distance between the end points of the model $x (t)$ then evolves toward the new equilibrium position. The distance $x (t)$ can be decomposed as $x (t) = x_{f} (t) + x_{s} (t)$ , where $x_{f} (t)$ and $x_{s} (t)$ are the time-dependent distances between end points of the two KV elements, representing their strains. The dynamics of $x (t)$ is then obtained by writing the force balance equation for the two KV elements

\begin{aligned} σ & = k_{f} x_{f} (t) + η_{f} \frac{d x_{f} (t)}{d t}, \end{aligned}

\begin{aligned} σ & = k_{s} x_{s} (t) + η_{s} \frac{d x_{s} (t)}{d t}, \end{aligned}

We solve for $x_{f} (t)$ and $x_{s} (t)$ to obtain

x_{f} (t) = \frac{σ}{k_{f}} (1 - e^{- t / τ_{f}}) + x_{f} (0) e^{- t / τ_{f}},

x_{s} (t) = \frac{σ}{k_{s}} (1 - e^{- t / τ_{s}}) + x_{s} (0) e^{- t / τ_{s}},

where

\begin{aligned} x_{f, s} (0) & = \frac{σ (1 - κ)}{k_{f, s}}, \end{aligned}

and $κ = k / (k + \bar{k})$ is the fraction of the overall model elasticity $k + \bar{k}$ destroyed by the ablation. The displacement relative to the initial configuration $Δ x (t) = x (t) - x (0)$ is therefore

\begin{aligned} Δ x (t) = X_{f} (1 - e^{- t / τ_{f}}) + X_{s} (1 - e^{- t / τ_{s}}), \end{aligned}

where we introduced the long time displacements associated with the two KV elements

\begin{aligned} X_{f, s} & = \frac{σ κ}{k_{f, s}} . \end{aligned}

For simplicity, in the main text we refer to the long time displacements $X_{f}$ and $X_{s}$ of the two KV elements simply as displacements.

Statistical analysis

Request a detailed protocol

Statistical analysis was done using R (R Development Core Team, 2020; RStudio Team, 2020). We first tested normality of the data using the Shapiro–Wilk test. When data were normal, we used Student’s t-test to test statistical significance between two groups and analysis of variance test for multiple groups. When data were not normally distributed, significance was tested using the Mann–Whitney U-test for two groups and Kruskal–Wallis test for multiple groups. Statistical test results are shown on the figure captions.

Appendix 1

During the course of this work, we identified a delay in the onset of blade elongation flows compared to previous work (Appendix 1—figure 1A; Etournay et al., 2015; Piscitello-Gómez et al., 2023). In the past, cells reached their maximum of cell elongation at 22.9 ± 0.4 hAPF, while now they reach it at 28 hAPF. Although we do not know the cause of this delay, we have ruled out differences in temperature (either on the microscope or during development), nutrition (plant vs yeast-based foods), genetic background, presence of the parasite Wolbachia, or circadian gating (Piscitello-Gómez et al., 2023). To deal with this variation and combine data acquired over the years, we present cell dynamics data aligned in time by the peaks of cell elongation, and we refer to this timepoint as 0 hRPCE (relative to peak cell elongation) (Appendix 1—figure 1A). We investigated the cell dynamics underlying blade elongation flows in the delayed flies and observed that the shear rates were comparable with the older flies (Appendix 1—figure 1B). Thus, it is reasonable to shift the curves by aligning them to a new reference time.

Appendix 1—figure 1

Download asset Open asset

Delay and time alignment of old and newer flies: (A) Left: Snapshots of the blade region of long-term timelapses of wt pupal wing morphogenesis acquired in different years.

Scale bar, 10 μm. Right: Cell elongation norm during blade elongation flows for old flies (orange palette, 2016 flies) and new flies (green curve, 2018 fly). Top plot: Cell elongation magnitude for each movie not aligned in time. The peak of cell elongation is delayed from around 23 to 28 hAPF. Bottom plot: Cell elongation magnitude after alignment in time to the peak of cell elongation. Time is now expressed in hours relative to peak cell elongation (hRPCE). (B) Cell dynamics underlying anisotropic tissue deformation for 2016 wt (n = 4, top left), 2018 wt flies (n = 1, top right), and 2 *fmi* flies imaged in 2018 (bottom left) and 2019 (bottom right). The vertical dashed line marks the timepoint where cell rearrangements flip from AP- to PD-oriented per movie. The two dotted lines mark the start and the end of the analyzed wt long-term timelapses acquired in 2016. The time is relative to the peak of cell elongation (hRPCE).

Data availability

Source data and code are provided for each figure.

References

Book
1. Adler PN
(2012)
The frizzled/Stan pathway and planar cell polarity in the Drosophila wing

In: Adler PN, editors. Planar Cell Polarity During Development. Elsevier. pp. 1–31.
- Google Scholar
1. Aigouy B
2. Farhadifar R
3. Staple DB
4. Sagner A
5. Röper JC
6. Jülicher F
7. Eaton S
(2010) Cell flow reorients the axis of planar polarity in the wing epithelium of Drosophila
Cell 142:773–786.
https://doi.org/10.1016/j.cell.2010.07.042
- PubMed
- Google Scholar
Book
(2016)
Segmentation and quantitative analysis of epithelial tissues

In: Aigouy B, Umetsu D, editors. Drosophila. Springer. pp. 227–239.
- Google Scholar
1. Berg S
2. Kutra D
3. Kroeger T
4. Straehle CN
5. Kausler BX
6. Haubold C
7. Schiegg M
8. Ales J
9. Beier T
10. Rudy M
11. Eren K
12. Cervantes JI
13. Xu B
14. Beuttenmueller F
15. Wolny A
16. Zhang C
17. Koethe U
18. Hamprecht FA
19. Kreshuk A
(2019) ilastik: interactive machine learning for (bio)image analysis
Nature Methods 16:1226–1232.
https://doi.org/10.1038/s41592-019-0582-9
- PubMed
- Google Scholar
1. Blasse C
2. Saalfeld S
3. Etournay R
4. Sagner A
5. Eaton S
6. Myers EW
(2017) PreMosa: extracting 2D surfaces from 3D microscopy mosaics
Bioinformatics 33:2563–2569.
https://doi.org/10.1093/bioinformatics/btx195
- PubMed
- Google Scholar
1. Bosveld F
2. Bonnet I
3. Guirao B
4. Tlili S
5. Wang Z
6. Petitalot A
7. Marchand R
8. Bardet P-L
9. Marcq P
10. Graner F
11. Bellaïche Y
(2012) Mechanical control of morphogenesis by Fat/Dachsous/Four-jointed planar cell polarity pathway
Science 336:724–727.
https://doi.org/10.1126/science.1221071
- PubMed
- Google Scholar
1. Butler MT
2. Wallingford JB
(2017) Planar cell polarity in development and disease
Nature Reviews. Molecular Cell Biology 18:375–388.
https://doi.org/10.1038/nrm.2017.11
- PubMed
- Google Scholar
1. Capilla A
2. Johnson R
3. Daniels M
4. Benavente M
5. Bray SJ
6. Galindo MI
(2012) Planar cell polarity controls directional Notch signaling in the Drosophila leg
Development 139:2584–2593.
https://doi.org/10.1242/dev.077446
- PubMed
- Google Scholar
1. Chung S
2. Vining MS
3. Bradley PL
4. Chan CC
5. Wharton KA
6. Andrew DJ
(2009) Serrano (sano) functions with the planar cell polarity genes to control tracheal tube length
PLOS Genetics 5:e1000746.
https://doi.org/10.1371/journal.pgen.1000746
- PubMed
- Google Scholar
(2005) Hexagonal packing of Drosophila wing epithelial cells by the planar cell polarity pathway
Developmental Cell 9:805–817.
https://doi.org/10.1016/j.devcel.2005.10.016
- PubMed
- Google Scholar
Book
(2008)
Imaging Drosophila Pupal wing Morphogenesis

In: Classen AK, Aigouy B, editors. Drosophila. Springer. pp. 265–275.
- Google Scholar
1. Deans MR
(2021) Conserved and divergent principles of planar polarity revealed by hair cell development and function
Frontiers in Neuroscience 15:742391.
https://doi.org/10.3389/fnins.2021.742391
- PubMed
- Google Scholar
1. Devenport D
(2014) The cell biology of planar cell polarity
The Journal of Cell Biology 207:171–179.
https://doi.org/10.1083/jcb.201408039
- PubMed
- Google Scholar
(2021) Self-organized patterning of cell morphology via mechanosensitive feedback
eLife 10:e57964.
https://doi.org/10.7554/eLife.57964
- PubMed
- Google Scholar
1. Etournay R
2. Popović M
3. Merkel M
4. Nandi A
5. Blasse C
6. Aigouy B
7. Brandl H
8. Myers G
9. Salbreux G
10. Jülicher F
11. Eaton S
(2015) Interplay of cell dynamics and epithelial tension during morphogenesis of the Drosophila pupal wing
eLife 4:e07090.
https://doi.org/10.7554/eLife.07090
- PubMed
- Google Scholar
1. Etournay R
2. Merkel M
3. Popović M
4. Brandl H
5. Dye NA
6. Aigouy B
7. Salbreux G
8. Eaton S
9. Jülicher F
(2016) TissueMiner: A multiscale analysis toolkit to quantify how cellular processes create tissue dynamics
eLife 5:e14334.
https://doi.org/10.7554/eLife.14334
- PubMed
- Google Scholar
(2007) The influence of cell mechanics, cell-cell interactions, and proliferation on epithelial packing
Current Biology 17:2095–2104.
https://doi.org/10.1016/j.cub.2007.11.049
- PubMed
- Google Scholar
(2021) Tissue fluidity mediated by adherens junction dynamics promotes planar cell polarity-driven ommatidial rotation
Nature Communications 12:6974.
https://doi.org/10.1038/s41467-021-27253-0
- PubMed
- Google Scholar
1. Gho M
2. Schweisguth F
(1998) Frizzled signalling controls orientation of asymmetric sense organ precursor cell divisions in Drosophila
Nature 393:178–181.
https://doi.org/10.1038/30265
- PubMed
- Google Scholar
(2001) Polarity controls forces governing asymmetric spindle positioning in the Caenorhabditis elegans embryo
Nature 409:630–633.
https://doi.org/10.1038/35054572
- PubMed
- Google Scholar
1. Gubb D
2. Green C
3. Huen D
4. Coulson D
5. Johnson G
6. Tree D
7. Collier S
8. Roote J
(1999) The balance between isoforms of the Prickle LIM domain protein is critical for planar polarity in Drosophila imaginal discs
Genes & Development 13:2315–2327.
https://doi.org/10.1101/gad.13.17.2315
- Google Scholar
1. Guirao B
2. Rigaud SU
3. Bosveld F
4. Bailles A
5. López-Gay J
6. Ishihara S
7. Sugimura K
8. Graner F
9. Bellaïche Y
(2015) Unified quantitative characterization of epithelial tissue development
eLife 4:e08519.
https://doi.org/10.7554/eLife.08519
- PubMed
- Google Scholar
1. Huang J
2. Zhou W
3. Dong W
4. Watson AM
5. Hong Y
(2009) From the Cover: Directed, efficient, and versatile modifications of the Drosophila genome by genomic engineering
PNAS 106:8284–8289.
https://doi.org/10.1073/pnas.0900641106
- PubMed
- Google Scholar
(2019) Epithelial Viscoelasticity Is Regulated by Mechanosensitive E-cadherin Turnover
Current Biology 29:578–591.
https://doi.org/10.1016/j.cub.2019.01.021
- PubMed
- Google Scholar
Book
1. Jenny A
(2010)
Planar cell polarity signaling in the Drosophila eye

In: Jenny A, editors. Current Topics in Developmental Biology. Elsevier. pp. 189–227.
- Google Scholar
1. Lavalou J
2. Mao Q
3. Harmansa S
4. Kerridge S
5. Lellouch AC
6. Philippe JM
7. Audebert S
8. Camoin L
9. Lecuit T
(2021) Formation of polarized contractile interfaces by self-organized Toll-8/Cirl GPCR asymmetry
Developmental Cell 56:1574–1588.
https://doi.org/10.1016/j.devcel.2021.03.030
- PubMed
- Google Scholar
1. Mao Y
2. Rauskolb C
3. Cho E
4. Hu WL
5. Hayter H
6. Minihan G
7. Katz FN
8. Irvine KD
(2006) Dachs: an unconventional myosin that functions downstream of Fat to regulate growth, affinity and gene expression in Drosophila
Development 133:2539–2551.
https://doi.org/10.1242/dev.02427
- PubMed
- Google Scholar
1. Matis M
2. Axelrod JD
(2013) Regulation of PCP by the Fat signaling pathway
Genes & Development 27:2207–2220.
https://doi.org/10.1101/gad.228098.113
- PubMed
- Google Scholar
Software
1. MATLAB
(2017) The Mathworks, Inc, version 9.2.0.556344
MATLAB (R2017a).

https://www.mathworks.com/company/newsroom/mathworks-announces-release-2017a-of-the-matlab-and-simulink-pro.html
1. Mayer M
2. Depken M
3. Bois JS
4. Jülicher F
5. Grill SW
(2010) Anisotropies in cortical tension reveal the physical basis of polarizing cortical flows
Nature 467:617–621.
https://doi.org/10.1038/nature09376
- PubMed
- Google Scholar
1. Merkel M
2. Sagner A
3. Gruber FS
4. Etournay R
5. Blasse C
6. Myers E
7. Eaton S
8. Jülicher F
(2014) The balance of prickle/spiny-legs isoforms controls the amount of coupling between core and fat PCP systems
Current Biology 24:2111–2123.
https://doi.org/10.1016/j.cub.2014.08.005
- PubMed
- Google Scholar
1. Merkel M
2. Etournay R
3. Popović M
4. Salbreux G
5. Eaton S
6. Jülicher F
(2017) Triangles bridge the scales: Quantifying cellular contributions to tissue deformation
Physical Review. E 95:032401.
https://doi.org/10.1103/PhysRevE.95.032401
- PubMed
- Google Scholar
1. Mirkovic I
2. Gault WJ
3. Rahnama M
4. Jenny A
5. Gaengel K
6. Bessette D
7. Gottardi CJ
8. Verheyen EM
9. Mlodzik M
(2011) Nemo kinase phosphorylates β-catenin to promote ommatidial rotation and connects core PCP factors to E-cadherin-β-catenin
Nature Structural & Molecular Biology 18:665–672.
https://doi.org/10.1038/nsmb.2049
- PubMed
- Google Scholar
Book
1. Nickoloff J
(2016)
Docker in Action

Manning Publications Co.
- Google Scholar
1. Paré AC
2. Vichas A
3. Fincher CT
4. Mirman Z
5. Farrell DL
6. Mainieri A
7. Zallen JA
(2014) A positional Toll receptor code directs convergent extension in Drosophila
Nature 515:523–527.
https://doi.org/10.1038/nature13953
- PubMed
- Google Scholar
Preprint
(2023) Sensitivity of the Timing of Drosophila Pupal Wing Morphogenesis to External Perturbations
bioRxiv.
https://doi.org/10.1101/2023.09.22.558969
- Google Scholar
(2015) Patterned anchorage to the apical extracellular matrix defines tissue shape in the developing appendages of Drosophila
Developmental Cell 34:310–322.
https://doi.org/10.1016/j.devcel.2015.06.019
- PubMed
- Google Scholar
Software
1. R Development Core Team
(2020) R: A language andEnvironment for statistical computing
R Foundation for Statistical Computing, Vienna, Austria.

https://www.R-project.org
Software
1. RStudio Team
(2020) Rstudio: integrated development environment for R
R Packages.

https://www.r-project.org/conferences/useR-2011/abstracts/180111-allairejj.pdf
1. Sagner A
2. Merkel M
3. Aigouy B
4. Gaebel J
5. Brankatschk M
6. Jülicher F
7. Eaton S
(2012) Establishment of global patterns of planar polarity during growth of the Drosophila wing epithelium
Current Biology 22:1296–1301.
https://doi.org/10.1016/j.cub.2012.04.066
- PubMed
- Google Scholar
1. Schindelin J
2. Arganda-Carreras I
3. Frise E
4. Kaynig V
5. Longair M
6. Pietzsch T
7. Preibisch S
8. Rueden C
9. Saalfeld S
10. Schmid B
11. Tinevez JY
12. White DJ
13. Hartenstein V
14. Eliceiri K
15. Tomancak P
16. Cardona A
(2012) Fiji: an open-source platform for biological-image analysis
Nature Methods 9:676–682.
https://doi.org/10.1038/nmeth.2019
- PubMed
- Google Scholar
1. Ségalen M
2. Johnston CA
3. Martin CA
4. Dumortier JG
5. Prehoda KE
6. David NB
7. Doe CQ
8. Bellaïche Y
(2010) The Fz-Dsh planar cell polarity pathway induces oriented cell division via Mud/NuMA in Drosophila and zebrafish
Developmental Cell 19:740–752.
https://doi.org/10.1016/j.devcel.2010.10.004
- PubMed
- Google Scholar
(2016) Measuring forces and stresses in situ in living tissues
Development 143:186–196.
https://doi.org/10.1242/dev.119776
- PubMed
- Google Scholar
1. Umetsu D
(2022) Cell mechanics and cell-cell recognition controls by Toll-like receptors in tissue morphogenesis and homeostasis
Fly 16:233–247.
https://doi.org/10.1080/19336934.2022.2074783
- PubMed
- Google Scholar
1. Wolff T
2. Rubin GM
(1998) Strabismus, a novel gene that regulates tissue polarity and cell fate decisions in Drosophila
Development 125:1149–1159.
https://doi.org/10.1242/dev.125.6.1149
- PubMed
- Google Scholar
(1995) Frizzled regulates mirror-symmetric pattern formation in the Drosophila eye
Development 121:3045–3055.
https://doi.org/10.1242/dev.121.9.3045
- PubMed
- Google Scholar

Article and author information

Author details

Romina Piscitello-Gómez
1. Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany
2. DFG Excellence Cluster Physics of Life, Technische Universität Dresden, Dresden, Germany
Contribution
Data curation, Formal analysis, Investigation, Software, Validation, Visualization, Writing – original draft, Writing – review and editing

Contributed equally with
Franz S Gruber

Competing interests
No competing interests declared
Franz S Gruber
1. Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany
2. National Phenotypic Screening Centre, University of Dundee, Dundee, United Kingdom
Contribution
Data curation, Formal analysis, Investigation, Software, Validation, Writing – review and editing

Contributed equally with
Romina Piscitello-Gómez

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-2008-8460
Abhijeet Krishna
1. Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany
2. DFG Excellence Cluster Physics of Life, Technische Universität Dresden, Dresden, Germany
3. Center for Systems Biology Dresden, Dresden, Germany
Contribution
Formal analysis, Methodology, Software, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-9291-500X
Charlie Duclut
1. Laboratoire Physico-Chimie Curie, CNRS UMR 168, Institut Curie, Université PSL, Sorbonne Université, Paris, France
2. Max Planck Institute for Physics of Complex Systems, Dresden, Germany
3. Université Paris Cité, Laboratoire Matière et Systèmes Complexes, Paris, France
Contribution
Formal analysis, Software

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-8595-6815
Carl D Modes
1. Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany
2. DFG Excellence Cluster Physics of Life, Technische Universität Dresden, Dresden, Germany
3. Center for Systems Biology Dresden, Dresden, Germany
Contribution
Validation, Investigation, Formal analysis, Software

Competing interests
No competing interests declared
Marko Popović
1. DFG Excellence Cluster Physics of Life, Technische Universität Dresden, Dresden, Germany
2. Center for Systems Biology Dresden, Dresden, Germany
3. Max Planck Institute for Physics of Complex Systems, Dresden, Germany
Contribution
Investigation, Formal analysis, Data curation, Software

Competing interests
No competing interests declared
Frank Jülicher
1. DFG Excellence Cluster Physics of Life, Technische Universität Dresden, Dresden, Germany
2. Center for Systems Biology Dresden, Dresden, Germany
3. Max Planck Institute for Physics of Complex Systems, Dresden, Germany
Contribution
Visualization, Validation, Investigation, Writing – original draft, Formal analysis, Writing – review and editing, Software

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-4731-9185
Natalie A Dye
1. Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany
2. DFG Excellence Cluster Physics of Life, Technische Universität Dresden, Dresden, Germany
3. Mildred Scheel Nachwuchszentrum P2, Medical Faculty, Technische Universität Dresden, Dresden, Germany
Contribution
Data curation, Supervision, Writing – original draft, Writing – review and editing

For correspondence
natalie_anne.dye@tu-dresden.de

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-4859-6670
Suzanne Eaton (deceased)
1. Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany
2. DFG Excellence Cluster Physics of Life, Technische Universität Dresden, Dresden, Germany
3. Center for Systems Biology Dresden, Dresden, Germany
4. Biotechnologisches Zentrum, Technische Universität Dresden, Dresden, Germany
Contribution
Visualization, Validation, Investigation, Writing – original draft, Formal analysis, Writing – review and editing

Competing interests
No competing interests declared

Funding

Max Planck Society

Franz S Gruber
Abhijeet Krishna
Charlie Duclut
Carl D Modes
Frank Jülicher
Natalie A Dye
Suzanne Eaton
Romina Piscitello-Gómez
Marko Popović

Deutsche Forschungsgemeinschaft (EXC-2068-390729961)

Natalie A Dye

Deutsche Forschungsgemeinschaft (SPP1782/EA4/10-1)

Suzanne Eaton
Natalie A Dye
Romina Piscitello-Gómez
Franz S Gruber

Deutsche Krebshilfe (MSNZ-P2 Dresden)

Natalie A Dye

Austrian Academy of Sciences (DOC Fellowship)

Franz S Gruber

Agence Nationale de la Recherche (ANR-11-LABX-0071)

Charlie Duclut

Agence Nationale de la Recherche (ANR-18-IDEX-0001)

Charlie Duclut

Deutsche Forschungsgemeinschaft (SPP1782/EA4/10-2)

Suzanne Eaton
Natalie A Dye
Romina Piscitello-Gómez

Open access funding provided by Max Planck Society. The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Stephan Grill for giving us access to the microscope used for laser ablation. We thank the Light Microscopy Facility, the Computer Department, and the Fly Keepers of the MPI-CBG for their support and expertise. We would like to thank Christian Dahmann and Jana Fuhrmann for comments on the manuscript prior to publication. This work was funded by Germany’s Excellence Strategy – EXC-2068-390729961– Cluster of Excellence Physics of Life of TU Dresden, as well as grants awarded to SE from the Deutsche Forschungsgemeinschaft (SPP1782, EA4/10-1, EA4/10-2) and core funding of the Max-Planck Society to SE and NAD. NAD additionally acknowledges funding from the Deutsche Krebshilfe (MSNZ P2 Dresden). AK and RPG were funded through the Elbe PhD program. FSG was supported by a DOC Fellowship of the Austrian Academy of Sciences. CD acknowledges the support of a postdoctoral fellowship from the LabEx 'Who Am I?' (ANR-11-LABX-0071) and the Université Paris Cité IdEx (ANR-18-IDEX-0001) funded by the French Government through its 'Investments for the Future'. We dedicate this work to our coauthor Prof. Dr. Suzanne Eaton, who tragically passed away before the finalization of the project.

Version history

Preprint posted: December 10, 2022 (view preprint)
Received: December 15, 2022
Accepted: December 18, 2023
Accepted Manuscript published: December 20, 2023 (version 1)
Version of Record published: February 2, 2024 (version 2)
Version of Record updated: February 5, 2024 (version 3)

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.