Skip to main content
Advertisement
  • Loading metrics

EGFR signaling coordinates patterning with cell survival during Drosophila epidermal development

Abstract

Extensive apoptosis is often seen in patterning mutants, suggesting that tissues can detect and eliminate potentially harmful mis-specified cells. Here, we show that the pattern of apoptosis in the embryonic epidermis of Drosophila is not a response to fate mis-specification but can instead be explained by the limiting availability of prosurvival signaling molecules released from locations determined by patterning information. In wild-type embryos, the segmentation cascade elicits the segmental production of several epidermal growth factor receptor (EGFR) ligands, including the transforming growth factor Spitz (TGFα), and the neuregulin, Vein. This leads to an undulating pattern of signaling activity, which prevents expression of the proapoptotic gene head involution defective (hid) throughout the epidermis. In segmentation mutants, where specific peaks of EGFR ligands fail to form, gaps in signaling activity appear, leading to coincident hid up-regulation and subsequent cell death. These data provide a mechanistic understanding of how cell survival, and thus appropriate tissue size, is made contingent on correct patterning.

Author summary

In many tissues, defective cells are eliminated by a process called apoptosis. This process prevents the emergence of rogue cells, which could be detrimental to normal physiology. Apoptosis is particularly apparent in developing embryos that lack appropriate positional information, and it has been suggested that in the absence of clear positional instructions, cells are unable to acquire a defined fate and commit suicide as a result. Here, we have used mutant fruit fly embryos lacking essential segmental determinants to identify the molecular signals that trigger apoptosis in response to mispatterning. We found that cells do not trigger apoptosis in response to conflicting fate determinants. Instead, mispatterning disrupts a tissue size control system that removes excess cells in oversized segments. Specifically, correct patterning information leads to the segmentally repeated production of survival signals, which activate the epidermal growth factor receptor, and this system is disrupted in patterning mutants leading to reproducible patterns of apoptosis. We propose that a similar, though less obvious, process also occurs in normal embryos. In such embryos, each segment would initially comprise a slight excess of cells and would then be trimmed down to a size specified by the pattern of survival signal production and the range of these signals. We suggest that a similar regulatory logic could ensure the coordination of tissue patterning and size in a variety of developing tissues.

Introduction

Defective cells are often eliminated by apoptosis during development and tissue homeostasis [14]. This has been particularly well studied during the process of cell competition, whereby unfit cells are eliminated when confronted with normal cells within a growing tissue [5]. Excess apoptosis is also seen in mutants that lack essential developmental determinants, a phenomenon that has been observed in a variety of model organisms, including zebrafish embryos lacking the signaling molecule Sonic Hedgehog [6], mice lacking the negative Wnt signaling regulator Adenomatous polyposis coli (APC) in the developing neural crest [7], and Drosophila segmentation mutants [812]. These observations have suggested the existence of a quality control system that detects conflicting or nonsense patterning inputs and, as a result, initiates apoptosis in response.

Although apoptosis was first observed in Drosophila segmentation mutants over 30 years ago [810], relatively little is known about the molecular basis of cell elimination in this context. Previous studies have revealed that apoptosis occurs within and around the areas where the missing developmental regulators are normally required. Thus, in the pair-rule mutant odd-skipped (odd), apoptosis is seen in alternate stripes that encompass the areas where odd is normally expressed [13,14]. Likewise, in embryos lacking the anterior determinant Bicoid, ectopic cell death occurs primarily at the anterior of the embryonic epidermis [13]. In Drosophila, apoptosis is initiated by a double inhibition mechanism: the activity of ubiquitously expressed inhibitor of apoptosis proteins (IAPs), which prevent caspase activation, is inhibited by the so-called IAP antagonists reaper (rpr), head involution defective (hid), grim, and sickle (skl) [15]. Most apoptotic events are initiated by the transcriptional up-regulation of one or more IAP antagonists. Indeed, hid is up-regulated in a pattern that prefigures apoptosis in segmentation mutants, i.e., at the anterior of bicoid mutants and in alternate stripes in odd mutants [13]. Therefore, it appears that cells within the segmental pattern “know” that they are missing key fate determinants and activate hid expression in response. Here, we take advantage of the genetic tools in Drosophila to investigate how apoptosis is triggered in mispatterned epidermal cells.

Results

Hid mediates apoptosis in Drosophila segmentation mutants

Consistent with an earlier report [13], extensive apoptosis was detected in Drosophila embryos lacking genes acting at various steps of the segmentation cascade, including mutants of the terminal gene tailless (tll), the gap gene krüppel (kr), the pair-rule gene fushi-tarazu (ftz), or the segment polarity gene hedgehog (hh) (Fig 1A–1E). In each instance, cleaved Death caspase-1 (Dcp1), a marker of apoptosis, was strongly enriched in the areas where the mutated patterning gene is known to be required during normal development. For example, in a null ftz mutant generated for this study (ftzΔ.attP), seven bands of apoptotic cells appeared in the segments where Ftz is normally required in the wild type (Fig 1E, see also [12]). Thus, in ftzΔ.attP mutants, alternate segments undergo massive apoptosis, while almost no cell death is detected in the remaining, normally patterned segments, which therefore serve as a useful control. For this reason, ftz mutants were selected as a prototypical condition to investigate mispatterning induced apoptosis.

thumbnail
Fig 1. Widespread apoptosis in segmentation mutants is mediated by hid.

(A–E) cDcp1 immunoreactivity (green), in control (A), tll (B), kr (C), hh (D), and ftz (E) stage 13/14 embryos. Segmental enrichment of cDcp1 is detected around the regions where the mutated gene is known to act during normal development. Arrowhead in B indicates the zone of elevated Dcp1 cleavage in the posterior epidermis of tll1 embryos. Anti-Engrailed (En, red) provides a positional reference along the A/P axis throughout. (F) Epidermal cDcp1 immunoreactivity in stage 13 hidΔ.attP ftzΔ.attP double homozygotes is strongly reduced. (G, H) Transcription of hid, as assayed with the hidΔ.GFP reporter, is detected at uniform low levels in a wild-type background (G) but is up-regulated in a banded pattern in hidΔ.GFP ftzΔ.attP double homozygotes. Embryos late stage 12. Scale bars 50 μm. (I) Quantification of cDcp1 levels in ftzΔ.attP single mutant and hidΔ.attP ftzΔ.attP double mutant embryos. (J) Quantification of mean epidermal GFP intensity values throughout the epidermis of hidΔ.GFP single mutant and hidΔ.GFP ftzΔ.attP double mutant embryos. In graphs, means are shown, and error bars display standard deviation (****p < 0.0001, ***p < 0.001, unpaired Student t test). The underlying data for (I) and (J) can be found in S1 Data. A/P, anterior/posterior; cDcp1, cleaved Death caspase-1; ftz, fushi-tarazu; GFP, green fluorescent protein; hh, hedgehog; hid, head involution defective; kr, krüppel; tll, tailless.

