Abstract
Epigenetic modifications can extend over long genomic regions to form domain-level chromatin states that play critical roles in gene regulation. The molecular mechanism for the establishment and maintenance of these states is not fully understood and remains challenging to study with existing experimental techniques. Here, we took a data-driven approach and parameterized an information-theoretic model to infer the formation mechanism of domain-level chromatin states from genome-wide epigenetic modification profiles. This model reproduces statistical correlations among histone modifications and identifies well-known states. Importantly, it predicts drastically different mechanisms and kinetic pathways for the formation of euchromatin and heterochromatin. In particular, long, strong enhancer and promoter states grow gradually from short but stable regulatory elements via a multistep process. On the other hand, the formation of heterochromatin states is highly cooperative, and no intermediate states are found along the transition path. This cooperativity can arise from a chromatin looping-mediated spreading of histone methylation mark and supports collapsed, globular three-dimensional conformations rather than regular fibril structures for heterochromatin. We further validated these predictions using changes of epigenetic profiles along cell differentiation. Our study demonstrates that information-theoretic models can go beyond statistical analysis to derive insightful kinetic information that is otherwise difficult to access.