Canonical nucleosome organization at promoters forms during genome activation

H3K4me2 distinguishes a distinct class of enhancers during the maternal-to-zygotic transition

After egg fertilization, an initially silent embryonic genome is transcriptionally activated during the maternal-to-zygotic transition. In zebrafish, maternal vertebrate pluripotency factors Nanog, Pou5f3 (OCT4 homolog), and Sox19b (SOX2 homolog) (NPS) play essential roles in orchestrating embryonic genome activation, acting as "pioneers" that open condensed chromatin and mediate acquisition of activating histone modifications. However, some embryonic gene transcription still occurs in the absence of these factors, suggesting the existence of other mechanisms regulating genome activation. To identify chromatin signatures of these unknown pathways, we profiled the histone modification landscape of zebrafish embryos using CUT&RUN. Our regulatory map revealed two subclasses of enhancers distinguished by presence or absence of H3K4me2. Enhancers lacking H3K4me2 tend to require NPS factors for de novo activation, while enhancers bearing H3K4me2 are epigenetically bookmarked by DNA hypomethylation to recapitulate gamete activity in the embryo, independent of NPS pioneering. Thus, parallel enhancer activation pathways combine to induce transcriptional reprogramming to pluripotency in the early embryo.

The developmental and evolutionary characteristics of transcription factor binding site clustered regions based on an explainable machine learning model

Gene expression is temporally and spatially regulated by the interaction of transcription factors (TFs) and cis-regulatory elements (CREs). The uneven distribution of TF binding sites across the genome poses challenges in understanding how this distribution evolves to regulate spatio-temporal gene expression and consequent heritable phenotypic variation. In this study, chromatin accessibility profiles and gene expression profiles were collected from several species including mammals (human, mouse, bovine), fish (zebrafish and medaka), and chicken. Transcription factor binding sites clustered regions (TFCRs) at different embryonic stages were characterized to investigate regulatory evolution. The study revealed dynamic changes in TFCR distribution during embryonic development and species evolution. The synchronization between TFCR complexity and gene expression was assessed across species using RegulatoryScore. Additionally, an explainable machine learning model highlighted the importance of the distance between TFCR and promoter in the coordinated regulation of TFCRs on gene expression. Our results revealed the developmental and evolutionary dynamics of TFCRs during embryonic development from fish, chicken to mammals. These data provide valuable resources for exploring the relationship between transcriptional regulation and phenotypic differences during embryonic development.

Chromatin conformation and histone modification profiling across human kidney anatomic regions

The three major anatomic regions of the human kidney include the cortex, medulla and papilla, with different functions and vulnerabilities to kidney diseases. Epigenetic mechanisms underlying these anatomic structures are incompletely understood. Here, we performed chromatin conformation capture with Hi-C and histone modification H3K4me3/H3K27me3 Cleavage Under Targets and Release Using Nuclease (CUT&RUN) sequencing on the kidney cortex, medulla and papilla dissected from one individual donor. Nuclear suspensions were generated from each region and split subjected to paired Hi-C and CUT&RUN sequencing. We evaluated the quality of next-generation sequencing data, Hi-C chromatin contact matrices and CUT&RUN peak calling. H3K4me3 and H3K27me3 histone modifications represent active and repressive gene transcription, respectively, and differences in chromatin conformation between kidney regions can be analyzed with this dataset. All raw and processed data files are publicly available, allowing researchers to survey the epigenetic landscape across regional human kidney anatomy.

Local nuclear to cytoplasmic ratio regulates chaperone-dependent H3 variant incorporation during zygotic genome activation

Early embryos often have relatively unstructured chromatin that lacks active and inactive domains typical of differentiated cells. In many species, these regulatory domains are established during zygotic genome activation (ZGA). In Drosophila, ZGA occurs after 13 fast, reductive, syncytial nuclear divisions during which the nuclear to cytoplasmic (N/C) ratio grows exponentially. These divisions incorporate maternally-loaded, cytoplasmic pools of histones into chromatin. Previous work found that chromatin incorporation of replication-coupled histone H3 decreases while its variant H3.3 increases in the cell cycles leading up to ZGA. In other cell types, H3.3 is associated with sites of active transcription as well as heterochromatin, suggesting a link between H3.3 incorporation and ZGA. Here, we examine the factors that contribute to H3.3 incorporation at ZGA. We identify a more rapid decrease in the nuclear availability of H3 than H3.3 over the final pre-ZGA cycles. We also observe an N/C ratio-dependent increase in H3.3 incorporation in mutant embryos with non-uniform local N/C ratios. We find that chaperone binding, not gene expression, controls incorporation patterns using H3/H3.3 chimeric proteins at the endogenous H3.3A locus. We test the specificity of the H3.3 chaperone pathways for H3.3 incorporation using Hira (H3.3 chaperone) mutant embryos. Overall, we propose a model in which local N/C ratios and specific chaperone binding regulate differential incorporation of H3.3 during ZGA.

PGE2 alters chromatin through H2A.Z-variant enhancer nucleosome modification to promote hematopoietic stem cell fate

Prostaglandin E2 (PGE2) and 16,16-dimethyl-PGE2 (dmPGE2) are important regulators of hematopoietic stem and progenitor cell (HSPC) fate and offer potential to enhance stem cell therapies [C. Cutler et al. Blood 122, 3074-3081(2013); W. Goessling et al. Cell Stem Cell 8, 445-458 (2011); W. Goessling et al. Cell 136, 1136-1147 (2009)]. Here, we report that PGE2-induced changes in chromatin at enhancer regions through histone-variant H2A.Z permit acute inflammatory gene induction to promote HSPC fate. We found that dmPGE2-inducible enhancers retain MNase-accessible, H2A.Z-variant nucleosomes permissive of CREB transcription factor (TF) binding. CREB binding to enhancer nucleosomes following dmPGE2 stimulation is concomitant with deposition of histone acetyltransferases p300 and Tip60 on chromatin. Subsequent H2A.Z acetylation improves chromatin accessibility at stimuli-responsive enhancers. Our findings support a model where histone-variant nucleosomes retained within inducible enhancers facilitate TF binding. Histone-variant acetylation by TF-associated nucleosome remodelers creates the accessible nucleosome landscape required for immediate enhancer activation and gene induction. Our work provides a mechanism through which inflammatory mediators, such as dmPGE2, lead to acute transcriptional changes and modify HSPC behavior to improve stem cell transplantation.

Epigenomic dynamics of early Xenopus Embryos

How the embryonic genome regulates accessibility to transcription factors is one of the major questions in understanding the spatial and temporal dynamics of gene expression during embryogenesis. Epigenomic analyses of embryonic chromatin provide molecular insights into cell-specific gene activities and genomic architectures. In recent years, significant advances have been made to elucidate the dynamic changes behind the activation of the zygotic genome in various model organisms. Here we provide an overview of the recent epigenomic studies pertaining to early Xenopus development.

