High-resolution Chromosome Conformation Research Articles

A chromosome-level genome assembly of a deep-sea symbiotic Aplacophora mollusc Chaetoderma sp.

The worm-shaped, shell-less Caudofoveata is one of the least known groups of molluscs. As early-branching molluscs, the lack of high-quality genomes hinders our understanding of their evolution and ecology. Here, we report a high-quality chromosome-scale genome of Chaetoderma sp. combining PacBio, Illumina, and high-resolution chromosome conformation capture sequencing. The final assembly has a size of 2.45 Gb, with a scaffold N50 length of 141.46 Mb, and is anchored to 17 chromosomes. Gene annotations showed a high level of accuracy and completeness, with 23,675 predicted protein-coding genes and 94.44% of the metazoan conserved genes by BUSCO assessment. We further present 16S rRNA gene amplicon sequencing of the gut microbiota in Chaetoderma sp., which was dominated by the chemoautotrophic bacteria (phylum Gammaproteobacteria). This chromosome-level genome assembly presents the first genome for the Caudofoveata, which constitutes an important resource for studies ranging from molluscan evolution, symposium, to deep-sea adaptation.

Chromosome-level genome assembly and population genomic analysis reveal evolution and local adaptation in common hairfin anchovy (Setipinna tenuifilis).

Understanding the genetic structure and the factors associated with adaptive diversity has significant implications for the effective management of wild populations under threat from overfishing and climate change. The common hairfin anchovy (Setipinna tenuifilis) is an economically and ecologically important pelagic fish species, spanning a broad latitudinal gradient along marginal seas of the Northwest Pacific. In this study, we constructed the first reference genome of S. tenuifilis using PacBio long reads and high-resolution chromosome conformation capture (Hi-C) technology. The assembled genome was 798.38 Mb with a contig N50 of 1.43 Mb and a scaffold N50 of 32.42 Mb, which were anchored onto 24 pseudochromosomes. A total of 22,019 genes were functionally annotated, which accounted for 95.27% of the predicted protein-coding genes. Chromosomal collinearity analysis revealed chromosome fusion or fission events in Clupeiformes species. Three genetic groups of S. tenuifilis were revealed along the Chinese coast using restriction site-associated DNA sequencing (RADseq). We investigated the influence of four bioclimatic variables as potential drivers of adaptive divergence in S. tenuifilis, suggesting that these environmental variables, especially sea surface temperature, may play important roles as drivers of spatially varying selection for S. tenuifilis. We also identified candidate functional genes underlying adaptive mechanisms and ecological tradeoffs using redundancy analysis (RDA) and BayeScan analysis. In summary, this study sheds light on the evolution and spatial patterns of genetic variation of S. tenuifilis, providing a valuable genomic resource for further biological and genetic studies on this species and other closely related Clupeiformes.

The Mediator complex regulates enhancer-promoter interactions

Enhancer-mediated gene activation generally requires physical proximity between enhancers and their target gene promoters. However, the molecular mechanisms by which interactions between enhancers and promoters are formed are not well understood. Here, we investigate the function of the Mediator complex in the regulation of enhancer-promoter interactions, by combining rapid protein depletion and high-resolution MNase-based chromosome conformation capture approaches. We show that depletion of Mediator leads to reduced enhancer-promoter interaction frequencies, which are associated with a strong decrease in gene expression. In addition, we find increased interactions between CTCF-binding sites upon Mediator depletion. These changes in chromatin architecture are associated with a redistribution of the Cohesin complex on chromatin and a reduction in Cohesin occupancy at enhancers. Together, our results indicate that the Mediator and Cohesin complexes contribute to enhancer-promoter interactions and provide insights into the molecular mechanisms by which communication between enhancers and promoters is regulated.

The first high-quality chromosome-level genome of the Sipuncula Sipunculus nudus using HiFi and Hi-C data

Sipuncula is a class of exocoelomic unsegmented animals whose evolutionary relationships are unresolved. The peanut worm Sipunculus nudus is a globally distributed, economically important species belonging to the class Sipuncula. Herein, we present the first high-quality chromosome-level assembly of S. nudus based on HiFi reads and high-resolution chromosome conformation capture (Hi-C) data. The assembled genome was 1,427 Mb, with a contig N50 length of 29.46 Mb and scaffold N50 length of 80.87 Mb. Approximately 97.91% of the genome sequence was anchored to 17 chromosomes. A BUSCO assessment showed that 97.7% of the expectedly conserved genes were present in the genome assembly. The genome was composed of 47.91% repetitive sequences, and 28,749 protein-coding genes were predicted. A phylogenetic tree demonstrated that Sipuncula belongs to Annelida and diverged from the common ancestor of Polychaeta. The high-quality chromosome-level genome of S. nudus will serve as a valuable reference for studies of the genetic diversity and evolution of Lophotrochozoa.

