STAT is an essential activator of the zygotic genome in the early Drosophila embryo

Exploring the versatility of zygotic genome regulators: A comparative and functional analysis

The activation of the zygotic genome constitutes an essential process during early embryogenesis that determines the overall progression of embryonic development. Zygotic genome activation (ZGA) is tightly regulated, involving a delicate interplay of activators and repressors, to precisely control the timing and spatial pattern of gene expression. While regulators of ZGA vary across species, they accomplish comparable outcomes. Recent studies have shed light on the unanticipated roles of ZGA components both during and after ZGA. Moreover, different ZGA regulators seem to have acquired unique functional modalities to manifest their regulatory potential. In this review, we explore these observations to assess whether these are simply anecdotal or contribute to a broader regulatory framework that employs a versatile means to arrive at the conserved outcome.

CBP and Gcn5 drive zygotic genome activation independently of their catalytic activity

Zygotic genome activation (ZGA) is a crucial step of embryonic development. So far, little is known about the role of chromatin factors during this process. Here, we used an in vivo RNA interference reverse genetic screen to identify chromatin factors necessary for embryonic development in Drosophila melanogaster. Our screen reveals that histone acetyltransferases (HATs) and histone deacetylases are crucial ZGA regulators. We demonstrate that Nejire (CBP/EP300 ortholog) is essential for the acetylation of histone H3 lysine-18 and lysine-27, whereas Gcn5 (GCN5/PCAF ortholog) for lysine-9 of H3 at ZGA, with these marks being enriched at all actively transcribed genes. Nonetheless, these HATs activate distinct sets of genes. Unexpectedly, individual catalytic dead mutants of either Nejire or Gcn5 can activate zygotic transcription (ZGA) and transactivate a reporter gene in vitro. Together, our data identify Nejire and Gcn5 as key regulators of ZGA.

Transcriptional activators in the early Drosophila embryo perform different kinetic roles

Combinatorial regulation of gene expression by transcription factors (TFs) may in part arise from kinetic synergy-wherein TFs regulate different steps in the transcription cycle. Kinetic synergy requires that TFs play distinguishable kinetic roles. Here, we used live imaging to determine the kinetic roles of three TFs that activate transcription in the Drosophila embryo-Zelda, Bicoid, and Stat92E-by introducing their binding sites into the even-skipped stripe 2 enhancer. These TFs influence different sets of kinetic parameters, and their influence can change over time. All three TFs increased the fraction of transcriptionally active nuclei; Zelda also shortened the first-passage time into transcription and regulated the interval between transcription events. Stat92E also increased the lifetimes of active transcription. Different TFs can therefore play distinct kinetic roles in activating the transcription. This has consequences for understanding the composition and flexibility of regulatory DNA sequences and the biochemical function of TFs. A record of this paper's transparent peer review process is included in the supplemental information.

The Drosophila AWP1 ortholog Doctor No regulates JAK/STAT signaling for left-right asymmetry in the gut by promoting receptor endocytosis

Many organs of Drosophila show stereotypical left-right (LR) asymmetry; however, the underlying mechanisms remain elusive. Here, we have identified an evolutionarily conserved ubiquitin-binding protein, AWP1/Doctor No (Drn), as a factor required for LR asymmetry in the embryonic anterior gut. We found that drn is essential in the circular visceral muscle cells of the midgut for JAK/STAT signaling, which contributes to the first known cue for anterior gut lateralization via LR asymmetric nuclear rearrangement. Embryos homozygous for drn and lacking its maternal contribution showed phenotypes similar to those with depleted JAK/STAT signaling, suggesting that Drn is a general component of JAK/STAT signaling. Absence of Drn resulted in specific accumulation of Domeless (Dome), the receptor for ligands in the JAK/STAT signaling pathway, in intracellular compartments, including ubiquitylated cargos. Dome colocalized with Drn in wild-type Drosophila. These results suggest that Drn is required for the endocytic trafficking of Dome, which is a crucial step for activation of JAK/STAT signaling and the subsequent degradation of Dome. The roles of AWP1/Drn in activating JAK/STAT signaling and in LR asymmetric development may be conserved in various organisms.

Mechanisms of Interaction between Enhancers and Promoters in Three Drosophila Model Systems

In higher eukaryotes, the regulation of developmental gene expression is determined by enhancers, which are often located at a large distance from the promoters they regulate. Therefore, the architecture of chromosomes and the mechanisms that determine the functional interaction between enhancers and promoters are of decisive importance in the development of organisms. Mammals and the model animal Drosophila have homologous key architectural proteins and similar mechanisms in the organization of chromosome architecture. This review describes the current progress in understanding the mechanisms of the formation and regulation of long-range interactions between enhancers and promoters at three well-studied key regulatory loci in Drosophila.

A timer gene network is spatially regulated by the terminal system in the Drosophila embryo

In insect embryos, anteroposterior patterning is coordinated by the sequential expression of the 'timer' genes caudal, Dichaete, and odd-paired, whose expression dynamics correlate with the mode of segmentation. In Drosophila, the timer genes are expressed broadly across much of the blastoderm, which segments simultaneously, but their expression is delayed in a small 'tail' region, just anterior to the hindgut, which segments during germband extension. Specification of the tail and the hindgut depends on the terminal gap gene tailless, but beyond this the regulation of the timer genes is poorly understood. We used a combination of multiplexed imaging, mutant analysis, and gene network modelling to resolve the regulation of the timer genes, identifying 11 new regulatory interactions and clarifying the mechanism of posterior terminal patterning. We propose that a dynamic Tailless expression gradient modulates the intrinsic dynamics of a timer gene cross-regulatory module, delineating the tail region and delaying its developmental maturation.

