Capturing Chromosome Conformation Across Length Scales

Summary

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.

Abstract

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.

Introduction

Chromosome conformation capture has been used since 2002¹. Fundamentally, every conformation capture variant relies on the fixation of DNA-protein and protein-protein interactions to preserve 3D chromatin organization. This is followed by DNA fragmentation, usually by restriction digestion, and, finally, religation of nearby DNA ends to convert spatially proximal loci into unique covalent DNA sequences. Initial 3C protocols used PCR to sample specific, "one-to-one" interactions. Subsequent 4C assays allowed the detection of "one-to-all" interactions², while 5C detected "many-to-many" interactions³. Chromosome conformation capture came to full fruition after implementing next-generation, high-throughput sequencing (NGS), which allowed detection of "all-to-all" genomic interactions using, genome-wide Hi-C⁴ and comparable techniques such as 3C-seq⁵, TCC⁶ and Micro-C^7,8 (see also review by Denker and De Laat⁹).

In Hi-C, biotinylated nucleotides are used to mark 5′ overhangs after digestion and before ligation (Figure 1). This allows for the selection of properly digested and religated fragments using streptavidin-coated beads, setting it apart from GCC¹⁰. An important update to the Hi-C protocol was implemented by Rao et al.¹¹, who performed the digestion and religation in intact nuclei (i.e., in situ) to reduce spurious ligation products. Moreover, substituting HindIII digestion with MboI (or DpnII) digestion reduced the fragment size and increased the resolution potential of Hi-C. This increase allowed the detection of relatively small-scale structures and a more precise genomic localization of points of contact, such as DNA loops between small cis-elements, e.g., loops between CTCF-bound sites generated by loop extrusion^11,12. However, this potential comes at a cost. First, a two-fold increase in resolution requires a four-fold (2²) increase in sequencing reads¹³. Second, the small fragment sizes increase the possibility of mistaking undigested neighboring fragments for digested and religated fragments¹⁴. As mentioned, in Hi-C, digested and religated fragments differ from undigested fragments by the presence of biotin at the ligation junction. However, proper biotin removal from unligated ends is required to assure that only ligation junctions are pulled down^14,15.

With the decreasing cost of NGS, it becomes feasible to study chromosome folding in greater detail. To decrease the size of DNA fragments, and thereby increase resolution, the Hi-C protocol can be adapted to use more frequently cutting restriction enzymes¹⁶ or to use combinations of restriction enzymes^17,18,19. Alternatively, MNase^7,8 in Micro-C and DNase in DNase Hi-C²⁰ can be titrated to achieve optimal digestion.

A recent systematic evaluation of the fundamentals of 3C methods showed that the detection of chromosome folding features at every length scale greatly improved with sequential crosslinking with 1% FA followed by 3 mM DSG¹⁷. Furthermore, Hi-C with HindIII digestion was the best option for detecting large-scale folding features, such as compartments, and that Micro-C was superior at detecting small-scale folding features such as DNA loops. These results led to the development of a single, high-resolution "Hi-C 3.0" strategy, that uses the combination of FA and DSG crosslinkers followed by double digestion with DpnII and DdeI endonucleases²¹. Hi-C 3.0 provides an effective strategy for general use because it accurately detects folding features across all length scales¹⁷. The experimental portion of the Hi-C 3.0 protocol is detailed here and typical results that can be expected after sequencing are shown.

Figure 1: Hi-C procedure in six steps. Cells are fixed first with FA, and then DSG (1). Then, lysis precedes a double digestion with DdeI and DpnII (2). Biotin is added by overhang fill-in and proximal blunted ends are ligated (3) before DNA purification (4). Biotin is removed from unligated ends before sonication and size selection (5). Finally, pull-down of biotin allows for adapter ligation and library amplification by PCR (6). Abbreviations: FA = formaldehyde; DSG = disuccinimidyl glutarate; B = Biotin. Please click here to view a larger version of this figure.

Protocol

1. Fixation by crosslinking

1. Formaldehyde fixation: starting from cells in monolayer
   1. Have cells seeded in appropriate medium to harvest 5 × 10⁶ cells per 150 mm plate.
      NOTE: Users can choose any preferred container that ensures optimal growth of any mammalian cell line. Additionally, cells can be isolated from tissue.
   2. Aspirate the medium with a Pasteur pipette coupled to a vacuum trap from 150 mm plate, wash 2x with ~10 mL HBSS.
   3. Immediately before crosslinking, prepare a 1% FA crosslinking solution in a 50 mL tube by combining 22.5 mL of HBSS and 625 µL of 37% FA to a final 1% concentration. Mix gently by rocking.
      CAUTION: Use fume hood; formaldehyde is toxic.
   4. To crosslink the cells, pour 23.125 mL of the 1% FA solution to each 15 cm plate.
   5. Incubate at room temperature for 10 min and gently rock the plates by hand every 2 min.
   6. Add 1.25 mL of 2.5 M glycine (128 mM final) and gently swirl the plate to quench the crosslinking reaction.
   7. Incubate at room temperature for 5 min and continue incubation on ice for at least 15 min to stop crosslinking.
   8. Scrape the cells from the plates with a cell scraper or rubber policeman.
   9. Transfer the cell suspension to a 50 mL conical tube with a pipet. Centrifuge at 1,000 × g for 10 min at room temperature and discard the supernatant by aspiration.
   10. Wash the cell pellet once with 10 mL of Dulbecco's phosphate-buffered saline (DPBS), using a pipet to resuspend. Then, centrifuge at 1,000 × g for 10 min at room temperature. Proceed immediately to DSG crosslinking.
       NOTE: Be careful when washing the cell pellet as cell pellets can be loose, and cells may be lost.