Reorganization of three-dimensional chromatin architecture in Medicago truncatula under phosphorus deficiency.

Emerging evidence reveals that the three-dimensional (3D) chromatin architecture plays a key regulatory role in various biological processes of plants. However, information on the 3D chromatin architecture of the legume model plant Medicago truncatula and its potential roles in the regulation of response to mineral nutrient deficiency are very limited. Using high-resolution chromosome conformation capture sequencing, we identified the 3D genome structure of M. truncatula in terms of A/B compartments, topologically associated domains (TADs) and chromatin loops. The gene density, expressional level, and active histone modification were higher in A compartments than in B compartments. Moreover, we analysed the 3D chromatin architecture reorganization in response to phosphorus (P) deficiency. The intra-chromosomal cis-interaction proportion was increased by P deficiency, and a total of 748 A/B compartment switch regions were detected. In these regions, density changes in H3K4me3 and H3K27ac modifications were associated with expression of P deficiency-responsive genes involved in root system architecture and hormonal responses. Furthermore, these genes enhanced P uptake and mobilization by increasing root surface area and strengthening signal transduction under P deficiency. These findings advance our understanding of the potential roles of 3D chromatin architecture in responses of plants in general, and in particular in M. truncatula, to P deficiency.

A review and performance evaluation of clustering frameworks for single-cell Hi-C data.

The three-dimensional genome structure plays a key role in cellular function and gene regulation. Single-cell Hi-C (high-resolution chromosome conformation capture) technology can capture genome structure information at the cell level, which provides the opportunity to study how genome structure varies among different cell types. Recently, a few methods are well designed for single-cell Hi-C clustering. In this manuscript, we perform an in-depth benchmark study of available single-cell Hi-C data clustering methods to implement an evaluation system for multiple clustering frameworks based on both human and mouse datasets. We compare eight methods in terms of visualization and clustering performance. Performance is evaluated using four benchmark metrics including adjusted rand index, normalized mutual information, homogeneity and Fowlkes-Mallows index. Furthermore, we also evaluate the eight methods for the task of separating cells at different stages of the cell cycle based on single-cell Hi-C data.

MCM complexes are barriers that restrict cohesin-mediated loop extrusion

Eukaryotic genomes are compacted into loops and topologically associating domains (TADs)1–3, which contribute to transcription, recombination and genomic stability4,5. Cohesin extrudes DNA into loops that are thought to lengthen until CTCF boundaries are encountered6–12. Little is known about whether loop extrusion is impeded by DNA-bound machines. Here we show that the minichromosome maintenance (MCM) complex is a barrier that restricts loop extrusion in G1 phase. Single-nucleus Hi-C (high-resolution chromosome conformation capture) of mouse zygotes reveals that MCM loading reduces CTCF-anchored loops and decreases TAD boundary insulation, which suggests that loop extrusion is impeded before reaching CTCF. This effect extends to HCT116 cells, in which MCMs affect the number of CTCF-anchored loops and gene expression. Simulations suggest that MCMs are abundant, randomly positioned and partially permeable barriers. Single-molecule imaging shows that MCMs are physical barriers that frequently constrain cohesin translocation in vitro. Notably, chimeric yeast MCMs that contain a cohesin-interaction motif from human MCM3 induce cohesin pausing, indicating that MCMs are ‘active’ barriers with binding sites. These findings raise the possibility that cohesin can arrive by loop extrusion at MCMs, which determine the genomic sites at which sister chromatid cohesion is established. On the basis of in vivo, in silico and in vitro data, we conclude that distinct loop extrusion barriers shape the three-dimensional genome.

The genome organization of Neurospora crassa at high resolution uncovers principles of fungal chromosome topology.

