Chromatin Architecture Emerges during Zygotic Genome Activation Independent of Transcription.

Chromatin architecture is fundamental in regulating gene expression. To investigate when spatial genome organization is first established during development, we examined chromatin conformation during Drosophila embryogenesis and observed the emergence of chromatin architecture within a tight time window that coincides with the onset of transcription activation in the zygote. Prior to zygotic genome activation, the genome is mostly unstructured. Early expressed genes serve as nucleation sites for topologically associating domain (TAD) boundaries. Activation of gene expression coincides with the establishment of TADs throughout the genome and co-localization of housekeeping gene clusters, which remain stable in subsequent stages of development. However, the appearance of TAD boundaries is independent of transcription and requires the transcription factor Zelda for locus-specific TAD boundary insulation. These results offer insight into when spatial organization of the genome emerges and identify a key factor that helps trigger this architecture.

Continued Activity of the Pioneer Factor Zelda Is Required to Drive Zygotic Genome Activation.

Reprogramming cell fate during the first stages of embryogenesis requires that transcriptional activators gain access to the genome and remodel the zygotic transcriptome. Nonetheless, it is not clear whether the continued activity of these pioneering factors is required throughout zygotic genome activation or whether they are only required early to establish cis-regulatory regions. To address this question, we developed an optogenetic strategy to rapidly and reversibly inactivate the master regulator of genome activation in Drosophila, Zelda. Using this strategy, we demonstrate that continued Zelda activity is required throughout genome activation. We show that Zelda binds DNA in the context of nucleosomes and suggest that this allows Zelda to occupy the genome despite the rapid division cycles in the early embryo. These data identify a powerful strategy to inactivate transcription factor function during development and suggest that reprogramming in the embryo may require specific, continuous pioneering functions to activate the genome.

Multi-level regulation of even-skipped stripes by the ubiquitous factor Zelda.

The zinc-finger protein Zelda (Zld) is a key activator of zygotic transcription in early Drosophila embryos. Here, we study Zld-dependent regulation of the seven-striped pattern of the pair-rule gene even-skipped (eve). Individual stripes are regulated by discrete enhancers that respond to broadly distributed activators; stripe boundaries are formed by localized repressors encoded by the gap genes. The strongest effects of Zld are on stripes 2, 3 and 7, which are regulated by two enhancers in a 3.8 kb genomic fragment that includes the eve basal promoter. We show that Zld facilitates binding of the activator Bicoid and the gap repressors to this fragment, consistent with its proposed role as a pioneer protein. To test whether the effects of Zld are direct, we mutated all canonical Zld sites in the 3.8 kb fragment, which reduced expression but failed to phenocopy the abolishment of stripes caused by removing Zld in trans. We show that Zld also indirectly regulates the eve stripes by establishing specific gap gene expression boundaries, which provides the embryonic spacing required for proper stripe activation.

Dense Bicoid hubs accentuate binding along the morphogen gradient.

Morphogen gradients direct the spatial patterning of developing embryos; however, the mechanisms by which these gradients are interpreted remain elusive. Here we used lattice light-sheet microscopy to perform in vivo single-molecule imaging in early Drosophila melanogaster embryos of the transcription factor Bicoid that forms a gradient and initiates patterning along the anteroposterior axis. In contrast to canonical models, we observed that Bicoid binds to DNA with a rapid off rate throughout the embryo such that its average occupancy at target loci is on-rate-dependent. We further observed Bicoid forming transient "hubs" of locally high density that facilitate binding as factor levels drop, including in the posterior, where we observed Bicoid binding despite vanishingly low protein levels. We propose that localized modulation of transcription factor on rates via clustering provides a general mechanism to facilitate binding to low-affinity targets and that this may be a prevalent feature of other developmental transcription factors.

Transcriptome analysis reveals the role of Zelda in the regulation of embryonic and wing development of Tribolium castaneum.

