Cohesin promotes stochastic domain intermingling to ensure proper regulation of boundary-proximal genes

The mammalian genome can be segmented into thousands of topologically associated domains (TADs) based on chromosome conformation capture studies, such as Hi-C. TADs have been proposed to act as insulated neighborhoods, spatially sequestering and insulating the enclosed genes and regulatory elements through chromatin looping and self-association. Recent results indicate that inter-TAD interactions can also occur, suggesting boundaries may be semi-permissible. However, the nature, extent, and function, if any, of these inter-TAD interactions remains unclear. Here, we combine super-and high-resolution microscopy with Oligopaint technology to precisely quantify the interaction frequency within and between neighboring domains in human cells. We find that intermingling across domain boundaries is a widespread feature of the human genome, with varying levels of interactions across different loci that correlate with their differing boundary strengths by Hi-C. Moreover, we find that cohesin depletion, which is known to abolish TADs at the population-average level, does not induce ectopic interactions but instead reduces both intra- and inter-domain interactions to a similar extent. Reduced chromatin intermixing due to cohesin loss affects domain incorporation and transcriptional bursting frequencies of genes close to architectural boundaries, potentially explaining the gene expression changes observed in the cohesinopathy Cornelia de Lange syndrome. Together, our results provide a mechanistic explanation for stochastic domain intermingling, arguing that cohesin partially bypasses boundaries to promote alternating incorporation of boundary-proximal genes into neighboring regulatory domains.