The genome is organised via CTCF/Cohesin binding sites, which partition chromosomes into 1-5Mb topologically associated domains (TADs), and further into smaller sub-domains (sub-TADs). Here we examined in vivo an ~80kb sub-TAD, containing the mouse α-globin gene cluster, lying within a ~1Mb TAD. We find that the sub-TAD is flanked by predominantly convergent CTCF/cohesin sites which are ubiquitously bound by CTCF but only interact during erythropoiesis, defining a self-interacting erythroid compartment. Whereas the α-globin regulatory elements normally act solely on promoters downstream of the enhancers, removal of a conserved upstream CTCF/cohesin boundary extends the sub-TAD to adjacent upstream CTCF/cohesin binding sites. The α-globin enhancers now interact with the flanking chromatin, upregulating expression of genes within this extended sub-TAD. Rather than acting solely as a barrier to chromatin modification, CTCF/cohesin boundaries in this sub-TAD delimit the region of chromatin to which enhancers have access and within which they interact with receptive promoters.

Complex gene regulation is one of the key requirements for the evolution of higher eukaryotes. 1 In these organisms, many genes are regulated by enhancers that are 10 4 -10 6 base pairs (bp) distant from the promoter. Enhancer sequences usually contain multiple small transcription factor binding sites (typically ~10bp), and physical contact between the promoter and enhancer is thought to be required to modulate gene expression. 2 Current methods have extensively defined chromatin architecture at scales above 1 kb but until now it has not been possible to define physical contacts at the scale of the key proteins determining gene expression. Here we define the interactions between different classes of regulatory elements (enhancers, promoters and boundary elements) in unprecedented detail, using a novel chromosome conformation capture method (Micro Capture-C (MCC)), which allows physical contacts to be determined at base-pair resolution. We find that highly punctate contacts occur between enhancers, promoters and CCCTC-binding factor (CTCF) sites and we show, using base pair resolution plots of ligation junctions, that transcription factors generate a key component of the contacts between enhancers and promoters. Our data show that contacts from CTCF sites highly correlate with cooccupancy of cohesin and that interactions between CTCF sites are increased when active promoters and enhancers are located within the intervening chromatin. We also find that promoters make the strongest contacts with both enhancers and CTCF sites and that while CTCF sites contact promoters strongly they only make weak contacts with enhancers. The highly punctate nature of the contacts is an unexpected finding because the current view is that physical contacts are constrained by much larger domains such as topological associated domains (TADs). 3 Our results support a model in which chromatin loop extrusion 4-6 is dependent on cohesin loading at active promoters and enhancers, explaining the formation of tissue-specific chromatin domains without changes in CTCF binding. The data suggest that a separate mechanism to loop extrusion underlies enhancer-/promoter contacts, which likely involves DNA binding proteins at enhancers and promoters. The unprecedented

Predicting the impact of non-coding genetic variation requires interpreting it in the context of 3D genome architecture. We have developed deepC, a transfer learning based deep neural network that accurately predicts genome folding from megabase-scale DNA sequence. DeepC predicts domain boundaries at high-resolution, learns the sequence determinants of genome folding and predicts the impact of both large-scale structural and single base pair variations.