Skip to main content

Insights about genome function from spatial organization of the genome

Abstract

Over the last 15 years, development of chromosome conformation capture (3C) and its subsequent high-throughput variants in conjunction with the fast development of sequencing technology has allowed investigators to generate large volumes of data giving insights into the spatial three-dimensional (3D) architecture of the genome. This huge data has been analyzed and validated using various statistical, mathematical, genomics, and biophysical tools in order to examine the chromosomal interaction patterns, understand the organization of the chromosome, and find out functional implications of the interactions. This review summarizes the data generated by several large-scale high-throughput chromosome conformation capture studies and the functional implications obtained from the data analyses. We also discuss emerging results on factors (both CCCTC binding factor (CTCF) related and CTCF independent) that could contribute to looping interactions.

Background

Recent advances in the field of chromosome conformation capture (3C) and its subsequent advancements like 4C (chromosome conformation capture-on-chip), 5C (carbon copy chromosome conformation capture), and high-throughput chromosome conformation capture (Hi-C) have greatly expanded the general understanding of the three-dimensional spatial arrangement of the genome. Consequently, the notion of a linear chromosomal distribution of various elements of the coding and the non-coding genome is being revisited. This has brought forth a new perspective on how genomic functions such as replication, transcription, and DNA repair would be investigated. Interestingly, the proximity of various regions of the genome, which were otherwise unintuitive, now opens new possibilities by which various aspects of gene regulation would be studied in future. Along with these, herein, we also discuss recent data suggesting the role of telomeres in looping interactions, both near and with interstitial regions of the genome. In addition, factors other than the chromatin-associated protein CCCTC binding factor (CTCF) that could mechanistically contribute to the formation of chromatin loops are also discussed.

The spatial organization of the genome from Hi-C data

All to all chromosome interaction matrices reveal compartments within genomes

Hi-C is a high-throughput NGS (next-generation sequencing)-based version of chromosome conformation capture assay commonly known as 3C [1]. It provides a global view of all chromosomal interactions across the genome. In such experiments, tagged nucleotides are used to capture ligated products after restriction digestion and appropriately diluted ligation. This is followed by sequencing of the captured DNA fragments. The resulting reads represent ligation of fragments both adjacent to each other as well as far apart in the linear genome. A genome-wide two-dimensional contact matrix is generated by dividing the genome or each chromosome into bins of equal distance and then aligning these reads into those bins. Such matrices depict the collective average of the interactions present within the population of cells used for Hi-C studies [2]. There are several other 3C-based assays to look at chromosomal interactions, but the advantage of Hi-C experiments is that they give a global view of all the chromosomal interactions across the genome.

Though we obtain vast knowledge on the biology of chromosomal interactions and genome organization from Hi-C data, there are some drawbacks of the method as well. The Hi-C data, for instance, does not provide much quantitative information on the stability and strength of chromosomal interactions. Most of the reported Hi-C studies so far used asynchronized cells and did not consider the variability that the different cell cycle stages (particularly mitotic phase) would most likely introduce into the chromatin conformation. Several reviews have further described the technical aspects of these 3C-based assays and complimentary microscopy techniques [3,4,5].

In 2009, Lieberman-Aiden et al. used karyotypically normal GM06990 human lymphoblastoid cells to generate 8.4 million reads that uniquely aligned to the human genome reference sequence in a Hi-C experiment; of these, 6.7 million reads corresponded to long-range contacts between segments of the genome > 20 Kb apart. The long-range contact reads were used to construct genome-wide contact matrices by dividing the genome into 1-Mb-sized bins. This revealed an interesting pattern comprising distinct intra- and inter-chromosomal genomic compartments with contacts that were primarily restricted within compartments [2].

A few years later, using Hi-C contact matrices of bin sizes ranging from 20 to 100 Kb, about 1.7 billion read-pairs of data were analyzed and genome-wide formation of compartments or topologically associated domains (TADs) in mouse and human-differentiated and embryonic stem (ES) cells was reported [6]. This was supported by work from Jin et al. who studied chromatin interactions in response to transient TNF-alpha signaling in IMR90 primary human fibroblasts by performing Hi-C before and after 1 h TNF-alpha treatment [7]. This gave about 3.4 billion uniquely mapped paired-end reads from multiple biological replicates, among which, approximately 1.4 billion were intra-chromosomal reads.

A slightly different approach termed as in situ Hi-C was used by Rao et al. where DNA-DNA proximity ligation was carried out inside intact nuclei—in situ ligation reduced the chances of spurious contacts due to random ligation in dilute solution [8]. This also yielded a higher resolution of up to 1 Kb bin size and 200- to 1000-fold more contacts. Using the high-resolution contact matrices, it was possible to see many smaller contact domains (sub-TADs) with higher intra-domain contact frequency like the topologically associated domains and other similar small domains discussed in several Hi-C studies (Fig. 1) [6, 9, 10].

Fig. 1
figure 1

Schematic representation of how smaller TAD-like structures (sub-TADs) emerge from Hi-C contact matrix and how they might be forming in three dimensions. Dynamic loop extrusion by factors like Cohesin might lead to looping and higher interactions in the extruded loci leading to emergence of sub-TADs

Contacts prevalent within chromosomes rather than across chromosomes

Chromosomal interactions across the whole genome noted that the average intra-chromosomal contact probability between pairs of loci in a chromosome decreased consistently with increasing genomic distance (Fig. 2) [2]. Interestingly, this suggested polymer-like behavior of the genome where the three-dimensional distance between pairs of loci increases with increasing genomic distance. Other studies on chromosomal organization via 3C and fluorescence in situ hybridization (FISH) also made similar observations [1, 11]. Even at distances greater than 200 Mb, the average intra-chromosomal contact probability was always more than the average contact probability between different chromosomes (inter-chromosomal), implying the existence of chromosomal territories.

