Abstract
Epithelial-to-mesenchymal transition (EMT) is a complex multistep process in which phenotype switches are mediated by a network of transcription factors (TFs). Systematic characterization of all dynamic TFs controlling EMT state transitions, especially for the intermediate partial-EMT state, represents a highly relevant yet largely unexplored task. Here, we performed a computational analysis that integrated time-course EMT transcriptomic data with public cistromic data and identified three synergistic master TFs (ETS2, HNF4A and JUNB) that regulate the transition through the partial-EMT state. Overexpression of these regulators predicted a poor clinical outcome, and their elimination readily abolished TGF-β-induced EMT. Importantly, these factors utilized a clique motif, physically interact and their cumulative binding generally characterized EMT-associated genes. Furthermore, analyses of H3K27ac ChIP-seq data revealed that ETS2, HNF4A and JUNB are associated with super-enhancers and the administration of BRD4 inhibitor readily abolished TGF-β-induced EMT. These findings have implications for systematic discovery of master EMT regulators and super-enhancers as novel targets for controlling metastasis.