https://doi.org/10.1371/journal.pbio.3000027.g001

Among the four known proapoptotic genes of Drosophila, hid is the key mediator of apoptosis in bicoid mutant embryos [13]. To assess the involvement of hid in ftz mutants, we created a clean hid allele (hidΔ.attP) by replacing the first coding exon, which encodes the IAP-binding motif, with an attP integrase site [1618]. This mutation was recombined with ftzΔ.attP, and the resulting double mutant embryos were stained with anti-cleaved Dcp1. Although some signal remained in the head region, immunoreactivity was much lower than in single ftzΔ.attP mutants (Fig 1F and 1I). In contrast, abundant apoptosis was detected in ftzΔ.attP mutants lacking the closely related IAP antagonists rpr and skl (S1 Fig). We conclude that, as in bicoid mutants, hid is the main mediator of apoptosis in ftz mutant embryos. Thus, in the absence of Hid, cells that would normally be eliminated are forced to contribute to the final pattern. To assess the fate of these “undead cells”, hidΔ.attP ftzΔ.attP double mutants were allowed to undergo terminal differentiation into first instar larvae, and the pattern of denticles was examined. Alternate denticle belts were missing while numerous nondescript ectopic denticles appeared at the posterior of each remaining belt (S2 Fig). In addition, the cuticle appeared particularly convoluted, suggesting that tissue size is at least partially restored compared to the single mutant. Therefore, as previously suggested for bicoid mutants [13], it appears that undead cells can contribute to the epidermis but are unable to adopt their normal fate.

To visualize the expression of hid in mispatterned embryos, we generated an authentic transcriptional reporter by integrating a cDNA encoding green fluorescent protein (GFP) into the attP site of hidΔ.attP [16]. These animals will be referred to as hidΔ.GFP to indicate that this genetic modification creates a null allele (Δ) as well as a fluorescent reporter (GFP). To increase signal intensity, the hidΔ.GFP sensor was examined in homozygous conditions throughout, which has the dual benefit of doubling the copy number of the GFP reporter in the genome whilst simultaneously preventing apoptosis, which could eliminate cells before the GFP signal has time to accumulate. In homozygous hidΔ.GFP embryos, a weak GFP signal could be detected throughout the epidermis, suggesting that low-level hid expression occurs at subapoptotic levels during normal development (Fig 1G). In hidΔ.GFP ftzΔ.attP double homozygotes, increased levels of GFP fluorescence were detected in broad stripes, resembling the bands of caspase immunoreactivity seen in ftzΔ.attP single mutants (Fig 1H and 1J). Time-lapse imaging showed that striped GFP fluorescence arises around stage 11 of embryogenesis (S1 Movie), approximately 4 hours after the time when pair-rule genes contribute to segmental patterning [19,20]. Likewise, cleaved Dcp1 immunoreactivity became only detectable at stage 11 in homozygous ftzΔ.attP embryos (S3A–S3C Fig). Therefore, loss of ftz appears to trigger hid expression and apoptosis with a delay, well after Ftz has fulfilled its patterning role in wild-type embryos. The same delay was noted in bicoid mutants [13], suggesting that it could be a general feature of segmentation mutants.

EGFR signaling and hid expression are anticorrelated in the embryonic epidermis

Mismatched or nonsense cell fates could conceivably be recognized through local cell interactions. With the aim of identifying the relevant mediator, we used the ubiquitous actin5C-Gal4 driver to overexpress modulators of major signaling pathways in a ftzΔ.attP mutant background and stained the resulting embryos with anti-cleaved Dcp1 (summarized in S1 Table). One manipulation, overexpression of an activated form of the Drosophila epidermal growth factor receptor (EGFRact, [21]), significantly reduced cleaved Dcp1 immunoreactivity (Fig 2A–2D). Since EGFR signaling is required for survival of the cells that secrete naked cuticle [22] and since EGFR signaling has been previously implicated in suppressing hid activity [23,24], we set out to determine whether the excess apoptosis and up-regulation of hid in ftz mutants could be accounted for by a loss of EGFR activity.

thumbnail
Fig 2. Loss of EGFR signaling and hid up-regulation in ftz mutants.

(A–D) cDcp1 immunoreactivity in control w1118 (A), ftzΔ.attP (B), and ftzΔ.attP actin5C-Gal4, UAS–EGFRact (C) stage 13/14 embryos. (D) Quantification of mean cDcp1 levels for genotypes shown in A–C. Error bars display standard deviation (***p < 0.001, **p < 0.01, unpaired Student t test). (E, F) dpERK/GFP immunoreactivity in stage 13 hidΔ.GFP single mutant (E, E′) and hidΔ.GFP ftzΔ.attP double mutant embryos (F, F′). Dashed line indicates representative region of interest for quantification in panels G and H. Scale bars 50 μm. (G, H) Mean dpERK (red) and GFP (green) fluorescence intensity profiles in control hidΔ.GFP (G) and hidΔ.GFP ftzΔ.attP (H) samples. Engrailed expression stripes were used as a spatial marker along the A/P axis (see material and methods). Shaded areas depict standard error of the mean. The underlying data for (D), (G) and (H) can be found in S1 Data. cDcp1, cleaved Death caspase-1; dpERK, phosphorylated Extracellular signal–regulated kinase; EGFR, epidermal growth factor receptor; ftz, fushi-tarazu; GFP, green fluorescent protein; hid, head involution defective; UAS, upstream activation sequence; w1118, white1118.

https://doi.org/10.1371/journal.pbio.3000027.g002

If EGFR signaling is needed to keep hid expression off, we would expect hid expression to be up-regulated in EGFR mutants. Indeed, strong ubiquitous GFP fluorescence was observed in embryos homozygous for both EGFRF3, an amorphic allele, and hidΔ.GFP (S4A and S4B Fig). This suggests that most, if not all, epidermal cells require EGFR signaling to keep hid expression low. Conversely, the prosurvival activity of EGFR signaling appears to be largely mediated by repression of hid activity since almost no apoptosis was detected in EGFRF3 hidΔ.attP double mutants (S4D Fig). In EGFRF3 single mutant embryos, cell death was widespread (S4C Fig) but did not occur until embryonic stage 11 (S3D–S3F Fig), the same time at which apoptosis becomes detectable in ftz mutants. We conclude that, in the absence of EGFR signaling, most cells of the epidermis up-regulate hid expression and undergo apoptosis around mid-embryogenesis.

