Abstract
Chromatin states are highly cell-type-specific, but the underlying mechanisms for the establishment and maintenance of their genome-wide patterns remain poorly understood. Here we present a computational approach for investigation of chromatin-state plasticity. We applied this approach to investigate an ENCODE ChIP-seq dataset profiling the genome-wide distributions of the H3K27me3 mark in 19 human cell lines. We found that the high plasticity regions (HPRs) can be divided into two functionally and mechanistically distinct subsets, which correspond to CpG island (CGI) proximal or distal regions, respectively. Although the CGI proximal HPRs are typically associated with continuous variation across different cell-types, the distal HPRs are associated with binary-like variations. We developed a computational approach to predict putative cell-type-specific modulators of H3K27me3 patterns and validated the predictions by comparing with public ChIP-seq data. Furthermore, we applied this approach to investigate mechanisms for poised enhancer establishment in primary human erythroid precursors. Importantly, we predicted and experimentally validated that the principal hematopoietic regulator T-cell acute lymphocytic leukemia-1 (TAL1) is involved in regulating H3K27me3 variations in collaboration with the transcription factor growth factor independent 1B (GFI1B), providing fresh insights into the context-specific role of TAL1 in erythropoiesis. Our approach is generally applicable to investigate the regulatory mechanisms of epigenetic pathways in establishing cellular identity.