Correct regionalization of a tissue primordium is essential for coordinated morphogenesis
Abstract
During organ development, tubular organs often form from flat epithelial primordia. In the placodes of the forming tubes of the salivary glands in the Drosophila embryo, we previously identified spatially defined cell behaviors of cell wedging, tilting, and cell intercalation that are key to the initial stages of tube formation. Here, we address what the requirements are that ensure the continuous formation of a narrow symmetrical tube from an initially asymmetrical primordium whilst overall tissue geometry is constantly changing. We are using live-imaging and quantitative methods to compare wild-type placodes and mutants that either show disrupted cell behaviors or an initial symmetrical placode organization, with both resulting in severe impairment of the invagination. We find that early transcriptional patterning of key morphogenetic transcription factors drives the selective activation of downstream morphogenetic modules, such as GPCR signaling that activates apical-medial actomyosin activity to drive cell wedging at the future asymmetrically placed invagination point. Over time, transcription of key factors expands across the rest of the placode and cells switch their behavior from predominantly intercalating to predominantly apically constricting as their position approaches the invagination pit. Misplacement or enlargement of the initial invagination pit leads to early problems in cell behaviors that eventually result in a defective organ shape. Our work illustrates that the dynamic patterning of the expression of transcription factors and downstream morphogenetic effectors ensures positionally fixed areas of cell behavior with regards to the invagination point. This patterning in combination with the asymmetric geometrical setup ensures functional organ formation.
Editor's evaluation
This paper addresses a fundamental question in developmental biology, that is, how morphogenetic movements driving tissue folding are patterned to occur with the correct spatiotemporal dynamics. By correlating dynamic patterns of transcription factor expression with rigorous, quantitative analyses of cell behaviors across the salivary gland primordium in Drosophila, their results suggest Hkb and Fkh transcription factor patterning induces switches in cell behaviors at fixed positions to promote continued morphogenesis of the tubular structure. This mechanism is likely to be more generally important for the development of complex tubular organs.
https://doi.org/10.7554/eLife.72369.sa0Introduction
Complex three-dimensional organs arise from simple tissue primordia, and in many cases, these primordia are flat polarised epithelial sheets. Early in development the expression of the first patterning genes broadly sets up embryonic regions. This is followed by the activation of gene regulatory networks that specify the fate of tissue primordia in defined locations (Gilmour et al., 2017; Sidor and Röper, 2016). The patterning and fate determination gene products do not directly affect morphogenetic changes, but rather instruct the expression of downstream morphogenetic effectors that drive a tissue primordium down a path of defined physical changes. In many well-studied cases of tissue morphogenesis, such as mesoderm invagination or germband extension in the fly, many morphogenetic effectors are induced evenly across the tissue primordium. Regional differences in the physical changes or ‘behaviors’ of cells observed in these primordia can arise due to physical effects and feedback or interference either from cells within the primordium or from surrounding tissues (Chanet et al., 2017; Collinet et al., 2015; Lye et al., 2015). A further layer of control and complexity is added as the primordia themselves can also be patterned to guide the differential behavior and changes of groups of cells within a primordium. Such pre-patterning of morphogenetic events within a single tissue primordium is much less understood.
We use the formation of a narrow-lumen tube in the Drosophila embryo as a model system to identify key requirements for successful organ formation. The symmetrical tubes of the embryonic salivary glands form from a flat, nearly circular epithelial primordium, the salivary gland placode (Figure 1A). We recently uncovered a patterning of key cell behaviors that drive the initiation and earliest stages of the tube budding process (Gillard et al., 2021; Sanchez-Corrales et al., 2018). Interestingly, the point of invagination from the circular placode primordium is not in the center of the placode, but in an asymmetric, eccentric position in the dorsal-posterior corner (Figure 1A). Using quantitative morphometric methods, we showed that the early apical constriction and associated cell wedging at the position of the invagination pit initiates the tissue bending (Booth et al., 2014; Sanchez-Corrales et al., 2018). Isotropic constriction driven by apical-medial actomyosin is concentrated near the forming invagination pit (Figure 1A’ and A’’), but at a distance to the pit a second cell behavior dominates cell intercalation. Driven by a polarised junctional accumulation of actomyosin, cell intercalation events such as T1 exchanges and rosette formation and resolution help to elongate the tissue radially toward the invagination pit and contract it circumferentially, via circumferential neighbor gains (Figure 1A’; Sanchez-Corrales et al., 2018).
This asymmetric setup of the placode raises two questions: first, is the asymmetry required for successful organ formation or could a symmetrical setup not lead to the same result? An indication that the asymmetry of the salivary gland placode is important for wild-type tube formation is given by the reports that in certain mutant situations where a more central invagination appears to form, the overall morphogenesis is disrupted and the invaginated structure that forms is sack-like and irregular. This is true for embryos lacking the transcription factor Hkb that is expressed early in the placode (Myat and Andrew, 2000b; Myat and Andrew, 2002), and could therefore potentially implicate Hkb in establishing early placode patterning.
Second, how is this asymmetry of the placode established prior to morphogenesis commencing and then maintained throughout the process? Previous studies have shown that several proteins known to be expressed in the placode can be detected at the mRNA level initially in the dorsal posterior region of the placode. These include the transcription factor Forkhead (Fkh) (Myat and Andrew, 2000a), the dynein-associated protein Klar (Myat and Andrew, 2002), as well as the kinase Btk29 (Chandrasekaran and Beckendorf, 2005). Defined expression of these and other factors at the dorsal-posterior corner could set up the initiation of cell shape changes in this region. In agreement with this, we previously found that fkh−/− mutants specifically fail to apically constrict in the dorsal-posterior corner, though they do show intercalation behavior, and these mutants do not invaginate a tube (Sanchez-Corrales et al., 2018). Following on from the early regionalized behaviors that drive the initial tissue bending and tube invagination (Sanchez-Corrales et al., 2018), it is thus far unclear what mechanisms continue to drive the sustained invagination until all secretory cells have internalized. This could be due to continued transcriptional patterning and activation of downstream effectors, and/or could be due to mechanical signaling and feedback.
Here, we investigate how the asymmetry is set up and maintained during salivary gland tube formation, in particular beyond the initial stages. We found that apical constriction is always strongest near the invagination pit at any moment in time. As more and more coronae of cells approach the region near the pit, they switch behavior from predominantly intercalating to predominantly apically constricting/wedging. This continued switch to apical constriction leads to a smooth continued invagination process akin to a ‘standing wave’ of apical constriction, through which the placodal cells flow. The continued near-pit constriction is driven by apical-medial myosin, whose peak intensity and peak pulsatile strength is tracked in a fixed position close to the pit over time. Using live-imaging of endogenous fluorescently tagged versions of Hkb and Fkh, we show that their expression levels pre-pattern where initial constriction is occurring. Both are upstream of a dynamic patterning of the GPCR ligand Fog, that in turn activates apical-medial myosin. Loss of the correct pre-patterning, such as in hkb−/− mutants or when Fog is overexpressed, leads to loss of the symmetrical final organ shape with a narrow lumen and instead gives rise to expanded and sack-like glands with widened lumena.
Our work establishes that after the tissue is patterned at the posterior corner, cells follow dynamic rules of behavior as they approach the pit. Symmetrical tube morphogenesis thus relies on a regionalization and patterning of the primordium long before cells show distinctive behaviors. Establishing general mechanisms that modulate the orderly and repetitive behavior of cells to form correctly shaped simple model organs such as the salivary glands will assist the understanding of the morphogenesis of other more complex tubular organs in mammals.
Results
Patterned apical constriction remains fixed around the pit over time
The initial morphogenetic change that occurs during the budding of the tube of the salivary glands in the Drosophila embryos is the apical constriction of cells at the point of the future site of invagination, the pit (Figure 1A and A’; Girdler and Röper, 2014; Myat and Andrew, 2000a; Myat and Andrew, 2000b; Sanchez-Corrales et al., 2018). We recently provided the first quantitative morphometric analysis of the very early stages of this process, and revealed that, first, isotropic apical constriction near the future pit commences long before changes at the tissue level are apparent, and that, second, the apical constriction is indicative of a 3D cell behavior of wedging of cells at the pit (Figure 1A’). Importantly, we uncovered a strong regionalization of cell behaviors during the early stages of tube budding: cells near the invagination point predominantly showed apical constriction behavior while cells at a distance from the pit predominantly showed intercalation behavior (Sanchez-Corrales et al., 2018; Figure 1A). What is unclear is how the early tube initiation via isotropic apical constriction of cells, combined with directional intercalation of the cells still positioned away from the forming pit, evolves over time. The process of tube formation is characterized by the continuous invagination of cells and the formation of a symmetrical narrow-lumen tube from an asymmetric placode primordium. We can imagine two scenarios: in one, cell behavior is fixed at early stages of budding, that is, apical cell constriction is restricted to the pit and only apparent during early stages. In a second scenario, cell behavior depends on a cell’s position relative to the invagination pit and thus cell behaviors change as cells move toward the invagination point.
To distinguish between these scenarios, we segmented and tracked placodal cells from time-lapse movies covering a –40 min to +70 min time interval of the invagination process, with t=0 min defined as the first occurrence of tissue bending (see Figure 1—figure supplement 1 for time intervals of analyzed movies and cell numbers analyzed) and performed quantitative analyses (Figure 1B and C). At the latest time point covered, the majority of secretory cells had completely internalized and formed the internal tube. As calculated previously (Sanchez-Corrales et al., 2018), we employ a radial coordinate system centered at the invagination point at t=0 min. We specify the ‘near to the pit’ region as the region located between the invagination point and up to 33% of the distance between the invagination point and the far edge of the placode (with this region usually being within 13.7±3 µm of the pit), and with the rest of the cells defined as ‘far from the pit.’
Cells positioned near the pit at t=0 min (Figure 1B, pink) constricted their apices isotropically and rapidly and were internalized by about t=+ 20 min (Figure 1B and C and dashed curve in Figure 1E). Cells positioned far from the pit at t=0 min (Figure 1B, gray) and moving into a position closer to the invagination point began constricting (Figure 1B and D and solid curve in Figure 1E; Figure 1—video 1), suggesting that cells changed their behavior depending on where they were positioned in the placode with respect to the invagination point. Cells near the pit displayed isotropic constriction and thus a negative cumulative apical area change from the start (Figure 1E, pink highlighted), whereas cells initially positioned far from the pit did not change their apical area over the first ~20 min (Figure 1C, gray highlighted), but then began to display isotropic apical constriction (Figure 1E, blue highlighted). The cell behaviors and changes of cell behaviors were very apparent when we analyzed individual cells: cells located in the near pit region (Figure 1F, yellow and green cells) constricted progressively, whereas cells located in the region far from the pit at t=0 min (Figure 1F, blue and magenta cells) retained their apical area initially and only began constricting from t>+20 min onwards. As described previously (Sanchez-Corrales et al., 2018), cells far from the pit initially showed predominantly a cell intercalation behavior that constricts the tissue primordium circumferentially whilst expanding it radially, thereby moving cells closer to the invagination pit (Figure 1G). As they approached the invagination pit, these cells also started to apically constrict (Figure 1E solid line, highlighted in blue and Figure 1F).
These data strongly suggest that an initially pre-patterned behavior, isotropic constriction at the pit and intercalation further away, is dynamically adjusted during tube budding from a flat epithelial sheet, so that cells within a similar distance to the invagination pit display the same behavior.
Sustained apical constriction near the pit is driven by sustained apical-medial myosin
Early apical constriction at the pit is driven by a highly dynamic apical-medial actomyosin network (Booth et al., 2014; Sanchez-Corrales et al., 2018). We next sought to determine whether the apical-medial actomyosin dynamics would differ among regions in the placode. We analyzed apical-medial actomyosin levels in cells across the placode, and firstly illustrate our findings in two snapshots: cells near the future invagination point at t=0 min showed strong accumulation of apical-medial myosin (Figure 2A, A’’ and C; see also Figure 2—video 1), whereas cells far from the pit and close to the anterior placode boundary (Figure 2A and A’) showed mainly junctional myosin II accumulation. At t=30 min, the cells now located right next to the invagination pit continued to display strong apical-medial myosin accumulation (Figure 2B and B’’), but in addition cells that were initially ‘far from the pit,’ but were now in a closer proximity to the pit though still near the anterior boundary, now also started to display apical-medial actomyosin (Figure 2B and B’). Thus, being moved into a position near the invagination pit appeared to determine levels of apical-medial myosin accumulation. Looking placode-wide over the whole time period analyzed, the peak of the apical-medial myosin fluorescence intensity tracked in a near-fixed position to the edge of the pit over time (Figure 2D, white dotted line). The location of the highest medial myosin intensity was also where cells showed the shortest cycle lengths of medial myosin II fluctuations (Figure 2E, white dotted line). Interestingly, cells near the pit at t=0 min showed an overall shorter myosin II cycle length (the time period elapsed between two peaks of highest myosin II intensity during the oscillation, see schematic in Figure 2E) compared to cells located far from the pit (Figure 2F). As longer cycle length correlates with unproductive myosin cycles (Booth et al., 2014), this was in line with cells further away from the pit showing less apical constriction. Furthermore, the strength of myosin oscillations can be expressed as the product of the amplitude and frequency of the oscillation (Booth et al., 2014). This myosin strength across the placode at t=0 min was much higher in cells near the invagination point when plotted against the distance to the pit, with a gradual decrease in myosin strength toward the placode boundary (Figure 2G). However, cells originally far from the pit at t=0 min increased their myosin strength over time (and concomitantly decreased their cycle length; see Figure 2—figure supplement 1) whilst moving into a position closer to the pit, so that by t=55 min their myosin strength had increased to levels previously shown by cells near the pit at t=0 min (Figure 2H).
Thus, our data show that the early asymmetric setup of high apical-medial myosin II intensity and activity that was clustered near the pit was maintained during the continued invagination of the tube, in turn driving the continued apical constriction near the pit.
hkb-/- mutants show a delayed symmetrical apical constriction
In order to uncover the role of the initial asymmetric patterning of cell behaviors within the placode, we turned our attention to mutant situations: we already knew that the transcription factor Fkh was important for invagination, as fkh−/− mutants do not show localized apical constriction or any invagination (Myat and Andrew, 2000a; Sanchez-Corrales et al., 2018). A possibly more interesting situation is present in a previously published mutant of the transcription factor Huckebein (Hkb) that was reported to show a central invagination pit, combined with malformed invaginated salivary glands at later stages (Myat and Andrew, 2000b; Myat and Andrew, 2002).
We collected and segmented time-lapse movies of hkb−/− mutant embryos over the same time period as the control wild-type embryos (Figure 3 and Figure 3—figure supplement 1; Figure 3—videos 1–3). Whereas in wild-type embryos, at t ~ +40 min, most secretory cells of the salivary gland placode have already invaginated (Figure 3A’’’, B’’’ and C’’’), in the hkb−/− mutants no invagination had occurred yet (Figure 3D’’’, E’’’ and F’’’). Eventually, these embryos invaginated a dilated tube through an enlarged pit (see Figure 3—figure supplement 2). We analyzed the spatial distribution of apical area changes over time in wild-type and hkb−/− placodes (Figure 3G and H). In the wild-type, the area of strongest apical area constriction was always located in a position near the pit, similar to the distribution of medial myosin intensity and activity (Figure 2D and E). In contrast, in hkb−/− mutant placodes cells began to constrict in an aberrant central position from about t=+20–30 min (Figure 3G and H). Thus, hkb−/− mutants show a delayed central invagination pit compared to the eccentric pit observed in wild-type placodes (Figure 3(I)1’). For our further quantitative analyses, we now took this change to a central pit position in the mutants into account. We plotted cumulative apical area change for cells in hkb−/− mutant placodes, split into cells near the (now central) pit and cells far from the pit (at the placode periphery). This analysis showed that cells in hkb−/− mutant placodes located near the central pit started to constrict at a much slower rate than cells in wild-type placodes, and only from about t=+20 min on (Figure 3J, wild-type data reproduced from Figure 1E). In fact, the rate of the constriction of the central cells in hkb−/− mutants was very similar to that of the cells located at a distance to the pit in the wild-type (compare blue dashed and magenta solid curves), suggesting that the central cells in hkb−/− mutants continue to behave identical to the wild-type cells occupying the same position. Cells far from the central pit in hkb−/− mutants showed a slight apical area expansion up to t=+20 min followed by a slight constriction, but did not invaginate over the time interval shown in the plots. In agreement with the delayed central constriction of cells in hkb−/− mutant placodes, we found no accumulation of apical-medial actomyosin at the eccentric position where a pit forms in the wild-type (compare Figure 3E and K to Figure 2D). Instead, only at t~+40 min did a central group of cells in hkb−/− mutant placodes display increased accumulation of apical-medial myosin II (Figure 3K), inducing a shallow tissue bending in this region at this point (see also Figure 3—figure supplement 2).
Thus, in comparison to wild-type placode with an eccentric invagination, in hkb−/− mutant placodes a central invagination formed with a significant delay in time.
hkb−/− mutants display disrupted cell behaviors across the placode
The above analyses strongly suggest that the initial asymmetrical setup of cell behaviors in the placode might be disrupted in hkb−/− embryos. We therefore performed a strain (deformation) rate analysis of placodes in hkb−/− mutants in comparison to wild-type embryos to assess this quantitatively. The strain rate analysis is based on the assumption that any change in tissue shape can be accounted for by changes in cell shape and cell intercalations (Blanchard et al., 2009; Sanchez-Corrales et al., 2018). In many tissues, cell division or cell death and delamination also contribute to strain rates, but these are absent in the salivary gland placode during invagination and can thus be neglected. The origin of the radial coordinate system used for this analysis was specified at t=0 min and positioned at the center of of the placode in the hkb−/− mutant embryos (where the central constriction eventually occurred) and in wild-type embryos at the dorsal-posterior corner of the placode (where the wild-type pit will form). In wild-type placodes, cells that were located near the eccentric invagination pit (Figure 4A; with the ‘near’ position defined at t=0 min; see also Figure 4—figure supplement 1) constricted isotropically, with near equal change in strain at the tissue and cell shape level, and very little intercalation strain (Figure 4A’ and A’’, solid curves). These cells in wild-type placodes could only be tracked up to t=20 min as at this point they had all invaginated into the embryo to form the initial part of the tube. In contrast, in the hkb−/− mutants, cells near the forming central pit (Figure 4A; defined to be in a ‘near’ position at t=0 min) showed a much reduced tissue and cell shape strain rate, with near identical radial and circumferential contributions, and constriction only commenced at about half the rate of the wild-type placodes from t=+20 min onwards (Figure 4A’ and A’’, dashed curves; see also Figure 4—figure supplement 1). In wild-type placodes, cells far from the eccentric pit (Figure 4B; defined to be ‘far’ at t=0 min), showed cell intercalations leading to circumferential convergence and radial extension of the tissue over the first approximately 20 min after invagination commenced, as we also reported previously (Figure 4B’ and B’’, solid curves; tissue and intercalation strain rate showing expansion radially and contraction circumferentially; Sanchez-Corrales et al., 2018). Interestingly, we then observed a clear switch in cell behavior of these cells to a second phase after t=+20 min. Once these cells were located closer to the pit position, they now also displayed apical contraction (beyond t=+20 min; Figure 4B’ and B’’, solid green curves). In the hkb−/− mutants, for cells at the placode periphery and hence far from the central constriction (Figure 4B), the tissue expanded slightly in both radial and circumferential direction (Figure 4B’ and B’’, dashed gray curves), with cell shapes narrowing radially (dashed green curve in Figure 4B’) and extending circumferentially (dashed green curve in Figure 4B’’). We could observe a similar rate of intercalation in these cells as in cells far from the pit in wild-type placodes, with some radial elongation and circumferential contraction (dashed orange curves in Figure 4B’ and B’’). As the tissue strain rate in the hkb−/− mutants showed no overall signature of directional changes with only a slight expansion both radially and circumferentially, we were curious to understand better what these tissue-scale signatures could represent at the individual cell level.
We therefore analyzed individual events of neighbor gains across the tissue in wild-type and hkb−/− mutant embryos. We compared two types of intercalation events, those leading to circumferential neighbor gains that in wild-type placodes explain the observed intercalation strain rates, and those leading to radial neighbor gains, and thus potentially opposing the changes seen at the tissue and intercalation strain rate level (Figure 4C). In wild-type placodes, circumferential neighbor gains outweigh radial ones by a large margin throughout (Figure 4D, dashed curves), leading to a steady increase in cumulative productive neighbor gains for cells across the placode (Figure 4E, yellow curves). In contrast, in hkb−/− mutant embryos circumferential and radial neighbor gains across the placode were overall much reduced in number and occurred with equal frequency (Figure 4D, solid curves). Furthermore, cells near the delayed central pit in hkb−/− mutants showed no productive neighbor gains (defined as the difference between circumferential and radial neighbor gains, corrected for the proportion of interfaces that could successfully be tracked; Figure 4E, purple dashed curve), whereas cells at the periphery away from the central pit showed productive neighbor gains (Figure 4E, solid purple curve). Active neighbor gains in the placode and other morphogenetic events are initiated and driven by junctional actomyosin (Sanchez-Corrales et al., 2018; Tetley et al., 2016), prompting us to analyze junctional myosin polarisation in hkb−/− mutant embryos. We previously reported on the unipolar and bipolar enrichment of junctional myosin in the wild-type at early stages of placode morphogenesis (Sanchez-Corrales et al., 2018) and expand it here to the whole process of secretory cell invagination in comparison to hkb−/− mutants (Figure 4F and G). Whereas in control placodes junctional myosin is clearly enriched in circumferential junctions by both measures, up to the point of most secretory cells having internalized by t=+40 min (Figure 4F and G; dashed green curve vs. dashed purple curve), in hkb−/− mutant placodes myosin unipolarity and bipolarity are strongly reduced, even though some bipolarity remained (Figure 4F and G, solid curves). Thus, the major driver of active intercalations appeared to be strongly diminished in the mutants, posing the question of what could drive the observed intercalations here. We also analyzed cell elongation in cells far from the pit, with cells in wild-type placodes being initially elongated circumferentially at t<+20 min, while early intercalation events are initiated by circumferential junctional myosin (Figure 4H, gray curve; Sanchez-Corrales et al., 2018). These cells then became increasingly elongated radially as they moved into a position nearer to the pit and were being actively funneled toward the pit. In hkb−/− mutant embryos, cells far from the pit at the periphery of the placode very quickly became circumferentially elongated and remained so throughout the time period analyzed. This aberrant circumferential elongation together with the mild bipolarity of junctional myosin polarisation that we observed prompted us to analyze individual events of cell intercalation in time-lapse movies (Figure 4I–L). In cells far from the central pit in hkb−/− mutants, we could observe groups of cells clearly undergoing intercalations (Figure 4I and J) with individual cells making circumferential neighbor gains and elongating (cells 1 and 2 in Figure 4K), whilst the overall shape of the local group of cells undergoing the intercalation event was not changing (Figure 4L). Thus, although cells far from the aberrant central pit in hkb−/− mutants still intercalated and gained circumferential neighbors, this behavior lacked several key signatures of wild-type intercalations, and overall in the hkb−/− mutants intercalations did not lead to convergence and extension of the tissue toward the misplaced pit.
These analyses show that in the absence of Hkb, the initial radial patterning of cell behaviors across the placode is highly aberrant. Cells located at the position corresponding to where the eccentric pit would be located in the wild-type clearly failed to constrict apically, suggesting that Hkb could be involved in the specification of this pit behavior observed in wild-type.
The salivary gland placodal primordium is asymmetrically patterned by Hkb and Fkh prior to morphogenesis
The loss of correct patterning of cell behaviors in the early placode in hkb−/− mutant embryos suggested that Hkb was involved in establishing these in wild-type embryos. Furthermore, hkb mRNA expression has been reported to show a dynamic pattern by in situ hybridization and β-gal reporter labeling (Myat and Andrew, 2000b). In order to be able to analyze Hkb protein levels dynamically and in comparison to cell shapes, we generated a Venus-tagged version of Hkb using CRISPR/Cas9 to tag the endogenous protein. Using this in vivo reporter, we detected Venus-Hkb at high levels concentrated in the dorsal-posterior region of the just forming placode, in the area where the invagination pit will form, already at t=–30 min prior to the first tissue bending (Figure 5A, A’, D and M and Figure 5—figure supplement 1B,B'' ), and even as early as t=–63 min prior to invagination (Figure 5—figure supplement 1A). The expression then expanded more broadly and anteriorly across the whole placode over time (Figure 5B–C’, E, F and M and Figure 5—video 1). At the moment of first tissue bending (t=0 min) Venus-Hkb was found across the placode, with elevated levels still at the invagination pit position, and a second increase toward the anterior. At late stages (t>>0 min) Venus-Hkb levels close to the invagination point were lower than further anterior (Figure 5M).
As mentioned above, embryos mutant in fkh also show a lack of early apical constriction at the dorsal-posterior corner of the placode (Myat and Andrew, 2000a; Sanchez-Corrales et al., 2018). Similar to Hkb, analysis of fkh mRNA and protein levels in fixed embryos indicated an early enrichment near the forming invagination pit (Myat and Andrew, 2000a; Sanchez-Corrales et al., 2018). Using an in vivo genomic GFP-Fkh reporter (Spokony and White, 2013), we followed Fkh protein levels and patterning over the same time period whilst also assessing apical cell shapes (Figure 5G–L). Already at t=–60 min prior to tissue bending, Fkh-GFP was visible in the cells near the future invagination pit (Figure 5G, G’, J and N and Figure 5—figure supplement 1C–C''), from where it spread anteriorly and ventrally across the whole placode (Figure 5—video 2). At t=0 min, Fkh-GFP was already detected in all cells of the placode, though with levels still highest near the pit and decreasing toward the anterior edge of the placode (Figure 5H, H’, K and N). At t=+30 min and beyond, levels of Fkh-GFP were more uniform across the placode (Figure 5I, I', L, N).
Thus, both Hkb and Fkh transcription factors showed a clear temporally and spatially graded pattern of expression within the early placode that developed dynamically. Changes were especially pronounced at the dorsal-posterior corner, with highest levels here already an hour before tissue bending commenced. Combined Hkb and Fkh action could thereby define the future position of the invagination pit. We therefore compared Venus-Hkb as well as Fkh-GFP fluorescence intensity levels from time-lapse movies that in parallel allowed us to determine the corresponding apical area of cells (Figure 5O–P’). Venus-Hkb intensity levels peaked in cells located at the site of the forming eccentric invagination pit at about t ~ –30 min (Figure 5O), with an appearance of smaller apical areas in this position following on about 30–40 min later (Figure 5O’). Similarly, Fkh-GFP intensity began to peak at about the same time, t ~ –30 min, and then tracked at high level at a fixed distance to the pit over time (Figure 5P). Again also here, the increased reduction in apical area appeared with a delay of about 30–40 min (Figure 5P’). Such a delay would be consistent with the estimated time required for the transcription of the fog locus, followed by translation, protein folding, secretion, binding to the receptor, and signaling leading to myosin phosphorylation (Garcia et al., 2013; Shamir et al., 2016).
Thus, Hkb and Fkh dynamic expression at the future eccentric pit could prefigure the start of apical constriction at this site. The continuous high levels of Fkh tracking at a fixed distance to the pit suggest that Fkh is key to the continuous apical constriction of cells moving into a near pit position. These dynamic changes might trigger the graded and patterned expression of targets of both transcription factors, leading to the observed patterned cell behaviors we uncovered, as well as the switch between intercalation and constriction.
Asymmetric Fog expression is controlled by Hkb and Fkh and is upstream of early differential behaviors in the placode
With the transcription factors Fkh and Hkb both displaying intriguing protein expression patterns across the early placode pre-morphogenesis, targets of both factors are likely to play key roles in instructing cell behaviors such as wedging and cell intercalation. Interestingly, an upstream activator of apical-medial actomyosin activity, the GPCR-ligand Folded gastrulation (Fog) was found to be dependent on Fkh (Chung et al., 2017; Dawes-Hoang et al., 2005; Figure 6) and suggested to be downstream of Hkb (Myat and Andrew, 2000b; though in this publication as ‘data not shown’). Fog is an apically secreted ligand that acts in an autocrine fashion, and thus the cells that switch on Fog expression will apically constrict (Dawes-Hoang et al., 2005). Previously published in situ hybridization of fog mRNA suggested increased levels near the forming pit, similar to the early Hkb and Fkh enrichment (Nikolaidou and Barrett, 2004). Fog expression and signaling in particular could not only be upstream of the isotropic constriction of cells near the forming pit, but could also be involved in maintaining this behavior in new coronas of cells moving into a position near the pit over time.
We therefore decided to carefully analyze Fog protein levels in comparison to apical area and apical-medial myosin intensity in salivary gland placodes from early to late invagination stages (Figure 6A–I). Fog levels were low at mid-stage 10 concomitant with placode specification (Figure 6A, A’ and C). Just before tissue bending commenced at t=0 min, Fog was already strongly enriched in cells near the forming invagination pit, and these cells were the ones already showing the decreased apical area (Figure 6D–F). Fog protein levels remained enriched near the invagination point once the tube had started to internalize (Figure 6G–I). In fact, Fog protein levels in cells were positively correlated with both smaller apical areas (Figure 6K) as well as higher average apical-medial myosin II intensity (Figure 6L; Figure 6—figure supplement 1).
hkb−/− mutant embryos, as detailed above, showed a very clear disruption of the early asymmetric placodal patterning. At late stage 11, Fog protein levels in hkb−/− mutant embryos appeared to be similar to wild-type embryos and enriched in the placode compared to the surrounding epidermis, though the pattern of Fog enrichment was altered. Cells at the posterior corner, where the invagination pit would have formed in wild-type placodes, lacked high Fog protein levels, whereas only cells in the center of the placode now showed high levels of Fog protein (Figure 6M–O). When quantitatively analyzing many hkb−/− mutant embryos in comparison to wild-type, the peak intensity of Fog protein was clearly shifted to a central, much more symmetrical, position (Figure 6P). In contrast, although we could detect some Fog protein in fkh−/− embryos at late stage 11, it was no longer enriched at the position where the invagination pit forms in the wild-type. The levels also did not appear raised anywhere within the placode compared to the surrounding epidermis (Figure 6Q-S and Figure 6—figure supplement 1 D-D’’), as the quantification of many placodes confirmed (Figure 6T). In both mutants, altered levels and patterning of Fog protein still correlated with the now altered levels and patterns of residual apical-medial myosin (Figure 6M and Q).
In order to understand the cause of the central and delayed constriction in hkb−/− mutant embryos upstream of the also only centrally localized Fog, we examined Fkh protein localization in hkb−/− mutants (Figure 6—figure supplement 1E-J). At late stage 10, when the Fkh protein expression pattern is just expanding from the eccentric pit position to the central part of the placode in the control, the pattern was indistinguishable in the hkb−/− mutants (Figure 6—figure supplement 1 E-F’ and H-I’). The pattern remained identical to wild-type placodes throughout stage 11 and beyond (Figure 6—figure supplement 1 G,G’, J, J’). Hence the natural expansion of Fkh protein expression to the cells in the center of the placode that is also observed in the hkb−/− mutant placodes might be the cause of the delayed constriction of the cells in this region, whereas the earlier initial constriction at the site of the eccentric pit was absent in the mutant because it requires both Fkh and Hkb. Both transcription factors appeared in fact to be independently regulated, as also Hkb levels and protein expression pattern were unchanged in fkh−/− mutant embryos (Figure 6—figure supplement 1 K-N’).
Thus, the asymmetrical setup of the salivary gland placode prior to initiation of tube budding relies on a correct spatial and temporal patterning of the key transcription factors, Hkb and Fkh. These transcription factors are in turn upstream of GPCR/Fog signaling that leads to patterned and sustained apical-medial myosin activity and hence contractile cell behaviors.
Changes to invagination pit size or position impact invagination and perturb final organ shape
Our data and previously published results show that the asymmetrical invagination point depends on patterning by transcription factors and downstream GPCRs/Fog signaling at the posterior corner, leading to a focussed and point-like invagination pit (Chung et al., 2017; Myat and Andrew, 2000a; Myat and Andrew, 2000b; Myat and Andrew, 2002; Sanchez-Corrales et al., 2018). We next sought to test the importance of this asymmetric setup for successful morphogenesis of a narrow-lumen and symmetrical tube. To do so, we overexpressed Fog in the salivary gland placode only, using UAS-Fog x fkhGal4 which led to placode-wide increased levels of Fog from early stages on (Figure 7 and Figure 7—figure supplement 1). Compared to the focused apical constriction near the forming invagination pit in wild-type placodes at early stage 11 (Figure 7A and A’), when Fog was overexpressed across the placode most secretory cells of the placode constricted simultaneously (Figure 7B, B’, C and C’). At late stage 11, when in wild-type placodes cells had begun to invaginate and a short symmetrical narrow-lumen tube had formed (Figure 7D and D’), in UAS-Fog x fkhGal4 embryos the apically constricting cells had formed a large shallow depression across the whole dorsal part of the placode, with the beginning of a deeper invagination at the dorsal-posterior corner (Figure 7F and F’). Following on from this, at stage 12, in contrast to the narrow opening to the invagination in wild-type placodes (Figure 7E), in UAS-Fog x fkhGal4 placodes this whole section usually formed a deepened very wide pit (Figure 7G). This widened tube in Fog-overexpressing embryos, at stages 15–16 when all cells had invaginated, displayed a significantly enlarged and flattened lumen compared to wild-type tubes (Figure 7H–J). Interestingly, overexpression of Fkh using UAS-Fkh x fkhGal4 led to a similarly enlarged area of apical constriction early on and invaginated tubes with aberrant and inflated lumens at late stages (Figure 7—figure supplement 2).
We set out to explain the link between an enlarged or widened invagination pit in the placode early during tube formation and the expanded lumen tube observed at late stages by considering why the invagination pit forms in its eccentric position in wild-type placodes. We formulated a simplistic theoretical consideration that assumes that following the specification of an initial group of cells that apically constrict and form the initial invagination (Figure 7K–M), neighboring coronae of cells will follow and invaginate (Figure 7K’–M’). Using an idealized version of apices of a wild-type placode with an eccentric pit, the number of cells invaginating, going from corona to corona, remains very steady and low (Figure 7K’ and green curve in Figure 7N). In contrast, a centrally located pit (Figure 7L) would lead to an ever-increasing number of cells per corona that would need to invaginate (Figure 7L’ and magenta curve in Figure 7N). An enlarged pit, a combination of the two previous scenarios and resembling aspects of the enlarged pit and invagination seen in UAS-Fog x fkhGal4 embryos (Figure 7M and M’), would lead to the invagination of a large number of cells initially, followed by similarly large numbers in each corona (Figure 7M’ and gray curve in Figure 7N). In fact, when counting invaginating coronae of cells in frames of time-lapse movies of wild-type (Figure 7O, green curve), hkb−/− (Figure 7O, magenta curve) curve and UAS-Fog x fkhGal4 (Figure 7O, gray curve) embryos, the actual numbers observed were in very good agreement with the theoretical considerations outlined above (Figure 7K–N).
In the hkb−/− mutant embryos, the delayed central area of constriction turned into a widened invagination and the enlarged pit at late stages often collapsed into an enlarged crescent of subduction at the boundary (Figure 3—figure supplement 2 and Figure 3—video 3; Figure 3—video 4), leading to malformed and enlarged lumen glands at the end of embryogenesis, as has been observed previously (Myat and Andrew, 2002).
Thus, the eccentrically placed invagination pit in the dorsal-posterior corner of the salivary gland placode is a key factor in ensuring the formation of a symmetrical narrow lumen tube during invagination, in part likely due to the control over the number of invaginating cells in each corona destined to internalize.
Discussion
In this work, we investigated how the sustained morphogenesis of a tissue is achieved, in this case, the continued invagination of epithelial cells forming a tubular organ. We previously identified cell behaviors that become established across the placode and drive the initial tissue bending (Sanchez-Corrales et al., 2018).
A number of mechanisms appear at play that in combination ensure that a symmetrical and narrow-lumen tube is formed from an initially flat epithelial primordium. A temporal and spatial component is provided by the dynamically changing pattern of the transcription factors Hkb and Fkh. These are expressed initially with the highest protein level at the point where the invagination pit will form, before expression expands further across the placode over time. This ensures that the cells in the dorsal-posterior corner are the first to experience high levels of activity of these transcription factors, and we show that one of the important targets, the GPCR-ligand Fog, is found at high levels in these cells in a Hkb- and Fkh-dependent manner with a slight temporal delay. A further spatial component is provided by the fact that, as cells continuously internalize to form the tube on the inside of the embryo, new coronas of cells are brought into a ‘near the pit’ position and adjust their behavior accordingly. This position-specific behavior of apical constriction is most likely triggered by the previous temporal expansion of Fkh and also Hkb expression to these cells, thereby priming them for their next step in the path to tube formation. Interestingly, though, the pattern of Fkh and Hkb expression from t=0 min onwards, though generally expanding across the placode, is not identical. This could suggest that downstream targets begin to diverge at later stages. The hkb−/− mutant analysis in particular suggests that the central cells in the placode that manage to constrict in the absence of Hkb might only require Fkh activity to initiate Fog expression and function. These cells therefore undergo their normal apical constriction once Fkh expression has reached the central position. As they now are the first cells to constrict and invaginate, the hkb−/− mutants show a central and delayed constriction and invagination.
Our work strongly suggests that the combined early patterning of Hkb and Fkh expression defines the early invagination pit and its eccentric position. How though is Hkb and Fkh expression initially limited to this point? Both are expressed in the placode downstream of the homeotic factor Sex Combs Reduced (Scr; Myat and Andrew, 2000b; Panzer et al., 1992). They are also not the only factors expressed near exclusively at this corner first. At the transcription factor level in particular one of the parasegmentally repeated stripes of Wingless (Wg) expression overlays the posterior part of the placode, its expression domain defined by previous parasegmental patterning (Clark, 2017). And at the time point when Hkb and Fkh expression at the future pit location commences, the Wg expression within the placode region is also confined to this corner only, rather than extending all the way to the ventral midline (data not shown). Future studies will show whether a combination of Scr and Wg transcription factor activity could activate early pit-defining expression of factors such as Fkh and Hkb. In addition to transcription factors, some potential downstream effectors also appear to be expressed near the future pit location first, these include Btk29/Tec29 (Chandrasekaran and Beckendorf, 2005), Crumbs (Myat and Andrew, 2002; Röper, 2012), 18-Wheeler (Kolesnikov and Beckendorf, 2007), and Klarsicht (Myat and Andrew, 2002). Cell-based genomics approaches such as single-cell RNAseq or ATAC-seq should in the future allow to precisely determine the genetic program that underlies the regionalization of the placode.
Activation of actomyosin contractility through GPCR signaling driven by the autocrine secretion of the ligand Fog is a module repeatedly used during tissue bending and invagination events in Drosophila embryogenesis. What differs are the upstream activating transcription factors: during mesoderm invagination Fog expression on the ventral side is controlled by Twist and Snail (Morize et al., 1998), whereas its expression during posterior midgut invagination is controlled by Tailless and Hkb (Weigel et al., 1990). The effect of Hkb on Fog expression is likely direct, as the Fog locus contains predicted Hkb-binding sites (using oPOSSUM, http://opossum.cisreg.ca/oPOSSUM3/). The identity of the GPCR that could mediate the Fog signal in the salivary gland placode is unclear thus far. Published and database (Flybase) in situs for the GPCRs previously reported to be involved in mesoderm invagination and germband extension, Mist and Smog (Kerridge et al., 2016; Manning et al., 2013), do not show an enriched expression in the placode. But with around 200 Drosophila GPCRs in existence, half of them orphan receptors (Hanlon and Andrew, 2015), the placode might well express yet another one, and we will focus on candidates previously identified in screens (Maybeck and Röper, 2009) and genomics approaches (our unpublished data) in the future.
The question of how cell behavior is modulated over time after the initial transcriptional patterning that defines primordium identity is key to understand general mechanisms that allow continuous organ formation. Different solutions seem to have been adopted in different tissues. Recent work analyzing mesoderm invagination in Drosophila has shown that the key morphogenetic effectors T48 and Fog show a gradient of expression due to differential timing of transcriptional activation, leading to different spatial accumulations of the transcripts (Lim et al., 2017) and proteins (Heer et al., 2017). Flattening the gradient experimentally (by expansion of the ventral domain) leads to a widened and perturbed invagination process. By contrast, the invagination of the posterior midgut or endoderm primordium in the fly embryo does not require sustained transcription, but instead a wave of apical constriction that is initiated by Fog expression in a subset of cells, appears to be propagated via mechanical feedback (Bailles et al., 2019). The salivary gland placode seems to follow a model of a temporal expansion or wave of transcriptional activity translating into positionally fixed behaviors. Nonetheless, our quantitative analysis of 3D cell behaviors strongly indicated that mechanical feedback does also play a role (Sanchez-Corrales et al., 2018). Which combination of mechanisms is ultimately at work to sustain a morphogenetic process could depend on a number of factors including the developmental stage, the speed of the process (consider ~15 min for mesoderm invagination versus ~30 min for endoderm invagination versus ~1.5 hr for secretory cell invagination in the salivary gland placode), as well as the influence of other nearby concomitant morphogenetic events.
The eccentric pit setup of the placode we suggest is key to wild-type tube formation with a narrow lumen, and we show that our simple geometrical considerations are in good agreement with the actual in vivo data of number of cells per corona invaginating. Whether the eccentricity only allows to control cell numbers of invaginating coronae, or whether such setup could also have further beneficial effects on for instance mechanical aspects of the tissue will be something we will actively explore both in vivo and in silico in the future. Curiously and maybe tellingly, the asymmetry of the placodal setup is retained in another tube formation process, the invagination of the posterior spiracles (Sidor et al., 2020; Simões et al., 2006). The posterior spiracle placode, like the salivary gland placode, is a flat epithelium with an eccentric invagination point. In this case, an initial tube is already present due to the previous invagination of the most posterior tracheal placode. Whether the asymmetrical setup is thus a conserved feature of further tube budding processes during development will be an exciting question to address. 3D mathematical modelling of tube morphogenesis assuming an eccentric versus centrally positioned invagination point will help to dissect the topological advantages of being asymmetric.
In summary, we have quantitatively unraveled a dynamic spatio-temporal patterning of transcription factors and switches in cell behaviors leading to positionally fixed behaviors during the morphogenesis of simple tubular organ. We foresee that such mechanisms will be important for establishing general mechanisms and morphogenetic modules at work during the morphogenesis of more complex tubular organs and potentially also for engineered tissues in a dish.
Materials and methods
Fly stocks and husbandry
Request a detailed protocolThe following transgenic fly lines were used: Armadillo-YFP (CPTI collection described in Lye et al., 2014), sqhAX3; sqh::sqhGFP42 (Royou et al., 2004) and fkhGal4 (Henderson and Andrew, 2000; Zhou et al., 2001) (kind gift of Debbie Andrew); UAS-Fog (kind gift of Thomas Lecuit); y1 w* cv1 sqhAX3; P[w+mC = sqh-GFP.RLC]C-42 M[w+mC = Ubi-TagRFP-T-CAAX]ZH-22A (Kyoto DGRC number 109822, referred to as sqhAX3;sqhGFP; UbiRFP); [sqhGFP42,UbiRFP; fkh[6]/TM3 Sb Twi-Gal4::UAS-GFP] (fkh[6] allele from Bloomington); [hkb2, p2, Scribble-GFP]; [sqhGFP42,UbiRFP; hkb2, p2/ TM3 Sb Twi-Gal4::UAS-GFP] hkb2 allele from Bloomington BDSC: 5457, [UbiRFP, Venus-HKB] (this study), FkhGFP (Bloomington BDSC:43951) and [FkhGFP; wgGal4, UAS-palmYFP] (generated from membrane Brainbow; Förster and Luschnig, 2012; Hampel et al., 2011). See Table 1 for details of genotypes used for individual figure panels.
The CRISPR Venus-HKB was created using the following gRNA targets flanking the N-terminus of Hkb:
CACCGCAAACCTACTCGCGACTT.
A linker sequence between the inserted fluorescent protein and the protein was also used:
GGAGGCGGAGGCTCGGGAGGCGGAGGCTCG.
Embryo immunofluorescence labeling, confocal, and time-lapse
Request a detailed protocolEmbryos were collected on apple juice-agar plates and processed for immunofluorescence using standard procedures. Briefly, embryos were dechorionated in 50% bleach, fixed in 4% formaldehyde, and stained with primary and secondary antibodies in PBT (PBS plus 0.5% bovine serum albumin and 0.3% Triton X-100). Anti-Hkb was a gift from Chris Doe (McDonald and Doe, 1997) anti-Fkh was a gift from Herbert Jäckle (Weigel et al., 1989) and anti-Fog was a gift from Naoyuki Fuse (Fuse et al., 2013), PY-20, (*P11120, Transduction Laboratories); anti-CrebA (DSHB). Secondary antibodies used were Alexa Fluor 488/Fluor 549/Fluor 649 coupled (Molecular Probes) and Cy3 coupled (Jackson ImmunoResearch Laboratories). Samples were embedded in Vectashield (Vector Labs).
Images of fixed samples were acquired on an Olympus FluoView 1200 or a Zeiss 780 Confocal Laser scanning system as z-stacks to cover the whole apical surface of cells in the placode. Z-stack projections were assembled in ImageJ or Imaris (Bitplane), 3D rendering was performed in Imaris.
For live time-lapse experiments (see Table 1), embryos were dechorionated in 50% bleach and extensively rinsed in water. Stage 10 embryos were manually aligned and attached to heptane-glue coated coverslips and mounted on custom-made metal slides; embryos were covered using halocarbon oil 27 (Sigma-Aldrich) and viability after imaging after 24 hr was controlled prior to further data analysis. Time-lapse sequences were imaged under a 40×/1.3NA oil objective on an inverted Zeiss 780 Laser scanning system, acquiring z-stacks every 0.58–3 min with a typical voxel xyz size of 0.22×0.22×1 μm3. Z-stack projections to generate movies in Supplementary Material were assembled in ImageJ or Imaris. The absence of fluorescent Twi-Gal4::UAS-GFP was used to identify homozygous hkb[2] and fkh[6] mutant embryos. The membrane channel images from time-lapse experiments were denoised using nd-safir software (Boulanger et al., 2010).
Cell segmentation and tracking
Request a detailed protocolCell segmentation and tracking were performed using custom software written in Blanchard et al., 2009. First, the curved surface of the embryonic epithelium was located by draping a ‘blanket’ over all image volumes over time, where the pixel-detailed blanket was caught by, and remained on top of binarised cortical fluorescence signal. Quasi-2D image layers were then extracted from image volumes at specified depths from the surface blanket. This step permitted to account for the curvature of the embryos in these quasi-2D projections. We took maximum intensity projections of a small number of near-surface image layers typically at 0–4 μm from membrane channels, for cell-apical analyses, apicomedial myosin and anti-Fog fluorescence channels and projections of image layers typically from 5 to 16 μm encompassing nuclear fluorescence for Venus-HKB and Fkh-GFP analyses. Membrane channels were filtered with median, top-hat, or high/low frequency filters as necessary to optimize subsequent cell tracking.
Cells in membrane channels were segmented using an adaptive watershedding algorithm as they were simultaneously linked in time. Manual correction of segmented cell outlines was performed for all fixed and time-lapse data. The segmentation of all the movies used in this study was manually corrected to ensure at least 90% tracking coverage of the placode at all times. Tracked cells were subjected to various quality filters (lineage length, area, aspect ratio, and relative velocity) so that incorrectly tracked cells were eliminated prior to further analysis. The number of embryos analyzed and the number of cells can be found in Figure 1—figure supplement 1, Figure 3—figure supplement 1.
Mobile radial coordinate system for the salivary gland placode
Request a detailed protocolWild-type movies were aligned in time using as t=0 min, the frame just before the first sign of invagination at the future tube pit was evident. hkb−/− mutants were aligned using as a reference of embryo development the level of invagination of the tracheal pits that are not affected in the hkb−/− mutant as well as other morphological markers such as appearance and depth of segmental grooves in the embryo. Cells belonging to the salivary placode (without the future duct cells that comprise the two most ventral rows of cells in the primordium) were then manually outlined at t=0 min using the surrounding myosin II cable as a guide and ramified forwards and backwards in time. Only cells of the salivary gland placode were included in subsequent analyses.
At t = 0 min, the center of the future tube pit was specified manually as the origin of a radial coordinate system, with radial distance (in µm) increasing away from the pit. Circumferential angle was set to zero toward posterior, proceeding anti-clockwise for the placode on the left-hand side of the embryo, and clockwise for the placode on the right so that data collected from different sides could be overlaid.
The radial coordinate system was ‘mobile,’ in the sense that its origin tracked the center of the pit, forwards and backwards in time, as the placode translated within the field of view due to embryo movement or to on-going morphogenesis.
As calculated previously (Sanchez-Corrales et al., 2018), we employ a radial coordinate system centered at the invagination point at t=0 min. We specify the ‘near to the pit’ region as the region located between the invagination point and up to 33% of the distance between the invagination point and the far edge of the placode (with this region usually being within 13.7±3 µm of the pit), and with the rest of the cells defined as ‘far from the pit.’
Morphogenetic strain rate analysis
Request a detailed protocolDetailed spatial patterns of the rates of deformation across the placode and over time quantify the net outcome of active stresses, viscoelastic material properties, and frictions both from within and outside the placode. We quantified strain (deformation) rates over small spatio-temporal domains composed of a focal cell and one corona of immediate neighbors over a ~3-min interval (Blanchard et al., 2009) and reviewed in Blanchard, 2017. On such 2D domains, strain rates are captured elliptically, as the strain rate in the orientation of greatest absolute strain rate, with a second strain rate perpendicular to this.
For the early morphogenesis of the salivary gland placode, in which there is no cell division or gain/loss of cells from the epithelium, three types of strain rates can be calculated. First, total tissue strain rates are calculated for all local domains using the relative movements of cell centroids, extracted from automated cell tracking. This captures the net effect of cell shape changes and cell rearrangements within the tissue, but these can also be separated out. Second, domain cell shape strain rates are calculated by approximating each cell with its best-fit ellipse and then finding the best mapping of a cell’s elliptical shape to its shape in the subsequent time point, and averaging over the cells of the domain. Third, intercalation strain rates that capture the continuous process of cells in a domain sliding past each other in a particular orientation, are calculated as the difference between the total tissue strain rates and the cell shape strain rates of cells. Strain rates were calculated using custom software written in L3Harris Geospatial IDL (code provided in Blanchard et al., 2009 or by email from G.B.B.).
The three types of elliptical strain rates were projected onto our radial coordinate system (Sanchez-Corrales et al., 2018) using the eccentric pit location in wild-type (dorsal-posterior corner) or the center of the misplaced invagination point in hkb−/− mutants. The radial coordinate system permitted us to analyze radial and circumferential contributions. Strain rates in units of proportional size change per minute can easily be averaged across space or accumulated over time. We present instantaneous strain rates over time for spatial subsets of cells in the placode (Figure 4—figure supplement 1), and cumulative strain ratios for the same regions over time. These plots were made from exported data as done previously in Sanchez-Corrales et al., 2018. To test for differences in instantaneous strain rates of wt and hkb−/− time-lapse movies we used a mixed-effects model implemented in R (lmer4 package as in Butler et al., 2009; Lye et al., 2015; Sanchez-Corrales et al., 2018) with a significance threshold of p<0.05. The phenotype (wt or hkb−/−) was considered a fixed effect while the variation between embryos from the same phenotype was considered a random effect.
Cell shape elongation (Figure 4H) was calculated by fitting an ellipse into each cell. The two axes of each ellipse were projected into radial coordinates and the ratio of them was calculated.
Neighbor exchange analysis
Request a detailed protocolWe used changes in neighbor connectivity in our tracked cell data to identify neighbor exchange events (T1 processes). Neighbor exchange events were defined by the identity of the pair of cells that lost connectivity in t and the pair that gained connectivity at t+1. The orientation of gain we defined as the orientation of the centroid-centroid line of the gaining pair at t+1. We further classified gains as either radially or circumferentially oriented, depending on which the gain axis was most closely aligned to locally. We did not distinguish between solitary T1s and T1s involved in rosette-like structures.
From visual inspection, we knew that some T1s were subsequently reversed, so we characterized not only the total number of gains in each orientation but also the net gain in the circumferential axis, by subtracting the number of radial gains. Furthermore, when comparing embryos and genotypes, we controlled for differences in numbers of tracked cells by expressing the net circumferential gain per time step as a proportion of half of the total number of tracked cell-cell interfaces in that time step. We accumulated numbers of gains, net gains, and proportional rate of gain over time for WT and hkb embryos (Figure 4D and E). Two sample Kolmogorov-Smirnov tests were used to determine significance at p<0.05 for data in Figure 3.
Automated medial myosin II quantification and oscillatory behavior and polarity
Request a detailed protocolWe extracted a quasi-2D layer image from the myosin II channel at a depth typically from 0 to 3 μm that maximized the capture of apical-medial myosin. We background-subtracted the myosin images and quantified the average intensity of apical-medial myosin as the fluorescence inside the segmented cell, excluding the fluorescence along cell-cell interfaces encompassing the cell edge pixels plus two pixels in a perpendicular direction on either side (considered junctional myosin).
The analysis of apical-medial fluorescence fluctuations was performed as described before in Blanchard et al., 2010 and Booth et al., 2014. Briefly, the apical-medial fluorescence values in the time series were detrended by subtracting a boxcar average with a window of 6 min (larger than the maximum expected fluctuation cycle length). All peaks and troughs were then identified in the fluorescence signal. Peaks or troughs that were associated with very small amplitude cycles were skipped. The remaining peaks and troughs were used to calculate the amplitude and the cycle length (peak to peak) and the amplitude (average peak to trough) of fluctuations.
The analysis of junctional myosin polarity was performed as described before in Sanchez-Corrales et al., 2018; Tetley et al., 2016. Briefly, the average junctional fluorescence intensity is considered a periodic signal around the cell and decomposed into individual components (modes) using Fourier analysis. The strength of myosin bipolarity corresponds to the amplitude of the second mode while the strength of unipolarity corresponds to the amplitude of the first mode. Similarly, the phase of modes 1 and 2 (uni-/bipolarity) represents the orientation of cell polarity. The uni- and bipolarity amplitudes are expressed as a proportion of the mean cell perimeter fluorescence. Cells at the border of the placode that are part of the supra-cellular actomyosin cable at the boundary were excluded from the analysis.
Hkb, Fkh, and Fog fluorescence intensity quantifications
Request a detailed protocolAverage fluorescence was calculated after background subtraction as the mean fluorescence inside a segmented cell excluding the cell edge pixels and two pixels in a perpendicular direction on either side. Fixed samples of Venus-Hkb and FkhGFP embryos of mid to late stage 10 (t<<0 min), early stage 11 (t+0 min) and late stage 11 (t>>0 min) were pooled together with relevant time points from live-imaging as shown in Table 1 and figure legends. Fluorescence was normalized in each image by dividing the average fluorescence of each cell by the 98th percentile value. Calculations were performed in R.
Data availability
All data generated and analysed in the manuscript are included in the manuscript and supporting files.
References
-
Taking the strain: quantifying the contributions of all cell behaviours to changes in epithelial shapePhilosophical Transactions of the Royal Society B 372:20150513.https://doi.org/10.1098/rstb.2015.0513
-
Patch-based nonlocal functional for denoising fluorescence microscopy image sequencesIEEE Transactions on Medical Imaging 29:442–454.https://doi.org/10.1109/TMI.2009.2033991
-
Actomyosin meshwork mechanosensing enables tissue shape to orient cell forceNature Communications 8:15014.https://doi.org/10.1038/ncomms15014
-
Local and tissue-scale forces drive oriented junction growth during tissue extensionNature Cell Biology 17:1247–1258.https://doi.org/10.1038/ncb3226
-
Controlling cell shape changes during salivary gland tube formation in DrosophilaSeminars in Cell & Developmental Biology 31:74–81.https://doi.org/10.1016/j.semcdb.2014.03.020
-
Outside-in signaling--a brief review of GPCR signaling with a focus on the Drosophila GPCR familyJournal of Cell Science 128:3533–3542.https://doi.org/10.1242/jcs.175158
-
Regulation and function of Scr, exd, and hth in the Drosophila salivary glandDevelopmental Biology 217:362–374.https://doi.org/10.1006/dbio.1999.9560
-
Reassessing the role and dynamics of nonmuscle myosin II during furrow formation in early Drosophila embryosMolecular Biology of the Cell 15:838–850.https://doi.org/10.1091/mbc.e03-06-0440
-
Genetic Control of Salivary Gland Tubulogenesis in DrosophilaOrganogenetic Gene Networks 1:125–149.https://doi.org/10.1007/978-3-319-42767-6_5
-
ConferencePersonal communication to FlyBaseSpokony insertions, Personal communication to FlyBase.
-
Primordium specific requirement of the homeotic gene fork head in the developing gut of the Drosophila embryoRoux’s Archives of Developmental Biology 198:201–210.https://doi.org/10.1007/BF00375906
Article and author information
Author details
Funding
Medical Research Council (U105178780)
- Yara E Sánchez-Corrales
- Katja Röper
Wellcome Trust (207553/Z/17/Z)
- Guy B Blanchard
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors would like to thank the following people; for help with an experiment: Ghislain Gillard; for reagents and fly stocks: Debbie Andrew, Thomas Lecuit, Stefan Luschnig, Herbert Jäckle, and Naoyuki Fuse. The authors thank members of the lab for input on the manuscript. KR and YSC were supported by the Medical Research Council (file reference number U105178780). GB was supported by a Wellcome Trust Investigator Award (207553/Z/17/Z) to Bénédicte Sanson, Department of PDN, Cambridge, UK.
Copyright
© 2021, Sánchez-Corrales et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,057
- views
-
- 156
- downloads
-
- 7
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
- Developmental Biology
Eukaryotic cells depend on exocytosis to direct intracellularly synthesized material toward the extracellular space or the plasma membrane, so exocytosis constitutes a basic function for cellular homeostasis and communication between cells. The secretory pathway includes biogenesis of secretory granules (SGs), their maturation and fusion with the plasma membrane (exocytosis), resulting in release of SG content to the extracellular space. The larval salivary gland of Drosophila melanogaster is an excellent model for studying exocytosis. This gland synthesizes mucins that are packaged in SGs that sprout from the trans-Golgi network and then undergo a maturation process that involves homotypic fusion, condensation, and acidification. Finally, mature SGs are directed to the apical domain of the plasma membrane with which they fuse, releasing their content into the gland lumen. The exocyst is a hetero-octameric complex that participates in tethering of vesicles to the plasma membrane during constitutive exocytosis. By precise temperature-dependent gradual activation of the Gal4-UAS expression system, we have induced different levels of silencing of exocyst complex subunits, and identified three temporarily distinctive steps of the regulated exocytic pathway where the exocyst is critically required: SG biogenesis, SG maturation, and SG exocytosis. Our results shed light on previously unidentified functions of the exocyst along the exocytic pathway. We propose that the exocyst acts as a general tethering factor in various steps of this cellular process.
-
- Cell Biology
- Developmental Biology
In most murine species, spermatozoa exhibit a falciform apical hook at the head end. The function of the sperm hook is not yet clearly understood. In this study, we investigate the role of the sperm hook in the migration of spermatozoa through the female reproductive tract in Mus musculus (C57BL/6), using a deep tissue imaging custom-built two-photon microscope. Through live reproductive tract imaging, we found evidence indicating that the sperm hook aids in the attachment of spermatozoa to the epithelium and facilitates interactions between spermatozoa and the epithelium during migration in the uterus and oviduct. We also observed synchronized sperm beating, which resulted from the spontaneous unidirectional rearrangement of spermatozoa in the uterus. Based on live imaging of spermatozoa-epithelium interaction dynamics, we propose that the sperm hook plays a crucial role in successful migration through the female reproductive tract by providing anchor-like mechanical support and facilitating interactions between spermatozoa and the female reproductive tract in the house mouse.