The eukaryotic genome must be precisely organized for its proper function, as genome topology impacts transcriptional regulation, cell division, replication, and repair, among other essential processes. Disruptions to human genome topology can lead to diseases, including cancer. The advent of chromosome conformation capture with high-throughput sequencing (Hi-C) to assess genome organization has revolutionized the study of nuclear genome topology; Hi-C has elucidated numerous genomic structures, including chromosomal territories, active/silent chromatin compartments, Topologically Associated Domains, and chromatin loops. While low-resolution heatmaps can provide important insights into chromosomal level contacts, high-resolution Hi-C datasets are required to reveal folding principles of individual genes. Of particular interest are high-resolution chromosome conformation datasets of organisms modeling the human genome. Here, we report the genome topology of the fungal model organism Neurospora crassa at a high resolution. Our composite Hi-C dataset, which merges 2 independent datasets generated with restriction enzymes that monitor euchromatin (DpnII) and heterochromatin (MseI), along with our DpnII/MseI double digest dataset, provide exquisite detail for both the conformation of entire chromosomes and the folding of chromatin at the resolution of individual genes. Within constitutive heterochromatin, we observe strong yet stochastic internal contacts, while euchromatin enriched with either activating or repressive histone post-translational modifications associates with constitutive heterochromatic regions, suggesting intercompartment contacts form to regulate transcription. Consistent with this, a strain with compromised heterochromatin experiences numerous changes in gene expression. Our high-resolution Neurospora Hi-C datasets are outstanding resources to the fungal community and provide valuable insights into higher organism genome topology.

Exploring the effects of genetic variation on gene regulation in cancer in the context of 3D genome structure

BackgroundNumerous genome-wide association studies (GWAS) conducted to date revealed genetic variants associated with various diseases, including breast and prostate cancers. Despite the availability of these large-scale data, relatively few variants have been functionally characterized, mainly because the majority of single-nucleotide polymorphisms (SNPs) map to the non-coding regions of the human genome. The functional characterization of these non-coding variants and the identification of their target genes remain challenging.ResultsIn this communication, we explore the potential functional mechanisms of non-coding SNPs by integrating GWAS with the high-resolution chromosome conformation capture (Hi-C) data for breast and prostate cancers. We show that more genetic variants map to regulatory elements through the 3D genome structure than the 1D linear genome lacking physical chromatin interactions. Importantly, the association of enhancers, transcription factors, and their target genes with breast and prostate cancers tends to be higher when these regulatory elements are mapped to high-risk SNPs through spatial interactions compared to simply using a linear proximity. Finally, we demonstrate that topologically associating domains (TADs) carrying high-risk SNPs also contain gene regulatory elements whose association with cancer is generally higher than those belonging to control TADs containing no high-risk variants.ConclusionsOur results suggest that many SNPs may contribute to the cancer development by affecting the expression of certain tumor-related genes through long-range chromatin interactions with gene regulatory elements. Integrating large-scale genetic datasets with the 3D genome structure offers an attractive and unique approach to systematically investigate the functional mechanisms of genetic variants in disease risk and progression.

Genome organization controls transcriptional dynamics during development.

Past studies offer contradictory claims for the role of genome organization in the regulation of gene activity. Here, we show through high-resolution chromosome conformation analysis that the Drosophila genome is organized by two independent classes of regulatory sequences, tethering elements and insulators. Quantitative live imaging and targeted genome editing demonstrate that this two-tiered organization is critical for the precise temporal dynamics of Hox gene transcription during development. Tethering elements mediate long-range enhancer-promoter interactions and foster fast activation kinetics. Conversely, the boundaries of topologically associating domains (TADs) prevent spurious interactions with enhancers and silencers located in neighboring TADs. These two levels of genome organization operate independently of one another to ensure precision of transcriptional dynamics and the reliability of complex patterning processes.

Capture-C: a modular and flexible approach for high-resolution chromosome conformation capture.

Chromosome conformation capture (3C) methods measure the spatial proximity between DNA elements in the cell nucleus. Many methods have been developed to sample 3C material, including the Capture-C family of protocols. Capture-C methods use oligonucleotides to enrich for interactions of interest from sequencing-ready 3C libraries. This approach is modular and has been adapted and optimised to work for sampling of disperse DNA elements (NuTi Capture-C), including from low-cell inputs (LI Capture-C), as well as to generate Hi-C like maps for specific regions of interest (Tiled-C) and to interrogate multi-way interactions (Tri-C). We present the design, experimental protocol and analysis pipeline for NuTi Capture-C in addition to the variations for generation of LI Capture-C, Tiled-C and Tri-C data. The entire procedure can be performed in three weeks and requires standard molecular biology skills and equipment, access to a next-generation sequencing platform, and basic bioinformatic skills. Implemented with other sequencing technologies, these methods can be used to identify regulatory interactions and to compare the structural organisation of the genome in different cell types and genetic models.

