To gain insights into coordinated lineage-specification and morphogenetic processes during early embryogenesis, here we report a systematic identification of transcriptional programs mediated by a key developmental regulator-Brachyury. High-resolution chromosomal localization mapping of Brachyury by ChIP sequencing and ChIP-exonuclease revealed distinct sequence signatures enriched in Brachyury-bound enhancers. A combination of genome-wide in vitro and in vivo perturbation analysis and cross-species evolutionary comparison unveiled a detailed Brachyury-dependent gene-regulatory network that directly links the function of Brachyury to diverse developmental pathways and cellular housekeeping programs. We also show that Brachyury functions primarily as a transcriptional activator genome-wide and that an unexpected gene-regulatory feedback loop consisting of Brachyury, Foxa2, and Sox17 directs proper stem-cell lineage commitment during streak formation. Target gene and mRNA-sequencing correlation analysis of the T(c) mouse model supports a crucial role of Brachyury in up-regulating multiple key hematopoietic and muscle-fate regulators. Our results thus chart a comprehensive map of the Brachyury-mediated gene-regulatory network and how it influences in vivo developmental homeostasis and coordination.

Enhancer-binding pluripotency regulators (Sox2 and Oct4) play a seminal role in embryonic stem (ES) cell-specific gene regulation. Here, we combine in vivo and in vitro single-molecule imaging, transcription factor (TF) mutagenesis, and ChIP-exo mapping to determine how TFs dynamically search for and assemble on their cognate DNA target sites. We find that enhanceosome assembly is hierarchically ordered with kinetically favored Sox2 engaging the target DNA first, followed by assisted binding of Oct4. Sox2/Oct4 follow a trial-and-error sampling mechanism involving 84-97 events of 3D diffusion (3.3-3.7 s) interspersed with brief nonspecific collisions (0.75-0.9 s) before acquiring and dwelling at specific target DNA (12.0-14.6 s). Sox2 employs a 3D diffusion-dominated search mode facilitated by 1D sliding along open DNA to efficiently locate targets. Our findings also reveal fundamental aspects of gene and developmental regulation by fine-tuning TF dynamics and influence of the epigenome on target search parameters.