We used recombineering (Warming et al., 2005) to modify a bacterial artificial chromosome (BAC) (Venken et al., 2006) containing the D. melanogaster eve gene and all of its enhancers and regulatory elements (Venken et al., 2009), replacing the coding region with an array of MS2 stem loops followed by the D. melanogaster yellow (y) gene (Figure 1D; Bothma et al., 2014). The 4329 base pair y gene, which is longer than the endogenous eve transcript, is phenotypically neutral and provides a means to increase the number of RNA Polymerase II (Pol II) molecules loaded onto the gene in order to amplify the signal (see Materials and methods for a discussion of how the structure of the reporter genes affects the fluorescence signal, analyses and inferences performed throughout this work). We inserted the engineered BAC into a targeted site on chromosome 3L using ΦC31 integrase-mediated recombination (Fish et al., 2007), and homozygosed the line, which showed no signs of adverse effects of the transgene.

We crossed males from this line with females carrying transgenes that express in embryos an RFP-labeled histone to visualize nuclei, and an MCP-GFP fusion that binds the MS2 RNA stem loops. The result is the direct visualization and quantification of nascent transcripts at the transgene locus as fluorescent puncta (Garcia et al., 2013). The temporal and spatial pattern of eve transgene transcription recapitulates the well-characterized dynamics of eve expression, most notably formation of the characteristic seven stripes in the late blastoderm (Figure 2; Video 1; Surkova et al., 2008; Fowlkes et al., 2008; Jiang et al., 2015; Ludwig et al., 2011). Further, as recently demonstrated in Lammers et al., 2020 this BAC reporter construct quantitatively recapitulates the cytoplasmic eve mRNA pattern as measured by FISH (Lammers et al., 2020; Luengo Hendriks et al., 2006; Lim et al., 2018).

We used laser-scanning confocal microscopy to record, with high temporal resolution and high magnification, two color (MCP-GFP and histone-RFP) movies of embryos from before nc14 through gastrulation. We optimized our data collection strategy to sample multiple stripes (3 to 5) in each movie, to obtain high temporal resolution (one Z-stack, corresponding to each time point of our movies, every 16.8 s) and to have optimal signal-to-noise with minimal bleaching. In total, we collected 11 movies (Videos 2–12), with every stripe imaged at least five times (see Table 1).

We used a custom image processing pipeline (Garcia et al., 2013; Lammers et al., 2020) to identify nuclei, track fluorescent puncta and extract fluorescence intensities in each nucleus over time. The resulting data (Supplementary file 1) contains fluorescence traces from 2961 nuclei at an interpolated time interval of 20 s, representative examples of which are shown in Figure 3A.

We first sought to reexamine the previously characterized temporal dynamics of stripe formation (Surkova et al., 2008; Fowlkes et al., 2008; Jiang et al., 2015; Ludwig et al., 2011) using the increased temporal resolution (relative to earlier analyses of fixed embryos and of slowly maturing fluorescent protein reporters) of these data (Figure 3B). Early imaging studies described eve as being expressed broadly in nc13 and early nc14 embryos before refining sequentially into four, then seven stripes (Macdonald et al., 1986; Frasch and Levine, 1987). Subsequent work with improved labeling and imaging techniques (Surkova et al., 2008; Fowlkes et al., 2008) revealed an initial phase with broad domains in the anterior and posterior, followed by the formation of stripes from within these broad domains and, eventually, amplification of the stripe pattern.

During nc14, we first observe eve transcription beginning approximately five minutes after the onset of anaphase. The initial transcription covers a broad swath of the embryo, from around 25 to 75% egg-length, with the highest activity in two domains roughly centered in the anterior and posterior halves of the embryo, respectively. The greatest fluorescence signal during the first 25 min of nc14, when stripes are not yet fully formed, is in the most anterior region of eve transcription, in an area in which stripe one will eventually form.