A chromosomal level genome sequence for Quasipaa spinosa (Dicroglossidae) reveals chromosomal evolution and population diversity.

Quasipaa spinosa is an Asian commercial Dicroglossidae species noted for its spiny chest found in adult males. Here, we report the first chromosomal level Q. spinosa genome employing PacBio long read sequencing and high-resolution chromosome conformation capture (Hi-C) technology. The total length of the final assembled genome was 2,839,292,578bp, with contig N50 of 3.79Mb and scaffold N50 of 327.44Mb. Approximately 99.30% of the length of the assembled genome sequences were anchored to 13 chromosomes with the assistance of Hi-C reads. A total of 26,173 protein-coding genes were predicted, and 95.98% of the genes were functionally annotated. The annotated genes covered a total of 92.10% of the complete vertebrate core gene set according to the BUSCO pipeline evaluation. Approximately 41million years ago, Q. spinosa began to diverge from its dicroglossid sister taxon Nanorana parkeri. The Q. spinosa genome revealed obvious chromosomal fissions compared with Xenopus tropicalis, which probably represented a specific chromosome evolutionary history within frogs. Population analysis showed that Chinese Q. spinosa could be divided into eastern and western genetic clusters, with the western population showing higher diversity than the eastern population. The effective population size of Q. spinosa showed a continuously decreasing trend from one million years ago to 10,000years ago. In summary, this study sheds light on Q. spinosa evolution and population differentiation, providing a valuable genomic resource for further biological and genetic studies on this species, and other closely related frog taxa.

The role of nuclear organization in trans-splicing based expression of heat shock protein 90 in Giardia lamblia.

Hsp90 gene of G. lamblia has a split nature comprising two ORFs separated by 777 kb on chromosome 5. The ORFs of the split gene on chromosome 5 undergo transcription to generate independent pre-mRNAs that join by a unique trans-splicing reaction that remains partially understood. The canonical cis-acting nucleotide elements such as 5'SS-GU, 3'SS-AG, polypyrimidine tract and branch point adenine are present in the independent pre-mRNAs and therefore trans-splicing of Hsp90 must be assisted by spliceosomes in vivo. Using an approach of RNA-protein pull down, we show that an RNA helicase selectively interacts with HspN pre-mRNA. Our experiments involving high resolution chromosome conformation capture technology as well as DNA FISH show that the trans-spliced genes of Giardia are in three-dimensional spatial proximity in the nucleus. Altogether our study provides a glimpse into the in vivo mechanisms involving protein factors as well as chromatin structure to facilitate the unique inter-molecular post-transcriptional stitching of split genes in G. lamblia.

Combinatorial CRISPR inversions of CTCF sites in HOXD cluster reveal complex insulator function

CTCF (CCCTC-binding factor) is a zinc-finger protein which plays a vital role in the three-dimensional (3D) genome architecture. A pair of forward-reverse convergent CTCF binding sites (CBS elements) mediates long-distance DNA interactions to form chromatin loops with the assistance of the cohesin complex, while CBS elements at the chromatin domain boundaries show reverse-forward divergent patterns and function as insulators to discriminate against DNA interaction between chromatin domains. However, there are still many unresolved problems regarding CTCF-mediated insulation function. In order to study the connections between chromatin loops and the insulation function of CBS elements, we combinatorically inverted CBS elements at the HOXD locus by using CRISPR/Cas9 DNA-fragment editing methods in the HEK293T cell line and obtained five different kinds of single-cell CRISPR clones. By performing quantitative high-resolution chromosome conformation capture copy (QHR-4C) experiments, we found that boundary CBS inversions abolish original chromatin loops and establish new loops from the opposite direction, thus shifting the insulator boundary to the new divergent CTCF sites. Furthermore, tandem CBS elements block cohesin permeated from the opposite orientation to function as insulators. RNA-seq experiments showed that alterations of local three-dimensional genome architecture would further influence gene expression of the HOXD cluster. In conclusion, a pair of divergent CBS elements function as insulators by forming chromatin loops within chromatin domains to block cohesin sliding.

RFX5 regulates gene expression of the Pcdhα cluster.