Zinc finger protein (Zelda) of Tribolium castaneum (TcZelda) has been showed to play pivotal roles in embryonic development and metamorphosis. However, the regulatory mechanism of TcZelda associated with these physiology processes is unclear. Herein, the developmental expression profile showed that Zelda of T. castaneum was highly expressed in early eggs. Tissue expression profiling revealed that TcZelda was mainly expressed in the larval head and adult ovary of late adults and late larvae. TcZelda knockdown led to a 95% mortality rate in adults. These results suggested that TcZelda is related to the activation of the zygote genome in early embryonic development. Furthermore, 592 differentially expressed genes were identified from the dsZelda treated group. Compared with the control group, altered disjunction (ALD) and AGAP005368-PA (GAP) in the dsZelda group were significantly down-regulated, while TGF-beta, propeptide (TGF) was significantly up-regulated, suggesting that TcZelda may be involved in insect embryonic development. In addition, the expression of Ubx ultrabithorax (UBX), Cx cephalothorax (CX), En engrailed (EN), and two Endocuticle structural glycoprotein sgabd (ABD) genes were significantly down-regulated, suggesting that they may cooperate with TcZelda to regulate the development of insect wings. Additionally, Elongation (ELO), fatty acid synthase (FAS), and fatty acyl-CoA desaturase (FAD) expression was inhibited in dsZelda insects, which could disturb the lipase signaling pathways, thus, disrupting the insect reproductive system and pheromone synthesis. These results may help reveal the function of TcZelda in insects and the role of certain genes in the gene regulatory network and provide new ideas for the prevention and control of T. castaneum.

Spatiotemporal control of gene expression boundaries using a feedforward loop.

BACKGROUND: A feedforward loop (FFL) is commonly observed in several biological networks. The FFL network motif has been mostly studied with respect to variation of the input signal in time, with only a few studies of FFL activity in a spatially distributed system such as morphogen-mediated tissue patterning. However, most morphogen gradients also evolve in time. RESULTS: We studied the spatiotemporal behavior of a coherent FFL in two contexts: (a) a generic, oscillating morphogen gradient and (b) the dorsal-ventral patterning of the early Drosophila embryo by a gradient of the NF-κB homolog dorsal with its early target Twist. In both models, we found features in the dynamics of the intermediate node-phase difference and noise filtering-that were largely independent of the parameterization of the models, and thus were functions of the structure of the FFL itself. In the dorsal gradient model, we also found that proper target gene expression was not possible without including the effect of maternal pioneer factor Zelda. CONCLUSIONS: An FFL buffers fluctuation to changes in the morphogen signal ensuring stable gene expression boundaries.

American Dreaming: Really Reading The Great Gatsby.

F. Scott Fitzgerald's The Great Gatsby (1925) is one of the best known and most widely read and taught novels in American literature. It is so familiar that even those who have not read it believe that they have and take for granted that they know about its main character and theme of the American Dream. We need to approach The Great Gatsby as if it were new and really read it, paying close attention to Fitzgerald's literary language. His novel gives us a vivid depiction of and insight into income inequality as it existed in the 1920s and, by extension, as it exists today, when the American Dream is even more limited to the fortunate few, not within reach of the many. When we really read The Great Gatsby, we perceive and understand the American dimension of the novel and appreciate, too, the global range and relevance that in it Fitzgerald has achieved. It is a great American book and a great book of world literature.

Initiation of diverse epigenetic states during nuclear programming of the Drosophila body plan.

Epigenetic patterns of histone modifications contribute to the maintenance of tissue-specific gene expression. Here, we show that such modifications also accompany the specification of cell identities by the NF-κB transcription factor Dorsal in the precellular Drosophila embryo. We provide evidence that the maternal pioneer factor, Zelda, is responsible for establishing poised RNA polymerase at Dorsal target genes before Dorsal-mediated zygotic activation. At the onset of cell specification, Dorsal recruits the CBP/p300 coactivator to the regulatory regions of defined target genes in the presumptive neuroectoderm, resulting in their histone acetylation and transcriptional activation. These genes are inactive in the mesoderm due to transcriptional quenching by the Snail repressor, which precludes recruitment of CBP and prevents histone acetylation. By contrast, inactivation of the same enhancers in the dorsal ectoderm is associated with Polycomb-repressed H3K27me3 chromatin. Thus, the Dorsal morphogen gradient produces three distinct histone signatures including two modes of transcriptional repression, active repression (hypoacetylation), and inactivity (H3K27me3). Whereas histone hypoacetylation is associated with a poised polymerase, H3K27me3 displaces polymerase from chromatin. Our results link different modes of RNA polymerase regulation to separate epigenetic patterns and demonstrate that developmental determinants orchestrate differential chromatin states, providing new insights into the link between epigenetics and developmental patterning.

