mRNA poly(A)-tail changes specified by deadenylation broadly reshape translation in Drosophila oocytes and early embryos

1The generation of a zygote through fertilization marks the onset of development for the offspring, and requires the proper completion of sperm and oocyte development in the respective parents. Spermatogenesis produces haploid sperm capable of penetrating the oocyte, whereas oogenesis produces differentiated oocytes that are stockpiled with maternal nutrients, proteins, and mRNAs, and have outer layers that protect the embryo and enable fertilization. In most animals the oocyte is arrested in meiosis, and fertilization leads to initiation of mitosis as the oocyte nucleus completes meiosis and fuses with the haploid sperm nucleus. Fertilization also leads to changes in mRNA translation and protein stability, which support a period of development driven off of maternal stockpiles. After this initial stage of maternal control, which lasts for 1–2 mitotic divisions in mammals and 13 mitotic divisions in Drosophila, widespread transcription begins from the zygotic nuclei.

2The production of viable offspring requires three key developmental events: oocyte maturation, the oocyte-to-embryo transition (OET), and the maternal-to-zygotic transition (MZT) (Figure 1A). Oocyte maturation involves the release of the primary meiotic arrest at prophase I and progression of the oocyte nucleus into meiotic divisions to produce a mature oocyte (egg) capable of being fertilized. The mature oocyte is arrested in metaphase II in most vertebrates and metaphase I in insects. The OET requires egg activation, which is coupled to fertilization in most animals. Egg activation releases the secondary meiotic arrest, allowing completion of meiosis. The MZT marks the transfer of control of development from the mother to the zygote as maternal mRNAs are degraded, transcription from the zygotic genome begins, and embryonic development becomes dependent on zygotic gene products. In Xenopus, zebrafish, and Drosophila the major activation of zygotic transcription occurs as the cell cycle lengthens and gastrulation begins, a developmental period referred to as the midblastula transition.

3In nearly all animals, oocyte maturation and the OET occur in the absence of transcription, and thus any changes in protein levels must occur through posttranscriptional regulatory mechanisms. Indeed, widespread changes in protein levels without corresponding changes in mRNA levels are observed at these developmental transitions. Posttranscriptional regulation of gene expression during oocyte maturation, the OET, and early embryogenesis faces unique challenges, including stably storing maternal mRNAs through the prolonged periods of oogenesis and maintaining some of these mRNAs in a translationally inactive state (referred to as masking) while selectively translating others at the proper times in development. During mouse oocyte maturation ~1300 mRNAs are translationally upregulated and ~1500 are translationally downregulated. Similarly, during Drosophila egg activation ~1000 mRNAs are translationally upregulated, and ~500 are translationally repressed.

4One long-appreciated mechanism of translational regulation during development acts through control of mRNA poly(A)-tail length. The observation of short poly(A) tails in oocytes from amphibians and marine invertebrates led to the proposal that short poly(A) tails help mask maternal mRNAs, preventing translation without triggering the destabilization that typically results from deadenylation. Specific mRNAs are then targeted for cytoplasmic polyadenylation, and the lengthened poly(A) tails in turn cause translational upregulation of these mRNAs, as demonstrated for the mRNAs for c-mos and several cyclin genes during oocyte maturation in Xenopus. The Cytoplasmic Polyadenylation Element Binding Protein (CPEB) plays a critical role in controlling translation of maternal mRNAs, both in recruiting Maskin to repress translation and later in promoting polyadenylation by binding cytoplasmic poly(A) polymerase together with Cleavage and Polyadenylation Specificity Factor (CPSF). In Drosophila, the cytoplasmic poly(A) polymerase encoded by the wispy gene promotes poly(A)-tail lengthening during both oocyte maturation and egg activation, and increased poly(A)-tail lengths of mRNAs for cyclin B, cortex, bicoid, torso, and Toll coincide with the translational upregulation of these mRNAs. Moreover, appending synthetic poly(A) tails of increasing lengths to the bicoid mRNA promotes the translation of this mRNA, which demonstrates a causal relationship between increased poly(A)-tail length and increased translation in early Drosophila embryos.

5Translational repressors, some of which have been identified through genetic analysis of pattern formation in the Drosophila embryo, also act during oogenesis and early embryogenesis. These repressors can impact either poly(A)-tail length, the translational machinery, or both. Brain Tumor and Pumilio (PUM) bind to target mRNAs and recruit the CCR4-NOT deadenylase complex, which shortens poly(A) tails. Bruno binds to target mRNAs during oogenesis and recruits the eIF4E-binding protein Cup, which represses translation by both blocking the interaction between eIF4E and eIF4G and promoting deadenylation while repressing decapping, which traps mRNAs in a translationally repressed state. Similarly, the RNA-binding protein Smaug (SMG) is reported to recruit either Cup or the CCR4-NOT deadenylase complex, with the consequent repression of translation. smg mRNA is translationally repressed in the mature oocyte by PUM, and at egg activation this repression is relieved by the PAN GU (PNG) protein kinase complex. PNG is also thought to regulate other translational regulators.

6mRNA expression arrays, RNA-seq and ribosome-footprint profiling have enabled global analyses of mRNA levels and relative translational efficiencies (TEs) in Drosophila oocytes and embryos. Measuring poly(A)-tail lengths has been more challenging due to the homopolymeric nature of the tail, and this technical hurdle has restricted studies to measuring tail lengths of a small number of Drosophila mRNAs or to estimating relative differences in tail lengths more broadly using hybridization-based methods. However, methods are now available to measure poly(A)-tails of essentially any length for millions of individual mRNA molecules. In zebrafish and Xenopus, high-throughput measurements demonstrate that poly(A)-tail length and TE are globally coupled during early embryonic development, as expected from single-gene studies. Surprisingly, however, this coupling disappears at gastrulation and is not observed in any of the post-embryonic samples examined. Thus, soon after a developmental shift in transcriptional control (the MZT), these vertebrate species undergo a developmental shift in translational control, in which TE abruptly changes from being strongly coupled to tail length to being uncoupled.

7Drosophila melanogaster affords the opportunity to define changes in poly(A)-tail length and TE during oocyte maturation, the OET, and the MZT, facilitated by the high quality of both the genome assembly and the coding-sequence (CDS) annotations, the wealth of knowledge on regulatory proteins active during early development, the availability of genetic mutants in these known regulatory proteins, and the ability to isolate sufficient quantities of oocytes and embryos at distinct developmental stages for high-throughput sequencing techniques. Here we systematically analyze changes in poly(A)-tail length and TE during oocyte maturation, the OET and the MZT, and delineate the roles of Wispy, PNG, and SMG in regulating translation during these crucial periods. These studies reveal fundamental insights into the regulatory mechanisms and regimes operating during these key developmental transitions.