2. Formaldehyde fixation: starting from cells in suspension
   1. Have cells seeded in appropriate medium to harvest 5 × 10⁶ cells per vessel.
      NOTE: Users can choose any preferred container that ensures optimal cell growth of any mammalian cell line.
   2. Immediately before harvest, count the cells and transfer 5 × 10⁶ cells to a 50 mL conical tube.
   3. Gently pellet the cells by centrifuging at 300 × g for 10 min at room temperature.
   4. Prepare 1% FA crosslinking solution by adding 1.25 mL of 37% FA to 45 mL of HBSS and mix by inverting the tube several times.
      NOTE: Add the entire 1.25 mL of FA without splitting the amount.
      CAUTION: Formaldehyde is highly toxic.
   5. Resuspend the cell pellet in the 46.25 mL of 1% FA crosslinking solution prepared in the previous step by pipetting up and down.
   6. Incubate at room temperature for exactly 10 min on rotator, rocker, or by gentle manual inversion of the tube every 1-2 min.
   7. Quench the crosslinking reaction by adding 2.5 mL of 2.5 M glycine (128 mM final) and mix well by inverting the tube.
   8. Incubate for 5 min at room temperature, and then on ice for at least 15 min to stop crosslinking completely.
   9. Centrifuge at room temperature to pellet the crosslinked cells at 1,000 × g for 10 min and discard the supernatant by aspiration.
   10. Wash the cells once with 10 mL of DPBS, and then centrifuge at 1,000 × g for 10 min at room temperature. Discard the supernatant completely using a pipet and immediately proceed to DSG crosslinking.
       NOTE: Be careful when washing the cell pellet as cell pellets can be loose, and cells may be lost.
3. Crosslinking with disuccinimidyl glutarate
   1. Resuspend the pelleted cells in 9.9 mL of DPBS before adding 100 µL of 300 mM DSG (3 mM final). Mix by inversion.
      NOTE: DSG is moisture-sensitive. It is important to prepare a fresh stock of 300 mM DSG in DMSO on the day of crosslinking.
      CAUTION: DSG in DMSO is highly toxic.
   2. Crosslink the cells at room temperature for 40 min on a rotator.
   3. Add 1.925 mL of 2.5 M glycine (400 mM final), invert to mix, and incubate at room temperature for 5 min.
   4. Centrifuge the cells at 2,000 × g for 15 min at room temperature.
      NOTE: Be careful when removing the supernatant from loose cell pellets.
   5. Resuspend the pellet in 1 mL of 0.05% bovine serum albumin (BSA)-DPBS and transfer to a 1.7 mL tube.
      NOTE: The addition of BSA can help to reduce cell clumping.
   6. Centrifuge the cells at 2,000 × g for 15 min at 4 °C and remove the supernatant by pipet.
      NOTE: To avoid losing the pellet, quickly and completely remove the supernatant.
   7. Snap-freeze the pellet in liquid nitrogen and store at -80 °C or immediately proceed to the next step.

2. Chromosome conformation capture

1. Cell lysis and chromatin digestion
   1. Resuspend (pipet) the crosslinked cell aliquots (~5 × 10⁶ cells) in 1 mL of ice-cold lysis buffer (recipe in Supplemental Table S1) containing 10 µL of protease inhibitor cocktail and transfer to a dounce homogenizer for a 15 min incubation on ice.
      NOTE: Add protease inhibitors to the lysis buffer immediately before use.
   2. Slowly move pestle A up and down 30 times to homogenize the cells on ice and incubate on ice for 1 min to allow the cells to cool down, before another 30 strokes.
   3. Transfer the lysate to a 1.7 mL microcentrifuge tube.
      NOTE: Keep the suspension moving because sometimes cells stick in the pipette tip.
   4. Centrifuge the lysed suspension at 2,500 × g for 5 min at room temperature.
   5. Discard the supernatant and flick or vortex the wet pellet to resuspend. Remove as much of the supernatant as possible to obtain a yoghurt-like substance with minimal clumps.
   6. Resuspend the pellet in 500 µL of ice-cold 1x Restriction Buffer (from 10x; see recipe in Supplemental Table S1) and centrifuge for 5 min at 2,500 × g. Repeat this step for a second wash.
      NOTE: The pellet in 1x Restriction Buffer is more granular than the previous pellet in lysis buffer.
   7. Resuspend the cells in a final volume of 360 µL of 1x Restriction Buffer by pipetting after adding ~340 µL to the carryover volume of the pellet, which depends on the cell size.
   8. Set aside 18 µL of each lysate for testing the Chromatin Integrity (CI). Store the CI samples at 4 °C.
   9. Add 38 µL of 1% sodium dodecyl sulfate (SDS) to each Hi-C tube (total volume of 380 µL) and mix carefully by pipetting without introducing bubbles.
      CAUTION: SDS is toxic.
   10. Incubate the samples at 65 °C without shaking for exactly 10 min to open up the chromatin.
   11. Immediately place the tubes on ice and prepare the Digestion Mixture with Triton X-100 to quench the SDS as described in Table 1.
   12. Add 107 µL of the Digestion Mixture to the Hi-C-tube (487 µL total) to digest the chromatin overnight (~16 h) at 37 °C in a thermomixer with interval shaking (e.g., 900 rpm, 30 s on, 4 min off).
       NOTE: The addition of Triton to a final concentration of 1% serves to quench the SDS.
2. Biotinylation of DNA ends
   1. After overnight digestion, transfer the samples to 65 °C for 20 min to deactivate the remaining endonuclease activity.
   2. During the incubation, prepare a fill-in master mix as shown in Table 2.
   3. After incubation, place the samples immediately on ice.
   4. Set aside a 10 µL Digestion Control (DC) for each sample and store at 4 °C.
   5. Remove the condensation from the lid with a pipet or by spinning. To each sample, add 58 µL of biotin fill-in mix (total sample volume 535 µL) and pipet gently without forming bubbles.
   6. Incubate the samples at 23 °C for 4 h in a thermomixer (e.g., 900 rpm, 30 s on, 4 min off).