Temporal Dynamic Analysis of Alternative Splicing During Embryonic Development in Zebrafish

Alternative splicing is pervasive in mammalian genomes and involved in embryo development, whereas research on crosstalk of alternative splicing and embryo development was largely restricted to mouse and human and the alternative splicing regulation during embryogenesis in zebrafish remained unclear. We constructed the alternative splicing atlas at 18 time-course stages covering maternal-to-zygotic transition, gastrulation, somitogenesis, pharyngula stages, and post-fertilization in zebrafish. The differential alternative splicing events between different developmental stages were detected. The results indicated that abundance alternative splicing and differential alternative splicing events are dynamically changed and remarkably abundant during the maternal-to-zygotic transition process. Based on gene expression profiles, we found splicing factors are expressed with specificity of developmental stage and largely expressed during the maternal-to-zygotic transition process. The better performance of cluster analysis was achieved based on the inclusion level of alternative splicing. The biological function analysis uncovered the important roles of alternative splicing during embryogenesis. The identification of isoform switches of alternative splicing provided a new insight into mining the regulated mechanism of transcript isoforms, which always is hidden by gene expression. In conclusion, we inferred that alternative splicing activation is synchronized with zygotic genome activation and discovered that alternative splicing is coupled with transcription during embryo development in zebrafish. We also unveiled that the temporal expression dynamics of splicing factors during embryo development, especially co-orthologous splicing factors. Furthermore, we proposed that the inclusion level of alternative splicing events can be employed for cluster analysis as a novel parameter. This work will provide a deeper insight into the regulation of alternative splicing during embryogenesis in zebrafish.

Epigenetic Reprogramming in Early Animal Development

Dramatic nuclear reorganization occurs during early development to convert terminally differentiated gametes to a totipotent zygote, which then gives rise to an embryo. Aberrant epigenome resetting severely impairs embryo development and even leads to lethality. How the epigenomes are inherited, reprogrammed, and reestablished in this critical developmental period has gradually been unveiled through the rapid development of technologies including ultrasensitive chromatin analysis methods. In this review, we summarize the latest findings on epigenetic reprogramming in gametogenesis and embryogenesis, and how it contributes to gamete maturation and parental-to-zygotic transition. Finally, we highlight the key questions that remain to be answered to fully understand chromatin regulation and nuclear reprogramming in early development.

Aberrant nucleosome organization in mouse SCNT embryos revealed by ULI-MNase-seq

Somatic cell nuclear transfer (SCNT) can reprogram terminally differentiated somatic cells into totipotent embryos, but with multiple defects. The nucleosome positioning, as an important epigenetic regulator for gene expression, is largely unexplored during SCNT embryonic development. Here, we mapped genome-wide nucleosome profiles in mouse SCNT embryos using ultra-low-input MNase-seq (ULI-MNase-seq). We found that the nucleosome-depleted regions (NDRs) around promoters underwent dramatic reestablishment, which is consistent with the cell cycle. Dynamics of nucleosome position in SCNT embryos were delayed compared to fertilized embryos. Subsequently, we found that the aberrant gene expression levels in inner cell mass (ICM) were positively correlated with promoter NDRs in donor cells, which indicated that the memory of nucleosome occupancy in donor cells was a potential barrier for SCNT-mediated reprogramming. We further confirmed that the histone acetylation level of donor cells was associated with the memory of promoter NDRs. Our study provides insight into nucleosome reconfiguration during SCNT preimplantation embryonic development.

Dynamic nucleosome organization after fertilization reveals regulatory factors for mouse zygotic genome activation

Chromatin remodeling is essential for epigenome reprogramming after fertilization. However, the underlying mechanisms of chromatin remodeling remain to be explored. Here, we investigated the dynamic changes in nucleosome occupancy and positioning in pronucleus-stage zygotes using ultra low-input MNase-seq. We observed distinct features of inheritance and reconstruction of nucleosome positioning in both paternal and maternal genomes. Genome-wide de novo nucleosome occupancy in the paternal genome was observed as early as 1 h after the injection of sperm into ooplasm. The nucleosome positioning pattern was continually rebuilt to form nucleosome-depleted regions (NDRs) at promoters and transcription factor (TF) binding sites with differential dynamics in paternal and maternal genomes. NDRs formed more quickly on the promoters of genes involved in zygotic genome activation (ZGA), and this formation is closely linked to histone acetylation, but not transcription elongation or DNA replication. Importantly, we found that NDR establishment on the binding motifs of specific TFs might be associated with their potential pioneer functions in ZGA. Further investigations suggested that the predicted factors MLX and RFX1 played important roles in regulating minor and major ZGA, respectively. Our data not only elucidate the nucleosome positioning dynamics in both male and female pronuclei following fertilization, but also provide an efficient method for identifying key transcription regulators during development.

The landscape of pioneer factor activity reveals the mechanisms of chromatin reprogramming and genome activation

Upon fertilization, embryos undergo chromatin reprogramming and genome activation; however, the mechanisms that regulate these processes are poorly understood. Here, we generated a triple mutant for Nanog, Pou5f3, and Sox19b (NPS) in zebrafish and found that NPS pioneer chromatin opening at >50% of active enhancers. NPS regulate acetylation across core histones at enhancers and promoters, and their function in gene activation can be bypassed by recruiting histone acetyltransferase to individual genes. NPS pioneer chromatin opening individually, redundantly, or additively depending on sequence context, and we show that high nucleosome occupancy facilitates NPS pioneering activity. Nucleosome position varies based on the input of different transcription factors (TFs), providing a flexible platform to modulate pioneering activity. Altogether, our results illuminate the sequence of events during genome activation and offer a conceptual framework to understand how pioneer factors interpret the genome and integrate different TF inputs across cell types and developmental transitions.

Establishment of developmental gene silencing by ordered polycomb complex recruitment in early zebrafish embryos

Vertebrate embryos achieve developmental competency during zygotic genome activation (ZGA) by establishing chromatin states that silence yet poise developmental genes for subsequent lineage-specific activation. Here, we reveal the order of chromatin states in establishing developmental gene poising in preZGA zebrafish embryos. Poising is established at promoters and enhancers that initially contain open/permissive chromatin with 'Placeholder' nucleosomes (bearing H2A.Z, H3K4me1, and H3K27ac), and DNA hypomethylation. Silencing is initiated by the recruitment of polycomb repressive complex 1 (PRC1), and H2Aub1 deposition by catalytic Rnf2 during preZGA and ZGA stages. During postZGA, H2Aub1 enables Aebp2-containing PRC2 recruitment and H3K27me3 deposition. Notably, preventing H2Aub1 (via Rnf2 inhibition) eliminates recruitment of Aebp2-PRC2 and H3K27me3, and elicits transcriptional upregulation of certain developmental genes during ZGA. However, upregulation is independent of H3K27me3 - establishing H2Aub1 as the critical silencing modification at ZGA. Taken together, we reveal the logic and mechanism for establishing poised/silent developmental genes in early vertebrate embryos.