8We performed poly(A)-tail length profiling by sequencing (PAL-seq) to measure poly(A)-tail lengths globally at eight developmental stages in Drosophila melanogaster, beginning with oocytes at the primary meiotic arrest point and ending with embryos late in gastrulation, at germ-band retraction (Figure 1A). PAL-seq determines the poly(A)-tail length of millions of mRNA molecules, providing for each molecule both a fluorescent measurement that reports on the poly(A)-tail length and a sequencing read that reports the cleavage and poly(A) site. Because alternative start sites or alternative splicing can generate different transcripts with the same poly(A) site, we considered our results with respect to a curated set of unique gene models. Moreover, because alternative cleavage and polyadenylation can generate different poly(A) sites for mRNAs from the same gene, we used our PAL-seq data to identify for each gene the alternative isoform with the longest 3' UTR and then combined its tail-length measurements with those of any shorter tandem 3'-UTR isoforms. Thus, the 'mRNAs' of our analyses each corresponded to a unique gene model representing a different protein-coding gene and the aggregated tail-length results from many (if not all) of the alternative transcripts from that gene. On the same samples as those analyzed by PAL-seq, we performed ribosome-footprint profiling and RNA-seq to determine TEs, which also report aggregate results for alternative isoforms, allowing comparison of mean poly(A)-tail length and TE for the thousands of mRNAs that passed our quantification cutoffs at each stage.

9In Drosophila, the first 13 division cycles occur in a nuclear syncytium. Cellularization (at 2–3 hr of embryogenesis) marks completion of the MZT and is immediately followed by the morphogenetic movements of gastrulation. Thus, early embryogenesis is controlled by maternal mRNAs, and based on the studies of mRNAs from a few genes, those with longer poly(A)-tails would be predicted to have higher TEs, as observed in vertebrate systems. Indeed, poly(A)-tail length and TE were correlated in cleavage-stage embryos (0–1 hr embryo), and the slope of this relationship was steep, such that small differences in poly(A)-tail length corresponded to large differences in TE (Figure 1B; e.g., the median TE of mRNAs with a poly(A)-tail length between 70–80 nt was >5 fold greater than that of mRNAs with a poly(A)-tail length between 30–40 nt). Likewise, when progressing through the OET, the changes in TE correlated with the changes in tail length (Figure 1B, 0–1 hr embryo / Stage 14 oocyte).

10Although a positive correlation between poly(A)-tail length and TE persisted in the 2–3 hr embryos, the slope of this relationship strongly diminished such that the median TE of mRNAs with a poly(A)-tail length between 70–80 nt was <2 fold greater than that of mRNAs with a poly(A)-tail length between 30–40 nt, as expected if poly(A)-tail length changes were no longer used to influence translation in these cellularized and gastrulating embryos. Indeed, the changes in poly(A)-tail length that occurred between cleavage-stage and cellularized embryos were not positively correlated with the changes in TE (Figure 1B, 2–3 hr embryo / 0–1 hr embryo). The positive correlation between poly(A)-tail length and TE disappeared in subsequent stages of gastrulation (3–4 hr embryos and 5–6 hr embryos), and although tail lengths changed between these stages, the TE changes spanned a very narrow range (Figure 1B, 3–4 hr embryo / 2–3 hr embryo and 5–6 hr embryo / 3–4 hr embryo). The tail-length changes observed between these stages presumably resulted from transcription and mRNA decay as well as deadenylation, whereas those in earlier stages resulted primarily from cytoplasmic polyadenylation and deadenylation.

11The narrower range of TE values and TE changes following the onset of gastrulation (Figure 1B) was consistent with a diminished role of translational regulation in specifying protein output after the embryo begins to modulate gene expression by altering mRNA abundance. Most importantly, these results showed that the developmental shift in the nature of translational control observed in vertebrates also occurs in flies; in Drosophila, just as in zebrafish and Xenopus, the coupling between poly(A)-tail length and TE observed in the early embryo disappeared around gastrulation.

12An analysis of the temporal patterns of translational regulation occurring from the OET through the onset of gastrulation identified several significantly enriched patterns, although all possible patterns were observed (Figure 1—figure supplement 1). Among the five enriched profiles were those for mRNAs that were translationally upregulated at the OET and then stayed elevated during early embryonic development, those for mRNAs that were upregulated at the OET and then downregulated during early embryonic development, and those for mRNAs that were downregulated before the OET and remained repressed (Figure 1—figure supplement 1). Because ribosome-profiling data reflect relative rather than absolute levels of translation of an mRNA across developmental stages, these results should be interpreted as reflecting changes in how well an mRNA was translated relative to all other mRNAs in the same stage. Thus, an mRNA translated at the same absolute efficiency in two stages might be perceived as downregulated if translation of many other mRNAs increased. Nonetheless, these dynamic and substantial changes in TE observed during the OET and early embryogenesis revealed a large-scale and rapid change in the translational capacity of mRNAs relative to each other during early development.

13Pairwise comparisons within the mRNA, translation, and tail-length datasets (from RNA-seq, ribosome-footprint profiling, and PAL-seq, respectively) for the different developmental stages further supported these overall conclusions (Supplementary file 1). As expected, mRNA abundances were highly correlated between all stages prior to the onset of transcription (Supplementary file 1). Over these same stages, poly(A)-tail lengths for many different mRNAs shortened or lengthened, but these changes seemed to have minimal impact on mRNA abundance measurements, consistent with stability of shorter-tail mRNAs during this time (and also indicating that the poly(A)-selection that was performed as part of the RNA-seq protocol did not bias our RNA abundance, and thus TE, measurements). In contrast to the highly correlated mRNA abundances, the number of ribosome-protected fragments (RPFs) differed substantially between stages prior to the onset of transcription (Supplementary file 1), consistent with dynamic and robust translational regulation during this time. Once transcription began, widespread changes in both mRNA levels and RPFs were observed between stages (Supplementary file 1). Poly(A)-tail lengths varied both prior to and after the onset of transcription, and generally mirrored the pattern of correlations observed for RPFs up until the MZT (Supplementary file 1). Little to no positive correlation between poly(A)-tail lengths was observed when comparing samples with no active transcription to those after the onset of robust transcription (Supplementary file 1), presumably because of the influence of nascent transcripts on mean poly(A)-tail lengths.