3. Ligation of proximal DNA fragments
   1. Prepare the ligation mix as shown in Table 3 while the biotin fill-in is incubating.
   2. Add 665 µL of the ligation mix to each sample (total sample volume 1,200 µL). Mix gently by pipetting.
   3. Incubate the samples at 16 °C for 4 h in a thermomixer with interval shaking (e.g., 900 rpm, 30 s on, 4 min off). Store these samples containing covalently linked chromatin at 4 °C for a few days.
4. Reversal of crosslinking
   1. Bring the volumes of the CI and DC samples to 50 µL with 1x Tris Low EDTA (TLE; see the recipe in Supplemental Table S1).
   2. Add 10 µL of 10 mg/mL Proteinase K to CI and DC samples.
   3. Incubate at 65 °C overnight with interval shaking (e.g., 900 rpm, 30 s on, 4 min off). Alternatively, perform a 30 min reversal of crosslinking for these controls during DNA purification of the Hi-C samples.
   4. To each Hi-C sample, add 50 µL of 10 mg/mL Proteinase K and incubate at 65 °C for at least 2 h with interval shaking (e.g., 900 rpm, 30 s on, 4 min off).
   5. Add another 50 µL of 10 mg/mL Proteinase K to each Hi-C tube (total sample volume 1,300 µL) and continue incubating at 65 °C overnight. Store at 4 °C until DNA purification.
      NOTE: Splitting up Proteinase K incubations ensures total protein digestion.
5. DNA purification
   1. Allow the tubes to cool from 65 °C down to room temperature.
   2. Transfer each sample to a 15 mL conical tube and add a 2.6 mL (2x volume) of phenol:chloroform:isoamyl alcohol to each tube.
      CAUTION: Phenol:chloroform:isoamyl alcohol is a very toxic irritant and is potentially carcinogenic.
   3. Vortex each tube for 1 min and then transfer its contents into a 15 mL phase-lock tube.
   4. Centrifuge the samples for 5 min at maximum speed (1,500-3,500 × g) in a benchtop centrifuge.
   5. Carefully pour the aqueous phase into a 35 mL ultracentrifuge tube and add ultrapure water to a final volume of 1,250 µL.
      NOTE: Use tubes to fit the available ultracentrifuge or split into multiple microcentrifuge tubes.
   6. Add a 1/10^th volume (~125 µL) of 3 M sodium acetate and mix well by inversion.
   7. Add a 2.5x volume (~3.4 mL) of ice-cold 100% ethanol to each sample, balance the tubes for ultracentrifugation by adding ice-cold 100% ethanol, and mix well by inversion.
   8. Incubate the tubes on dry ice for ~15 min (avoid solidification).
   9. Centrifuge the tubes at 18,000 × g for 30 min at 4 °C.
      NOTE: For angled rotors: mark the tubes for where the pellet will be.
   10. Using a pipette, remove and discard the supernatant completely from the non-pellet side.
       NOTE: At this point, the pellet should become visible and can be marked on the tube, because it may not be clearly visible after drying in the next step.
   11. Air-dry the samples for approximately 10 min or until they are visibly dry.
   12. Solubilize each pellet in 450 µL of 1x TLE by pipetting or swirling and transfer to 0.5 mL Centrifugal Filter Unit (CFU) with a 3 kDa molecular weight cut-off.
   13. Centrifuge the CFU at maximum speed for 10 min and discard the flowthrough. Wash each ultracentrifuge tube with an additional 450 µL of 1x TLE and transfer to its CFU for another wash.
       NOTE: Washing the CFU this way limits DNA loss while reducing salt concentration.
   14. Centrifuge the CFU at maximum speed for 10 min and discard the flowthrough.
   15. Add 80 µL of 1x TLE to the column and turn the column around into a new collection tube before centrifuging for 2 min at maximum speed to obtain a final volume of ~100 µL.
   16. Add 1 µL of RnaseA (1 mg/mL; 10-fold dilution of 10 mg/mL stock) to each sample and incubate at 37 °C on a heat block, in a water bath, or a thermomixer for at least 30 min.
   17. After the Rnase treatment, remove the samples from 37 °C and store at 4 °C until the quality control step.
6. Checking chromatin quality, enzyme digestion, and sample ligation
   1. Cool the CI and DC samples to room temperature after reversing the crosslinks in step 2.4.5. Then, transfer to a prespun 2 mL phase-lock tube.
      NOTE: Make sure the phase lock content is centrifuged to a pellet (maximum speed for 2 min).
   2. Add 200 µL of phenol:chloroform:isoamyl alcohol and mix the samples by vortexing for 1 min.
   3. Centrifuge the tubes for 5 min at maximum speed.
   4. Transfer the aqueous phase from each sample (~50 µL) to a new 1.7 mL microfuge tube.
   5. Add 1 µL of Rnase A (from 1 mg/mL) and incubate at 37 °C for at least 30 min.
   6. Load the samples onto a 0.8% agarose gel as recommended in Table 4.
      NOTE: The results expected from a quality control are shown in Figure 2.
   7. Quantify the DNA by densitometry from the gel or using Qubit or Nanodrop.
      NOTE: Accurate quantification ensures the correct input volume in the next part of the protocol. Use several standards with a known quantity to construct a standard curve.

Table 1: Digestion reagents. Please click here to download this Table.

Table 2: Biotin fill-in reagents. *Note that changing enzymes may require different buffers and biotinylated dNTPs. Please click here to download this Table.

Table 3: Ligation mix reagents. Abbreviation: BSA = bovine serum albumin. Please click here to download this Table.

Table 4: Gel loading parameters for quality and size selection assessment. Please click here to download this Table.

Figure 2: Agarose gel showing typical postDNA purification quality control results. (A) The CI control should indicate a band of high molecular weight DNA. (B) The DC and Hi-C samples show a range of DNA sizes. The Hi-C sample, having been combined into larger fragments, should be of higher molecular weight than the DC. The concentration range of markers allows for generating a standard curve. Note, in this example, the CI was loaded on a separate gel, but it is advised to load and run all samples and controls together. Abbreviations: CI = chromatin integrity; DC = digestion control; Hi-C = proximity-ligated. Please click here to view a larger version of this figure.