RNA polymerase II clusters form in line with surface condensation on regulatory chromatin

It is essential for cells to control which genes are transcribed into RNA. In eukaryotes, two major control points are recruitment of RNA polymerase II (Pol II) into a paused state, and subsequent pause release toward transcription. Pol II recruitment and pause release occur in association with macromolecular clusters, which were proposed to be formed by a liquid-liquid phase separation mechanism. How such a phase separation mechanism relates to the interaction of Pol II with DNA during recruitment and transcription, however, remains poorly understood. Here, we use live and super-resolution microscopy in zebrafish embryos to reveal Pol II clusters with a large variety of shapes, which can be explained by a theoretical model in which regulatory chromatin regions provide surfaces for liquid-phase condensation at concentrations that are too low for canonical liquid-liquid phase separation. Model simulations and chemical perturbation experiments indicate that recruited Pol II contributes to the formation of these surface-associated condensates, whereas elongating Pol II is excluded from these condensates and thereby drives their unfolding.

Antisense RNAs during early vertebrate development are divided in groups with distinct features

Long noncoding RNAs or lncRNAs are a class of non-protein-coding RNAs that are >200 nt in length. Almost 50% of lncRNAs during zebrafish development are transcribed in an antisense direction to a protein-coding gene. However, the role of these natural antisense transcripts (NATs) during development remains enigmatic. To understand NATs in early vertebrate development, we took a computational biology approach and analyzed existing as well as novel data sets. Our analysis indicates that zebrafish NATs can be divided into two major classes based on their coexpression patterns with respect to the overlapping protein-coding genes. Group 1 NATs have characteristics similar to maternally deposited RNAs in that their levels decrease as development progresses. Group 1 NAT levels are negatively correlated with that of overlapping sense-strand protein-coding genes. Conversely, Group 2 NATs are coexpressed with overlapping protein-coding genes. In contrast to Group 1, which is enriched in genes involved in developmental pathways, Group 2 protein-coding genes are enriched in housekeeping functions. Group 1 NATs also show larger overlap and higher complementarity with the sense-strand mRNAs compared to other NATs. In addition, our transcriptomics data, quantifying RNA levels from cytoplasmic and nuclear compartments, indicates that Group 1 NATs are more abundant in the cytosol. Based on their expression pattern, cytosolic nature, and their higher complementarity to the overlapping developmental mRNAs, we speculate that Group 1 NATs function post-transcriptionally to silence spurious expression of developmental genes.

Control of zygotic genome activation in Xenopus

The fertilized frog egg contains all the materials needed to initiate development of a new organism, including stored RNAs and proteins deposited during oogenesis, thus the earliest stages of development do not require transcription. The onset of transcription from the zygotic genome marks the first genetic switch activating the gene regulatory network that programs embryonic development. Zygotic genome activation occurs after an initial phase of transcriptional quiescence that continues until the midblastula stage, a period called the midblastula transition, which was first identified in Xenopus. Activation of transcription is programmed by maternally supplied factors and is regulated at multiple levels. A similar switch exists in most animals and is of great interest both to developmental biologists and to those interested in understanding nuclear reprogramming. Here we review in detail our knowledge on this major switch in transcription in Xenopus and place recent discoveries in the context of a decades old problem.

Nucleosome Positioning and Spacing: From Mechanism to Function

Eukaryotes associate their genomes with histone proteins, forming nucleosome particles. Nucleosomes regulate and protect the genetic information. They often assemble into evenly spaced arrays of nucleosomes. These regular nucleosome arrays cover significant portions of the genome, in particular over genes. The presence of these evenly spaced nucleosome arrays is highly conserved throughout the entire eukaryotic domain. Here, we review the mechanisms behind the establishment of this primary structure of chromatin with special emphasis on the biogenesis of evenly spaced nucleosome arrays. We highlight the roles that transcription, nucleosome remodelers, DNA sequence, and histone density play towards the formation of evenly spaced nucleosome arrays and summarize our current understanding of their cellular functions. We end with key unanswered questions that remain to be explored to obtain an in-depth understanding of the biogenesis and function of the nucleosome landscape.

Nucleosome Positioning around Transcription Start Site Correlates with Gene Expression Only for Active Chromatin State in Drosophila Interphase Chromosomes

We analyzed the whole-genome experimental maps of nucleosomes in Drosophila melanogaster and classified genes by the expression level in S2 cells (RPKM value, reads per kilobase million) as well as the number of tissues in which a gene was expressed (breadth of expression, BoE). Chromatin in 5′-regions of genes we classified on four states according to the hidden Markov model (4HMM). Only the Aquamarine chromatin state we considered as Active, while the rest three states we defined as Non-Active. Surprisingly, about 20/40% of genes with 5′-regions mapped to Active/Non-Active chromatin possessed the minimal/at least modest RPKM and BoE. We found that regardless of RPKM/BoE the genes of Active chromatin possessed the regular nucleosome arrangement in 5′-regions, while genes of Non-Active chromatin did not show respective specificity. Only for genes of Active chromatin the RPKM/BoE positively correlates with the number of nucleosome sites upstream/around TSS and negatively with that downstream TSS. We propose that for genes of Active chromatin, regardless of RPKM value and BoE the nucleosome arrangement in 5′-regions potentiates transcription, while for genes of Non-Active chromatin, the transcription machinery does not require the substantial support from nucleosome arrangement to influence gene expression.

Histones, Their Variants and Post-translational Modifications in Zebrafish Development

Complex multi-cellular organisms are shaped starting from a single-celled zygote, owing to elaborate developmental programs. These programs involve several layers of regulation to orchestrate the establishment of progressively diverging cell type-specific gene expression patterns. In this scenario, epigenetic modifications of chromatin are central in influencing spatiotemporal patterns of gene transcription. In fact, it is generally recognized that epigenetic changes of chromatin states impact on the accessibility of genomic DNA to regulatory proteins. Several lines of evidence highlighted that zebrafish is an excellent vertebrate model for research purposes in the field of developmental epigenetics. In this review, I focus on the dynamic roles recently emerged for histone post-translational modifications (PTMs), histone modifying enzymes, histone variants and histone themselves in the coordination between the precise execution of transcriptional programs and developmental progression in zebrafish. In particular, I first outline a synopsis of the current state of knowledge in this field during early embryogenesis. Then, I present a survey of histone-based epigenetic mechanisms occurring throughout morphogenesis, with a stronger emphasis on cardiac formation. Undoubtedly, the issues addressed in this review take on particular importance in the emerging field of comparative biology of epigenetics, as well as in translational research.

Chromatin dynamics at the maternal to zygotic transition: recent advances from the zebrafish model

