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 border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<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&#39;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/&#181;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>&#946;-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/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/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/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/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/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/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/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/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

Research

JoVE Journal

Biology

3.1K 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

High-resolution Measurement of Odor-Driven Behavior in Drosophila Larvae

Olfactory responses in Drosophila larvae have been traditionally studied in Petri dishes comprising a single peripheral odor source. In this behavioral paradigm, the experimenter usually assumes that the rapid diffusion of odorant molecules from the source leads to the creation of a stable gradient in the dish. To establish a quantitative correlation between sensory inputs and behavioral responses, it is necessary to achieve a more thorough characterization of the odorant stimulus conditions. In this video article, we describe a new method allowing the construction of odorant gradients with stable and controllable geometries. We briefly illustrate how these gradients can be used to screen for olfactory defects (full and partial anosmia) and to study more subtle features of chemotaxis behavior.

In this video article, we describe a new method allowing the construction of odorant gradients with stable and controllable geometries. We briefly illustrate how these gradients can be used to screen for olfactory defects (full and partial anosmia) and to study more subtle features of chemotaxis behavior.

Olfactory responses in Drosophila larvae have been traditionally studied in Petri dishes comprising a single peripheral odor source. In this behavioral ...

Neuronal Cell Cultures from Aplysia for High-Resolution Imaging of Growth Cones

Neuronal growth cones are the highly motile structures at the tip of axons that can detect guidance cues in the environment and transduce this information into directional movement towards the appropriate target cell. To fully understand how guidance information is transmitted from the cell surface to the underlying dynamic cytoskeletal networks, one needs a model system suitable for live cell imaging of protein dynamics at high temporal and spatial resolution. Typical vertebrate growth cones are too small to quantitatively analyze F-actin and microtubule dynamics. Neurons from the sea hare Aplysia californica are 5-10 times larger than vertebrate neurons, can easily be kept at room temperature and are very robust cells for micromanipulation and biophysical measurements. Their growth cones have very defined cytoplasmic regions and a well-described cytoskeletal system. The neuronal cell bodies can be microinjected with a variety of probes for studying growth cone motility and guidance. In the present protocol we demonstrate a procedure for dissection of the abdominal ganglion, culture of bag cell neurons and setting up an imaging chamber for live cell imaging of growth cones.

Aplysia californica neurons develop large growth cones in culture that are excellent for high-resolution imaging of growth cone motility and guidance. Here, we present a protocol for dissection and plating of Aplysia bag cell neurons as well as for setting up a chamber for live cell imaging.

Neuronal growth cones are the highly motile structures at the tip of axons that can detect guidance cues in the environment and transduce this information ...

High-Resolution Video Tracking of Locomotion in Adult Drosophila Melanogaster

Flies provide an important model for studying complex behavior due to the plethora of genetic tools available to researchers in this field. Studying locomotor behavior in Drosophila melanogaster relies on the ability to be able to quantify changes in motion during or in response to a given task. For this reason, a high-resolution video tracking system, such as the one we describe in this paper, is a valuable tool for measuring locomotion in real-time. Our protocol involves the use of an initial air pulse to break the flies momentum, followed by a thirty second filming period in a square chamber. A tracking program is then used to calculate the instantaneous speed of each fly within the chamber in 10 msec increments. Analysis software then compiles this data, and outputs a variety of parameters such as average speed, max speed, time spent in motion, acceleration, etc. This protocol will discuss proper feeding and management of flies for behavioral tasks, handling flies without anesthetization or immobilization, setting up a controlled environment, and running the assay from start to finish.

The study of complex locomotor behavior in Drosophila melanogaster is dependent upon the ability to quantify changes in a given fly's movement. This article demonstrates how to do this using a high-resolution tracking system.

Flies provide an important model for studying complex behavior due to the plethora of genetic tools available to researchers in this field.  Studying ...

Zebrafish Whole Mount High-Resolution Double Fluorescent In Situ Hybridization