To visualize EGFR signaling activity in mispatterned ftz embryos, we stained embryos with an antibody that recognizes phosphorylated extracellular signal–regulated kinase (dpERK) [25]. As validation, we first stained EGFRF3 homozygotes. No immunoreactivity was detectable above background levels (S5 Fig), a strong indication that EGFR is the principal contributor of ERK phosphorylation in the embryonic epidermis. Anti-dpERK was next used to compare the patterns of EGFR signaling activity in ftzΔ.attP mutant and control embryos. The hidΔ.GFP allele (homozygous) was included to allow simultaneous assessment of hid expression and to avoid possible confounding effects of apoptosis on signaling activity. In control hidΔ.GFP embryos, dpERK immunoreactivity was observed in a wide range of tissues, including the epidermis but also the peripheral nervous system, oenocytes, and the tracheal precursors (Fig 2E) [2629]. In the epidermis, dpERK immunoreactivity was detectable in nearly all cells, though not uniformly so. Signal intensity followed a broadly undulating pattern along the anterior/posterior (A/P) axis, with smooth peaks near segment boundaries and intervening shallow troughs (Fig 2G). Little GFP fluorescence was detectable in the epidermis (Fig 2E′ and 2G), confirming that hid transcription (as reported by hidΔ.GFP) remains low in normally patterned embryos (see also Figs 1G and S3A). In the absence of ftz activity, dpERK immunoreactivity dropped significantly in seven broad bands (Fig 2F and 2H), leaving only half the number of signaling peaks intact. At the same time, a complementary pattern of GFP fluorescence (from the hidΔ.GFP reporter) appeared (Fig 2F and 2H). We conclude that, in ftz mutant embryos, the transcriptional activation of hid increases in regions where EGFR signaling activity fails to reach the threshold level required for epidermal cell survival.

Other segmentation mutants, including tll, kr, and hh, also displayed anticorrelated hidΔ.GFP reporter activity and dpERK immunoreactivity (Fig 3A–3C). Importantly, the pattern of hid expression in these mutants corresponded to that of apoptosis (as shown in Fig 1A–1D). Embryos lacking another segmentation gene, patched, were particularly informative, as they had increased dpERK immunoreactivity and did not up-regulate hidΔ.GFP (Fig 3D). Crucially, they also showed no excess apoptosis, as assayed with anti-cleaved Dcp1 (S6 Fig). This observation confirms the tight correlation between EGFR signaling and the absence of hid expression. It also shows that not all instances of patterning error are associated with apoptotic cell death.

thumbnail
Fig 3. hid and dpERK are anticorrelated in various patterning mutants.

(A–D) GFP and dpERK immunoreactivity in hidΔ.GFP tll1 (A), hidΔ.GFP hhΔ.attP (B), kr1; hidΔ.GFP (C), and ptc9; hidΔ.GFP (D) stage 13 double mutants. In each instance, dpERK intensity and hidΔ.GFP expression are anticorrelated. Scale bars 50 μm. dpERK, phosphorylated extracellular signal–regulated kinase; GFP, green fluorescent protein; hh, hedgehog; hid, head involution defective; kr, krüppel; ptc, patched; tll, tailless.

https://doi.org/10.1371/journal.pbio.3000027.g003

Segmentally expressed EGF ligands coordinate apoptosis with the patterning cascade

EGFR signaling has previously been widely implicated in cell survival in the embryo and other tissues [22,30,31]. Known Drosophila EGFR ligands include the TGF-α homologues Spitz, Keren, and Gurken, which require activation by Rhomboid proteases [3235], as well as the Rhomboid-independent neuregulin-like ligand Vein [36]. To assess the nature of the ligands involved in the survival of embryonic epidermal cells, we first examined embryos homozygous for a deficiency (Df rho-1,3) spanning the rhomboid (rho-1) and roughoid (rho-3) loci. In this background, a relatively mild apoptotic phenotype was observed (S7B Fig). Likewise, homozygous vein (vn) mutants were also characterized by a relatively minor increase in the number of apoptotic cells compared to wild-type embryos (S7C Fig). In contrast, simultaneous removal of both ligand systems, in Df rho-1,3 vnL6 double homozygotes, led to widespread apoptosis (S7D Fig), consistent with the severe cuticle phenotype previously reported [22]. We conclude that multiple ligands act redundantly to elicit the EGFR signaling landscape that, at mid-embryogenesis, prevents apoptosis.

In agreement with previously described expression patterns [36,37], vn and rho-1 were found to be expressed in a generally segmental manner in wild-type embryos (Fig 4A and 4C). Thus, peaks of ligand production correspond roughly to the areas of high-dpERK immunoreactivity seen around segmental borders (Fig 2G). We therefore considered whether changes in ligand production could account for the observed changes in EGFR signaling in patterning mutants. The pattern of rho-1 expression is expected to depend on the segmentation network since it is determined by segment polarity genes [37,38], which are themselves regulated by upstream pair-rule and gap genes. It is likely that the segmentation network similarly controls vn expression. Therefore, as expected, peaks of vn and rho-1 expression failed to appear in alternate segments of ftzΔ.attP mutants (Fig 4B and 4D), foreshadowing the loss of dpERK immunoreactivity (Fig 2F), hid up-regulation (Fig 1H), and apoptosis (Fig 1E). Importantly, the same logic can explain the pattern of apoptosis in other segmentation mutants, which are all expected to affect the pattern of EGFR ligand production within their realm of action. It also neatly explains the absence of apoptosis in patched (ptc) mutants, in which ectopic peaks of rho-1 expression arise [37].

thumbnail
Fig 4. Patterned expression of vn and rho-1 ensure epidermal cell survival.