Early animal development is characterized by intense reorganization of the embryonic genome, including large-scale changes in chromatin structure and in the DNA and histone modifications that help shape this structure. Particularly profound shifts in the chromatin landscape are associated with the maternal-to-zygotic transition, when the zygotic genome is first transcribed and maternally loaded transcripts are degraded. The accessibility of the early zebrafish embryo facilitates the interrogation of chromatin during this critical window of development, making it an important model for early chromatin regulation. Here, we review our current understanding of chromatin dynamics during early zebrafish development, highlighting new advances as well as similarities and differences between early chromatin regulation in zebrafish and other species.

Inhibition of methyltransferase activity of enhancer of zeste 2 leads to enhanced lipid accumulation and altered chromatin status in zebrafish

Background

Recent studies indicate that exposure to environmental chemicals may increase susceptibility to developing metabolic diseases. This susceptibility may in part be caused by changes to the epigenetic landscape which consequently affect gene expression and lead to changes in lipid metabolism. The epigenetic modifier enhancer of zeste 2 (Ezh2) is a histone H3K27 methyltransferase implicated to play a role in lipid metabolism and adipogenesis. In this study, we used the zebrafish (Danio rerio) to investigate the role of Ezh2 on lipid metabolism and chromatin status following developmental exposure to the Ezh1/2 inhibitor PF-06726304 acetate. We used the environmental chemical tributyltin (TBT) as a positive control, as this chemical is known to act on lipid metabolism via EZH-mediated pathways in mammals.

Results

Zebrafish embryos (0–5 days post-fertilization, dpf) exposed to non-toxic concentrations of PF-06726304 acetate (5 μM) and TBT (1 nM) exhibited increased lipid accumulation. Changes in chromatin were analyzed by the assay for transposase-accessible chromatin sequencing (ATAC-seq) at 50% epiboly (5.5 hpf). We observed 349 altered chromatin regions, predominantly located at H3K27me3 loci and mostly more open chromatin in the exposed samples. Genes associated to these loci were linked to metabolic pathways. In addition, a selection of genes involved in lipid homeostasis, adipogenesis and genes specifically targeted by PF-06726304 acetate via altered chromatin accessibility were differentially expressed after TBT and PF-06726304 acetate exposure at 5 dpf, but not at 50% epiboly stage. One gene, cebpa, did not show a change in chromatin, but did show a change in gene expression at 5 dpf. Interestingly, underlying H3K27me3 marks were significantly decreased at this locus at 50% epiboly.

Conclusions

Here, we show for the first time the applicability of ATAC-seq as a tool to investigate toxicological responses in zebrafish. Our analysis indicates that Ezh2 inhibition leads to a partial primed state of chromatin linked to metabolic pathways which results in gene expression changes later in development, leading to enhanced lipid accumulation. Although ATAC-seq seems promising, our in-depth assessment of the cebpa locus indicates that we need to consider underlying epigenetic marks as well.

Brd4 and P300 Confer Transcriptional Competency during Zygotic Genome Activation

The awakening of the genome after fertilization is a cornerstone of animal development. However, the mechanisms that activate the silent genome after fertilization are poorly understood. Here, we show that transcriptional competency is regulated by Brd4- and P300-dependent histone acetylation in zebrafish. Live imaging of transcription revealed that genome activation, beginning at the miR-430 locus, is gradual and stochastic. We show that genome activation does not require slowdown of the cell cycle and is regulated through the translation of maternally inherited mRNAs. Among these, the enhancer regulators P300 and Brd4 can prematurely activate transcription and restore transcriptional competency when maternal mRNA translation is blocked, whereas inhibition of histone acetylation blocks genome activation. We conclude that P300 and Brd4 are sufficient to trigger genome-wide transcriptional competency by regulating histone acetylation on the first zygotic genes in zebrafish. This mechanism is critical for initiating zygotic development and developmental reprogramming.

Chromatin accessibility established by Pou5f3, Sox19b and Nanog primes genes for activity during zebrafish genome activation

In many organisms, early embryonic development is driven by maternally provided factors until the controlled onset of transcription during zygotic genome activation. The regulation of chromatin accessibility and its relationship to gene activity during this transition remain poorly understood. Here, we generated chromatin accessibility maps with ATAC-seq from genome activation until the onset of lineage specification. During this period, chromatin accessibility increases at regulatory elements. This increase is independent of RNA polymerase II-mediated transcription, with the exception of the hypertranscribed miR-430 locus. Instead, accessibility often precedes the transcription of associated genes. Loss of the maternal transcription factors Pou5f3, Sox19b, and Nanog, which are known to be required for zebrafish genome activation, results in decreased accessibility at regulatory elements. Importantly, the accessibility of regulatory regions, especially when established by Pou5f3, Sox19b and Nanog, is predictive for future transcription. Our results show that the maternally provided transcription factors Pou5f3, Sox19b, and Nanog open up chromatin and prime genes for activity during zygotic genome activation in zebrafish.

Combinatorial action of NF-Y and TALE at embryonic enhancers defines distinct gene expression programs during zygotic genome activation in zebrafish

Animal embryogenesis is initiated by maternal factors, but zygotic genome activation (ZGA) shifts regulatory control to the embryo during blastula stages. ZGA is thought to be mediated by maternally provided transcription factors (TFs), but few such TFs have been identified in vertebrates. Here we report that NF-Y and TALE TFs bind zebrafish genomic elements associated with developmental control genes already at ZGA. In particular, co-regulation by NF-Y and TALE is associated with broadly acting genes involved in transcriptional control, while regulation by either NF-Y or TALE defines genes in specific developmental processes, such that NF-Y controls a cilia gene expression program while TALE controls expression of hox genes. We also demonstrate that NF-Y and TALE-occupied genomic elements function as enhancers during embryogenesis. We conclude that combinatorial use of NF-Y and TALE at developmental enhancers permits the establishment of distinct gene expression programs at zebrafish ZGA.

Systemic Loss and Gain of Chromatin Architecture throughout Zebrafish Development

The spatial organization of chromosomes is critical in establishing gene expression programs. We generated in situ Hi-C maps throughout zebrafish development to gain insight into higher-order chromatin organization and dynamics. Zebrafish chromosomes segregate in active and inactive chromatin (A/B compartments), which are further organized into topologically associating domains (TADs). Zebrafish A/B compartments and TADs have genomic features similar to those of their mammalian counterparts, including evolutionary conservation and enrichment of CTCF binding sites at TAD borders. At the earliest time point, when there is no zygotic transcription, the genome is highly structured. After zygotic genome activation (ZGA), the genome loses structural features, which are re-established throughout early development. Despite the absence of structural features, we see clustering of super-enhancers in the 3D genome. Our results provide insight into vertebrate genome organization and demonstrate that the developing zebrafish embryo is a powerful model system to study the dynamics of nuclear organization.