Whole mount in situ hybridization is one of the most widely used techniques in developmental biology. Here, we present a high-resolution double fluorescent in situ hybridization protocol for analyzing the precise expression pattern of a single gene and for determining the overlap of the expression domains of two genes. The protocol is a modified version of the standard in situ hybridization using alkaline phosphatase and substrates such as NBT/BCIP and Fast Red 1,2. This protocol utilizes standard digoxygenin and fluorescein labeled probes along with tyramide signal amplification (TSA) 3. The commercially available TSA kits allow flexible experimental design as fluorescence emission from green to far-red can be used in combination with various nuclear stains, such as propidium iodide, or fluorescence immunohistochemistry for proteins. TSA produces a reactive fluorescent substrate that quickly covalently binds to moieties, typically tyrosine residues, in the immediate vicinity of the labeled antisense riboprobe. The resulting staining patterns are high resolution in that subcellular localization of the mRNA can be observed using laser scanning confocal microscopy 3,4. One can observe nascent transcripts at the chromosomal loci, distinguish nuclear and cytoplasmic staining and visualize other patterns such as cortical localization of mRNA. Studies in Drosophila indicate that roughly 70% of mRNAs exhibit specific patterns of subcellular localization that frequently correlate with the function of the encoded protein 5. When combined with computer-aided reconstruction of 3D confocal datasets, our protocol allows the detailed analysis of mRNA distribution with sub-cellular resolution in whole vertebrate embryos.

Zebrafish Whole Mount High-Resolution Double Fluorescent In Situ Hybridization

Whole mount in situ hybridization is one of the most widely used techniques in developmental biology. Here, we present a high-resolution double fluorescent in situ hybridization protocol for analyzing the precise expression pattern of a single gene and for determining the overlap of the expression domains of two genes. We include a propidium iodide nuclear counter-stain to highlight tissue organization.

Whole mount in situ hybridization is one of the most widely used techniques in developmental biology.  Here, we present a high-resolution double ...

Hi-C: A Method to Study the Three-dimensional Architecture of Genomes.

The three-dimensional folding of chromosomes compartmentalizes the genome and and can bring distant functional elements, such as promoters and enhancers, into close spatial proximity 2-6. Deciphering the relationship between chromosome organization and genome activity will aid in understanding genomic processes, like transcription and replication. However, little is known about how chromosomes fold. Microscopy is unable to distinguish large numbers of loci simultaneously or at high resolution. To date, the detection of chromosomal interactions using chromosome conformation capture (3C) and its subsequent adaptations required the choice of a set of target loci, making genome-wide studies impossible 7-10.
We developed Hi-C, an extension of 3C that is capable of identifying long range interactions in an unbiased, genome-wide fashion. In Hi-C, cells are fixed with formaldehyde, causing interacting loci to be bound to one another by means of covalent DNA-protein cross-links. When the DNA is subsequently fragmented with a restriction enzyme, these loci remain linked. A biotinylated residue is incorporated as the 5' overhangs are filled in. Next, blunt-end ligation is performed under dilute conditions that favor ligation events between cross-linked DNA fragments. This results in a genome-wide library of ligation products, corresponding to pairs of fragments that were originally in close proximity to each other in the nucleus. Each ligation product is marked with biotin at the site of the junction. The library is sheared, and the junctions are pulled-down with streptavidin beads. The purified junctions can subsequently be analyzed using a high-throughput sequencer, resulting in a catalog of interacting fragments.
Direct analysis of the resulting contact matrix reveals numerous features of genomic organization, such as the presence of chromosome territories and the preferential association of small gene-rich chromosomes. Correlation analysis can be applied to the contact matrix, demonstrating that the human genome is segregated into two compartments: a less densely packed compartment containing open, accessible, and active chromatin and a more dense compartment containing closed, inaccessible, and inactive chromatin regions. Finally, ensemble analysis of the contact matrix, coupled with theoretical derivations and computational simulations, revealed that at the megabase scale Hi-C reveals features consistent with a fractal globule conformation.