Although the full seven stripe pattern is not fully formed until around 25 min, the three anterior-most stripes are already apparent as locally high areas of fluorescence intensity as early as 10 min. By 20 min stripes 1, 2, and 3 have clearly separated from background, stripes 4 and 6 appear to split off from a large posterior domain, and stripe 7 forms de novo in the posterior. Stripe 5 appears as a distinct stripe last, emerging in an area of low transcriptional activity left behind following the splitting of stripes 4 and 6. The stripes persist for the remainder of nc14, gradually increasing in fluorescence intensity until they reach a peak at around 35 min into nc14.

The positions of stripes 1–3 along the anterior-posterior (AP) axis are largely stable after they form, while stripes 4–6 show small anterior shifts. Stripe 7 makes a more dramatic movement towards the anterior, moving approximately 8% of egg-length, or around 40 μm from its initial location. The quantitative characterization of this stripe movement, the decoupling between stripes and nuclei, and the quantification of transcriptional bursting dynamics in each nucleus necessitated the development of a method, described below, to dynamically define the position of stripes throughout each movie.

As expected, the fluorescence traces from individual nuclei show clear hallmarks of transcriptional bursting, with apparent stochastic fluctuations between states with low and high fluorescence output (Figure 3A). Following previous work in the field (Golding et al., 2005; Chubb et al., 2006; Zenklusen et al., 2008; Lionnet et al., 2011; Muramoto et al., 2012; Little et al., 2013; Xu et al., 2015; Lenstra et al., 2015; Fukaya et al., 2016; Desponds et al., 2016; Hendy et al., 2017; Zoller et al., 2018; Bothma et al., 2014; Paré et al., 2009; Lim et al., 2018), we model bursting as a simple Markovian process in which a promoter switches stochastically between an OFF and an ON state with rates k_on and k_off. When the promoter is in the ON state, we assume it loads polymerases continuously with a constant rate r (Figure 1A).

In our implementation of the MS2 system, once a polymerase molecule transcribes the stem loops located at the 5’ end of the gene, the MCP-GFP molecules bound to the stem loops produce a constant fluorescent signal at the locus that persists until this polymerase completes its traversal of the gene. Building off of the method presented in Lammers et al., 2020, we estimated this polymerase transit time as the displacement that gives the minimum value in the autocorrelation of the single frame differences in the fluorescent signal (see Materials and methods). The rationale for this approach was that every increase in signal due to polymerase loading at time t should be accompanied by a corresponding decrease in signal at time t+t_elong due to the completion of a transcriptional elongation cycle with a delay equal to the elongation time (Coulon and Larson, 2016; Desponds et al., 2016). We arrived at an estimate of 140 s (Figure 4A), consistent with a direct measurement of the rate of polymerase elongation of ~2700 bp/min from Fukaya et al., 2017 and the length of the construct (6,563 bp).

We model the bursting process at each promoter in discrete time steps of Δt = 20 s, set by the time resolution of our imaging. Under our model, in each time window a promoter is either OFF and not loading polymerases, or ON and loading polymerases at a fixed rate. A promoter that is in the ON state loads Δt x r polymerases, producing a single pulse of fluorescence proportional to Δt x r (with the proportionality factor determined by the fluorescence of GFP and the fraction of MS2 loops bound by MCP-GFP). This pulse lasts at the locus for 140 s, at which point all polymerase molecules loaded during the original time window have terminated transcribing (Figure 4B). Since we do not calibrate the fluorescence signal to the number of polymerase molecules for this construct, in practice we fold the proportionality factor into r altering its units from polymerases loaded per unit time to fluorescence signal produced per unit time. Since many transcriptional bursts last for longer than 20 s, the fluorescence output of a single burst is a sum of the pulses generated during each time window.

In the embryos we imaged here, the MS2 BAC is heterozygous, contributed only by the father, while the mother contributes the MCP-GFP. However, DNA replication occurs within an average of 10 min for loci in nc14 (McKnight and Miller, 1977; Rabinowitz, 1941; Shermoen et al., 2010), meaning that there are actually two sister chromatids with the MS2-containing transgene in every nucleus. Because of sister chromatid cohesion, we cannot, in general, discriminate both copies (Fung et al., 1998; Wilkie et al., 1999; Little et al., 2013). Thus we model the locus as having two promoters, which implies that the system can be in one of three distinct states: OFF, one promoter ON, and two promoters ON. For ease of exposition, however, we will frame our discussion in terms of an effective 2-state model in which the locus is ON so long as at least one promoter is in the ON state. As it is still unclear how the sister chromatids influence each other’s transcription (McKnight and Miller, 1977; Lammers et al., 2020; Zoller et al., 2018), we make no prior assumptions about the nature or degree of correlations between sister chromatids (though, for illustration purposes, in Figure 4B we entertain the case where sister chromatids transcribe independently).