Mechanisms regulating zygotic genome activation

Following fertilization, the two specified gametes must unite to create an entirely new organism. The genome is initially transcriptionally quiescent, allowing the zygote to be reprogrammed into a totipotent state. Gradually, the genome is activated through a process known as the maternal-to-zygotic transition, which enables zygotic gene products to replace the maternal supply that initiated development. This essential transition has been broadly characterized through decades of research in several model organisms. However, we still lack a full mechanistic understanding of how genome activation is executed and how this activation relates to the reprogramming of the zygotic chromatin architecture. Recent work highlights the central role of transcriptional activators and suggests that these factors may coordinate transcriptional activation with other developmental changes.

The maternal-to-zygotic transition revisited

The development of animal embryos is initially directed by maternal gene products. Then, during the maternal-to-zygotic transition (MZT), developmental control is handed to the zygotic genome. Extensive research in both vertebrate and invertebrate model organisms has revealed that the MZT can be subdivided into two phases, during which very different modes of gene regulation are implemented: initially, regulation is exclusively post-transcriptional and post-translational, following which gradual activation of the zygotic genome leads to predominance of transcriptional regulation. These changes in the gene expression program of embryos are precisely controlled and highly interconnected. Here, we review current understanding of the mechanisms that underlie handover of developmental control during the MZT.

Widespread Enhancer Dememorization and Promoter Priming during Parental-to-Zygotic Transition

The epigenome plays critical roles in controlling gene expression and development. However, how the parental epigenomes transit to the zygotic epigenome in early development remains elusive. Here we show that parental-to-zygotic transition in zebrafish involves extensive erasure of parental epigenetic memory, starting with methylating gametic enhancers. Surprisingly, this occurs even prior to fertilization for sperm. Both parental enhancers lose histone marks by the 4-cell stage, and zygotic enhancers are not activated until around zygotic genome activation (ZGA). By contrast, many promoters remain hypomethylated and, unexpectedly, acquire histone acetylation before ZGA at as early as the 4-cell stage. They then resolve into either activated or repressed promoters upon ZGA. Maternal depletion of histone acetyltransferases results in aberrant ZGA and early embryonic lethality. Finally, such reprogramming is largely driven by maternal factors, with zygotic products mainly contributing to embryonic enhancer activation. These data reveal widespread enhancer dememorization and promoter priming during parental-to-zygotic transition.

Packaging development: how chromatin controls transcription in zebrafish embryogenesis

How developmental gene expression is activated, co-ordinated and maintained is one of the biggest questions in developmental biology. While transcription factors lead the way in directing developmental gene expression, their accessibility to the correct repertoire of genes can depend on other factors such as DNA methylation, the presence of particular histone variants and post-translational modifications of histones. Collectively, factors that modify DNA or affect its packaging and accessibility contribute to a chromatin landscape that helps to control the timely expression of developmental genes. Zebrafish, perhaps better known for their strength as a model of embryology and organogenesis during development, are coming to the fore as a powerful model for interpreting the role played by chromatin in gene expression. Several recent advances have shown that zebrafish exhibit both similarities and differences to other models (and humans) in the way that they employ chromatin mechanisms of gene regulation. Here, I review how chromatin influences developmental transcriptional programmes during early zebrafish development, patterning and organogenesis. Lastly, I briefly highlight the importance of zebrafish chromatin research towards the understanding of human disease and transgenerational inheritance.

Pou5f3, SoxB1, and Nanog remodel chromatin on high nucleosome affinity regions at zygotic genome activation

The zebrafish embryo is transcriptionally mostly quiescent during the first 10 cell cycles, until the main wave of zygotic genome activation (ZGA) occurs, accompanied by fast chromatin remodeling. At ZGA, homologs of the mammalian stem cell transcription factors (TFs) Pou5f3, Nanog, and Sox19b bind to thousands of developmental enhancers to initiate transcription. So far, how these TFs influence chromatin dynamics at ZGA has remained unresolved. To address this question, we analyzed nucleosome positions in wild-type and maternal-zygotic (MZ) mutants for pou5f3 and nanog by MNase-seq. We show that Nanog, Sox19b, and Pou5f3 bind to the high nucleosome affinity regions (HNARs). HNARs are spanning over 600 bp, featuring high in vivo and predicted in vitro nucleosome occupancy and high predicted propeller twist DNA shape value. We suggest a two-step nucleosome destabilization-depletion model, in which the same intrinsic DNA properties of HNAR promote both high nucleosome occupancy and differential binding of TFs. In the first step, already before ZGA, Pou5f3 and Nanog destabilize nucleosomes at HNAR centers genome-wide. In the second step, post-ZGA, Nanog, Pou5f3, and SoxB1 maintain open chromatin state on the subset of HNARs, acting synergistically. Nanog binds to the HNAR center, whereas the Pou5f3 stabilizes the flanks. The HNAR model will provide a useful tool for genome regulatory studies in a variety of biological systems.

Ionizing radiation induces transgenerational effects of DNA methylation in zebrafish

Ionizing radiation is known to cause DNA damage, yet the mechanisms underlying potential transgenerational effects of exposure have been scarcely studied. Previously, we observed effects in offspring of zebrafish exposed to gamma radiation during gametogenesis. Here, we hypothesize that these effects are accompanied by changes of DNA methylation possibly inherited by subsequent generations. We assessed DNA methylation in F1 embryos (5.5 hours post fertilization) with whole genome bisulfite sequencing following parental exposure to 8.7 mGy/h for 27 days and found 5658 differentially methylated regions (DMRs). DMRs were predominantly located at known regulatory regions, such as gene promoters and enhancers. Pathway analysis indicated the involvement of DMRs related to similar pathways found with gene expression analysis, such as development, apoptosis and cancers, which could be linked to previous observed developmental defects and genomic instability in the offspring. Follow up of 19 F1 DMRs in F2 and F3 embryos revealed persistent effects up to the F3 generation at 5 regions. These results indicate that ionizing radiation related effects in offspring can be linked to DNA methylation changes that partly can persist over generations. Monitoring DNA methylation could serve as a biomarker to provide an indication of ancestral exposures to ionizing radiation.

Transcriptional defects and reprogramming barriers in somatic cell nuclear reprogramming as revealed by single-embryo RNA sequencing

Background

Nuclear reprogramming reinstates totipotency or pluripotency in somatic cells by changing their gene transcription profile. This technology is widely used in medicine, animal husbandry and other industries. However, certain deficiencies severely restrict the applications of this technology.

Results