(A, B) Expression of vn (green) in control hidΔ.attP (A) and hidΔ.attP ftzΔ.attP (B) stage 12 homozygotes, as detected by fluorescent in situ hybridisation. hidΔ.attP is included in the background to avoid any confounding effects of apoptosis. (C, D) Expression of rho-1 (green) in hidΔ.attP (C) and hidΔ.attP ftzΔ.attP (D) stage 12 homozygotes, as detected by fluorescent in situ hybridisation. Scale bars 50 μm. (E) Representative raw mSPIM image of a stage 12 w1118 embryo (control) stained with anti-cDcp1 and anti-En (red). Anterior is to the left and posterior to the right, with the ventral midline running horizontally through the middle. (F) Processed version of the image shown in E. cDcp1-positive cells were segmented and reduced to individual pixels (green). Similarly, stripes of en expression were reduced to a skeleton (red) to allow the position of each apoptotic cell to be mapped along the A/P axis of the segment and thus to derive a relative position applicable to all segments. (G) Histogram displaying the frequency of apoptosis along the A/P axis of the segment. Illustration of the en, rho-1, and vn expression domains is included to indicate the approximate sites of ligand production. Error bars show standard deviation. (H) Model explaining the pattern of apoptosis in ftz and other segmentation mutants. In normal development, the segmental expression of EGFR ligands maintains signaling activity above the threshold level required for cell survival (green dashed line). When patterning errors disrupt the landscape of EGFR activity, signaling falls below the survival threshold, triggering hid up-regulation and patterned apoptosis. The underlying data for (G) can be found in S1 Data. A/P, anterior/posterior; cDcp1, cleaved Death caspase-1; EGFR, epidermal growth factor receptor; en, engrailed; ftz, fushi-tarazu; hid, head involution defective; mSPIM, multiview single plane illumination microscopy; rho-1, rhomboid; vn, vein; w1118, white1118.

https://doi.org/10.1371/journal.pbio.3000027.g004

During normal development, EGFR signaling is lowest in cells located around the middle of each segment, which are likely to be exposed to the lowest amounts of ligand (Fig 2G). We asked whether these cells are particularly susceptible to apoptosis. To map the pattern of apoptosis along the A/P axis, we performed in toto multiview single plane illumination microscopy (mSPIM) on wild-type embryos stained with anti-cleaved Dcp1 and anti-Engrailed, which marks the anterior edge of each segment boundary. This approach allowed us to visualize, at a given time point, all epidermal apoptotic events, which are relatively rare during normal development (see Fig 1A and Fig 2A). Raw image files were generated (Fig 4E), segmented (Fig 4F), and the position of each caspase-positive cell was then mapped and normalized relative to the nearest anterior and posterior segment boundary. These data were then tabulated in a histogram displaying the frequency of apoptosis at various positions along the A/P axis within each segment (Fig 4G). Cell death was relatively infrequent near segment boundaries, where vn and rho-1 are most highly expressed (Fig 4A–4D), while apoptosis was more abundant in the middle of the intervening areas. We conclude that, in the cells that are most distant from a source of ligand, EGFR signaling activity falls near, and occasionally below, the threshold level needed to protect against apoptosis. Cells in these regions are therefore more likely to commit to an apoptotic cell fate.

Despite the observed correlation between the distance to segment boundaries and the rate of apoptosis in wild-type embryos, it is clear that cells located near the sources of EGFR ligand production can still undergo apoptosis. In this context, it is worth pointing out that only a subset of embryonic apoptotic events are regulated by the EGFR–Hid axis, since residual apoptosis is still detected in hid mutants [39], and Hid-independent apoptotic processes have been characterized in the Drosophila embryo [40]. Therefore, it is likely that additional factors, besides EGFR signaling, contribute to the control of apoptosis in the embryonic epidermis.

Discussion

Excess apoptosis in patterning mutants gives the impression that mispatterned cells are recognized as deleterious and eliminated during development. We have shown that excess apoptosis in embryonic patterning mutants can be explained by the requirement of EGFR signaling for cell survival and the fact that the segmentation cascade culminates in the segmental production of EGF ligands (model summarized in Fig 4H). Thus, in segmentation mutants, sources of EGF ligand production are lost, leaving remaining sources unable to reach all intervening cells. This results in tissue loss until all cells are brought back within signaling range, restoring a “normal” distance between segment boundaries. EGFR signaling has previously been linked to compartment size regulation in the embryo [41]. Specifically, it was proposed that epidermal cells in the posterior regions of each segment compete for a finite amount of Spitz emanating from an anterior source. Our results are consistent with an extended version of this model, whereby multiple EGFR ligands (including Spitz and Vein) produced primarily at either side of the segment boundary control the size of entire segments, not just posterior compartments. The pattern of apoptosis in wild-type embryos (Fig 4G) suggests that segments are slightly oversized upon completion of patterning and that the extra cells are eliminated as a result of subthreshold EGFR signaling and hid activation. Accordingly, there is no need to invoke the existence of a mechanism that recognizes and eliminates mis-specified cells, and mispatterning-induced cell death could be the byproduct of a size control system.

As we have shown, the lack of EGFR signaling triggers hid expression only at stage 11, when proliferation is largely complete [42]. This would serve as a particularly suitable time for a size control checkpoint and pruning of excess cells would be an effective means of ensuring reproducible segment dimensions in a fast-growing embryo with no compensatory proliferation. Similar prosurvival signaling by limited ligand diffusion from specific locations could form the basis of tissue size sensing and explain mispatterning-induced apoptosis in a variety of developmental contexts. It is intriguing that, while EGFR signaling over a relatively low threshold is required for cell survival, higher signaling activity (as occurs near the segment boundary) contributes to the specification of pattern elements (i.e., denticles) [22,37,4345]. Thus, the same ligands contribute to both patterning and cell survival, perhaps ensuring a tight coordination of these two essential developmental activities.

Materials and methods

Drosophila stocks and husbandry

Flies were maintained at 25°C on standard fly food containing yeast, agar, and cornmeal. w1118 was used as a wild-type control throughout. Fly stocks used in this study are listed in Table 1.

To identify regulators of apoptosis in mispatterned cells, a collection of UAS lines was assembled. Each expressed a gene that modulates the activity of a signaling pathway in the presence of Gal4 (S1 Table). These so-called UAS lines were crossed to the null ftz allele ftzΔ.attP, and male progeny containing both the ftz allele and the UAS transgene were crossed to females bearing the ubiquitous actin5C-Gal4 driver recombined to ftzΔ.attP. The embryos produced from this cross were then stained for cleaved Dcp1.

Generation of transgenic flies

