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.
Historically, microscopy was the primary method of investigating nuclear organization, [6] which can be dated back to 1590. [7]
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]
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.
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 ]
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.
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 adaptations of these methods, such as ChIP-seq and Hi-C can be used to investigate the interactions occurring in individual cells. [41] [42]
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]
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 ]
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]
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]
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]
Adenocarcinoma of the lung can be caused by a duplication of enhancer element for MYC gene. [58]
T-cell acute lymphoblastic leukemia is caused by an introduction of a new enhancer. [59]
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]
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 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.
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]
Chromatin is a complex of DNA and protein found in eukaryotic cells. The primary function is to package long DNA molecules into more compact, denser structures. This prevents the strands from becoming tangled and also plays important roles in reinforcing the DNA during cell division, preventing DNA damage, and regulating gene expression and DNA replication. During mitosis and meiosis, chromatin facilitates proper segregation of the chromosomes in anaphase; the characteristic shapes of chromosomes visible during this stage are the result of DNA being coiled into highly condensed chromatin.
In genetics, an enhancer is a short region of DNA that can be bound by proteins (activators) to increase the likelihood that transcription of a particular gene will occur. These proteins are usually referred to as transcription factors. Enhancers are cis-acting. They can be located up to 1 Mbp away from the gene, upstream or downstream from the start site. There are hundreds of thousands of enhancers in the human genome. They are found in both prokaryotes and eukaryotes. Active enhancers typically get transcribed as enhancer or regulatory non-coding RNA, whose expression levels correlate with mRNA levels of target genes.
A regulatory sequence is a segment of a nucleic acid molecule which is capable of increasing or decreasing the expression of specific genes within an organism. Regulation of gene expression is an essential feature of all living organisms and viruses.
In molecular biology and genetics, transcriptional regulation is the means by which a cell regulates the conversion of DNA to RNA (transcription), thereby orchestrating gene activity. A single gene can be regulated in a range of ways, from altering the number of copies of RNA that are transcribed, to the temporal control of when the gene is transcribed. This control allows the cell or organism to respond to a variety of intra- and extracellular signals and thus mount a response. Some examples of this include producing the mRNA that encode enzymes to adapt to a change in a food source, producing the gene products involved in cell cycle specific activities, and producing the gene products responsible for cellular differentiation in multicellular eukaryotes, as studied in evolutionary developmental biology.
In biology, the epigenome of an organism is the collection of chemical changes to its DNA and histone proteins that affects when, where, and how the DNA is expressed; these changes can be passed down to an organism's offspring via transgenerational epigenetic inheritance. Changes to the epigenome can result in changes to the structure of chromatin and changes to the function of the genome. The human epigenome, including DNA methylation and histone modification, is maintained through cell division. The epigenome is essential for normal development and cellular differentiation, enabling cells with the same genetic code to perform different functions. The human epigenome is dynamic and can be influenced by environmental factors such as diet, stress, and toxins.
An insulator is a type of cis-regulatory element known as a long-range regulatory element. Found in multicellular eukaryotes and working over distances from the promoter element of the target gene, an insulator is typically 300 bp to 2000 bp in length. Insulators contain clustered binding sites for sequence specific DNA-binding proteins and mediate intra- and inter-chromosomal interactions.
Transcriptional repressor CTCF also known as 11-zinc finger protein or CCCTC-binding factor is a transcription factor that in humans is encoded by the CTCF gene. CTCF is involved in many cellular processes, including transcriptional regulation, insulator activity, V(D)J recombination and regulation of chromatin architecture.
Upstream stimulatory factor 1 is a protein that in humans is encoded by the USF1 gene.
DNA adenine methyltransferase identification, often abbreviated DamID, is a molecular biology protocol used to map the binding sites of DNA- and chromatin-binding proteins in eukaryotes. DamID identifies binding sites by expressing the proposed DNA-binding protein as a fusion protein with DNA methyltransferase. Binding of the protein of interest to DNA localizes the methyltransferase in the region of the binding site. Adenine methylation does not occur naturally in eukaryotes and therefore adenine methylation in any region can be concluded to have been caused by the fusion protein, implying the region is located near a binding site. DamID is an alternate method to ChIP-on-chip or ChIP-seq.
Chromatin Interaction Analysis by Paired-End Tag Sequencing is a technique that incorporates chromatin immunoprecipitation (ChIP)-based enrichment, chromatin proximity ligation, Paired-End Tags, and High-throughput sequencing to determine de novo long-range chromatin interactions genome-wide.
Cas9 is a 160 kilodalton protein which plays a vital role in the immunological defense of certain bacteria against DNA viruses and plasmids, and is heavily utilized in genetic engineering applications. Its main function is to cut DNA and thereby alter a cell's genome. The CRISPR-Cas9 genome editing technique was a significant contributor to the Nobel Prize in Chemistry in 2020 being awarded to Emmanuelle Charpentier and Jennifer Doudna.
STARR-seq is a method to assay enhancer activity for millions of candidates from arbitrary sources of DNA. It is used to identify the sequences that act as transcriptional enhancers in a direct, quantitative, and genome-wide manner.
A topologically associating domain (TAD) is a self-interacting genomic region, meaning that DNA sequences within a TAD physically interact with each other more frequently than with sequences outside the TAD. The average size of a topologically associating domain (TAD) is 1000 kb in humans, 880 kb in mouse cells, and 140 kb in fruit flies. Boundaries at both side of these domains are conserved between different mammalian cell types and even across species and are highly enriched with CCCTC-binding factor (CTCF) and cohesin. In addition, some types of genes appear near TAD boundaries more often than would be expected by chance.
Gene deserts are regions of the genome that are devoid of protein-coding genes. Gene deserts constitute an estimated 25% of the entire genome, leading to the recent interest in their true functions. Originally believed to contain inessential and "Junk DNA" due to their inability to create proteins, gene deserts have since been linked to several vital regulatory functions, including distal enhancing and conservatory inheritance. Thus, an increasing number of risks that lead to several major diseases, including a handful of cancers, have been attributed to irregularities found in gene deserts.
Nuclear organization refers to the spatial distribution of chromatin within a cell nucleus. There are many different levels and scales of nuclear organisation. Chromatin is a higher order structure of DNA.
In mammalian biology, insulated neighborhoods are chromosomal loop structures formed by the physical interaction of two DNA loci bound by the transcription factor CTCF and co-occupied by cohesin. Insulated neighborhoods are thought to be structural and functional units of gene control because their integrity is important for normal gene regulation. Current evidence suggests that these structures form the mechanistic underpinnings of higher-order chromosome structures, including topologically associating domains (TADs). Insulated neighborhoods are functionally important in understanding gene regulation in normal cells and dysregulated gene expression in disease.
Single cell epigenomics is the study of epigenomics in individual cells by single cell sequencing. Since 2013, methods have been created including whole-genome single-cell bisulfite sequencing to measure DNA methylation, whole-genome ChIP-sequencing to measure histone modifications, whole-genome ATAC-seq to measure chromatin accessibility and chromosome conformation capture.
Hi-C is a high-throughput genomic and epigenomic technique to capture chromatin conformation (3C). In general, Hi-C is considered as a derivative of a series of chromosome conformation capture technologies, including but not limited to 3C, 4C, and 5C. Hi-C comprehensively detects genome-wide chromatin interactions in the cell nucleus by combining 3C and next-generation sequencing (NGS) approaches and has been considered as a qualitative leap in C-technology development and the beginning of 3D genomics.
Proximity ligation-assisted chromatin immunoprecipitation sequencing (PLAC-seq) is a chromatin conformation capture(3C)-based technique to detect and quantify genomic chromatin structure from a protein-centric approach. PLAC-seq combines in situ Hi-C and chromatin immunoprecipitation (ChIP), which allows for the identification of long-range chromatin interactions at a high resolution with low sequencing costs. Mapping long-range 3-dimensional(3D) chromatin interactions is important in identifying transcription enhancers and non-coding variants that can be linked to human diseases.
Pore-C is a genomic technique which utilizes chromatin conformation capture (3C) and Oxford Nanopore Technologies' (ONT) long-read sequencing to characterize three-dimensional (3D) chromatin structure. To characterize concatemers, the originators of Pore-C developed an algorithm to identify alignments that are assigned to a restriction fragment; concatemers with greater than two associated fragments are deemed high order. Pore-C attempts to improve on previous 3C technologies, such as Hi-C and SPRITE, by not requiring DNA amplification prior to sequencing. This technology was developed as a simpler and more easily scalable method of capturing higher-order chromatin structure and mapping regions of chromatin contact. In addition, Pore-C can be used to visualize epigenomic interactions due to the capability of ONT long-read sequencing to detect DNA methylation. Applications of this technology include analysis of combinatorial chromatin interactions, the generation of de novo chromosome scale assemblies, visualization of regions associated with multi-locus histone bodies, and detection and resolution of structural variants.