Decoding Nucleosome-Depleted Regions: Insights from Epigenetic Marks, Nucleosome Size, and Thermodynamic Modelling

Elucidating the global and local rules that govern genome-wide nucleosome organisation and chromatin architecture remains a critical challenge. Thermodynamic modelling based on DNA elastic properties predicts the presence of sequence-encoded nucleosome-inhibiting energy barriers (NIEBs) along vertebrate genomes. They delineate in vivo nucleosome-depleted regions (NDRs) flanked by 2-3 well positioned nucleosomes. Here, we compared mouse NIEBs to NDRs observed at CTCF binding sites and active TSSs to reveal specific chromatin organizations. We uncover in MNase-seq chromatin profiles the presence of particles of subnucleosomal length specifically positioned at the border of NIEBs with an enrichment of H3.3 and its modification H3.3 S31Ph, whereas the positioning of nucleosomes bearing H3K27ac appears insensitive to NIEBs. Surprisingly, post-translational modifications affect the size distribution of nucleosomes as seen by MNase digestion and so likely their breathing capability. We implemented an extension of our thermodynamic model allowing for variable particle size and suggest that subnucleomes at NIEB borders would result from the recruitment of chromatin remodellers at NIEBs. Our findings provide new insights into the mechanisms by which the DNA sequence and epigenetic marks shape the nucleosome positioning and breathing.