Using single-embryo RNA-seq, our study provides complete transcriptome blueprints of embryos generated by cumulus cell (CC) donor nuclear transfer (NT), embryos generated by mouse embryonic fibroblast (MEF) donor NT and in vivo embryos at each stage (zygote, 2-cell, 4-cell, 8-cell, morula, and blastocyst). According to the results from further analyses, NT embryos exhibit RNA processing and translation initiation defects during the zygotic genome activation (ZGA) period, and protein kinase activity and protein phosphorylation are defective during blastocyst formation. Two thousand three constant genes are not able to be reprogrammed in CCs and MEFs. Among these constant genes, 136 genes are continuously mis-transcribed throughout all developmental stages. These 136 differential genes may be reprogramming barrier genes (RBGs) and more studies are needed to identify.

Conclusions

These embryonic transcriptome blueprints provide new data for further mechanistic studies of somatic nuclear reprogramming. These findings may improve the efficiency of somatic cell nuclear transfer.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5091-1) contains supplementary material, which is available to authorized users.

TALE factors use two distinct functional modes to control an essential zebrafish gene expression program

TALE factors are broadly expressed embryonically and known to function in complexes with transcription factors (TFs) like Hox proteins at gastrula/segmentation stages, but it is unclear if such generally expressed factors act by the same mechanism throughout embryogenesis. We identify a TALE-dependent gene regulatory network (GRN) required for anterior development and detect TALE occupancy associated with this GRN throughout embryogenesis. At blastula stages, we uncover a novel functional mode for TALE factors, where they occupy genomic DECA motifs with nearby NF-Y sites. We demonstrate that TALE and NF-Y form complexes and regulate chromatin state at genes of this GRN. At segmentation stages, GRN-associated TALE occupancy expands to include HEXA motifs near PBX:HOX sites. Hence, TALE factors control a key GRN, but utilize distinct DNA motifs and protein partners at different stages – a strategy that may also explain their oncogenic potential and may be employed by other broadly expressed TFs.

Link of Zygotic Genome Activation and Cell Cycle Control

The activation of the zygotic genome and onset of transcription in blastula embryos is linked to changes in cell behavior and remodeling of the cell cycle and constitutes a transition from exclusive maternal to zygotic control of development. This step in development is referred to as mid-blastula transition and has served as a paradigm for the link between developmental program and cell behavior and morphology. Here, we discuss the mechanism and functional relationships between the zygotic genome activation and cell cycle control during mid-blastula transition with a focus on Drosophila embryos.

Transcriptional Pause Sites Delineate Stable Nucleosome-Associated Premature Polyadenylation Suppressed by U1 snRNP

Regulation of RNA polymerase II (Pol II) elongation is a critical step in gene regulation. Here, we report that U1 snRNP recognition and transcription pausing at stable nucleosomes are linked through premature polyadenylation signal (PAS) termination. By generating RNA exosome conditional deletion mouse embryonic stem cells, we identified a large class of polyadenylated short transcripts in the sense direction destabilized by the RNA exosome. These PAS termination events are enriched at the first few stable nucleosomes flanking CpG islands and suppressed by U1 snRNP. Thus, promoter-proximal Pol II pausing consists of two processes: TSS-proximal and +1 stable nucleosome pausing, with PAS termination coinciding with the latter. While pausing factors NELF/DSIF only function in the former step, flavopiridol-sensitive mechanism(s) and Myc modulate both steps. We propose that premature PAS termination near the nucleosome-associated pause site represents a common transcriptional elongation checkpoint regulated by U1 snRNP recognition, nucleosome stability, and Myc activity.

Cohesin facilitates zygotic genome activation in zebrafish

At zygotic genome activation (ZGA), changes in chromatin structure are associated with new transcription immediately following the maternal-to-zygotic transition (MZT). The nuclear architectural proteins cohesin and CCCTC-binding factor (CTCF) contribute to chromatin structure and gene regulation. We show here that normal cohesin function is important for ZGA in zebrafish. Depletion of the cohesin subunit Rad21 delays ZGA without affecting cell cycle progression. In contrast, CTCF depletion has little effect on ZGA, whereas complete abrogation is lethal. Genome-wide analysis of Rad21 binding reveals a change in distribution from pericentromeric satellite DNA and other locations, including the miR-430 locus (the products of which are responsible for maternal transcript degradation), to genes, as embryos progress through the MZT. After MZT, a subset of Rad21 binding overlaps the pioneer factor Pou5f3, which activates early expressed genes. Rad21 depletion disrupts the formation of nucleoli and RNA polymerase II foci, suggestive of global defects in chromosome architecture. We propose that Rad21/cohesin redistribution to active areas of the genome is key to the establishment of chromosome organization and the embryonic developmental program.

Hypomethylated domain-enriched DNA motifs prepattern the accessible nucleosome organization in teleosts

Background

Gene promoters in vertebrate genomes show distinct chromatin features such as stably positioned nucleosome array and DNA hypomethylation. The nucleosomes are known to have certain sequence preferences, and the prediction of nucleosome positioning from DNA sequence has been successful in some organisms such as yeast. However, at gene promoters where nucleosomes are much more stably positioned than in other regions, the sequence-based model has failed to work well, and sequence-independent mechanisms have been proposed.

Results

Using DNase I-seq in medaka embryos, we demonstrated that hypomethylated domains (HMDs) specifically possess accessible nucleosome organization with longer linkers, and we reassessed the DNA sequence preference for nucleosome positioning in these specific regions. Remarkably, we found with a supervised machine learning algorithm, k-mer SVM, that nucleosome positioning in HMDs is accurately predictable from DNA sequence alone. Specific short sequences (6-mers) that contribute to the prediction are specifically enriched in HMDs and distribute periodically with approximately 200-bp intervals which prepattern the position of accessible linkers. Surprisingly, the sequence preference of the nucleosome and linker in HMDs is opposite from that reported previously. Furthermore, the periodicity of specific motifs at hypomethylated promoters was conserved in zebrafish.

Conclusion

This study reveals strong link between nucleosome positioning and DNA sequence at vertebrate promoters, and we propose hypomethylated DNA-specific regulation of nucleosome positioning.

Electronic supplementary material

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

Uniform gene expression in embryos is achieved by temporal averaging of transcription noise

Here, Stapel et al. investigated how stochastic transcription activation can result in uniform expression during development of cellularized embryos. Using smFISH and mathematical modeling, the authors show that during zebrafish embryogenesis, transcription activation is stochastic due to (1) genes acquiring transcriptional competence at different times in different cells, (2) differences in cell cycle stage between cells, and (3) the stochastic nature of transcription.

Transcription is often stochastic. This is seemingly incompatible with the importance of gene expression during development. Here we show that during zebrafish embryogenesis, transcription activation is stochastic due to (1) genes acquiring transcriptional competence at different times in different cells, (2) differences in cell cycle stage between cells, and (3) the stochastic nature of transcription. Initially, stochastic transcription causes large cell-to-cell differences in transcript levels. However, variability is reduced by lengthening cell cycles and the accumulation of transcription events in each cell. Temporal averaging might provide a general context in which to understand how embryos deal with stochastic transcription.

