The TAGteam motif facilitates binding of 21 sequence-specific transcription factors in the Drosophila embryo

Zelda is dispensable for Drosophila melanogaster histone gene regulation

To ensure that the embryo can package exponentially increasing amounts of DNA, replication-dependent histones are some of the earliest transcribed genes from the zygotic genome. However, how the histone genes are identified is not known. The Drosophila melanogaster pioneer factor CLAMP regulates the embryonic histone genes and helps establish the histone locus body, a suite of factors that controls histone mRNA biosynthesis, but CLAMP is not unique to the histone genes. Zelda collaborates with CLAMP across the genome to regulate zygotic genome activation and target early activated genes. We hypothesized that Zelda helps identify histone genes for early embryonic expression. We found that Zelda targets the histone gene locus early during embryogenesis, prior to histone gene expression. However, depletion of zelda in the early embryo does not affect histone mRNA levels or prevent the recruitment of other factors. These results suggest the earliest events responsible for specifying the zygotic histone genes remain undiscovered.

HIRA and dPCIF1 coordinately establish totipotent chromatin and control orderly ZGA in Drosophila embryos

Early embryos undergo profound changes in their genomic architecture to establish the totipotent state, enabling pioneer factors to access chromatin and drive zygotic genome activation (ZGA). However, the mechanisms by which the totipotent state is established and properly interpreted by pioneer factors to allow orderly ZGA remain unknown. Here, we identify the H3.3-specific chaperone HIRA as a factor involving establishing totipotent-state chromatin in Drosophila early embryos. Through cophase separation with HIRA, the pioneer factor GAGA factor (GAF) efficiently binds to H3.3-marked nucleosomes to activate major-wave zygotic genes. Importantly, dPCIF1, a chromatin-associated protein, antagonized the GAF-HIRA interaction by competitively binding to HIRA, thereby restricting GAF on earlier chromatin and avoiding premature ZGA. Hence, the coordinated action of HIRA and dPCIF1 ensures sequential ZGA from the minor to major wave in early embryos. This study provides insights into understanding how a totipotent state is established and properly controlled during ZGA.

An oncogenic enhancer promotes melanoma progression via regulating ETV4 expression

BACKGROUND

Enhancers are important gene regulatory elements that promote the expression of critical genes in development and disease. Aberrant enhancer can modulate cancer risk and activate oncogenes that lead to the occurrence of various cancers. However, the underlying mechanism of most enhancers in cancer remains unclear. Here, we aim to explore the function and mechanism of a crucial enhancer in melanoma.

METHODS

Multi-omics data were applied to identify an enhancer (enh17) involved in melanoma progression. To evaluate the function of enh17, CRISPR/Cas9 technology were applied to knockout enh17 in melanoma cell line A375. RNA-seq, ChIP-seq and Hi-C data analysis integrated with luciferase reporter assay were performed to identify the potential target gene of enh17. Functional experiments were conducted to further validate the function of the target gene ETV4. Multi-omics data integrated with CUT&Tag sequencing were performed to validate the binding profile of the inferred transcription factor STAT3.

RESULTS

An enhancer, named enh17 here, was found to be aberrantly activated and involved in melanoma progression. CRISPR/Cas9-mediated deletion of enh17 inhibited cell proliferation, migration, and tumor growth of melanoma both in vitro and in vivo. Mechanistically, we identified ETV4 as a target gene regulated by enh17, and functional experiments further support ETV4 as a target gene that is involved in cancer-associated phenotypes. In addition, STAT3 acts as a transcription factor binding with enh17 to regulate the transcription of ETV4.

CONCLUSIONS

Our findings revealed that enh17 plays an oncogenic role and promotes tumor progression in melanoma, and its transcriptional regulatory mechanisms were fully elucidated, which may open a promising window for melanoma prevention and treatment.

Zelda is dispensable for Drosophila melanogaster histone gene regulation

To ensure that the embryo can package exponentially increasing amounts of DNA, replication-dependent histones are some of the earliest transcribed genes from the zygotic genome. However, how the histone genes are identified is not known. The pioneer factors Zelda and CLAMP collaborate at a subset of genes to regulate zygotic genome activation in Drosophila melanogaster and target early activated genes to induce transcription. CLAMP also regulates the embryonic histone genes and helps establish the histone locus body, a suite of factors that controls histone mRNA biosynthesis. The relationship between Zelda and CLAMP led us to hypothesize that Zelda helps identify histone genes for early embryonic expression. We found that Zelda targets the histone locus early during embryogenesis, prior to histone gene expression. However, depletion of zelda in the early embryo does not affect histone mRNA levels or histone locus body formation. While surprising, these results concur with other investigations into Zelda's role in the early embryo, suggesting the earliest factors responsible for specifying the zygotic histone genes remain undiscovered.

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.

Setting the stage for development: the maternal-to-zygotic transition in Drosophila

