- Methodology article
- Open Access
A streamlined tethered chromosome conformation capture protocol
BMC Genomics volume 17, Article number: 274 (2016)
Identification of locus-locus contacts at the chromatin level provides a valuable foundation for understanding of nuclear architecture and function and a valuable tool for inferring long-range linkage relationships. As one approach to this, chromatin conformation capture-based techniques allow creation of genome spatial organization maps. While such approaches have been available for some time, methodological advances will be of considerable use in minimizing both time and input material required for successful application.
Here we report a modified tethered conformation capture protocol that utilizes a series of rapid and efficient molecular manipulations. We applied the method to Caenorhabditis elegans, obtaining chromatin interaction maps that provide a sequence-anchored delineation of salient aspects of Caenorhabditis elegans chromosome structure, demonstrating a high level of consistency in overall chromosome organization between biological samples collected under different conditions. In addition to the application of the method to defining nuclear architecture, we found the resulting chromatin interaction maps to be of sufficient resolution and sensitivity to enable detection of large-scale structural variants such as inversions or translocations.
Our streamlined protocol provides an accelerated, robust, and broadly applicable means of generating chromatin spatial organization maps and detecting genome rearrangements without a need for cellular or chromatin fractionation.
The spatial organization of the eukaryotic genome is now accessible through techniques involving massive parallel high-throughput sequencing ([1–3]). An understanding of how chromosomes fold can provide insight into complex relationships between chromatin structure, genetic activity and functional state of the cell ([4, 5]). In addition, genome wide chromatin interaction data sets can reveal long-range information about the grouping and linear organization of sequences along entire chromosomes, enabling high quality chromosome-scale de novo genome assembly .
Chromosome conformation capture (3C)-based techniques have emerged as powerful tools for mapping chromatin contacts ([1, 3, 7–13]). One recently-described technique, tethered conformation capture (TCC) , was developed to improve signal to noise ratio over previously published techniques, allowing in-depth analysis of both intra and inter-chromosomal contacts. 3C-based techniques use proximity ligation and massively parallel sequencing to probe the three-dimensional architecture of chromosomes within the nucleus, with closely interacting regions captured via the ligation step and identified through sequence analysis. In the resulting data sets, the probability of intra-chromosomal contacts is on average much higher than that of inter-chromosomal contacts, as expected if chromosomal territories are at least partially distinct. These data have supported a model in which chromosomal territories are indeed distinct entities: although the probability of interaction decays with linear distance, even loci separated by megabases on the same chromosome are more likely to interact than loci on different chromosomes.
We have developed a rapid tethered conformation capture (RTCC) technique to allow fast application of a TCC-based protocol in C. elegans. Our method allows detection of chromatin contacts using unfractionated whole tissues or whole organisms as starting material, avoiding extensive cellular and molecular fractionation steps. Applied to C. elegans, this protocol illuminates both the large scale structural partitioning (e.g., [15–21]) and fine resolution genome architecture dynamics (i.e. Crane et al. ).
Results and discussion
Detecting genome-wide chromatin contacts using RTCC
RTCC differs from previous chromosome conformation capture protocols in two key respects: (i) the lack of a need for cellular or chromatin fractionation, (ii) the application of efficient transposon tagging approach toward capture of potential ligation junctions.
Figure 1b outlines the RTCC protocol with approximations of time for each step. Of note are the avoidance of cellular or chromatin fractionation steps in standard Hi-C, and the adaptation of the protocol for low input volumes in sequencing library preparation by adapting Nextera tagmentase  to sequencing of biotin-labeled junctions.
To preserve native features of genomic organization, whole worms were flash frozen and finely ground under liquid nitrogen. These frozen samples were quickly resuspended and subjected to formaldehyde treatment to chemically crosslink DNA and proteins. This material could then be used directly for molecular manipulations and analysis (Fig. 1), avoiding any need for an intervening nuclear or chromatin isolation step. In the cross-linked lysates, DNA was digested with a restriction enzyme, proteins were biotinylated non-specifically (cysteine biotinylation) , protein-DNA cross-linked complexes were immobilized at a low surface density on streptavidin-coated beads, and physically juxtaposed free DNA ends were filled in with a nucleotide mixture containing biotin-14-dCTP. Ligation of ends was then performed while fragments remained tethered to the surface of the beads. Following reversal of crosslinking and release of DNA, the Nextera tagmentation (transposon tagging/fragmentation) protocol was used to fragment the DNA and add linkers for sequencing in a single step . The ligation junctions were purified by selection for biotinylated DNA fragments, which were subjected to massively parallel sequencing. Mapping of the sequenced junctions allowed detection of genomic locations for pairs of contacting loci.
We have applied RTCC to a variety of different tissue samples derived from C. elegans, using the DpnII, AvaII, and HindIII restriction enzymes, which cut respectively with 4, 4.5, and 6 base recognition specificities (appropriate restriction enzymes for this protocol must have the properties of leaving a 5' overhang, and of having an overhang for which one base is a “C”.) Exemplary chromatin interaction intensity maps are shown in Fig. 2. We then applied several approaches in evaluating the consistency and accuracy of this large-scale dataset.
First, we evaluated the consistency of results from this approach by comparing results from different experimental replicates. The replicates were carried out with slight deviations in the underlying protocol, distinct stages and tissue distributions, and different choices of restriction enzymes; hence their consistency becomes a test for both biological and technical reproducibility. Indeed we observed that experimental replicates (Additional file 1: Figure S1) were highly correlated (Pearson’s correlation coefficients vary between 0.76 and 0.97 for 50KB binned data).
Second we compared the chromatin interaction matrices obtained with RTCC with a recent Hi-C analysis performed using C. elegans embryos, acquired by Crane et al.  towards understanding of X chromosome topology remodeling during dosage compensation. This comparison likewise demonstrates experiment-to-experiment correlation, both by inspection (Additional file 2: Figure S2) and from calculating a formal Pearson correlation coefficient (>0.64).
Further validation of the RTCC data that we obtained comes from analysis of distances between interacting regions captured by RTCC ligation. With any proximity ligation protocol (even when optimized), we expect a subset of artefactual ligation events that will join DNA segments with no true association (in some cases, even DNA sequences from two different cells). Such events might be expected to lack a clear dependence on inter-locus distance, while bone-fide contacts that were captured from a physiological chromosome configuration would be expected to be much more frequent for closely linked loci.
Our data shows an expected and dramatic inverse association between the distance between the interacting loci and the number of intra-chromosomal contacts (Fig. 3). Dependence of the chromatin contacts on the distance between interacting regions is consistent with DNA polymer-like behavior in which three-dimensional distance between loci increases with increasing genomic distance . Analysis of frequencies of different types of contacts shows greater frequency of intra-chromosomal contacts than inter-chromosomal (Fig. 4) supporting the idea of chromosome territoriality .
Chromatin interaction maps
Using the interaction information we obtained from the mapping of the sequenced DNA library (see Methods), we created a genome-wide raw chromatin interaction matrix. In order to create this matrix, the genome was divided into segments, the size of which depends on the depth of sequencing analysis and level of resolution required (Fig. 2 shows a series of such maps with a 50KB segment length). Each cell in the matrix m i,j corresponds to the number of contacts (proximity ligation products) between segment i and segment j of the genome. The interaction matrix can be depicted visually with a heat map, in which the color intensity correlates with contact frequency (Fig. 2a).
A whole genome raw contacts map for C. elegans shows 6 distinct squares aligned on the diagonal, each corresponding to an individual chromosome, supporting the idea that the probability of intra-chromosomal contacts is on average much higher than that of inter-chromosomal contacts and that chromosomes occupy distinct territories.
The whole genome chromatin interaction matrix was normalized using the expected number of contacts assuming each region has an equal chance of interacting with every other region in the genome, essentially normalizing to read coverage at each region. The resulting normalized whole genome matrix, presented as a heat map in Fig. 2b, exhibits more clear separation into 6 chromosomes and shows some of the intra-chromosomal organization features as well.
Zooming in to the single chromosome level (Fig. 2 c-f) allows study of intra-chromosomal structural organization features ([25, 26]). Tendencies for arm-arm and center-center interaction on the autosomes are evident from two dimensional heat maps that associate each combination of genome positions (x and y) with a color indicating the degree of over- or under-representation in the population of novel junctions. This observation is indicative of a physical basis for the observation that C. elegans autosomes show distinct features in arm and center regions (although both can house active genes, there is a tendency for constraint of contacts between chromosome arms and centers). These findings are consistent with results from linear analyses of chromosomal features ([15–21]), from cytological studies  and from another recently communicated chromosome capture analysis by Crane et al. . Additional file 3: Figure S3 shows a remarkable association between center-arm positioning within each autosome for our data and data from Crane et al. . One physical correlate of the association is suggested by alignment with a dataset derived from immunoprecipitation  using an antibody against the nuclear envelope component LEM-2 (Additional file 3: Figure S3, cyan). As noted also by Crane et al.  our data point to a more complex organization on the X chromosome, with evidence for a domain organization that differs from the end-center-end organization observed on the five autosomes.
Several years ago, it was noted that a subset of DNA segments on C. elegans autosomal arms exhibit a strongly periodic sequence character, with phased runs of A and T residues. Known as the “PATCs” regions, these genomic features are characteristic of introns and other noncoding sequences for a subset of genes expressed in the C. elegans germ line . We found that PATC character (for which we use the quantitative measure defined in ) was strongly associated with partitioning of contacts between central and arm regions of chromosomes (Table 1) This tendency appears to be a general property of this unusual genomic partition, as all five autosomes show similar arm-association-with-arm enrichment on both arms (Table 1). This correspondence suggests a strong tendency for PATC-rich regions to inhabit their own subdomain of the connectivity network, and is consistent with proposed roles of such extended DNA sequence features in long range chromosomal organization .
Modeling-based evaluation suggests utility of RTCC in identification of chromosome structural variation
Information on genomic architecture from Hi-C and similar approaches has been useful in diverse genome-structure applications (e.g., [30, 31]). To evaluate the utility of RTCC in investigating structural variation, we produced three conceptual “model reference genomes” with substantial structural variation, one where the standard N2 reference genome is related by a chromosomal fusion, one where the N2 genome is related by a reciprocal translocation, and one where the N2 genome is related by a large inversion. To model “real world” structural variants, we constructed the translocation with breakpoints at a homologous point in an extended DNA repeat (the transposon Tn5). This rearrangement would have remained undetected with paired end or split read sequencing. As shown in Fig. 5 all of the simulated structural variants could be readily detected using RTCC analysis.
In this study we modified and optimized TCC protocol to work with limited amounts of unfractionated tissue, avoiding subcellular fractionation and chromatin isolation.
The ability to measure proximity between genomic loci in linear space can be utilized both in structural studies of chromosomes and in applications such as genome assembly or phasing.
The RTCC protocol that we present here should provide an efficient and robust means of performing Hi-C experiments in general and in nematodes in particular. Our method allows detection of long-range contacts between genomic loci, supplying useful information to study three-dimensional organization of the genome. Further, Hi-C data sets generated using our RTCC protocol can provide a measurement of proximity between genomic elements in linear sense, allowing this knowledge to be used in bioinformatics applications such as LACHESIS for high quality genome assembly .
Experiments were carried out with two C. elegans cultivars, each derived from the original wild type “N2” isolate used by Brenner , PD7052 is an N2 stock obtained from the Caenorhabditis Genetics Center (Minnesota, USA) in 2001, while PD1074 is a clonal isolate from the genomically defined N2-derivative VC2010 obtained from M. Edgley and colleagues . Germ cell deficient animals were obtained as described below from the temperature sensitive mutant strain glp-1(e2141ts)III .
L1 stage larvae starved and fed animals
To obtain synchronized wild-type L1 stage larvae, animals from N2 (PD7052) strain worms were grown in liquid culture in S-complete media  supplemented with Escherichia coli HB101 bacteria at 20 °C shaking at 180 rounds per minute. Embryos were obtained by standard bleaching protocol and hatched in sterile S-complete liquid media. The animals were starved for 24 h for population synchronization. Half of the synchronized starved L1 stage larvae animals were harvested and frozen in liquid nitrogen [“starved” sample], while the remainder [“fed” sample] were fed on HB101 bacteria for 3 h before harvesting and freezing in liquid nitrogen.
Wild-type N2 and mutant glp-1(e2142ts)III young adult animals
Young adult populations were grown on enriched nematode growth medium plates with Escherichia coli OP50 bacteria  at 23 °C, the worm population was synchronized by standard bleaching protocol  and starvation for 24 h on unseeded nematode growth media plates at 16 °C. Synchronized L1 stage larvae animals were transferred to enriched nematode growth media seeded with Escherichia coli OP50 bacteria  and grown at 23 °C until reaching young adulthood.
In order to obtain germline-depleted populations, we used a temperature sensitive mutant, glp-1(e2141ts)III, which produces a gonad with approximately 20 sperm, in contrast to the thousand or more germ cells present in wild type animals [33, 35]. glp-1(e2141ts)III animals were grown on enriched nematode growth medium plates with Escherichia coli OP50 bacteria at permissive temperature (16 °C). The worm population was synchronized by standard bleaching protocol and starvation for 24 h on unseeded nematode growth media plates at 16 °C. Synchronized L1 stage larvae animals were transferred to enriched nematode growth media seeded with Escherichia coli OP50 bacteria at restrictive temperature (23 °C) . The worms were grown at 23 °C until reaching young adulthood.
Animals were harvested by chilling on ice, centrifugation at 950g for two minutes and washing several times with cold M9. Pellets consisting each of approximately 100μl of closely packed worms were flash frozen in liquid nitrogen and stored in -80 °C. The flash frozen worm pellets were ground to fine powder in liquid nitrogen. The grinding was performed either using mortar and pestle or using an electric drill with “Cellcrusher” drill-bit and “Cellcrusher” base held at liquid nitrogen temperatures [http://cellcrusher.com/tissuepulverizer/, http://cellcrusher.com/drill-bit-2/]. The grinding was done for several minutes, until reaching a fine powder. The powder was stored at -80 °C.
Approximately 100 μl of liquid-nitrogen-ground tissue were used in L1 experiments, corresponding to ~2*106 L1 animals, ~7.4*109 haploid genomic copies, and yielding 5-10 μg of DNA in the final steps. For N2/glp-1 young adult experiments, we also used about 100 μl of worm powder, in this case corresponding to ~104 adult animals, ~7*107 haploid genomic copies  in N2 and about half of that in glp-1. In the final steps of the adult experiments we had 100-350 ng of DNA for Nextera tagmentation.
For both L1 and adult protocols, freeze-ground tissue was processed using procedures modified from . Ground tissue was directly resuspended in 1 ml ice cold buffer A (15 mM Hepes-Na, pH 7.5, 60 mM KCl, 15 mM NaCl, 0.15 mM beta-mercaptoethanol, 0.15 mM spermine, 0.15 mM spermidine, 0.34 M sucrose) containing 1/100 dilution of HALT protease and phosphatase inhibitor cocktail (Thermo Scientific). 16 % Formaldehyde Solution (Thermo Scientific) was added to the final concentration of 1 %. Samples were incubated for 20 min at room temperature while rocking. Formaldehyde was quenched with 120 μl of 2 M Glycine stock solution, samples were incubated at room temperature for 15 min while rocking. Samples were spun for 3 min at maximum speed (15,000 RPM) at room temperature. Pellets were washed twice with 500 μl of ice cold phosphate buffered saline (PBS) containing 1/100 dilution of HALT protease and phosphatase inhibitor cocktail (Thermo Scientific). Worm pellets were resuspended in 500 μl Hi-C lysis buffer (10 mM HEPES pH = 8.0, 10 mM NaCl, 0.2 % IGEPAL CA-630, and 1/100 HALT) and incubated on ice for 15 min. The resulting lysate was spun at maximum speed for 5 min at 4 °C. The supernatant was discarded and pellets were washed with 500 μl of ice-cold wash buffer (50 mM Tris-HCl at pH 8.0, 50 mM NaCl, 1 mM EDTA).
Pellets were resuspended in 500 μl Hi-C nuclear lysis buffer (50 mM Hepes pH = 7.3, 150 mM NaCl, 1 % Triton X-100, 1 mM EDTA, 1 % SDS, 0.1 % Sodium deoxycholate) and rotated at 4 °C for 20 min. Samples were spun at maximum speed for 5 min at 4 °C and supernatant was discarded. Pellets were washed twice with 500 μl ice-cold wash buffer (50 mM Tris-HCl at pH 8.0, 50 mM NaCl, 1mM EDTA) and resuspended in the same buffer to a final volume of 250 μl.
In order to solubilize cross-linked chromatin, samples were mixed with 95 μl of 2 % SDS and incubated at 65 °C for 10 min. Suspensions were cooled down to room temperature before they were mixed with 105 μl 25 mM EZ-Link Iodoacetyl-PEG2-Biotin (IPB) (Thermo Fisher Scientific) to biotinylate proteins. After incubating for 1 h at room temperature while rotating, the SDS was neutralized by adding 1.3 ml 1× NEBuffer 2 (New England Biolabs [NEB], Ipswich, MA, USA). Samples were mixed with 225 μl 10 % Triton X-100 to a final concentration of 1 % and incubated for 10 min on ice, followed by 10 min at 37 °C.
5 μl 1 M DTT, 100 μl 10× NEBuffer 2, 415 μl water and 100 μl of DpnII (or AvaII) restriction enzyme (NEB) (10 U/μl) was added to digest the DNA overnight at 37 °C in a total volume of 2530 μl. After digestion, samples were loaded into a Slide-A-Lyzer Dialysis Cassette G2 (Thermo Fisher Scientific) and dialyzed for 4 h at room temperature against 1 L of dialysis buffer (10 mM Tris-HCl at pH 8.0, 1 mM EDTA) to eliminate excess IPB remaining from the biotinylation step. Dialysis buffer was renewed after 3 h.
400 μl MyOne Streptavidin T1 beads (Life Technologies) were washed 3 times with PBS + 0.01 % Tween-20 (PBST) and beads were resuspended in 2 ml PBST. Dialyzed samples were divided into 5 equal aliquots of 500 μl in 1.5 ml Eppendorf Protein LoBind tubes. 400 μl beads were added to each tube and samples were incubated for 30 min at room temperature while rotating. To prevent interference of unbound streptavidin on the beads with later steps (adding biotinylated dCTP) 5 μl neutralized IPB was added to each tube. IPB was neutralized by adding an equimolar amount of 2-mercaptoethanol. Samples were incubated for an additional 15 min at room temperature while rotating. Non-biotinylated chromatin and non-cross-linked DNA were removed by washing the magnetic T1 beads once with 600 μl PBST and once with 600 μl wash buffer (10 mM Tris-HCl at pH 8.0, 50 mM NaCl, 0.4 % Triton X-100). Beads were resuspended in 100 μl of the same wash buffer. Restriction enzyme (DpnII or AvaII) generated 5’ overhangs were filled in by adding 63 μl water, 1 μl 1 M MgCl, 10 μl 10× NEBuffer 2, 0.7 μl 10 mM dATP, 0.7 μl 10 mM dTTP, 0.7 μl 10 mM 2’-Deoxyguanosine-5’-O-(1-thiotriphosphate), sodium salt, Sp-isomer (Axxora, San Diego, CA, USA), 15 μl 0.4 mM Biotin-14-dCTP (Life Technologies), 4 μl 10 % Triton X-100 and 5 μl 5 U/μl DNA Polymerase I, Large (Klenow) Fragment (NEB).
Samples were incubated for 40 min at room temperature while rotating. Reaction was stopped by adding 5 μl 0.5 M EDTA to the suspension. After 2 min of incubation at room temperature while rotating, beads were washed twice with 600 μl buffer (50 mM Tris-HCl at pH 7.4, 0.4 % Triton X-100, 0.1 mM EDTA) and resuspended in 500 μl of the same buffer. Each sample was transferred into a 15 ml centrifuge tube. For blunt-end ligation under dilute conditions 500 μl sample was mixed with 4 ml water, 250 μl 10× Ligase Buffer (NEB), 100 μl 1 M Tris-HCl at pH 7.4, 90 μl 20 % Triton X-100, 50 μl 100× BSA and 2 μl 2000 U/μl T4 DNA Ligase (NEB), and incubated overnight at 16 °C.
The overnight ligation reaction was stopped by adding 200 μl 0.5M EDTA to each tube. The magnetic T1 beads were collected on the wall of the tube using a magnet and the solution was aspirated out of the tube. The beads were resuspended in 400 μl extraction buffer (50 mM Tris-HCl at pH 8.0, 0.2 % SDS, 1 mM EDTA, 500 mM NaCl) and the mix was transferred into a Eppendorf Protein LoBind microcentrifuge tube. Samples were treated with 5 μl RNase A (20 mg/ml) (Life Technologies) for 45 min at 37 °C and with 20 μl Proteinase K (20 mg/ml) (NEB) overnight at 45 °C.
After overnight incubation, an additional 5 μl Proteinase K were added and samples were incubated for another 2 h at 45 °C. Beads were collected on the wall of the tube and DNA was extracted from the supernatant once with an equal volume of phenol:chloroform (1:1) and once with an equal volume of chloroform. The aqueous phase was mixed with 100 μl of 5 M Ammonium Acetate and 4 μl of 15 mg/ml of glycoblue (Ambion). DNA was precipitated by adding 2.5 volumes of pure ethanol. Precipitated DNA was pelleted by centrifugation at maximal speed (15000RPM) for 30 min at 4 °C. Pellets were washed with ice-cold 70 % ethanol and resuspended in 20 μl 10 mM Tris-HCl at pH 8.0.
100 ng of DNA were subjected to 5 μl Nextera tagmentase (TDE1) at 55 °C for 10 min and purified and concentrated using the DNA Clean & Concentrator™-5 kit (Zymo Research). The purified DNA was eluted in 50 μl of elution buffer. 10 μl of MyOne Streptavidin C1 magnetic beads (Invitrogen) were washed twice with 500 μl 1× Bind & Wash (B&W) buffer (5 mM Tris-HCl at pH 7.4, 0.5 mM EDTA, 1 M NaCl) and resuspended in 50 μl 2× B&W buffer. The purified DNA sample and the C1 beads were mixed and incubated at room temperature for 30 min. The beads were washed once with 500 μl 1× B&W buffer with 0.1 % Triton, once with 500 μl 10 mM Tris-HCl at pH = 8.0 and were resuspended in 20 μl of Resuspension buffer (Nextera DNA Library Prep Kit, Illumina [http://www.illumina.com/products/nextera_dna_library_prep_kit.html]).
The solution with the beads was used directly for PCR amplification according to the Nextera DNA Library Kit protocol, with 12 rounds of PCR. Size selected fragments from a 1 % agarose gel (~500 base pairs) were used for sequencing. The libraries were sequenced using various paired end read lengths with MiSeq, NextSeq and HiSeq Illumina instruments.
Read pairs obtained by massive parallel sequencing were aligned to the C. elegans reference genome (ce10) using an iterative mapping approach utilizing ICE software as described in  and available for download from [https://bitbucket.org/mirnylab/hiclib]. As a starting point for this approach, a first portion of the read is aligned to the reference genome, while the read is truncated to a certain length; subsequently the algorithm aggregates alignments over increasing truncation lengths. The mapping utilizes Bowtie2-2.2.5 software  and allows detection of several types of double sided mapped reads.
The genome was binned into 50KB non-overlapping intervals, and uniquely mapped read pairs were used to create matrices of contacts between 50KB intervals. Binned data for observed contacts between any two regions was normalized to products of sequencing coverage for the two regions, using HOMER v4.7 software .
Raw read counts and bulk read properties from each individual experiment are presented in Tables 2, 3, 4, 5 and 6, which also provides NCBI-GEO accession numbers. All data are available at the NCBI Gene Expression Omnibus (GEO) repository, accession number GSE76930.
chromosome conformation capture
gene expression omnibus
periodic A/T cluster
rapid tethered conformation capture
tethered conformation capture
Lieberman-Aiden E, van Berkum NL, Williams L, Imakaev M, Ragoczy T, Telling A, Amit I, Lajoie BR, Sabo PJ, Dorschner MO, Sandstrom R, Bernstein B, Bender MA, Groudine M, Gnirke A, Stamatoyannopoulos J, Mirny LA, Lander ES, Dekker J. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009;326:289–93.
van Berkum NL, Lieberman-Aiden E, Williams L, Imakaev M, Gnirke A, Mirny LA, Dekker J, Lander ES. Hi-C: a method to study the three-dimensional architecture of genomes. J Vis Exp. 2010;39:1869.
Dekker J, Rippe K, Dekker M, Kleckner N. Capturing chromosome conformation. Science. 2002;295:1306–11.
Misteli T. Beyond the sequence: cellular organization of genome function. Cell. 2007;128:787–800.
Branco MR, Pombo A. Chromosome organization: new facts, new models. Trends Cell Biol. 2007;17:127–34.
Burton JN, Adey A, Patwardhan RP, Qiu R, Kitzman JO, Shendure J. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat Biotechnol. 2013;31:1119–25.
Tolhuis B, Palstra RJ, Splinter E, Grosveld F, de Laat W. Looping and interaction between hypersensitive sites in the active beta-globin locus. Mol Cell. 2002;10:1453–65.
Duan Z, Andronescu M, Schutz K, McIlwain S, Kim YJ, Lee C, Shendure J, Fields S, Blau CA, Noble WS. A three-dimensional model of the yeast genome. Nature. 2010;465:363–7.
Spilianakis CG, Flavell RA. Long-range intrachromosomal interactions in the T helper type 2 cytokine locus. Nat Immunol. 2004;5:1017–27.
Würtele H, Chartrand P. Genome-wide scanning of HoxB1-associated loci in mouse ES cells using an open-ended Chromosome Conformation Capture methodology. Chromosome Res. 2006;14:477–95.
Zhao Z, Tavoosidana G, Sjölinder M, Göndör A, Mariano P, Wang S, Kanduri C, Lezcano M, Sandhu KS, Singh U, Pant V, Tiwari V, Kurukuti S, Ohlsson R. Circular chromosome conformation capture (4C) uncovers extensive networks of epigenetically regulated intra- and interchromosomal interactions. Nat Genet. 2006;38:1341–7.
van Steensel B, Dekker J. Genomics tools for unraveling chromosome architecture. Nat Biotechnol. 2010;28:1089–95.
Simonis M, Klous P, Splinter E, Moshkin Y, Willemsen R, de Wit E, van Steensel B, de Laat W. Nuclear organization of active and inactive chromatin domains uncovered by chromosome conformation capture-on-chip (4C). Nat Genet. 2006;38:1348–54.
Kalhor R, Tjong H, Jayathilaka N, Alber F, Chen L. Genome architectures revealed by tethered chromosome conformation capture and population-based modeling. Nat Biotechnol. 2012;30:90–8.
Brenner S. The genetics of Caenorhabditis elegans. Genetics. 1974;77:71–94.
Gu SG, Fire A. Partitioning the C. elegans genome by nucleosome modification, occupancy, and positioning. Chromosoma. 2010;119:73–87.
Liu T, Rechtsteiner A, Egelhofer TA, Vielle A, Latorre I, Cheung M-S, Ercan S, Ikegami K, Jensen M, Kolasinska-Zwierz P, Rosenbaum H, Shin H, Taing S, Takasaki T, Iniguez AL, Desai A, Dernburg AF, Kimura H, Lieb JD, Ahringer J, Strome S, Liu XS. Broad chromosomal domains of histone modification patterns in C. elegans. Genome Res. 2011;21:227–36.
Barnes TM, Kohara Y, Coulson A, Hekimi S. Meiotic recombination, noncoding DNA and genomic organization in Caenorhabditis elegans. Genetics. 1995;141:159–79.
Rockman MV, Kruglyak L. Recombinational landscape and population genomics of Caenorhabditis elegans. PLoS Genet. 2009;5:e1000419.
Ross JA, Koboldt DC, Staisch JE, Chamberlin HM, Gupta BP, Miller RD, Baird SE, Haag ES. Caenorhabditis briggsae recombinant inbred line genotypes reveal inter-strain incompatibility and the evolution of recombination. PLoS Genet. 2011;7:e1002174.
Garrigues JM, Sidoli S, Garcia BA, Strome S. Defining heterochromatin in C. elegans through genome-wide analysis of the heterochromatin protein 1 homolog HPL-2. Genome Res. 2015;25:76–88.
Crane E, Bian Q, McCord RP, Lajoie BR, Wheeler BS, Ralston EJ, Uzawa S, Dekker J, Meyer BJ. Condensin-driven remodelling of X chromosome topology during dosage compensation. Nature. 2015;523:240–4.
Adey A, Morrison HG, Asan, Xun X, Kitzman JO, Turner EH, Stackhouse B, MacKenzie AP, Caruccio NC, Zhang X, Shendure J. Rapid, low-input, low-bias construction of shotgun fragment libraries by high-density in vitro transposition. Genome Biol. 2010, 11:R119.
Molnar M, Kleckner N. Examination of interchromosomal interactions in vegetatively growing diploid Schizosaccharomyces pombe cells by Cre/loxP site-specific recombination. Genetics. 2008;178:99–112.
Nora EP, Lajoie BR, Schulz EG, Giorgetti L, Okamoto I, Servant N, Piolot T, van Berkum NL, Meisig J, Sedat J, Gribnau J, Barillot E, Blüthgen N, Dekker J, Heard E. Spatial partitioning of the regulatory landscape of the X-inactivation centre. Nature. 2012;485:381–5.
Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, Hu M, Liu JS, Ren B. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012;485:376–80.
Towbin BD, González-Aguilera C, Sack R, Gaidatzis D, Kalck V, Meister P, Askjaer P, Gasser SM. Step-wise methylation of histone H3K9 positions heterochromatin at the nuclear periphery. Cell. 2012;150:934–47.
Ikegami K, Egelhofer TA, Strome S, Lieb JD. Caenorhabditis elegans chromosome arms are anchored to the nuclear membrane via discontinuous association with LEM-2. Genome Biol. 2010;11:R120.
Fire A, Alcazar R, Tan F. Unusual DNA structures associated with germline genetic activity in Caenorhabditis elegans. Genetics. 2006;173:1259–73.
Fudenberg G, Getz G, Meyerson M, Mirny LA. High order chromatin architecture shapes the landscape of chromosomal alterations in cancer. Nat Biotechnol. 2011;29:1109–13.
Simonis N, Rual J-F, Carvunis A-R, Tasan M, Lemmens I, Hirozane-Kishikawa T, Hao T, Sahalie JM, Venkatesan K, Gebreab F, Cevik S, Klitgord N, Fan C, Braun P, Li N, Ayivi-Guedehoussou N, Dann E, Bertin N, Szeto D, Dricot A, Yildirim MA, Lin C, de Smet A-S, Kao H-L, Simon C, Smolyar A, Ahn JS, Tewari M, Boxem M, Milstein S, et al. Empirically controlled mapping of the Caenorhabditis elegans protein-protein interactome network. Nat Methods. 2009;6:47–54.
Thompson O, Edgley M, Strasbourger P, Flibotte S, Ewing B, Adair R, Au V, Chaudhry I, Fernando L, Hutter H, Kieffer A, Lau J, Lee N, Miller A, Raymant G, Shen B, Shendure J, Taylor J, Turner EH, Hillier LW, Moerman DG, Waterston RH. The million mutation project: a new approach to genetics in Caenorhabditis elegans. Genome Res. 2013;23:1749–62.
Priess JR, Schnabel H, Schnabel R. The glp-1 locus and cellular interactions in early C. elegans embryos. Cell. 1987;51:601–11.
Riddle DL, Blumenthal T, Meyer BJ, Priess JR. C. Elegans Ii. Cold Spring Harbor (NY): Cold Spring Harbor Laboratory Press; 1997.
Austin J, Kimble J. glp-1 is required in the germ line for regulation of the decision between mitosis and meiosis in C. elegans. Cell. 1987;51:589–99.
Kodoyianni V, Maine EM, Kimble J. Molecular basis of loss-of-function mutations in the glp-1 gene of Caenorhabditis elegans. Mol Biol Cell. 1992;3:1199–213.
Crittenden SL, Leonhard KA, Byrd DT, Kimble J. Cellular analyses of the mitotic region in the Caenorhabditis elegans adult germ line. Mol Biol Cell. 2006;17:3051–61.
Imakaev M, Fudenberg G, McCord RP, Naumova N, Goloborodko A, Lajoie BR, Dekker J, Mirny LA. Iterative correction of Hi-C data reveals hallmarks of chromosome organization. Nat Methods. 2012;9:999–1003.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, Cheng JX, Murre C, Singh H, Glass CK. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010;38:576–89.
We are grateful to M. Stadler, K. Artiles, J. Pak, L. Hansen, R. Bell, E. Sarinay Cenik, M Shoura, B. Fu, N. Jain, S. Silas, L. Wahba, C. Frokjaer-Jensen and J. Arribere for help and discussion.
Some strains were provided by the CGC, which is funded by NIH Office of Research Infrastructure Programs (P40 OD010440). This work was supported by a Human Frontiers Science Program Fellowship (LT000517/2011) to IG and by NIH grants R01GM37706 to AF and R01GM53804 to AV.
The authors declare that they have no competing interests.
IG conceived of the study, participated in the study design, carried out chromatin conformation capture experiments and data analysis and wrote the manuscript. AZF conceived of the study, and participated in its design and coordination and helped to draft the manuscript. AMV participated in the design of the study and helped to draft the manuscript. SR participated in chromatin conformation capture experiments and data analysis and helped to draft the manuscript. All authors read and approved the final manuscript.
Correlation between different experiments. The genome was divided into 50KB non-overlapping segments and chromatin contact matrices were generated. In order to evaluate the correlation between experiments we calculated the normalized levels of contacts for all matrix locations (dividing by the total number of contacts detected). Then we calculated correlations between normalized values of contacts of different experiments. x-axis and y-axis stand for normalized contacts values calculated for different experiments. The correlation analysis was performed between experiments done with four C. elegans populations:(i) N2 L1 starved animals, (ii) N2 L1 fed animals, (iii) N2 young adult animals,(iv) glp-1(e2141ts)III [33, 35] young adults (populations grown at the permissive temperature (16 °C) to L1 stage, then shifted to the restrictive temperature (23 °C)  to adulthood). Analysis was performed with DpnII and AvaII restriction enzymes as noted. Any contacts between any location on chromosome I and region containing rRNA on chromosome I (the bin 15,050,000-end of chromosome I) are colored in green. (PDF 9816 kb)
Correlation between N2 DpnII experimental data and data from Crane et al. . The 50KB chromatin contact matrix constructed using N2 young adults treated with DpnII restriction enzyme data (GSM2041038- SRR3105476) was compared with 50KB resolution chromatin contacts matrix constructed using the Crane et al. data (GSM1556154 - SRR1665087) from . The total number of paired-ended 37X2 reads in our dataset was 88,466,514, while Crane et al. provide a total of 115,983,178 paired-ended 100X2 reads was. In order to build the chromatin contacts matrix we used the ICE pipeline  iterative mapping implementation from [https://bitbucket.org/mirnylab/hiclib], starting from 21nt up to 37nt in increments of 8. The number of detected Hi-C valid pairs in our dataset was 18,779,498 , consisting of 4,542,078 inter-chromosomal contacts and 14,237,420 intra-chromosomal contacts. In Crane’s dataset the number of valid Hi-C pairs was 59,200,047, consisting of 6,457,271 inter-chromosomal contacts and 52,742,776 intra-chromosomal contacts. Similarly to Additional file 1: Figure S1, any contacts between any location on chromosome I and region containing rRNA on chromosome I (the bin 15,050,000-end of chromosome I) are colored in green. Contacts between genomic loci in adjacent regions (up to 100KB apart are colored in purple). Versions of software used for analysis are as follows: Bowtie2-2.2.6 , and mirnylib/hiclib [https://bitbucket.org/mirnylab/hiclib] downloaded on December 1, 2015. Slight differences in aligned read counts from Tables 2 and 3 reflect updates in alignment software in the concerted package compared to the legacy versions used in Tables 2 and 3. (PDF 1132 kb)
Comparison between RTCC experimental data (this work; black lines), Hi-C data from Crane et al. (2015; magenta lines); and relative representation in anti-LEM2 ChIP-chip data . The curves were obtained by running the ICE pipeline  on our N2 dataset (N2 DpnII GSM2041038- SRR3105476) and on Crane et al. dataset (GSM1556154 - SRR1665087), as implemented in [https://bitbucket.org/mirnylab/hiclib], downloaded on Dec 1, 2015. To obtain the first Eigen Vector values, representing compartments along the chromosome axis, we have followed the tutorial from [https://bitbucket.org/mirnylab/hiclib], using the binnedData class function doEig(numPCs = 1). To inspect the correlation to LEM-2 binding compartments we added LEM-2 binding data  (MA2C normalized log2 ratio of ChIP signal over control), lifted from the ce4 genome assembly to the ce10 assembly, and averaged in 50KB bins. (PDF 46 kb)
About this article
Cite this article
Gabdank, I., Ramakrishnan, S., Villeneuve, A.M. et al. A streamlined tethered chromosome conformation capture protocol. BMC Genomics 17, 274 (2016) doi:10.1186/s12864-016-2596-3
- Caenorhabditis elegans