Accelerated “ends out” homologous recombination was used to generate the hidΔ.attP and ftzΔ.attP alleles (as described in [18]). Homology arms of approximately 2 kb were PCR amplified from genomic DNA. These fragments were cloned into a targeting vector (pTVCherry) containing a GMR > white eye marker and an attP landing site flanked with FRT sites. UAS-rpr was included downstream of the 3′ homology arm. Targeting vectors were integrated into the genome via P element–mediated transformation to generate so-called donor lines. To release the homology cassette from its genomic location, donor flies were crossed to hs-FLP and heat shocked at 24-hour intervals throughout larval development. Adult female progeny with mottled eye color (indicating the presence of both the targeting vector and hs-FLP constructs) were crossed to the strong ubi-Gal4 driver to eliminate individuals that had either failed to excise the homology cassette or had undergone an unsuccessful recombination event. Resultant white+ candidates were then isolated and verified with PCR.

To generate the hidΔ.GFP reporter, a GFP cDNA was cloned into the RIVCherry vector [18], which contains an attB sequence and a 3xPax3-cherry selection marker. This construct was injected into hidΔ.attP embryos along with a source of ΦC31 integrase, and progeny were screened for Cherry expression. Rainbow Transgenic Flies, Inc. (Camarillo, California, United States) and Bestgene, Inc. (Chino Hills, California, United States) were used for embryo injection services.

Embryo collection

Embryos were collected overnight, transferred to baskets, and washed with water. Embryos were then dechorionated using 75% bleach for 3 minutes and washed again. After drying on a paper towel, embryos were transferred with a brush to microcentrifuge tubes containing equal volumes of heptane and 6% formaldehyde in phosphate buffered saline (PBS). Fixation was performed for 25 minutes on a rotating platform. The fixative was then removed and replaced with 1 ml methanol. Embryos were shaken vigorously for 45 seconds to remove the vitelline membrane, and the heptane was removed. Embryos were washed twice more in methanol and moved to −20°C for long-term storage.

Immunofluorescence

Fixed embryos were rehydrated in 0.3% triton in PBS (PBT) and blocked for 30 minutes at room temperature in 4% bovine serum albumin in PBT. Embryos were incubated with primary antibody at 4°C overnight, washed for 1 hour in PBT at room temperature, and then incubated in a solution of Alexa Fluor 488, 568, or 633 conjugated species appropriate secondary antibodies (Thermo Fisher Scientific, 1:1,000) for 2 hours. Further washes (1 hour) were conducted before samples were mounted in Vectashield containing DAPI (Vector Laboratories). Primary antibodies used in this study are listed in Table 2.

Confocal microscopy

Confocal microscopy was conducted with a Leica SP5 confocal microscope using a 20x (Leica, NA 0.7) or 40x (Leica, NA 1.25) oil immersion objective. Typically, 8 confocal planes were imaged at 1.8 μm intervals and processed using Fiji (NIH) and Photoshop (Adobe).

Live confocal imaging

Samples were prepared using the hanging drop method described in [47]. A small drop of oxygen permeable Voltalef 10S oil was spotted onto a cover slip. A dechorionated embryo was then placed inside the oil and oriented with a 27-G hypodermic needle. The sample was then inverted so that the oil containing the sample hung below the cover slip. Images were acquired from above at 10 minute intervals using a Leica SP5 confocal microscope and processed using imageJ (NIH) and Photoshop (Adobe) as described above.

mSPIM imaging

Fluorescently labeled samples were imaged on a multiview light sheet microscope [48]. The optics consisted of two detection and illumination arms. Here, each detection arm formed an epifluorescence microscope using in sequence: a water-dipping lens (Apo LWD 25x, NA 1.1, Nikon Instruments, Inc.), a filter wheel (HS-1032, Finger Lakes Instrumentation LLC) equipped with emission filters (BLP01-488R-25, BLP02-561R-25, Semrock, Inc.), a tube lens (200 mm, Nikon Instruments, Inc.), and sCMOS camera (Hamamatsu Flash 4.0 v2.0). The imaging produced an effective pixel size of 0.26 μm. The illumination arms each consisted in sequence: a water-dipping objective (CPI Plan Fluor 10x, NA 0.3, Nikon Instruments, Inc.), a tube lens (200 mm, Nikon Instruments, Inc.) a scan lens (S4LFT0061/065, Sill optics GmbH & Co. KG), and a galvanometric scanner (6215h, Cambridge Technology, Inc.). The illumination arms were fed by a combination of lasers (06-MLD 488 nm, Cobolt AB, and 561 LS OBIS 561 nm, Coherent, Inc.). Samples were translated through the resulting light sheet using a linear piezo stage (P-629.1CD together with E-753 controller). Multiple rotation views were achieved using a rotational piezo stage (U-628.03 with C-867 controller) in combination with a linear actuator (M-231.17 with C-863 controller, all Physik Instrumente GmbH and Co. KG). Microscope operation was done using Micro Manager [49], running on a Super Micro 7047GR-TF Server, with 12 Core Intel Xeon 2.5 GHz, 64-GB PC3 RAM, and hardware Raid 0 with 7 2.0 TB SATA drives. Four views with a rotational offset of 90° were recorded at 1 μm sectioning. Fusion of individual views is based on a diagnostic specimen used to determine an initial guess for an affine transformation that maps each view into a common frame of reference. Together with the data, this initial guess was used with a rigid image registration algorithm [50] to fuse individual views and resulted in isotropic resolution of 0.26 μm in the final registered image. ImSAnE [51] was used to obtain tissue cartographic projections showing global maps of the embryo surface of interest.

Fluorescent in situ hybridisation

Digoxigenin labeled antisense RNA probes were synthesized from Gold Collection (Berkeley Drosophila Genome Project) clones containing full length versions of the vn or rho-1 cDNAs (clone number LP21849 and LD06131, respectively). Vectors were linearized with the EcoRI restriction enzyme and used as a template for in vitro transcription as per manufacturer instructions (Roche, SP6/T7 DIG RNA-labeling kit). Fixed embryos were bathed in 3% hydrogen peroxide in methanol for 15 minutes to remove endogenous peroxidase activity and rinsed extensively in PBT. Embryos were then incubated for 3 minutes in 10 μg/ml proteinase K in PBT and washed in 2 mg/ml Glycine in PBT. Samples were next washed in a hybridization buffer before the relevant probe was added, and the samples incubated overnight at 55°C. The next day, the embryos were washed once more in PBT and incubated overnight at 4°C with a sheep anti-DIG antibody diluted at 1:2,000. On the third day, embryos were washed and incubated with a biotinylated anti-sheep secondary antibody, and the final signal was developed using a Tyramide-FITC signal amplification kit (Thermo Fisher Scientific). Standard immunostaining was then performed on these samples, using the protocol outlined above.

Quantification and statistical analysis