3. Hi-C sequencing library preparation

1. Removal of biotin from unligated ends
   1. Prepare biotin removal reactions as shown in Table 5.
      NOTE: Typically, 10 µg of DNA is sufficient, but up to 30 µg can be used.
   2. Distribute 2 x 65 µL aliquots from each 130 µL reaction into two PCR tubes.
   3. Transfer to a thermocycler or PCR machine and incubate as described in Table 5.
      NOTE: Samples can be stored here at 4 °C (days to weeks), -20 °C (long term), or immediately carried over to sonication.
2. Sonication
   1. Pool the sample duplicates from the biotin removal step (130 µL total sample volume) into a 130 µL sonicator tube for sonication.
   2. Sonicate the samples using the parameters given in Table 6 to achieve a tight narrow distribution below 500 bp.
      NOTE: Different sonicator types can be used, but for a narrow fragment distribution (100-500 bp), sonicator settings may require optimization.
3. Size selection with magnetic beads
   1. Pipet the sonicated DNA from the sonicator tube(s) into a 1.7 mL low binding tube.
   2. Bring each sample to a total volume of 500 µL with 1x TLE. Try to get the volume as close as possible to 500 µL, since the ratio of sample to magnetic beads mixture is essential for size selection.
   3. Add 400 µL of magnetic beads mixture to each tube to attain a ratio of magnetic beads mixture to the sample volume of 0.8.
      NOTE: Under these conditions, beads capture DNA fragments >300 bp, which will be the Upper Fraction. The supernatant will contain fragments <300 bp, which will be the Lower Fraction.
   4. Mix the tubes by vortexing and incubate for 10 min at room temperature on a rotator. For this and other size selection ratios, make sure that the full volume of the sample mixes well. Look for the "erratic mode" that some rotors have, which works well for smaller volumes.
   5. Incubate for 5 min at room temperature on a magnetic particle separator (MPS).
   6. While incubating, add 500 µL of the magnetic beads mixture to a fresh 1.7 µL low binding tube for each sample.
      NOTE: These tubes will be used for the next size selection of the lower fraction by generating a magnetic beads mixture to sample ratio of 1.1:1.
   7. Leave the tubes on the MPS for 5 min.
   8. Remove the supernatant from the beads and resuspend with 150 µL of the magnetic beads mixture.
      NOTE: This step avoids saturation of the beads with DNA by increasing the number of beads without increasing the volume.
   9. Transfer the supernatant from step 3.3.5 to the labeled tube prepared to select the Lower Fraction (step 3.3.8).
      NOTE: 150 µL of magnetic beads mixture + 400 µL of 0.8x magnetic beads mixture (550 µL total) divided by 500 µL of initial sample = 1.1x magnetic beads mixture to sample ratio.
   10. Mix the Lower Fraction tubes by vortexing and incubate for 10 min at room temperature on a rotator.
       NOTE: The beads will bind DNA fragments >100 bp, resulting in a final bead-bound fraction of 100-300 bp.
   11. Put the Lower Fraction tubes on the MPS for 5 min (room temperature).
   12. Remove the supernatant and centrifuge the tubes briefly to further remove the supernatant as much as possible.
   13. Wash the beads from both Fractions twice using 200 µL of 70% ethanol and reclaim the beads for 5 min on the MPS each time.
   14. After a quick spin in a centrifuge, remove the ethanol completely and further dry the beads on the MPS.
       NOTE: Dry until the alcohol is completely evaporated. The pellet should look like dark chocolate without cracking (can take ~10 min).
   15. Resuspend both fractions in 50 µL of 1x TLE buffer. Incubate at room temperature for 10 min and tap or flick the tubes every other minute to stimulate mixing and elution.
   16. Separate the beads from the supernatant on the MPS for 5 min for both fractions.
   17. Keep the supernatant from each sample. Pipet the supernatant into a 1.7 mL low binding tube.
       NOTE: The samples can be kept at 4 °C for a few days or at -20 °C for long term.
   18. Run a 2% agarose gel as in Table 4 to determine the sample quality and quantity. See Figure 3 for an example of such a gel.
       NOTE: If from experience, sonication is highly reproducible, one can skip this gel and proceed to end repair immediately. It is recommended that users hold on to the Upper Fractions until after titration PCR. Suboptimal DNA quantities that would require a lot of PCR amplification can be rescued from material in the Upper Fraction.
   19. Quantify the amount of DNA from the gel directly after generating a standard curve from the known DNA ladder input or by using a Qubit or Nanodrop.
4. End repair
   1. Prepare the end repair mix as in Table 7 (amount per reaction given).
   2. Transfer the remaining 46 µL of eluted DNA from the Lower Fraction to PCR tubes and add 24 µL of the prepared end repair mix. Incubate in a PCR machine, as proposed in Table 7.
   3. Once the program is completed, keep the samples at 4 °C until pull-down.
5. Pull-down of biotinylated ligation products with streptavidin coated-beads
   1. Determine the amount of streptavidin-coated beads for each library from the quantified size selection (step 3.3.19).
      NOTE: These streptavidin-coated beads (10 mg/mL solution) can bind 20 µg of double-stranded DNA per mg beads (= 20 µg / 100 µL beads). Use 2 µL for each 1 µg of Hi-C DNA but no less than 10 µL.
   2. Mix the streptavidin-coated beads and pipet the volume of beads needed for each library (calculated in the previous step) into individual 1.7 mL low binding tubes.
   3. Resuspend the beads in 400 µL of Tween wash buffer (TWB; see recipe in Supplemental Table S1) and incubate for ~3 min at room temperature on a rotator (see instructions in step 3.3.4).
   4. Separate the beads from the supernatant on the MPS for 1 min and remove the supernatant.
   5. Wash the beads by pipetting another 400 µL of TWB.
   6. Separate the beads from the supernatant on the MPS for 1 min and remove the supernatant.
   7. Add 400 µL of 2x Binding Buffer (BB) (recipe in Supplemental Table S1) to the beads and resuspend. Additionally, add 330 µL of 1x TLE and the solution from End-repair (from step 3.4.3).
   8. Incubate the samples for 15 min at room temperature while mixing on a rotator.
   9. Separate the beads from the supernatant on the MPS for 1 min and remove the supernatant.
   10. Add 400 µL of 1x BB to the beads and resuspend.
   11. Separate the beads from the supernatant on the MPS for 1 min and remove the supernatant.
   12. Add 100 µL of 1x TLE to wash the beads.
   13. Separate the beads from the supernatant on the MPS for 1 min and remove the supernatant.
   14. Finally, add 41 µL of 1x TLE to resuspend the beads.