14The striking increases in TE observed soon after fertilization presumably generated factors needed for later developmental transitions. For example, the TE of mRNAs for zelda and stat92E, two transcription factors required for the onset of zygotic transcription at the MZT, increased markedly during the OET (Figure 1C), as expected if translational activation of these mRNAs produces the proteins necessary to initiate the MZT, and accounting for the previously noted increase in Zelda protein in early embryos. Likewise, the TE of smaug and bicoid mRNA increased markedly during the OET (Figure 1C), yielding factors essential for embryonic development and patterning, respectively. Although the TE of stat92E and smaug began to increase at stage 14, their TEs were still relatively low until after the OET (Figure 1C), consistent with the roles of their translation products during early embryonic development.

15Having established that Drosophila embryos resemble those of fish and frog with respect to the coupling between TE and poly(A)-tail length, we sought to extend these findings to oocytes. The ability to isolate individual oocytes at distinct developmental stages in Drosophila enabled us to measure poly(A)-tail lengths and TEs and examine their relationship at key stages of oogenesis. Stage 11 oocytes are arrested in prophase I. Nuclear envelope breakdown, a hallmark of the release of this arrest, occurs between stages 12 and 13, and the oocytes arrest in metaphase I during stage 14, the final stage of oogenesis. Thus, comparison of oocytes from stages 11 through 14 captures the events of meiotic progression and oocyte differentiation that occur over the course of oocyte maturation.

16Poly(A)-tail lengths and TEs were well correlated at stages 11, 12, and 13 (Rs ≥0.51), but poorly correlated at stage 14 (Rs = 0.20) (Figure 2). Previously, the relationship between poly(A)-tail length and translation had been examined for the c-mos mRNA and several cyclin mRNAs during Xenopus oocyte maturation. Our results extended these analyses to an invertebrate and to a global scale, showing that, as for the c-mos mRNA and several cyclin mRNAs in the Xenopus oocyte, tail-length increases correspond to translational activation in the maturing Drosophila oocyte. Thus, the global coupling between poly(A)-tail length and TE, which had previously been examined only after fertilization, extends well before fertilization to potentially influence key developmental transitions in the oocyte.

17Of the mRNAs exceeding our read cutoffs, nearly half (46%) had net TE changes of ≥2 fold over the course of oocyte maturation, with mRNAs for 1264 and 1037 genes up- or downregulated, respectively (Figure 2—figure supplement 1A). In contrast, far fewer mRNAs (1.5%) changed in abundance ≥2 fold (25 mRNAs increased, and 52 mRNAs decreased), indicating the posttranscriptional nature of this regulation. More detailed comparison of TEs for each mRNA from stages 11 through 14 revealed unexpectedly dynamic and complex temporal patterns of translational control (Figure 3A). The relative TEs for some mRNAs increased progressively at each stage of oogenesis (Figure 3A cluster 1), and this set of mRNAs included cyclin B (Figure 3B), whose translational control via polyadenylation is reported to be crucial for oocyte maturation. Some of the downregulated mRNAs similarly showed progressive decreases in TE through the stages of oocyte maturation (Figure 3A cluster 26). However, more complex patterns of developmental changes in TE were more commonly observed (Figure 3A). Moreover, changes in TE correlated poorly with previously measured changes in protein abundance, potentially due to a large role of posttranslational regulation during this time (Figure 2—figure supplement 1B, Rs = 0.15).

18The net changes in TE and poly(A)-tail length over the course of oocyte maturation were largely concordant for both translationally up- and downregulated mRNAs, indicating strong coupling between TE and poly(A)-tail length during oocyte maturation (Figure 4A). Nonetheless, a potential deviation from the typical coupled behavior occurred between oocyte stages 13 and 14, when a set of mRNAs that underwent substantial tail-length reductions did not undergo concordant TE reductions (Figure 2, Stage 14 / Stage 13, and Figure 2—figure supplement 2A). These mRNAs tended to have both long poly(A)-tails and efficient translation in stage 13 oocytes, and then remained well translated compared to other mRNAs in stage 14 oocytes, despite undergoing substantial tail shortening (Figure 2—figure supplement 2B). Notable among this set were mRNAs for 22 of the 30 genes encoding subunits of the proteasome (Figure 2—figure supplement 2A and B).

19To examine the relationship between tail-length changes and TE changes for mRNAs of key factors in oocyte maturation, we focused on those of three of the most highly upregulated cell-cycle regulators involved in meiotic progression, Cyclins B and B3 and the Cdc20 ortholog Fizzy. Between oocyte stages 11 through 14, the TEs of these mRNAs increased 103, 248, and 319 fold (compared to a median TE change of only 1.3 fold), respectively, whereas poly(A)-tail lengths increased 1.6, 2.4, and 2.0 fold (compared to the median of 1.2 fold) (Figure 4A). In addition to these examples, mRNAs from 121 other genes were very poorly translated in stage 11 oocytes (TE ≤ 0.1) and became ≥10-fold better translated during oocyte maturation (Figure 2—figure supplement 1C). As with cyclin B, cyclin B3 and fizzy mRNAs, these also had modest increases in tail length (Figure 2—figure supplement 1C; median tail-length fold change = 1.7 for mRNAs with TE ≤ 0.1 in stage 11 oocytes and a ≥10-fold increase in TE in stage 14 oocytes), raising the question of the degree to which modest tail-length increases might explain the large TE increases observed during oocyte maturation.

20One possibility is that a subset of the mRNAs for a gene underwent extensive polyadenylation, and this change in a subset of mRNAs did not cause much change in the mean tail-length value despite causing a large change in TE. However, for cyclin B, cyclin B3 and fizzy mRNAs, the distribution of poly(A)-tail lengths remained unimodal as overall lengths increased, which indicated that the focus on the mean tail-length change in our analyses did not miss specific regulation of a subset of mRNAs for a gene (Figure 4B). Also helpful for considering the effects of increasing tail lengths during oocyte maturation is the steep relationship between poly(A)-tail length and TE during this developmental period (Figure 2). The steepness of this relationship in stage 11, 12, and 13 oocytes implied a 6- to 7-fold increase in the median TE for mRNAs with a poly(A)-tail length between 70–80 nt relative to those with a poly(A)-tail length between 30–40 nt. In stage 14 oocytes, the poor correlation between tail length and TE resulted in only a 1.8-fold increase in median TE when comparing mRNAs with such tail lengths. Thus during oocyte maturation, the modest increases observed in mean poly(A)-tail lengths presumably explain part, but not all, of the robust translational activation of cyclin B, cyclin B3, fizzy, and mRNAs from many other genes. Moreover, a broad range of TE values was observed for mRNAs with the same mean poly(A)-tail length, presumably as a consequence of additional regulatory processes that are independent of poly(A)-tail length. To the extent that these additional processes act concordantly with tail-length–dependent processes, they would increase the apparent effect of tail lengthening.

