Targeted Chromatin Capture (T2C): a novel high resolution high throughput method to detect genomic interactions and regulatory elements
© Kolovos et al.; licensee BioMed Central Ltd. 2014
Received: 28 April 2014
Accepted: 28 May 2014
Published: 16 June 2014
Significant efforts have recently been put into the investigation of the spatial organization and the chromatin-interaction networks of genomes. Chromosome conformation capture (3C) technology and its derivatives are important tools used in this effort. However, many of these have limitations, such as being limited to one viewpoint, expensive with moderate to low resolution, and/or requiring a large sequencing effort. Techniques like Hi-C provide a genome-wide analysis. However, it requires massive sequencing effort with considerable costs. Here we describe a new technique termed Targeted Chromatin Capture (T2C), to interrogate large selected regions of the genome. T2C provides an unbiased view of the spatial organization of selected loci at superior resolution (single restriction fragment resolution, from 2 to 6 kbp) at much lower costs than Hi-C due to the lower sequencing effort.
We applied T2C on well-known model regions, the mouse β-globin locus and the human H19/IGF2 locus. In both cases we identified all known chromatin interactions. Furthermore, we compared the human H19/IGF2 locus data obtained from different chromatin conformation capturing methods with T2C data. We observed the same compartmentalization of the locus, but at a much higher resolution (single restriction fragments vs. the common 40 kbp bins) and higher coverage. Moreover, we compared the β-globin locus in two different biological samples (mouse primary erythroid cells and mouse fetal brain), where it is either actively transcribed or not, to identify possible transcriptional dependent interactions. We identified the known interactions in the β-globin locus and the same topological domains in both mouse primary erythroid cells and in mouse fetal brain with the latter having fewer interactions probably due to the inactivity of the locus. Furthermore, we show that interactions due to the important chromatin proteins, Ldb1 and Ctcf, in both tissues can be analyzed easily to reveal their role on transcriptional interactions and genome folding.
T2C is an efficient, easy, and affordable with high (restriction fragment) resolution tool to address both genome compartmentalization and chromatin-interaction networks for specific genomic regions at high resolution for both clinical and non-clinical research.
KeywordsChromatin conformation capture Long range interactions Enhancers Promoters
A number of recent studies have shown that the genome is organized in self-associating domains  that are separated by linker regions. These so-called ‘topological domains’ or ‘topological associated domains’ generally range from 300 kilobasepairs (kbp) to 1 megabasepairs (1 Mb) and consist of a series of different types of chromatin loops in agreement with earlier models of the genome ( and references therein).
One loop is defined as two distant chromatin regions coming, spatially, into close proximity (interact with each other), thereby creating DNA loops. Such ‘long-range interactions’ have been first observed between promoters and distant enhancers ([3, 4] and references therein) and can bring DNA elements together that are separated by a large distance on the linear DNA strand ([5, 6] and references therein). These regulatory elements (enhancers or silencers) are short sequences containing several binding sites for transcription factors, which regulate the activation (reviewed in ) repression (reviewed in ) genes and their subsequent transcription (reviewed in ). In the linear genome the distance between enhancer(s) and gene can be quite large, for example, the sonic hedgehog (shh) enhancer is located about 1 Mb away from its target gene Shh. Changes or differences within these elements and their interaction with genes can be responsible for changes in gene expression , causing intrinsic differences between individuals, disease susceptibility, and disease progression.
A number of chromatin loops are thought to be purely structural, that is, to enable the folding of the genome creating distinct topological domains, while other loops have a function in the expression of genes. Loops of the latter type are frequently found within topological domains, but are less frequently observed between different topological domains [1, 12]. These regulatory chromatin loops change and depend on a large number of proteins including Ctcf , cohesin , and a series of transcription factors [15–18], which are mostly involved in the transcriptional regulation of genes within the domain.
Comparison between different chromatin conformation capturing techniques (adopted and modified from)
Laborious, requires knowledge of the locus and proper controls
Good resolution, good signal-to-noise ratio
Restricted to single viewpoint per experiment when multiplexing several viewpoints, analysis requires extra bioinformatics expertise, not an all-to-all genome-wide method
Relatively simple data analysis
Poor signal-to-noise ratio, difficult to obtain genome-wide coverage
Identifies interactions between many individual fragments
Very laborious, no genome-wide coverage, primer design can be challenging. Analysis requires advanced bioinformatics expertise
Explores the genome-wide interactions between all individual fragments
Very expensive, requires a large sequence effort to obtain sufficient coverage, approximately 10 to 40 kbp resolution, requires advanced bioinformatics expertise
Explores the interactome of a selected region in cis but also in trans, high (restriction fragment) resolution, cheaper than Hi-C and 5C, requiring only half a lane of Illumina HiSeq2000
Is restricted to the selected regions of the genome, requires advanced bioinformatics expertise
The analysis of the interactions of several viewpoints with the aforementioned techniques in 3C and 4C is possible, but the choice for several viewpoints will increase the costs and work effort linearly. However, the number of viewpoints can also be limited due to the (often) limiting amount of available cell material. 5C is demanding in primer design and allows the analysis of interactions only among the primer designed fragments. Furthermore, genome-wide coverage is not possible. Hi-C is very expensive as it requires extremely deep sequencing in order to cover the whole genome, even at a relatively low resolution of 40 kbp. The most recent Hi-C data analysis has used a new algorithm and provided a genome-wide interaction map of 10 kbp resolution. However, an enormous amount of sequencing is required (3.4 billion mapped paired-end reads from six biological replicates) . Such effort is not affordable for most research groups and, in addition, the scientific interest is most of the time focused on a specific question involving a limited set of specific loci or domains. Hence, there is a need for a technique which eliminates most of the aforementioned limitations.
Here we present Targeted Chromatin Capture (T2C), a new 3C method, which does not involve a massive sequencing effort, but which results in a high resolution map of interactions for particular loci of interest. We used the well-studied human H19/IGF2 locus and compared the results of our new method with data from other chromatin conformation capturing techniques. Using the mouse β-globin locus we demonstrated that the method can reliably identify chromatin structural changes between different tissues and also allows the study of the role of individual transcription factors in the chromatin architecture.
Overview of the procedure
To overcome the aforementioned problems of the 5C and Hi-C techniques we have developed the novel method T2C. The method has the advantage that it allows the analysis of the structure of the genome and all the interactions of selected regions of the genome at high resolution (single restriction fragments) without a massive sequencing effort and associated costs.
T2C employs a selective enrichment of the 3C ligation products in preselected regions of interest in order to identify their interactions within a domain as well as the compartmentalization of one or several specific regions of the genome. These regions can be continuous Mb sized genomic regions, but could also be a collection of smaller regions (a few kbp each). Every captured restriction fragment can be used as a single ‘4C-seq viewpoint’ and analyzed accordingly. The results of T2C provide a local interaction map at a restriction fragment-level resolution accompanied with a lower sequencing effort and less intricate bioinformatics analysis than Hi-C. T2C also overcomes the limits of 5C since it identifies not only interactions within the targeted region(s), but also interactions between the targeted region(s) and with regions outside of them.
In brief, we have designed sets of unique oligonucleotide probes (ranging from 62 to 90 nucleotides) specific for all the restriction fragments and as close as possible to the end of the first restriction site (Mm - HindIII + NlaIII digest, Hs - BglII + NlaIII digest) in our regions of interest, the mouse β-globin locus and the human H19/IGF2 locus (see Methods). Alternative to continuous regions, separate genomic regions within one (or more) chromosomes could be analyzed simultaneously. The oligonucleotides are spotted on an array or can alternatively be captured on beads. Some fragment ends cannot be captured by a designed oligonucleotide due to the presence of repeat elements or the insufficient size of the restriction fragment end. Repetitive sequences are a general problem in all 3C-based methods, including Hi-C. The size limitation of the fragment end can be circumvented if necessary by a backup procedure with different enzymes (changing either the first or the second restriction enzyme or both), which generates a new set of end fragments or by mechanically shearing of the chromatin (instead of the second restriction enzyme digestion) which can result in fragment sizes of different length (see Discussion).
Summary of information about the different experiments
Genome assembly version
Coordinates oligo-nucleotide positions
Size of area of interest (Mb)
Median resolution (kbp)
Raw paired reads (n)
Paired reads that could be mapped to the whole genome (n)
Mapped paired reads between the region of interest and the whole genome (n)
Uniquely mapped paired-reads in the whole genome without self-ligation and and non-digestion (n)
Uniquely mapped paired-reads between the region of interest and the whole genome without self-ligation and non-digestion (n)
Uniquely mapped paired-reads inside the region of interest without self-ligation and non-digestion (n)
‘Interactions’ inside the region of interest (n)
Average number of reads/interaction in the region of interest (n)
Mouse fetal liver
Mouse fetal brain
chr11: 1100646 - 3173091
T2C identifies known long-range interactions
We selected unique oligonucleotides mapping near the ends of 344 BglII generated fragments spanning 2.1 Mb around the H19/IGF2 locus (Table 2). This set of 525 oligonucleotides was spotted on a capture array. A ligation fragment library was generated from the breast endothelial cell line 1-7HB2 (abbreviated HB2) after digestion with BglII and NlaIII according to the 3C-seq protocol  (see also Figure 1). The library was subsequently hybridized to the capture array. After elution from the capture array the captured DNA fragments were amplified by a PCR with low cycle number (12 cycles) and sequenced by paired-end Illumina sequencing (see Methods).
To demonstrate that T2C reveals a similar overall interaction pattern and compartmentalization of the locus as observed by Hi-C in IMR90 cells  we first binned the paired-reads into 40 kbp bins (Figure 2A, B). The interaction patterns at this level of resolution show that the topological domain is maintained between different cell types, HB2 versus IRM90  with a Spearman’s rank correlation coefficient r s = 0.64 (P <2.2 × 10-16).
However, with T2C we obtained a chromatin interaction map at restriction fragment resolution (Figure 2C, each block represents one restriction fragment), revealing significantly more detail with respect to the general chromatin organization of the region when visualized by a logarithmic and rainbow-like colored interaction frequency. To first validate T2C in comparison to 3C and 4C-seq we compared the interactions of a single restriction fragment (CTCF AD viewpoint)  to interactions detected for this fragment by 3C  and 4C-seq  (Figure 2D, E, F). Although there are some variations in the read coverage of the individual interactions, similar interactions can be observed by both 4C-seq and T2C. Moreover, both methods detect interactions which we previously observed with 3C . It should be noted that an important difference between 4C-seq and T2C is the number of PCR amplification cycles. For T2C this is on average 12 cycles (only after capture) whereas for 4C-seq it is 30 cycles. The lower number of cycles will give less PCR bias of the different fragments relative to each other, because fragments have different PCR efficiencies.
We conclude that the T2C method yields interaction data at a resolution identical to 4C-seq for the individual restriction fragments (median approximately 4 kbp resolution) and that when T2C is performed for a continuous region over 2 Mb it can reproduce the overall topological domain structure that was observed by Hi-C.
T2C identifies different interaction networks based on different biological materials
The mouse β-globin locus undergoes structural changes upon activation in erythroid tissue [20, 34, 35], but is surrounded by silent olfactory receptor genes, which are only expressed in the olfactory epithelium. The major difference between the H19/IGF2 locus and the β-globin locus is that the β-globin locus is embedded in a large area of inactive genes. Thus two patterns of interactions may be expected in erythroid cells, those important for the globin locus and those present in inactive chromatin. We selected a region of 2.1 Mb around the locus (Table 2) containing 719 restriction fragments of the restriction enzyme HindIII (6 bp recognition site). About 800 oligonucleotide probes were designed close to the ends of the fragments. To analyze the locus in its active state we used primary erythroid cells from fetal liver which were compared to fetal brain cells as a model of inactive loci. Based on results from previous 3C studies of the locus [20, 35] we expected in primary erythroid cells a higher number of interactions around the β-globin gene and between the β-globin gene and its regulatory elements.
The analysis of the hybridized fragments shows that almost the entire 2.1 Mb appears to be part of one topological domain (with two possible subdomains, one of which contains the β-globin locus) with the next domain starting near the end of the selected sequences (due to the repetitive sequences and the borders of the region of interest, that topological domain cannot be depicted clearly, in agreement with Dixon et al.) both in mouse primary erythroid cells (Figure 3A, right hand side) and mouse fetal brain cells (Figure 3B) with many interactions within the topological domain (Figure 3C and 3D). Although the topological domain structure between the different biological materials is similar, there appear to be fewer interactions in mouse fetal brain cells relative to mouse primary erythroid cells due to the inactivity of the locus in the brain (Figure 3). Focusing on the β-globin region, all the well-known interactions in the β-globin locus are detected in the primary erythroid cells. The known interactions, such as between the β-globin promoter and Locus Control Region (LCR) (Figure 4B, adapted and modified from Drissen et al., with blue line depicting the interactions for primary erythroid cells and with grey the interactions for mouse fetal brain cells) and between the LCR-3′HS1 are clearly visualized [16, 20, 35] (Figure 4A). These interactions are absent from the fetal brain sample (Figure 4C). Furthermore, the main regulatory region (HS1-6) shows the well-known interaction with the β-globin genes and HS1 at the 3′end of the locus in fetal liver cells but not in brain [16, 20]. In addition, for the β-globin promoter we identify a few additional interactions further away than the ones previously reported. These are located even approximately 1 Mb far from the β-globin promoter (Figure 3A). It is unknown whether these interactions are related to the functioning of the β-globin genes or whether these DNA elements are in close proximity due to the folding of the domain, although their absence in the fetal brain suggests they have a role in the regulation of the globin β-globin. In addition to the interactions in cis, the β-globin (Hbb-b1) gene and the LCR also contact a number of positions on other chromosomes.
T2C in combination with ChIP-seq identifies factor specific interactions
The importance of the role of chromatin interactions in the regulation of the gene transcription is well established [9, 39–42]. However, there is still an increasing need for a quick, easy, and affordable technique to provide the information on chromatin interactions and the compartmentalization of the genome. T2C is affordable to most scientific groups and will meet in a satisfactory manner their needs for detecting high resolution chromatin organization of selected loci. Every restriction fragment can serve as a ‘viewpoint’ and all their interactions, either short or long or to other chromosomes (data not shown), can be identified. Thus multiple 3C-seq, 4C-seq or 5C experiments do not have to be performed. Moreover, with T2C the compartmentalization of the genome can be identified in the regions of interest without requiring the large sequencing effort of Hi-C, which would increase the costs tremendously. Furthermore, due to the T2C design, a better coverage and resolution of the locus is obtained when compared to other genome wide techniques (like Hi-C and 3C with its derivatives) using a 6 bp cutter as first restriction enzyme. Here we multiplexed two samples, but by multiplexing more than two samples the costs are likely to be reduced significantly without sacrificing the quality of the output. We have recently successfully used 13 samples per sequencing lane, including the β- globin locus which showed the same interactions (data not shown).
The resolution of T2C is based on the restriction enzyme used. Digesting cross-linked chromatin from primary erythroid cells and HB2 cells with HindIII or BglII, resulted in a median resolution of 2 kbp and 4.1 kbp, respectively (Table 2). That provides a significantly better resolution than the usual 40 kbp bins obtained with Hi-C. Furthermore, comparing T2C with 4C-seq  and Hi-C  for the H19/IGF2 locus (Figure 2) and with already published 3C-qPCR data for the β-globin locus [16, 20, 35], the same topological domains and chromatin interaction networks were identified. Taken together, all these results reveal the strengths of the T2C as a tool to identify all the interactions and the compartmentalization of specific regions of the genome.
In addition, the T2C interactions are easily connected to the factors that play a role in these interactions or the type of elements (promoters/enhancers) involved in the interactions. Ldb1 and Ctcf are important proteins which mediate chromatin interactions. Ldb1 is an important transcription factor necessary for primitive mouse hematopoiesis and for the development of megakaryocytes [43, 44] and controls essential hematopoietic pathways in mouse early development . Depletion of Ldb1 is lethal for mouse embryos after E9.5 with severe effects such as impairment of hematopoietic and vascular development . It is well established that the LCR has higher interaction frequencies with the β-globin locus in mouse primary erythroid cells comparing to mouse brain cells [16, 20, 35] and that Ldb1 is significantly enriched in the LCR region in mouse primary erythroid cells relative to mouse fetal brain cells  (Figure 5E vs. Figure 5F). Furthermore, Ctcf is an insulator binding protein known to be involved in chromatin conformation  and is enriched at the boundaries of topological domains . Ctcf mediates long range interactions in the β-globin locus  (Figure 5C vs. Figure 5D and Figure 5G vs. Figure 5H). Hence, it is no surprise that for Ldb1 and Ctcf occupied restriction fragments we observe a higher number of interacting fragments at larger linear distances of fragments that interact in mouse primary erythroid cells than in mouse brain cells (Figure 6). This effect can be explained by the fact that the β-globin locus is active in mouse primary erythroid cells. Furthermore, we observe that the boundaries of the topological domain, which contains the β-globin locus, are easily observed in mouse erythroid cells (Figure 3A). That is prominent when depicting only the Ctcf interacting fragments (Figure 5C vs. Figure 5D). Furthermore, the number of interactions within that topological domain, appear higher in the erythroid cells comparing to fetal brain cells (Figure 3A vs. Figure 3B, Figure 6A, B). We hypothesize that this is due to the fact that the β -globin locus is active with open chromatin in mouse primary erythroid cells. Hence, the chromatin has a different conformation by enabling the interaction between many different elements necessary for the regulation of the gene . However, in mouse fetal brain cells, where β-globin locus is not active, that is not necessary and there are no important elements that need to spatially be in close proximity.
Because T2C is focused on particular regions of interest, it would be easy to design a set of oligonucleotides for a number of loci that are known to be associated with a particular disease and design a diagnostic kit on that basis that could handle many samples at the same time. Since SNPs are often linked to diseases, dedicated oligonucleotides for them can be designed in order to assess their effect in long range interactions and the regulation of the gene transcription. For non-clinical research purposes the size of the region used in our experiments is sufficient (more than 2 Mb) to extract safe conclusions about the local chromatin interactome and the compartmentalization of the genome.
We conclude that T2C can be used as an affordable, cost-effective, diagnostic tool with single restriction fragment resolution to explore the local spatial organization of the genome and chromatin interactions without requiring laborious procedures or massive sequencing efforts.
A microarray for the β-globin locus was designed containing unique oligonucleotides and physically as close as possible to the HindIII restriction sites spanning 2.1 Mb around the gene (chr7: 109876329-111966581, mm9). For the H19/IGF2 locus unique oligonucleotides were designed close to BglII restriction sites (chr11:1100646-3173091, hg18) spanning an area of 2.1 Mb (Table 2).
The oligonucleotides were designed with the following criteria, they should be: (1) as close as possible to the first restriction site; (2) a unique DNA sequence within the area of interest and preferably in the entire genome; (3) similar melting temperatures, but with different base composition and the length; (4) oligonucleotides which exceed the second restriction site due to very small end fragments, were trimmed keeping in mind to stay close to the same melting temperature.
A custom-made NimbleGen Sequence Capture 2.1 M capture array is produced separately for the H19/IGF2 locus and for the β-globin locus containing for each one the oligonucleotides which satisfy the aforementioned criteria. The oligonucleotides, 525 for the H19/IGF2 locus and 800 for the β-globin locus, were replicated proportionally and equally up to 2.1 M in total for each design, that is, for the β-globin locus each of the 800 oligonucleotides was spotted in 2,625 spots.
Chromatin isolation and library preparation
Nuclei from approximately 107 mouse primary erythroid cells from mouse fetal liver E12.5, mouse fetal brain cells E12.5, and a human breast endothelial cell line (HB2) were isolated, cross-linked (in 2% formaldehyde at room temperature) quenched with 1 M glycine and were re-suspended in lysis buffer (10 mM Tris–HCl (pH 8.0), 10 mM NaCl, 0.2% (vol/vol) NP-40 and 1× protease inhibitor solution). The chromatin was digested with a 6-cutter (400 units of Hin dIII for mouse cells and Bgl II for the HB2 cells) and ligated using 100 units of T4 DNA ligase (Promega) under conditions favoring intramolecular ligation events. After reversing the cross-link at 65°C overnight, 50 μg of the resulting DNA chromatin library were digested with a frequent 4-cutter (Dpn II or Nla III for the mouse cells, Nla III for the HB2 cells, at a DNA concentration of 100 ng/μL, using 1 unit of enzyme per μg of DNA). All these steps were performed according to the initial steps of 3C-seq protocol, as described previously .
The final library is prepared for analysis on the Illumina Cluster Station and HiSeq 2000 Sequencer according to the Illumina TruSeq DNA protocol with modifications (http://www.illumina.com). In short, the digested library is purified using AMPure XP beads (Beckman Coulter), end-repaired, and cleaned using AMPure XP beads. The now blunt-ended fragments were A-tailed using the Klenow exo enzyme in the presence of ATP and purified again using AMPure XP beads. Then indexed adapters provided by Illumina were ligated to the A-tailed DNA fragments with subsequent purification using AMPure XP beads.
The resulting adapter-modified DNA library (300 to 500 ng) was hybridized in 35 μL for 64 h at 42°C on a custom-made NimbleGen Sequence Capture 2.1 M capture array according to NimbleGen Sequence Capture array protocol (http://www.nimblegen.com/seqcapez) on the NimbleGen Hybridization System. The captured DNA fragments are eluted from the capture array and purified using MinElute columns (Qiagen). The yield for a positive region (a fragment inside the region of interest) and a negative region (a fragment outside the region of interest) differ by >30-fold on average. The captured DNA fragments are amplified by 12 PCR cycles. PCR products are purified using AMPure XP beads and eluted in 30 μL of re-suspension buffer. One microliter is loaded on an Agilent Technologies 2100 Bioanalyzer using a DNA 1000 assay to determine the library concentration and to check for quality.
Cluster generation and high throughput sequencing
Cluster generation is performed according to the Illumina Cluster Reagents preparation protocol (http://www.illumina.com). Briefly, 1 μL of a 10 nM TruSeq DNA library stock DNA is denatured with NaOH, diluted to 9-10 pM and hybridized onto the flowcell. The hybridized fragments are sequentially amplified, linearized, and end-blocked according to the Illumina Paired-end Sequencing user guide protocol. After hybridization of the sequencing primer, sequencing by synthesis is performed using the HiSeq 2000 sequencer with a 101 cycle protocol according to manufacturer’s instructions. The sequenced fragments were denatured with NaOH using the HiSeq 2000 and the index-primer was hybridized onto the fragments. The index was sequenced with a seven-cycle protocol. The fragments are denatured with NaOH, sequentially amplified, linearized, and end-blocked. After hybridization of the sequencing primer, sequencing by synthesis of the third read is performed using the HiSeq 2000 sequencer with a 101-cycle protocol.
Targeted Chromatin Capture data analysis
The generated HiSeq 2000 sequencing reads were trimmed if the reads contained the first enzyme restriction recognition site (Hin dIII for the mouse derived reads and Bgl II for the human derived reads) For each read with one or more enzyme recognition sites, the DNA sequence after the 3′ end of the first site was removed, that is, after the trimming procedure the trimmed reads contained and ended with a single restriction recognition site. Subsequently, consecutive bases with a quality score lower than 10 were cut off from the ends of all the reads and the reads that contained less than 12 bases were omitted using Trimmomatic . We used the Burrows-Wheeler Alignment tool (BWA, version 0.6.1) to the whole genome NCBI36/hg18 assembly for the human derived reads and to NCBI37/mm9 assembly for the mouse derived reads, using default settings . Aligned reads that localized between two second enzyme recognition sites that did not contain a first enzyme recognition site, that is, all Nla III-Nla III restriction fragments were removed using BEDtools .
In the alignment, paired reads were removed if one of the reads was not uniquely mapped. Furthermore, paired reads that were a result of a self-ligation event, non-digestion/re-ligation event, or a ligation of identical ends were removed from the analysis, since these paired reads introduce a common bias in chromosome conformation capture techniques [50, 51]. The alignments were further processed with SAMtools  to generate paired-end Binary Alignment/Map (BAM) files. BEDtools  was used to remove reads that overlapped more than one restriction fragment. Interaction matrices were generated from the alignments at a resolution of the restriction fragments and at 40 kb resolution (using BEDtools on a 40 kbp binned genome). In addition, the human T2C 40 kb binned data were compared to IMR90 40 kb Hi-C data of the combined replicates . The T2C interaction plots were normalized for capture efficiency of the fragments. For each interaction the number reads of each interaction was normalized through dividing it by the sum of the reads of both fragments involved in the interaction. Similarly, the T2C plots of the 40 kb bins were normalized after all the fragments were divided into 40 kb bins along each chromosome. ChIP-seq and T2C interaction-intersection plots were generated from normalized T2C interaction plots and intersected with fragments that contained a ChIP-seq peak signal of the protein of interest. The statistical software package R (version 3.1.0) was used to generate the interaction plots and to conduct the statistical calculations .
Published ChIP-seq datasets [37, 38] were obtained and analyzed. MACS  was used to identify peaks (fdr ≤0.01, peak height ≥20 overlapping reads) to intersect their positions with the interacting fragments obtained from T2C.
We thank the members of the FG Laboratory, Argyris Papantonis, Robert-Jan Palstra, and Danny Huylebroeck for discussions and reading the manuscript. PK was supported by grants from EpiGenSys/ERASysBio +/FP7 (NL: NWO, UK: BSRC, D: BMBF). JZ was supported by an NWO ALW grant and KSW by E-RARE/TARGET-CdLS (NL: ZonMW). NK and the grid infrastructure were supported by the BMBF (grant #01AK803A (German MediGRID), and #01IG07015G (Services@MediGRID). HJGvdW was supported by Zenith (93511036) grant from the Netherlands Genomics Initiative (NGI). The work was supported by EpiGenSys/ERASysBio +/FP7 (NL: NWO, UK: BSRC, D: BMBF), the Bluescript EU Integrated Project, the Netherlands Institute for Regenerative Medicine (NIRM), the MEC Booster grant by the Netherlands Genomics Institute (MEC Booster grant) and the KNAW (professorship to FG).
The accession number is SRP042002.
- Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, Hu M, Liu JS, Ren B: Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012, 485: 376-380. 10.1038/nature11082.PubMed CentralView ArticlePubMedGoogle Scholar
- Cremer T, Cremer C: Chromosome territories, nuclear architecture and gene regulation in mammalian cells. Nat Rev Genet. 2001, 2: 292-301. 10.1038/35066075.View ArticlePubMedGoogle Scholar
- Dillon N, Trimborn T, Strouboulis J, Fraser P, Grosveld F: The effect of distance on long-range chromatin interactions. Mol Cell. 1997, 1: 131-139. 10.1016/S1097-2765(00)80014-3.View ArticlePubMedGoogle Scholar
- Mueller-Storm HP, Sogo JM, Schaffner W: An enhancer stimulates transcription in trans when attached to the promoter via a protein bridge. Cell. 1989, 58: 767-777. 10.1016/0092-8674(89)90110-4.View ArticlePubMedGoogle Scholar
- Jhunjhunwala S, van Zelm MC, Peak MM, Cutchin S, Riblet R, van Dongen JJ, Grosveld FG, Knoch TA, Murre C: The 3D structure of the immunoglobulin heavy-chain locus: implications for long-range genomic interactions. Cell. 2008, 133: 265-279. 10.1016/j.cell.2008.03.024.PubMed CentralView ArticlePubMedGoogle Scholar
- Medvedovic J, Ebert A, Tagoh H, Tamir IM, Schwickert TA, Novatchkova M, Sun Q, Veld PJHI’t, Guo C, Yoon HS, Denizot Y, Holwerda SJ, de Laat W, Cogne M, Alt FW, Busslinger M: Flexible long-range loops in the VH gene region of the Igh locus facilitate the generation of a diverse antibody repertoire. Immunity. 2013, 39: 229-244. 10.1016/j.immuni.2013.08.011.View ArticlePubMedGoogle Scholar
- Bulger M, Groudine M: Functional and mechanistic diversity of distal transcription enhancers. Cell. 2011, 144: 327-339. 10.1016/j.cell.2011.01.024.PubMed CentralView ArticlePubMedGoogle Scholar
- Maeda RK, Karch F: Gene expression in time and space: additive vs hierarchical organization of cis-regulatory regions. Curr Opin Genet Dev. 2011, 21: 187-193. 10.1016/j.gde.2011.01.021.View ArticlePubMedGoogle Scholar
- Kolovos P, Knoch TA, Grosveld FG, Cook PR, Papantonis A: Enhancers and silencers: an integrated and simple model for their function. Epigenetics Chromatin. 2012, 5: 1-10.1186/1756-8935-5-1.PubMed CentralView ArticlePubMedGoogle Scholar
- Maas SA, Fallon JF: Single base pair change in the long-range Sonic hedgehog limb-specific enhancer is a genetic basis for preaxial polydactyly. Dev Dyn. 2005, 232: 345-348. 10.1002/dvdy.20254.View ArticlePubMedGoogle Scholar
- Lin YC, Murre C: Nuclear location and the control of developmental progression. Curr Opin Genet Dev. 2013, 23: 104-108. 10.1016/j.gde.2012.11.010.PubMed CentralView ArticlePubMedGoogle Scholar
- Splinter E, de Wit E, Nora EP, Klous P, van de Werken HJ, Zhu Y, Kaaij LJ, van Ijcken W, Gribnau J, Heard E, de Laat W: The inactive X chromosome adopts a unique three-dimensional conformation that is dependent on Xist RNA. Genes Dev. 2011, 25: 1371-1383. 10.1101/gad.633311.PubMed CentralView ArticlePubMedGoogle Scholar
- Splinter E, Heath H, Kooren J, Palstra RJ, Klous P, Grosveld F, Galjart N, de Laat W: CTCF mediates long-range chromatin looping and local histone modification in the beta-globin locus. Genes Dev. 2006, 20: 2349-2354. 10.1101/gad.399506.PubMed CentralView ArticlePubMedGoogle Scholar
- Zuin J, Dixon JR, van der Reijden MI, Ye Z, Kolovos P, Brouwer RW, van de Corput MP, van de Werken HJ, Knoch TA, van Ijcken WF, Grosveld FG, Ren B, Wendt KS: Cohesin and CTCF differentially affect chromatin architecture and gene expression in human cells. Proc Natl Acad Sci U S A. 2014, 111: 996-1001. 10.1073/pnas.1317788111.PubMed CentralView ArticlePubMedGoogle Scholar
- Deng W, Lee J, Wang H, Miller J, Reik A, Gregory PD, Dean A, Blobel GA: Controlling long-range genomic interactions at a native locus by targeted tethering of a looping factor. Cell. 2012, 149: 1233-1244. 10.1016/j.cell.2012.03.051.PubMed CentralView ArticlePubMedGoogle Scholar
- Drissen R, Palstra RJ, Gillemans N, Splinter E, Grosveld F, Philipsen S, de Laat W: The active spatial organization of the beta-globin locus requires the transcription factor EKLF. Genes Dev. 2004, 18: 2485-2490. 10.1101/gad.317004.PubMed CentralView ArticlePubMedGoogle Scholar
- Hou C, Dale R, Dean A: Cell type specificity of chromatin organization mediated by CTCF and cohesin. Proc Natl Acad Sci U S A. 2010, 107: 3651-3656. 10.1073/pnas.0912087107.PubMed CentralView ArticlePubMedGoogle Scholar
- Lin YC, Benner C, Mansson R, Heinz S, Miyazaki K, Miyazaki M, Chandra V, Bossen C, Glass CK, Murre C: Global changes in the nuclear positioning of genes and intra- and interdomain genomic interactions that orchestrate B cell fate. Nat Immunol. 2012, 13: 1196-1204. 10.1038/ni.2432.PubMed CentralView ArticlePubMedGoogle Scholar
- Dekker J, Rippe K, Dekker M, Kleckner N: Capturing chromosome conformation. Science. 2002, 295: 1306-1311. 10.1126/science.1067799.View ArticlePubMedGoogle Scholar
- Tolhuis B, Palstra RJ, Splinter E, Grosveld F, de Laat W: Looping and interaction between hypersensitive sites in the active beta-globin locus. Mol Cell. 2002, 10: 1453-1465. 10.1016/S1097-2765(02)00781-5.View ArticlePubMedGoogle Scholar
- Hagege H, Klous P, Braem C, Splinter E, Dekker J, Cathala G, de Laat W, Forne T: Quantitative analysis of chromosome conformation capture assays (3C-qPCR). Nat Protoc. 2007, 2: 1722-1733. 10.1038/nprot.2007.243.View ArticlePubMedGoogle Scholar
- Naumova N, Smith EM, Zhan Y, Dekker J: Analysis of long-range chromatin interactions using Chromosome Conformation Capture. Methods. 2012, 58: 192-203. 10.1016/j.ymeth.2012.07.022.View ArticlePubMedGoogle Scholar
- Stadhouders R, Kolovos P, Brouwer R, Zuin J, van den Heuvel A, Kockx C, Palstra RJ, Wendt KS, Grosveld F, van Ijcken W, Soler E: Multiplexed chromosome conformation capture sequencing for rapid genome-scale high-resolution detection of long-range chromatin interactions. Nat Protoc. 2013, 8: 509-524. 10.1038/nprot.2013.018.View ArticlePubMedGoogle Scholar
- van de Werken HJ, Landan G, Holwerda SJ, Hoichman M, Klous P, Chachik R, Splinter E, Valdes-Quezada C, Oz Y, Bouwman BA, Verstegen MJ, de Wit E, Tanay A, de Laat W: Robust 4C-seq data analysis to screen for regulatory DNA interactions. Nat Methods. 2012, 9: 969-972. 10.1038/nmeth.2173.View ArticlePubMedGoogle Scholar
- Simonis M, Klous P, Splinter E, Moshkin Y, Willemsen R, de Wit E, van Steensel B, de Laat W: Nuclear organization of active and inactive chromatin domains uncovered by chromosome conformation capture-on-chip (4C). Nat Genet. 2006, 38: 1348-1354. 10.1038/ng1896.View ArticlePubMedGoogle Scholar
- Sexton T, Kurukuti S, Mitchell JA, Umlauf D, Nagano T, Fraser P: Sensitive detection of chromatin coassociations using enhanced chromosome conformation capture on chip. Nat Protoc. 2012, 7: 1335-1350. 10.1038/nprot.2012.071.View ArticlePubMedGoogle Scholar
- Gondor A, Rougier C, Ohlsson R: High-resolution circular chromosome conformation capture assay. Nat Protoc. 2008, 3: 303-313. 10.1038/nprot.2007.540.View ArticlePubMedGoogle Scholar
- Fullwood MJ, Liu MH, Pan YF, Liu J, Xu H, Mohamed YB, Orlov YL, Velkov S, Ho A, Mei PH, Chew EG, Huang PY, Welboren WJ, Han Y, Ooi HS, Ariyaratne PN, Vega VB, Luo Y, Tan PY, Choy PY, Wansa KD, Zhao B, Lim KS, Leow SC, Yow JS, Joseph R, Li H, Desai KV, Thomsen JS, Lee YK: An oestrogen-receptor-alpha-bound human chromatin interactome. Nature. 2009, 462: 58-64. 10.1038/nature08497.PubMed CentralView ArticlePubMedGoogle Scholar
- Dostie J, Dekker J: Mapping networks of physical interactions between genomic elements using 5C technology. Nat Protoc. 2007, 2: 988-1002. 10.1038/nprot.2007.116.View ArticlePubMedGoogle Scholar
- Lieberman-Aiden E, van Berkum NL, Williams L, Imakaev M, Ragoczy T, Telling A, Amit I, Lajoie BR, Sabo PJ, Dorschner MO, Sandstrom R, Bernstein B, Bender MA, Groudine M, Gnirke A, Stamatoyannapoulos J, Mirny LA, Lander ES, Dekker J: Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009, 326: 289-293. 10.1126/science.1181369.PubMed CentralView ArticlePubMedGoogle Scholar
- Dostie J, Richmond TA, Arnaout RA, Selzer RR, Lee WL, Honan TA, Rubio ED, Krumm A, Lamb J, Nusbaum C, Green RD, Dekker J: Chromosome Conformation Capture Carbon Copy (5C): a massively parallel solution for mapping interactions between genomic elements. Genome Res. 2006, 16: 1299-1309. 10.1101/gr.5571506.PubMed CentralView ArticlePubMedGoogle Scholar
- Jin F, Li Y, Dixon JR, Selvaraj S, Ye Z, Lee AY, Yen CA, Schmitt AD, Espinoza CA, Ren B: A high-resolution map of the three-dimensional chromatin interactome in human cells. Nature. 2013, 503: 290-294.PubMedGoogle Scholar
- Nativio R, Wendt KS, Ito Y, Huddleston JE, Uribe-Lewis S, Woodfine K, Krueger C, Reik W, Peters JM, Murrell A: Cohesin is required for higher-order chromatin conformation at the imprinted IGF2-H19 locus. PLoS Genet. 2009, 5: e1000739-10.1371/journal.pgen.1000739.PubMed CentralView ArticlePubMedGoogle Scholar
- van de Corput MP, de Boer E, Knoch TA, van Cappellen WA, Quintanilla A, Ferrand L, Grosveld FG: Super-resolution imaging reveals three-dimensional folding dynamics of the beta-globin locus upon gene activation. J Cell Sci. 2012, 125: 4630-4639. 10.1242/jcs.108522.View ArticlePubMedGoogle Scholar
- Palstra RJ, Tolhuis B, Splinter E, Nijmeijer R, Grosveld F, de Laat W: The beta-globin nuclear compartment in development and erythroid differentiation. Nat Genet. 2003, 35: 190-194. 10.1038/ng1244.View ArticlePubMedGoogle Scholar
- Song SH, Hou C, Dean A: A positive role for NLI/Ldb1 in long-range beta-globin locus control region function. Mol Cell. 2007, 28: 810-822. 10.1016/j.molcel.2007.09.025.PubMed CentralView ArticlePubMedGoogle Scholar
- Bernstein BE, Birney E, Dunham I, Green ED, Gunter C, Snyder M, ENCODE Project Consortium: An integrated encyclopedia of DNA elements in the human genome. Nature. 2012, 489: 57-74. 10.1038/nature11247.View ArticleGoogle Scholar
- Soler E, Andrieu-Soler C, de Boer E, Bryne JC, Thongjuea S, Stadhouders R, Palstra RJ, Stevens M, Kockx C, van Ijcken W, Hou J, Steinhoff C, Rijkers E, Lenhard B, Grosveld F: The genome-wide dynamics of the binding of Ldb1 complexes during erythroid differentiation. Genes Dev. 2010, 24: 277-289. 10.1101/gad.551810.PubMed CentralView ArticlePubMedGoogle Scholar
- Fraser P, Bickmore W: Nuclear organization of the genome and the potential for gene regulation. Nature. 2007, 447: 413-417. 10.1038/nature05916.View ArticlePubMedGoogle Scholar
- Misteli T: Beyond the sequence: cellular organization of genome function. Cell. 2007, 128: 787-800. 10.1016/j.cell.2007.01.028.View ArticlePubMedGoogle Scholar
- Simonis M, Klous P, Homminga I, Galjaard RJ, Rijkers EJ, Grosveld F, Meijerink JP, de Laat W: High-resolution identification of balanced and complex chromosomal rearrangements by 4C technology. Nat Methods. 2009, 6: 837-842. 10.1038/nmeth.1391.View ArticlePubMedGoogle Scholar
- Stadhouders R, van den Heuvel A, Kolovos P, Jorna R, Leslie K, Grosveld F, Soler E: Transcription regulation by distal enhancers: who’s in the loop?. Transcription. 2012, 3: 181-186. 10.4161/trns.20720.PubMed CentralView ArticlePubMedGoogle Scholar
- Li L, Lee JY, Gross J, Song SH, Dean A, Love PE: A requirement for Lim domain binding protein 1 in erythropoiesis. J Exp Med. 2010, 207: 2543-2550. 10.1084/jem.20100504.PubMed CentralView ArticlePubMedGoogle Scholar
- Li L, Jothi R, Cui K, Lee JY, Cohen T, Gorivodsky M, Tzchori I, Zhao Y, Hayes SM, Bresnick EH, Zhao K, Westphal H, Love PE: Nuclear adaptor Ldb1 regulates a transcriptional program essential for the maintenance of hematopoietic stem cells. Nat Immunol. 2011, 12: 129-136. 10.1038/ni.1978.PubMed CentralView ArticlePubMedGoogle Scholar
- Mylona A, Andrieu-Soler C, Thongjuea S, Martella A, Soler E, Jorna R, Hou J, Kockx C, van Ijcken W, Lenhard B, Grosveld F: Genome-wide analysis shows that Ldb1 controls essential hematopoietic genes/pathways in mouse early development and reveals novel players in hematopoiesis. Blood. 2013, 121: 2902-2913. 10.1182/blood-2012-11-467654.View ArticlePubMedGoogle Scholar
- Mukhopadhyay M, Teufel A, Yamashita T, Agulnick AD, Chen L, Downs KM, Schindler A, Grinberg A, Huang SP, Dorward D, Westphal H: Functional ablation of the mouse Ldb1 gene results in severe patterning defects during gastrulation. Development. 2003, 130: 495-505. 10.1242/dev.00225.View ArticlePubMedGoogle Scholar
- Bolger AM, Lohse M, Usadel B: Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014, [Epub ahead of print]Google Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.PubMed CentralView ArticlePubMedGoogle Scholar
- Quinlan AR, Hall IM: BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010, 26: 841-842. 10.1093/bioinformatics/btq033.PubMed CentralView ArticlePubMedGoogle Scholar
- Nagano T, Lubling Y, Stevens TJ, Schoenfelder S, Yaffe E, Dean W, Laue ED, Tanay A, Fraser P: Single-cell Hi-C reveals cell-to-cell variability in chromosome structure. Nature. 2013, 502: 59-64. 10.1038/nature12593.View ArticlePubMedGoogle Scholar
- van de Werken HJ, de Vree PJ, Splinter E, Holwerda SJ, Klous P, de Wit E, de Laat W: 4C technology: protocols and data analysis. Methods Enzymol. 2012, 513: 89-112.View ArticlePubMedGoogle Scholar
- R-Core-Team: R: A language and environment for statistical computing. 2013, Vienna, Austria: R Foundation for Statistical Computing, [http://www.r-project.org/]Google Scholar
- Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS: Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9: R137-10.1186/gb-2008-9-9-r137.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.