Co-activation of microRNAs by Zelda is essential for early Drosophila development.

Transcription factors and microRNAs (miRNAs) are two important classes of trans-regulators in differential gene expression. Transcription factors occupy cis-regulatory motifs in DNA to activate or repress gene transcription, whereas miRNAs specifically pair with seed sites in target mRNAs to trigger mRNA decay or inhibit translation. Dynamic spatiotemporal expression patterns of transcription factors and miRNAs during development point to their stage- and tissue-specific functions. Recent studies have focused on miRNA functions during development; however, much remains to explore regarding how the expression of miRNAs is initiated and how dynamic miRNA expression patterns are achieved by transcriptional regulatory networks at different developmental stages. Here, we focused on the identification, regulation and function of miRNAs during the earliest stage of Drosophila development, when the maternal-to-zygotic transition (MZT) takes place. Eleven miRNA clusters comprise the first set of miRNAs activated in the blastoderm embryo. The transcriptional activator Zelda is required for their proper activation and regulation, and Zelda binding observed in genome-wide binding profiles is predictive of enhancer activity. In addition, other blastoderm transcription factors, comprising both activators and repressors, the activities of which are potentiated and coordinated by Zelda, contribute to the accurate temporal and spatial expression of these miRNAs, which are known to function in diverse developmental processes. Although previous genetic studies showed no early phenotypes upon loss of individual miRNAs, our analysis of the miR-1; miR-9a double mutant revealed defects in gastrulation, demonstrating the importance of co-activation of miRNAs by Zelda during the MZT.

Comparing mRNA levels using in situ hybridization of a target gene and co-stain.

In situ hybridization is an important technique for measuring the spatial expression patterns of mRNA in cells, tissues, and whole animals. However, mRNA levels cannot be compared across experiments using typical protocols. Here we present a semi-quantitative method to compare mRNA levels of a gene across multiple samples. This method yields an estimate of the error in the measurement to allow statistical comparison. Our method uses a typical in situ hybridization protocol to stain for a target gene and an internal standard, which we refer to as a co-stain. As a proof of concept, we apply this method to multiple lines of transgenic Drosophila embryos, harboring constructs that express reporter genes to different levels. We generated this test set by mutating enhancer sequences to contain different numbers of binding sites for Zelda, a transcriptional activator. We demonstrate that using a co-stain with in situ hybridization is an effective method to compare mRNA levels across samples. This method requires only minor modifications to existing in situ hybridization protocols and uses straightforward analysis techniques. This strategy can be broadly applied to detect quantitative, spatially resolved changes in mRNA levels.

The embryonic transcription factor Zelda of Drosophila melanogaster is also expressed in larvae and may regulate developmentally important genes.

The transcription factor Zelda plays a pivotal role in promoting the maternal to zygotic transition during embryogenesis in Drosophila melanogaster. However, little is known about its role later in development. Here we are showing that Zelda is essential for proper wing development through gain and loss of function experiments. Zelda's transcript variants RB, RC and RD are present in imaginal wing discs of third instar larvae and the production of 2 protein isoforms of ∼180 and ∼70kD was detected in the same tissue. In ChIP experiments using larval wing discs, Zelda was found to bind to a region of the optomotor-blind gene, suggesting an interaction with a Dpp target that promotes wing growth and patterning.

Chromatin accessibility in the Drosophila embryo is determined by transcription factor pioneering and enhancer activation.

