Modelling genome-wide topological associating domains in mouse embryonic stem cells

Chromosome conformation capture (3C)-based techniques such as chromosome conformation capture carbon copy (5C) and Hi-C revealed that the folding of mammalian chromosomes is highly hierarchical. A fundamental structural unit in the hierarchy is represented by topologically associating domains (TADs), sub-megabase regions of the genome within which the chromatin fibre preferentially interacts. 3C-based methods provide the mean contact probabilities between chromosomal loci, averaged over a large number of cells, and do not give immediate access to the single-cell conformations of the chromatin fibre. However, coarse-grained polymer models based on 5C data can be used to extract the single-cell conformations of single TADs. Here, we extend this approach to analyse around 2500 TADs in murine embryonic stem cells based on high-resolution Hi-C data. This allowed to predict the cell-to-cell variability in single contacts within genome-wide TADs and correlations between them. Based on these results, we predict that TADs are more similar to ideal chains than to globules in terms of their physical size and three-dimensional shape distribution. Furthermore, we show that their physical size and the degree of structural anisotropy of single TADs are correlated with the level of transcriptional activity of the genes that it harbours. Finally, we show that a large number of multiplets of genomic loci co-localize more often than expected by random, and these loci are particularly enriched in promoters, enhancers and CTCF-bound sites. These results provide the first genome-wide structural reconstruction of TADs using polymeric models obeying the laws of thermodynamics and reveal important universal trends in the correlation between chromosome structure and transcription.