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

Summary

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

Abstract

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

Introduction

The identification of regulatory elements for gene expression has been facilitated by the Encyclopedia of DNA Elements (ENCODE) Project that comprehensively annotated functional activity for 80% of the human genome^1,2. The identification of the sites for in vivo transcription factor binding, DNaseI hypersensitivity, and epigenetic histone and DNA methylation modifications in individual cell types paved the way for the functional analyses of candidate regulatory elements for target gene expression. Armed with these findings, we are faced with the challenge of determining the functional interconnectivity between regulatory elements and genes. Specifically, what is the relationship between a given target gene and its enhancer(s)? The chromatin conformation capture (3C) method directly addresses this question by identifying physical, and likely functional, interactions between a region of interest and candidate interacting sequences through captured events in fixed chromatin³. As our understanding of chromatin interactions has increased, however, it is clear that the investigation of preselected candidate loci is insufficient to provide a complete understanding of gene-enhancer interactions. For example, ENCODE used the high-throughput chromosomal conformation capture carbon copy (5C) method to examine a small portion of the human genome (1%, pilot set of 44 loci) and reported complex interconnectivity of the loci. Genes and enhancers with identified interactions averaged 2–4 different interacting partners, many of which were hundreds of kilobases away in linear space⁴. Further, Li et al. used Chromatin Interaction Analysis by Paired-End Tag Sequencing (ChIA-PET) to analyze whole-genome promoter interactions and found that 65% of RNA polymerase II binding sites were involved in chromatin interactions. Some of these interactions resulted in large, multi-gene complexes spanning hundreds of kilobases of genomic distance and containing, on average, 8-9 genes each⁵. Together, these findings highlight the need for unbiased whole-genome methods for interrogating chromatin interactions. Some of these methods are reviewed in Schmitt et al.⁶.

More recent methods for chromatin conformation capture studies coupled with next-generation sequencing (Hi-C and 4C-seq) enable the discovery of unknown sequences interacting with a region of interest⁶. Specifically, circular chromosome conformation capture with next-generation sequencing (4C-seq) was developed to identify loci interacting with a sequence of interest in an unbiased manner⁷ by sequencing DNA from captured chromatin proximal to the region of interest in 3D space. Briefly, chromatin is fixed to preserve its native protein-DNA interactions, cleaved with a restriction enzyme, and subsequently ligated under dilute conditions to capture biologically relevant "tangles" of interacting loci (Figure 1). The cross links are reversed to remove the protein, thus leaving the DNA available for additional cleavage with a second restriction enzyme. A final ligation generates smaller circles of interacting loci. Primers to the sequence of interest are then used to generate an amplified library of unknown sequences from the circularized fragments, followed by downstream next generation sequencing.

The protocol presented here, which focuses on sample preparation, makes two major alterations to existing 4C-seq methods^8,9,10,11,12. First, it uses a qPCR-based method to empirically determine the optimal number of amplification cycles for 4C-seq library preparation steps and thus mitigates the potential for PCR bias stemming from over-amplification of libraries. Second, it uses an additional restriction digest step in an effort to reduce the uniformity of known "bait" sequences that hinders accurate base-calling by the sequencing instrument and, hence, maximizes the unique, informative sequence in each read. Other protocols circumvent this issue by pooling many (12-15)⁸ 4C-seq libraries with different bait sequences and/or restriction sites, a volume of experiments which may not be achievable by other laboratories. The modifications presented here allow a small number of experiments, samples, and/or replicates to be indexed and pooled into a single lane.

Protocol

1. Restriction Enzyme Selection