Optogenetic control of the Bicoid morphogen reveals fast and slow modes of gap gene regulation

Developmental patterning networks are regulated by multiple inputs and feedback connections that rapidly reshape gene expression, limiting the information that can be gained solely from slow genetic perturbations. Here we show that fast optogenetic stimuli, real-time transcriptional reporters, and a simplified genetic background can be combined to reveal the kinetics of gene expression downstream of a developmental transcription factor in vivo. We engineer light-controlled versions of the Bicoid transcription factor and study their effects on downstream gap genes in embryos. Our results recapitulate known relationships, including rapid Bicoid-dependent transcription of giant and hunchback and delayed repression of Krüppel. In addition, we find that the posterior pattern of knirps exhibits a quick but inverted response to Bicoid perturbation, suggesting a noncanonical role for Bicoid in directly suppressing knirps transcription. Acute modulation of transcription factor concentration while recording output gene activity represents a powerful approach for studying developmental gene networks in vivo.

The Dynamic Changes of Transcription Factors During the Development Processes of Human Biparental and Uniparental Embryos

Previous studies have revealed that transcription factors (TFs) play important roles in biparental (BI) early human embryogenesis. However, the contribution of TFs during early uniparental embryo development is still largely unknown. Here we systematically studied the expression profiles of transcription factors in early embryonic development and revealed the dynamic changes of TFs in human biparental and uniparental embryogenesis by single-cell RNA sequencing (scRNA-seq). In general, the TF expression model of uniparental embryos showed a high degree of conformity with biparental embryos. The detailed network analysis of three different types of embryos identified that 10 out of 17 hub TFs were shared or specifically owned, such as ZNF480, ZNF581, PHB, and POU5F1, were four shared TFs, ZFN534, GTF3A, ZNF771, TEAD4, and LIN28A, were androgenic (AG) specific TFs, and ZFP42 was the only one parthenogenetic (PG) specific TF. All the four shared TFs were validated using human embryonic stem cell (hESC) differentiation experiments; most of their target genes are responsible for stem cell maintenance and differentiation. We also found that Zf-C2H2, HMG, and MYB were three dominant transcription factor families that appeared in early embryogenesis. Altogether, our work provides a comprehensive regulatory framework and better understanding of TF function in human biparental and uniparental embryogenesis.

Molecular competition can shape enhancer activity in the Drosophila embryo

Transgenic reporters allow the measurement of regulatory DNA activity in vivo and consequently have long been useful tools for studying enhancers. Despite their utility, few studies have investigated the effects these reporters may have on the expression of other genes. Understanding these effects is required to accurately interpret reporter data and characterize gene regulatory mechanisms. By measuring the expression of Kruppel (Kr) enhancer reporters in live Drosophila embryos, we find reporters inhibit one another's expression and that of a nearby endogenous gene. Using synthetic transcription factor (TF) binding site arrays, we present evidence that competition for TFs is partially responsible for the observed transcriptional inhibition. We develop a simple thermodynamic model that predicts competition of the measured magnitude specifically when TF binding is restricted to distinct nuclear subregions. Our findings underline an unexpected role of the non-homogenous nature of the nucleus in regulating gene expression.

Regulatory principles governing the maternal-to-zygotic transition: insights from Drosophila melanogaster

The onset of metazoan development requires that two terminally differentiated germ cells, a sperm and an oocyte, become reprogrammed to the totipotent embryo, which can subsequently give rise to all the cell types of the adult organism. In nearly all animals, maternal gene products regulate the initial events of embryogenesis while the zygotic genome remains transcriptionally silent. Developmental control is then passed from mother to zygote through a process known as the maternal-to-zygotic transition (MZT). The MZT comprises an intimately connected set of molecular events that mediate degradation of maternally deposited mRNAs and transcriptional activation of the zygotic genome. This essential developmental transition is conserved among metazoans but is perhaps best understood in the fruit fly, Drosophila melanogaster. In this article, we will review our understanding of the events that drive the MZT in Drosophila embryos and highlight parallel mechanisms driving this transition in other animals.

Synthetic enhancer design by in silico compensatory evolution reveals flexibility and constraint in cis-regulation

BACKGROUND

Models that incorporate specific chemical mechanisms have been successful in describing the activity of Drosophila developmental enhancers as a function of underlying transcription factor binding motifs. Despite this, the minimum set of mechanisms required to reconstruct an enhancer from its constituent parts is not known. Synthetic biology offers the potential to test the sufficiency of known mechanisms to describe the activity of enhancers, as well as to uncover constraints on the number, order, and spacing of motifs.

RESULTS