To measure fluorescence intensity along the A/P axis, rectangular regions of interest with a fixed height of 75 μm were drawn over four contiguous segments. In each instance, efforts were made to ensure the same embryonic regions were analyzed between samples. GFP (from hidΔ.GFP), dpERK, and Engrailed intensity profiles were generated across the regions of interest in Fiji [52], and these data were imported into MATLAB (Mathworks) for further analysis. GFP and dpERK intensity profiles were fitted to the peaks of engrailed expression (i.e., segment boundaries) along the A/P axis to normalize data to segment width, which can vary between samples. This was achieved by creating 50 evenly spaced intervals along the A/P axis of each segment before extracting the raw data values and averaging the resulting numbers across all samples analyzed. standard error of the mean was calculated for each point, and data was plotted using the boundedline.m function developed by Kelly Kearney and downloaded from the MATLAB stack exchange.

To measure the frequency of apoptosis along the segmental A/P axis, the Engrailed, and cleaved Dcp1 channels were separated from the raw mSPIM image files and then segmented using the open source iLastik software (University of Heidelberg). Segmented images were imported into ImageJ and reduced to single pixel lines and points using the software’s inbuilt skeletonization algorithms. From these skeletons, individual segment boundaries were manually classified and masks were generated that covered each of the segments to be analyzed. Processed files were then transferred to MATLAB for further analysis. Distances were calculated from each Dcp-1 positive pixel to the nearest anterior and posterior Engrailed positive pixel using a k-nearest neighbors algorithm, and these measurements were used to determine the position of each apoptotic cell along the segmental A/P axis as a percentage of segment width. Cells located in the immediate vicinity of the ventral midline (defined as 20% of linear Dorsal to Ventral distance) were excluded from the analysis to avoid apoptotic figures in the developing nervous system. The retained data were collated and tabulated as a histogram.

Supporting information

S1 Fig. Apoptosis in ftz mutants does not require rpr.

(A) Schematic representation of the Drosophila H99 locus showing the genes encoding the four major IAP antagonists. The Df(3L)XR38 deficiency, which removes the rpr and skl genes, is highlighted in blue. (B, C) Cleaved Dcp1 immunoreactivity in stage 12 ftzΔ.attP (B) and ftzΔ.attP Df(3L)XR38 homozygotes (C). Scale bars 50 μm. Df, deficiency; ftz, fushi-tarazu; IAP, inhibitor of apoptosis protein; rpr, reaper; skl, sickle.

https://doi.org/10.1371/journal.pbio.3000027.s001

(TIF)

S2 Fig. “Undead cells” undergo terminal differentiation in ftz mutants that cannot undergo apoptosis.

(A, B) Cuticle preparations of ftzΔ.attP (A) and hidΔ.attP ftzΔ.attP (B) embryos. In ftzΔ.attP homozygotes, remaining denticle belts have near wild-type morphology. In hidΔ.attP ftzΔ.attP homozygotes, the remaining denticle belts are expanded posteriorly with multiple rows of nondescript ectopic denticles (B′). ftz, fushi-tarazu; hid, head involution defective.

https://doi.org/10.1371/journal.pbio.3000027.s002

(TIF)

S3 Fig. Apoptosis in ftz and EGFR mutants is first detectable at embryonic stage 11.

(A–C) Cleaved Dcp1 immunoreactivity in ftzΔ.attP mutants at embryonic stages 9 (A), 11 (B), and 13+ (C). (D–F) Cleaved Dcp1 immunoreactivity in EGFRF3 mutants at embryonic stages 9 (D), 11 (E), and 13+ (F). In both genotypes, cleaved Dcp1 is first detected in stage 11 and persists throughout the remainder of embryonic development. Scale bars 50 μm. Dcp1, Death caspase-1; EGFR, epidermal growth factor receptor; ftz, fushi-tarazu.

https://doi.org/10.1371/journal.pbio.3000027.s003

(TIF)

S4 Fig. EGFR signaling maintains epidermal cell survival through repression of hid.

(A, B) Transcription of hid, as assayed with the hidΔ.GFP reporter, in control (A) and EGFRF3 mutant (B) stage 12 embryos. hid is up-regulated in most epidermal cells upon loss of EGFR signaling. (C, D) Cleaved Dcp1 immunoreactivity is strongly up-regulated throughout the epidermis in stage 12 EGFRF3 single mutants (C) and this up-regulation is lost in stage 12 EGFRF3; hidΔ.attP double homozygotes (D). Scale bars 50 μm. Dcp1, Death caspase-1; EGFR, epidermal growth factor receptor; GFP, green fluorescent protein; hid, head involution defective.

https://doi.org/10.1371/journal.pbio.3000027.s004

(TIF)

S5 Fig. dpERK immunoreactivity is markedly reduced upon loss of EGFR signaling.

(A, B) dpERK immunoreactivity in control (EGFR+/+) hidΔ.attP (A) and EGFRF3; hidΔ.attP double homozygotes. Extensive dpERK immunoreactivity is detected in wild-type control embryos (A) and this signal is largely lost in EGFRF3; hidΔ.attP double mutants (B). We take this as evidence that EGFR signaling is the major source of ERK phosphorylation in the embryonic epidermis. Scale bars 50 μm. dpERK, phosphorylated extracellular signal–regulated kinase; EGFR, epidermal growth factor receptor; hid, head involution defective.

https://doi.org/10.1371/journal.pbio.3000027.s005

(TIF)

S6 Fig. Absence of ectopic apoptosis in ptc mutants.

(A, B) Cleaved Dcp1 immunoreactivity in control w1118 (A) and ptc9 mutant embryos (B) at embryonic stage 13. No increase in Dcp1 cleavage was detected in ptc mutants, despite disruption to the segmental pattern. Scale bars 50 μm. Dcp1, Death caspase-1; ptc, patched; w1118, white1118.

https://doi.org/10.1371/journal.pbio.3000027.s006

(TIF)

S7 Fig. Vein and Rhomboid act redundantly to maintain epidermal cell survival.

(A–D) Cleaved Dcp1 in control w1118 (A), Df rho-1,3 (B), vnL6 (C), and Df rho-1,3 vnL6 double homozygotes at embryonic stage 13 (D). A mild increase in Dcp1 immunoreactivity is seen in vn and rho-1,3 single mutants (compared to w1118 controls). This signal is strongly enhanced in the double mutants. Scale bars 50 μm. Dcp1, Death caspase-1; Df, deficiency; rho-1, rhomboid; rho-3, roughoid; vn, vein; w1118, white1118.

https://doi.org/10.1371/journal.pbio.3000027.s007