As Illustrated in Figure 4C for the case of independent sister chromatids, if we know the state of the promoter over time, we can reconstruct its expected fluorescence output by summing 140 s pulses beginning at each point where the promoter is ON and having height r if one promoter is ON or height 2 x r if two promoters are ON. Traces modeled from hypothetical promoter state sequences (Figure 4C) have the features of the observed fluorescence signal: linear increases in intensity (corresponding to periods when the promoter is ON); plateaus (corresponding to periods when transcriptional initiation is matched with previously initiated polymerases completing their transit of the gene); and linear signal decays (corresponding to periods when the promoter is OFF but previously initiated polymerases are still transiting the gene) (Bothma et al., 2014; Garcia et al., 2013).

However, when given a fluorescence trace, it is not trivial to infer the promoter state sequence that generated it, owing to the time convolution between promoter state and fluorescence output. To solve this problem, we developed a compound-state hidden Markov Model (cpHMM, described in Lammers et al., 2020) that estimates global parameters k_on, k_off, and r for a set of traces, and allows us to identify the maximum-likelihood promoter state sequence under these parameters for every trace via the Viterbi algorithm.

The cpHMM thus accomplishes two aims central to treating these data in a more rigorous and biologically meaningful manner. First, it allows us to describe the bursting behavior of any set of nuclei in quantitative terms. Across all seven stripes, the model infers approximate k_on k_off values of 0.60 events per minute and an r of 67 AU per minute. And second, by providing a means to fit a sequence of ON and OFF states to the data from each nucleus, the cpHMM allows us to shift the focus in the analysis of individual traces from fluorescence, which only indirectly reflects the temporal behavior of a promoter, to the instantaneous promoter state (Figure 4D–F; see also Supplementary file 1 which the inferred promoter state for each nucleus at every time point and the corresponding modeled fluorescence intensity, and Videos 13–23).

Before analyzing the data further we had to solve two practical problems. To compare the kinetic behavior of individual stripes, we had to determine which nuclei were in each stripe at every time point, a process complicated by the movement of stripes relative to both the embryo and nuclei. Further, to analyze the data in toto, we also had to register the 11 movies relative to each other and to the embryo.

To address these challenges, we used a Gaussian mixture model to cluster bursting nuclei in each movie in a series of overlapping six-minute time windows based on their x and y positions in the image (Figure 5A). This clustering reliably separates nuclei into individual stripes. We next determined the orientation of each stripe to the AP and imaging axes by fitting a line to coordinates of all nuclei assigned to that stripe in each movie (Figure 5B). We fit a line with this slope to bursting nuclei from each time window (Figure 5C and D), and use these fits to generate a linear model of the position of each stripe in each image over time, which we use to reorient the stripe so that it is perpendicular to the image x-axis (Figure 5E).

We next use the known coordinates of the anterior and posterior poles of the embryo to convert the image x-axis to AP position, and register the examples of each stripe from different movies by setting the AP position of the center of each stripe at 35 min in nc14 to the mean AP position of all examples of that stripe at 35 min, adjusting the position of the stripe at other time points by the same correction (Figure 5F). As the stripes are not all present until after 25 min in nc14, we assign and register nuclei before that point based on the stripe position at 25 min. The stripe assignment is invariant over bootstrapping of movies, and the standard deviation of the AP displacement over bootstrapping of movies is 0.16% of embryo length.

Collectively these data represent an easy way to visualize and interpret kinetic fingerprint of stripe formation: a record of every transcriptional burst that occurred during the formation of eve stripes in these embryos (Figure 6; Video 24).