Using a functional model and in silico compensatory evolution, we generated putative synthetic even-skipped stripe 2 enhancers with varying degrees of similarity to the natural enhancer. These elements represent the evolutionary trajectories of the natural stripe 2 enhancer towards two synthetic enhancers designed ab initio. In the first trajectory, spatially regulated expression was maintained, even after more than a third of binding sites were lost. In the second, sequences with high similarity to the natural element did not drive expression, but a highly diverged sequence about half the length of the minimal stripe 2 enhancer drove ten times greater expression. Additionally, homotypic clusters of Zelda or Stat92E motifs, but not Bicoid, drove expression in developing embryos.

CONCLUSIONS

Here, we present a functional model of gene regulation to test the degree to which the known transcription factors and their interactions explain the activity of the Drosophila even-skipped stripe 2 enhancer. Initial success in the first trajectory showed that the gene regulation model explains much of the function of the stripe 2 enhancer. Cases where expression deviated from prediction indicates that undescribed factors likely act to modulate expression. We also showed that activation driven Bicoid and Hunchback is highly sensitive to spatial organization of binding motifs. In contrast, Zelda and Stat92E drive expression from simple homotypic clusters, suggesting that activation driven by these factors is less constrained. Collectively, the 40 sequences generated in this work provides a powerful training set for building future models of gene regulation.

Canonical and non-canonical JAK/STAT transcriptional targets may be involved in distinct and overlapping cellular processes

Background

The Janus kinase-signal transducer and activator of transcription (JAK/STAT) pathway has been well-characterized as a crucial signal transduction cascade that regulates vital biological responses including development, immunity and oncogenesis. Additionally to its canonical pathway that uses the phosphorylated form of the STAT transcription factor, recently the non-canonical pathway involving heterochromatin formation by unphosphorylated STAT was recently uncovered. Considering the significant role of the JAK/STAT pathway, we used the simple Drosophila system in which the non-canonical pathway was initially characterized, to compare putative canonical versus non-canonical transcriptional targets across the genome. We analyzed microarray expression patterns of wildtype, Jak gain- and loss-of-function mutants, as well as the Stat loss-of-function mutant during embryogenesis, since the contribution of the canonical signal transduction pathway has been well-characterized in these contexts. Previous studies have also demonstrated that Jak gain-of-function and Stat mutants counter heterochromatin silencing to de-repress target genes by the non-canonical pathway.

Results

Compared to canonical target genomic loci, non-canonical targets were significantly more associated with sites enriched with heterochromatin-related factors (p = 0.004). Furthermore, putative canonical and non-canonical transcriptional targets identified displayed some differences in biological pathways they regulate, as determined by Gene Ontology (GO) enrichment analyses. Canonical targets were enriched mainly with genes relevant to development and immunity, as expected, whereas the non-canonical target gene set mainly showed enrichment of genes for various metabolic responses and stress response, highlighting the possibility that some differences may exist between the two loci.

Conclusions

Canonical and non-canonical JAK/STAT genes may regulate distinct and overlapping sets of genes and may perform specific overall functions in physiology. Further studies at different developmental stages, or using distinct tissues may identify additional targets and provide insight into which gene targets are unique to the canonical or non-canonical pathway.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4058-y) contains supplementary material, which is available to authorized users.

A sequence level model of an intact locus predicts the location and function of nonadditive enhancers

Metazoan gene expression is controlled through the action of long stretches of noncoding DNA that contain enhancers-shorter sequences responsible for controlling a single aspect of a gene's expression pattern. Models built on thermodynamics have shown how enhancers interpret protein concentration in order to determine specific levels of gene expression, but the emergent regulatory logic of a complete regulatory locus shows qualitative and quantitative differences from isolated enhancers. Such differences may arise from steric competition limiting the quantity of DNA that can simultaneously influence the transcription machinery. We incorporated this competition into a mechanistic model of gene regulation, generated efficient algorithms for this computation, and applied it to the regulation of Drosophila even-skipped (eve). This model finds the location of enhancers and identifies which factors control the boundaries of eve expression. This model predicts a new enhancer that, when assayed in vivo, drives expression in a non-eve pattern. Incorporation of chromatin accessibility eliminates this inconsistency.

The human PRD-like homeobox gene LEUTX has a central role in embryo genome activation

Leucine twenty homeobox (LEUTX) is a paired (PRD)-like homeobox gene that is expressed almost exclusively in human embryos during preimplantation development. We previously identified a novel transcription start site for the predicted human LEUTX gene based on the transcriptional analysis of human preimplantation embryos. The novel variant encodes a protein with a complete homeodomain. Here, we provide a detailed description of the molecular cloning of the complete homeodomain-containing LEUTX Using a human embryonic stem cell overexpression model we show that the complete homeodomain isoform is functional and sufficient to activate the transcription of a large proportion of the genes that are upregulated in human embryo genome activation (EGA), whereas the previously predicted partial homeodomain isoform is largely inactive. Another PRD-like transcription factor, DPRX, is then upregulated as a powerful repressor of transcription. We propose a two-stage model of human EGA in which LEUTX acts as a transcriptional activator at the 4-cell stage, and DPRX as a balancing repressor at the 8-cell stage. We conclude that LEUTX is a candidate regulator of human EGA.

lolal Is an Evolutionarily New Epigenetic Regulator of dpp Transcription during Dorsal-Ventral Axis Formation

