Work in the Heard laboratory was supported by a European Research Council Advanced Investigator award (XPRESS - AdG671027). A.L. is supported by a European Union Marie Sk&#322;odowska-Curie Actions Individual Fellowship (IF-838408). A.H. is supported by the ITN Innovative and Interdisciplinary Network ChromDesign, under the Marie Sk&#322;odowska-Curie Grant agreement 813327. The authors are thankful to Daniel Ibrahim (MPI for Molecular Genetics, Berlin) for helpful technical advice, to the NGS platform at Institut Curie (Paris), and to Vladimir Benes and the Genomics Core Facility at EMBL (Heidelberg), for support and assistance.

The mouse embryonic stem cells (mESCs) used in this study were derived from a cross of a TX/TX R26rtTA/rtTA female50 with a Mus musculus castaneus male according to the animal care guidelines of Institut Curie (Paris)51.
1. Probe design
<ol>
	<li>Design an array of biotinylated probes (120-mer RNA oligonucleotides) covering the target region of interest.
	<ol>
		<li>Tile the region of interest with overlapping oligonucleotides so that on average each sequence within the target is covered by two unique probes (2x coverage) (Figure 1).</li>
		<li>Exclude repetitive sequences from probe coverage to avoid enrichment of unspecific interactions. 
		&#8203;NOTE: To maximize the enrichment of informative ligation fragments, regions spanning 300 bp upstream and downstream of each DpnII restriction site across the target were defined (ChrX: 102,475,000-105,475,000), and 28,913 biotinylated probes were designed according to the SureSelect DNA target enrichment technology through the Sure Design platform52. According to this strategy, up to a maximum of 40 bases of repetitive sequences are allowed in each oligonucleotide to minimize the enrichment of unspecific interactions. The probe array was synthesized by Agilent. Here, DpnII is used as a restriction enzyme for two reasons: (1) it is a four-cutter routinely used in several 3C-based methods53; and (2) it maximizes the chances of capturing informative single nucleotide polymorphisms (SNPs) in the proximity of the cutting sites in comparison with other restriction enzymes that were tested in silico in F1 hybrid lines used in this study (C57BL/6J x CASTEi/J).</li>
	</ol>
	</li>
