Dynamic combinatorial interactions of RUNX1 and cooperating partners regulates megakaryocytic differentiation in cell line models

Specific interactions of transcription factors (TFs) with their targets are crucial for specifying gene expression programs during cell differentiation. How specificity is maintained despite limited selectivity of individual TF-DNA interactions is not fully understood. RUNX1 TF is among the most frequently mutated genes in human leukemia and an important regulator of megakaryopoiesis. We used megakaryocytic cell lines to characterize the network of RUNX1 targets and cooperating TFs in differentiating megakaryocytes and demonstrated how dynamic partnerships between RUNX1 and cooperating TFs facilitated regulatory plasticity and specificity during this process. After differentiation onset, RUNX1 directly activated a large number of genes through interaction with preexisting and de novo binding sites. Recruitment of RUNX1 to de novo occupied sites occurred at H3K4me1-marked preprogrammed enhancers. A significant number of these de novo bound sites lacked RUNX motif but were occupied by AP-1 TFs. Reciprocally, AP-1 TFs were up-regulated by RUNX1 after 12-O-tetradecanoylphorbol-13-acetate induction and recruited to RUNX1-occupied sites lacking AP-1 motifs. At other differentiation stages, additional combinatorial interactions occurred between RUNX1 and its coregulators, GATA1 and ETS. The findings suggest that in differentiating megakaryocytic cell lines, RUNX1 cooperates with GATA1, AP-1, and ETS to orchestrate cell-specific transcription programs through dynamic TF partnerships.