The zygote has a daunting task ahead of itself; it must develop from a single cell (fertilized egg) into a fully functioning adult with a multitude of different cell types. In the beginning, the zygote has help from its mother, in the form of gene products deposited into the egg, but eventually, it must rely on its own resources to proceed through development. The transfer of developmental control from the mother to the embryo is called the maternal-to-zygotic transition (MZT). All animals undergo this transition, which is defined by two main processes-the degradation of maternal RNAs and the synthesis of new RNAs from the zygote's own genome. Here, we review the regulation of the MZT in Drosophila, but given the broad conservation of this essential process, much of the regulation is shared among metazoans.

Transcriptional Enhancers in Drosophila

Key discoveries in Drosophila have shaped our understanding of cellular "enhancers." With a special focus on the fly, this chapter surveys properties of these adaptable cis-regulatory elements, whose actions are critical for the complex spatial/temporal transcriptional regulation of gene expression in metazoa. The powerful combination of genetics, molecular biology, and genomics available in Drosophila has provided an arena in which the developmental role of enhancers can be explored. Enhancers are characterized by diverse low- or high-throughput assays, which are challenging to interpret, as not all of these methods of identifying enhancers produce concordant results. As a model metazoan, the fly offers important advantages to comprehensive analysis of the central functions that enhancers play in gene expression, and their critical role in mediating the production of phenotypes from genotype and environmental inputs. A major challenge moving forward will be obtaining a quantitative understanding of how these cis-regulatory elements operate in development and disease.

Bicoid-Dependent Activation of the Target Gene hunchback Requires a Two-Motif Sequence Code in a Specific Basal Promoter

In complex genetic loci, individual enhancers interact most often with specific basal promoters. Here we investigate the activation of the Bicoid target gene hunchback (hb), which contains two basal promoters (P1 and P2). Early in embryogenesis, P1 is silent, while P2 is strongly activated. In vivo deletion of P2 does not cause activation of P1, suggesting that P2 contains intrinsic sequence motifs required for activation. We show that a two-motif code (a Zelda binding site plus TATA) is required and sufficient for P2 activation. Zelda sites are present in the promoters of many embryonically expressed genes, but the combination of Zelda plus TATA does not seem to be a general code for early activation or Bicoid-specific activation per se. Because Zelda sites are also found in Bicoid-dependent enhancers, we propose that simultaneous binding to both enhancers and promoters independently synchronizes chromatin accessibility and facilitates correct enhancer-promoter interactions.

ATAC-seq reveals regional differences in enhancer accessibility during the establishment of spatial coordinates in the Drosophila blastoderm

Establishment of spatial coordinates during Drosophila embryogenesis relies on differential regulatory activity of axis patterning enhancers. Concentration gradients of activator and repressor transcription factors (TFs) provide positional information to each enhancer, which in turn promotes transcription of a target gene in a specific spatial pattern. However, the interplay between an enhancer regulatory activity and its accessibility as determined by local chromatin organization is not well understood. We profiled chromatin accessibility with ATAC-seq in narrow, genetically tagged domains along the antero-posterior axis in the Drosophila blastoderm. We demonstrate that one-quarter of the accessible genome displays significant regional variation in its ATAC-seq signal immediately after zygotic genome activation. Axis patterning enhancers are enriched among the most variable intervals, and their accessibility changes correlate with their regulatory activity. In an embryonic domain where an enhancer receives a net activating TF input and promotes transcription, it displays elevated accessibility in comparison to a domain where it receives a net repressive input. We propose that differential accessibility is a signature of patterning cis-regulatory elements in the Drosophila blastoderm and discuss potential mechanisms by which accessibility of enhancers may be modulated by activator and repressor TFs.

Pure early zygotic genes in the Asian malaria mosquito Anopheles stephensi

BACKGROUND

The Asian malaria mosquito, Anopheles stephensi, is a major urban malaria vector in the Middle East and on the Indian subcontinent. Early zygotic transcription, which marks the maternal-to-zygotic transition, has not been systematically studied in An. stephensi or any other Anopheles mosquitoes. Improved understanding of early embryonic gene expression in An. stephensi will facilitate genetic and evolutionary studies and help with the development of novel control strategies for this important disease vector.

RESULTS

We obtained RNA-seq data in biological triplicates from four early An. stephensi embryonic time points. Using these data, we identified 70 and 153 pure early zygotic genes (pEZGs) under stringent and relaxed conditions, respectively. We show that these pEZGs are enriched in functional groups related to DNA-binding transcription regulators, cell cycle modulators, proteases, transport, and cellular metabolism. On average these pEZGs are shorter and have less introns than other An. stephensi genes. Some of the pEZGs may arise de novo while others have clear non-pEZG paralogs. There is no or very limited overlap between An. stephensi pEZGs and Drosophila melanogaster or Aedes aegypti pEZGs. Interestingly, the upstream region of An. stephensi pEZGs lack significant enrichment of a previously reported TAGteam/VBRGGTA motif found in the regulatory region of pEZGs in D. melanogaster and Ae. aegypti. However, a GT-rich motif was found in An. stephensi pEZGs instead.