Chromatin accessibility is integral to the process by which transcription factors (TFs) read out cis-regulatory DNA sequences, but it is difficult to differentiate between TFs that drive accessibility and those that do not. Deep learning models that learn complex sequence rules provide an unprecedented opportunity to dissect this problem. Using zygotic genome activation in Drosophila as a model, we analyzed high-resolution TF binding and chromatin accessibility data with interpretable deep learning and performed genetic validation experiments. We identify a hierarchical relationship between the pioneer TF Zelda and the TFs involved in axis patterning. Zelda consistently pioneers chromatin accessibility proportional to motif affinity, whereas patterning TFs augment chromatin accessibility in sequence contexts where they mediate enhancer activation. We conclude that chromatin accessibility occurs in two tiers: one through pioneering, which makes enhancers accessible but not necessarily active, and the second when the correct combination of TFs leads to enhancer activation.

CLAMP regulates zygotic genome activation in Drosophila embryos.

Embryonic patterning is critically dependent on zygotic genome activation (ZGA). In Drosophila melanogaster embryos, the pioneer factor Zelda directs ZGA, possibly in conjunction with other factors. Here, we have explored the novel involvement of Chromatin-Linked Adapter for MSL Proteins (CLAMP) during ZGA. CLAMP binds thousands of sites genome-wide throughout early embryogenesis. Interestingly, CLAMP relocates to target promoter sequences across the genome when ZGA is initiated. Although there is a considerable overlap between CLAMP and Zelda binding sites, the proteins display distinct temporal dynamics. To assess whether CLAMP occupancy affects gene expression, we analyzed transcriptomes of embryos zygotically compromised for either clamp or zelda and found that transcript levels of many zygotically activated genes are similarly affected. Importantly, compromising either clamp or zelda disrupted the expression of critical segmentation and sex determination genes bound by CLAMP (and Zelda). Furthermore, clamp knockdown embryos recapitulate other phenotypes observed in Zelda-depleted embryos, including nuclear division defects, centrosome aberrations, and a disorganized actomyosin network. Based on these data, we propose that CLAMP acts in concert with Zelda to regulate early zygotic transcription.

Spatial organization of transcribing loci during early genome activation in Drosophila.

The early Drosophila embryo provides unique experimental advantages for addressing fundamental questions of gene regulation at multiple levels of organization, from individual gene loci to the entire genome. Using 1.5-h-old Drosophila embryos undergoing the first wave of genome activation,1 we detected ∼110 discrete "speckles" of RNA polymerase II (RNA Pol II) per nucleus, two of which were larger and localized to the histone locus bodies (HLBs).2,3 In the absence of the primary driver of Drosophila genome activation, the pioneer factor Zelda (Zld),1,4,5 70% fewer speckles were present; however, the HLBs tended to be larger than wild-type (WT) HLBs, indicating that RNA Pol II accumulates at the HLBs in the absence of robust early-gene transcription. We observed a uniform distribution of distances between active genes in the nuclei of both WT and zld mutant embryos, indicating that early co-regulated genes do not cluster into nuclear sub-domains. However, in instances whereby transcribing genes did come into close 3D proximity (within 400 nm), they were found to have distinct RNA Pol II speckles. In contrast to the emerging model whereby active genes are clustered to facilitate co-regulation and sharing of transcriptional resources, our data support an "individualist" model of gene control at early genome activation in Drosophila. This model is in contrast to a "collectivist" model, where active genes are spatially clustered and share transcriptional resources, motivating rigorous tests of both models in other experimental systems.

CLAMP and Zelda function together to promote Drosophila zygotic genome activation.

During the essential and conserved process of zygotic genome activation (ZGA), chromatin accessibility must increase to promote transcription. Drosophila is a well-established model for defining mechanisms that drive ZGA. Zelda (ZLD) is a key pioneer transcription factor (TF) that promotes ZGA in the Drosophila embryo. However, many genomic loci that contain GA-rich motifs become accessible during ZGA independent of ZLD. Therefore, we hypothesized that other early TFs that function with ZLD have not yet been identified, especially those that are capable of binding to GA-rich motifs such as chromatin-linked adaptor for male-specific lethal (MSL) proteins (CLAMP). Here, we demonstrate that Drosophila embryonic development requires maternal CLAMP to (1) activate zygotic transcription; (2) increase chromatin accessibility at promoters of specific genes that often encode other essential TFs; and (3) enhance chromatin accessibility and facilitate ZLD occupancy at a subset of key embryonic promoters. Thus, CLAMP functions as a pioneer factor that plays a targeted yet essential role in ZGA.