</ol>
2. Experimental procedure
<ol>
	<li>Cell preparation
	<ol>
		<li>Seed the appropriate number of cells on one or multiple cell culture plates to achieve a total cell number of &#8805; 5 x 107 cells on the day of fixation. 
		NOTE: Mouse embryonic stem cells (mESCs) were used in this study. mESCs are plated on gelatinized (0.1% gelatine in 1x PBS - o/n at 37 &#176;C, 5% CO2 incubator) cell culture plates in mESCs medium containing 2i + LIF and batch-tested fetal calf serum (DMEM, 15% FBS, 0.1 mM &#946;-mercaptoethanol, 1,000 U/mL&#8722;1 leukaemia inhibitory factor (LIF), CHIR99021 (3 &#181;M), and PD0325901 (1 &#181;M)). For this cell type, one 80% confluent 10 cm plate contains approximately 2 x 107 cells.</li>
		<li>Prepare one additional cell culture plate for cell counting. 
		NOTE: A smaller cell culture plate can be used to reduce media usage. In this case, the number of cells to seed on the smaller plate needs to be adjusted accordingly (e.g., 3x fewer cells on a 10 cm plate compared to a 15 cm plate).</li>
	</ol>
	</li>
	<li>Formaldehyde fixation
	<ol>
		<li>Estimate the total number of cells to be cross-linked.
		<ol>
			<li>Before starting the cross-linking reaction, trypsinize and count cells from the control plate prepared specifically for cell counting using an automated cell counter according to the manufacturer&#39;s instructions.</li>
			<li>Include a viability staining (e.g., Trypan Blue) to determine the percentage of viable cells54. From this cell count, estimate the total number of cells in the plate(s) prepared for cross-linking.</li>
		</ol>
		</li>
		<li>Remove the culture medium from the plates prepared for cross-linking and replace it with the appropriate amount of fixation solution (2% formaldehyde in cell culture medium). Use 10 mL on a 10 cm plate (e.g., ~20 mL for a 15 cm plate). 
		NOTE: Add an exact volume of fixation solution. If fixing adherent cells is not possible, this step can be adapted to trypsinized cells and performed in 30 mL of fixation solution in 50 mL conical centrifuge tubes. Formaldehyde must not be more than 1 year old. It is preferred to use single-use vials. The fixation solution must be brought to room temperature (RT) before use. 
		CAUTION: Formaldehyde is hazardous and needs to be handled according to appropriate health and safety regulations.</li>
		<li>Fix for 10 min at RT under gentle mixing on a shaker.</li>
		<li>Quench the fixation reaction by addition of 2.5 M glycine-1x PBS to a final concentration of 0.125 M. Add 530 &#181;L of 2.5 M glycine-1x PBS to 10 mL on a 10 cm plate (e.g., 1060 &#181;L to 20 mL on a 15 cm plate). 
		NOTE: If the cells were fixed in solution, quench the fixation reaction with 1590 &#181;L of 2.5 M glycine-1x PBS.</li>
		<li>Incubate for 5 min at RT, gently mixing on a shaker.</li>
		<li>Transfer the plates to ice and incubate for an additional 15 min on ice while gently mixing on a shaker. 
		NOTE: From now on, the cells must be kept on ice, and buffers must be pre-chilled to avoid further cross-linking. Move to a cold room if many plates need to be processed.</li>
		<li>Remove the fixation solution from the cells by pouring it into a beaker to ensure fast handling. 
		NOTE: Make sure to dispose of the formaldehyde-containing liquid waste according to the appropriate health and safety regulations.</li>
		<li>Rinse the 10 cm plate quickly two times with 5 mL of cold 0.125 M glycine-1x PBS (8 mL for a 15 cm plate) to wash off the debris and dead cells. Remove the liquid from the plate by pouring it into a beaker to ensure fast handling.</li>
		<li>Add 5 mL of cold 0.125 M glycine-1x PBS to the 10 cm plate (10 mL for a 15 cm plate) and quickly scrape the cells from the plate using a plastic cell scraper.</li>
		<li>Transfer the cell suspension into a pre-cooled 50 mL conical centrifuge tube using a serological pipette.</li>
		<li>Rinse the plate twice with 5 mL of cold 0.125 M glycine-1x PBS and add the cell suspension to the conical centrifuge tube.</li>
		<li>Spin down at 480 x g for 10 min at 4 &#176;C. 
		NOTE: If the cells were fixed in solution, transfer the cell to a pre-chilled conical centrifuge tube and spin down at 480 x g for 10 min at 4 &#176;C. Remove the fixation solution by pouring it into a beaker and wash three times in 10 mL of cold 0.125 M glycine-1x PBS. Make sure to resuspend the cells at each washing step.</li>
		<li>Remove the supernatant by aspirating with a benchtop aspiration system. Resuspend the cells in 500 &#181;L of 1x PBS per 1 x 107 cells by carefully pipetting up and down with a P1000 pipette. To resuspend cells in the accurate volume, refer to the total cell number estimation obtained in 2.2.1.</li>
		<li>Aliquot 500 &#181;L of the cell suspension into the calculated number of 1.5 mL microcentrifuge tubes (1 x 107 cells/tube).</li>
		<li>Spin down at 480 x g for 10 min at 4 &#176;C.</li>
		<li>Remove the supernatant with a benchtop aspiration system and snap freeze the cell pellets in liquid nitrogen. Store the dry cell pellets at -80 &#176;C. 
		NOTE: Samples can be stored for at least 1 year.</li>
	</ol>
	</li>
	<li>Cell lysis
	<ol>
		<li>Thaw the frozen pellet(s) on ice.</li>
		<li>Prepare 1.5 mL of lysis buffer in H20 per sample: Add 10 mM Tris-HCl, pH 8.0, 10 mM NaCl, and 0.2% NP40.</li>
		<li>Add 600 &#181;L of the cold lysis buffer and resuspend well on ice.</li>
		<li>Incubate on ice for 15 min to let the cells swell.</li>
		<li>Spin down at 2655 x g for 5 min at 4 &#176;C and remove the supernatant using a benchtop aspiration system.</li>
		<li>To remove debris, resuspend the pellet in 1 mL of the cold lysis buffer, spin down at 2655 x g for 5 min at 4 &#176;C, and remove the supernatant.</li>
		<li>Spin again shortly at 2655 x g and 4 &#176;C and remove as much of the remaining supernatant as possible using a benchtop aspiration system fitted with a P200 tip.</li>
		<li>Resuspend in 100 &#181;L of 0.5% (vol/vol) SDS.</li>
		<li>Incubate in a thermomixer at 62 &#176;C, swirling at 1400 rpm for 10 min.</li>
		<li>Add 290 &#181;L of H2O + 50 &#181;L of 10% TritonX-100 and mix well, avoiding air bubbles.</li>
		<li>Incubate in a thermomixer at 37 &#176;C, swirling at 1400 rpm for 15 min.</li>
		<li>Add 50 &#181;L of 10x Dpnll buffer and invert the tube to mix.</li>
		<li>Take 50 &#181;L of undigested DNA for quality control into a separate tube. Do not forget to take the undigested control sample.</li>
	</ol>
	</li>
	<li>DpnII digestion
	<ol>
		<li>Add 10 &#181;L of Dpnll high concentration (500 U total) and mix by inverting.</li>
		<li>Incubate the samples and the undigested control in a thermomixer at 37 &#176;C, swirling at 1400 rpm for &#62;4 h.</li>
		<li>Add 10 &#181;L of Dpnll high concentration (500 U total) at the end of the day.</li>
		<li>Incubate the samples and the undigested control at 37 &#176;C, swirling at 1400 rpm overnight.</li>
		<li>Add 10 &#181;L of Dpnll high concentration (500 U total) at the beginning of the next day to the samples.</li>
		<li>Incubate the samples and the undigested control in a thermomixer at 37 &#176;C, swirling at 1400 rpm for 4 h.</li>
	</ol>
	</li>
	<li>Ligation and reversal of cross-linking
	<ol>
		<li>Incubate the tubes at 65 &#176;C for 20 min at 1400 rpm. 
		NOTE: Do not add SDS at this point. The idea is to preserve nuclear integrity, so the ligation is carried out inside the nuclei, circumventing the need for extreme dilution.</li>
		<li>Cool down the samples on ice for a maximum of 5-10 min. To avoid SDS precipitation, do not leave the samples on ice longer than this.</li>
		<li>Take 50 &#181;L of the unligated digested DNA for quality control in a separate tube. Store the undigested and unligated controls at -20 &#176;C. 
		NOTE: Do not forget to take the unligated control sample.</li>
		<li>Add 800 &#181;L of ligation cocktail: 122 &#181;L of 10x ligase buffer, 8 &#181;L of T4 ligase (30 U/&#181;L), and 670 &#181;L of H20.</li>
		<li>Incubate at 16 &#176;C, swirling at 1000 rpm overnight.</li>
		<li>Add 7.5 &#181;L of Proteinase K (20 mg/mL) to the samples and 2 &#181;L to the controls.</li>
		<li>Incubate at 65 &#176;C for 4 h at 1000 rpm.</li>
	</ol>
	</li>
	<li>DNA purification
	<ol>
		<li>Transfer the samples on ice to pre-cooled 15 mL conical centrifuge tubes and add 2 mL of water, 10.5 mL of ice-cold EtOH, and 583 &#181;L of 3 M NaAC. 
		NOTE: Additional water aims to prevent carry-over of DTT into the pellet.</li>
		<li>Add 200 &#181;L of ice cold EtOH, 10.8 &#181;L of NaAC, and 1 &#181;L of the coprecipitant to the undigested and unligated quality controls.</li>
		<li>Incubate at -80 &#176;C for at least 4 h up to overnight.</li>
		<li>Spin the 15 mL tubes at 2200 x g at 4 &#176;C for 45 min.</li>
		<li>Spin the 1.5 mL control tubes at 20,500 x g at 4 &#176;C for 30 min.</li>
		<li>Wash once with 3 mL (samples) and 1 mL (controls) of ice-cold 70% EtOH.</li>
		<li>Spin at 2200 x g (samples) or 20,500 x g (controls) at 4 &#176;C for 10 min.</li>
		<li>Carefully remove EtOH and air dry at RT for 10-15 min; do not overdry.</li>
		<li>Resuspend the samples and controls in 100 &#181;L and 20 &#181;L of H20, respectively.</li>
		<li>Add 1 &#181;L of RNAseA and incubate at 37 &#176;C, swirling at 1400 rpm for 30 min.</li>
	</ol>
	</li>
	<li>Quality control of 3C template preparation
	<ol>
		<li>Quantify each sample and control using a fluorometer kit for high-sensitivity DNA concentration measurements.</li>
		<li>Load 100-200 ng of each sample and each control on a 1% agarose/1x TBE gel.</li>
		<li>Verify that the gel image shows the expected result by comparing the differences in DNA fragment sizes of the controls and the 3C template as shown in Figure 2A.</li>
		<li>Store the samples and controls at -20 &#176;C.</li>
	</ol>
	</li>
	<li>Hybridization, capture, and sample processing for multiplexed sequencing
	<ol>
		<li>To hybridize the array of biotinylated RNA probes to the 3C template, capture the targeted ligation fragments, and prepare the samples for multiplexed sequencing according to the target enrichment system used in this study for paired-end multiplexed sequencing (see Table of Materials). Follow the protocol according to the manufacturer&#39;s instructions while introducing the following minor modifications:
		<ol>
			<li>Section 2 of the manufacturer&#39;s protocol: Sample preparation
			<ol>
				<li>Follow the instructions for target enrichment starting from 3 &#181;g of gDNA input.</li>
				<li>Shear the DNA in a sonicator using the following specifications: 10% duty cycle, 4 intensity, 200 cyc/burst, and 130 s. Start with 4 &#181;g of 3C template resuspended in 130 &#181;L of water for each capture reaction to ensure enough material to continue the sample preparation with 3 &#181;g of the sheared DNA.</li>
				<li>Assess the quality of sheared DNA. Run 1 &#181;L of the sheared DNA on a DNA bioanalyzer according to the high-sensitivity protocol. Expect a distribution of fragment size between 150-700 bp (Figure 2).</li>
				<li>Purify the sample using solid phase reversible immobilization (SPRI) beads. Add 124 &#181;L of SPRI beads to 124 &#181;L of the DNA sample to perform a 1:1 left side-size selection according to the manufacturer&#39;s instructions and elute in 25 &#181;L of nuclease-free water. This purification step will remove shorter fragments to enrich fragments of around 300 bp (Figure 2). 
				NOTE: The amount of samples and SPRI beads used at this step take into account the volume loss that occurred while transferring the samples to new tubes and running the quality controls at the Bioanalyzer. All subsequent size selection steps are performed according to the ratios recommended by the manufacturer&#39;s protocol. DNA elution from SPRI beads is performed at RT throughout the protocol.</li>
				<li>Assess the quality of size-selected sheared DNA. Run 1 &#181;L of the sheared DNA on the DNA bioanalyzer according to the high-sensitivity (HS) protocol. Expect a distribution of fragment sizes with the highest enrichment at 300 bp (Figure 2). Go ahead with the quantification of the sheared DNA if shearing was successful.</li>
				<li>Quantify the sheared DNA with a fluorometer kit for HS DNA concentration measurements. 
				NOTE: If DNA shearing results in a DNA yield of &#60;3 &#181;g, perform a second round of DNA shearing with another 4 &#181;g of DNA and combine the sheared DNA samples after the first SPRI bead purification step to achieve a total of 3 &#181;g of sheared DNA.</li>
				<li>Add nuclease-free water to the size-selected cleaned DNA sample (3 &#181;g total) to a final volume of 48 &#181;L and proceed with the end repair reaction according to the manufacturer&#39;s protocol.</li>
				<li>After ligation of paired-end adaptors, amplify the library by performing five cycles of pre-capture PCR according to the manufacturer&#39;s instructions (the conditions for PCR and the primers are provided in the kit).</li>
			</ol>
			</li>
			<li>Section 4 of the manufacturer&#39;s protocol: Hybridization and capture
			<ol>
				<li>To hybridize the prepared DNA samples to the target-specific RNA probes, dilute 750 ng of DNA samples in a final volume of 3.4 &#181;L, resulting in an initial concentration of 221 ng/&#181;L. For DNA samples diluted in bigger volumes, use a speed-vacuum concentrator to reduce to the final volume. A speed-vacuum concentration (250 x g; &#8804;45 &#176;C) for 15-20 min is normally sufficient for the samples resuspended in 10 &#181;L. Make sure to have the same input volume for each sample before starting the speed-vacuum concentrator.</li>
				<li>Incubate the hybridization mixture for 16-18 h at 65 &#176;C with a heated lid at 105 &#176;C according to the manufacturer&#39;s instructions.</li>
			</ol>
			</li>
			<li>Section 5 of the manufacturer&#39;s protocol: Indexing and sample processing for multiplexed sequencing
			<ol>
				<li>To amplify the captured libraries with indexing primers, perform 12 cycles of post-capture PCR according to the manufacturer&#39;s instructions (the conditions for PCR and the primers are provided in the kit).</li>
			</ol>
			</li>
		</ol>
		</li>
	</ol>
	</li>
	<li>Next-generation sequencing
	<ol>
		<li>To run multiple capture Hi-C libraries on the same flow cell, prepare an equimolar mixture of the capture libraries and sequence 100-120 M reads per library.</li>
		<li>If the allele-specific analysis is needed, sequence 150 bp paired-end to ensure sufficient SNP coverage.</li>
	</ol>
	</li>
</ol>
3. Data analysis
<ol>
	<li>Apply the HiC-Pro pipeline to perform Capture Hi-C data analysis55. HiC-Pro provides quality controls at each step of the processing, including (Figure 3): 
	(i)The alignment rate on the reference genome specifying the fraction of reads spanning a ligation site, as well as the number of pairs and singletons. 
	(ii)The fraction of valid ligation products and non-informative read pairs (dangling-end, self-ligation, etc). 
	(iii)The fraction of short/long-range and intra/inter-chromosomal contacts. 
	(iv)The fraction of on-target contacts for Capture Hi-C. 
	(v)The fraction of allele-specific reads if specified. 
	NOTE: HiC-Pro supports a wide range of protocols including in situ Hi-C and Capture Hi-C. In the latter case, the user simply needs to specify the targeted region (BED format) in the configuration file. Once the data are processed, the HiC-Pro outputs can easily be converted into a cooler object for downstream analysis56. At this step, the contact maps at various resolutions are normalized using the ICE method previously described by Imakaev and colleagues57. Several analyses can then be run to call chromosome compartments, TADs, or chromatin loops (for review58). The workflow of the protocol is shown in Figure 4. Here, the &#39;cooltools&#39; suite is applied to calculate the insulation score and TADs boundaries, as illustrated in Figure 5 and Figure 659.</li>
</ol>

Kai Hauschulz is Field Application Scientist at Agilent Technologies - Diagnostic and Genomics Group. All other authors declare no competing interests.

Here we describe a relatively quick and easy Capture Hi-C protocol to characterize the higher-order organization of megabase-sized genomic regions at 5-10 kb resolution. Capture Hi-C belongs to the family of Capture-C technologies that are designed to enrich targeted chromatin interactions from genome-wide 3C or Hi-C templates. To date, the large majority of Capture-C applications have been exploited to map chromatin contacts of relatively small regulatory elements scattered across the entire genome. In the first Capture-C protocol, multiple overlapping RNA biotinylated probes were used to capture &#62;400 pre-selected promoters in 3C libraries prepared from erythroid cells31. The same strategy was subsequently improved in Next Generation (NG) and Nuclear Titrated (NuTi) Capture-C to achieve high-resolution interaction profiles of &#62;8,000 promoters by using single 120 bp DNA baits spanning single restriction sites and two sequential rounds of Capture to maximize the enrichment of informative ligation fragments32,40. These strategies led to the functional dissection of cis-acting elements in many different contexts, including mouse embryonic development, cell differentiation, X-chromosome inactivation, and gene mis-regulation in pathological conditions46,63,65,66,67,68,69,70,71.
In Promoter Capture Hi-C (PCHi-C), &#62;22,000 annotated promoters containing restriction fragments were pulled down from Hi-C libraries by hybridization of single RNA 120-mer biotinylated probes at either or both ends of the restriction fragment34,72. This method allowed dissection of the interactome of thousands of promoters in a rapidly increasing number of cell types, including mouse embryonic stem cells, fetal liver cells, and adipocytes34,35,72,73, but also human lymphoblastoid lines, hematopoietic progenitors, epidermal keratinocytes, and pluripotent cells37,74,75,76,77.
In comparison with these target enrichment technologies, Capture Hi-C targets contiguous genomic regions up to the megabase scale, thereby spanning one or more TADs and encompassing regulatory landscapes of genes. The entire region of interest must be tiled with an array of biotinylated probes encompassing each DpnII restriction site within the target. The hybridization of the biotinylated array to the 3C template, its subsequent streptavidin-based capture, and processing for multiplexed sequencing is performed using a target enrichment system for Illumina Paired-End multiplexed sequencing. The entire protocol is fast, as it can be performed in 1 week from 3C library preparation through to NGS sequencing, and it requires only minor adaptations and/or custom-specific troubleshooting.
The protocol also provides advantages in comparison with other 3C-based methods. To obtain interaction maps at a resolution of 5-10 kb, we sequenced 100-120 M paired-end reads. As a comparison, we used here a Hi-C dataset of 571 M reads to reach a 20 kb resolution64 (GSM2053973), and at least 1 billion reads would be required to reach a 5 kb resolution with chromosome-wide Hi-C22.
Capture Hi-C as used in the present study reaches a much higher resolution than the previously published 5C based on a 6-bp cutter restriction enzyme47 (Supplementary Table 1). Importantly, the strategy designed to enrich and amplify targeted interactions in 5C does not allow for allele-specific analysis of chromatin interactions. On the contrary, Capture Hi-C data can be mapped allele-specifically, allowing the dissection of the 3D structural landscapes of pairs of homologous chromosomes, for example in human cells or in F1 hybrid cell lines derived by crossing genetically different mouse strains78. To generate allele-specific Capture Hi-C interaction maps at 5 kb resolution, we sequenced 150 bp paired-end reads to increase SNP coverage. Similar allele-specific approaches can be applied to human cell lines, for which the annotation of SNPs is available22.
Importantly, although Capture Hi-C generally ensures high resolution while improving the affordability of sequencing costs, the production of custom-tailored biotinylated oligonucleotides does have an impact on the overall cost of this method. Therefore, the choice of the most suitable 3C method will differ for different applications, and will depend on the biological question that is being addressed and the resolution required, as well as the size of the region of interest. Other Capture Hi-C protocols developed share key features with the protocol described here. For example, a Capture Hi-C strategy was applied to characterize ~50 kb to 1 Mb genomic regions spanning noncoding variants associated with breast and colorectal cancer risk; in this protocol, target regions were pulled down from Hi-C libraries by hybridizing 120-mer RNA baits tiling the target regions at a 3x coverage33,38,79. Similarly, HYbrid Capture Hi-C (Hi-C2) was used to target interactions within regions of interest up to 2 Mb80. In both protocols, the use of a Hi-C template enriched for biotin pulled-down ligation fragments increased the percentage of total informative reads compared to our protocol. For example, in the Hi-C dataset we used here for comparison64 (GSM2053973), the percentage of valid pairs after the removal of duplicates is 4.8-fold higher than the valid pairs obtained in Capture Hi-C as described in Figure 3 and Supplementary Table 1. However, the consecutive pull-down of biotinylated ligated fragments and hybridized probes makes the protocol significantly more complex and time consuming while possibly decreasing the complexity of the captured region.
Another available method to enrich 3C templates with tiling probes is Tiled-C, which was applied to study chromatin architecture at high spatial and temporal resolution during mouse erythroid differentiation43. In Tiled-C, a panel of 70 bp biotinylated probes is used to enrich contacts within large-scale regions in two consecutive rounds of capture to generate very high-resolution maps of targeted interactions43,81. The double capture enrichment also makes the protocol longer and more complex when compared to Capture Hi-C. However, different to the Capture-C strategies targeting single restriction sites, in Tiled-C the second round of capture does not seem to significantly increase the capture efficiency, and therefore can probably be omitted43. Finally, a similar tiling approach based on the same target enrichment strategy used in this study was applied to the dissection of regulatory landscapes encompassing structural variants described in patients with congenital malformations and re-engineered in transgenic mice41,42. In this case, the tiling array of probes was designed across the entire target rather than in the proximity of DpnII restriction sites41. Nevertheless, this work was seminal in highlighting the sensitivity and power of this strategy to achieve high-resolution characterization of large genomic regions in different contexts41,42,48.
In conclusion, the protocol described here represents an easy, robust, and powerful strategy for the high-resolution 3D characterization of any genomic regions of interest. The application of this approach to different model systems, cell types, developmentally-regulated chromatin landscapes, and gene regulation in healthy and pathological conditions is likely to facilitate our understanding of the interplay and causality between genome topology and gene regulation, one of the fundamental open questions in the epigenetics field. Furthermore, applying Capture Hi-C to map long-range interactions and higher-order chromatin folding of risk variants identified by GWAS studies has the potential to reveal the functional relevance of noncoding genomic loci associated with human diseases in different contexts, thereby providing novel insights into the processes potentially underlying pathogenesis.

The linear genome holds all the information necessary for an organism to undergo embryonic development and survive throughout adulthood. However, instructing genetically identical cells to perform different functions is fundamental for accurately controlling which information is used in specific contexts, including different tissues and/or developmental stages. The three-dimensional organization of the genome is thought to participate in this accurate spatio-temporal regulation of gene activity by facilitating or preventing the physical interaction between regulatory elements that can be separated by several hundred kilobases in the linear genome (for reviews1,2,3). In the last 20 years, our understanding of the interplay between genome folding and activity has rapidly increased, largely owing to the development of chromosome conformation capture technologies (3C) (for review4,5,6,7). These methods measure the frequency of interactions between any regions of the genome and rely on the ligation of DNA sequences that are in close 3D proximity within the nucleus. The most common 3C protocols start with the fixation of cell populations with a cross-linking agent such as formaldehyde. The cross-linked chromatin is then digested with a restriction enzyme, although MNase digestion has also been used8,9. After digestion, free DNA ends in close spatial proximity are re-ligated, and cross-linking is reversed. This step gives rise to the 3C &#39;library&#39; or &#39;template&#39;, a mixed pool of hybrid fragments in which sequences that were in 3D proximity to the nucleus have higher chances of getting ligated in the same DNA fragment. The downstream quantification of these hybrid fragments enables inferring the 3D conformation of genomic regions that are located thousands of base pairs apart in the linear genome but might interact in the 3D space.
Many different approaches have been developed to characterize the 3C library, differing both in terms of which subsets of ligation fragments are analyzed and which technology is used for their downstream quantification. The original 3C protocol relied on the selection of two regions of interest and the quantification of their &#39;one versus one&#39; interaction frequency by PCR10,11. The 4C approach (circular chromosome conformation capture) measures the interactions between a single locus of interest (i.e., the &#39;view-point&#39;) and the rest of the genome (&#39;one versus all&#39;)12,13,14. In 4C, the 3C library undergoes a second round of digestion and re-ligation to generate small circular DNA molecules that are PCR amplified by view-point specific primers15. 5C (chromosome conformation capture carbon copy) enables the characterization of 3D interactions across larger regions of interest, providing insights into higher-order chromatin folding within that region (&#39;many versus many&#39;)16. In 5C, the 3C library is hybridized to a pool of oligonucleotides overlapping restriction sites that can be subsequently amplified by multiplex PCR with universal primers15. In both 4C and 5C, the informative DNA fragments were initially quantified by microarrays and later by next-generation sequencing (NGS)17,18,19. These strategies characterize targeted regions of interest but cannot be applied to map genome-wide interactions. This latter goal is achieved with Hi-C, a 3C-based high-throughput strategy in which massively parallel sequencing of the 3C template allows the unbiased characterization of chromatin folding at the genome-wide level (&#39;all versus all&#39;)20. The Hi-C protocol includes the incorporation of a biotinylated residue at the digested fragments&#39; ends, which is followed by pull-down of ligation fragments with streptavidin beads to increase the recovery of ligated fragments20.
Hi-C revealed that mammalian genomes are structurally organized at multiple scales in the 3D nucleus. At the megabase scale, the genome is divided into regions of active and inactive chromatin, the A and B compartments, respectively20,21. The existence of further subcompartments represented by different chromatin and activity states was also subsequently shown22. At higher resolution, the genome is further partitioned into sub-megabase self-interacting domains called topologically associating domains (TADs), first revealed by Hi-C and 5C analysis of the human and mouse genomes23,24. Unlike compartments which vary in a tissue-specific manner, TADs tend to be constant (although there are many exceptions). Importantly, TAD boundaries are conserved across species25. In mammalian cells, TADs frequently encompass genes sharing the same regulatory landscape and have been shown to represent a structural framework that facilitates gene co-regulation while limiting the interactions with neighboring regulatory domains (for review3,26,27,28). Furthermore, within TADs, interactions due to CTCF sites at the base of cohesin-extruded loops may increase the probability of promoter-enhancer or enhancer-enhancer interactions (for review29).
In Hi-C, compartments and TADs can be detected at 1 Mb to 40 kb resolution, but higher resolution can be achieved to characterize smaller scale contacts such as looping interactions between distal elements at the scale of 5-10 kb. However, increasing the resolution to be able to detect such loops efficiently by HiC requires a significant increase in sequencing depth and, therefore, sequencing costs. This is exacerbated if the analysis needs to be allele-specific. Indeed, an X-fold increase in resolution requires an X2 increase in sequencing depth, meaning that high-resolution and allele-specific genome-wide approaches can be prohibitively expensive30.
To improve cost-effectiveness and affordability while maintaining high-resolution, target regions of interest can be physically pulled down from genome-wide 3C or Hi-C libraries following their hybridization with complementary biotin-labeled oligonucleotide probes before downstream sequencing. These target enrichment strategies are referred to as Capture-C methods and allow the interrogation of interactions of hundreds of target loci scattered across the genome (i.e., Promoter Capture (PC) Hi-C; Next Generation (NG) Capture-C; Low Input (LI) Capture-C; Nuclear Titrated (NuTi) Capture-C; Tri-C)31,32,33,34,35,36,37,38,39,40, or across regions spanning up to several megabases (i.e., Capture HiC; HYbrid Capture Hi-C (Hi-C2); Tiled-C)41,42,43. Two aspects can vary in capture-based methods: (1) the nature and design of biotinylated oligonucleotides (i.e., RNA or DNA, single oligos capturing dispersed genomic targets or multiple oligos tiling a region of interest); and (2) the template that is used for pulling down targets which can be the 3C or Hi-C library, the latter consisting of biotinylated restriction fragments pulled down from the 3C library.
Here, a Capture Hi-C protocol based on the enrichment of target contacts from the 3C library is described. The protocol relies on the design of a custom-tailored tiling array of biotinylated RNA probes and can be performed in 1 week from the 3C library preparation to the NGS sequencing. The protocol is fast, simple, and allows characterizing the higher order 3D organization of megabase-sized regions of interest at 5 kb resolution while improving time effectiveness and affordability in comparison with other 3C methods. The Capture Hi-C protocol was applied to the master regulatory locus of X-chromosome inactivation (XCI), the X-inactivation center (Xic), which hosts the Xist noncoding RNA. The Xic has previously been the subject of extensive structural and functional analyses (for review44,45). In mammals, XCI compensates for the dosage of X-linked genes between females (XX) and males (XY) and involves the transcriptional silencing of almost the entirety of one of the two X chromosomes in female cells. The Xic has represented a powerful, gold standard locus for studies in 3D genome topology and the interplay with gene regulation44. 5C analysis of the Xic in mouse embryonic stem cells (mESCs) led to the discovery and naming of TADs, providing the first insights into the functional relevance of topological partitioning and gene co-regulation24. The topological organization of the Xic was subsequently shown to be critically involved in the appropriate developmental timing of Xist upregulation and XCI46, and unsuspected cis-regulatory elements that can influence gene activity within and between TADs were also recently discovered within the Xic47,48,49. Applying Capture Hi-C to 3 Mb of the mouse X chromosome spanning the Xic demonstrates the power of this approach in dissecting large-scale chromatin folding at high resolution. A detailed and easy-to-follow protocol is provided, starting from the design of the array of biotinylated probes across every DpnII restriction site within the region of interest to the generation of the genome-wide 3C library, the hybridization and capture of target contacts, and downstream data analysis. An overview of the appropriate quality controls and expected results is also included, and both the strengths and limitations of the approach are discussed in light of similar existing methods.

<table><tbody><tr><td>10x PBS pH 7.4</td><td>Gibco</td><td>10010-023</td><td></td></tr><tr><td>37% (vol/vol) paraformaldehyde solution</td><td>Electron Microscopy Sciences</td><td>15686</td><td>single use glass-vials; do not reuse</td></tr><tr><td>50 mL PP conical tube</td><td>Falcon</td><td>352070</td><td></td></tr><tr><td>Agarose</td><td>Sigma</td><td>A9539-500g</td><td></td></tr><tr><td>Bioanalyzer</td><td>Agilent</td><td>G2939BA</td><td></td></tr><tr><td>Cell Scrapers - 25 cm Handle and 3.0 cm Blade</td><td>Falcon</td><td>353089</td><td></td></tr><tr><td>CHIR99021</td><td>Axon Medchem BV</td><td>Axon 1386</td><td></td></tr><tr><td>cOmplete Mini, Protease inhibitor cocktail (EDTA-free)</td><td>Merck</td><td>11836170001</td><td></td></tr><tr><td>Countess Cell Counting Chamber Slides</td><td>Invitrogen</td><td>C10228</td><td></td></tr><tr><td>Countess II FL</td><td>Invitrogen</td><td>ZGEXSCCOUNTESS2FL</td><td>Automated cell counter</td></tr><tr><td>Covaris S2</td><td>Covaris</td><td>500217</td><td>Sonicator</td></tr><tr><td>DNA LoBind tube, 1.5 mL</td><td>Eppendorf</td><td>30108051</td><td></td></tr><tr><td>DpnII (50000 units/mL)</td><td>New England Biolabs</td><td>R0543M</td><td></td></tr><tr><td>Dulbecco&#039;s Modified Eagle Medium (DMEM)</td><td>Merck</td><td>D6429</td><td></td></tr><tr><td>Ethanol (100%)</td><td>Merck</td><td>1.00983.2500</td><td></td></tr><tr><td>Fetal Bovine Serum (FBS)</td><td>Thermo Scientific</td><td>10270106</td><td></td></tr><tr><td>gelatine from porcine skin</td><td>Sigma</td><td>G1890</td><td></td></tr><tr><td>GeneRuler 1 kb Plus DNA Ladder</td><td>Thermo Scientific</td><td>SM0313</td><td></td></tr><tr><td>GlycoBlue</td><td>Thermo Scientific</td><td>AM9516</td><td>Coprecipitant</td></tr><tr><td>High-Sensitivity Bioanlayzer chips</td><td>Agilent</td><td>5067-4626</td><td></td></tr><tr><td>Large Cooling Centrifuge 5920 R</td><td>Eppendorf</td><td>5948000018</td><td></td></tr><tr><td>leukaemia inhibitory factor (LIF)</td><td>Merck</td><td>ESG1107</td><td></td></tr><tr><td>Liquiport</td><td>KNF</td><td>NF300</td><td>Benchtop aspiration system</td></tr><tr><td>Low-binding filter tips</td><td>Biozym</td><td>VT0260U, VT0240, VT0220, VT0200U</td><td></td></tr><tr><td>Molecular biology grade water</td><td>Merck</td><td>W3500-6x500ML</td><td></td></tr><tr><td>Next Seq 500</td><td>Illumina</td><td>SY-415-1001</td><td></td></tr><tr><td>Next Seq 500 High Output v2 Kit (300 cycles)</td><td>Illumina</td><td>FC-404-2004</td><td></td></tr><tr><td>Nonidet P40 Substitute (NP40)</td><td>Merck</td><td>11332473001</td><td></td></tr><tr><td>PD0325901</td><td>Axon Medchem BV</td><td>Axon 1408</td><td></td></tr><tr><td>Protease inhibitor cocktail (EDTA-free)</td><td>Merck</td><td>11873580001</td><td></td></tr><tr><td>Proteinase K - recombinant, PCR-grade (20 mg/mL)</td><td>Thermo Scientific</td><td>EO0491</td><td></td></tr><tr><td>Qubit 2.0</td><td>Thermo Scientific</td><td>Q32871</td><td></td></tr><tr><td>Qubit assay tubes</td><td>Thermo Scientific</td><td>Q32856</td><td></td></tr><tr><td>Qubit dsDNA High Sensitivity kit</td><td>Thermo Scientific</td><td>Q32851</td><td></td></tr><tr><td>RNase A (10 mg/mL)</td><td>Thermo Scientific</td><td>EN0531</td><td></td></tr><tr><td>Sodium acetate pH 5.2 (3M)</td><td>Merck</td><td>S7899</td><td></td></tr><tr><td>speed vacuum concentrator</td><td>Eppendorf</td><td>EP5305000100-1EA</td><td></td></tr><tr><td>Agencourt AMPureXP</td><td>Beckman Coulter</td><td>A63881</td><td>SPRI beads</td></tr><tr><td>SureSelect Target Enrichment Box 1</td><td>Agilent</td><td>5190-8645</td><td></td></tr><tr><td>SureSelect Target Enrichment Kit ILM Indexing Hyb Module Box 2</td><td>Agilent</td><td>5190-4455</td><td></td></tr><tr><td>SureSelect XT Library Prep Kit ILM</td><td>Agilent</td><td>5500-0132</td><td></td></tr><tr><td>T4 ligase (30 units/&amp;micro;L)</td><td>Thermo Scientific</td><td>EL0013</td><td></td></tr><tr><td>table-top Centrifuge 5427 R</td><td>Eppendorf</td><td>5409000012</td><td></td></tr><tr><td>Triton-X-100 (500 mL)</td><td>Merck</td><td>X100-500ML</td><td></td></tr><tr><td>Trypan Blue</td><td>Invitrogen</td><td>T10282</td><td></td></tr><tr><td>Trypsine</td><td>Thermo Scientific</td><td>25300054</td><td></td></tr><tr><td>UltraPure Glycine</td><td>Thermo Scientific</td><td>15527013</td><td></td></tr><tr><td>&amp;beta;-mercaptoethanol</td><td>Thermo Scientific</td><td>31350010</td><td></td></tr></tbody></table>

deciphering high-resolution 3d chromatin organization via capture hi-c

The spatial organization of the genome contributes to its function and regulation in many contexts, including transcription, replication, recombination, and repair. Understanding the exact causality between genome topology and function is therefore crucial and increasingly the subject of intensive research. Chromosome conformation capture technologies (3C) allow inferring the 3D structure of chromatin by measuring the frequency of interactions between any region of the genome. Here we describe a fast and simple protocol to perform Capture Hi-C, a 3C-based target enrichment method that characterizes the allele-specific 3D organization of megabased-sized genomic targets at high-resolution. In Capture Hi-C, target regions are captured by an array of biotinylated probes before downstream high-throughput sequencing. Thus, higher resolution and allele-specificity are achieved while improving the time-effectiveness and affordability of the technology. To demonstrate its strengths, the Capture Hi-C protocol was applied to the mouse X-inactivation center (Xic), the master regulatory locus of X-chromosome inactivation (XCI).

The described Capture Hi-C protocol is based on the preparation of the genome-wide 3C-template using a four-base cutter (DpnII). The subsequent enrichment of ligation fragments across the genomic region of interest is obtained by hybridization of an array of tiling RNA probes and their streptavidin-based capture according to the target enrichment system used in this study (Figure 1). Biotinylated RNA probes were selected as they show tighter binding affinity to their targets compared to DNA probes52,60. Captured libraries are then indexed and pooled for multiplexed high-throughput sequencing. Capture Hi-C data can be visualized as high-resolution Hi-C interaction maps but also as 4C-like single view-point contact maps to specifically visualize the interactions of smaller sequences such as promoters or enhancers within the entire captured region. The workflow of the protocol is shown in Figure 4. Pre-sequencing quality controls are shown in Figure 2 and include the assessment of proper digestion and re-ligation of the 3C template and its efficient shearing and purification across the different steps of the protocol. The sheared 3C template DNA is expected to run between 150 to 700 bp, and no enrichment of fragments &#62;2 kb should be detected. During the following steps, several bead-based DNA cleanup and size selection steps are performed, first after the shearing, then after the pre-capture and post-capture PCRs. Cleaned libraries show a distinct fragment enrichment profile as visualized on a high sensitivity DNA bioanalyzer (Figure 2). The mean fragment size increases over the course of the library preparation due to the ligation of adaptors, sequencing, and indexing primers. Post-sequencing quality controls are obtained via Hi-C Pro and shown in Figure 3. Many different bioinformatics software applications have been proposed for 3C-like data processing and analysis. Among them, the HiC-Pro pipeline is one of the most popular solutions, allowing the processing of raw sequencing data to the final contact maps at various resolutions55. HiC-Pro uses a two-step mapping strategy to align the sequencing reads on the reference genome. The 3C products are then reconstructed and filtered out to remove non-informative pairs of contact and to generate the contact maps. In addition, it is able to use a list of known polymorphisms to perform allele-specific analysis and to separate the contacts coming from the two parental alleles in distinct contact maps. More recently, HiC-Pro has been included and extended into the nf-core framework (nf-core-hic), providing a highly scalable and reproducible community-driven pipeline61,62.
To capture the mouse Xic, an array of 28,913 RNA probes tiling 3 Mb of the X chromosome was designed. This region includes the key player in XCI, the long noncoding gene Xist, and its known ~800 kb regulatory landscape (Figure 5). This ~800 kb region is partitioned into two TADs: one including the Xist promoter and its known positive regulators (i.e., the noncoding transcripts Ftx, Jpx, and Xert and the protein-coding gene Rnf12), and the neighboring TAD encompassing the negative cis-regulators of Xist (i.e., its antisense transcript Tsix, the enhancer element Xite, and the noncoding transcript Linx) (for review44,45).
By applying the described Capture Hi-C protocol to the Xic, the topological organization of this locus was obtained at unprecedented resolution (Figure 6 and Figure 7). This is particularly clear when comparing the Capture Hi-C profile to previously published 5C47 (Figure 6 and Figure 7; Supplementary Table 1) and Hi-C61 (Figure 6 and Figure 7; Supplementary Table 1) profiles. For instance, sub-TAD structures are more evident &#8212; the TAD containing the Xist promoter (Xist-TAD) is clearly subdivided into two smaller domains (Figure 6A, blue arrowhead). Previously, this could only be visually &#34;guessed&#34; from the 5C profile (Figure 6B), albeit the detection of a boundary in this region using the insulation score algorithm. Likewise, the resolution of the Capture Hi-C profile allows the identification of two smaller domains in the neighboring TAD (Figure 6A, B), which contains the promoter of the Tsix locus (Tsix-TAD); this was not previously achieved with 5C (Figure 6B). Of note, topological boundaries determined by the insulation score from the Capture Hi-C and 5C data are generally detected at slightly different locations and with different relative strengths.
Moreover, other sub-TAD structures such as contact loops are clearly visible from the Capture Hi-C data, such as the loop between Xist and Ftx (Figure 7A), previously identified with Capture-C63, and the loop between Xist and Xert (Figure 7B), recently identified using a similar protocol for Capture Hi-C48. Other contacts can also be mapped more precisely due to the increased resolution of the Capture Hi-C profiles, such as those forming the known contact hotspots within the Tsix-TAD between the Linx, Chic1, and Xite loci (Figure 7A).
In comparison with the Hi-C data shown in Figure 7, Capture Hi-C allowed for a fourfold increase in resolution, yet it required only one-fourth of sequencing depth (i.e., 126 M reads versus 571 M) (Supplementary Table 1). This increase in resolution allows for the detection of subTADs and looping interactions that could not be detected by Hi-C at the sequencing depth shown in Figure 6 and Figure 7. The described protocol for Capture Hi-C thus allows for a much more detailed, high-resolution characterization of a large genomic region of interest, when compared to previous approaches.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/64166/64166fig01.jpg" /> 
Figure 1: Probe design. Schematic representation of the strategy used for probe design. Regions of 300 bp upstream and downstream of each DpnII restriction site across the 3 Mb target region were selected and tiled with overlapping biotinylated RNA probes. One of these selected regions is shown, chrX: 102,474,805-102,475,500. No more than 40 bases of repetitive sequences are allowed in each probe. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64166/64166fig01large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/64166/64166fig02.jpg" /> 
Figure 2: Capture Hi-C pre-sequencing quality controls. (A) Representative example of 3C template quality controls. 200 ng of DNA were loaded on a 1% agarose gel. Lane 1: 1 kb ladder. Lane 2: Undigested, cross-linked, and intact chromatin runs as a sharp band at &#62;10 kb. Lane 3: DpnII-digested cross-linked chromatin runs as a smear between 1 kb to 3 kb in size. Lane 4: Final 3C library or template; free ends of digested cross-linked DNA fragments are re-ligated. The DNA smear of lower molecular size is almost undetectable, and the ligation product is detected as a band of &#62;10 kb. (B) Representative examples of high sensitivity bioanalyzer DNA profiles. Top left: successfully sheared 3C library showing a distribution of fragment size between 150 bp and 700 bp. Top right: unsatisfactory sheared 3C library. Unsheared DNA is detected as broad enrichment of fragments &#62;2 kb. (C) Bottom left: sheared DNA sample following a 1:1 left side-size selection using SPRI beads. Fragments of ~300 bp are enriched. Bottom middle: Pre-capture PCR profile after ligation of paired-end adaptors according to the manufacturer&#39;s protocol. Bottom right: final Capture Hi-C library including adaptors, sequencing, and indexing primers for multiplexed sequencing. Abbreviations: bp = base pairs, FU = arbitrary fluorescence unit. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64166/64166fig02large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/64166/64166fig03.jpg" /> 
Figure 3: Capture Hi-C post-sequencing quality controls with HiC-Pro. (A) Example of mapping rate on the reference genome for the first mate of the sequencing pairs. The light blue fraction represents the reads aligned by HiC-Pro and spanning a ligation junction. This metric can thus be used to validate the experimental ligation step. (B) Once sequencing mates are aligned on the genome, only uniquely aligned read pairs are kept for analysis. (C) Non-valid pairs (in red) such as dangling-end, self-circle, or re-ligation are discarded from the analysis. The fraction of valid pairs is a good indicator of the ligation and pull-down efficiency. (D) The valid pairs can be further divided into intra/inter-chromosomal and short/long-range contacts. Duplicated read pairs that are likely to represent PCR artifacts are discarded from the analysis. (E) For allele-specific analysis, HiC-Pro reports the number of allelic reads supported by either one or two mates for each parental genome (i.e., C57BL/6J x CASTEi/J). The same fraction of reads assigned to the maternal and paternal allele are expected. (F) Finally, only valid pairs overlapping the capture region are selected to build the contact maps. Capture-capture pairs represent contacts within the targeted region, while capture-reporter pairs involve interaction between the targeted region and an off-target one. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64166/64166fig03large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/64166/64166fig04.jpg" /> 
Figure 4: Workflow of Capture Hi-C protocol. Schematic representation of different protocol steps. To generate the genome-wide 3C template, chromatin is first cross-linked with formaldehyde and then digested with the DpnII restriction enzyme. Free DNA ends are then re-ligated, cross-linking is reversed, and DNA is purified. To enrich fragments encompassing the target region, an array of biotinylated RNA probes is hybridized to the 3C template and captured by streptavidin-mediated pull-down. Capture libraries are processed for multiplexed sequencing, and valid ligation fragments are quantified to infer the frequency of chromatin contacts across the target, which are visualized as high-resolution interaction maps. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64166/64166fig04large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/64166/64166fig05.jpg" /> 
Figure 5: Overview of the region encompassing the Xic on the mouse X chromosome. Schematic representation of the mouse X chromosome and zoom in of the 3 Mb captured region (ChrX: 102,475,000-105,475,000). The targeted region includes ~800 kb of DNA corresponding to the Xic, the master regulatory locus of XCI. The Xic includes the long noncoding genes, Xist, a key player of XCI, and its regulatory landscape. Positive regulators of Xist are shown in green, and negative regulators in purple. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64166/64166fig05large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 6" class="xfigimg" src="/files/ftp_upload/64166/64166fig06.jpg" /> 
Figure 6: Capture Hi-C, 5C, and Hi-C interaction maps across the 3 Mb captured region. (A) Capture Hi-C interaction map of the 3 Mb target encompassing the mouse Xic at 10 kb resolution (this study). (B)&#160;5C interaction map of the same target region as in A at 6 kb resolution (data reprocessed from47). Repetitive regions not included in the analyses are masked in white. The 5C data require their own bioinformatics processing (see47). After cleaning and alignment, the 5C maps at the primer resolution are binned using a running median (window = 30 kb, step = 5) to reach a final resolution of 6 kb. (C) Hi-C interaction map of the same genomic region as in A and&#160;B at 40 kb resolution (data reprocessed from64). All interaction maps were generated from mouse ESCs. The insulation score was calculated using cooltools and is represented as histograms with insulation minimas at TAD boundaries. TAD boundaries are shown as vertical lines below the map. The height of each line indicates boundary strength. Genes are shown as arrows pointing in the direction of transcription. Sub-TAD boundaries that are detected exclusively or more precisely in Capture Hi-C maps are indicated by magenta and blue arrowheads for sub-TADs in the Tsix and Xist TADs, respectively. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64166/64166fig06large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 7" class="xfigimg" src="/files/ftp_upload/64166/64166fig07.jpg" /> 
Figure 7: Capture Hi-C, 5C, and Hi-C interaction maps across 1 Mb within the captured region. (A) Capture Hi-C interaction map of the 1 Mb genomic region encompassing the mouse Xic at 5 kb resolution (this study). (B) 5C interaction map of the same genomic region as in A. at 6 kb resolution (data reprocessed from47). Repetitive regions not included in the analyses are masked in white. Of note, the 5C data require their own bioinformatics processing (see47). After cleaning and alignment, the 5C maps at the primer resolution are binned using a running median (window = 30 kb, step = 5) to reach a final resolution of 6 kb. (C) Hi-C interaction map of the same genomic region as in A and B of Hi-C at 20 kb resolution (data reprocessed from64). All interaction maps were generated from mESCs. The insulation score was calculated using cooltools and is represented as histograms with insulation minimas at TAD boundaries. TAD boundaries are shown as vertical lines below the map. The height of each line indicates boundary strength. Genes are shown as arrows pointing to the direction of transcription. Contact loops that are detected exclusively or more precisely in Capture Hi-C are indicated by magenta and blue asterisks for loops in the Tsix and Xist TADs, respectively. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64166/64166fig07large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
Supplementary Table 1: Post-sequencing statistics for the datasets used in this manuscript: Capture Hi-C (this study), Hi-C64, and 5C47. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64166/Supplementary Table 1-corrected.xlsx" target="_blank">Please click here to download this File.</a>

Watch this Scientific Journal Video about Deciphering High-Resolution 3D Chromatin Organization via Capture Hi-C at JoVE.com

Deciphering High-Resolution 3D Chromatin Organization via Capture Hi-C

This protocol describes the Capture Hi-C method used to characterize the 3D organization of megabased-sized targeted genomic regions at high-resolution, including boundaries of topologically associating domains (TADs) and long-range chromatin interactions between regulatory and other DNA sequence elements.

Deciphering High-Resolution 3D Chromatin Organization via Capture Hi-C

The spatial organization of the genome contributes to its function and regulation in many contexts, including transcription, replication, recombination, ...

deciphering-high-resolution-3d-chromatin-organization-via-capture-hi-c

Nanyang Technological University

Research

JoVE Journal

Biology

3.2K Views.  EMBL: European Molecular Biology Laboratory. This protocol describes the Capture Hi-C method used to characterize the 3D organization of megabased-sized targeted genomic regions at high-resolution, including boundaries of topologically associating domains (TADs) and long-range chromatin interactions between regulatory and other DNA sequence elements. 

Video: Deciphering High-Resolution 3D Chromatin Organization via Capture Hi-C

Generation of Genome-wide Chromatin Conformation Capture Libraries from Tightly Staged Early Drosophila Embryos

Investigating the three-dimensional architecture of chromatin offers invaluable insight into the mechanisms of gene regulation. Here, we describe a protocol for performing the chromatin conformation capture technique in situ Hi-C on staged Drosophila melanogaster embryo populations. The result is a sequencing library that allows the mapping of all chromatin interactions that occur in the nucleus in a single experiment. Embryo sorting is done manually using a fluorescent stereo microscope and a transgenic fly line containing a nuclear marker. Using this technique, embryo populations from each nuclear division cycle, and with defined cell cycle status, can be obtained with very high purity. The protocol may also be adapted to sort older embryos beyond gastrulation. Sorted embryos are used as inputs for in situ Hi-C. All experiments, including sequencing library preparation, can be completed in five days. The protocol has low input requirements and works reliably using 20 blastoderm stage embryos as input material. The end result is a sequencing library for next generation sequencing. After sequencing, the data can be processed into genome-wide chromatin interaction maps that can be analyzed using a wide range of available tools to gain information about topologically associating domain (TAD) structure, chromatin loops, and chromatin compartments during Drosophila development.

Generation of Genome-wide Chromatin Conformation Capture Libraries from Tightly Staged Early Drosophila Embryos

This work describes a protocol for the generation of high resolution in situ Hi-C libraries from tightly staged pre-gastrulation Drosophila melanogaster embryos.

Investigating the three-dimensional architecture of chromatin offers invaluable insight into the mechanisms of gene regulation. Here, we describe a ...

Genetics

Promoter Capture Hi-C: High-resolution, Genome-wide Profiling of Promoter Interactions

The three-dimensional organization of the genome is linked to its function. For example, regulatory elements such as transcriptional enhancers control the spatio-temporal expression of their target genes through physical contact, often bridging considerable (in some cases hundreds of kilobases) genomic distances and bypassing nearby genes. The human genome harbors an estimated one million enhancers, the vast majority of which have unknown gene targets. Assigning distal regulatory regions to their target genes is thus crucial to understand gene expression control. We developed Promoter Capture Hi-C (PCHi-C) to enable the genome-wide detection of distal promoter-interacting regions (PIRs), for all promoters in a single experiment. In PCHi-C, highly complex Hi-C libraries are specifically enriched for promoter sequences through in-solution hybrid selection with thousands of biotinylated RNA baits complementary to the ends of all promoter-containing restriction fragments. The aim is to then pull-down promoter sequences and their frequent interaction partners such as enhancers and other potential regulatory elements. After high-throughput paired-end sequencing, a statistical test is applied to each promoter-ligated restriction fragment to identify significant PIRs at the restriction fragment level. We have used PCHi-C to generate an atlas of long-range promoter interactions in dozens of human and mouse cell types. These promoter interactome maps have contributed to a greater understanding of mammalian gene expression control by assigning putative regulatory regions to their target genes and revealing preferential spatial promoter-promoter interaction networks. This information also has high relevance to understanding human genetic disease and the identification of potential disease genes, by linking non-coding disease-associated sequence variants in or near control sequences to their target genes.

DNA regulatory elements, such as enhancers, control gene expression by physically contacting target gene promoters, often through long-range chromosomal interactions spanning large genomic distances. Promoter Capture Hi-C (PCHi-C) identifies significant interactions between promoters and distal regions, enabling the assignment of potential regulatory sequences to their target genes.

The three-dimensional organization of the genome is linked to its function. For example, regulatory elements such as transcriptional enhancers control the ...

CRISPR-Mediated Reorganization of Chromatin Loop Structure

Recent studies have clearly shown that long-range, three-dimensional chromatin looping interactions play a significant role in the regulation of gene expression, but whether looping is responsible for or a result of alterations in gene expression is still unknown. Until recently, how chromatin looping affects the regulation of gene activity and cellular function has been relatively ambiguous, and limitations in existing methods to manipulate these structures prevented in-depth exploration of these interactions. To resolve this uncertainty, we engineered a method for selective and reversible chromatin loop re-organization using CRISPR-dCas9 (CLOuD9). The dynamism of the CLOuD9 system has been demonstrated by successful localization of CLOuD9 constructs to target genomic loci to modulate local chromatin conformation. Importantly, the ability to reverse the induced contact and restore the endogenous chromatin conformation has also been confirmed. Modulation of gene expression with this method establishes the capacity to regulate cellular gene expression and underscores the great potential for applications of this technology in creating stable de novo chromatin loops that markedly affect gene expression in the contexts of cancer and development.

Chromatin looping plays a significant role in gene regulation; however, there have been no technological advances that allow for selective and reversible modification of chromatin loops. Here we describe a powerful system for chromatin loop re-organization using CRISPR-dCas9 (CLOuD9), demonstrated to selectively and reversibly modulate gene expression at targeted loci.

Recent studies have clearly shown that long-range, three-dimensional chromatin looping interactions play a significant role in the regulation of gene ...

High-throughput Identification of Gene Regulatory Sequences Using Next-generation Sequencing of Circular Chromosome Conformation Capture (4C-seq)

The identification of regulatory elements for a given target gene poses a significant technical challenge owing to the variability in the positioning and effect sizes of regulatory elements to a target gene. Some progress has been made with the bioinformatic prediction of the existence and function of proximal epigenetic modifications associated with activated gene expression using conserved transcription factor binding sites. Chromatin conformation capture studies have revolutionized our ability to discover physical chromatin contacts between sequences and even within an entire genome. Circular chromatin conformation capture coupled with next-generation sequencing (4C-seq), in particular, is designed to discover all possible physical chromatin interactions for a given sequence of interest (viewpoint), such as a target gene or a regulatory enhancer. Current 4C-seq strategies directly sequence from within the viewpoint but require numerous and diverse viewpoints to be simultaneously sequenced to avoid the technical challenges of uniform base calling (imaging) with next generation sequencing platforms. This volume of experiments may not be practical for many laboratories. Here, we report a modified approach to the 4C-seq protocol that incorporates both an additional restriction enzyme digest and qPCR-based amplification steps that are designed to facilitate a greater capture of diverse sequence reads and mitigate the potential for PCR bias, respectively. Our modified 4C method is amenable to the standard molecular biology lab for assessing chromatin architecture.

High-throughput Identification of Gene Regulatory Sequences Using Next-generation Sequencing of Circular Chromosome Conformation Capture 4C-seq

The identification of physical interactions between genes and regulatory elements is challenging but has been facilitated by chromosome conformation capture methods. This modification to the 4C-seq protocol mitigates PCR bias by minimizing over-amplification of PCR templates and maximizes the mappability of reads by incorporating an addition restriction enzyme digest step.

The identification of regulatory elements for a given target gene poses a significant technical challenge owing to the variability in the positioning and ...

Identification of Enhancer-Promoter Contacts in Embryoid Bodies by Quantitative Chromosome Conformation Capture (4C)

During mammalian development, cell fates are determined through the establishment of regulatory networks that define the specificity, timing, and spatial patterns of gene expression. Embryoid bodies (EBs) derived from pluripotent stem cells have been a popular model to study the differentiation of the main three germ layers and to define regulatory circuits during cell fate specification. Although it is well-known that tissue-specific enhancers play an important role in these networks by interacting with promoters, assigning them to their relevant target genes still remains challenging. To make this possible, quantitative approaches are needed to study enhancer-promoter contacts and their dynamics during development. Here, we adapted a 4C method to define enhancers and their contacts with cognate promoters in the EB differentiation model. The method uses frequently cutting restriction enzymes, sonication, and a nested-ligation-mediated PCR protocol compatible with commercial DNA library preparation kits. Subsequently, the 4C libraries are subjected to high-throughput sequencing and analyzed bioinformatically, allowing detection and quantification of all sequences that have contacts with a chosen promoter. The resulting sequencing data can also be used to gain information about the dynamics of enhancer-promoter contacts during differentiation. The technique described for the EB differentiation model is easy to implement.

Identification of Enhancer-Promoter Contacts in Embryoid Bodies by Quantitative Chromosome Conformation Capture 4C

We report the application of quantitative chromosome conformation capture followed by high-throughput sequencing in embryoid bodies generated from embryonic stem cells. This technique allows to identify and quantitate the contacts between putative enhancers and promoter regions of a given gene during embryonic stem cell differentiation.

During mammalian development, cell fates are determined through the establishment of regulatory networks that define the specificity, timing, and spatial ...

In-Nucleus Hi-C in Drosophila Cells

The genome is organized into topologically associating domains (TADs) delimited by boundaries that isolate interactions between domains. In Drosophila, the mechanisms underlying TAD formation and boundaries are still under investigation. The application of the in-nucleus Hi-C method described here helped to dissect the function of architectural protein (AP)-binding sites at TAD boundaries isolating the Notch gene. Genetic modification of domain boundaries that cause loss of APs results in TAD fusion, transcriptional defects, and long-range topological alterations. These results provided evidence demonstrating the contribution of genetic elements to domain boundary formation and gene expression control in Drosophila. Here, the in-nucleus Hi-C method has been described in detail, which provides important checkpoints to assess the quality of the experiment along with the protocol. Also shown are the required numbers of sequencing reads and valid Hi-C pairs to analyze genomic interactions at different genomic scales. CRISPR/Cas9-mediated genetic editing of regulatory elements and high-resolution profiling of genomic interactions using this in-nucleus Hi-C protocol could be a powerful combination for the investigation of the structural function of genetic elements.

The genome is organized in the nuclear space into different structures that can be revealed through chromosome conformation capture technologies. The in-nucleus Hi-C method provides a genome-wide collection of chromatin interactions in Drosophila cell lines, which generates contact maps that can be explored at megabase resolution at restriction fragment level.

The genome is organized into topologically associating domains (TADs) delimited by boundaries that isolate interactions between domains. In Drosophila, ...

Capturing Chromosome Conformation Across Length Scales

Chromosome conformation capture (3C) is used to detect three-dimensional chromatin interactions. Typically, chemical crosslinking with formaldehyde (FA) is used to fix chromatin interactions. Then, chromatin digestion with a restriction enzyme and subsequent religation of fragment ends converts three-dimensional (3D) proximity into unique ligation products. Finally, after reversal of crosslinks, protein removal, and DNA isolation, DNA is sheared and prepared for high-throughput sequencing. The frequency of proximity ligation of pairs of loci is a measure of the frequency of their colocalization in three-dimensional space in a cell population.
A sequenced Hi-C library provides genome-wide information on interaction frequencies between all pairs of loci. The resolution and precision of Hi-C relies on efficient crosslinking that maintains chromatin contacts and frequent and uniform fragmentation of the chromatin. This paper describes an improved in situ Hi-C protocol, Hi-C 3.0, that increases the efficiency of crosslinking by combining two crosslinkers (formaldehyde [FA] and disuccinimidyl glutarate [DSG]), followed by finer digestion using two restriction enzymes (DpnII and DdeI). Hi-C 3.0 is a single protocol for the accurate quantification of genome folding features at smaller scales such as loops and topologically associating domains (TADs), as well as features at larger nucleus-wide scales such as compartments.

Hi-C 3.0 is an improved Hi-C protocol that combines formaldehyde and disuccinimidyl glutarate crosslinkers with a cocktail of DpnII and DdeI restriction enzymes to increase the signal-to-noise ratio and the resolution of chromatin interaction detection.

Chromosome conformation capture (3C) is used to detect three-dimensional chromatin interactions. Typically, chemical crosslinking with formaldehyde (FA) ...

A Comprehensive Approach to Studying the Promoter Interactome in Scarce Cells

An Integrated Workflow to Study the Promoter-Centric Spatio-Temporal Genome Architecture in Scarce Cell Populations

Spatiotemporal gene transcription is tightly regulated by distal regulatory elements, such as enhancers and silencers, which rely on physical proximity with their target gene promoters to control transcription. Although these regulatory elements are easy to identify, their target genes are difficult to predict, since most of them are cell-type specific and may be separated by hundreds of kilobases in the linear genome sequence, skipping over other non-target genes. For several years, Promoter Capture Hi-C (PCHi-C) has been the gold standard for the association of distal regulatory elements to their target genes. However, PCHi-C relies on the availability of millions of cells, prohibiting the study of rare cell populations such as those commonly obtained from primary tissues. To overcome this limitation, low input Capture Hi-C (liCHi-C), a cost-effective and customizable method to identify the repertoire of distal regulatory elements controlling each gene of the genome, has been developed. liCHi-C relies on a similar experimental and computational framework as PCHi-C, but by employing minimal tube changes, modifying&#160;the reagent concentration and volumes, and swapping or eliminating steps, it accounts&#160;for minimal material loss during library construction. Collectively, liCHi-C enables the study of gene regulation and spatiotemporal genome organization in the context of developmental biology and cellular function.

Author Spotlight: An Integrated Workflow to Study the Promoter-Centric Spatio-Temporal Genome Architecture in Scarce Cell Populations  

Gene expression is regulated by interactions of gene promoters with distal regulatory elements. Here, we descirbe how low input Capture Hi-C (liCHi-C) allows the identification of these interactions in rare cell types, which were previously unmeasurable.

Spatiotemporal gene transcription is tightly regulated by distal regulatory elements, such as enhancers and silencers, which rely on physical proximity ...

Chromosome Conformation Capture (3C) for Undergraduate Research and Teaching

Getting an A with the 3Cs: Chromosome Conformation Capture for Undergraduates

Chromosome conformation capture (3C) is a powerful tool that has spawned a family of similar techniques (e.g., Hi-C, 4C, and 5C, referred to here as 3C techniques) that provide detailed information of the three-dimensional organization of chromatin. The 3C techniques have been used in a wide range of studies, from monitoring the changes in chromatin organization in cancer cells to identifying enhancer contacts made with gene promoters. While many of the studies using these techniques are asking big genome-wide questions with intricate sample types (i.e., single-cell analysis), what is often lost is that the 3C techniques are grounded in basic molecular biology methods that are applicable to a broad range of studies. By addressing tightly focused questions of chromatin organization, this cutting-edge technique can be used to enhance the undergraduate research and teaching lab experience. This paper presents a 3C protocol and provides adaptations and points of emphasis for implementation at primarily undergraduate institutions in undergraduate research and teaching experiences.

Author Spotlight: Getting an A with the 3Cs: Chromosome Conformation Capture for Undergraduates  

Here, we present an adaption of the chromosome conformation capture (3C) technique in detail with an emphasis on undergraduate involvement and learning.

Chromosome conformation capture (3C) is a powerful tool that has spawned a family of similar techniques (e.g., Hi-C, 4C, and 5C, referred to here as 3C ...

Mapping Mammalian 3D Genome Interactions with Micro-C-XL

Three-dimensional (3D) chromosome organization is a major factor in genome regulation and cell-type specification. For example, cis-regulatory elements, known as enhancers, are thought to regulate the activity of distal promoters via interaction in 3D space. Genome-wide chromosome conformation capture (3C)-technologies, such as Hi-C, have transformed our understanding of how genomes are organized in cells. The current understanding of 3D genome organization is limited by the resolution with which the topological organization of chromosomes in 3D space can be resolved. Micro-C-XL measures chromosome folding with resolution at the level of the nucleosome, the basic unit of chromatin, by utilizing micrococcal nuclease (MNase) to fragment genomes during the chromosome conformation capture protocol. This results in an improved signal-to-noise ratio in the measurements, thus facilitating the better detection of insulation sites and chromosome loops compared to other genome-wide 3D technologies. A visually supported, detailed, step-by-step protocol for preparing high-quality Micro-C-XL samples from mammalian cells is presented in this article.

A protocol for mapping the three-dimensional genome organization with nucleosome resolution using the genome-wide chromosome conformation capture method Micro-C-XL is presented here.

Three-dimensional (3D) chromosome organization is a major factor in genome regulation and cell-type specification. For example, cis-regulatory elements, ...

Deciphering High-Resolution 3D Chromatin Organization via Capture Hi-C

Podsumowanie

Przeglądaj więcej filmów

Rozdziały w tym wideo

Generation of Genome-wide Chromatin Conformation Capture Libraries from Tightly Staged Early Drosophila Embryos

Promoter Capture Hi-C: High-resolution, Genome-wide Profiling of Promoter Interactions

CRISPR-Mediated Reorganization of Chromatin Loop Structure

High-throughput Identification of Gene Regulatory Sequences Using Next-generation Sequencing of Circular Chromosome Conformation Capture (4C-seq)

Identification of Enhancer-Promoter Contacts in Embryoid Bodies by Quantitative Chromosome Conformation Capture (4C)

In-Nucleus Hi-C in Drosophila Cells

Capturing Chromosome Conformation Across Length Scales

An Integrated Workflow to Study the Promoter-Centric Spatio-Temporal Genome Architecture in Scarce Cell Populations

Getting an A with the 3Cs: Chromosome Conformation Capture for Undergraduates

Mapping Mammalian 3D Genome Interactions with Micro-C-XL