Secreted ligands in the Dpp/BMP family drive dorsal–ventral (D/V) axis formation in all Bilaterian species. However, maternal factors regulating Dpp/BMP transcription in this process are largely unknown. We identified the BTB domain protein longitudinals lacking-like (lolal) as a modifier of decapentaplegic (dpp) mutations. We show that Lolal is evolutionarily related to the Trithorax group of chromatin regulators and that lolal interacts genetically with the epigenetic factor Trithorax-like during Dpp D/V signaling. Maternally driven Lolal^HA is found in oocytes and translocates to zygotic nuclei prior to the point at which dpp transcription begins. lolal maternal and zygotic mutant embryos display significant reductions in dpp, pMad, and zerknullt expression, but they are never absent. The data suggest that lolal is required to maintain dpp transcription during D/V patterning. Phylogenetic data revealed that lolal is an evolutionarily new gene present only in insects and crustaceans. We conclude that Lolal is the first maternal protein identified with a role in dpp D/V transcriptional maintenance, that Lolal and the epigenetic protein Trithorax-like are essential for Dpp D/V signaling and that the architecture of the Dpp D/V pathway evolved in the arthropod lineage after the separation from vertebrates via the incorporation of new genes such as lolal.

Stepwise Progression of Embryonic Patterning

It is long established that the graded distribution of Dorsal transcription factor influences spatial domains of gene expression along the dorsoventral (DV) axis of Drosophila melanogaster embryos. However, the more recent realization that Dorsal levels also change with time raises the question of whether these dynamics are instructive. An overview of DV axis patterning is provided, focusing on new insights identified through quantitative analysis of temporal changes in Dorsal target gene expression from one nuclear cycle to the next ('steps'). Possible roles for the stepwise progression of this patterning program are discussed including (i) tight temporal regulation of signaling pathway activation, (ii) control of gene expression cohorts, and (iii) ensuring the irreversibility of the patterning and cell fate specification process.

Neural Progenitors Adopt Specific Identities by Directly Repressing All Alternative Progenitor Transcriptional Programs

In the vertebrate neural tube, a morphogen-induced transcriptional network produces multiple molecularly distinct progenitor domains, each generating different neuronal subtypes. Using an in vitro differentiation system, we defined gene expression signatures of distinct progenitor populations and identified direct gene-regulatory inputs corresponding to locations of specific transcription factor binding. Combined with targeted perturbations of the network, this revealed a mechanism in which a progenitor identity is installed by active repression of the entire transcriptional programs of other neural progenitor fates. In the ventral neural tube, sonic hedgehog (Shh) signaling, together with broadly expressed transcriptional activators, concurrently activates the gene expression programs of several domains. The specific outcome is selected by repressive input provided by Shh-induced transcription factors that act as the key nodes in the network, enabling progenitors to adopt a single definitive identity from several initially permitted options. Together, the data suggest design principles relevant to many developing tissues.

•

Specific vertebrate neural progenitor populations generated in vitro

•

Gene expression dynamics, transcription factor binding assessed in neural progenitors

•

Progenitor fate selected by repressors blocking entire programs of other identities

•

Repressors counteract non-selective morphogen and pan-neural activatory inputs

Krüppel Expression Levels Are Maintained through Compensatory Evolution of Shadow Enhancers

Many developmental genes are controlled by shadow enhancers—pairs of enhancers that drive overlapping expression patterns. We hypothesized that compensatory evolution can maintain the total expression of a gene, while individual shadow enhancers diverge between species. To test this hypothesis, we analyzed expression driven by orthologous pairs of shadow enhancers from Drosophila melanogaster, Drosophila yakuba, and Drosophila pseudoobscura that control expression of Krüppel, a transcription factor that patterns the anterior-posterior axis of blastoderm embryos. We found that the expression driven by the pair of enhancers is conserved between these three species, but expression levels driven by the individual enhancers are not. Using sequence analysis and experimental perturbation, we show that each shadow enhancer is regulated by different transcription factors. These results support the hypothesis that compensatory evolution can occur between shadow enhancers, which has implications for mechanistic and evolutionary studies of gene regulation.

Regulation and Function of Maternal Gene Products During the Maternal-to-Zygotic Transition in Drosophila

Drosophila late-stage oocytes and early embryos are transcriptionally silent. Thus, control of gene expression during these developmental periods is posttranscriptional and posttranslational. Global changes in the transcriptome and proteome occur during oocyte maturation, after egg activation and fertilization, and upon zygotic genome activation. We review the scale, content, and dynamics of these global changes; the factors that regulate these changes; and the mechanisms by which they are accomplished. We highlight the intimate relationship between the clearance of maternal gene products and the activation of the embryo's own genome, and discuss the fact that each of these complementary components of the maternal-to-zygotic transition can be subdivided into several phases that serve different biological roles and are regulated by distinct factors.

What does it take to evolve an enhancer? A simulation-based study of factors influencing the emergence of combinatorial regulation