1. Identify a region of interest (e.g., gene promoter, single nucleotide polymorphism (SNP), enhancer) and obtain the DNA sequence from repositories such as National Center for Biotechnology Information (NCBI).
2. Identify the candidate restriction enzymes (REs) for the first restriction enzyme digestion (RE1) that do not cut within the sequence of interest, which produce sticky ends (DNA termini with overhangs) after RE digestion, and whose activities are not inhibited by CpG methylation.
   Note: The resolution of the data is determined by RE1. An enzyme with a 6 base pair (bp) recognition sequence produces, on average, fragments approximately 4 kilobases (kb) in length, while an enzyme with a 4 bp recognition sequence produces fragments approximately 250 bp in length, allowing more precise identification of interacting sequences.
3. Select a RE1 from the candidate enzymes identified in Step 1.2 that produces a “viewpoint” restriction fragment at least 500 bp long that includes the region of interest or, alternatively, a neighboring region.
   NOTE: The enzyme selected must be amenable to primer design (see Step 2).
4. For the second digestion, identify a restriction enzyme (RE2) with a 4 bp recognition sequence that produces sticky ends (DNA termini with overhangs) after RE digestion, whose activity is not inhibited by CpG methylation, and that cuts within the restriction fragment generated by RE1 to produce a 250–500 bp bait fragment (Figure 1).
   NOTE: The enzyme selected must be amenable to primer design (see Step 2).

2. Design and Test Primers for Inverse PCR and Digestion Efficiency qPCR

