Minimal vertex model explains how the amnioserosa avoids fluidization during Drosophila dorsal closure

Abstract

Dorsal closure is a process that occurs during embryogenesis of Drosophila melanogaster. During dorsal closure, the amnioserosa (AS), a one-cell thick epithelial tissue that fills the dorsal opening, shrinks as the lateral epidermis sheets converge and eventually merge. During this process, the aspect ratio of amnioserosa cells increases markedly. The standard 2-dimensional vertex model, which successfully describes tissue sheet mechanics in multiple contexts, would in this case predict that the tissue should fluidize via cell neighbor changes. Surprisingly, however, the amnioserosa remains an elastic solid with no such events. We here present a minimal extension to the vertex model that explains how the amnioserosa can achieve this unexpected behavior. We show that continuous shrink-age of the preferred cell perimeter and cell perimeter polydispersity lead to the retention of the solid state of the amnioserosa. Our model accurately captures measured cell shape and orientation changes and predicts non-monotonic junction tension that we confirm with laser ablation experiments.

Keywords: Drosophila dorsal closure; amnioserosa; epithelial tissue; morphogenesis; tissue rigidity; vertex model.

Fig. 1.. Experiment and vertex model for dorsal closure.

(A) The geometry of the dorsal hole during early (left), middle (center), and late (right) dorsal closure. Enlargements show tissue with selected cells, several of which ingress (highlighted by triangles). (B) We model the dorsal closure process as a quasistatic uni-axial deformation. The geometry of the model is shown at the beginning (left), in the middle (center, at 45% closure), and towards the end (right, 80% closure) of the process. $\Delta A\left(t\right)=\frac{A_0-A\left(t\right)}{A_0}$ is the fractional change in total AS area of the closure process, where $A_0$ is the AS area at the onset of dorsal closure. (C) An initial normal distribution of the preferred shape index of the model tissue (dashed red) with the standard deviation adjusted to be 0.45, leads to a distribution of the actual shape index after minimization (solid blue) that is in excellent agreement with the distribution of the experimentally observed shape index (solid black) at the beginning of dorsal closure. (D) Sketch of AS tissue regions included in model comparison (white center), with edge regions excluded (gray regions). (E) In the model, we reduce the preferred cell perimeter at a linear rate (blue) to capture the experimentally observed decrease of junction lengths (black). For comparison, we normalize the average perimeter by its value at the onset of the process. Inset: schematic representation of the reduction of cellular junction length and apical area during dorsal closure.

Fig. 2.. Results from experiment (black solid) and model (blue dashed).

(A) A comparison of average shape index $\overline{q}$ as a function of $\Delta A\left(t\right)=\frac{A_0-A\left(t\right)}{A_0}$. Here $A_0$ is the AS area at the onset of dorsal closure and $A\left(t\right)$ is the area as it shrinks during dorsal closure, so that $\Delta A\left(t\right)=0$ at onset. Lines show the mean, and shaded regions indicate the standard deviation among 12 embryos (experiment) or 10 different initial configurations (model). (B) Comparison of cell to cell standard deviation of the shape index $({\sigma }_q)$ during dorsal closure. (C) Orientational order parameter $\left(\overline{Q}\right)$ of the cells during dorsal closure. (D) Experimental initial junction recoil velocity (left y-axis) of the vertices after performing laser ablation of the junction, and predicted average cellular cortical tension $\left({\overline{\tau }}_J\right)$ (right y-axis) of the model during dorsal closure. The boxplots represent data across three intervals of $\Delta A (\Delta A<0.4,0.4\le \Delta A<0.7,\Delta A\ge 0.7)$. Whiskers extend to the 5th and 95th percentiles, while the boxes delineate the interquartile range, and the horizontal lines within the boxes indicate the median values. An ANOVA followed by a post-hoc Tukey’s HSD test was conducted to assess statistical significance (*: $p<0.1$, **: $p<0.05$). We performed and evaluated cuts of $N=97$ junctions. (E) Average initial recoil velocity of vertices after laser cutting as a function of junction straightness (ratio of the inter-vertex distance $\left(d_v\right)$ to the junction length $\left(L\right)$, see inset) immediately before cutting. Junction recoil velocity is independent of junction straightness (fitted with the red dashed line) until $S=d_v/L\gtrsim 0.93$. The crossover point at $d_v/L\approx 0.93$ marks the intersection of the red and blue dashed lines; the latter fits the data points in the gray-shaded region, indicating that the recoil velocity increases strongly and approximately linearly with junction straightness in this regime. (F) Comparison of experimental junction straightness (left y-axis) and model cellular junction tension (right y-axis) during dorsal closure.

Fig. 3.. Phase diagram of amnioserosa during dorsal closure.

Phase diagram in $\overline{q}$ vs. $\Delta A$ space, with solid (open black circles) and fluid (open gray squares) states as shown. The solid-fluid transition is marked by the black dashed line. Both experiment (green dashed line) and model (blue dashed line) are within the solid phase throughout dorsal closure.