The Hi-C method allows unbiased, genome-wide identification of chromatin interactions (1). Hi-C couples proximity ligation and massively parallel sequencing. The resulting data can be used to study genomic architecture at multiple scales: initial results identified features such as chromosome territories, segregation of open and closed chromatin, and chromatin structure at the megabase scale.

The three-dimensional folding of chromosomes compartmentalizes the genome and and can bring distant functional elements, such as promoters and enhancers, ...

Combined Immunofluorescence and DNA FISH on 3D-preserved Interphase Nuclei to Study Changes in 3D Nuclear Organization

Fluorescent in situ hybridization using DNA probes on 3-dimensionally preserved nuclei followed by 3D confocal microscopy (3D DNA FISH) represents the most direct way to visualize the location of gene loci, chromosomal sub-regions or entire territories in individual cells. This type of analysis provides insight into the global architecture of the nucleus as well as the behavior of specific genomic loci and regions within the nuclear space. Immunofluorescence, on the other hand, permits the detection of nuclear proteins (modified histones, histone variants and modifiers, transcription machinery and factors, nuclear sub-compartments, etc). The major challenge in combining immunofluorescence and 3D DNA FISH is, on the one hand to preserve the epitope detected by the antibody as well as the 3D architecture of the nucleus, and on the other hand, to allow the penetration of the DNA probe to detect gene loci or chromosome territories 1-5. Here we provide a protocol that combines visualization of chromatin modifications with genomic loci in 3D preserved nuclei.

Here we describe a protocol for simultaneous detection of histone modifications by immunofluorescence and DNA sequences by DNA FISH followed by 3D microscopy and analyses (3D immuno-DNA FISH).

Fluorescent in situ hybridization using DNA probes  on 3-dimensionally preserved nuclei followed by 3D confocal microscopy (3D DNA  FISH) represents the ...

Isolation of Specific Genomic Regions and Identification of Associated Molecules by enChIP

The identification of molecules associated with specific genomic regions of interest is required to understand the mechanisms of regulation of the functions of these regions. To enable the non-biased identification of molecules interacting with a specific genomic region of interest, we recently developed the engineered DNA-binding molecule-mediated chromatin immunoprecipitation (enChIP) technique. Here, we describe how to use enChIP to isolate specific genomic regions and identify the associated proteins and RNAs. First, a genomic region of interest is tagged with a transcription activator-like (TAL) protein or a clustered regularly interspaced short palindromic repeats (CRISPR) complex consisting of a catalytically inactive form of Cas9 and a guide RNA. Subsequently, the chromatin is crosslinked and fragmented by sonication. The tagged locus is then immunoprecipitated and the crosslinking is reversed. Finally, the proteins or RNAs that are associated with the isolated chromatin are subjected to mass spectrometric or RNA sequencing analyses, respectively. This approach allows the successful identification of proteins and RNAs associated with a genomic region of interest.

The identification of molecules associated with specific genomic regions of interest is required to understand the mechanisms of regulation of the functions of these regions. This protocol describes procedures to perform engineered DNA-binding molecule-mediated chromatin imunoprecipitation (enChIP) for identification of proteins and RNAs associated with a specific genomic region.

The identification of molecules associated with specific genomic regions of interest is required to understand the mechanisms of regulation of the ...

A High Resolution Method to Monitor Phosphorylation-dependent Activation of IRF3

The IRF3 transcription factor is critical for the first line of defense against pathogens mainly through interferon &#946; and antiviral gene expression. A detailed analysis of IRF3 activation is essential to understand how pathogens induce or evade the innate antiviral response. Distinct activated forms of IRF3 can be distinguished based on their phosphorylation and monomer vs dimer status. In vivo discrimination between the different activated species of IRF3 can be achieved through the separation of IRF3 phosphorylated forms based on their mobility shifts on SDS-PAGE. Additionally, the levels of IRF3 monomer and dimer can be monitored using non-denaturing electrophoresis. Here, we detail a procedure to reach the highest resolution to gain the most information regarding IRF3 activation status. This is achieved through the combination of a high resolution SDS-PAGE and a native-PAGE coupled to immunoblots using multiple total and phosphospecific antibodies. This experimental strategy constitutes an affordable and sensitive approach to acquire all the necessary information for a complete analysis of the phosphorylation-mediated activation of IRF3.