There is widespread interest today in understanding enhancers, which are regulatory elements typically harboring several transcription factor binding sites and mediating the combinatorial effect of transcription factors on gene expression. The evolution of enhancers poses interesting unanswered questions, for example, the evolutionary time taken for a typical enhancer to emerge or the factors shaping its evolution. Existing approaches to cis-regulatory evolution have often ignored the combinatorial nature and varied biochemical mechanisms of gene regulation encoded in enhancers. We report on our investigation of enhancer evolution through the use of PEBCRES, a framework for evolutionary simulation of enhancers that employs a mechanistic and well-supported sequence-to-expression model to assign fitness to the evolving enhancer genotype. We estimated the time necessary to evolve, from genomic background, enhancers capable of driving complex gene expression patterns similar to those involved in early development in Drosophila. We found the time-to-evolve to range between 0.5 and 10 Myr, and to vary greatly with the target expression pattern, complexity of the real enhancer known to encode that pattern, and the strength of input from specific transcription factors. To our knowledge, this is the first estimate of waiting times for realistic enhancers to evolve. The in silico evolved enhancers had, with a few interesting exceptions, site compositions similar to those seen in real enhancers for the same patterns. Our simulations also revealed that certain features of an enhancer might evolve not due to their biological function but as aids to the evolutionary process itself.

Discovery of transcription factors and regulatory regions driving in vivo tumor development by ATAC-seq and FAIRE-seq open chromatin profiling

Genomic enhancers regulate spatio-temporal gene expression by recruiting specific combinations of transcription factors (TFs). When TFs are bound to active regulatory regions, they displace canonical nucleosomes, making these regions biochemically detectable as nucleosome-depleted regions or accessible/open chromatin. Here we ask whether open chromatin profiling can be used to identify the entire repertoire of active promoters and enhancers underlying tissue-specific gene expression during normal development and oncogenesis in vivo. To this end, we first compare two different approaches to detect open chromatin in vivo using the Drosophila eye primordium as a model system: FAIRE-seq, based on physical separation of open versus closed chromatin; and ATAC-seq, based on preferential integration of a transposon into open chromatin. We find that both methods reproducibly capture the tissue-specific chromatin activity of regulatory regions, including promoters, enhancers, and insulators. Using both techniques, we screened for regulatory regions that become ectopically active during Ras-dependent oncogenesis, and identified 3778 regions that become (over-)activated during tumor development. Next, we applied motif discovery to search for candidate transcription factors that could bind these regions and identified AP-1 and Stat92E as key regulators. We validated the importance of Stat92E in the development of the tumors by introducing a loss of function Stat92E mutant, which was sufficient to rescue the tumor phenotype. Additionally we tested if the predicted Stat92E responsive regulatory regions are genuine, using ectopic induction of JAK/STAT signaling in developing eye discs, and observed that similar chromatin changes indeed occurred. Finally, we determine that these are functionally significant regulatory changes, as nearby target genes are up- or down-regulated. In conclusion, we show that FAIRE-seq and ATAC-seq based open chromatin profiling, combined with motif discovery, is a straightforward approach to identify functional genomic regulatory regions, master regulators, and gene regulatory networks controlling complex in vivo processes.

Coordinate regulation of stem cell competition by Slit-Robo and JAK-STAT signaling in the Drosophila testis

Stem cells in tissues reside in and receive signals from local microenvironments called niches. Understanding how multiple signals within niches integrate to control stem cell function is challenging. The Drosophila testis stem cell niche consists of somatic hub cells that maintain both germline stem cells and somatic cyst stem cells (CySCs). Here, we show a role for the axon guidance pathway Slit-Roundabout (Robo) in the testis niche. The ligand Slit is expressed specifically in hub cells while its receptor, Roundabout 2 (Robo2), is required in CySCs in order for them to compete for occupancy in the niche. CySCs also require the Slit-Robo effector Abelson tyrosine kinase (Abl) to prevent over-adhesion of CySCs to the niche, and CySCs mutant for Abl outcompete wild type CySCs for niche occupancy. Both Robo2 and Abl phenotypes can be rescued through modulation of adherens junction components, suggesting that the two work together to balance CySC adhesion levels. Interestingly, expression of Robo2 requires JAK-STAT signaling, an important maintenance pathway for both germline and cyst stem cells in the testis. Our work indicates that Slit-Robo signaling affects stem cell function downstream of the JAK-STAT pathway by controlling the ability of stem cells to compete for occupancy in their niche.

Zygotic genome activation during the maternal-to-zygotic transition

Embryogenesis depends on a highly coordinated cascade of genetically encoded events. In animals, maternal factors contributed by the egg cytoplasm initially control development, whereas the zygotic nuclear genome is quiescent. Subsequently, the genome is activated, embryonic gene products are mobilized, and maternal factors are cleared. This transfer of developmental control is called the maternal-to-zygotic transition (MZT). In this review, we discuss recent advances toward understanding the scope, timing, and mechanisms that underlie zygotic genome activation at the MZT in animals. We describe high-throughput techniques to measure the embryonic transcriptome and explore how regulation of the cell cycle, chromatin, and transcription factors together elicits specific patterns of embryonic gene expression. Finally, we illustrate the interplay between zygotic transcription and maternal clearance and show how these two activities combine to reprogram two terminally differentiated gametes into a totipotent embryo.

The role of Dichaete in transcriptional regulation during Drosophila embryonic development

Background

Group B Sox domain transcription factors play conserved roles in the specification and development of the nervous system in higher metazoans. However, we know comparatively little about how these transcription factors regulate gene expression, and the analysis of Sox gene function in vertebrates is confounded by functional compensation between three closely related family members. In Drosophila, only two group B Sox genes, Dichaete and SoxN, have been shown to function during embryonic CNS development, providing a simpler system for understanding the functions of this important class of regulators.