21Although egg activation occurs rapidly (within ~20 min), this developmental transition induced TE changes for mRNAs from many genes, some of which exceeded 30 fold (Figure 2, lower right). As in the preceding period of oocyte maturation, tail-length changes corresponded to TE changes (Figure 4A and C). Indeed, during egg activation the correlation between tail-length change and TE change exceeded that observed for any other transition (Figure 1B and Figure 2, lower panels; Figure 4A). Thus during egg activation, changes in poly(A)-tail length can account for a large proportion of the translational changes.

22mRNAs for the embryonic patterning genes bicoid, torso, and Toll are translationally upregulated during egg activation in a manner reported to depend on changes in poly(A)-tail length. In our analyses, translation of bicoid mRNA was undetectable or detected with very few reads in oocytes from stages 11 through 14, but it became much more highly translated following egg activation, increasing by 36 fold in the activated egg relative to the stage 14 oocyte (median TE change, 1.2 fold). mRNAs for torso and Toll were somewhat better translated in stage 14 oocytes, but nonetheless their TEs increased 13 and 46 fold, respectively (Figure 4C). During this time, the Toll poly(A)-tail length increased 3.3 fold (median tail length change, 1.1 fold), which corresponded well to the trend observed globally for tail-length and TE changes (Figure 4C), suggesting that its increased tail-length might fully account for its increased TE. However, concordantly acting tail-length–independent processes that increase the apparent effect of tail-length changes might also be influencing Toll translation. Indeed, our subsequent experiments with wispy-mutant animals showed that in the absence of cytoplasmic polyadenylation, Toll underwent substantial translational upregulation during egg activation and remained one of the most efficiently translated mRNAs in the wispy-mutant laid egg, despite having a mean tail length of only 14 nt. The poly(A)-tail lengths of bicoid and torso increased by only 1.7 and 1.3 fold, respectively, putting them outside the global trend (Figure 4C), and the unimodal tail-length distributions showed again that our focus on the mean change did not miss a unique behavior of a subpopulation of transcripts (Figure 4D). Thus, the translational regulation of these mRNAs seems to also involve tail-length–independent mechanisms.

23The strong translational upregulation of bicoid mRNA, from essentially untranslated in oocytes from stages 11 through 14 to relatively well translated in the activated egg, implied the influence of either a translational repressor acting on bicoid mRNA in the oocyte that was relieved at egg activation, or a translational activator acting on it during egg activation. Because bicoid mRNA had an average tail length of 42 nt in stage 14 oocytes, which was similar to that of many well translated mRNAs at this stage, we conclude that there was likely a translational repressor acting on it in oocytes to prevent its translation prior to egg activation. The same is presumably true for other mRNAs that had little to no detectable translation in oocytes but then became efficiently translated following egg activation. For mRNAs such as torso, which had somewhat greater translation in stage 14 oocytes but subsequently experienced a TE increase that was out of proportion to its tail-length increase, either possibility might account for its translational upregulation at egg activation. Likewise, the 15-fold and 4-fold translational activation observed between stage 14 oocytes and 0–1 hr embryos for zelda and stat92E, respectively (Figure 1C), were accompanied by only 1.2-fold changes in poly(A)-tail length for both mRNAs (median TE change, 1.2 fold; median tail length change, 1.1 fold), again implying that a component of their activation may have occurred through a tail-length independent mechanism (note that the poly(A)-tail length of zelda was quantified by only 31 and 72 tags in stage 14 oocytes and 0–1 hr embryos, respectively).

24During oocyte maturation and egg activation, changes in poly(A)-tail length occur through the opposing activities of deadenylases and cytoplasmic poly(A)-polymerases, which shorten and lengthen poly(A) tails, respectively. Wispy is the cytoplasmic poly(A)-polymerase responsible for poly(A)-tail lengthening at the end of oocyte maturation and continuing through egg activation. Indeed, within this developmental window, most maternal transcripts bind to oligo(dT) with less affinity in wispy mutants, indicating that many poly(A) tails shorten in the absence of Wispy, but the degree to which tail-lengths shorten and the impact on TE had not been previously measured.

25Our PAL-seq measurements of poly(A)-tail lengths in wispy-mutant stage 13 oocytes and laid eggs confirmed the widespread shortening of poly(A) tails in the absence of Wispy (Figure 5A). Of the mRNAs from 3571 genes measured in stage 13 oocytes, 1801 had tails that were at least 50% shorter in wispy-mutant oocytes relative to wild-type, and for all but a few mRNAs the poly(A) tail was at least somewhat shorter in the absence of Wispy (Figure 5A). Although wispy-mutant eggs do not complete meiosis and thus do not begin nuclear divisions, most of the molecular changes characteristic of the OET are independent of these divisions, motivating a comparison of wispy-mutant laid eggs with wild-type early embryos. In wispy-mutant laid eggs, 1219 of the 1836 measured mRNAs had poly(A) tails at least 50% shorter than those in the analogous wild-type sample (0–1 hr embryos), and again for all but a few mRNAs the poly(A) tail was shorter in the absence of Wispy (Figure 5A). Examination of the mRNAs with tails that were least affected by the absence of Wispy revealed a remarkable enrichment for mRNAs of ribosomal proteins, both in oocytes and laid eggs (Figure 5—figure supplement 1A; p<10–41 and p<10–44, respectively, one-tailed Wilcoxon rank sum test).