Fig. 2
figure 2

Frequency distribution of interactions (above cutoff of 10 units of interaction value for a pair of loci) with distance between the interacting loci by analyzing intra-chromosomal Hi-C contact matrix of human chromosome 5 (analysis was done using normalized interaction values from Hi-C data given in Rao et al. [8])

Intriguingly, probabilities of inter-chromosomal contacts show that small, gene-rich chromosomes (chromosomes 16, 17, 19, 20, 21, and 22) preferentially interact with each other. Furthermore, FISH studies also observed these chromosomes to frequently localize to the center of the nucleus [12, 13]. On the other hand, another small but gene-poor chromosome, i.e., 18, was found to interact less with other chromosomes, and FISH studies showed chromosome 18 tends to be located near the nuclear periphery [14].

Genome topologies constructed from interaction matrices

Using principal component analysis, each chromosome could be partitioned into two compartments (termed A and B) such that loci within the same compartment had correlated contact profiles. Through this correlation, even loci belonging to different chromosomes were assigned the same compartment, resulting in the whole genome being divided into two spatial compartments such that greater interactions occurred within than across compartments [15].

It was observed that loci in compartment B had a higher tendency for close spatial localization, suggesting a relatively compact state of chromatin. On the other hand, loci within compartment A showed significant correlation with the presence of genes, higher expression, and chromatin accessibility (DNase I hypersensitivity) and were enriched for activating chromatin marks (H3K36me3). Thus, compartment A could be associated with open, accessible, and actively transcribed chromatin. Together, these suggested open and closed chromatin domains throughout the genome occupy different spatial compartments in the nucleus.

At higher resolution (bin sizes less than 100 Kb), highly self-interacting regions were found to emerge, seen as triangles in the interaction matrix heat map [6]. These regions were termed as topological domains. Topological domains were found to be bound by narrow segments where the chromatin interactions appear to end abruptly. Using a statistic termed directionality index (DI), authors identified 2200 topological domains in mouse ES cells with a median size of 880 Kb that covered ~ 91% of the genome. DI models the difference between a number of upstream and downstream interactions at a given locus along a chromosome, and thereby boundaries of topological domains were detected where there was a significant shift from contact points oriented with upstream vis-à-vis downstream bias. Also, as expected, the frequency of intra-domain interactions was noted to be higher than that of inter-domain interactions.

Boundary demarcations between distinct genome topologies

Consistent with this, fluorescent in situ hybridization (FISH) experiments revealed pairs of loci within a particular topological domain were closer in space than pairs of loci in different topological domains in spite of similar genomic distances between the loci [16]. The genomic regions between the topological domains were defined as either “topological boundary regions” or “unorganized chromatin,” depending on their sizes (topological boundary regions: median ~ 0 Kb, 76.3% of the regions < 50 Kb; or unorganized chromatin: median ~ 560 Kb). Moreover, the topological domains correlated with other described components of the genome-like compartments A and B [2], replication time zones [17, 18], and large organized chromatin K9 modification (LOCK) domains [19]. A large subset of the identified domain boundaries also appeared to mark the transition between LAD (lamina-associated domain) and non-LAD regions in the genome [20, 21].

An in situ Hi-C study using high-resolution contact matrices revealed numerous relatively small contact domains with higher intra-domain contact frequency [8] like topologically associated domains or other small domains discussed in other Hi-C studies [6, 9, 10].

Interestingly, the in situ Hi-C data revealed six nuclear sub-compartments based on long-range interaction patterns, both intra-chromosomal and inter-chromosomal, using different approaches to clustering. On comparison with compartment A/B [2], two of the six interaction patterns correlated with loci in compartment A—termed as sub-compartments A1 and A2. These loci were gene dense, harbored highly expressed genes, enriched in activating chromatin marks, and were depleted at the nuclear lamina- and nucleolus-associated domains (NADs). Rest of interaction patterns correlated with loci in compartment B with very different properties than A1 and A2.

The DI data used to identify TADs can vary substantially depending on the sliding window size selected with small window sizes giving smaller TADs and larger ones yielding larger TADs which often nest groups of smaller domains [22]. In fact, reanalysis of the original Hi-C data from which megabase-sized TADs were identified [6] with a different algorithm using smaller window sizes led to the observation that larger conserved TADs tend to consist entirely of smaller domains. These domains were found to be stable across cell lines and persistent across resolutions, with their boundaries having high enrichment in CTCF binding and activating histone marks [22].

Genome compartments: fractal globule versus topologically associated domain architecture

Chromosomal regions have been perceived as an “equilibrium globule”—a compact and densely knotted configuration originally used to describe a polymer in a poor solvent at equilibrium [23, 24]. An alternative model proposes that polymers, including interphase DNA, can self-organize into a long-lived, non-equilibrium conformation described as “fractal globule” [25, 26]. This dense, compact state is adopted by an untangled polymer as it crumbles into a series of small globules in a “beads-on-a-string” configuration. These beads act as monomers in further rounds of spontaneous crumpling until only a single globule of globules of globules remain.