Formation, interpretation, and regulation of the Drosophila Dorsal/NF-κB gradient.

The morphogen gradient of the transcription factor Dorsal in the early Drosophila embryo has become one of the most widely studied tissue patterning systems. Dorsal is a Drosophila homolog of mammalian NF-κB and patterns the dorsal-ventral axis of the blastoderm embryo into several tissue types by spatially regulating upwards of 100 zygotic genes. Recent studies using fluorescence microscopy and live imaging have quantified the Dorsal gradient and its target genes, which has paved the way for mechanistic modeling of the gradient. In this review, we describe the mechanisms behind the initiation of the Dorsal gradient and its regulation of target genes. The main focus of the review is a discussion of quantitative and computational studies of the Dl gradient system, including regulation of the Dl gradient. We conclude with a discussion of potential future directions.

Quantitative dissection of transcription in development yields evidence for transcription-factor-driven chromatin accessibility.

Thermodynamic models of gene regulation can predict transcriptional regulation in bacteria, but in eukaryotes, chromatin accessibility and energy expenditure may call for a different framework. Here, we systematically tested the predictive power of models of DNA accessibility based on the Monod-Wyman-Changeux (MWC) model of allostery, which posits that chromatin fluctuates between accessible and inaccessible states. We dissected the regulatory dynamics of hunchback by the activator Bicoid and the pioneer-like transcription factor Zelda in living Drosophila embryos and showed that no thermodynamic or non-equilibrium MWC model can recapitulate hunchback transcription. Therefore, we explored a model where DNA accessibility is not the result of thermal fluctuations but is catalyzed by Bicoid and Zelda, possibly through histone acetylation, and found that this model can predict hunchback dynamics. Thus, our theory-experiment dialogue uncovered potential molecular mechanisms of transcriptional regulatory dynamics, a key step toward reaching a predictive understanding of developmental decision-making.

Cells in the brain, liver and skin, as well as many other organs, all contain the same DNA, yet behave in very different ways. This is because before a gene can produce its corresponding protein, it must first be transcribed into messenger RNA. As an organism grows, the transcription of certain genes is switched on or off by regulatory molecules called transcription factors, which guide cells towards a specific 'fate'. These molecules bind to specific locations within the regulatory regions of DNA, and for decades biologist have tried to use the arrangement of these sites to predict which proteins a cell will make. Theoretical models known as thermodynamic models have been able to successfully predict transcription in bacteria. However, this has proved more challenging to do in eukaryotes, such as yeast, fruit flies and humans. One of the key differences is that DNA in eukaryotes is typically tightly wound into bundles called nucleosomes, which must be disentangled in order for transcription factors to access the DNA. Previous thermodynamic models have suggested that DNA in eukaryotes randomly switches between being in a wound and unwound state. The models assume that once unwound, regulatory proteins stabilize the DNA in this form, making it easier for other transcription factors to bind to the DNA. Now, Eck, Liu et al. have tested some of these models by studying the transcription of a gene involved in the development of fruit flies. The experiments showed that no thermodynamic model could accurately mimic how this gene is regulated in the embryos of fruit flies. This led Eck, Liu et al. to identify a model that is better at predicting the activation pattern of this developmental gene. In this model, instead of just 'locking' DNA into an unwound shape, transcription factors can also actively speed up the unwinding of DNA. This improved understanding builds towards the goal of predicting gene regulation, where DNA sequences can be used to tell where and when cell decisions will be made. In the future, this could allow the development of new types of therapies that can regulate transcription in different diseases.

Odd-paired is a pioneer-like factor that coordinates with Zelda to control gene expression in embryos.