26Strikingly, the global changes in poly(A)-tail length that occurred in the absence of Wispy did not affect the coupling between poly(A)-tail length and TE. In wispy-mutant stage 13 oocytes and laid eggs, poly(A)-tail length and TE were correlated to a similar extent as in wild-type oocytes and embryos (Figure 1B, 2, and Figure 5—figure supplement 1B). These correlations were observed even though the RNA-seq measurements used to calculate TEs were made using poly(A)-selected RNA, which in the wispy-mutant samples likely resulted in inflated TE estimates for mRNAs with very short poly(A)-tails, as these would be poorly recovered during the poly(A) selection. Indeed, RNA-seq measurements correlated with poly(A)-tail lengths in the wispy-mutant samples but not the wild-type samples (Figure 5—figure supplement 1C), indicating a poly(A)-selection bias in the wispy-mutant RNA-seq data but not in the wild-type data. Because neither developmental stage being considered had appreciable transcription or RNA degradation, and because of the poly(A)-selection bias in the wispy-mutant RNA-seq data, we reasoned that the RNA-seq measurements from corresponding wild-type samples should more accurately reflect mRNA abundances in the wispy-mutant samples. Indeed, results obtained from determining TEs for the wispy-mutant samples using the wispy-mutant ribosome profiling and the wild-type RNA-seq (Figure 5B) were similar to those obtained when only wispy-mutant measurements were used (Figure 5—figure supplement 1B), except they had somewhat stronger correlations between poly(A)-tail length and TE, which supported the idea that these modified TEs more accurately reflected the true TEs at these stages.

27Importantly, the strong correlation between changes in poly(A)-tail length and changes in TE that occurred at egg activation in wild-type flies was retained in wispy-mutant flies (Figures 1B and 5C, respectively), indicating that changes in the length of the poly(A) tail—not the act of cytoplasmic polyadenylation per se—is what influences TE. Moreover, the changes in TE that occurred during egg activation in wispy-mutant samples strongly correlated with those in the wild-type samples, as did the relative changes in poly(A)-tail length (Figure 5D). Thus, although Wispy lengthens the poly(A)-tails of most mRNAs during egg activation, our results indicated that this lengthening was largely dispensable for reshaping TEs during egg activation, as the rank order of tail-length and TE changes at this time were largely preserved in the wispy-mutant samples.

28The observation that cytoplasmic polyadenylation was not needed to generate the rank order of tail-length and TE changes for most mRNAs warranted further investigation, particularly when considering both the critical influence of cytoplasmic polyadenylation for generating these differences during vertebrate development and the importance of Wispy for relieving meiotic arrest during the OET. Perhaps some mRNAs most affected by the loss of Wispy dropped out of our tail-length analysis because they had very short poly(A)-tails in the wispy mutant and thus were poorly recovered in our PAL-seq libraries. Although the tail lengths of such mRNAs might have been difficult to measure by PAL-seq, the RPF measurements of mRNAs for all genes can be analyzed without consideration of the corresponding tail-length or mRNA abundance measurements, either of which might selectively lose the set of mRNAs most dependent on Wispy for their regulation. The RPFs in wild-type and wispy-mutant stage 13 oocytes were well correlated, and the same was true for 0–1 hr wild-type embryos compared with wispy-mutant laid eggs (Rs = 0.76 for both, Figure 5—figure supplement 1D). Although globally the RPFs were largely preserved in the wispy-mutant samples, some mRNAs had many RPFs in the wild-type samples but none in the wispy-mutant samples (Figure 5—figure supplement 1D). Included among these was cortex, a known target of Wispy, the upregulation of which is essential for the OET.

29We conclude that during egg activation in Drosophila, cytoplasmic polyadenylation extends the tails of nearly all mRNAs but does not play the major role in generating the relative differences in tail lengths and TEs observed between mRNAs, which implies that poly(A)-tail shortening instead plays the major role in generating these differences. Nonetheless, the mRNAs responsible for the wispy-mutant phenotype might be among those that are most strongly dependent on Wispy for their translational state.

30The PNG kinase, which is required for nuclear divisions in the early embryo, influences the translational up- or downregulation of hundreds of different mRNAs at egg activation (Figure 6A), presumably through its regulation of translational regulators, such as SMG. To characterize the extent to which these downstream regulators might use tail-length changes to influence TE, we profiled poly(A)-tail lengths and TEs in png-mutant stage 14 oocytes and activated eggs. Poly(A)-tail lengths and TEs were highly correlated between png-mutant and wild-type stage 14 oocytes, whereas the correlation was substantially lower between png-mutant and wild-type activated eggs (Figure 6—figure supplement 1A). Additionally, the poly(A)-tail lengths and TEs in png-mutant activated eggs were quite similar to those in png-mutant stage 14 oocytes, whereas this was not the case for the same comparisons among wild-type samples (Figure 6—figure supplement 1B). The requirement of PNG for achieving the gene-expression profile of the wild-type activated egg but not the stage 14 oocyte was consistent with PNG kinase becoming active during egg activation to serve as a master regulator of the OET.

31The set of mRNAs dependent on PNG for their translational downregulation was predominantly those that underwent tail-length shortening during egg activation in wild-type samples (median TE log2 fold change of –3.03 after median centering and median tail-length log2 fold change of –0.35) (Figure 6B). In png-mutant embryos, these mRNAs failed to be deadenylated, as expected if tail-length shortening mediated their translational repression (median TE log2 fold change of 0.10 after median centering and median tail-length log2 fold change of 0.32, which closely resembled the median value for all mRNAs) (Figure 6B). Correspondingly, the set of mRNAs dependent on PNG for their translational upregulation was predominantly those that underwent tail-length extension during wild-type egg activation (median TE log2 fold change of 2.98 after median centering and median tail-length log2 fold change of 0.62), and these mRNAs failed to undergo tail lengthening in png mutants, again as expected if tail-lengthening mediated their translational activation (median TE log2 fold change of 0.36 after median centering and median tail-length log2 fold change of 0.32, which closely resembled the median value for all mRNAs) (Figure 6C). mRNAs that were translationally up- or downregulated by PNG also typically increased or decreased in poly(A)-tail length, respectively (Figure 6D), and the tail-length changes for either the up- or downregulated sets of mRNAs were largely dependent on PNG (Figure 6B and C), which showed that PNG plays a central role in regulating polyadenylation and deadenylation for the hundreds of mRNAs present at egg activation.