Here we describe a procedure allowing a detailed analysis of the phosphorylation-dependent activation of the IRF3 transcription factor. This is achieved through the combination of a high resolution SDS-PAGE and a native-PAGE coupled to immunoblots using multiple phosphospecific antibodies.

The IRF3 transcription factor is critical for the first line of defense against pathogens mainly through interferon &#946; and antiviral gene expression. ...

Amplification, Next-generation Sequencing, and Genomic DNA Mapping of Retroviral Integration Sites

Retroviruses exhibit signature integration preferences on both the local and global scales. Here, we present a detailed protocol for (1) generation of diverse libraries of retroviral integration sites using ligation-mediated PCR (LM-PCR) amplification and next-generation sequencing (NGS), (2) mapping the genomic location of each virus-host junction using BEDTools, and (3) analyzing the data for statistical relevance. Genomic DNA extracted from infected cells is fragmented by digestion with restriction enzymes or by sonication. After suitable DNA end-repair, double-stranded linkers are ligated onto the DNA ends, and semi-nested PCR is conducted using primers complementary to both the long terminal repeat (LTR) end of the virus and the ligated linker DNA. The PCR primers carry sequences required for DNA clustering during NGS, negating the requirement for separate adapter ligation. Quality control (QC) is conducted to assess DNA fragment size distribution and adapter DNA incorporation prior to NGS. Sequence output files are filtered for LTR-containing reads, and the sequences defining the LTR and the linker are cropped away. Trimmed host cell sequences are mapped to a reference genome using BLAT and are filtered for minimally 97% identity to a unique point in the reference genome. Unique integration sites are scrutinized for adjacent nucleotide (nt) sequence and distribution relative to various genomic features. Using this protocol, integration site libraries of high complexity can be constructed from genomic DNA in three days. The entire protocol that encompasses exogenous viral infection of susceptible tissue culture cells to integration site analysis can therefore be conducted in approximately one to two weeks. Recent applications of this technology pertain to longitudinal analysis of integration sites from HIV-infected patients.

We describe a protocol for amplifying retroviral integration sites from the genomic DNA of infected cells, sequencing the amplified virus-host junctions, and then mapping these sequences to a reference genome. We also describe techniques to quantify the distribution of integration sites relative to various genomic annotations using BEDTools.

Retroviruses exhibit signature integration preferences on both the local and global scales. Here, we present a detailed protocol for (1) generation of ...

High-resolution Optical Mapping of the Mouse Sino-atrial Node

Sino-atrial node (SAN) dysfunctions and associated complications constitute important causes of morbidity in patients with cardiac diseases. The development of novel pharmacological therapies to cure these patients relies on the thorough understanding of both normal physiology and pathophysiology of the SAN. Among the studies of cardiac pacemaking, the mouse SAN is widely used due to its feasibility for modifications in the expression of different genes that encode SAN ion channels or calcium handling proteins. Emerging evidence from electrophysiological and histological studies has also proved the representativeness and similarity of the mouse SAN structure and functions to larger mammals, including the presence of specialized conduction pathways from the SAN to the atrium and a complex pacemakers' hierarchy within the SAN. Recently, the technique of optical mapping has greatly facilitated the exploration and investigation of the origin of excitation and conduction within and from the mouse SAN, which in turn has extended the understanding of the SAN and benefited clinical treatments of SAN dysfunction associated diseases. In this manuscript, we have described in detail how to perform the optical mapping of the mouse SAN from the intact, Langendorff-perfused heart and from the isolated atrial preparation. This protocol is a useful tool to enhance the understanding of mouse SAN physiology and pathophysiology.

Here, we present a protocol for optical mapping of electrical activity from the mouse right atrium and especially the sino-atrial node, at a high spatial and temporal resolution.