6. A-tailing
   1. Prepare the A-tailing mix as in Table 8.
   2. Pipet the reactions into PCR tubes and incubate as in Table 8.
   3. Place the PCR tubes on ice immediately after removal from thermocycler and transfer the contents to 1.7 mL low binding tubes.
   4. Separate the beads from the supernatant on the MPS for 1 min and discard the supernatant.
   5. Add 400 µL of 1x ligation buffer, diluted from 5x T4 DNA ligase buffer with ultrapure water.
   6. Separate the beads from the supernatant on the MPS for 1 min, and then discard the supernatant.
   7. Add 1x ligation buffer to a final volume of 40 µL.
7. Annealing adapter oligos
   1. Prepare adapter oligo stocks of 100 µM (Table 9).
      NOTE: Order 250 nmoles of HPLC-purified oligos.
   2. Anneal the adapters in PCR tubes as described in Table 9.
   3. Use a PCR thermocycler to gradually ramp up the temperature at 0.5 °C/s to 97.5 °C. Hold at 97.5 °C for 2.5 min.
   4. Use a PCR thermocycler to gradually ramp up temperature at 0.1 °C/s for 775 cycles (reaching 20 °C). Hold the temperature at 4 °C until further use.
   5. Add 83 µL of 1x Annealing Buffer (recipe in Supplemental Table S1) to dilute the adapters to 15 µM. Store the adapters at -20 °C.
8. Sequencing adaptor ligation
   1. Prepare the adapter ligation mix in a 1.7 mL low binding tube (Table 10).
   2. Ligate for 2 h at room temperature.
   3. Separate the beads from the supernatant on the MPS for 1 min and discard the supernatant.
   4. Add 400 µL of TWB and pipet the beads up and down carefully before incubation on the rotator for 5 min at room temperature. Separate the beads from the supernatant on the MPS and repeat this step once more.
   5. Separate the beads from the supernatant on the MPS (~1 min), discard the supernatant, and add 200 µL of 1x BB.
   6. Separate the beads from the supernatant on the MPS (~1 min) and discard the supernatant.
   7. Add 200 µL of 1x pre-PCR Buffer (from 10x; recipe in Supplemental Table S1) and transfer to a new 1.7 mL low binding tube.
   8. Separate the beads from the supernatant on the MPS (~1 min) and discard the supernatant.
   9. Add 20 µL of 1x pre-PCR Buffer and mix by pipetting.
   10. Keep the tubes on ice while in use or store at 4 °C.
9. Optimization the PCR cycle number by titration
   1. Set up 30 µL of master mix reactions per sample as in Table 11 (amount per reaction given).
   2. Run the lowest number of cycles and take a 5 µL aliquot. Run additional 2-3 cycles on the remaining reaction before taking the next 5 µL aliquot. Repeat to collect four aliquots.
   3. Use PCR parameters from Table 11 for each aliquot.
   4. Add 5 µL of water and 2 µL of 6x dye to each 5 µL sample. Run on a 2% agarose TBE gel (recipe in Supplemental Table S1) with 25-150 ng of Low Molecular weight ladder. See Figure 4 for expected results.
      NOTE: An optimal number of cycles for the final library PCR amplification [step 3.10] is the lowest number of cycles to get visible product on the gel minus one cycle.
10. Final library PCR amplification
    1. Set up 12 x 30 µL reactions to amplify each final library for sequencing as in Table 11.
    2. Cycle the PCR reactions according to Table 11 after determining the number of cycles after PCR titration (step 3.9.3).
    3. Once the PCR is complete, pool the replicate samples into a 1.7 mL low binding microfuge tube.
    4. Place the tubes on the MPS and transfer the supernatant to a new 1.7 mL low binding microfuge tube.
    5. Resuspend the leftover streptavidin-coated beads in 20 µL of 1x pre-PCR Buffer.
       NOTE: This template can be reused when stored at 4 °C for days to weeks or long-term at -20 °C.
11. Removing primers with magnetic beads mixture
    1. Use 1x TLE buffer (recipe in Supplemental Table S1) to adjust the volume from step 3.10.4 to exactly 360 µL.
    2. To each sample, add 360 µL of the magnetic beads mixture and pipette up and down to mix.
    3. On a rotator, mix the samples for 10 min at room temperature.
    4. Separate the beads from the supernatant on the MPS at room temperature (3-5 min).
    5. Wash the beads twice using 200 µL of 70% ethanol and reclaim the beads for 5 min on the MPS each time.
    6. Quick-spin in a centrifuge and completely pipette off the ethanol. Dry the beads in air on the MPS to further evaporate the ethanol.
       NOTE: The pellet should look like dark chocolate without cracking (can take ~10 min).
    7. Add 30 µL of ultrapure water and resuspend to elute the DNA for 10 min at room temperature. Flick the tubes every 2 min to help mixing.
    8. Separate the beads from the supernatant on the MPS for 5 min.
    9. Collect the supernatant from each sample in a fresh 1.7 mL tube.
    10. Run 1 µL of the library on a 2% agarose TBE (recipe in Supplemental Table S1) gel to obtain a fragment size distribution and to quantify the final library (Figure 5).
        NOTE: ClaI digestion can only occur for DpnII-DpnII junctions and serves as a positive ligation control that should result in a lower fragment size distribution for the final library. Final libraries can be stored for a few days at 4 °C, long term at -20 °C, or immediately diluted and submitted for sequencing.