We used the output of the cpHMM and registration process to examine the locations of transcriptional bursts along the AP axis and over time (Figure 6). The most striking feature is the almost complete lack of observable transcriptional bursts in the regions between stripes from 25 min into nc14, with the exception of the 5–6 interstripe which is discussed below (note that this is not an artifact of the movie alignment and orientation process, as this effect is seen clearly in individual movies). We took advantage of the fact that we were tracking bursts in individual nuclei in order to analyze the relationship between this absence of bursting in interstripe regions and the single-nucleus bursting behavior within stripes.

Stripes are defined by sharp spatial boundaries, with the transition between the low-bursting (quiescent) state and the frequently bursting (active) state occurring from one column of nuclei to the next (Figure 6), consistent with the classical descriptions of eve stripe patterns (Small et al., 1991; Frasch and Levine, 1987; Fujioka et al., 1999; Small et al., 1992; Warrior and Levine, 1990; Clyde et al., 2003). They also have sharp temporal boundaries: all of the interstripe regions, save that between stripes 6 and 7, form in regions where there was appreciable bursting early in nc14 that disappears at around 25 min into the nuclear cycle (Figure 6).

To better understand how the low-bursting state in interstripes is established, we looked at the bursting history of the nuclei in these regions (Figure 7). The first feature we noticed was that most of the nuclei that ultimately form the interstripe were never detected to burst at any point in nc14 (Figure 7A,B). With the exception of the 5–6 interstripe, these never-ON nuclei effectively form the boundaries between stripes, as essentially every nucleus within each stripe bursts at some point during nc14 (Figure 7A,B).

The contrast in bursting history between stripes and interstripes is less pronounced in the posterior, where there are fewer such never-ON nuclei in the interstripe region (Figure 7B, notice the lower density of red single-nuclei tracks corresponding to never-ON nuclei). In order to reveal the source of this reduced number of never-ON nuclei in posterior interstripes, we analyzed their bursting history. Figure 7C shows the AP positions of the nuclei in one movie covering stripe seven as a function of time, with the period in which they are part of the stripe highlighted. Although the stripe is clearly present throughout this period, no nuclei remain a part of the stripe for the entirety of this 25 min period. As time progresses, nuclei at the posterior edge of stripe seven shift from an active state, in which the promoter stochastically alternates between the ON and OFF transcriptional states, to a quiescent state in which we observe no appreciable bursting. In contrast, nuclei just off the anterior edge of the stripe switch from a quiescent to an active state at roughly the same rate. This leads to a net overall anterior movement of the stripe, akin to treadmilling, at a velocity of approximately one percent of embryo length every three minutes.

Consistent with Lim et al., 2018, the other stripes exhibit smaller and varied anterior shifts (Figure 7—figure supplement 1), but in every case the shift is associated with a similar coupled gain of active nuclei along the anterior edge and loss along the posterior edge. This effect is most clearly seen in Figure 7D, which shows, for each time point where a nucleus initiates a new burst, the difference in activity (defined as the difference between the fraction of the time the nucleus is in the ON state in the subsequent 10 min minus the fraction of the time the nucleus is in the ON state in the preceding 10 min). For all seven stripes there is a clear spatial pattern, with nuclei along the anterior edge of the stripe entering a bursting state and nuclei along the posterior edge becoming quiescent, indicating a movement of stripes relative to nuclei. Hence, stripe movement is associated with the dynamic switching of nuclei between active and quiescent states, and not just with the movement of nuclei themselves.

We next turned to the questions of how the spatial pattern of nuclear transcriptional activity described above is produced by regulating bursting kinetics, and whether this regulation differs among the seven eve stripes. In principle, any pattern of transcriptional activity could be achieved by modulating the duration, separation and/or amplitude of bursts across space and time. For example, a stripe could be created by varying burst separation along the anterior-posterior axis, with nuclei in the stripe center having lower burst separation, and those outside the stripe having long periods without bursts, or no bursts at all. Alternatively, the same stripe could be created with uniform burst separation across nuclei, but elevated burst duration or amplitude within the stripe, or by modulating multiple parameters simultaneously.