Gene expression and three-dimensional (3D) genome organization are thought to be closely related. The protocadherin (Pcdh) clusters are central for neuronal self-avoidance in brain development and have been used as model genes to explore the role of 3D genome structures in gene regulation. Transcription factor RFX5 (regulatory factor x 5) is a member of the winged-helix subfamily of the helix-turn-helix superfamily proteins. The RFX5 protein contains four domains: oligomerization, DNA-binding, helical, and transactivation domains. RFX5 plays an important role in regulating expression of the major histocompatibility complex class II (MHC II) gene complex. Here we report a role of RFX5 in the regulation of clustered Pcdh genes. We first knocked out the RFX5 gene by using CRISPR/Cas9 DNA-fragment editing in HEC-1-B cell line. By RNA-seq experiments, we found that RFX5 deletion results in a significant increase of expression levels of the Pcdhα6, Pcdhα12 and Pcdhαc2 genes. By ChIP-nexus experiments, we found that RFX5 deletion leads to the increased binding of CTCF and RAD21 in the Pcdhα cluster. Finally, through quantitative high-resolution chromosome conformation capture copy (QHR-4C) experiments, we found that RFX5 deletion results in a significant increase of long-distance chromatin interactions between the HS5-1 enhancer and its target promoters in the Pcdhα cluster. In conclusion, RFX5 regulates gene expression of the Pcdhα cluster through higher-order chromatin structure.

Computational methods for predicting 3D genomic organization from high-resolution chromosome conformation capture data.

The advent of high-resolution chromosome conformation capture assays (such as 5C, Hi-C and Pore-C) has allowed for unprecedented sequence-level investigations into the structure–function relationship of the genome. In order to comprehensively understand this relationship, computational tools are required that utilize data generated from these assays to predict 3D genome organization (the 3D genome reconstruction problem). Many computational tools have been developed that answer this need, but a comprehensive comparison of their underlying algorithmic approaches has not been conducted. This manuscript provides a comprehensive review of the existing computational tools (from November 2006 to September 2019, inclusive) that can be used to predict 3D genome organizations from high-resolution chromosome conformation capture data. Overall, existing tools were found to use a relatively small set of algorithms from one or more of the following categories: dimensionality reduction, graph/network theory, maximum likelihood estimation (MLE) and statistical modeling. Solutions in each category are far from maturity, and the breadth and depth of various algorithmic categories have not been fully explored. While the tools for predicting 3D structure for a genomic region or single chromosome are diverse, there is a general lack of algorithmic diversity among computational tools for predicting the complete 3D genome organization from high-resolution chromosome conformation capture data.

Charting the cis-regulome of activated B cells by coupling structural and functional genomics.

Cis-regulomes underlying immune-cell-specific genomic states have been extensively analyzed by structure-based chromatin profiling. By coupling such approaches with a high-throughput enhancer screen (self-transcribing active regulatory region sequencing (STARR-seq)), we assembled a functional cis-regulome for lipopolysaccharide-activated B cells. Functional enhancers, in contrast with accessible chromatin regions that lack enhancer activity, were enriched for enhancer RNAs (eRNAs) and preferentially interacted in vivo with B cell lineage-determining transcription factors. Interestingly, preferential combinatorial binding by these transcription factors was not associated with differential enrichment of their sites. Instead, active enhancers were resolved by principal component analysis (PCA) from all accessible regions by co-varying transcription factor motif scores involving a distinct set of signaling-induced transcription factors. High-resolution chromosome conformation capture (Hi-C) analysis revealed multiplex, activated enhancer-promoter configurations encompassing numerous multi-enhancer genes and multi-genic enhancers engaged in the control of divergent molecular pathways. Motif analysis of pathway-specific enhancers provides a catalog of diverse transcription factor codes for biological processes encompassing B cell activation, cycling and differentiation.

PF153 MLL‐AF4 CAUSES ABERRANT UPREGULATION OF PROM1 (CD133) IN ACUTE LYMPHOBLASTIC LEUKEMIA BY CONTROLLING ENHANCER‐PROMOTER INTERACTIONS