Table 5: Biotin removal reagents and temperatures Please click here to download this Table.

Table 6: Parameters for sonication. Please click here to download this Table.

Table 7: End repair reagents and temperatures. Please click here to download this Table.

Table 8: A-tailing reagents and temperatures. Please click here to download this Table.

Table 9: PCR Primers and Paired-End-adapter oligos with reagents annealing for annealing. Abbreviation: 5PHOS = 5' phosphate. Asterisks indicate phosphorothioated DNA bases. ^#Combine an indexed oligo with the Universal oligo to anneal into an indexed adapter. Please click here to download this Table.

Table 10: Adapter ligation reagents. Please click here to download this Table.

Table 11: PCR reagents and cycling parameters. Please click here to download this Table.

Figure 3: Agarose gel showing typical post-size selection results. The Upper and Lower Fractions for four samples (numbered 1-4) of DpnII-DdeI Hi-C are shown. The first lane for each sample contains the Upper fraction, derived from a 0.8x magnetic beads mixture, and the second and third lanes contain a dilution of the Lower Fraction derived from a 1.1x magnetic beads mixture. Please click here to view a larger version of this figure.

Figure 4: Agarose gel with PCR titration results. Starting from 5 cycles of PCR, samples are taken after every 2 cycles (5, 7, 9, and 11 cycles) for each of four libraries. Based on this figure, 6 cycles was chosen as the optimal cycle for each sample. Please click here to view a larger version of this figure.

Figure 5: Final PCR products. After cleaning and size-selection, PCR products (Hi-C) were loaded next to a ClaI-digested fraction of the same library (ClaI). ClaI-digested fragments indicate the presence of sought-after DpnII-DpnII ligations. Note that ClaI does not digest DpnII-DdeI junctions and, therefore, not all ligations will contribute to a size reduction from this restriction. Please click here to view a larger version of this figure.

Results