Results

Using a combination of transcriptional profiling and genome-wide binding analysis we conservatively identify over 1000 high confidence direct Dichaete target genes in the Drosophila genome. We show that Dichaete plays key roles in CNS development, regulating aspects of the temporal transcription factor sequence that confer neuroblast identity. Dichaete also shows a complex interaction with Prospero in the pathway controlling the switch from stem cell self-renewal to neural differentiation. Dichaete potentially regulates many more genes in the Drosophila genome and was found to be associated with over 2000 mapped regulatory elements.

Conclusions

Our analysis suggests that Dichaete acts as a transcriptional hub, controlling multiple regulatory pathways during CNS development. These include a set of core CNS expressed genes that are also bound by the related Sox2 gene during mammalian CNS development. Furthermore, we identify Dichaete as one of the transcription factors involved in the neural stem cell transcriptional network, with evidence supporting the view that Dichaete is involved in controlling the temporal series of divisions regulating neuroblast identity.

JAK-STAT pathway in Drosophila morphogenesis: From organ selector to cell behavior regulator

One of the main contributions of Drosophila to the JAK-STAT field is the study of morphogenesis. JAK-STAT signaling controls the formation of many different structures through surprisingly different morphogenetic behaviors that include induction of cell rearrangements, invagination, folding of tissues, modulation of cell shape, and migration. This variability may be explained by the many transcription factors and signaling molecules STAT regulates at early stages of development. But is STAT just acting as an upstream inducer of morphogenesis or does it have a more direct role in controlling cell behaviors? Here we review what is known about how the canonical phosphorylation of STAT contributes to shaping the embryonic and imaginal structures.

Extensive polymerase pausing during Drosophila axis patterning enables high-level and pliable transcription

Cascades of zygotic gene expression pattern the anterior-posterior (AP) and dorsal-ventral (DV) axes of the early Drosophila embryo. Here, we used the global run-on sequencing assay (GRO-seq) to map the genome-wide RNA polymerase distribution during early Drosophila embryogenesis, thus providing insights into how genes are regulated. We identify widespread promoter-proximal pausing yet show that the presence of paused polymerase does not necessarily equate to direct regulation through pause release to productive elongation. Our data reveal that a subset of early Zelda-activated genes is regulated at the level of polymerase recruitment, whereas other Zelda target and axis patterning genes are predominantly regulated through pause release. In contrast to other signaling pathways, we found that bone morphogenetic protein (BMP) target genes are collectively more highly paused than BMP pathway components and show that BMP target gene expression requires the pause-inducing negative elongation factor (NELF) complex. Our data also suggest that polymerase pausing allows plasticity in gene activation throughout embryogenesis, as transiently repressed and transcriptionally silenced genes maintain and lose promoter polymerases, respectively. Finally, we provide evidence that the major effect of pausing is on the levels, rather than timing, of transcription. These data are discussed in terms of the efficiency of transcriptional activation required across cell populations during developmental time constraints.

Context-dependent transcriptional interpretation of mitogen activated protein kinase signaling in the Drosophila embryo

Terminal regions of the Drosophila embryo are patterned by the localized activation of Mitogen Activated Protein Kinase (MAPK), which induces zygotic genes through relief of their repression by transcriptional repressor Capicua. The levels of MAPK activation at the anterior and posterior termini are close to each other, but the expression patterns of MAPK-target genes, such as zerknüllt (zen) and tailless (tll), display strong anterior-posterior (AP) asymmetry. This region-specific response to MAPK activation provides a clear example of context-dependent interpretation of inductive signaling, a common developmental effect that remains poorly understood. In the past, the AP asymmetry of zen expression was attributed to a mechanism that depends on MAPK substrate competition. We present data suggesting that the asymmetric expression of tll is generated by a different mechanism, based on feedforward control and multiple enhancers of the tll gene. A simple mathematical model of this mechanism correctly predicts how the wild-type expression pattern of tll changes in mutants affecting the anterior, dorsoventral, and terminal patterning systems and some of their direct targets.

Uncovering a dynamic feature of the transcriptional regulatory network for anterior-posterior patterning in the Drosophila embryo

Anterior-posterior (AP) patterning in the Drosophila embryo is dependent on the Bicoid (Bcd) morphogen gradient. However, most target genes of Bcd also require additional inputs to establish their expression domains, reflective of the operation of a cross-regulatory network and contributions of other maternal signals. This is in contrast to hunchback (hb), which has an anterior expression domain driven by an enhancer that appears to respond primarily to the Bcd input. To gain a better understanding of the regulatory logic of the AP patterning network, we perform quantitative studies that specifically investigate the dynamics of hb transcription during development. We show that Bcd-dependent hb transcription, monitored by the intron-containing nascent transcripts near the P2 promoter, is turned off quickly–on the order of a few minutes–upon entering the interphase of nuclear cycle 14A. This shutdown contrasts with earlier cycles during which active hb transcription can persist until the moment when the nucleus enters mitosis. The shutdown takes place at a time when the nuclear Bcd gradient profile in the embryo remains largely intact, suggesting that this is a process likely subject to control of a currently unknown regulatory mechanism. We suggest that this dynamic feature offers a window of opportunity for hb to faithfully interpret, and directly benefit from, Bcd gradient properties, including its scaling properties, to help craft a robust AP patterning outcome.

