Abstract

In interphase nuclei, chromatin is organized into interspersed dense domains with characteristic sizes, both in the nuclear interior and periphery. However, the quantitative impact of transcription and histone modifications on the size and distribution of these domains remains unclear. Here, we introduce a mesoscale theoretical model that investigates the relationship between heterochromatic domain sizes and loop extrusion rates from these domains. The model considers chromatin-chromatin and chromatin-lamina interactions, methylation and acetylation kinetics, and diffusion of epigenetic marks and nucleoplasm. Our model generates testable predictions that help reveal the biophysics underlying chromatin organization in the presence of transcription-driven loop extrusion. This process is kinetically captured through the conversion of heterochromatin to euchromatin in response to RNAPII activity. We discovered that a balance between diffusive and reactive fluxes governs the steady-state sizes of heterochromatin domains. Using theory and simulations, we predicted that a loss of transcription results in increased chromatin compaction and larger heterochromatin domain sizes. To validate our predictions, we employed complementary super-resolution and nano-imaging techniques on five different cell lines with impaired transcription. We quantitatively assessed how domain sizes scale with loop extrusion rates at the hetero-euchromatin interfaces. Our analysis of previously obtained super-resolution images of nuclei revealed that excessive loop extrusion leads to smaller heterochromatin domains. The model successfully recapitulated these observations, explaining how transcription loss can counteract the effects of cohesin overloading. As the general biophysical mechanisms regulating heterochromatin domain sizes are independent of cell type, our findings have significant implications for understanding the role of transcription in global genome organization.