1. Use a primer design tool such as Primer3 (http://primer3.ut.ee) to design inverse primers that are directed outwards towards the ends of the bait fragment (i.e., the 5’ primer is a “reverse” primer, complementary to the plus strand, while the 3’ primer is a “forward” primer, complementary to the minus strand) and as close to the restriction sites as possible to minimize the amplification of non-informative bait sequence and maximize PCR efficiency (Figure 2A).
   NOTE: Primers should be predicted to have minimal non-specific amplification from genomic DNA by in silico PCR predictions and should not align elsewhere in the genome except their intended target with more than 16/18 identity⁹. As sequencing adapters are not incorporated into the inverse PCR primers, the site of primer binding is more flexible than in other 4C preparation methods; primers should optimally anneal within 50 bp of the end of the corresponding restriction site.
   1. Confirm the specificity of inverse PCR primers by amplification using purified genomic DNA (gDNA) as template.
   2. Optimal primers should yield few products when using gDNA as a template (see Step 11.2 for the description of expected 4C products). If substantial amplification from gDNA occurs, design new primers. If no acceptable primer sets are identified, return to Step 1 and select a new bait by selecting new restriction enzymes.
2. Design qPCR primers to amplify the fragments 70–200 nucleotide (nt) in length to monitor the restriction digestion efficiency (Figure 2B).
   1. Design one pair of primers to amplify across each restriction site (selected in Step 1) that define the bait sequence. Design the primers for both RE1 (see Step 6.5.6) and RE2 sites (see Step 9.2).
   2. Design a set of primers that amplifies a region not containing the sites for either restriction enzyme as a normalization control (uncut DNA) for RE1 and RE2 (see also Step 6.5.6).

3. Collection of Cells

1. Using a method that is appropriate for the tissue or cell culture, obtain a single-cell suspension. Pellet the cells for 5 min at 200 x g, discard the supernatant and resuspend the pellet in 500 µL of 1x phosphate buffered saline (PBS) per 10⁷ cells.

4. Formaldehyde Cross-linking of Cells to Preserve Chromatin Interactions

1. Per 10⁷ cells, add 9.5 mL of 1% electron microscopy (EM)-grade formaldehyde (methanol-free) in PBS and incubate, while tumbling on rocker (or similar), for 10 min at room temperature (RT, 18–22 °C).
   NOTE: Fixation conditions may have to be optimized for each cell type. Prior published work in mouse cells reported that optimal fixation results in at least 40% of mapped reads aligning to the chromosome containing the viewpoint⁹ assuming a non-telomeric region for an average sized chromosome.
2. Transfer the reaction tubes to ice and add ice-cold 1 M glycine to a final concentration of 0.125 M (1.425 mL) to quench the crosslinking reaction. Mix by gentle inversion.
3. Centrifuge for 5 min at 200 x g, 4 °C and carefully remove all the supernatant. Store the cell pellet at -80 °C or proceed immediately to cell lysis.

5. Cell Lysis

1. Resuspend the cell pellet from Step 4.3 in 500 µL of 5 mM ethylenediaminetetraacetic acid (EDTA) to wash. Centrifuge at 200 x g for 5 min at RT. Discard the supernatant.
2. Resuspend the cell pellet in 125 µL of 5 mM EDTA with 0.5% sodium dodecyl sulfate (SDS) and 1x protease inhibitors, such as a 100x cocktail containing 5 mg/mL antipain, 10 mg/mL chymostatin, 10 mg/mL leupeptin, and 10 mg/mL pepstatin A.
   NOTE: Detergent concentration may have to be optimized for each cell type. Insufficient detergent concentrations will result in incomplete cell lysis, leaving chromatin inaccessible to restriction enzymes. If restriction digestion is persistently low in Step 6.5.6 and the activity of the enzyme has been confirmed, additional detergent may be required. Increase SDS concentration in 0.1% increments.
3. Incubate the cell suspension on ice for 10 min.
   NOTE: For primary human keratinocytes, optimal lysis was achieved using a 10 min incubation at 65 °C followed by a 37 °C overnight incubation in a shaking heating block with agitation (900 rpm).
4. Ensure that cell lysis is complete.
   1. Mix 6 µL of the cells with 6 µL of Trypan Blue on a slide and cover with a coverslip. View under a microscope.
      Note: Upon successful lysis, the interior of the cells should appear blue and un-lysed cells will appear white.
   2. If the cell lysis appears insufficient, pellet the cells at 200 x g for 5 min and save the supernatant. Resuspend the pellet and repeat Steps 5.2–5.4 and/or alternatively with different incubation temperatures (see the Note in Step 5.3).
   3. Combine the re-lysed pellet with the saved supernatant and proceed.
      Note: Visual inspection is not a guarantee of sufficient lysis. Digestion efficiency should be objectively determined as in Step 6.5.6.

6. First Restriction Digestion

1. Add 30 μL of 10x restriction enzyme buffer (specified by the manufacturer) and 27 μL of 20% Triton X-100 (final 1.8%) to the suspension from Step 5.4. Bring the total volume to 300 µL with H₂O.
   NOTE: The composition of the restriction enzyme buffer will depend on the RE chosen in Step 1.
2. Remove a 15 μL aliquot as an “Undigested control” and store at 4 °C.
3. Add 200 U RE1 to the remaining reaction mixture and incubate overnight at the temperature appropriate to the enzyme in a shaking heating block with agitation (900 rpm). The next day, add an additional 200 U RE1 and continue the incubation overnight.
4. Remove a 15 μL aliquot as a “Digested control” and store at 4 °C.
5. Determine digestion efficiency:
   1. Add 82.5 µL of 10 mM Tris-HCl pH 7.5 to the 15 µL samples from Steps 6.2 and 6.4. Add 2.5 µL of proteinase K (20 mg/mL) and incubate for 1 h at 65 °C.
   2. Add 100 µL of phenol-chloroform and mix vigorously by inversion to remove residual protein contamination. Centrifuge for 5 min, 16,100 x g at room temperature.
   3. Transfer the aqueous phase to a new tube. Add 6.66 µL of 3 M sodium acetate pH 5.2, 1 µL of 20 mg/mL glycogen, and 300 µL of 100% ethanol (EtOH). Mix gently by inversion and place at -80 °C for 1 h.
   4. Centrifuge for 20 min at 16,100 x g at 4 °C. Remove the supernatant, add 500 µL of 70% EtOH, and centrifuge for 5 min at 16,100 x g at RT.
   5. Remove the supernatant and dry the pellet at room temperature for 2 min. Resuspend the pellet in 50 µL of nuclease-free H₂O.
   6. Determine digest efficiency by qPCR^13,14 using the ∆∆Ct method. Use the primer set not flanking a restriction site (see Step 2.2.2) as the normalization control. Proceed if the digestion efficiency is >85%. Otherwise, pellet the cells and repeat Steps 5.2–6.5, omitting the incubation at 65 °C in Step 5.3.

7. First Ligation

1. Heat-inactivate the restriction enzyme by incubating 20 min at 65 °C. Alternatively, if the enzyme cannot be heat inactivated, phenol-chloroform extract and ethanol precipitate the sample.
   NOTE: The higher inactivation temperatures recommended for some restriction enzymes denature the proteins in the chromatin, and this may negatively affect the sample quality. Additionally, phenol:chloroform extraction is not ideal, as it results in the loss of sample.
2. Transfer to a 50 mL conical tube and add 6 mL of nuclease-free H₂O, 700 µL of 10x ligase buffer (660 mM Tris-HCl pH 7.5, 50 mM magnesium chloride (MgCl₂), 10 mM dithiothreitol (DTT), 10 mM adenosine triphosphate (ATP)), and 50 U T4 DNA Ligase. Mix gently by swirling and incubate overnight at 16 °C.
3. Remove a 100 µL aliquot of the sample as the ‘‘Ligation control’’.
4. Determine ligation efficiency:
   1. Add 2.5 µL of Proteinase K (20mg/mL) and incubate 1 h at 65 °C.
   2. Add 100 µL of phenol-chloroform and mix vigorously by inversion. Centrifuge for 5 min, 16,100 x g at RT.
   3. Transfer the aqueous phase to a new tube. Add 6.66 µL of 3 M sodium acetate pH 5.2, 1 µL of glycogen, and 300 µL of 100% EtOH. Mix gently by inversion and place at -80 °C about 1 h.
   4. Centrifuge for 20 min at 16,100 x g at 4 °C. Remove the supernatant, add 500 µL of 70% EtOH, and centrifuge for 5 min at 16,100 x g at 4 °C.
   5. Remove the supernatant and dry the pellet at room temperature. Resuspend the pellet in 20 µL of water and load on a 0.6% agarose gel next to the “digestion controls” from Step 6.5.
      NOTE: A well-ligated sample should appear as a relatively tight, high molecular weight band (Figure 3).
   6. If the ligation is sufficient, proceed with Step 8. Otherwise, add fresh ATP (1 mM final concentration) and new ligase; incubate overnight at 16 °C and repeat Steps 7.3–7.4.

8. Reverse Cross-linking and Isolate Chromatin

1. Add 15 µL of Proteinase K (20 mg/mL) and incubate overnight at 65 °C to reverse cross-links.
2. Add 30 µL of RNase A (10 mg/mL) and incubate 45 min at 37 °C.
3. Add 7 mL of phenol-chloroform and mix vigorously by inversion. Centrifuge 15 min, 3,300 x g at RT.
4. Transfer the aqueous phase to a new 50 mL tube and add 7.5 mL of nuclease-free H₂O (to dilute the DTT present in the ligase buffer, which otherwise precipitates with the DNA), 1 mL of 3 M sodium acetate pH 5.6, 7 µL of glycogen (20 mg/mL), and 35 mL of 100% EtOH. Mix and incubate at -80 °C for 1 h.
5. Centrifuge 20 min, 3,900 x g at 4 °C. Remove the supernatant (pellet may be difficult to see), wash the pellet with 10 mL of ice-cold 70% EtOH, and centrifuge 15 min, 3,300 x g at 4 °C.
6. Remove the supernatant and briefly dry the pellet at RT. Dissolve the pellet in 150 µL of 10 mM Tris-HCl pH 7.5 at 37 °C. Store at -20 °C or continue with Step 9.

9. Second Restriction Digestion: Trimming the Circles

NOTE: This step creates smaller circles to minimize the overrepresentation of smaller captured fragments due to PCR bias in downstream amplification steps.

1. To the sample from Step 8.6, add 50 µL of 10x restriction enzyme buffer (specified by the manufacturer), 300 µL of nuclease-free H₂O, and 50 U restriction enzyme RE2. Incubate overnight at the temperature appropriate for the chosen enzyme.
2. Remove a 15 µL aliquot of the sample as the “Digestion control”. Determine digestion efficiency as described in Step 6.5.

10. Second Ligation and DNA Purification

1. Inactivate the restriction enzyme as recommended by the manufacturer. If the enzyme cannot be heat-inactivated, remove the enzyme using a column-based purification kit.
   NOTE: As this results in the loss of sample, column purification is not ideal.
2. Transfer the sample to a 50 mL tube and add 12.1 mL of nuclease-free H₂O, 1.4 mL of 10x ligation buffer (660 mM Tris-HCl, pH 7.5; 50 mM MgCl₂; 10 mM DTT; 10 mM ATP), and 100 U T4 DNA ligase. Incubate overnight at 16 °C.
3. Add 467 µL of 3 M sodium acetate pH 5.6, 233 µL nuclease-free H₂O, 7 µL glycogen (20 mg/mL), and 35 mL 100% EtOH. Mix well and incubate at -80 °C 1 h.
4. Centrifuge 45 min, 3,900 x g at 4 °C. Remove the supernatant, add 10 mL of cold 70% EtOH, and centrifuge 15 min, 3,300 x g at 4 °C.
5. Remove the supernatant and briefly dry the pellet at room temperature. Add 150 µL of 10 mM Tris-HCl pH 7.5 and incubate at 37 °C to dissolve the pellet.
6. Purify the samples with a silica column-based PCR purification kit, following the manufacturer’s instructions. Use 1 column per 3 x 10⁶ cells, based on the initial number of the cells. Elute the columns with 50 µL of 10 mM Tris-HCl, pH 7.5 and pool the samples.
7. Measure the concentration by fluorimetry or spectrophotometry using the absorbance at 260 nm. The 4C template is now ready for inverse PCR. Store at -20 °C or proceed directly to Step 11.

11. PCR Amplify Unknown Interacting Sequences by Inverse PCR

1. Determine the linear range of amplification by performing a PCR using template dilutions of 12.5, 25, 50, and 100 ng 4C template. If desired, amplify from gDNA in parallel to directly compare products in order to identify non-specific amplification. Run PCR reactions using an initial denaturation at 94 °C for 2 min; 30 cycles consisting of a denaturation step at 94 °C for 10 s, an annealing step at 55 °C for 1 min, and an extension step at 68 °C for 3 min; and a final extension at 68 °C for 5 min.
2. Separate 15 µL of each PCR product on a 1.5% agarose gel to confirm linear amplification and assess template quality (Figure 4). Amplification from 4C template should yield discrete banding at low DNA concentrations and a smear at higher concentrations. The presence of the smear indicates increased complexity of the amplified 4C ligations.
3. When satisfied with the quality and quantity of the inverse PCR product generated, set up a qPCR to determine the optimal number of cycles to use for amplification:
   1. Set up the reactions containing SYBR and ROX dyes using the reaction mixture in Table 2. Unless the amplification is not linear at high concentrations in Step 11.2, use 100 ng of the template per reaction.
      NOTE: The addition of ROX facilitates the normalization of fluorescent signal from well to well and cycle to cycle.
   2. Run the PCR reactions using an initial denaturation at 94 °C for 2 min; 40 cycles consisting of a denaturation step at 94 °C for 10 s, an annealing step at 55 °C for 1 min, and an extension step at 68 °C for 3 min; and a final extension at 68 °C for 5 min.
   3. Determine the peak (endpoint) fluorescence of the reactions using the amplification plot. Determine the cycle at which reactions reach 25% of peak fluorescence (Figure 5).
      NOTE: This is the number of cycles that will be used to amplify the 4C libraries (see Step 11.4).
4. Set up inverse PCR as in Table 3 to amplify unknown sequences ligated to the bait from the 4C template. Divide into 16 reactions of 50 µL before running. Run PCR reactions using an initial denaturation at 94 °C for 2 min and cycles consisting of a denaturation step at 94 °C for 10 s, an annealing step at 55 °C for 1 min, and an extension step at 68 °C for 3 min. Use the number of cycles determined in Step 11.3.3.
5. Collect and pool the reactions. Purify using a silica column-based PCR purification kit. Use at least 2 columns per 16 reactions. Pool the purified PCR products.
6. Determine the sample quantity and purity by spectrophotometry. Typical yields are between 10 and 20 μg with A260/A280 ~1.85. If the absorption ratios are sub-optimal, re-purify in order to prevent the problems during sequencing.
7. Assess the complexity of library by separating 300 ng of purified PCR product on a 1.5% agarose gel.
   NOTE: Amplified product should resemble that from Step 11.2.

12. Third Restriction Digestion: Trim off Bait Sequences

NOTE: This step removes non-informative bait sequences from the inverse PCR products to maximize informative captured sequences in the downstream sequencing steps. To monitor digest efficiency, a “digest monitor”¹⁵ is digested in parallel using equivalent DNA and enzyme concentration. If RE1 and RE2 are incompatible for simultaneous digestion, for example, due to different optimal incubation temperatures or reaction buffers, this must be done as a sequential digest (this is not ideal).

1. Obtain a RE digest monitor and test the digestion in a 50 µL reaction.
   NOTE: This is a molecule of dsDNA (for example, a plasmid, PCR amplicon, or synthetic DNA) that contains the RE site(s). The only requirement is that it should be easy to distinguish cut from uncut monitor on an agarose gel. Optimize RE monitor mass and enzyme concentration if needed.
2. Digest 1 µg of purified inverse PCR product from Step 11.6 and, in parallel, the RE monitor from Step 12.1.
   NOTE: DNA and enzyme concentration, as well as incubation time, should be the same as for the test digest in Step 12.1 for both the inverse PCR product and the monitor digests. Adjust the reaction volume as needed.
3. Run the digested RE monitor on an agarose gel of concentration appropriate for the expected fragments. Digestion is considered sufficient when <10% of the DNA remains uncut (Figure 6).
4. Purify the digested inverse PCR product on a silica column-based DNA purification kit. If sequential digests are required, repeat Steps 12.1–12.4 with the second enzyme.

13. Preparation of Sequencing Library

1. Ligate the adapters compatible with the respective next generation sequencing platform.
   1. Design the adapters for each RE1 and RE2 so that, after annealing the oligos, each RE adaptor has the respective 1-sided overhang.
      NOTE: For example, if HindIII is used, the annealed adapter should contain a 5’ phosphorylated “AGCT” overhang. Alternatively, if standard T-overhang adapters are used, end-repair and A-tail the library prior to adapter ligation.
   2. Anneal the respective adaptor RE oligos. Resuspend the oligos in annealing buffer (10 mM Tris, pH 7.5, 50 mM sodium chloride (NaCl), 1 mM EDTA), mix in equimolar quantities to 50 µM, heat to 95 °C, and allow to cool slowly to 25 °C.
   3. Perform the ligation reaction. Mix RE-digested 4C library with a total 5- to 10-fold molar excess of annealed adaptors (Step 13.1.2; 50:50 ratio for each RE adaptor used) in 1x ligase buffer and add 6U DNA ligase. Incubate according to manufacturer instructions.
   4. Remove excess adapters by purifying using a low-elution volume silica-based column kit.
2. Size-select.
   NOTE: Size selection can also be performed using a bead-based cleanup kit and is recommended even after column purification.
   1. Cast a 2% agarose gel using a high-resolution agarose, adding a non-UV stain prior to pouring. Ensure that any water lost due to evaporation during heating is replaced by weighing the bottle or flask before and after heating and adding water to recover the lost mass.
   2. Run the adapter-ligated libraries on the gel, leaving empty lanes between the samples.
   3. Excise the gel slices corresponding to the size range 150 bp to 1 kb using a clean scalpel or razor blade.
   4. Purify the libraries from the gel using a gel extraction kit. To minimize GC bias in the recovery of DNA, dissolve the gel slices at room temperature with rotation.
3. Determine the number of cycles for PCR amplification by qPCR using the reaction mixture in Table 4. Run PCR reactions using an initial denaturation at 98 °C for 2 min and 40 cycles consisting of a denaturation step at 98 °C for 10 s, an annealing step at 60 °C for 1 min, and an extension step at 72 °C for 3 min.
   NOTE: The cycle at which fluorescence reaches 25% of the maximum is the number of cycles that will be used to amplify the library, as in Step 11.3.
4. Amplify the libraries using the reaction mixture in Table 5. Run PCR reactions using an initial denaturation at 98 °C for 2 min and cycles consisting of a denaturation step at 98 °C for 10 s, an annealing step at 60 °C for 1 min, and an extension step at 72 °C for 3 min. The number of amplification cycles to be used for each library was determined in Step 13.3.
5. Purify amplified libraries using a low-elution volume silica-based column kit.

14. Sequencing and Analysis of Sequencing Data

1. Determine the concentrations of amplified libraries using a fluorimetric assay. Dilute the libraries to the appropriate concentration for sequencing (check with the sequencing service or core for their recommendation).
2. Determine the quality of the libraries using a microfluidic nucleic acid analysis platform.
   NOTE: The size profile of the library should mirror that seen on the size selection gel in Step 13.2, and the concentration of the sample determined in this step is to be used for sequencing.
3. Pool indexed libraries to obtain at least 3 million reads (at least 50 bp read; single [1X50] or paired end [2X50]) per sample. A high-quality library requires at least 1 million mapped reads⁹, but there will be some attrition during mapping. For example, a sequencing platform that yields approximately 150 million reads per lane can accommodate up to 50 samples in a single lane.

15. Analysis of Sequence Data

1. Demultiplex the data by using index sequences to assign the reads to their appropriate samples.
   NOTE: Depending on their familiarity, users may perform downstream computational analyses of raw FASTA/FASTQ sequence files using basic command line interface or, alternatively, the user-friendly, web-based Galaxy graphical user interface (GUI)¹⁶.
2. Trim adapter sequences and any bait sequence arising from incomplete RE digestion in Step 12, for example, using the FAST-X Toolkit (http://hannonlab.cshl.edu/fastx_toolkit/).
3. Map demultiplexed reads to the appropriate genome of interest (e.g., UCSC mm10, hg19) using the Burrows-Wheeler Aligner (BWA)¹⁷ or other alignment software resulting in SAM output files. If needed, convert SAM files to BAM via samtools¹⁸ (see Step 15.4).
4. Use a 4C-seq analysis pipeline such as Basic4Cseq¹⁹, 4C-ker²⁰, FourCSeq²¹, or fourSig²² to analyze the data and identify the interactions.
   NOTE: Data analysis and the assessment of the quality thereof should be conducted according to the selected software package’s reference manual. Packages generate BED output files except where noted. Briefly, Basic4CSeq¹⁹ uses an input SAM file (Step 15.3) to visualize the interactions (txt, tiff, BED, and wig output files) and assesses the quality of data based on the criteria set forth in van de Werken et al.⁹ 4C-ker²⁰ (input trimmed FASTQ, step 15.2) and FourCSeq²¹ (input BAM, Step 15.3) identify differential interactions between conditions. fourSig²² (input SAM) also identifies significant interactions and prioritizes those that are likely to be reproducible. Tools such as the Genomic Regulatory Enrichment of Annotations Tool (GREAT)²³, or integration with ENCODE data sets may also be used to predict the biological function of interactions discovered by 4C-seq.

Results

Primary human keratinocytes were isolated from 2–3 discarded neonatal foreskins, pooled, and cultured in KSFM supplemented with 30 µg/mL bovine pituitary extract, 0.26 ng/mL recombinant human epidermal growth factor, and 0.09 mM calcium chloride (CaCl₂) at 37 °C, 5% carbon dioxide. The cells were split into two flasks, and one flask was differentiated by the addition of CaCl₂ to a final concentration of 1.2 mM for 72 h. 10⁷ cells each from proliferating an...

Discussion

4C results have the potential to reveal chromatin interactions that can identify previously unknown regulatory elements and/or target genes that are important in a specific biological context^24,25,26. However, technical hurdles may limit the data obtained from these experiments. PCR bias stemming from over-amplification of template in 4C protocols is likely. This protocol addresses this issue by utilizing qPCR to determine the o...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
HindIII	NEB	R0104S
CviQI	NEB	R0639S
DNA oligonucleotide primers	IDT		To be designed by the reader
50 mL conical centrifuge tubes	Fisher Scientific	06-443-19
1.7 mL microcentrifuge tubes	MidSci	AVSS1700
Phosphate buffered saline	Thermo Fisher	14190-136
Formaldehyde, methanol free	Electron Microscopy Sciences	15710
Nutator	VWR	15172-203
Glycine	JT Baker	4059-00
Benchtop centrifuge
Refrigerated microcentrifuge
Ethylenediaminetetraacetic acid (EDTA)	Sigma-Aldrich	ED2SS
20% SDS solution	Sigma-Aldrich	05030
Trypan Blue	Thermo Fisher	15250061
Glass slides	Fisher Scientific	12-550-143
Cover slips	VWR	16004-094
Light microscope
Triton X-100	Alfa Aesar	A16046
Shaking heat block
2 M Tris-HCl	Quality Biological	351-048-101
Proteinase K	NEB	P8107S
Phenol:chloroform:isoamyl alcohol (25:24:1)	Sigma-Aldrich	P2069
Sodium acetate	Sigma-Aldrich	M5661
20 mg/mL glycogen	Thermo Fisher	R0561
Ethanol	Fisher Scientific	04-355-223
Nuclease-free water	Fisher Scientific	MT-46-000-CM
qPCR cycler	Thermo Fisher	4453536
qPCR plates	Thermo Fisher	4309849
Thermocycler	Thermo Fisher	4375786
PCR strip tubes	MidSci	AVSST-FL
1 M Magnesium chloride	Quality Biological	351-033-721
Dithiothreotol	Sigma-Aldrich	43815
Adenosine triphosphate	Sigma-Aldrich	A2383
T4 DNA Ligase	NEB	M0202S
Agarose	Sigma-Aldrich	A6013
RNase A	Thermo Fisher	EN0531
Qiaquick PCR purification kit	Qiagen	28104
MinElute PCR Purification kit	Qiagen	28004
Spectrophotometer
Expand Long Template PCR System	Sigma-Aldrich	11681834001
dNTP mix	Thermo Fisher	R0191
SYBR Green I	Sigma-Aldrich	S9430
ROX	BioRad	172-5858
Sodium chloride	Sigma-Aldrich	S5886
End-It DNA End-Repair Kit	Lucigen	ER0720
LigaFast Rapid DNA Ligation System	Promega	M8221
SYBR Safe	Thermo Fisher	S33102
Taq Polymerase	NEB	M0267S
UltraSieve Agarose	IBI Scientific	IB70054
Qiaquick Gel Extraction Kit	Qiagen	28704