Transcriptional and epigenetic signatures of zygotic genome activation during early drosophila embryogenesis
1 Technological Advances for Genomics and Clinics (TAGC), INSERM U1090, Université de la Méditerranée, Campus de Luminy, 13288 Marseille Cedex 9, France
2 Institut de Biologie du Développement de Marseille-Luminy (IBDML), Aix Marseille Université (AMU), UMR 7288 Case 907 - Parc Scientifique de Luminy, 13288 Marseille Cedex 9, France
3 Institut de Biologie de l’Ecole Normale Supérieure (IBENS) - UMR ENS and CNRS 8197 and INSERM 1024, 46 rue d’Ulm, 75005 Paris, France
4 Laboratoire de Bioinformatique des Génomes et des Réseaux (BiGRe), Université Libre de Bruxelles, Campus Plaine, CP 263, Bld du Triomphe, B-1050 Bruxelles, Belgium
BMC Genomics 2013, 14:226 doi:10.1186/1471-2164-14-226Published: 5 April 2013
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.
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.
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.