Abstract:The zebrafish embryo remains transcriptionally quiescent during the first 10 cell cycles. Only then Zygotic Genome Activation (ZGA) occurs and is accompanied by fast chromatin remodeling. At ZGA, homologs of mammalian stem cell transcription factors (TFs) Pou5f3/Oct4, Nanog and Sox19b bind to thousands of developmental enhancers to initiate 5 transcription. So far, how these three ZGA TFs influence chromatin dynamics at ZGA has remained unresolved. To address this question, we analyzed before and after ZGA nucleosome positions in wild-type and Maternal-Zygotic (MZ) mutants for pou5f3 and nanog. 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 10 and high propeller twist DNA shape value. Just prior to ZGA, Pou5f3 and Nanog nonspecifically compete with histones and reduce nucleosome occupancy on strong nucleosome positioning sequences genome-wide. After ZGA, specific binding of Nanog and Pou5f3/SoxB1 complex is necessary to maintain open chromatin state on more than 6,000 HNARs. Nanog binds to the strongest nucleosome positioning sequence within HNAR, while the Pou5f3/SoxB1 15 complex binds to the flanks of it. We suggest a model, where the same intrinsic DNA properties of HNAR promote both high nucleosome occupancy and differential binding of TFs. We hope that our model will provide a useful framework for genome regulatory studies in the variety of biological systems.