(TIF)

S1 Table. A genetic survey to identify regulators of apoptosis in patterning mutants.

List of pathways targeted and transgenes used to assess the rescue of apoptosis in ftz mutants. In each instance, the transgene of interest was expressed globally using the actin5C-Gal4 driver and a positive score was attributed if a reduction in the levels of cleaved Dcp1 occurred. Dcp1, Death caspase-1; ftz, fushi-tarazu.

https://doi.org/10.1371/journal.pbio.3000027.s008

(XLSX)

S1 Movie. hid is up-regulated at mid-embryogenesis in ftz mutants.

Activity of the hidΔ.GFP reporter in a ftzΔ.attP mutant embryo. Negligible fluorescence is detected during the early stages of embryogenesis but around embryonic stage 11 (approximately 7 hours after egg laying) bands of fluorescence appear. Hours after egg lay are displayed in the lower right corner. ftz, fushi-tarazu; GFP, green fluorescent protein; hid, head involution defective.

https://doi.org/10.1371/journal.pbio.3000027.s009

(AVI)

Acknowledgments

We thank Alberto Baena-Lopez for sharing fly stocks and reagents and Donald Bell for imaging support. We also acknowledge the advice and help of Ingrid Poernbacher, Kristina Stapornwongkul, Alberto Baena-Lopez, Cyrille Alexandre, Nic Tapon, Iris Salecker, and Jorge Beira during the early phase of this research. We finally thank Barry Thompson, Ingrid Poernbacher, and Ian McGough for comments on the manuscript.

