Chromosome conformation capture techniques (often abbreviated to 3C technologies or 3C-based methods[1]) are a set of molecular biology methods used to analyze the spatial organization of chromatin in a cell. These methods quantify the number of interactions between genomic loci that are nearby in 3-D space, but may be separated by many nucleotides in the linear genome.[2] Such interactions may result from biological functions, such as promoter-enhancer interactions, or from random polymer looping, where undirected physical motion of chromatin causes loci to collide.[3] Interaction frequencies may be analyzed directly,[4] or they may be converted to distances and used to reconstruct 3-D structures.[5]
The chief difference between 3C-based methods is their scope. For example, when using PCR to detect interaction in a 3C experiment, the interactions between two specific fragments are quantified. In contrast, Hi-C quantifies interactions between all possible pairs of fragments simultaneously. Deep sequencing of material produced by 3C also produces genome-wide interactions maps.
In 1993, the Nuclear Ligation Assay was published, a method that could determine circularization frequencies of DNA in solution. This assay was used to show that estrogen induces an interaction between the prolactin gene promoter and a nearby enhancer.[16]
In 2002, Job Dekker introduced the new idea that dense matrices of interaction frequencies between loci could be used to infer the spatial organization of genomes. This idea was the basis for his development of the chromosome conformation capture (3C) assay, published in 2002 by Job Dekker and colleagues in the Kleckner lab at Harvard University.[17][18]
In 2013, Takashi Nagano and Peter Fraser introduced in-nuclei ligation for Hi-C and single-cell Hi-C.[26]
In 2014, Suhas Rao, Miriam Huntley, et al. developed in-situ Hi-C and the use of 4-cutter restriction enzymes, and released the first high-resolution datasets down to kilobase resolution for several human cell lines. They also identified the first clear evidence of CTCF-Cohesin looping in Hi-C maps and identified the convergent CTCF motif rule underlying these loops.[27]
All 3C methods start with a similar set of steps, performed on a sample of cells.
First, the cell genomes are cross-linked with formaldehyde,[28] which introduces bonds that "freeze" interactions between genomic loci. Treatment of cells with 1-3% formaldehyde, for 10-30min at room temperature is most common, however, standardization for preventing high protein-DNA cross linking is necessary, as this may negatively affect the efficiency of restriction digestion in the subsequent step.[29] The genome is then cut into fragments with a restriction endonuclease. The size of restriction fragments determines the resolution of interaction mapping. Restriction enzymes (REs) that make cuts on 6bp recognition sequences, such as EcoR1 or HindIII, are used for this purpose, as they cut the genome once every 4000bp, giving ~ 1 million fragments in the human genome.[29][30] For more precise interaction mapping, a 4bp recognizing RE may also be used. The next step is, proximity based ligation. This takes place at low DNA concentrations or within intact, permeabilized nuclei[26] in the presence of T4 DNA ligase,[31] such that ligation between cross-linked interacting fragments is favored over ligation between fragments that are not cross-linked. Subsequently, interacting loci are quantified by amplifying ligated junctions by PCR methods.[29][31]
Original methods
3C (one-vs-one)
The chromosome conformation capture (3C) experiment quantifies interactions between a single pair of genomic loci. For example, 3C can be used to test a candidate promoter-enhancer interaction. Ligated fragments are detected using PCR with known primers.[2][17] That is why this technique requires the prior knowledge of the interacting regions.
4C (one-vs-all)
Chromosome conformation capture-on-chip (4C) (also known as circular chromosome conformation capture) captures interactions between one locus and all other genomic loci. It involves a second ligation step, to create self-circularized DNA fragments, which are used to perform inverse PCR. Inverse PCR allows the known sequence to be used to amplify the unknown sequence ligated to it.[32][2][19] In contrast to 3C and 5C, the 4C technique does not require the prior knowledge of both interacting chromosomal regions. Results obtained using 4C are highly reproducible with most of the interactions that are detected between regions proximal to one another. On a single microarray, approximately a million interactions can be analyzed. [citation needed]
5C (many-vs-many)
Chromosome conformation capture carbon copy (5C) detects interactions between all restriction fragments within a given region, with this region's size typically no greater than a megabase.[2][20] This is done by ligating universal primers to all fragments. However, 5C has relatively low coverage. The 5C technique overcomes the junctional problems at the intramolecular ligation step and is useful for constructing complex interactions of specific loci of interest. This approach is unsuitable for conducting genome-wide complex interactions since that will require millions of 5C primers to be used.[citation needed]
Hi-C uses high-throughput sequencing to find the nucleotide sequence of fragments[2][22] and uses paired end sequencing, which retrieves a short sequence from each end of each ligated fragment. As such, for a given ligated fragment, the two sequences obtained should represent two different restriction fragments that were ligated together in the proximity based ligation step. The pair of sequences are individually aligned to the genome, thus determining the fragments involved in that ligation event. Hence, all possible pairwise interactions between fragments are tested.
Sequence capture-based methods
A number of methods use oligonucleotide capture to enrich 3C and Hi-C libraries for specific loci of interest.[33][34] These methods include Capture-C,[35] NG Capture-C,[36] Capture-3C,[35] HiCap,[33][37] Capture Hi-C.[38] and Micro Capture-C.[39] These methods are able to produce higher resolution and sensitivity than 4C based methods,[40] Micro Capture-C provides the highest resolution of the available 3C techniques and it is possible to generate base pair resolution data.[39]
Single-cell methods
Single-cell adaptations of these methods, such as ChIP-seq and Hi-C can be used to investigate the interactions occurring in individual cells.[41][42]
Multi-interaction methods
A number of methods sequence multiple ligation junctions simultaneously to detect higher-order structures where multiple regions of chromatin may be interacting. These methods include Tri-C,[43] 3way 4C/C-walks,[44] and multi-contact 4C (MC-4C).[45]
Immunoprecipitation-based methods
ChIP-loop
ChIP-loop combines 3C with ChIP-seq to detect interactions between two loci of interest mediated by a protein of interest.[2][46] The ChIP-loop may be useful in identifying long-range cis-interactions and trans interaction mediated through proteins since frequent DNA collisions will not occur. [citation needed]
Genome wide methods
ChIA-PET combines Hi-C with ChIP-seq to detect all interactions mediated by a protein of interest.[2][23] HiChIP was designed to allow similar analysis as ChIA-PET with less input material.[47]
Biological impact
3C methods have led to a number of biological insights, including the discovery of new structural features of chromosomes, the cataloguing of chromatin loops, and increased understanding of transcriptional regulation mechanisms (the disruption of which can lead to disease).[6]
3C methods have demonstrated the importance of spatial proximity of regulatory elements to the genes that they regulate. For example, in tissues that express globin genes, the β-globin locus control region forms a loop with these genes. This loop is not found in tissues where the gene is not expressed.[48] This technology has further aided the genetic and epigenetic study of chromosomes both in model organisms and in humans.[not verified in body]
These methods have revealed large-scale organization of the genome into topologically associating domains (TADs), which correlate with epigenetic markers. Some TADs are transcriptionally active, while others are repressed.[49] Many TADs have been found in D. melanogaster, mouse and human.[50] Moreover, CTCF and cohesin play important roles in determining TADs and enhancer-promoter interactions. The result shows that the orientation of CTCF binding motifs in an enhancer-promoter loop should be facing to each other in order for the enhancer to find its correct target.[51]
Human disease
There are several diseases caused by defects in promoter-enhancer interactions, which are reviewed in this paper.[52]
Beta thalassemia is a certain type of blood disorder caused by a deletion of LCR enhancer element.[53][54]
Holoprosencephaly is cephalic disorder caused by a mutation in the SBE2 enhancer element, which in turn weakened the production of SHH gene.[55]
PPD2 (polydactyly of a triphalangeal thumb) is caused by a mutation of ZRS enhancer, which in turn strengthened the production of SHH gene.[56][57]
The different 3C-style experiments produce data with very different structures and statistical properties. As such, specific analysis packages exist for each experiment type.[34]
Hi-C data is often used to analyze genome-wide chromatin organization, such as topologically associating domains (TADs), linearly contiguous regions of the genome that are associated in 3-D space.[49] Several algorithms have been developed to identify TADs from Hi-C data.[4][64]
Hi-C and its subsequent analyses are evolving. Fit-Hi-C [3] is a method based on a discrete binning approach with modifications of adding distance of interaction (initial spline fitting, aka spline-1) and refining the null model (spline-2). The result of Fit-Hi-C is a list of pairwise intra-chromosomal interactions with their p-values and q-values.[63]
The 3-D organization of the genome can also be analyzed via eigendecomposition of the contact matrix. Each eigenvector corresponds to a set of loci, which are not necessarily linearly contiguous, that share structural features.[65]
A significant confounding factor in 3C technologies is the frequent non-specific interactions between genomic loci that occur due to random polymer behavior. An interaction between two loci must be confirmed as specific through statistical significance testing.[3]
Normalization of Hi-C contact map
There are two major ways of normalizing raw Hi-C contact heat maps. The first way is to assume equal visibility, meaning there is an equal chance for each chromosomal position to have an interaction. Therefore, the true signal of a Hi-C contact map should be a balanced matrix (Balanced matrix has constant row sums and column sums). An example of algorithms that assumes equal visibility is Sinkhorn-Knopp algorithm, which scales the raw Hi-C contact map into a balanced matrix.
The other way is to assume there is a bias associated with each chromosomal position. The contact map value at each coordinate will be the true signal at that position times bias associated with the two contact positions. An example of algorithms that aim to solve this model of bias is iterative correction, which iteratively regressed out row and column bias from the raw Hi-C contact map. There are a number of software tools available for analysis of Hi-C data.[66]
DNA motif analysis
DNA motifs are specific short DNA sequences, often 8-20 nucleotides in length[67] which are statistically overrepresented in a set of sequences with a common biological function. Currently, regulatory motifs on the long-range chromatin interactions have not been studied extensively. Several studies have focused on elucidating the impact of DNA motifs in promoter-enhancer interactions.
Bailey et al. has identified that ZNF143 motif in the promoter regions provides sequence specificity for promoter-enhancer interactions.[68] Mutation of ZNF143 motif decreased the frequency of promoter-enhancer interactions suggesting that ZNF143 is a novel chromatin-looping factor.
For genome-scale motif analysis, in 2016, Wong et al. reported a list of 19,491 DNA motif pairs for K562 cell line on the promoter-enhancer interactions.[69] As a result, they proposed that motif pairing multiplicity (number of motifs that are paired with a given motif) is linked to interaction distance and regulatory region type. In the next year, Wong published another article reporting 18,879 motif pairs in 6 human cell lines.[70] A novel contribution of this work is MotifHyades, a motif discovery tool that can be directly applied to paired sequences.
Cancer genome analysis
The 3C-based techniques can provide insights into the chromosomal rearrangements in the cancer genomes.[71] Moreover, they can show changes of spatial proximity for regulatory elements and their target genes, which bring deeper understanding of the structural and functional basis of the genome.[72]
^ abSimonis M, Klous P, Splinter E, Moshkin Y, Willemsen R, de Wit E, van Steensel B, de Laat W (November 2006). "Nuclear organization of active and inactive chromatin domains uncovered by chromosome conformation capture-on-chip (4C)". Nature Genetics. 38 (11): 1348–54. doi:10.1038/ng1896. PMID17033623. S2CID22787572.
^Gavrilov A, Eivazova E, Priozhkova I, Lipinski M, Razin S, Vassetzky Y (2009). "Chromosome Conformation Capture (From 3C to 5C) and Its ChIP-Based Modification". Chromatin Immunoprecipitation Assays. review. Methods in Molecular Biology. Vol. 567. pp. 171–88. doi:10.1007/978-1-60327-414-2_12. ISBN978-1-60327-413-5. PMID19588093.
^Horike S, Cai S, Miyano M, Cheng JF, Kohwi-Shigematsu T (January 2005). "Loss of silent-chromatin looping and impaired imprinting of DLX5 in Rett syndrome". Nature Genetics. 37 (1): 31–40. doi:10.1038/ng1491. PMID15608638. S2CID2884412.
^Krijger PH, de Laat W (December 2016). "Regulation of disease-associated gene expression in the 3D genome". Nature Reviews Molecular Cell Biology. 17 (12): 771–782. doi:10.1038/nrm.2016.138. PMID27826147. S2CID11484886.
^Van der Ploeg LH, Konings A, Oort M, Roos D, Bernini L, Flavell RA (February 1980). "gamma-beta-Thalassaemia studies showing that deletion of the gamma- and delta-genes influences beta-globin gene expression in man". Nature. 283 (5748): 637–42. Bibcode:1980Natur.283..637V. doi:10.1038/283637a0. PMID6153459. S2CID4371542.
^Wieczorek D, Pawlik B, Li Y, Akarsu NA, Caliebe A, May KJ, et al. (January 2010). "A specific mutation in the distant sonic hedgehog (SHH) cis-regulator (ZRS) causes Werner mesomelic syndrome (WMS) while complete ZRS duplications underlie Haas type polysyndactyly and preaxial polydactyly (PPD) with or without triphalangeal thumb". Human Mutation. 31 (1): 81–9. doi:10.1002/humu.21142. PMID19847792. S2CID1715146.
^Bailey, S. D., Zhang, X., Desai, K., Aid, M., Corradin, O., Cowper-Sal·lari, R., … Lupien, M. (2015). ZNF143 provides sequence specificity to secure chromatin interactions at gene promoters. Nature Communications, 2, 6186. Retrieved from https://doi.org/10.1038/ncomms7186
^K. Wong, Y. Li, and C. Peng, “Identification of coupling DNA motif pairs on long-range chromatin interactions in human,” vol. 32, no. September 2015, pp. 321–324, 2016.
^Ka-Chun Wong; MotifHyades: expectation maximization for de novo DNA motif pair discovery on paired sequences, Bioinformatics, Volume 33, Issue 19, 1 October 2017, Pages 3028–3035, https://doi.org/10.1093/bioinformatics/btx381