Background:MLL rearrangements (MLLr) are a common cause of aggressive, incurable acute lymphoblastic leukaemias (ALL) in infants and children, most of which originate in utero. The most common MLLr produces an MLL‐AF4 fusion protein. MLL‐AF4 is thought to promote leukaemogenesis by activating key target genes, mainly by altering the epigenetic profile of the cell through recruitment of DOT1L and increasing histone H3 lysine 79 methylation (H3K79me2/3). One of these key genes is PROM1 which encodes the cell surface protein CD133 (Prominin‐1; Mak2016). CD133 has previously been shown to be expressed on ALL blasts and is of potential therapeutic interest. However, the mechanism by which MLL‐AF4 dysregulates PROM1 expression and the impact of PROM1 expression on leukaemia cell survival remains unclear.Aims:To determine the mechanism by which MLL‐AF4 causes aberrant expression of PROM1 in leukaemia.Methods:We characterised normal fetal haematopoietic stem and progenitor cells (HSPC), blasts from patient leukaemia samples, and SEM cells (a MLL‐AF4+ cell line) by immunophenotyping and transcriptome analysis. ChIP‐Seq and Capture C (a high resolution chromosome conformation capture technique) was used to elucidate enhancer‐promoter structure at the PROM1 locus.Results:We find that the MLL‐AF4 complex aberrantly upregulates PROM1 transcription by controlling enhancer promoter‐interactions.Immunophenotyping and transcriptome analysis: Immunophenotypic analysis of normal human fetal bone marrow HSPC confirmed that CD133 expression is confined to the stem cell compartment and is downregulated through B cell commitment. In contrast, MLLr patient blasts and SEM cells showed high levels of surface CD133 expression. Analysis of RNA‐Sequencing data showed that PROM1 gene expression in normal fetal B progenitors, patient blasts and SEM cells reflected the pattern of surface CD133 expressionFunctional assays: Upon inhibition of MLL‐AF4 binding by a DOT1L inhibitor, CD133 expression is downregulated in SEM cells and nascent RNA‐Seq shows down‐regulation of PROM1 at the transcriptional level. Direct knock down of PROM1 in SEM cells confers slower growth rates/less colony formation and increased cell death.ChIP‐Seq and Capture‐C: ChIP‐sequencing identified MLL‐AF4 and associated complex components binding to the PROM1 promoter and within the gene body at putative enhancer elements. Capture C profiles identified intragenic enhancer elements within the PROM1 gene body and within the nearby TAPT1 gene. Inhibition of MLL‐AF4 complex activity upon treatment with DOT1L inhibitors causing loss of H34K79me2/3 disrupted gene expression as well as enhancer‐promoter interactions.Summary/Conclusion:Together these results indicate that MLL‐AF4 directly up‐regulates PROM1 expression in ALL cells and promotes their proliferation and survival. Therefore, CD133 expression should not only be a useful marker for diagnosis and residual disease monitoring but might also be exploited as a therapeutic target in MLL re‐arranged ALL.

Long-range interactions between topologically associating domains shape the four-dimensional genome during differentiation.

Genomic information is selectively used to direct spatial and temporal gene expression during differentiation. Interactions between topologically associating domains (TADs) and between chromatin and the nuclear lamina organize and position chromosomes in the nucleus. However, how these genomic organizers together shape genome architecture is unclear. Here, using a dual-lineage differentiation system, we report long-range TAD-TAD interactions that form constitutive and variable TAD cliques. A differentiation-coupled relationship between TAD cliques and lamina-associated domains suggests that TAD cliques stabilize heterochromatin at the nuclear periphery. We also provide evidence of dynamic TAD cliques during mouse embryonic stem-cell differentiation and somatic cell reprogramming and of inter-TAD associations in single-cell high-resolution chromosome conformation capture (Hi-C) data. TAD cliques represent a level of four-dimensional genome conformation that reinforces the silencing of repressed developmental genes.

A tissue-specific self-interacting chromatin domain forms independently of enhancer-promoter interactions

Self-interacting chromatin domains encompass genes and their cis-regulatory elements; however, the three-dimensional form a domain takes, whether this relies on enhancer–promoter interactions, and the processes necessary to mediate the formation and maintenance of such domains, remain unclear. To examine these questions, here we use a combination of high-resolution chromosome conformation capture, a non-denaturing form of fluorescence in situ hybridisation and super-resolution imaging to study a 70 kb domain encompassing the mouse α-globin regulatory locus. We show that this region forms an erythroid-specific, decompacted, self-interacting domain, delimited by frequently apposed CTCF/cohesin binding sites early in terminal erythroid differentiation, and does not require transcriptional elongation for maintenance of the domain structure. Formation of this domain does not rely on interactions between the α-globin genes and their major enhancers, suggesting a transcription-independent mechanism for establishment of the domain. However, absence of the major enhancers does alter internal domain interactions. Formation of a loop domain therefore appears to be a mechanistic process that occurs irrespective of the specific interactions within.