32Because much of the PNG-mediated translational repression is attributed to its derepression of smg translation, we examined the impact of SMG on the set of mRNAs that depended on PNG for tail-length shortening and translational repression. SMG is reported to downregulate translation of its target mRNAs through either direct recruitment of the deadenylation machinery or recruitment of CUP, which directly represses translation initiation and can also promote deadenylation. Our results showed that in the absence of SMG, the poly(A)-tail lengths for the set of PNG-dependent downregulated mRNAs were less shortened during egg activation (median tail-length log2 fold changes of –0.22 and –0.04 for the wild-type and smg mutant, respectively, relative to the median of all mRNAs, p<10–27, one-tailed Wilcoxon rank sum test), and this coincided with their partial translational derepression (median TE log2 fold change of –2.2 and –1.2, respectively, after median centering, p<10–32, one-tailed Wilcoxon rank sum test) (Figure 7A). In contrast, no significant difference was observed for either poly(A)-tail length or TE changes for the PNG-dependent upregulated mRNAs (median tail-length log2 fold changes of 0.27 for both the wild-type and smg mutant relative to the median of all mRNAs, and median TE log2 fold changes of 2.8 and 2.7, respectively, after median centering for the wild-type and smg mutant, p=0.89 and 0.37 for tail-length or TE changes, respectively, two-tailed Wilcoxon rank sum test), consistent with SMG functioning only as a repressor (Figure 7B). In general, mRNAs that were translationally repressed by SMG also decreased in poly-(A)-tail length, which indicated that recruitment of deadenylases was the dominant mechanism of SMG-dependent translational repression occurring during egg activation and early embryogenesis (Figure 7C). Additionally, in the absence of SMG, a set of mRNAs reported to be bound by SMG in the early embryo had longer poly(A) tails and were better translated following the OET (median tail-length log2 fold changes of –0.21 and –0.04 for the wild-type and smg mutant, respectively, relative to the median of all mRNAs, and median TE log2 fold changes of –1.0 and –0.15, respectively, after median centering, p<10–16 and 10–11 for tail length or TE changes, respectively, one-tailed Wilcoxon rank sum test) (Figure 7—figure supplement 1).

33Taken together, our results indicate that SMG mediates much of the PNG-dependent translational repression occurring during egg activation and early embryogenesis, with most of this repression attributable to poly(A)-tail shortening. Nevertheless, the residual repression observed in the absence of SMG (Figure 7A) indicates that one or more additional PNG-activated translational repressors also operate at this time.

34Our results from Drosophila oocytes and embryos increase the number of cellular environments known to possess global coupling between poly(A)-tail length and TE, which now include several stages of Drosophila oocytes as well as Drosophila, zebrafish, and Xenopus embryos prior to the onset of gastrulation. At least as early as the primary meiotic arrest point, changes in poly(A)-tail length largely explain changes in TE, and this regulatory strategy continues through early embryonic development. A switch around the onset of gastrulation uncouples poly(A)-tail length and TE in Drosophila, mirroring the developmental timing of the same uncoupling in zebrafish and Xenopus embryos. The uncoupling of TE from poly(A)-tail length is thus a broadly conserved developmental switch in the nature of translational control, occurring around the onset of gastrulation in both vertebrate and invertebrate embryos.

35Changing the length of a poly(A) tail would influence both the number of poly(A)-binding proteins (PABPs) that could bind the tail and the probability that at least one PABP would be bound, with the former being more relevant in conditions of saturating PABP and the later being more relevant in conditions in which tails are competing for PABP. If the number of PABPs that bind the tail were most relevant, then TE changes might have correlated better with absolute tail-length changes than with relative tail-length changes. However, repeating our analyses using absolute tail-length changes rather than tail-length log2 fold changes yielded the same results for the respective comparisons, with nearly identical Rs values (Figure 7—figure supplement 2), suggesting that PABP is limiting in the developmental stages in which tail length and TE are coupled, such that mRNAs are competing for PABP.

36We presume that the correlation between poly(A)-tail length and TE that we observed in oocytes and pre-gastrulation embryos reflects a causal relationship, i.e., that changes in poly(A)-tail length at least partially cause the corresponding changes in TE. The other possibilities are that changes in TE cause changes in tail length or that these two measurements merely correlate without having a causal relationship. Because previous single-gene studies clearly demonstrate a causal role for cytoplasmic polyadenylation, and poly(A)-tail length itself, in regulating the translation of mRNAs during these developmental times, we interpret our observation of a global correlation between poly(A)-tail length and TE as reflecting a causal relationship between the two wherein tail-length changes cause TE changes. Given this interpretation, our results indicate that widespread changes in poly(A)-tail length broadly reshape the translational profiles of oocytes and early embryos.

37Despite this strong qualitative conclusion, we cannot provide quantitative estimates of the proportion of each TE change that is attributable to each tail-length change, as additional translational regulatory mechanisms presumably act to reinforce or oppose the effects of tail-length–dependent translational regulation. Indeed, we suspect that tail-length–independent processes often act concordantly with tail-length–dependent processes and thereby amplify the apparent effects of tail-length changes. Such amplification would explain why some mRNAs, including cyclin B, cyclin B3, fizzy, zelda, stat92E, bicoid, and torso, undergo translational activation that seems out of proportion to their tail-length changes.

38A primary motivation for using Drosophila in our study was the opportunity to examine the impact of Wispy, PNG, and SMG on poly(A)-tail length and TE. Of particular interest was Wispy, the cytoplasmic poly(A) polymerase that acts during late oogenesis and early embryogenesis, and is essential for the OET. Analysis of wispy mutants revealed key insights into both the mechanism of translational activation and the importance of tail shortening (as opposed to tail lengthening) in specifying translational control in Drosophila maturing oocytes and early embryos. Although poly(A)-tail lengths were typically shorter in wispy-mutant stage 13 oocytes and laid eggs compared to the corresponding wild-type samples, tail length and TE were coupled in both of these wispy-mutant contexts, and the relative changes in poly(A)-tail lengths during the OET in wild-type samples were largely preserved in the wispy-mutant samples, as were the changes in TE. Thus, the act of cytoplasmic polyadenylation is not necessary for the coupling of poly(A)-tail length to translation. Moreover, cytoplasmic polyadenylation is largely dispensable for specifying which mRNAs are to be translationally activated or repressed, implying that in Drosophila, selective deadenylation is more important for imparting this specificity. This greater importance for deadenylation does not preclude a role for selective cytoplasmic polyadenylation in translational control in Drosophila—differences were observed when comparing TE changes in wild-type and wispy-mutant samples (Figure 5D), and when comparing RPF measurements between wild-type and wispy-mutant samples (Figure 5—figure supplement 1D), including for cortex, a critical OET regulator. However, the greater importance of selective deadenylation for specifying translational control in Drosophila differs from the prevailing paradigm of vertebrates, which centers on selective cytoplasmic polyadenylation of mRNAs with cytoplasmic polyadenylation elements. This suggests either a fundamental difference between invertebrates and vertebrates or that an underappreciated component of translational control in vertebrates occurs through more global cytoplasmic polyadenylation with selective deadenylation.