Lieberman-Aiden et al. analyzed the scaling of contact probability of fractal globule and found that it is close to the contact probability observed from the Hi-C data [2]. The predicted scaling of three-dimensional (3D) distance between pairs of loci based on the fractal globule model is also close to the scaling reported by 3D FISH for genomic distances between 500 Kb and 2 Mb [24]. At a scale of several megabases, the data is consistent with a fractal globule model for chromatin organization. Fractal globule is an attractive model to define chromatin organization since they are free of knots [27] and in principle consistent with unfolding and refolding, for instance, during gene regulatory events like activation and repression or processes such as replication and recombination.

Functional implications of domain formation in the genome

Domain boundaries associated with gene promoters and transcription

A strong enrichment of CTCF binding sites was observed at the TAD boundary regions [6], a property also common with many known insulator or barrier elements [20, 28, 29]. Like a classical boundary element is known to stop the spread of heterochromatin, a clear segregation of the heterochromatin marker H3K9me3 modification was observed within the TAD boundary regions (Fig. 3).

Fig. 3
figure 3

Schematic representation to show that TAD boundary restricts the spread of heterochromatinization

Several studies found TADs to be majorly conserved across cell types in a given organism, suggesting TADs to be stable structures of the 3D genome organization [6, 9, 30]. On the other hand, the smaller scale structures, like sub-TADs, loops, and insulation neighborhoods, all show at least partial variation between different cell lineages, with the variations in their organization appearing to be related with cell type [8, 31,32,33]. The dynamic chromatin interactions varying across cell types were also enriched for differentially expressed genes [6, 34, 35].

These observations suggest that chromatin organization as TADs is mostly stable across cell types, within which specific structures and dynamic interactions can form to play lineage and context-specific regulatory roles contributing to molecular events associated with differentiation [36].

Interestingly, the topological domains do not seem to be the consequence of heterochromatin formation as the detected boundaries were present in both pluripotent cells and their differentiated progeny, i.e., before and after heterochromatinization associated with cellular differentiation. Importantly, this implied that the topological domains along with boundaries delineate the endpoints of heterochromatic spreading [6].

Furthermore, enrichment of chromatin marks associated with promoters and gene bodies and depletion of repressive chromatin marks were detected in the TAD boundaries along with enrichment of housekeeping genes, transcription start sites, and global run on sequencing signal in the boundaries. Together, these suggest a high level of transcriptional activity associated with boundary formation; however, whether boundary formation was a cause or consequence of transcriptional activity was not clear [6].

Point-to-point direct interactions—looping of chromatin affects gene transcription

Interestingly, Hi-C data also revealed pairs of loci that had significantly stronger interaction than any loci lying between them. These were designated as loops, and can appear within topology domains (discussed above), independently and/or across TADs. Interestingly, chromatin loops were not only conserved among human cell lines but also found to be conserved between mouse and human cells [8]; chromatin looping interactions were significantly enriched within cis-regulatory elements like active promoters and enhancers while being depleted at inactive TSS or regions with repressive chromatin marks [6,7,8]. About 30% of chromatin loops brought promoters and enhancers together (versus 7% expected by chance), and genes with promoters associated with loops were more expressed than genes whose promoters were not associated with any loop (sixfold) [8]. Hi-C data analysis also revealed 55% of distal enhancers interact with at least one active promoter, confirming previous observations that promoters and enhancers often form complex networks to regulate transcription [7, 37]. Interestingly, a particular case study showed many genes without any NF-kappaB (p65) binding site in promoters were induced simultaneously by TNF-alpha (which is known to trigger NF-kappaB signaling), possibly due to sharing of overlapping distal interacting regions containing multiple NF-kappaB binding sites [7].

Somewhat intriguingly, there was little or no change in promoter-enhancer looping interactions at a vast majority of TNF-alpha-responsive enhancers on TNF-alpha treatment. This suggested that in general, promoter-enhancer contacts in untreated cells, which are the existing DNA loops, did not alter upon transient activation or repression of enhancers following treatment. Interestingly, chromatin interactions involving cell type-specific enhancers are variable between cell types indicating context-specific interaction structures. This discrepancy between signal-dependent and cell type-specific enhancers correlates with H3K4me1 chromatin marks, which unlike H3K27me3 remain unchanged upon TNF-alpha treatment [7]. Other studies have also observed pre-existing looping interactions at several loci induced by p53, FOXO3, and glucocorticoid receptors [33, 38, 39].

Chromatin loops marked through CTCF binding

Analysis of ENCODE ChIP-seq data revealed loci with chromatin loops were typically bound by the insulator protein CTCF (86%) and the Cohesin subunits RAD21 (86%) and SMC3 (87%) [8]. This was consistent with many reports which, using a variety of experimental approaches, suggest a role for CTCF and cohesion in mediating DNA loops [28, 40, 41]. As many of these loops demarcate domains, this observation was also concordant with studies showing CTCF delimits structural and regulatory domains [6, 42, 43]. Furthermore, most peak loci possessed a unique DNA site containing a CTCF binding motif to which all the three proteins (CTCF, Rad21, and SMC3) bind. A vast majority of these motif pairs present in the peak loci were oriented in a convergent manner suggesting that a pair of CTCF motifs in the convergent orientation might be required for the formation of a loop [8].

Looping through telomere ends

Studies by Shay and Wright’s groups have demonstrated that telomeres loop to specific loci (within 10 Mb)—a phenomenon called TPE-OLD (telomere positioning effect—over long distances). These studies showed that genes close to the telomere were silenced in young primary cells with long telomeres, but were activated when telomeres became short with cellular aging, an effect that was reversed by re-elongation of telomeres upon exogenous expression of the hTERT (human telomerase) gene. Interactions of the telomere and sub-telomere with chromosomal regions up to 10 Mb away from the telomere end were revealed using modified Hi-C and 3C experiments [44,45,46,47]. An important function of such looping showed looping of chromosome 5 sub-telomere resulted in heterochromatic silencing of the telomerase gene in young primary cells [44]. Reports also suggest that several telomere binding proteins like TRF2 and TIN2 can associate with interstitial telomere-like repeats [48]. Such regions of repetitive DNA (often referred to as interstitial telomeric repeats or ITS) may be crucial in forming sub-telomeric loops by recruitment of telomeric factors and structural proteins like Lamins to mediate interaction with telomeres [49, 50].