Ideally, we would like to have a measure of the bursting parameters governing the behavior of every nucleus. However, individual MS2 traces have too few time points to allow for accurate cpHMM inference of burst parameters at the single trace level. We therefore used the cpHMM to infer k_on, k_off, and r for groups of nuclei binned on their mean fluorescence output and stripe. The logic of the fluorescence binning was that, given that ⟨ fluorescence ⟩∝rkonkon+koff (Lammers et al., 2020), nuclei that have similar k_on, k_off, and r will have similar fluorescence outputs. Our inference shows that k_on is very strongly regulated as a function of average fluorescence in a consistent manner across stripes (Figure 8A). In contrast, only a weak drop in k_off is observed (Figure 8B). Finally, r also featured a strong upregulation as a function of average fluorescence across stripes (Figure 8C).

As shown in Figure 8D, each stripe contains nuclei with a relatively wide range of average fluorescence values. In order to reveal the bursting parameters across the AP axis for each stripe, we averaged the single-cell bursting parameters determined in each stripe (Figure 8A–C) weighted by the relative number of nuclei in each fluorescence bin present at each position along the AP axis (Figure 8D). We find that the variation in bursting parameters observed as a function of average fluorescence largely echoes the modulation of fluorescence in space (Figure 8E). Specifically, while there is a subtle downregulation of k_off within stripes, k_on and r are significantly upregulated in the center of each stripe.

Thus, not only do the five eve enhancers employ a common regulatory strategy for modulating the fluorescence output of nuclei to create a stripe, decreasing burst separation and increasing burst amplitude with a constant burst duration, the precise quantitative relationship among these bursting parameters is maintained across a wide range of molecular inputs and fluorescence outputs.

Just as gene expression patterns are dynamic in time (Bothma et al., 2014), they are dynamic in space, resulting in the movement of expression domains throughout the embryo during development (Jaeger et al., 2004b; Keränen et al., 2006). The anterior movement of eve stripes during nc14 has been previously described (Keränen et al., 2006; Lim et al., 2018), and proposed to arise from a combination of nuclear movement (nuclear flow) and movement in the pattern of regulators (expression flow), especially repressors, which are known to shift anteriorly during nc14 as well (Jaeger et al., 2004b; Jaeger et al., 2004a). While Keränen et al., 2006 concluded that the relative contributions of these two forces were roughly equal, our data suggest that, especially in the posterior, expression flow dominates the anterior shift of the eve stripes.

A typical nucleus in stripe seven moves around one percent of embryo length in the final 25 min of nc14. The stripe, however, moves around five percent of embryo length during that time (see Figure 7C). Because we are tracking both the position and activity of individual nuclei, we can visualize expression flow in action. We see nuclei transition from low activity in the anterior interstripe to high activity in the stripe, from high activity in the stripe to low activity off the posterior flank of eve expression, and in some cases both.

This effect is most pronounced for the posterior stripes, but is observed for the more anterior stripes as well, although the magnitude of the shift decreases for more anterior stripes (Figure 7—figure supplement 1). The difference in the amount of the effect we and Keränen et al. attribute to expression flow is likely an effect of differences in the data used. Because we are looking at instantaneous transcription rates while they looked at accumulated mRNA, there is a considerable temporal lag and integration of the transcriptional activity over the life time of eve mRNA in their data, which has the effect of underestimating the extent to which the stripes actually move.

We also note that the extent to which nuclear flow by itself would be expected to shift output patterns measured as a function of position in the embryo is unclear, as it would depend on the extent to which the repositioning of regulators drives movement of nuclei (which it is believed to do [Blankenship and Wieschaus, 2001]), and the corresponding effect that nuclear movement has on the positioning of regulators, which remains largely unknown.

One open question relates to the temporal relationship between changes in the position of the repressor array that drives stripe position and the transcriptional output of the stripes. For example, the anterior shift of the stripes of eve as well as fushi tarazu has been proposed to originate, in part, from cross-repression between these two genes (Lim et al., 2018). Recent advances in the simultaneous monitoring of protein concentration and transcriptional output in living embryos should help answer this question in the near future (Bothma et al., 2018; Lammers et al., 2020).