Zygotic Genome Activation in Vertebrates

The first major developmental transition in vertebrate embryos is the maternal-to-zygotic transition (MZT) when maternal mRNAs are degraded and zygotic transcription begins. During the MZT, the embryo takes charge of gene expression to control cell differentiation and further development. This spectacular organismal transition requires nuclear reprogramming and the initiation of RNAPII at thousands of promoters. Zygotic genome activation (ZGA) is mechanistically coordinated with other embryonic events, including changes in the cell cycle, chromatin state, and nuclear-to-cytoplasmic component ratios. Here, we review progress in understanding vertebrate ZGA dynamics in frogs, fish, mice, and humans to explore differences and emphasize common features.

Using local chromatin structure to improve CRISPR/Cas9 efficiency in zebrafish

Although the CRISPR/Cas9 has been successfully applied in zebrafish, considerable variations in efficiency have been observed for different gRNAs. The workload and cost of zebrafish mutant screening is largely dependent on the mutation rate of injected embryos; therefore, selecting more effective gRNAs is especially important for zebrafish mutant construction. Besides the sequence features, local chromatin structures may have effects on CRISPR/Cas9 efficiency, which remain largely unexplored. In the only related study in zebrafish, nucleosome organization was not found to have an effect on CRISPR/Cas9 efficiency, which is inconsistent with recent studies in vitro and in mammalian cell lines. To understand the effects of local chromatin structure on CRISPR/Cas9 efficiency in zebrafish, we first determined that CRISPR/Cas9 introduced genome editing mainly before the dome stage. Based on this observation, we reanalyzed our published nucleosome organization profiles and generated chromatin accessibility profiles in the 256-cell and dome stages using ATAC-seq technology. Our study demonstrated that chromatin accessibility showed positive correlation with CRISPR/Cas9 efficiency, but we did not observe a clear correlation between nucleosome organization and CRISPR/Cas9 efficiency. We constructed an online database for zebrafish gRNA selection based on local chromatin structure features that could prove beneficial to zebrafish homozygous mutant construction via CRISPR/Cas9.

CAM: A quality control pipeline for MNase-seq data

Nucleosome organization affects the accessibility of cis-elements to trans-acting factors. Micrococcal nuclease digestion followed by high-throughput sequencing (MNase-seq) is the most popular technology used to profile nucleosome organization on a genome-wide scale. Evaluating the data quality of MNase-seq data remains challenging, especially in mammalian. There is a strong need for a convenient and comprehensive approach to obtain dedicated quality control (QC) for MNase-seq data analysis. Here we developed CAM, which is a comprehensive QC pipeline for MNase-seq data. The CAM pipeline provides multiple informative QC measurements and nucleosome organization profiles on different potentially functional regions for given MNase-seq data. CAM also includes 268 historical MNase-seq datasets from human and mouse as a reference atlas for unbiased assessment. CAM is freely available at: http://www.tongji.edu.cn/~zhanglab/CAM.

Chromatin Control of Developmental Dynamics and Plasticity

Chromatin structure is intimately connected with gene expression and cell identity. Here we review recent advances in the field and discuss how establishment of cell identity during development is accompanied by large-scale remodeling of the epigenetic landscape and how this remodeling drives and supports lineage specification and maintenance. We discuss maternal control of the early embryonic epigenetic landscape, selective usage of enhancer clusters via 3D chromatin contacts leading to activation of transcription factor networks, and conserved regulation of developmental pathways by specific DNA demethylation of key regulatory regions. Together, these processes establish an epigenetic framework regulating different phases of embryonic development.

A comparison of nucleosome organization in Drosophila cell lines

Changes in the distribution of nucleosomes along the genome influence chromatin structure and impact gene expression by modulating the accessibility of DNA to transcriptional machinery. However, the role of genome-wide nucleosome positioning in gene expression and in maintaining differentiated cell states remains poorly understood. Drosophila melanogaster cell lines represent distinct tissue types and exhibit cell-type specific gene expression profiles. They thus could provide a useful tool for investigating cell-type specific nucleosome organization of an organism's genome. To evaluate this possibility, we compared genome-wide nucleosome positioning and occupancy in five different Drosophila tissue-specific cell lines, and in reconstituted chromatin, and then tested for correlations between nucleosome positioning, transcription factor binding motifs, and gene expression. Nucleosomes in all cell lines were positioned in accordance with previously known DNA-nucleosome interactions, with helically repeating A/T di-nucleotide pairs arranged within nucleosomal DNAs and AT-rich pentamers generally excluded from nucleosomal DNA. Nucleosome organization in all cell lines differed markedly from in vitro reconstituted chromatin, with highly expressed genes showing strong nucleosome organization around transcriptional start sites. Importantly, comparative analysis identified genomic regions that exhibited cell line-specific nucleosome enrichment or depletion. Further analysis of these regions identified 91 out of 16,384 possible heptamer sequences that showed differential nucleosomal occupation between cell lines, and 49 of the heptamers matched one or more known transcription factor binding sites. These results demonstrate that there is differential nucleosome positioning between these Drosophila cell lines and therefore identify a system that could be used to investigate the functional significance of differential nucleosomal positioning in cell type specification.

Competition between histone and transcription factor binding regulates the onset of transcription in zebrafish embryos

Upon fertilization, the genome of animal embryos remains transcriptionally inactive until the maternal-to-zygotic transition. At this time, the embryo takes control of its development and transcription begins. How the onset of zygotic transcription is regulated remains unclear. Here, we show that a dynamic competition for DNA binding between nucleosome-forming histones and transcription factors regulates zebrafish genome activation. Taking a quantitative approach, we found that the concentration of non-DNA-bound core histones sets the time for the onset of transcription. The reduction in nuclear histone concentration that coincides with genome activation does not affect nucleosome density on DNA, but allows transcription factors to compete successfully for DNA binding. In agreement with this, transcription factor binding is sensitive to histone levels and the concentration of transcription factors also affects the time of transcription. Our results demonstrate that the relative levels of histones and transcription factors regulate the onset of transcription in the embryo.

DOI:http://dx.doi.org/10.7554/eLife.23326.001

The DNA in a fertilized egg contains all the information required to form an animal’s body. In order for the animal to develop properly, particular genes encoded in the DNA are only active at specific times. The DNA is wrapped around proteins called histones, which allows the DNA to be tightly packed inside the cell. However, histones can block other proteins called transcription factors from binding to the DNA to activate the genes. Young embryos initially develop with all of their genes switched off, relying on the nutrients and other molecules provided by their mother. After some time, the embryo starts to switch on its own genes to take control of its own development, but it was not clear how this happens.

Joseph et al. investigated how genes are activated in zebrafish embryos, which are often used as models to study how animals develop. The experiments show that competition between histones and transcription factors for binding to DNA controls when genes are switched on. In young fish embryos, there are so many histones present that transcription factors have no opportunity to bind to DNA. Over time, however, the numbers of histones decrease, allowing transcription factors to bind to DNA and switch on genes.