Single cell Hi-C

Though Hi-C has given insights into the functional chromosomal organization, one can argue that this is an average, probabilistic view of chromosomal interactions with much cell-to-cell variation, such that observed domain organization and chromosomal interactions might represent just a fraction of the cells. Some recent studies attempt to address this with Hi-C at the single cell level. In the first such report, pooled data recapitulated the formation of TAD-like structures indicating these domains are robust and form the basis of chromosome conformation in each cell [51]. However, the observed variability in inter-domain contacts suggested significant differences might be possible in the higher order folding of chromosomes. Another single cell study concluded that, though structures of individual TADs and chromosome loops vary substantially from cell to cell, the higher order organizational signatures like the A/B compartments are mostly retained [15, 52]. In addition, it was also noted that LADs and active enhancers and promoters are consistently organized genome-wide in every cell [52].

Other factors that nucleate chromatin interactions in the 3D genome

The finding that a large fraction of loop boundaries are bound by CTCF and the Cohesin subunits [8] led to popularization of the extrusion model of loop formation. In this model, a pair of factors (Cohesin), possibly with motor function, can dynamically bind to DNA and move along the DNA in opposite directions extruding the chromatin to form loops in between until they dissociate or stall at a boundary element (CTCF motif) [53,54,55]. This model largely helps to understand the nested nature of TADs and loops as well as the consequences of CTCF motif deletion and inversion [36, 56]. However, many convergent CTCF motifs exist that do not delineate loops and, interestingly, a considerable fraction of identified loops did not appear to possess any CTCF or cohesion binding sites. Together, these suggest other factors that could be involved in loop formation.

The fact that CTCF motifs can act as boundary elements only when they are oriented in convergent fashion indicate dimerization of CTCF is required for stalling extruding cohesion loops. This indicates that interactions between DNA binding proteins can possibly play a role in mediating loop formation and chromatin interactions by bringing distant genomic loci together. Indeed, DNA binding proteins, YY1 and ZNF143, were found to be enriched in loop loci [8]. Both homodimer and heterodimer formation by proteins bound to distant loci can lead to chromosomal looping interactions. In a somewhat similar context, the telomere binding factor TRF2 has already been implicated in mediating telomeric looping into the extra-telomeric interstitial regions including the TERT loci [44]. The fact that TRF2 can bind to many extra-telomeric sites throughout the human genome implicates possibility of it mediating other looping interactions as well [49, 50, 57, 58].

Apart from proteins, nucleic acids like lncRNAs and secondary DNA structures might also contribute to mediating chromosomal looping interactions. lncRNAs called as activating/enhancer RNAs were found to bind with mediator proteins (MED12 and MED1), and their knockdown led to diminished looping between the ncRNA loci and their target genes and a decrease in mediator binding to the genes [59, 60].

Another interesting context that could support loop formation in a CTCF-independent way comes from the possible involvement of DNA secondary structures such as G quadruplexes. G quadruplexes are typically formed by a stretch of DNA with 3-guanines repeated at least four times at close interval and are reported to be involved in various biological functions across life forms [61,62,63,64]. The formation of inter-molecular G quadruplexes raises the possibility that such secondary structure DNA motifs can bring together two distant genomic loci. G quadruplexes have been found to be enriched in DNAse hypersensitive (DHS) promoters as well as DHS cis-regulatory elements [65, 66]. Half G quadruplexes (i.e., two runs of guanines, both containing at least three consecutive Gs) were also found to be enriched at the boundaries of these DHS promoters and cis-elements but were depleted in the vicinity of these sites. Computational analyses showed such half G quadruplexes, one from DHS promoter and one from DHS enhancer, could come together to form G quadruplexes thereby bringing together promoters and enhancers elements via looping [65]. Several studies have reported enrichment of potential G quadruplex forming sequences in fragile genomic regions associated with pathogenic, cancer-causing breakpoints, and structural variations [67,68,69,70]. The possible formation of chromosomal loops mediated by G quadruplexes could be one of the contributing factors in increasing the fragility of such chromosomal regions. Disruption of these DNA structure motifs could give a clear idea regarding their roles in mediating functional chromosomal interactions.

Disruption of TAD boundaries has deleterious effects: clinical implications

Deletion of a TAD boundary, on the other hand, was found to result in spreading of contacts across the deleted region and transcriptional misregulation [9]. Several reports of TAD disruptions leading to pathological outcomes have further consolidated TADs as functional units. Structural variants disrupting CTCF-associated TAD boundaries were noted to allow de novo promoter-enhancer interactions and ectopic gene expression causing limb malformation [71] and oncogenic transformation [72]. Oncogenic chromosomal rearrangements were also shown to induce aberrant oncogenic expression in AML and medulloblastoma by causing enhancer trafficking where an enhancer acts on a gene other than its normal target due to genomic rearrangements like TAD disruptions [73, 74]. Moreover, it was also reported from Hi-C in prostate cancer cells that generation of smaller TADs due to establishment of additional boundaries and particular cancer-specific interactions within TADs could be associated with oncogenic transformation [75], and perturbation of CTCF binding by aberrant DNA methylation could cause oncogenic gene expression leading to gliomas by a mechanism similar to “enhancer hijacking” [76]. These pathogenic outcomes based on deregulated enhancer function along with another reporter construct insertion study [77] implicates that cis-regulatory elements can act in a non-specific manner but within a given TAD in the genome (Fig. 4).