The figures in this manuscript were generated from a separate, replicate experiment of the one published previously by Lafontaine et al.²¹. After obtaining high-throughput sequencing data, the Open Chromatin Collective (Open2C: https://github.com/open2c) was used to process the Hi-C data. A similar pipeline can be found on the data portal of the 4D Nucleome project (https://data.4dnucleome.org/resources/data-analysis/hi_c-processing-pipeline). Briefly, the Nextflow pipeline distiller (https://github.com/open2c/distiller-nf) was implemented to (1) align the sequences of Hi-C molecules to the reference genome, (2) parse .sam alignment and form files with Hi-C pairs, (3) filter PCR duplicates, and (4) aggregate pairs into binned matrices of Hi-C interactions. These HDF5 formatted matrices, called coolers, can then be (1) viewed on a HiGlass server (https://higlass.io/) and (2) analyzed using a large set of open source computational tools present in the “cooltools” collection maintained by the Open Chromatin Collective (https://github.com/open2c/cooltools) to extract and quantify folding features such as compartments, TADs, and loops.

Some quality indicators of the Hi-C3.0 libraries can be assessed immediately after mapping read pairs to a reference genome, using a few simple metrics/indicators. First, typically ~50% of sequenced read pairs can be uniquely mapped for human cells. Due to the polymeric nature of chromosomes, most of these mapped reads (~60%-90%) represent interactions within a chromosome (cis), with interaction frequencies rapidly decaying with increasing genomic distance (distance-dependent decay). The distance-dependent decay can be visualized best in a "scaling plot", which shows the contact probability (per chromosome arm) as a function of genomic distance. We found that the use of different crosslinkers and enzymes can alter the distance-dependent decay at long- and short-range distances¹⁷. The addition of DSG crosslinking increases the detectability of interactions at short distances when combined with enzymes such as Mnase and combinations of DpnII-DdeI that produce smaller fragments (Figure 6A).

Distance-dependent decay can also be observed directly from 2D interaction matrices: interactions become more infrequent when located further away from the central diagonal (Figure 6B). Additionally, genomic folding features, such as compartments, TADs, and loops can be identified from Hi-C matrices and scaling plots as deviations from the general genome-wide average distance-dependent decay. Importantly, crosslinking with DSG in addition to FA decreases random ligations, which are unrestricted due to the polymer nature of chromosomes and, therefore, more likely to occur between chromosomes (in trans) (Figure 6C). Reducing random ligation leads to increased signal-to noise ratios, especially for inter-chromosomal and very long range (>10-50 Mb) intrachromosomal interactions.

Figure 6: Representative results of mapped and filtered Hi-C libraries. (A) Scaling plots with contact probability and its derivative for various enzymes, ordered by fragment length (top) and crosslinking with either FA or FA + DSG (bottom). Digestion with MNAse (microC) or DpnII-DdeI (Hi-C 3.0) significantly increases short range contacts (top) as does adding DSG to FA (bottom). (B) Columns show Hi-C heatmaps of DpnII digestion after just FA crosslinking and DpnII or DdeI digestion after FA+DSG crosslinking. White arrows show increasing strength of "dots" after DSG crosslinking and DdeI digestion, implying better detection of DNA loops. Rows show different parts of chromosome 3 at increasing resolution, aligning with panel C: top row: entire chromosome 3 (0-198,295,559 Mb); middle row: 186-196 Mb; bottom row: 191.0-191.5 Mb. (C) Coverage graphs for the regions depicted in A. Black arrows show the lower coverage (%cis reads) for FA-only crosslinking. Abbreviations: FA = formaldehyde; DSG = disuccinimidyl glutarate; chr = chromosome. Please click here to view a larger version of this figure.

Not all mapped reads are useful. A second quality indicator is the number of PCR duplicates. Exact duplicate reads are highly unlikely to occur by chance after ligation and sonication. Thus, such reads likely resulted from PCR amplification and need to be filtered out. Duplicates often arise when too many PCR cycles are required to amplify low complexity libraries. Generally, for Hi-C, most libraries only need 5-8 cycles of final PCR amplification, as determined by titration PCR (see step 3.9; Figure 4). However, libraries with sufficient complexity can be obtained even after 14 cycles of PCR amplification.

Another category of duplicate reads, so-called optical duplicates, can arise from the amplification process on Illumina sequencing platforms that use patterned flow cells (such as HiSeq4000). Optical duplicates are found from either overloading the flow cell, causing (large) clusters to be called two separate clusters, or from local reclustering of the original paired-end molecule after a first round of PCR. Because both types of optical duplicates are local, they can be identified and distinguished from PCR duplicates by their location on the flow cell. Whereas libraries with >15% PCR duplicates would need regeneration, libraries with optical duplicates can be reloaded after optimizing the loading process.

Supplemental Table S1: Buffers and solutions. Please click here to download this Table.

Discussion

Critical steps for cell handling
Although it is possible to use a lower number of input cells, this protocol has been optimized for ~5 × 10⁶ cells per sequencing lane (~400 M reads) to ensure proper complexity after deep sequencing. Cells are best counted prior to fixation. For the generation of ultradeep libraries, we generally multiply the number of lanes (and cells) until the desired read-depth is reached. For optimal fixation, serum-containing medium should be replaced with PBS prior to FA fixation, and fixative solutions should be added immediately and without concentration gradients^15,22. For cell harvesting, scraping is preferred over trypsinization, because the transition from a flatter to a spherical shape after trypsinization could affect the nuclear conformation. After the addition of DSG, loose and clumpy cell pellets are easily lost. Be careful when handling cells at this stage and add up to 0.05% BSA to decrease clumping.

Modifications to the method
This protocol was developed using human cells¹⁷. Yet, based on experience with chromosome conformation capture, this protocol should work for most eukaryotic cells. For a significantly lower input (~1 × 10⁶ cells), we advise using half the volumes for the lysis and conformation capture procedures [steps 2.1-2.4]. This would also allow DNA isolation [step 2.5] to be performed in a tabletop centrifuge with 1.7 mL tubes, which could improve pelleting for low DNA concentrations. The quantification of DNA (step 2.6.6) will indicate how to proceed. For low amounts of isolated DNA (1-5 µg), we suggest skipping the size selection (step 3.3) and proceed with biotin removal after reducing the volume from 130 µL to ~45 µL with a CFU.

This protocol was developed specifically to ensure high-quality data after subsequent crosslinking with FA and DSG and digestion with DpnII and DdeI. However, alternative crosslinking strategies such as FA followed by EGS (ethylene glycol bis(succinimidyl succinate)), which is also used in ChIP-seq²³ and ChIA-PET²⁴, might work equally well¹⁷. Similarly, different enzyme combinations, such as DpnII and HinfI¹⁸ or MboI, MseI, and NlaIII¹⁹ can be used for digestion. When adapting enzyme combinations, be sure to use biotinylated nucleotides that can fill in the specific 5' overhangs and use the most optimal buffers for each cocktail. DpnII comes with its own buffer and the enzyme manufacturer recommends a specific buffer for DdeI digestion. Yet, for the double-digestion with DpnII and DdeI in this protocol, Restriction Buffer is recommended because it is rated at 100% activity for both enzymes.

Troubleshooting conformation capture
The three key steps in chromosome conformation capture: crosslinking, digestion, and religation have all been performed before the results can be visualized on gel. To determine the quality of each of these three steps and discern where problems could have arisen, aliquots before (CI) and after digestion (DC) are taken and loaded on the gel along with the ligated Hi-C sample (Figure 2). This gel is used to determine the quality of the Hi-C sample and whether it will be worth continuing the protocol. Without the CI and DC, it is difficult to pinpoint potential suboptimal step(s). It is worth noting that suboptimal ligation could be due to a problem in the ligation itself, the fill-in, or a problem with crosslinking. To troubleshoot crosslinking, be sure not to use more than 1 × 10⁷ cells per library and start with fresh crosslinking reagents and clean cells (i.e., rinsed with PBS). For ligation, make sure cells and ligation mixture are kept on ice. Add T4 DNA ligase just before the 4 h incubation at 16 °C and mix well.

Troubleshooting library preparation
If more than 10 PCR cycles are needed or no PCR product can be seen on gel after PCR titration (Figure 4), there are a few options to save the Hi-C sample. Working back from the PCR titration, the first option is to try the PCR again. If there is still not enough product, it is possible to attempt another round of A-tailing and adapter ligation (step 3.6) after washing the beads twice with 1x TLE buffer. After this additional A-tailing and adapter ligation, one can proceed to the PCR titration as before. If there is still no product, the last option is to resonicate the 0.8x fraction from step 3.3 and proceed from there.

Limitations and advantages of Hi-C3.0
It is important to realize that Hi-C is a population-based method that captures the average frequency of interactions between pairs of loci in the cell population. Some computational analyses are designed to disentangle combinations of conformations from a population²⁵, but in principle, Hi-C is blind to differences between cells. Although it is possible to perform single-cell Hi-C^26,27 and computational inferences can be made²⁸, single-cell Hi-C is not suitable for obtaining ultrahigh-resolution 3C information. An additional limitation of Hi-C is that it only detects pairwise interactions. To detect multicontact interactions, one can either use frequent cutters combined with short-read sequencing (Illumina)¹⁶ or perform multicontact 3C²⁹ or 4C³⁰, using long-read sequencing from PacBio or Oxford Nanopore platforms. Hi-C derivatives to specifically detect contacts between and along sister chromatids have also been developed^31,32.

Although Hi-C¹⁹ and Micro-C³³ can be used to generate contact maps at subkilobase resolutions, both require a large amount of sequencing reads and this can become a costly undertaking. To get to similar or even higher resolution without the costs, enrichment for specific genomic regions (capture-C³⁴) or specific protein interactions (ChiA-PET³⁵, PLAC-seq³⁶, Hi-ChIP³⁷) can be applied. The strength and downside of these enrichment applications is that only a limited number of interactions are sampled. With such enrichments, the global aspect of Hi-C (and the option of global normalization) is lost.

Importance and potential applications of Hi-C3.0
This protocol was designed to enable high-resolution, ultradeep 3C while simultaneously detecting large-scale folding features such as TADs and compartments¹⁷ (Figure 6). This protocol starts with 5 × 10⁶ cells per tube for each Hi-C library, which should be more than enough material to sequence one or two lanes on a flow cell to obtain up to 1 billion paired-end reads. For ultradeep sequencing, multiple tubes of 5 × 10⁶ cells should be prepared, depending on the number of mapped reads and PCR duplicates. At the highest resolution (<1 kb), looping interactions are mostly found between CTCF sites, but promoter-enhancer interactions can also be detected. Readers can refer to Akgol Oksuz et al.¹⁷ for a detailed description of the data analysis.

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
1 kb Ladder	New England Biolabs	N3232L
Agarose	Invitrogen	16500100
Agencourt AMPure XP magnetic beads , 60 mL	Beckman Coulter	A63881
Amicon Ultra-0.5 Centrifugal Filter Unit (CFU)	EMD Millipore	UFC500396
Annealing Buffer (5x)	See recipe in supplemental materials
ATP 10 mM	ThermoFisher	R0441
Avanti J-25i High Speed Refrigerated ultra-centrifuge	Beckman Coulter
beckman ultracentrifuge tube 35 mL	Beckman Coulter	357002
Binding Buffer (2x)	See recipe in supplemental materials
biotin-14-dATP 0.4 mM	Invitrogen	19524-016
BSA 10 mg/mL	New England Biolabs	B9000S	dilute from 20 mg/mL
Cell scraper	Falcon	353089
Cell scraper	Corning	3008
Conical polypropylene tubes 50 mL	Denville	C1062-P
Conical tube 15 mL	Denville	C1017-P
Covaris micro tube AFA fiber with snap-cap 130 µL	Covaris	520045/520077
Covaris Sonicator	Covaris	E220/E220evolution/M220
Culture flask 175 cm2	Falcon	353112
Culture plates 150 mm x 25 mm	Corning	430599
dATP 1 mM	Invitrogen	56172
dATP 10 mM	Invitrogen	56172
dCTP 10 mM	Invitrogen	56173
DdeI	New England Biolabs	R0175L
dGTP 10 mM	Invitrogen	56174
DMSO	Sigma	D2650-5x10ML
dNTP mix 25 mM	Invitrogen	10297117
Dounce homogenizer	DWK Life Sciences	8853010002/8853030002
DPBS	Gibco	14190-144
DpnII	New England Biolabs	R0543M
DSG	ThermoScientific	20593
dTTP 10 mM	Invitrogen	56175
Ethanol 70%	Fisher	A409-4	Diluted from 100%
Ethidium Bromide	Fisher	BP1302-10
Formaldehyde (37%)	Fisher	BP531-500
Gel loading dye (6x )	New England Biolabs	B7024S
Glycine in ultrapure water 2.5 M	Sigma	G8898-1KG
HBSS	Gibco	14025-092
Igepal CA-630 detergent	MP Biomedicals	198596
Klenow DNA polymerase 5 U/µL	New England Biolabs	M0210L
Klenow Fragment 3-->5’ exo-, 5 U/µL	New England Biolabs	M0212L
ligation buffer (10x)	New England Biolabs	B7203S
Liquid nitrogen
LoBind microcentrifuge tube 1.7 mL	Eppendorf	22431021
Low Molecular Weight DNA Ladder	New England Biolabs	N3233L
Lysis buffer	See recipe in supplemental materials
Magnetic Particle separator	ThermoFisher	12321D
Microfuge tubes 1.7 mL	Axygen	MCT-175-C
MyOne Streptavidin C1 beads	Invitrogen	65001
NEBuffer 2.1 (10x)	New England Biolabs	B7002S
NEBuffer 3.1 (10x)	New England Biolabs	B7203S
PBS	Gibco	70013-032
PCR (strip) tubes	Biorad	TBS0201/ TCS0803
PCR thermocycler	Biorad	T100
Pfu Ultra II Buffer (10x)	Agilent		Comes with Pfu Ultra
PfuUltra II Fusion HS DNA Polymerase	Agilent	600674
Phase lock tube 15 mL	Qiagen	129065
Phase lock tubes 2 mL	Qiagen	129056
Phenol:chloroform:isoamyl alcohol	Invitrogen	15593-049
Protease inhibitor cocktail	ThermoFisher	78440
Proteinase K in ultrapure water 10 mg/mL	Invitrogen	25530-031
Refrigerated Centrifuge	Eppendorf	5810R
RNase A, DNase and protease-free 10 mg/mL	Thermo Scientific	EN0531
Rotator	Argos technologies	EW-04397-40 or rocking platform
SDS 1%	Fisher	BP13111
Sodium acetate pH = 5.2, 3 M	Sigma
Sub-Cell GT Horizontal Electrophoresis System	Biorad	1704401
T4 DNA ligase 1 U/µL	Invitrogen	100004817
T4 DNA polymerase	New England Biolabs	M0203L
T4 DNA polymerase	New England Biolabs	M0203L
T4 DNA polymerase 3 U/µL	New England Biolabs	M0203L
T4 ligation buffer (5x)	Invitrogen	Y90001
T4 polynucleotide kinase 10 U/µL	New England Biolabs	M0201L
Tabletop centrifuge	Eppendorf	5425
TBE buffer	See recipe in supplemental materials
Tris Low EDTA Buffer (TLE)	See recipe in supplemental materials
Triton X-100 (10%)	Sigma	93443
Truseq adapter oligos	Integrated DNA Technologies (IDT))	https://www.idtdna.com/site/order/oligoentry	250 nmole and HPLC purified
Tween 20 detergent	Fisher	9005-64-5
Tween Wash Buffer	See recipe in supplemental materials
Vortex	Scientific Industries	(G560)SI-0236