That gene expression is a fundamentally dynamic process is not new information. However, the tools we have had at our disposal to study gene expression so far have tended to emphasize its static features, down to the language we use to describe the transcriptional output of a gene. In textbooks and the scientific literature, eve has a gene expression pattern consisting of seven stripes. But, as some earlier work emphasized (Janssens et al., 2006), and we have directly visualized here, the transcriptional output of eve, rather than a single ‘pattern’ is a rapidly changing function of time and space: it is dynamic at many time scales and across space and nuclear positions. Indeed, at no point does eve approach anything even remotely like a steady state.

We are at the dawn of a new period in the study of transcription, as new experimental techniques and advanced microscopy allow us to monitor transcriptional regulators, observe their behavior at the single-molecule level, and track the transcriptional output of a gene in living, developing animals. We have only barely begun to understand this new data and what it can tell us about biology. While the focus in this paper was on a single gene in a single species, we hope that this and our accompanying work (Lammers et al., 2020) will have a broader impact by beginning to establish rigorous frameworks for quantifying, characterizing and visualizing the dynamics of transcription at the single-cell level during development that will be required in the era of live imaging of transcription in development.

We used bacterial recombineering (Warming et al., 2005) to modify a bacterial artificial chromosome (BAC) (Venken et al., 2006) containing the D. melanogaster eve gene and all of its enhancers and regulatory elements (BAC CH322-103K22) (Venken et al., 2009). We replaced the coding region with an array of 24 MS2 stem loops fused to the D. melanogaster yellow gene (Figure 1B; Bothma et al., 2014) as described below. We inserted our eve::MS2::yellow BAC-based construct in the D. melanogaster genome at chromosome 3L through ΦC31 integrase-mediated recombination (see Generation of fly lines), and generated a viable homozygous fly line (w-; +; eve::MS2::yellow) as detailed below.

In principle the length of the reporter should not limit our ability to estimate burst parameters. However, in practice a reporter construct that is too short will have insufficient signal. Further, one that is too long will increase the dwell time of each RNA polymerase molecule on the gene and, as a result, our cpHMM inference will require too many computational resources. Our choice of reporter construct structure strikes a balance between these two limitations and is ideally suited for inferring bursting parameters in the time range where eve resides, as well as for boosting the signal-to-noise ratio. See Lammers et al., 2020 for a more detailed discussion of reporter length-related tradeoffs.

We modified a CHORI BAC CH322-103K22 derived from Venken et al., 2009, which contained the entire eve locus and a GFP reporter instead of the eve coding sequence (CH322-103K22-GFP). We replaced the GFP reporter with MS2::yellow (total insert size, 6665 bp) through a two step, scarless, galK cassette-mediated bacterial recombineering (Warming et al., 2005). Briefly, we transformed our starting CH322-103K22-GFP BAC into E. coli recombineering strain SW102. We then electroporated the strain with a galK cassette flanked by 50bp-long DNA homology arms homologous to the MS2::yellow (6665 bp) reporter. Upon electroporation, we selected transformants on M63 minimal media plates with galactose as a single carbon source. We achieved a correct replacement of GFP sequence by galK cassette in the BAC context (CH322-103K22-galK), validated by observing the digestion patterns produced by ApaLI restriction enzyme.

We next purified the CH322-103K22-galK BAC and transformed it into fresh E. coli SW102 cells. We electroporated these cells with the purified MS2::yellow insert and used M63 minimal media plates with 2-deoxy-galactose to select against bacteria with a functional galK gene. We used colony PCR to screen for colonies with a correct MS2::yellow insertion (CH322-103K22-MS2) replacing the galK cassette. We validated this insertion by observing ApaLI, XhoI, SmaI, and EcoRI restriction digestion patterns and through PCR and Sanger sequencing of the insertion junctions. We transformed our CH322-103K22-MS2 BAC in E.coli EPI300 cells to induce high copy numbers and purified it with a Qiagen plasmid Midiprep kit.