Fig. 4
figure 4

Schematic representation showing coordinated expression of genes mediated by an enhancer interacting with multiple promoters within the same TAD and a TAD boundary restricts enhancer interaction and activity only to target genes within the same TAD

Apart from how abrogation of TADs could lead to pathogenic phenotypes, chromatin looping interactions, even when not associated with a particular TAD, might have significant clinical implications. Recently, an interaction has been reported that might be crucial for TERT reactivation across different cancer types. The two cytosine to thymidine single-point mutations in the TERT proximal promoter (− 146 and − 124 bp from the translation start site), recurrent in several cancers, were found to create a de novo consensus binding motif for ETS transcription factor family [78,79,80,81,82,83]. GABPA (an ETS transcription factor) was shown to bind to the mutant promoter leading to a long-range chromatin interaction with a region 300 Kb upstream of the TERT promoter. As an effect of this binding, the promoter locus was changed into an open, active chromatin region, eventually enhancing TERT expression which has been widely associated with cancer development [84,85,86]. In another study, it was shown that telomere looping at the 4q35 locus regulates expression of SORBS2, which is disrupted in the age-associated genetic disease facioscapulohumeral muscular dystrophy [45].

Conclusion

Many large-scale studies with vast Hi-C data have given us important insights into possible mechanisms of looping, subsequent higher order chromatin organization, and functional significance of domain formation and chromosomal interactions in different processes including transcriptional regulation. Given that these observations encourage a shift in perspective regarding how 3D organization of the genome impacts biological processes, it is of much interest to understand the underlying mechanistic rules governing causation of the folded architecture. Though cohesion and CTCF have been implicated in the formation of chromosomal loops, these proteins are absent in a fraction of loop loci (~ 10–14%) suggesting a role of other factors in determining how chromosomal interactions arise and how these contribute to specific cellular and context-specific functions.

Abbreviations

3C:

Chromosome conformation capture

3D:

Three-dimensional

4C:

Chromosome conformation capture-o-chip

5C:

Carbon copy chromosome conformation capture

AML:

Acute myeloid leukemia

ChIP-seq:

Chromatin immunoprecipitation-sequencing

CTCF:

CCCTC binding factor

DI:

Directionality index

ENCODE:

Encyclopedia of DNA elements

ES:

Embryonic stem cells

FISH:

Fluorescence in situ hybridization

hTERT:

Human telomerase reverse transcriptase

LAD:

Lamina-associated domains

lncRNA:

Long non-coding RNA

LOCK:

Large organized chromatin K9 modifications

NAD:

Nucleolus-associated domains

NGS:

Next-generation sequencing

TAD:

Topologically associated domain

TIN2:

TERF1-interacting nuclear factor 2

TNF:

Tumor necrosis factor

TPE-OLD:

Telomere position effect-over long distances

TRF2:

Telomere repeat binding factor 2

TSS:

Transcription start site