39Despite the surprisingly modest perturbation of translational regulation that we observed in wispy-mutant laid eggs during egg activation, these laid eggs are blocked at the OET because of defective progression through meiosis and failure of pronuclei to fuse. Why does this developmental arrest occur if translational regulation is largely preserved during egg activation in wispy mutants? One possibility is that the absolute level of translation in wispy-mutant eggs is reduced such that they make an inadequate amount of one or more essential protein. The other possibility rests on the observation that although TE changes and RPFs were remarkably concordant in the wispy-mutant and wild-type samples, they were not identical, and thus an aberrant change in expression for one or more key factor, such as cortex, might be sufficient to trigger the developmental arrest of wispy-mutant laid eggs.

40Activity of the PNG kinase during egg activation affects most mRNAs that undergo translational activity changes at this developmental transition, both those that are up- and downregulated. Although the direct targets of PNG kinase remain to be identified, PNG is required to relieve the repressive effects of PUM and thereby derepresses translation of smg and other mRNAs. In wild-type activated eggs the mRNAs dependent on PNG for translational up- or downregulation exhibited corresponding increases or decreases in poly(A)-tail length, and these tail-length changes did not occur in png mutants. Thus, for both translationally activated and translationally repressed mRNAs, the downstream effectors of PNG appear to regulate translation primarily through changes in poly(A)-tail length.

41Translational repressors in addition to SMG have been proposed to contribute to PNG-dependent regulation. Our results are consistent with the idea that PNG activates at least one additional translational repressor, as the set of mRNAs dependent on PNG for their translational downregulation were still somewhat translationally repressed and deadenylated relative to all other mRNAs in smg-mutant embryos, although less than in wild-type embryos. Additionally, although SMG is reported to act partly through a mechanism that does not involve deadenylation, essentially all the derepression observed in the smg-mutant embryos was accompanied by tail-length increases. Perhaps at other times SMG causes direct translational repression, but in 0–1 hr embryos it seems to act predominantly through deadenylation.

42Within the overall framework in which tail-length changes cause the up- or downregulation of mRNAs for a vast repertoire of different genes, we discovered unanticipated complexities and exceptions. One unanticipated complexity was the evidence for tail-length–independent processes that appear to amplify the tail-length–dependent effects, as already discussed. Additional unanticipated complexity was observed during oocyte maturation, in which surprisingly diverse patterns of regulatory dynamics produce the ultimate changes in TE. Single-gene studies in Xenopus, which show that c-mos mRNA becomes polyadenylated and translated prior to the mRNAs for several cyclin B genes, hint that the complex landscape of translational control observed in Drosophila extends to vertebrate oocyte maturation. The diverse temporal behaviors of translational regulation during oocyte maturation might reflect the evolving demands on the maturing oocyte, or perhaps simply differences in the timing of expression of RNA-binding proteins that interact with these different sets of mRNAs. Regardless, the observation of different temporal patterns of regulation indicates that rather than being a sudden switch, oocyte maturation in Drosophila is a progressive process, whereby mRNAs are translationally activated or repressed at different times throughout this developmental transition.

43Another surprise was the identification of a set of mRNAs for which tail-length changes did not correspond to TE changes during oocyte maturation. This exception to the general pattern occurred for the mRNAs that had undergone deadenylation between oocyte stages 13 and 14 yet continued to be translated with relatively high efficiency (Figure 2—figure supplement 2). To escape the general regulatory environment of the cytoplasm, these mRNAs might exploit unique interactions with the translational machinery analogous to the translational activation of histone mRNAs through stem-loop binding protein, or they might localize to a subcellular compartment that lacks coupling. The deadenylation of these mRNAs during oogenesis might prime them for rapid degradation during embryogenesis. Most mRNAs from genes encoding proteasome subunits fell into this class. Previous comparisons of the proteome and translation changes suggest that extensive proteolysis occurs at egg activation. Upregulation of the proteasome in late oogenesis might necessitate a mechanism to downregulate it rapidly during embryogenesis, perhaps facilitated by pre-deadenylating the corresponding mRNAs before egg activation.

44The molecular mechanism that couples TE with poly(A)-tail length is unknown, as is the mechanism of the uncoupling that occurs at gastrulation. Nevertheless, our results, which show that widespread changes in poly(A)-tail length broadly reshape translational activity during oocyte and embryo development, demonstrate the importance of this coupling mechanism and provide a resource that documents the thousands of affected genes with information on the timing, magnitude, and inferred consequences of the tail-length changes for mRNAs from each of these genes. Moreover, the identification of unanticipated complexities and exceptions reveals the rich posttranscriptional regulatory landscape of oocytes and early embryos and points to additional mechanisms that operate for some mRNAs in specific settings, which can now be targeted for investigation.

45All flies were kept at 18, 22 or 25°C according to standard procedure. Oregon R (OrR) flies were used as a wild-type control. The null png1058 allele was previously described. The wispy hemizygote mutants were derived from a cross between wisp41/FM6 (kindly provided by Mariana Wolfner, Cornell) and Df(1)RA47/FM7c flies (BL961, obtained from the Bloomington Drosophila Stock Center). The smaug hemizygote mutants resulted from a cross between smg1/TM3, Sb[1] (BL5930, obtained from the Bloomington Drosophila Stock Center) and Df(3L)Scf-R6, Diap1[1] st[1] cu[1] sr[1] e[s] ca[1]/TM3, Sb[1] flies (BL4500, obtained from the Bloomington Drosophila Stock Center). Unfertilized eggs were collected from crosses to twineHB5mutant males, which fail to make sperm.

46Egg chambers were hand-dissected in Grace’s Unsupplemented Insect Media (Gibco) from three-day old flies that had been fattened for two days with wet yeast at 22°C. Approximately 300 egg chambers were collected per oogenesis stage. Oocyte stages 11 through 14 were distinguished using morphological criteria. Activated eggs were collected as previously described. Samples were transferred to RPF lysis buffer (10 mM Tris-HCl, pH 7.4, 5 mM MgCl2, 100 mM KCl, 2 mM DTT, 100 µg ml–1 cycloheximide, 1% Triton X-100, 500 U ml–1 RNasin Plus, and protease inhibitor (1x complete, EDTA-free, Roche); wild-type stage 14 oocyte and activated egg samples were transferred to RPF lysis buffer lacking cycloheximide), homogenized and then centrifuged at 10000 rpm for 5 min.The supernatants were collected and flash frozen in liquid nitrogen to be stored at –80°C. Ultimately, the yield was approximately 40–50 μg of total RNA in 70–120 μl lysate.