Transcriptional and epigenetic signatures of zygotic genome activation during early Drosophila embryogenesis

Background

In all Metazoa, transcription is inactive during the first mitotic cycles after fertilisation. In Drosophila melanogaster, Zygotic Genome Activation (ZGA) occurs in two waves, starting respectively at mitotic cycles 8 (approximately 60 genes) and 14 (over a thousand genes). The regulatory mechanisms underlying these drastic transcriptional changes remain largely unknown.

Results

We developed an original gene clustering method based on discretized transition profiles, and applied it to datasets from three landmark early embryonic transcriptome studies. We identified 417 genes significantly up-regulated during ZGA. De novo motif discovery returned nine motifs over-represented in their non-coding sequences (upstream, introns, UTR), three of which correspond to previously known transcription factors: Zelda, Tramtrack and Trithorax-like (Trl). The nine discovered motifs were combined to scan ZGA-associated regions and predict about 1300 putative cis-regulatory modules. The fact that Trl is known to act as chromatin remodelling factor suggests that epigenetic regulation might play an important role in zygotic genome activation. We thus systematically compared the locations of predicted CRMs with ChIP-seq profiles for various transcription factors, 38 epigenetic marks from ModENCODE, and DNAse1 accessibility profiles. This analysis highlighted a strong and specific enrichment of predicted ZGA-associated CRMs for Zelda, CBP, Trl binding sites, as well as for histone marks associated with active enhancers (H3K4me1) and for open chromatin regions.

Conclusion

Based on the results of our computational analyses, we suggest a temporal model explaining the onset of zygotic genome activation by the combined action of transcription factors and epigenetic signals. Although this study is mainly based on the analysis of publicly available transcriptome and ChiP-seq datasets, the resulting model suggests novel mechanisms that underly the coordinated activation of several hundreds genes at a precise time point during embryonic development.

Chromatin-linked determinants of zygotic genome activation

The merging of the maternal and paternal genomes into a single pronucleus after fertilization is accompanied by a remarkable reconfiguration of chromatin in the newly formed zygote. The first stages of embryonic chromatin remodeling take place in the absence of ongoing transcription, during a species-specific developmental time-frame. Once post-fertilization chromatin states are organized, zygotic genome activation (ZGA) is initiated, and embryonic transcripts gradually take control of development. We review here transitions in chromatin modifications associated with the onset of ZGA, and the role of transcription factors and DNA motifs in the regulation of ZGA. We propose a model of sequential chromatin remodeling events preceding ZGA, leading to the onset of embryonic transcription.

Temporal dynamics, spatial range, and transcriptional interpretation of the Dorsal morphogen gradient

Dorsoventral pattern of Drosophila embryo is specified by the nuclear localization gradient of the transcription factor Dorsal. Genetic and genomic studies of this morphogen gradient provided important insights into spatial control of gene expression in development. Recent live imaging experiments revealed hitherto unappreciated dynamics of the Dorsal gradient and posed new questions about the mechanisms of its transcriptional interpretation. Some of these questions can be answered by models in which the morphogenetic capacity of the Dorsal gradient is potentiated by spatially uniform factors, such as Zelda, a transcription factor that plays a key role in the activation of zygotic transcription. Combinatorial effects of uniform and graded factors play an important role in the transcriptional and signaling cascades initiated by Dorsal and may explain differential positioning of gene expression borders by other morphogen gradients.

i-cisTarget: an integrative genomics method for the prediction of regulatory features and cis-regulatory modules

The field of regulatory genomics today is characterized by the generation of high-throughput data sets that capture genome-wide transcription factor (TF) binding, histone modifications, or DNAseI hypersensitive regions across many cell types and conditions. In this context, a critical question is how to make optimal use of these publicly available datasets when studying transcriptional regulation. Here, we address this question in Drosophila melanogaster for which a large number of high-throughput regulatory datasets are available. We developed i-cisTarget (where the 'i' stands for integrative), for the first time enabling the discovery of different types of enriched 'regulatory features' in a set of co-regulated sequences in one analysis, being either TF motifs or 'in vivo' chromatin features, or combinations thereof. We have validated our approach on 15 co-expressed gene sets, 21 ChIP data sets, 628 curated gene sets and multiple individual case studies, and show that meaningful regulatory features can be confidently discovered; that bona fide enhancers can be identified, both by in vivo events and by TF motifs; and that combinations of in vivo events and TF motifs further increase the performance of enhancer prediction.

Identification of early zygotic genes in the yellow fever mosquito Aedes aegypti and discovery of a motif involved in early zygotic genome activation