References

  1. Dekker J. Capturing chromosome conformation. Science. 2002;295:1306–11.

    Article  CAS  PubMed  Google Scholar 

  2. Lieberman-Aiden E, van Berkum NL, Williams L, Imakaev M, Ragoczy T, Telling A, et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009;326:289–93.

  3. Simonis M, Kooren J, de Laat W. An evaluation of 3C-based methods to capture DNA interactions. Nat Methods. 2007;4:895–901.

    Article  CAS  PubMed  Google Scholar 

  4. De Wit E, De Laat W. A decade of 3C technologies-insights into nuclear organization. Genes Dev. 2012;26:11–24.

  5. Fraser J, Williamson I, Bickmore WA, Dostie J. An overview of genome organization and how we got there: from FISH to Hi-C. Microbiol Mol Biol Rev. 2015;79:347–72.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, et al. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012;485:376–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Jin F, Li Y, Dixon JR, Selvaraj S, Ye Z, Lee AY, et al. A high-resolution map of the three-dimensional chromatin interactome in human cells. Nature. 2013;503:290–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Rao SSP, Huntley MH, Durand NC, Stamenova EK, Bochkov ID, Robinson JT, et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell. 2014;159:1665–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Nora EP, Lajoie BR, Schulz EG, Giorgetti L, Okamoto I, Servant N, et al. Spatial partitioning of the regulatory landscape of the X-inactivation centre. Nature. 2012;485:381–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Sexton T. Three-dimensional folding and functional organization principles of the Drosophila genome. Cell. 2012;148:458–72.

    Article  CAS  PubMed  Google Scholar 

  11. Yokota H, van den Engh G, Hearst JE, Sachs RK, Trask BJ. Evidence for the organization of chromatin in megabase pair-sized loops arranged along a random walk path in the human G0/G1 interphase nucleus. J Cell Biol. 1995;130:1239–49.

    Article  CAS  PubMed  Google Scholar 

  12. Tanabe H, Habermann FA, Solovei I, Cremer M, Cremer T. Non-random radial arrangements of interphase chromosome territories: evolutionary considerations and functional implications. Mutat Res - Fundam Mol Mech Mutagen. 2002;504:37–45.

    Article  CAS  Google Scholar 

  13. Boyle S, Gilchrist S, Bridger JM, Mahy NL, Ellis J a, Bickmore W a. The spatial organization of human chromosomes within the nuclei of normal and emerin-mutant cells. Hum Mol Genet. 2001;10:211–9.

    Article  CAS  PubMed  Google Scholar 

  14. Croft JA, Bridger JM, Boyle S, Perry P, Teague P, Bickmore WA. Differences in the localization and morphology of chromosomes in the human nucleus. J Cell Biol. 1999;145:1119–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Lieberman-Aiden E. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009;326:289–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Eskeland R. Ring1B compacts chromatin structure and represses gene expression independent of histone ubiquitination. Mol Cell. 2010;38:452–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Hiratani I. Genome-wide dynamics of replication timing revealed by in vitro models of mouse embryogenesis. Genome Res. 2010;20:155–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Ryba T. Evolutionarily conserved replication timing profiles predict long-range chromatin interactions and distinguish closely related cell types. Genome Res. 2010;20:761–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Wen B, Wu H, Shinkai Y, Irizarry RA, Feinberg AP. Large histone H3 lysine 9 dimethylated chromatin blocks distinguish differentiated from embryonic stem cells. Nat Genet. 2009;41:246–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Guelen L. Domain organization of human chromosomes revealed by mapping of nuclear lamina interactions. Nature. 2008;453:948–51.

    Article  CAS  PubMed  Google Scholar 

  21. Peric-Hupkes D. Molecular maps of the reorganization of genome-nuclear lamina interactions during differentiation. Mol Cell. 2010;38:603–13.

    Article  CAS  PubMed  Google Scholar 

  22. Filippova D, Patro R, Duggal G, Kingsford C. Identification of alternative topological domains in chromatin. Algorithms Mol Biol. 2014;9:14.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Munkel C, Langowski J. Chromosome structure predicted by a polymer model. Phys Rev E. 1998;57:5888–96.

    Article  CAS  Google Scholar 

  24. Mateos-Langerak J, Bohn M, de Leeuw W, Giromus O, Manders EMM, Verschure PJ, et al. Spatially confined folding of chromatin in the interphase nucleus. Proc Natl Acad Sci U S A. 2009;106:3812–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Grosberg AY, Nechaev SK, Shakhnovich EI. The role of topological constraints in the kinetics of collapse of macromolecules. J Phys Fr. 1988;49:2095–100.

    Article  CAS  Google Scholar 

  26. Grosberg a, Rabin Y, Havlin S, Neer a. Crumpled globule model of the three-dimensional structure of DNA. Europhys Lett. 1993;23:373–8.

    Article  CAS  Google Scholar 

  27. Vasilyev OA, Nechaev SK. Topological correlations in trivial knots: new arguments in favor of the representation of a crumpled polymer globule. Theor Math Phys. 2003;134:142–59.

    Article  Google Scholar 

  28. Phillips JE, Corces VG. CTCF: master weaver of the genome. Cell. 2009;137:1194–211.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Handoko L. CTCF-mediated functional chromatin interactome in pluripotent cells. Nat Genet. 2011;43:630–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Dixon JR, Jung I, Selvaraj S, Shen Y, Antosiewicz-Bourget JE, Lee AY, et al. Chromatin architecture reorganization during stem cell differentiation. Nature. 2015;518:331–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Dowen JM, Fan ZP, Hnisz D, Ren G, Abraham BJ, Zhang LN, et al. Control of cell identity genes occurs in insulated neighborhoods in mammalian chromosomes. Cell. 2014;159:374–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Ji X, Dadon DB, Powell BE, Fan ZP, Borges-Rivera D, Shachar S, et al. 3D chromosome regulatory landscape of human pluripotent cells. Cell Stem Cell. 2016;18:262–75.

    Article  CAS  PubMed  Google Scholar 

  33. Phillips-Cremins JE, Sauria MEG, Sanyal A, Gerasimova TI, Lajoie BR, Bell JSK, et al. Architectural protein subclasses shape 3D organization of genomes during lineage commitment. Cell. 2013;153:1281–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Kagey MH. Mediator and cohesin connect gene expression and chromatin architecture. Nature. 2010;467:430–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Noordermeer D. The dynamic architecture of Hox gene clusters. Science. 2011;334:222–5.

    Article  CAS  PubMed  Google Scholar 

  36. Dixon JR, Gorkin DU, Ren B. Chromatin domains: the unit of chromosome organization. Mol Cell. 2016;62:668–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Smallwood A, Ren B. Genome organization and long-range regulation of gene expression by enhancers. Curr Opin Cell Biol. 2013;25:387–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Melo CA. eRNAs are required for p53-dependent enhancer activity and gene transcription. Mol Cell. 2013;49:524–35.

    Article  CAS  PubMed  Google Scholar 

  39. Tan PY. Integration of regulatory networks by NKX3-1 promotes androgen-dependent prostate cancer survival. Mol Cell Biol. 2012;32:399–414.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Splinter E, Heath H, Kooren J, Palstra RJ, Klous P, Grosveld F, et al. CTCF mediates long-range chromatin looping and local histone modification in the beta-globin locus. Genes Dev. 2006;20:2349–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Hou C, Zhao H, Tanimoto K, Dean A. CTCF-dependent enhancer-blocking by alternative chromatin loop formation. Proc Natl Acad Sci U S A. 2008;105:20398–403.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Xie X, Mikkelsen TS, Gnirke A, Lindblad-Toh K, Kellis M, Lander ES. Systematic discovery of regulatory motifs in conserved regions of the human genome, including thousands of CTCF insulator sites. Proc Natl Acad Sci U S A. 2007;104:7145–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Cuddapah S, Jothi R, Schones DE, Roh TY, Cui K, Zhao K. Global analysis of the insulator binding protein CTCF in chromatin barrier regions reveals demarcation of active and repressive domains. Genome Res. 2009;19:24–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Kim W, Ludlow AT, Min J, Robin JD, Stadler G, Mender I, et al. Regulation of the human telomerase gene TERT by telomere position effect—over long distances (TPE-OLD): implications for aging and cancer. PLoS Biol. 2016;14:e2000016.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Robin JD, Ludlow AT, Batten K, Gaillard MC, Stadler G, Magdinier F, et al. SORBS2 transcription is activated by telomere position effect-over long distance upon telomere shortening in muscle cells from patients with facioscapulohumeral dystrophy. Genome Res. 2015;25:1781–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Robin JD, Ludlow AT, Batten K, Magdinier F, Stadler G, Wagner KR, et al. Telomere position effect: regulation of gene expression with progressive telomere shortening over long distances. Genes Dev. 2014;28:2464–76.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Stadler G, Rahimov F, King OD, Chen JC, Robin JD, Wagner KR, et al. Telomere position effect regulates DUX4 in human facioscapulohumeral muscular dystrophy. Nat Struct Mol Biol. 2013;20:671–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Mignon-Ravix C, Depetris D, Delobel B, Croquette M-F, Mattei M-G. A human interstitial telomere associates in vivo with specific TRF2 and TIN2 proteins. Eur J Hum Genet. 2002;10:107.

    Article  CAS  PubMed  Google Scholar 

  49. Wood AM, Danielsen JMR, Lucas CA, Rice EL, Scalzo D, Shimi T, et al. TRF2 and lamin A/C interact to facilitate the functional organization of chromosome ends. Nat Commun. 2014;5:5467.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Wood AM, Laster K, Rice EL, Kosak ST. A beginning of the end: new insights into the functional organization of telomeres. Nucleus. Taylor & Francis. 2015;6:172–8.

    CAS  Google Scholar 

  51. Nagano T, Lubling Y, Stevens TJ, Schoenfelder S, Yaffe E, Dean W, et al. Single-cell Hi-C reveals cell-to-cell variability in chromosome structure. Nature. 2013;502:59–64.

    Article  CAS  PubMed  Google Scholar 

  52. Stevens TJ, Lando D, Basu S, Liam P, Cao Y, Lee SF, et al. 3D structures of individual mammalian genomes studied by single-cell Hi-C. Nature. 2017;544:59–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Alipour E, Marko JF, Marko JF, Mirny L, Alipour E, Marko J, et al. Self-organization of domain structures by DNA-loop-extruding enzymes. Nucleic Acids Res. 2012;40:11202–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Fudenberg G, Imakaev M, Lu C, Goloborodko A, Abdennur N, Mirny LA. Formation of chromosomal domains by loop extrusion. Cell Rep. 2016;15:2038–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Sanborn AL, Rao SSP, Huang S-C, Durand NC, Huntley MH, Jewett AI, et al. Chromatin extrusion explains key features of loop and domain formation in wild-type and engineered genomes. Proc Natl Acad Sci. 2015;112:201518552.

    Article  Google Scholar 

  56. Boyan B, Giacomo C. Organization and function of the 3D genome. Nat Rev Genet. 2016;17:661–78.

    Article  Google Scholar 

  57. Yang D, Xiong Y, Kim H, He Q, Li Y, Chen R, et al. Human telomeric proteins occupy selective interstitial sites. Cell Res. 2011;21:1013–27.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Simonet T, Zaragosi L-E, Philippe C, Lebrigand K, Schouteden C, Augereau A, et al. The human TTAGGG repeat factors 1 and 2 bind to a subset of interstitial telomeric sequences and satellite repeats. Cell Res. 2011;21:1028–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Trimarchi T, Bilal E, Ntziachristos P, Fabbri G, Dalla-Favera R, Tsirigos A, et al. Genome-wide mapping and characterization of notch-regulated long noncoding RNAs in acute leukemia. Cell. 2014;158:593–606.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Lai F, Orom UA, Cesaroni M, Beringer M, Taatjes DJ, Blobel GA, et al. Activating RNAs associate with Mediator to enhance chromatin architecture and transcription. Nature. 2013;494:497–501.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Rawal P, Kummarasetti VBR, Ravindran J, Kumar N, Halder K, Sharma R, et al. Genome-wide prediction of G4 DNA as regulatory motifs: role in Escherichia coli global regulation. Genome Res. 2006;16:644–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Yadav VK, Abraham JK, Mani P, Kulshrestha R, Chowdhury S. QuadBase: genome-wide database of G4 DNA—occurrence and conservation in human, chimpanzee, mouse and rat promoters and 146 microbes. Nucleic Acids Res. 2008;36:381–5.

    Article  Google Scholar 

  63. Lipps HJ, Rhodes D. G-quadruplex structures: in vivo evidence and function. Trends Cell Biol. 2009;19:414–22.

    Article  CAS  PubMed  Google Scholar 

  64. Rhodes D, Lipps HJ. Survey and summary G-quadruplexes and their regulatory roles in biology. Nucleic Acids Res. 2015;43:8627–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Hegyi H. Enhancer-promoter interaction facilitated by transiently forming G-quadruplexes. Sci Rep. 2015;5:1–6.

    Article  Google Scholar 

  66. Thurman RE, Rynes E, Humbert R, Vierstra J, Maurano MT, Haugen E, et al. The accessible chromatin landscape of the human genome. Nature. 2012;489:75–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Bose P, Hermetz KE, Conneely KN, Rudd MK. Tandem repeats and G-rich sequences are enriched at human CNV breakpoints. PLoS One. 2014;9:1–8.

    Google Scholar 

  68. Murat P, Balasubramanian S. Existence and consequences of G-quadruplex structures in DNA. Curr Opin Genet Dev. 2014;25:22–9.

    Article  CAS  PubMed  Google Scholar 

  69. De S, Michor F. DNA secondary structures and epigenetic determinants of cancer genome evolution. Nat Struct Mol Biol. 2011;18:950–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Katapadi VK, Nambiar M, Raghavan SC. Potential G-quadruplex formation at breakpoint regions of chromosomal translocations in cancer may explain their fragility. Genomics. 2012;100:72–80.

    Article  CAS  PubMed  Google Scholar 

  71. Lupiáñez DG, Kraft K, Heinrich V, Krawitz P, Brancati F, Klopocki E, et al. Disruptions of topological chromatin domains cause pathogenic rewiring of gene-enhancer interactions. Cell. 2015;161:1012–25.

    Article  PubMed  PubMed Central  Google Scholar 

  72. Hnisz D, Weintraub AS, Day DS, Valton A-L, Bak RO, Li CH, et al. Activation of proto-oncogenes by disruption of chromosome neighborhoods. Science. 2016;351:1454–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Northcott PA, Lee C, Zichner T, Stütz AM, Erkek S, Kawauchi D, et al. Enhancer hijacking activates GFI1 family oncogenes in medulloblastoma. Nature. 2014;511:428–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Gröschel S, Sanders MA, Hoogenboezem R, De Wit E, Bouwman BAM, Erpelinck C, et al. A single oncogenic enhancer rearrangement causes concomitant EVI1 and GATA2 deregulation in leukemia. Cell. 2014;157:369–81.

    Article  PubMed  Google Scholar 

  75. Taberlay PC, Achinger-Kawecka J, Lun ATL, Fabian A, Bauer DC, Smyth GK, et al. Three-dimensional disorganisation of the cancer genome occurs coincident with long range genetic and epigenetic alterations. Genome Res. 2016;26:719–31.

  76. Flavahan WA, Drier Y, Liau BB, Gillespie SM, Venteicher AS, Stemmer-Rachamimov AO, et al. Insulator dysfunction and oncogene activation in IDH mutant gliomas. Nature. 2015;529:110–4.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Symmons O, Uslu VV, Tsujimura T, Ruf S, Nassari S, Schwarzer W, et al. Functional and topological characteristics of mammalian regulatory domains. Genome Res. 2014;24:390–400.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Horn S, Figl A, Rachakonda PS, Fischer C, Sucker A, Gast A, et al. TERT promoter mutations in familial and sporadic melanoma. Science. 2013;339:959–61.

    Article  CAS  PubMed  Google Scholar 

  79. Borah S, Xi L, Zaug AJ, Powell NM, Dancik GM, Cohen SB, et al. TERT promoter mutations and telomerase reactivation in urothelial cancer. Science. 2015;347:1006–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Huang FW, Hodis E, Xu MJ, Kryukov GV, Chin L, Garraway LA. Highly recurrent TERT promoter mutations in human melanoma. Science. 2013;339:957–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Vinagre J, Almeida A, Pópulo H, Batista R, Lyra J, Pinto V, et al. Frequency of TERT promoter mutations in human cancers. Nat Commun. 2013;4:2185.

  82. Bell RJA, Rube HT, Kreig A, Mancini A, Fouse SD, Nagarajan RP, et al. The transcription factor GABP selectively binds and activates the mutant TERT promoter in cancer. Science. 2015;348:1036–39.

  83. Li Y, Zhou QL, Sun W, Chandrasekharan P, Cheng HS, Ying Z, et al. Non-canonical NF-κB signalling and ETS1/2 cooperatively drive C250T mutant TERT promoter activation. Nat Cell Biol. 2015;17:1327–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Stern JL, Theodorescu D, Vogelstein B, Papadopoulos N, Cech TR. Mutation of the TERTpromoter, switch to active chromatin, and monoallelic TERTexpression in multiple cancers. Genes Dev. 2015;29:2219–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Xi L, Schmidt JC, Zaug AJ, Ascarrunz DR, Cech TR. A novel two-step genome editing strategy with CRISPR-Cas9 provides new insights into telomerase action and TERT gene expression. Genome Biol. 2015;16:1–17.

    Article  Google Scholar 

  86. Akıncılar SC, Khattar E, Boon PLS, Unal B, Fullwood MJ, Tergaonkar V. Long-range chromatin interactions drive mutant TERT promoter activation. Cancer Discov. 2016;6:1276–92.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Research fellowships from CSIR (to SSR and AKM) are acknowledged. SC is a recipient of Wellcome Trust/DBT India Alliance Fellowship [grant number 500127/Z/09/Z].

Funding

No experimental expenses had been incurred in preparation of this manuscript.

Availability of data and materials

Data sharing is not applicable to this article as no datasets were generated during the current study. All data that was used to generate Fig. 1 can be found in Supplementary information section of Rao et al. [8].

Author information

Authors and Affiliations

Authors

Contributions

SSR contributed to the literature review, Hi-C data analysis (for Fig. 2), making of figures, and manuscript writing; AKM contributed to the making of figures and manuscript writing; and SC contributed to the overall conceptualization of the review article and manuscript writing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Shantanu Chowdhury.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Roy, S.S., Mukherjee, A.K. & Chowdhury, S. Insights about genome function from spatial organization of the genome. Hum Genomics 12, 8 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s40246-018-0140-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s40246-018-0140-z

Keywords