We sent a sample of our purified CH322-103K22-MS2 BAC to Rainbow Transgenic Flies, Inc for injection in D. melanogaster embryos bearing a ΦC31 AttP insertion site in chromosome 3L (Bloomington stock #24871; landing site VK00033; cytological location 65B2). We received the flies that resulted from that injection and used a balancer fly line (w-; +; +/TM3sb) to obtain a viable MS2 homozygous line (w-; +; MS2::yellow). We used line (yw; His::RFP; MCP::GFP) as the maternal source of Histone-RFP and MCP-GFP (Garcia et al., 2013).

Embryo collection and mounting was done as specified in Garcia and Gregor, 2018. In short, we set fly crosses between ~30 males (w-; +; eve::MS2::yellow) and ~80 females (yw; His::RFP; MCP::GFP) in a plastic cage capped with a grape juice agar plate. We collected embryos from cages two to ten days old by adding a fresh plate for 30 min and aging for 60 min to target embryos 90 min or younger.

Embryos were mounted on a gas-permeable Lumox Film (Sarstedt - Catalog # 94.6077.317) embedded on a microscope slide hollowed on the center. Then, we coated the hydrophobic side of the Lumox film with heptane glue and let it dry. The film allows for the oxygenation of embryos during the 2–3 hr long imaging sessions while heptane immobilizes them.

We soaked an agar plate with Halocarbon 27 oil, picked embryos with forceps, and laid them down on a 3 × 3 cm piece of paper tissue. We dechorionated embryos by adding 2 drops of bleach diluted in water (5.25%) on the paper tissue and incubating for 1.5 min. We removed bleach with a clean tissue and rinsed with ~4 drops of distilled water. We then placed the tissue paper with dechorionated embryos in water, and picked buoyant embryos with a brush.

We lined ~ 30 apparently healthy embryos on the Lumox film slide and added 2–3 drops of Halocarbon 27 oil to avoid desiccation, and covered the embryos with a cover slip (Corning Cover Glass, No. 1, 18 × 18 mm) for live imaging.

Movies of embryonic development were recorded on a Zeiss-800 confocal laser-scanning microscope in two channels, (EGFP: 488 nm; TagRFP: 561 nm). We imaged embryos on a wide field of view, along their anterior-posterior axis, of 1024 × 256 pixels (202.8 x 50.7 µm), encompassing 3–5 stripes per movie. We tuned laser power, scanning parameters, master gain, pinhole size and laser power to optimize signal-to-noise ratio without significant photobleaching and phototoxicity.

For imaging, the following microscope settings were used: 63x oil-objective, scan mode ‘frame’, pixel size of 0.2 µm, 16 bits per pixel, bidirectional scanning at a speed of 7, line step of 1, laser scanner dwelling per pixel of 1.03 µs, laser scanner averaging of 2, averaging method Mean, averaging mode Line, 488 nm laser power of 30 µW (EGFP), 561 nm laser power of 7.5µW (TagRFP) (both powers were measured with a 10x air-objective), Master Gain in EGFP detector of 550V, Master Gain in TagRFP detector of 650V, Digital Offset in both detectors of 0, Digital Gain in both detectors of 1.0, and a pinhole size of 1 airy unit under the imaging conditions mentioned above (44, 0.7 µm/section), laser filters EGFP:SP545 and TagRFP:LBF640. This resulted in an imaging time of 633 ms per frame and a full Z-stack of 21 frames in intervals of 0.5 µm every 16.8 s. Following (Bothma et al., 2014; Bothma et al., 2015; Bothma et al., 2018; Lammers et al., 2020), the imaging conditions were determined not to affect normal development as reported by the timing of the nuclear cycles in early development. We stopped imaging after 50 min into nuclear cycle 14, and took mid-sagittal and surface pictures of the whole embryo for localization of the recorded field of view along the embryo’s AP axis.

We used a Matlab computational pipeline based on Garcia et al., 2013; Lammers et al., 2020 to segment and extract numeric data from our raw movies. Briefly, this software segments and processes the images from the two channels (channel 1: MCP::GFP, channel 2: Histone::RFP) on which we collected our data. For segmentation of channel 1, we used Fiji-Weka Segmentation 3D software; this machine-learning-based method relies on the manual segmentation of a variety of MCP::GFP labeled transcriptional foci in a given 21 frame Z-stack from a single dataset (EVE_D11) to produce a model for the segmentation of all datasets recorded under the same imaging conditions. Next, we segmented and tracked the Histone::RFP-labeled nuclei on channel 2. Subsequently, we assigned MCP::GFP labeled transcriptional foci to their corresponding Histone::RFP-labeled nuclei. Since we collected whole embryo pictures of each of our datasets, we were able to match and locate the recorded fields of view to their right position along the body of their corresponding embryos. Finally, we extracted position and fluorescence values over time of all transcriptional foci to generate data structures ready to use in further analyses.

To estimate the transit time of the polymerase along the construct (which is used to determine the persistence of the fluorescence signal from a single transcript at the locus) we first calculated, for each nucleus, the difference in fluorescence signal between adjacent timepoints D_n,t = F_n,t+1 - F_n,t where F_n,t is the fluorescence signal for nucleus n at time point t and then calculated the Pearson correlation coefficient of the vectors [..., D_n,t, D_n,t+1, D_n,t+2,. ..] and [..., D_n,t+d, D_n,_t+d+1, D_n,t+d+2,. ..] over values of d from 1 to 20 representing time displacements of 20 to 400 s. The minimum correlation occurred at 140 s.

For this work we employed a statistical method that utilizes a compound-state hidden Markov Model to infer bursting parameters from experimental fluorescence traces. The theory and implementation of this method are described in detail in Lammers et al., 2020. Briefly, parameters were inferred using a standard version of the Expectation Maximization Algorithm implemented using custom-written scripts in Matlab. Our inference is carried over the full duration of activity of each active nucleus during nuclear cycle 14. Bootstrap sampling was used to estimate the standard error in our parameter estimates. Subsets of 3000 data points were used to generate time-averaged parameter estimates. Inference was not conducted for groups for which fewer than 1000 time points were available.

All data were analyzed in Python using a Jupyter notebook with custom code to process raw data and generate figures. The Jupyter notebook and all data required to run it is available in Supplementary file 1 and at (https://github.com/mbeisen/Berrocal_2020; Berrocal, 2020; copy archived at swh:1:rev:d983098bd5183f9907d633c425f80b2cb5282a8b).

We first filtered the raw data to remove data with observations spanning less than 2,000 s, as well as nuclei that were poorly tracked over time defined as nuclei that moved across the movies at an average rate of over 4.5 pixels per minute. This left 430,073 observations from 2959 nuclei.

We used the Gaussian mixture model module of the Python library scikit-learn (Pedregosa et al., 2011) to cluster all nuclei time points in each movie in each of a series of overlapping 428 s time windows beginning at 25 min in nc14, specifying the number of components equal to the number of stripes captured in the movie and using the setting covariance_type='tied'. We preliminarily assigned nuclei time points to a stripe if they were consistently clustered in that stripe in the relevant time windows. We then pooled all nuclei time points assigned to the same stripe and fit a line to the median x and y positions in the bottom (y < 128) and top (y > 128) halves of the image. We considered the slope of this line to represent the orientation of the stripe to the image x-axis. We then went back to each time window and fit the nuclei assigned to the stripe with a line with the previously computed slope fixed. This produced an association of time with stripe position, from which we derived a linear model that describes the position of each stripe in each movie at every time point.

We assigned all nuclei time points (not just bursting ones) to stripes by identifying the stripe whose predicted position at the relevant time was closest (along the x-axis) to the nucleus being analyzed, and assigned a nucleus to the most common stripe assignment for its individual time points. We then corrected the reorientation of the stripe at each time point to be perpendicular to the image x-axis (to enable projection along the AP axis) by setting its new image x-axis position to be the x position of the stripe in the middle of the y-axis (y = 128) plus the offset of the nucleus to the unoriented stripe along the x-axis. Finally, we used the positions of the anterior and posterior poles of the embryo to map image x coordinates to AP position. We then adjusted the AP position of each stripe in each movie such that the center of the stripe at 35 min in nc14 had the same AP position.