Histones and transcription factors regulate the activity of genes throughout the life of the animal. Therefore, competition between these two types of protein may also control gene activity in other situations. A better understanding of how gene activity is controlled could allow researchers to more easily grow different types of cell in the laboratory or to reprogram specific cells in the body. As such, these new findings may aid the development of therapies to regenerate organs or tissues that have been damaged by injury or disease.

DOI:http://dx.doi.org/10.7554/eLife.23326.002

Transcriptome dynamics in early zebrafish embryogenesis determined by high-resolution time course analysis of 180 successive, individual zebrafish embryos

Background

Recently, much progress has been made in the field of gene-expression in early embryogenesis. However, the dynamic behaviour of transcriptomes in individual embryos has hardly been studied yet and the time points at which pools of embryos are collected are usually still quite far apart. Here, we present a high-resolution gene-expression time series with 180 individual zebrafish embryos, obtained from nine different spawns, developmentally ordered and profiled from late blastula to mid-gastrula stage. On average one embryo per minute was analysed. The focus was on identification and description of the transcriptome dynamics of the expressed genes in this embryonic stage, rather than to biologically interpret profiles in cellular processes and pathways.

Results

In the late blastula to mid-gastrula stage, we found 6,734 genes being expressed with low variability and rather gradual changes. Ten types of dynamic behaviour were defined, such as genes with continuously increasing or decreasing expression, and all expressed genes were grouped into these types. Also, the exact expression starting and stopping points of several hundred genes during this developmental period could be pinpointed. Although the resolution of the experiment was so high, that we were able to clearly identify four known oscillating genes, no genes were observed with a peaking expression. Additionally, several genes showed expression at two or three distinct levels that strongly related to the spawn an embryo originated from.

Conclusion

Our unique experimental set-up of whole-transcriptome analysis of 180 individual embryos, provided an unparalleled in-depth insight into the dynamics of early zebrafish embryogenesis. The existence of a tightly regulated embryonic transcriptome program, even between individuals from different spawns is shown. We have made the expression profile of all genes available for domain experts. The fact that we were able to separate the different spawns by their gene-expression variance over all expressed genes, underlines the importance of spawn specificity, as well as the unexpectedly tight gene-expression regulation in early zebrafish embryogenesis.

Electronic supplementary material

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

Establishing nucleosome architecture and stability at promoters: Roles of pioneer transcription factors and the RSC chromatin remodeler

Improvements in deep sequencing, together with methods to rapidly deplete essential transcription factors (TFs) and chromatin remodelers, have recently led to a more detailed picture of promoter nucleosome architecture in yeast and its relationship to transcriptional regulation. These studies revealed that ∼40% of all budding yeast protein-coding genes possess a unique promoter structure, where we propose that an unusually unstable nucleosome forms immediately upstream of the transcription start site (TSS). This "fragile" nucleosome (FN) promoter architecture relies on the combined action of the essential RSC (Remodels Structure of Chromatin) nucleosome remodeler and pioneer transcription factors (PTFs). FNs are associated with genes whose expression is high, coupled to cell growth, and characterized by low cell-to-cell variability (noise), suggesting that they may promote these features. Recent studies in metazoans suggest that the presence of dynamic nucleosomes upstream of the TSS at highly expressed genes may be conserved throughout evolution.

NucTools: analysis of chromatin feature occupancy profiles from high-throughput sequencing data

Background

Biomedical applications of high-throughput sequencing methods generate a vast amount of data in which numerous chromatin features are mapped along the genome. The results are frequently analysed by creating binary data sets that link the presence/absence of a given feature to specific genomic loci. However, the nucleosome occupancy or chromatin accessibility landscape is essentially continuous. It is currently a challenge in the field to cope with continuous distributions of deep sequencing chromatin readouts and to integrate the different types of discrete chromatin features to reveal linkages between them.

Results

Here we introduce the NucTools suite of Perl scripts as well as MATLAB- and R-based visualization programs for a nucleosome-centred downstream analysis of deep sequencing data. NucTools accounts for the continuous distribution of nucleosome occupancy. It allows calculations of nucleosome occupancy profiles averaged over several replicates, comparisons of nucleosome occupancy landscapes between different experimental conditions, and the estimation of the changes of integral chromatin properties such as the nucleosome repeat length. Furthermore, NucTools facilitates the annotation of nucleosome occupancy with other chromatin features like binding of transcription factors or architectural proteins, and epigenetic marks like histone modifications or DNA methylation. The applications of NucTools are demonstrated for the comparison of several datasets for nucleosome occupancy in mouse embryonic stem cells (ESCs) and mouse embryonic fibroblasts (MEFs).

Conclusions

The typical workflows of data processing and integrative analysis with NucTools reveal information on the interplay of nucleosome positioning with other features such as for example binding of a transcription factor CTCF, regions with stable and unstable nucleosomes, and domains of large organized chromatin K9me2 modifications (LOCKs). As potential limitations and problems we discuss how inter-replicate variability of MNase-seq experiments can be addressed.

Electronic supplementary material

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

Genomes of Multicellular Organisms Have Evolved to Attract Nucleosomes to Promoter Regions

Sequences that influence nucleosome positioning in promoter regions, and their relation to gene regulation, have been the topic of much research over the last decade. In yeast, significant nucleosome-depleted regions are found, which facilitate transcription. With the arrival of nucleosome positioning maps for the human genome, it was discovered that in our genome, unlike in that of yeast, promoters encode for high nucleosome occupancy. In this work, we look at the genomes of a range of different organisms, to provide a catalog of nucleosome positioning signals in promoters across the tree of life. We utilize a computational model of the nucleosome, based on crystallographic analyses of the structure and elasticity of the nucleosome, to predict the nucleosome positioning signals in promoter regions. To be able to apply our model to large genomic datasets, we introduce an approximative scheme that makes use of the limited range of correlations in nucleosomal sequence preferences to create a computationally efficient approximation of the full biophysical model. Our predictions show that a clear distinction between unicellular and multicellular life is visible in the intrinsically encoded nucleosome affinity. Furthermore, the strength of the nucleosome positioning signals correlates with the complexity of the organism. We conclude that encoding for high nucleosome occupancy, as in the human genome, is in fact a universal feature of multicellular life.

The timing of zygotic genome activation

After fertilization, the embryonic genome is inactive until transcription is initiated during the maternal-to-zygotic transition. How the onset of transcription is regulated in a precisely timed manner, however, is a long standing question in biology. Several mechanisms have been shown to contribute to the temporal regulation of genome activation but none of them can fully explain the general absence of transcription as well the gene specific onset that follows. Here we review the work that has been done toward elucidating the mechanisms underlying the temporal regulation of transcription in embryos.