CONCLUSIONS

We have identified a number of pEZGs whose predicted functions and structures are consistent with their collective roles in the degradation of maternally deposited components, activation of the zygotic genome, cell division, and metabolism. The pEZGs appear to rapidly turn over within the Dipteran order and even within the Culicidae family. These pEZGs, and the shared regulatory motif, could provide the promoter or regulatory sequences to drive gene expression in the syncytial or early cellular blastoderm, a period when the developing embryo is accessible to genetic manipulation. In addition, these molecular resources may be used to achieve sex separation of mosquitoes for sterile insect technique.

The dynamics of energy metabolism in the tick embryo

The cattle tick Rhipicephalus (Boophilus) microplus is an ectoparasite capable of transmitting a large number of pathogens, causing considerable losses in the cattle industry, with substantial damage to livestock. Over the years, important stages of its life cycle, such as the embryo, have been largely ignored by researchers. Tick embryogenesis has been typically described as an energy-consuming process, sustaining cell proliferation, differentiation, and growth. During the embryonic stage of arthropods, there is mobilization of metabolites of maternal origin for the development of organs and tissues of the embryo. Glycogen resynthesis in late embryogenesis is considered as an effective indicator of embryonic integrity. In the cattle tick R.(B. (B.) microplus, glycogen resynthesis is sustained by protein degradation through the gluconeogenesis pathway at the end of the embryonic period. Despite recent advancements in research on tick energy metabolism at the molecular level, the dynamics of nutrient utilization during R. (B.) microplus embryogenesis is still poorly understood. The present review aims to describe the regulatory mechanisms of carbohydrate metabolism during maternal-zygotic transition and identify possible new targets for the development of novel drugs and other control measures against R. (B.) microplus infestations.

TFIIH localization is highly dynamic during zygotic genome activation in Drosophila, and its depletion causes catastrophic mitosis

In Drosophila, zygotic genome activation occurs in pre-blastoderm embryos during rapid mitotic divisions. How the transcription machinery is coordinated to achieve this goal in a very brief time span is still poorly understood. Transcription factor II H (TFIIH) is fundamental for transcription initiation by RNA polymerase II (RNAPII). Herein, we show the in vivo dynamics of TFIIH at the onset of transcription in Drosophila embryos. TFIIH shows an oscillatory behaviour between the nucleus and cytoplasm. TFIIH foci are observed from interphase to metaphase, and colocalize with those for RNAPII phosphorylated at serine 5 (RNAPIIS5P) at prophase, suggesting that transcription occurs during the first mitotic phases. Furthermore, embryos with defects in subunits of either the CAK or the core subcomplexes of TFIIH show catastrophic mitosis. Although, transcriptome analyses show altered expression of several maternal genes that participate in mitosis, the global level of RNAPIIS5P in TFIIH mutant embryos is similar to that in the wild type, therefore, a direct role for TFIIH in mitosis cannot be ruled out. These results provide important insights regarding the role of a basal transcription machinery component when the zygotic genome is activated.

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.

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.

Genome-wide search for Zelda-like chromatin signatures identifies GAF as a pioneer factor in early fly development

BACKGROUND

The protein Zelda was shown to play a key role in early Drosophila development, binding thousands of promoters and enhancers prior to maternal-to-zygotic transition (MZT), and marking them for transcriptional activation. Recently, we showed that Zelda acts through specific chromatin patterns of histone modifications to mark developmental enhancers and active promoters. Intriguingly, some Zelda sites still maintain these chromatin patterns in Drosophila embryos lacking maternal Zelda protein. This suggests that additional Zelda-like pioneer factors may act in early fly embryos.

RESULTS

We developed a computational method to analyze and refine the chromatin landscape surrounding early Zelda peaks, using a multichannel spectral clustering. This allowed us to characterize their chromatin patterns through MZT (mitotic cycles 8-14). Specifically, we focused on H3K4me1, H3K4me3, H3K18ac, H3K27ac, and H3K27me3 and identified three different classes of chromatin signatures, matching "promoters," "enhancers" and "transiently bound" Zelda peaks. We then further scanned the genome using these chromatin patterns and identified additional loci-with no Zelda binding-that show similar chromatin patterns, resulting with hundreds of Zelda-independent putative enhancers. These regions were found to be enriched with GAGA factor (GAF, Trl) and are typically located near early developmental zygotic genes. Overall our analysis suggests that GAF, together with Zelda, plays an important role in activating the zygotic genome.

CONCLUSIONS

As we show, our computational approach offers an efficient algorithm for characterizing chromatin signatures around some loci of interest and allows a genome-wide identification of additional loci with similar chromatin patterns.

Regulation of DNA Replication in Early Embryonic Cleavages

Early embryonic cleavages are characterized by short and highly synchronous cell cycles made of alternating S- and M-phases with virtually absent gap phases. In this contracted cell cycle, the duration of DNA synthesis can be extraordinarily short. Depending on the organism, the whole genome of an embryo is replicated at a speed that is between 20 to 60 times faster than that of a somatic cell. Because transcription in the early embryo is repressed, DNA synthesis relies on a large stockpile of maternally supplied proteins stored in the egg representing most, if not all, cellular genes. In addition, in early embryonic cell cycles, both replication and DNA damage checkpoints are inefficient. In this article, we will review current knowledge on how DNA synthesis is regulated in early embryos and discuss possible consequences of replicating chromosomes with little or no quality control.

Genome-wide measurement of spatial expression in patterning mutants of Drosophila melanogaster

Patterning in the Drosophila melanogaster embryo is affected by multiple maternal factors, but the effect of these factors on spatial gene expression has not been systematically analyzed. Here we characterize the effect of the maternal factors Zelda, Hunchback and Bicoid by cryosectioning wildtype and mutant blastoderm stage embryos and sequencing mRNA from each slice. The resulting atlas of spatial gene expression highlights the intersecting roles of these factors in regulating spatial patterns, and serves as a resource for researchers studying spatial patterning in the early embryo. We identify a large number of genes with both expected and unexpected patterning changes, and through integrated analysis of transcription factor binding data identify common themes in genes with complex dependence on these transcription factors.

Regulation of Zygotic Genome and Cellular Pluripotency

Events, manifesting transition from maternal to zygotic period of development are studied for more than 100 years, but underlying mechanisms are not yet clear. We provide a brief historical overview of development of concepts and explain the specific terminology used in the field. We further discuss differences and similarities between the zygotic genome activation and in vitro reprogramming process. Finally, we envision the future research directions within the field, where biochemical methods will play increasingly important role.

Zygotic Genome Activators, Developmental Timing, and Pluripotency

The transcription factors Pou5f1, Sox2, and Nanog are central regulators of pluripotency in mammalian ES and iPS cells. In vertebrate embryos, Pou5f1/3, SoxB1, and Nanog control zygotic genome activation and participate in lineage decisions. We review the current knowledge of the roles of these genes in developing vertebrate embryos from fish to mammals and suggest a model for pluripotency gene regulatory network functions in early development.

Enhancer divergence and cis-regulatory evolution in the human and chimp neural crest

cis-regulatory changes play a central role in morphological divergence, yet the regulatory principles underlying emergence of human traits remain poorly understood. Here, we use epigenomic profiling from human and chimpanzee cranial neural crest cells to systematically and quantitatively annotate divergence of craniofacial cis-regulatory landscapes. Epigenomic divergence is often attributable to genetic variation within TF motifs at orthologous enhancers, with a novel motif being most predictive of activity biases. We explore properties of this cis-regulatory change, revealing the role of particular retroelements, uncovering broad clusters of species-biased enhancers near genes associated with human facial variation, and demonstrating that cis-regulatory divergence is linked to quantitative expression differences of crucial neural crest regulators. Our work provides a wealth of candidates for future evolutionary studies and demonstrates the value of "cellular anthropology," a strategy of using in-vitro-derived embryonic cell types to elucidate both fundamental and evolving mechanisms underlying morphological variation in higher primates.

Zelda overcomes the high intrinsic nucleosome barrier at enhancers during Drosophila zygotic genome activation

The Drosophila genome activator Vielfaltig (Vfl), also known as Zelda (Zld), is thought to prime enhancers for activation by patterning transcription factors (TFs). Such priming is accompanied by increased chromatin accessibility, but the mechanisms by which this occurs are poorly understood. Here, we analyze the effect of Zld on genome-wide nucleosome occupancy and binding of the patterning TF Dorsal (Dl). Our results show that early enhancers are characterized by an intrinsically high nucleosome barrier. Zld tackles this nucleosome barrier through local depletion of nucleosomes with the effect being dependent on the number and position of Zld motifs. Without Zld, Dl binding decreases at enhancers and redistributes to open regions devoid of enhancer activity. We propose that Zld primes enhancers by lowering the high nucleosome barrier just enough to assist TFs in accessing their binding motifs and promoting spatially controlled enhancer activation if the right patterning TFs are present. We envision that genome activators in general will utilize this mechanism to activate the zygotic genome in a robust and precise manner.

Zelda is differentially required for chromatin accessibility, transcription factor binding, and gene expression in the early Drosophila embryo

The transition from a specified germ cell to a population of pluripotent cells occurs rapidly following fertilization. During this developmental transition, the zygotic genome is largely transcriptionally quiescent and undergoes significant chromatin remodeling. In Drosophila, the DNA-binding protein Zelda (also known as Vielfaltig) is required for this transition and for transcriptional activation of the zygotic genome. Open chromatin is associated with Zelda-bound loci, as well as more generally with regions of active transcription. Nonetheless, the extent to which Zelda influences chromatin accessibility across the genome is largely unknown. Here we used formaldehyde-assisted isolation of regulatory elements to determine the role of Zelda in regulating regions of open chromatin in the early embryo. We demonstrate that Zelda is essential for hundreds of regions of open chromatin. This Zelda-mediated chromatin accessibility facilitates transcription-factor recruitment and early gene expression. Thus, Zelda possesses some key characteristics of a pioneer factor. Unexpectedly, chromatin at a large subset of Zelda-bound regions remains open even in the absence of Zelda. The GAGA factor-binding motif and embryonic GAGA factor binding are specifically enriched in these regions. We propose that both Zelda and GAGA factor function to specify sites of open chromatin and together facilitate the remodeling of the early embryonic genome.

How gene expression in fast-proliferating cells keeps pace

The development of living organisms requires a precise coordination of all basic cellular processes, in space and time. Early embryogenesis of most species with externally deposited eggs starts with a series of extremely fast cleavage cycles. These divisions have a strong influence on gene expression as mitosis represses transcription and pre-mRNA processing. In this review, we will describe the distinct adaptations for efficient gene expression and discuss the emerging role of the multifunctional NineTeen Complex (NTC) in gene expression and genomic stability during fast proliferation.

BMPs regulate msx gene expression in the dorsal neuroectoderm of Drosophila and vertebrates by distinct mechanisms

In a broad variety of bilaterian species the trunk central nervous system (CNS) derives from three primary rows of neuroblasts. The fates of these neural progenitor cells are determined in part by three conserved transcription factors: vnd/nkx2.2, ind/gsh and msh/msx in Drosophila melanogaster/vertebrates, which are expressed in corresponding non-overlapping patterns along the dorsal-ventral axis. While this conserved suite of “neural identity” gene expression strongly suggests a common ancestral origin for the patterning systems, it is unclear whether the original regulatory mechanisms establishing these patterns have been similarly conserved during evolution. In Drosophila, genetic evidence suggests that Bone Morphogenetic Proteins (BMPs) act in a dosage-dependent fashion to repress expression of neural identity genes. BMPs also play a dose-dependent role in patterning the dorsal and lateral regions of the vertebrate CNS, however, the mechanism by which they achieve such patterning has not yet been clearly established. In this report, we examine the mechanisms by which BMPs act on cis-regulatory modules (CRMs) that control localized expression of the Drosophila msh and zebrafish (Danio rerio) msxB in the dorsal central nervous system (CNS). Our analysis suggests that BMPs act differently in these organisms to regulate similar patterns of gene expression in the neuroectoderm: repressing msh expression in Drosophila, while activating msxB expression in the zebrafish. These findings suggest that the mechanisms by which the BMP gradient patterns the dorsal neuroectoderm have reversed since the divergence of these two ancient lineages.

The trunk nervous system of both vertebrates and invertebrates develops from three primary rows of neural stem cells whose fate is determined by neural identity genes expressed in an evolutionarily conserved dorso-ventral pattern. Establishment of this pattern requires a shared signaling pathway in both groups of animals. Previous studies suggested that a shared signaling pathway functions in opposite ways in vertebrates and invertebrates, despite the final patterning outcomes having remained the same. Here, we employ bioinformatics, biochemistry, and transgenic animal technology to elucidate the genetic mechanism by which this pathway can engage the same components to generate opposite instructions and yet arrive at similar outcomes in patterning of the nervous system. Our findings highlight how natural selection can act to conserve a particular output pattern despite changes during evolution in the genetic mechanisms underlying the formation of this pattern.

From egg to gastrula: how the cell cycle is remodeled during the Drosophila mid-blastula transition

Many, if not most, embryos begin development with extremely short cell cycles that exhibit unusually rapid DNA replication and no gap phases. The commitment to the cell cycle in the early embryo appears to preclude many other cellular processes that only emerge as the cell cycle slows just prior to gastrulation at a major embryonic transition known as the mid-blastula transition (MBT). As reviewed here, genetic and molecular studies in Drosophila have identified changes that extend S phase and introduce a post-replicative gap phase, G2, to slow the cell cycle. Although many mysteries remain about the upstream regulators of these changes, we review the core mechanisms of the change in cell cycle regulation and discuss advances in our understanding of how these might be timed and triggered. Finally, we consider how the elements of this program may be conserved or changed in other organisms.

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.

Hit-and-run transcriptional control by bZIP1 mediates rapid nutrient signaling in Arabidopsis

The dynamic nature of gene regulatory networks allows cells to rapidly respond to environmental change. However, the underlying temporal connections are missed, even in kinetic studies, as transcription factor (TF) binding within at least one time point is required to identify primary targets. The TF-regulated but unbound genes are dismissed as secondary targets. Instead, we report that these genes comprise transient TF-target interactions most relevant to rapid signal transduction. We temporally perturbed a master TF (Basic Leucine Zipper 1, bZIP1) and the nitrogen (N) signal it transduces and integrated TF regulation and binding data from the same cell samples. Our enabling approach could identify primary TF targets based solely on gene regulation, in the absence of TF binding. We uncovered three classes of primary TF targets: (i) poised (TF-bound but not TF-regulated), (ii) stable (TF-bound and TF-regulated), and (iii) transient (TF-regulated but not TF-bound), the largest class. Unexpectedly, the transient bZIP1 targets are uniquely relevant to rapid N signaling in planta, enriched in dynamic N-responsive genes, and regulated by TF and N signal interactions. These transient targets include early N responders nitrate transporter 2.1 and NIN-like protein 3, bound by bZIP1 at 1-5 min, but not at later time points following TF perturbation. Moreover, promoters of these transient targets are uniquely enriched with cis-regulatory motifs coinherited with bZIP1 binding sites, suggesting a recruitment role for bZIP1. This transient mode of TF action supports a classic, but forgotten, "hit-and-run" transcription model, which enables a "catalyst TF" to activate a large set of targets within minutes of signal perturbation.

Zelda potentiates morphogen activity by increasing chromatin accessibility

Zygotic genome activation (ZGA) is a major genome programming event whereby the cells of the embryo begin to adopt specified fates. Experiments in Drosophila and zebrafish have revealed that ZGA depends on transcription factors that provide large-scale control of gene expression by direct and specific binding to gene regulatory sequences. Zelda (Zld) plays such a role in the Drosophila embryo, where it has been shown to control the action of patterning signals; however, the mechanisms underlying this effect remain largely unclear. A recent model proposed that Zld binding sites act as quantitative regulators of the spatiotemporal expression of genes activated by Dorsal (Dl), the morphogen that patterns the dorsoventral axis. Here we tested this model experimentally, using enhancers of brinker (brk) and short gastrulation (sog), both of which are directly activated by Dl, but at different concentration thresholds. In agreement with the model, we show that there is a clear positive correlation between the number of Zld binding sites and the spatial domain of enhancer activity. Likewise, the timing of expression could be advanced or delayed. We present evidence that Zld facilitates binding of Dl to regulatory DNA, and that this is associated with increased chromatin accessibility. Importantly, the change in chromatin accessibility is strongly correlated with the change in Zld binding, but not Dl. We propose that the ability of genome activators to facilitate readout of transcriptional input is key to widespread transcriptional induction during ZGA.

Impacts of the ubiquitous factor Zelda on Bicoid-dependent DNA binding and transcription in Drosophila

The Drosophila transcription factor Bicoid (Bcd) binds thousands of genomic sites during early embryogenesis, but it is unclear how many of these binding events are functionally important. Here, Small and colleagues test the role of the maternal factor Zelda (Zld) in Bcd-mediated binding and transcription. Embryos lacking Zld show enhanced Bcd binding to a subset of genomic locations, causing early activation of target genes normally silent until later stages. This study demonstrates a critical role for Zld in controlling Bcd binding and target gene activation in the early embryo.

In vivo cross-linking studies suggest that the Drosophila transcription factor Bicoid (Bcd) binds to several thousand sites during early embryogenesis, but it is not clear how many of these binding events are functionally important. In contrast, reporter gene studies have identified >60 Bcd-dependent enhancers, all of which contain clusters of the consensus binding sequence TAATCC. These studies also identified clusters of TAATCC motifs (inactive fragments) that failed to drive Bcd-dependent activation. In general, active fragments showed higher levels of Bcd binding in vivo and were enriched in predicted binding sites for the ubiquitous maternal protein Zelda (Zld). Here we tested the role of Zld in Bcd-mediated binding and transcription. Removal of Zld function and mutations in Zld sites caused significant reductions in Bcd binding to known enhancers and variable effects on the activation and spatial positioning of Bcd-dependent expression patterns. Also, insertion of Zld sites converted one of six inactive fragments into a Bcd-responsive enhancer. Genome-wide binding experiments in zld mutants showed variable effects on Bcd-binding peaks, ranging from strong reductions to significantly enhanced levels of binding. Increases in Bcd binding caused the precocious Bcd-dependent activation of genes that are normally not expressed in early embryos, suggesting that Zld controls the genome-wide binding profile of Bcd at the qualitative level and is critical for selecting target genes for activation in the early embryo. These results underscore the importance of combinatorial binding in enhancer function and provide data that will help predict regulatory activities based on DNA sequence.

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.

Evolution of the sex-lethal gene in insects and origin of the sex-determination system in Drosophila

Sex-lethal (Sxl) functions as the switch gene for sex-determination in Drosophila melanogaster by engaging a regulatory cascade. Thus far the origin and evolution of both the regulatory system and SXL protein's sex-determination function have remained largely unknown. In this study, we explore systematically the Sxl homologs in a wide range of insects, including the 12 sequenced Drosophila species, medfly, blowflies, housefly, Megaselia scalaris, mosquitoes, butterfly, beetle, honeybee, ant, and aphid. We find that both the male-specific and embryo-specific exons exist in all Drosophila species. The homologous male-specific exon is also present in Scaptodrosophila lebanonensis, but it does not have in-frame stop codons, suggesting the exon's functional divergence between Drosophila and Scaptodrosophila after acquiring it in their common ancestor. Two motifs closely related to the exons' functions, the SXL binding site poly(U) and the transcription-activating motif TAGteam, surprisingly exhibit broader phylogenetic distributions than the exons. Some previously unknown motifs that are restricted to or more abundant in Drosophila and S. lebanonensis than in other insects are also identified. Finally, phylogenetic analysis suggests that the SXL's novel sex-determination function in Drosophila is more likely attributed to the changes in the N- and C-termini rather than in the RNA-binding region. Thus, our results provide a clearer picture of the phylogeny of the Sxl's cis-regulatory elements and protein sequence changes, and so lead to a better understanding of the origin of sex-determination in Drosophila and also raise some new questions regarding the evolution of Sxl.

Identification of transcription factor binding sites from ChIP-seq data at high resolution

MOTIVATION

Chromatin immunoprecipitation coupled to next-generation sequencing (ChIP-seq) is widely used to study the in vivo binding sites of transcription factors (TFs) and their regulatory targets. Recent improvements to ChIP-seq, such as increased resolution, promise deeper insights into transcriptional regulation, yet require novel computational tools to fully leverage their advantages.

RESULTS

To this aim, we have developed peakzilla, which can identify closely spaced TF binding sites at high resolution (i.e. resolves individual binding sites even if spaced closely), as we demonstrate using semisynthetic datasets, performing ChIP-seq for the TF Twist in Drosophila embryos with different experimental fragment sizes, and analyzing ChIP-exo datasets. We show that the increased resolution reached by peakzilla is highly relevant, as closely spaced Twist binding sites are strongly enriched in transcriptional enhancers, suggesting a signature to discriminate functional from abundant non-functional or neutral TF binding. Peakzilla is easy to use, as it estimates all the necessary parameters from the data and is freely available.

AVAILABILITY AND IMPLEMENTATION

The peakzilla program is available from https://github.com/steinmann/peakzilla or http://www.starklab.org/data/peakzilla/.

CONTACT

stark@starklab.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pou5f1 transcription factor controls zygotic gene activation in vertebrates

The development of multicellular animals is initially controlled by maternal gene products deposited in the oocyte. During the maternal-to-zygotic transition, transcription of zygotic genes commences, and developmental control starts to be regulated by zygotic gene products. In Drosophila, the transcription factor Zelda specifically binds to promoters of the earliest zygotic genes and primes them for activation. It is unknown whether a similar regulation exists in other animals. We found that zebrafish Pou5f1, a homolog of the mammalian pluripotency transcription factor Oct4, occupies SOX-POU binding sites before the onset of zygotic transcription and activates the earliest zygotic genes. Our data position Pou5f1 and SOX-POU sites at the center of the zygotic gene activation network of vertebrates and provide a link between zygotic gene activation and pluripotency control.

Computational identification of diverse mechanisms underlying transcription factor-DNA occupancy

ChIP-based genome-wide assays of transcription factor (TF) occupancy have emerged as a powerful, high-throughput method to understand transcriptional regulation, especially on a global scale. This has led to great interest in the underlying biochemical mechanisms that direct TF-DNA binding, with the ultimate goal of computationally predicting a TF's occupancy profile in any cellular condition. In this study, we examined the influence of various potential determinants of TF-DNA binding on a much larger scale than previously undertaken. We used a thermodynamics-based model of TF-DNA binding, called “STAP,” to analyze 45 TF-ChIP data sets from Drosophila embryonic development. We built a cross-validation framework that compares a baseline model, based on the ChIP'ed (“primary”) TF's motif, to more complex models where binding by secondary TFs is hypothesized to influence the primary TF's occupancy. Candidates interacting TFs were chosen based on RNA-SEQ expression data from the time point of the ChIP experiment. We found widespread evidence of both cooperative and antagonistic effects by secondary TFs, and explicitly quantified these effects. We were able to identify multiple classes of interactions, including (1) long-range interactions between primary and secondary motifs (separated by ≤150 bp), suggestive of indirect effects such as chromatin remodeling, (2) short-range interactions with specific inter-site spacing biases, suggestive of direct physical interactions, and (3) overlapping binding sites suggesting competitive binding. Furthermore, by factoring out the previously reported strong correlation between TF occupancy and DNA accessibility, we were able to categorize the effects into those that are likely to be mediated by the secondary TF's effect on local accessibility and those that utilize accessibility-independent mechanisms. Finally, we conducted in vitro pull-down assays to test model-based predictions of short-range cooperative interactions, and found that seven of the eight TF pairs tested physically interact and that some of these interactions mediate cooperative binding to DNA.

Chromatin Immunoprecipitation (ChIP)-based genome-wide assays of transcription factor (TF) occupancy have emerged as a powerful, high throughput method to understand transcriptional regulation, especially on a global scale. Here, we utilize 45 ChIP-chip and ChIP-SEQ data sets from Drosophila to explore the underlying mechanisms of TF-DNA binding. For this, we employ a biophysically motivated computational model, in conjunction with over 300 TF motifs (binding specificities) as well as gene expression and DNA accessibility data from different developmental stages in Drosophila embryos. Our findings provide robust statistical evidence of the role played by TF-TF interactions in shaping genome-wide TF-DNA binding profiles, and thus in directing gene regulation. Our method allows us to go beyond simply recognizing the existence of such interactions, to quantifying their effects on TF occupancy. We are able to categorize the probable mechanisms of these effects as involving direct physical interactions versus accessibility-mediated indirect interactions, long-range versus short-range interactions, and cooperative versus antagonistic interactions. Our analysis reveals widespread evidence of combinatorial regulation present in recently generated ChIP data sets, and sets the stage for rich integrative models of the future that will predict cell type-specific TF occupancy values from sequence and expression data.

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.

Interpreting the regulatory genome: the genomics of transcription factor function in Drosophila melanogaster

Researchers have now had access to the fully sequenced Drosophila melanogaster genome for over a decade, and the sequenced genomes of 11 additional Drosophila species have been available for almost 5 years, with more species' genomes becoming available every year [Adams MD, Celniker SE, Holt RA, et al. The genome sequence of Drosophila melanogaster. Science 2000;287:2185-95; Clark AG, Eisen MB, Smith DR, et al. Evolution of genes and genomes on the Drosophila phylogeny. Nature 2007;450:203-18]. Although the best studied of the D. melanogaster transcription factors (TFs) were cloned before sequencing of the genome, the availability of sequence data promised to transform our understanding of TFs and gene regulatory networks. Sequenced genomes have allowed researchers to generate tools for high-throughput characterization of gene expression levels, genome-wide TF localization and analyses of evolutionary constraints on DNA elements across multiple species. With an estimated 700 DNA-binding proteins in the Drosophila genome, it will be many years before each potential sequence-specific TF is studied in detail, yet the last decade of functional genomics research has already impacted our view of gene regulatory networks and TF DNA recognition.

Disentangling the many layers of eukaryotic transcriptional regulation

Regulation of gene expression in eukaryotes is an extremely complex process. In this review, we break down several critical steps, emphasizing new data and techniques that have expanded current gene regulatory models. We begin at the level of DNA sequence where cis-regulatory modules (CRMs) provide important regulatory information in the form of transcription factor (TF) binding sites. In this respect, CRMs function as instructional platforms for the assembly of gene regulatory complexes. We discuss multiple mechanisms controlling complex assembly, including cooperative DNA binding, combinatorial codes, and CRM architecture. The second section of this review places CRM assembly in the context of nucleosomes and condensed chromatin. We discuss how DNA accessibility and histone modifications contribute to TF function. Lastly, new advances in chromosomal mapping techniques have provided increased understanding of intra- and interchromosomal interactions. We discuss how these topological maps influence gene regulatory models.

HOT DNAs: a novel class of developmental enhancers

Enhancers mediate localized patterns of gene expression during development. A common feature of "traditional" enhancers is the presence of clustered binding motifs for sequence-specific transcription factors (TFs). In this issue of Genes & Development, Kvon and colleagues (pp. 908-913) present new evidence that HOT (highly occupied transcription) DNAs direct specific patterns of gene expression, despite being depleted for TF-binding motifs.

HOT regions function as patterned developmental enhancers and have a distinct cis-regulatory signature

HOT (highly occupied target) regions bound by many transcription factors are considered to be one of the most intriguing findings of the recent modENCODE reports, yet their functions have remained unclear. We tested 108 Drosophila melanogaster HOT regions in transgenic embryos with site-specifically integrated transcriptional reporters. In contrast to prior expectations, we found 102 (94%) to be active enhancers during embryogenesis and to display diverse spatial and temporal patterns, reminiscent of expression patterns for important developmental genes. Remarkably, HOT regions strongly activate nearby genes and are required for endogenous gene expression, as we show using bacterial artificial chromosome (BAC) transgenesis. HOT enhancers have a distinct cis-regulatory signature with enriched sequence motifs for the global activators Vielfaltig, also known as Zelda, and Trithorax-like, also known as GAGA. This signature allows the prediction of HOT versus control regions from the DNA sequence alone.