Pioneer factors such as Zelda (Zld) help initiate zygotic transcription in Drosophila early embryos, but whether other factors support this dynamic process is unclear. Odd-paired (Opa), a zinc-finger transcription factor expressed at cellularization, controls the transition of genes from pair-rule to segmental patterns along the anterior-posterior axis. Finding that Opa also regulates expression through enhancer sog_Distal along the dorso-ventral axis, we hypothesized Opa's role is more general. Chromatin-immunoprecipitation (ChIP-seq) confirmed its in vivo binding to sog_Distal but also identified widespread binding throughout the genome, comparable to Zld. Furthermore, chromatin assays (ATAC-seq) demonstrate that Opa, like Zld, influences chromatin accessibility genome-wide at cellularization, suggesting both are pioneer factors with common as well as distinct targets. Lastly, embryos lacking opa exhibit widespread, late patterning defects spanning both axes. Collectively, these data suggest Opa is a general timing factor and likely late-acting pioneer factor that drives a secondary wave of zygotic gene expression.

The Drosophila Pioneer Factor Zelda Modulates the Nuclear Microenvironment of a Dorsal Target Enhancer to Potentiate Transcriptional Output.

Connecting the developmental patterning of tissues to the mechanistic control of RNA polymerase II remains a long-term goal of developmental biology. Many key elements have been identified in the establishment of spatial-temporal control of transcription in the early Drosophila embryo, a model system for transcriptional regulation. The dorsal-ventral axis of the Drosophila embryo is determined by the graded distribution of Dorsal (Dl), a homolog of the nuclear factor κB (NF-κB) family of transcriptional activators found in humans [1, 2]. A second maternally deposited factor, Zelda (Zld), is uniformly distributed in the embryo and is thought to act as a pioneer factor, increasing enhancer accessibility for transcription factors, such as Dl [3-9]. Here, we utilized the MS2 live imaging system to evaluate the expression of the Dl target gene short gastrulation (sog) to better understand how a pioneer factor affects the kinetic parameters of transcription. Our experiments indicate that Zld modifies probability of activation, the timing of this activation, and the rate at which transcription occurs. Our results further show that this effective rate increase is due to an increased accumulation of Dl at the site of transcription, suggesting that transcription factor "hubs" induced by Zld [10] functionally regulate transcription.

Zelda and the maternal-to-zygotic transition in cockroaches.

In the endopterygote Drosophila melanogaster, Zelda is an activator of the zygotic genome during the maternal-to-zygotic transition (MZT). Zelda binds cis-regulatory elements (TAGteam heptamers), making chromatin accessible for gene transcription. Zelda has been studied in other endopterygotes: Apis mellifera and Tribolium castaneum, and the paraneopteran Rhodnius prolixus. We studied Zelda in the cockroach Blattella germanica, a hemimetabolan, short germ-band, and polyneopteran species. B. germanica Zelda has the complete set of functional domains, which is typical of species displaying ancestral features concerning embryogenesis. Interestingly, we found D. melanogaster TAGteam heptamers in the B. germanica genome. The canonical one, CAGGTAG, is present at a similar proportion in the genome of these two species and in the genome of other insects, suggesting that the genome admits as many CAGGTAG motifs as its length allows. Zelda-depleted embryos of B. germanica show defects involving blastoderm formation and abdomen development, and genes contributing to these processes are down-regulated. We conclude that in B. germanica, Zelda strictly activates the zygotic genome, within the MZT, a role conserved in more derived endopterygote insects. In B. germanica, zelda is expressed during MZT, whereas in D. melanogaster and T. castaneum it is expressed beyond this transition. In these species and A. mellifera, Zelda has functions even in postembryonic development. The expansion of zelda expression beyond the MZT in endopterygotes might be related with the evolutionary innovation of holometabolan metamorphosis. DATABASES: The RNA-seq datasets of B. germanica, D. melanogaster, and T. castaneum are accessible at the GEO databases GSE99785, GSE18068, GSE63770, and GSE84253. In addition, the RNA-seq library from T. castaneum adult females is available at SRA: SRX021963. The B. germanica reference genome is available as BioProject PRJNA203136.