47Embryos were collected at 25°C by discarding two initial one-hour collections to avoid collecting embryos held within females for prolonged time. Only the subsequent 0–1 hr collection was then processed for sequencing analysis. The same embryo collection strategy was applied to smg or wispy mutants, though for the wispy mutant the collection was extended to 0–2 hr as one hour was too short to gather a sufficient number of laid eggs for sequencing. To obtain OrR 2–3 h, 3–4 hr and 5–6 hr embryo samples, 0–1 hr embryos (collected after two previous one-hour collections were discarded) were incubated in a “humid chamber” at 25°C for additional 2, 3 or 5 hr, respectively.

48Drosophila females can hold developing embryos or lay unfertilized eggs, and so to ensure that the vast majority of the embryos processed for downstream analyses were in the expected developmental stage, a fraction of each embryo collection (approximately 10–20%) was stained with DAPI. This approach was also used to confirm that the embryos laid by mutant mothers (png, smg, or wispy) were developmentally arrested as previously described. Briefly, collected embryos were dechorionated in 50% bleach for 3 min at room temperature, extensively washed with water, dried and then transferred into a scintillation vial that contained methanol. To remove the vitelline membrane from the embryos, an equal volume of heptane was added and the vial was vigorously shaken by hand for 2 min. Next, embryos were fixed overnight in methanol at 4°C. Embryos were then sequentially rehydrated in a mix of methanol:PBS (9:1, 3:1, 1:1, 1:3) for at least 5 min per step. After a wash in PBS, embryos were incubated for 15 min in a 1 μg ml–1 solution of DAPI in PBS, followed by 30 min wash in PBS containing 0.1% Triton X-100. Finally, after removing all the washing solution, samples were mounted in Vectashield.

49Reference transcript annotations were downloaded from the UCSC Genome browser (dm6 in refFlat format), and for each gene the longest transcript was chosen as the initial representative transcript model. Non-coding genes and overlapping genes on the same strand were removed from this analysis. PAL-seq tags from all wild-type samples were combined, filtered to remove tags with a poly(A) tail <20 nt to ensure they represent genuine cleavage and polyadenylation sites, and then aligned with STAR. 3' UTRs were annotated by identifying positions where ≥2 PAL-seq tags aligned at the same position between the stop codon of one gene and the transcription start site of the neighboring gene on the same strand. Isoforms were annotated as described previously, using a 60 nt window and a maximum 3'-UTR length of 4100 nt (so as to not exceed the 99th percentile of previously annotated fly 3' UTRs). If a new 3' terminus was annotated that was distal to the initial model, the 3' end of the model was extended to the distal isoform. If not, the initial representative transcript model was used. The resulting set of representative transcript models (referred to as 'mRNAs') used for our analyses is available at http://www.ncbi.nlm.nih.gov/geo under accession number GSE83616. Models for which the selected 3'-end isoform (either the initial refFlat 3' end or an extension based on our re-annotation) overlapped with a snoRNA or snRNA were included in the analysis of sequencing data but excluded from figures.

50Samples were split and prepared for ribosome-footprint profiling and RNA-seq as described, as well as for PAL-seq as described. Detailed protocols are available at http://bartellab.wi.mit.edu/protocols.html. The wild-type stage 14 oocyte and activated egg samples were prepared using RPF lysis buffer and sucrose gradients that did not contain cycloheximide, all other samples were prepared using RPF lysis buffer and sucrose gradients that contained cycloheximide. Because TE changes between wild-type stage 14 oocytes and activated eggs for samples prepared with and without cycloheximide correlated very well, the omission of cycloheximide did not affect our observations pertaining to these stages. RPF and RNA-seq reads were mapped to ORFs of representative transcript models as described, which excluded the first 50 nt of each ORF to eliminate signal from ribosomes that initiated after adding cycloheximide. PAL-seq tags were mapped to 3' UTRs of representative transcript models, and PAL-seq fluorescence intensities were converted to tail lengths as described.

51TE values (log2) for stage 11, 12, 13, and 14 oocytes, or for stage 14 oocytes, 0–1 hr embryos, and 2–3 hr embryos were normalized by the median TE value for the mRNAs in the corresponding sample from Figure 2 or Figure 1B, respectively, and were then entered into the Short Time-series Expression Miner (STEM), which normalized the data relative to the first sample and then used the STEM clustering method to identify significantly enriched profiles (p value < 0.05, following Bonferonni correction). The default settings were used, except a maximum unit change of one in model profiles was allowed between time points, and only mRNAs with ≥0.5 or ≤–0.5 log2 fold change between any two samples were placed in a cluster. A heatmap of all expression patterns was generated in Matlab using the normalized data.

52The RNA-seq and ribosome-footprint profiling data for wild-type stage 14 oocytes and activated eggs, and png-mutant activated eggs were previously published and were reanalyzed for this study. The raw data for these previously published datasets are available at the Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo) under accession number GSE52799, and the processed data files used for this study are under accession number GSE83616. The PAL-seq data for those samples had not been published and are available along with all of the other RNA-seq, ribosome profiling, and PAL-seq data and processed data files under accession number GSE83616. All of the processed data are also available in Supplementary file 2.

53All correlations report the Spearman correlation coefficient, and all statistical tests were two-sided unless noted otherwise. The TE changes observed between biological replicates for wild-type stage 14 oocytes and activated eggs correlated very well (Rs = 0.94), even though cycloheximide was omitted from one of the replicates. Likewise, TE of activated egg highly correlated with that of 0–1 hr embryo (Rs = 0.88). Although we did not analyze replicates of PAL-seq measurements in the current study, previous analysis of biological replicates demonstrates the reproducibility of mean tail lengths measured using this technique (Rs = 0.83). Similar correlations were observed in the current study when comparing stages for which little if any differences in tail lengths were expected, including activated egg and 0–1 hr embryo (Rs = 0.86), stage 11 and stage 12 oocyte (Rs = 0.94), stage 12 and stage 13 oocyte (Rs = 0.88), as well as png and wild-type oocytes at stage 14 (Rs = 0.86), which was prior to PNG activation (Supplementary file 1 and Figure 6—figure supplement 1A).