During early embryogenesis the zygotic genome is transcriptionally silent and all mRNAs present are of maternal origin. The maternal-zygotic transition marks the time over which embryogenesis changes its dependence from maternal RNAs to zygotically transcribed RNAs. Here we present the first systematic investigation of early zygotic genes (EZGs) in a mosquito species and focus on genes involved in the onset of transcription during 2–4 hr. We used transcriptome sequencing to identify the “pure” (without maternal expression) EZGs by analyzing transcripts from four embryonic time ranges of 0–2, 2–4, 4–8, and 8–12 hr, which includes the time of cellular blastoderm formation and up to the start of gastrulation. Blast of 16,789 annotated transcripts vs. the transcriptome reads revealed evidence for 63 (P<0.001) and 143 (P<0.05) nonmaternally derived transcripts having a significant increase in expression at 2–4 hr. One third of the 63 EZG transcripts do not have predicted introns compared to 10% of all Ae. aegypti genes. We have confirmed by RT-PCR that zygotic transcription starts as early as 2–3 hours. A degenerate motif VBRGGTA was found to be overrepresented in the upstream sequences of the identified EZGs using a motif identification software called SCOPE. We find evidence for homology between this motif and the TAGteam motif found in Drosophila that has been implicated in EZG activation. A 38 bp sequence in the proximal upstream sequence of a kinesin light chain EZG (KLC2.1) contains two copies of the mosquito motif. This sequence was shown to support EZG transcription by luciferase reporter assays performed on injected early embryos, and confers early zygotic activity to a heterologous promoter from a divergent mosquito species. The results of these studies are consistent with the model of early zygotic genome activation via transcriptional activators, similar to what has been found recently in Drosophila.

Pattern formation by graded and uniform signals in the early Drosophila embryo

The early Drosophila embryo is patterned by graded distributions of maternal transcription factors. Recent studies revealed that pattern formation by these graded signals depends on uniformly expressed transcriptional activators, such as Zelda. Removal of Zelda influences both the timing and the spatial expression domains for most of the genes controlled by maternal gradients. We demonstrate that some of these patterning defects, which range from temporal delay to loss of expression, can be rationalized with the use of a mathematical model based on cooperative binding of graded and uniform factors. This model makes a number of predictions, which we confirm experimentally by analyzing the expression of short gastrulation (sog), a gene that is controlled by a combination of the Dorsal morphogen gradient and Zelda. The proposed model suggests a general mechanism for the formation of nested gene expression domains, which is a hallmark of tissue patterning by morphogen gradients. According to this mechanism, the differential effects of a morphogen on its target genes can depend on their differential sensitivity to uniform factors.

Transcription co-activator SAYP mediates the action of STAT activator

Jak/STAT is an important signaling pathway mediating multiple events in development. We describe participation of metazoan co-activator SAYP/PHF10 in this pathway downstream of STAT. The latter, via its activation domain, interacts with the conserved core of SAYP. STAT is associated with the SAYP-containing co-activator complex BTFly and recruits BTFly onto genes. SAYP is necessary for stimulating STAT-driven transcription of numerous genes. Mutation of SAYP leads to maldevelopments similar to those observed in STAT mutants. Thus, SAYP is a novel co-activator mediating the action of STAT.

Zelda binding in the early Drosophila melanogaster embryo marks regions subsequently activated at the maternal-to-zygotic transition

The earliest stages of development in most metazoans are driven by maternally deposited proteins and mRNAs, with widespread transcriptional activation of the zygotic genome occurring hours after fertilization, at a period known as the maternal-to-zygotic transition (MZT). In Drosophila, the MZT is preceded by the transcription of a small number of genes that initiate sex determination, patterning, and other early developmental processes; and the zinc-finger protein Zelda (ZLD) plays a key role in their transcriptional activation. To better understand the mechanisms of ZLD activation and the range of its targets, we used chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-Seq) to map regions bound by ZLD before (mitotic cycle 8), during (mitotic cycle 13), and after (late mitotic cycle 14) the MZT. Although only a handful of genes are transcribed prior to mitotic cycle 10, we identified thousands of regions bound by ZLD in cycle 8 embryos, most of which remain bound through mitotic cycle 14. As expected, early ZLD-bound regions include the promoters and enhancers of genes transcribed at this early stage. However, we also observed ZLD bound at cycle 8 to the promoters of roughly a thousand genes whose first transcription does not occur until the MZT and to virtually all of the thousands of known and presumed enhancers bound at cycle 14 by transcription factors that regulate patterned gene activation during the MZT. The association between early ZLD binding and MZT activity is so strong that ZLD binding alone can be used to identify active promoters and regulatory sequences with high specificity and selectivity. This strong early association of ZLD with regions not active until the MZT suggests that ZLD is not only required for the earliest wave of transcription but also plays a major role in activating the genome at the MZT.

The newly fertilized eggs of most animal species begin development with a series of rapid cell divisions. During this time of rapid DNA replication, there is little or no transcription of the embryo's genome, with the synthesis of new proteins being directed by a store of maternally deposited mRNAs. Several hours after fertilization, at a period known as the maternal-to-zygotic transition (MZT), transcription of the embryo's genome begins in earnest, but little is known about how this process is initiated. In this paper we investigate the role of a protein known as Zelda (or ZLD) at the MZT in the laboratory model insect Drosophila melanogaster. ZLD had been previously shown to control the activation of a small number of genes expressed prior to the MZT. Here, using an experimental technique (ChIP-Seq) that allowed us to visualize where on the genome a protein is bound, we show that, approximately an hour prior to the MZT, ZLD is bound to most of the genomic regions active at the MZT. This suggests that ZLD may act as a kind of an “on switch” for the zygotic genome, poising regions where it binds for activation at the MZT, and this raises the possibility that similar master regulators of the MZT exist in other species.