Bimodal specificity of TF–DNA recognition in embryonic stem cells

Transcription factors (TFs) bind genomic DNA regulating gene expression and developmental programs in embryonic stem cells (ESCs). Even though comprehensive genome-wide molecular maps for TF-DNA binding are experimentally available for key pluripotency-associated TFs, the understanding of molecular design principles responsible for TF-DNA recognition remains incomplete. Here, we show that binding preferences of key pluripotency TFs, such as Pou5f1 (Oct4), Smad1, Otx2, Srf, and Nanog, exhibit bimodality in the local GC-content distribution. Sequence-dependent binding specificity of these TFs is distributed across three major contributions. First, local GC-content is dominant in high-GC-content regions. Second, recognition of specific k-mers is predominant in low-GC-content regions. Third, short tandem repeats (STRs) are highly predictive in both low- and high-GC-content regions. In sharp contrast, the binding preferences of c-Myc are exclusively dominated by local GC-content and STRs in high-GC-content genomic regions. We demonstrate that the transition in the TF-DNA binding landscape upon ESC differentiation is regulated by the concentration of c-Myc, which forms a bivalent c-Myc-Max heterotetramer upon promoter binding, competing with key pluripotency factors such as Smad1. Finally, a direct interaction between c-Myc and key pluripotency factors is not required to achieve this transition.