References

  1. 1. Merino Marisa M, Rhiner C, Lopez-Gay Jesus M, Buechel D, Hauert B, et al. (2015) Elimination of Unfit Cells Maintains Tissue Health and Prolongs Lifespan. Cell 160: 461–476. pmid:25601460
  2. 2. Koto A, Kuranaga E, Miura M (2011) Apoptosis ensures spacing pattern formation of Drosophila sensory organs. Curr Biol 21: 278–287. pmid:21276725
  3. 3. Fuchs Y, Steller H (2011) Programmed Cell Death in Animal Development and Disease. Cell 147: 742–758. pmid:22078876
  4. 4. Derry WB, Putzke AP, Rothman JH (2001) Caenorhabditis elegans p53: role in apoptosis, meiosis, and stress resistance. Science 294: 591–595. pmid:11557844
  5. 5. Merino MM, Levayer R, Moreno E (2016) Survival of the Fittest: Essential Roles of Cell Competition in Development, Aging, and Cancer. Trends Cell Biol 26: 776–788. pmid:27319281
  6. 6. Prykhozhij SV (2010) In the Absence of Sonic Hedgehog, p53 Induces Apoptosis and Inhibits Retinal Cell Proliferation, Cell-Cycle Exit and Differentiation in Zebrafish. PLoS ONE 5: e13549. pmid:21042410
  7. 7. Hasegawa S, Sato T, Akazawa H, Okada H, Maeno A, et al. (2002) Apoptosis in neural crest cells by functional loss of APC tumor suppressor gene. Proceedings of the National Academy of Sciences 99: 297.
  8. 8. Ingham PW, Howard KR, Ish-Horowicz D (1985) Transcription pattern of the Drosophila segmentation gene hairy. Nature 318: 439.
  9. 9. Jackle H, Rosenberg UB, Preiss A, Seifert E, Knipple DC, et al. (1985) Molecular analysis of Kruppel, a segmentation gene of Drosophila melanogaster. Cold Spring Harb Symp Quant Biol 50: 465–473. pmid:3007002
  10. 10. Martinez-Arias A (1985) The development of fused- embryos of Drosophila melanogaster. J Embryol Exp Morphol 87: 99–114. pmid:3928800
  11. 11. White RA, Lehmann R (1986) A gap gene, hunchback, regulates the spatial expression of Ultrabithorax. Cell 47: 311–321. pmid:2876779
  12. 12. Magrassi L, Lawrence PA (1988) The pattern of cell death in fushi tarazu, a segmentation gene of Drosophila. Development 104: 447–451. pmid:3256471
  13. 13. Werz C, Lee TV, Lee PL, Lackey M, Bolduc C, et al. (2005) Mis-specified cells die by an active gene-directed process, and inhibition of this death results in cell fate transformation in Drosophila. Development 132: 5343–5352. pmid:16280349
  14. 14. DiNardo S, O'Farrell PH (1987) Establishment and refinement of segmental pattern in the Drosophila embryo: spatial control of engrailed expression by pair-rule genes. Genes Dev 1: 1212–1225. pmid:3123316
  15. 15. Goyal L, McCall K, Agapite J, Hartwieg E, Steller H (2000) Induction of apoptosis by Drosophila reaper, hid and grim through inhibition of IAP function. The EMBO Journal 19: 589–597. pmid:10675328
  16. 16. Bateman JR, Lee AM, Wu CT (2006) Site-specific transformation of Drosophila via phiC31 integrase-mediated cassette exchange. Genetics 173: 769–777. pmid:16547094
  17. 17. Wright CW, Clem RJ (2002) Sequence requirements for Hid binding and apoptosis regulation in the baculovirus inhibitor of apoptosis Op-IAP. Hid binds Op-IAP in a manner similar to Smac binding of XIAP. J Biol Chem 277: 2454–2462. pmid:11717313
  18. 18. Baena-Lopez LA, Alexandre C, Mitchell A, Pasakarnis L, Vincent JP (2013) Accelerated homologous recombination and subsequent genome modification in Drosophila. Development 140: 4818–4825. pmid:24154526
  19. 19. Clark E, Akam M (2016) Odd-paired controls frequency doubling in Drosophila segmentation by altering the pair-rule gene regulatory network. eLife 5: e18215. pmid:27525481
  20. 20. DiNardo S, O'Farrell PH (1987) Establishment and refinement of segmental pattern in the Drosophila embryo: spatial control of engrailed expression by pair-rule genes. Genes & development 1: 1212–1225.
  21. 21. Lesokhin AM, Yu S-Y, Katz J, Baker NE (1999) Several Levels of EGF Receptor Signaling during Photoreceptor Specification in Wild-Type,Ellipse,and Null MutantDrosophila. Developmental Biology 205: 129–144. pmid:9882502
  22. 22. Urban S, Brown G, Freeman M (2004) EGF receptor signalling protects smooth-cuticle cells from apoptosis during Drosophila ventral epidermis development. Development 131: 1835–1845. pmid:15084467
  23. 23. Bergmann A, Agapite J, McCall K, Steller H (1998) The Drosophila Gene hid Is a Direct Molecular Target of Ras-Dependent Survival Signaling. Cell 95: 331–341. pmid:9814704
  24. 24. Kurada P, White K (1998) Ras Promotes Cell Survival in Drosophila by Downregulating hid Expression. Cell 95: 319–329. pmid:9814703
  25. 25. Gabay L, Seger R, Shilo BZ (1997) In situ activation pattern of Drosophila EGF receptor pathway during development. Science 277: 1103–1106. pmid:9262480
  26. 26. Reichman-Fried M, Dickson B, Hafen E, Shilo BZ (1994) Elucidation of the role of breathless, a Drosophila FGF receptor homolog, in tracheal cell migration. Genes & Development 8: 428–439.
  27. 27. Sutherland D, Samakovlis C, Krasnow MA (1996) branchless encodes a Drosophila FGF homolog that controls tracheal cell migration and the pattern of branching. Cell 87: 1091–1101. pmid:8978613
  28. 28. Elstob PR, Brodu V, Gould AP (2001) spalt-dependent switching between two cell fates that are induced by the Drosophila EGF receptor. Development 128: 723–732. pmid:11171397
  29. 29. Bier E, Jan LY, Jan YN (1990) rhomboid, a gene required for dorsoventral axis establishment and peripheral nervous system development in Drosophila melanogaster. Genes & Development 4: 190–203.
  30. 30. Bergmann A, Tugentman M, Shilo BZ, Steller H (2002) Regulation of cell number by MAPK-dependent control of apoptosis: a mechanism for trophic survival signaling. Dev Cell 2: 159–170. pmid:11832242
  31. 31. Gilboa L, Lehmann R (2006) Soma-germline interactions coordinate homeostasis and growth in the Drosophila gonad. Nature 443: 97–100. pmid:16936717
  32. 32. Urban S, Lee JR, Freeman M (2002) A family of Rhomboid intramembrane proteases activates all Drosophila membrane-tethered EGF ligands. EMBO J 21: 4277–4286. pmid:12169630
  33. 33. Rutledge BJ, Zhang K, Bier E, Jan YN, Perrimon N (1992) The Drosophila spitz gene encodes a putative EGF-like growth factor involved in dorsal-ventral axis formation and neurogenesis. Genes Dev 6: 1503–1517. pmid:1644292
  34. 34. Reich A, Shilo BZ (2002) Keren, a new ligand of the Drosophila epidermal growth factor receptor, undergoes two modes of cleavage. EMBO J 21: 4287–4296. pmid:12169631
  35. 35. Shmueli A, Cohen-Gazala O, Neuman-Silberberg FS (2002) Gurken, a TGF-alpha-like protein involved in axis determination in Drosophila, directly binds to the EGF-receptor homolog Egfr. Biochem Biophys Res Commun 291: 732–737. pmid:11866425
  36. 36. Schnepp B, Grumbling G, Donaldson T, Simcox A (1996) Vein is a novel component in the Drosophila epidermal growth factor receptor pathway with similarity to the neuregulins. Genes Dev 10: 2302–2313. pmid:8824589
  37. 37. Alexandre C, Lecourtois M, Vincent J (1999) Wingless and Hedgehog pattern Drosophila denticle belts by regulating the production of short-range signals. Development 126: 5689–5698. pmid:10572045
  38. 38. Gritzan U, Hatini V, DiNardo S (1999) Mutual antagonism between signals secreted by adjacent wingless and engrailed cells leads to specification of complementary regions of the Drosophila parasegment. Development 126: 4107–4115. pmid:10457019
  39. 39. Grether ME, Abrams JM, Agapite J, White K, Steller H (1995) The head involution defective gene of Drosophila melanogaster functions in programmed cell death. Genes & Development 9: 1694–1708.
  40. 40. Lohmann I, McGinnis N, Bodmer M, McGinnis W (2002) The Drosophila Hox Gene Deformed Sculpts Head Morphology via Direct Regulation of the Apoptosis Activator reaper. Cell 110: 457–466. pmid:12202035
  41. 41. Parker J (2006) Control of Compartment Size by an EGF Ligand from Neighboring Cells. Current Biology 16: 2058–2065. pmid:17055987
  42. 42. Edgar BA, O’Farrell PH (1990) The Three Postblastoderm Cell Cycles of Drosophila Embryogenesis Are Regulated in G2 by string. Cell 62: 469–480. pmid:2199063
  43. 43. O'Keefe L, Dougan ST, Gabay L, Raz E, Shilo BZ, et al. (1997) Spitz and Wingless, emanating from distinct borders, cooperate to establish cell fate across the Engrailed domain in the Drosophila epidermis. Development 124: 4837–4845. pmid:9428420
  44. 44. Szuts D, Freeman M, Bienz M (1997) Antagonism between EGFR and Wingless signalling in the larval cuticle of Drosophila. Development 124: 3209–3219. pmid:9272961
  45. 45. Payre F, Vincent A, Carreno S (1999) ovo/svb integrates Wingless and DER pathways to control epidermis differentiation. Nature 400: 271. pmid:10421370
  46. 46. Zeidler MP, Perrimon N, Strutt DI (1999) Polarity determination in the Drosophila eye: a novel role for unpaired and JAK/STAT signaling. Genes Dev 13: 1342–1353. pmid:10346822
  47. 47. Reed BH, McMillan SC, Chaudhary R (2009) The preparation of Drosophila embryos for live-imaging using the hanging drop protocol. J Vis Exp.
  48. 48. Krzic U, Gunther S, Saunders TE, Streichan SJ, Hufnagel L (2012) Multiview light-sheet microscope for rapid in toto imaging. Nat Methods 9: 730–733. pmid:22660739
  49. 49. Edelstein AD, Tsuchida MA, Amodaj N, Pinkard H, Vale RD, et al. (2014) Advanced methods of microscope control using μManager software. Journal of biological methods 1: e10. pmid:25606571
  50. 50. Klein S, Staring M, Murphy K, Viergever MA, Pluim JP (2010) elastix: a toolbox for intensity-based medical image registration. IEEE Trans Med Imaging 29: 196–205. pmid:19923044
  51. 51. Heemskerk I, Streichan SJ (2015) Tissue cartography: compressing bio-image data by dimensional reduction. Nat Methods 12: 1139–1142. pmid:26524242
  52. 52. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, et al. (2012) Fiji: an open-source platform for biological-image analysis. Nat Methods 9: 676–682. pmid:22743772