Skip to main content

An Overview of Methods for Reconstructing 3-D Chromosome and Genome Structures from Hi-C Data


Over the past decade, methods for predicting three-dimensional (3-D) chromosome and genome structures have proliferated. This has been primarily due to the development of high-throughput, next-generation chromosome conformation capture (3C) technologies, which have provided next-generation sequencing data about chromosome conformations in order to map the 3-D genome structure. The introduction of the Hi-C technique—a variant of the 3C method—has allowed researchers to extract the interaction frequency (IF) for all loci of a genome at high-throughput and at a genome-wide scale. In this review we describe, categorize, and compare the various methods developed to map chromosome and genome structures from 3C data—particularly Hi-C data. We summarize the improvements introduced by these methods, describe the approach used for method evaluation, and discuss how these advancements shape the future of genome structure construction.


After decades of research about the organization of the nucleus of the eukaryotic cell, there exists substantial evidence that the genome architecture plays a key role in nuclear functions. [1,2,3,4,5,6,7,8]. For instance, the spatial arrangement and proximity of genes has been linked to biological functions such as gene replication, regulation and transcription. [6, 9,10,11].

The impact of genome architecture on nuclear processes spans multiple hierarchical levels including the spatial compartmentalization of the process, the higher-order organization of chromatin and the arrangement of the genome within the nucleus. Despite the dynamic nature of their process components, processes such as transcription and DNA repair have been shown to be constrained to specific spatial locations rather than randomly dispersed throughout the nucleus. Genes tend to be more active in sparse euchromatin than dense heterochromatin, purportedly due to the impact of folding density on regulatory factor availability. The homogeneous topology of chromatin has potential to capture nuclear proteins, affecting their probability of interaction with binding sites. Small, kilo-base sized chromatin loops can localize promoters with upstream elements while larger mega-base sized loops can spatially segregate nuclear regions imposing independence on different processes.

Understanding the 3-D organization of the eukaryotic genome is essential to explain the important chromosomal activities within the cell. Hence, a fundamental question in genome and biological studies is how the spatial conformation of the chromosome in the nucleus affects a number of genetic and biological functions such as gene regulation [12, 13], gene expression [14], transcription regulation [15], DNA repair, and DNA replication [16, 17].

Early studies of chromosome conformation relied on the use of cytogenetic techniques. An example of the which is fluorescence in situ hybridization (FISH), employed to detect the presence of a specific chromosome region and the proximity between two regions in a genome sequence [18, 19]. Fluorescence in situ hybridization uses fluorescent probes that bind to specific regions of a chromosome with a high degree of sequence complementarity. Using fluorescence microscopy, the location of the loci or DNA sequence with which a probe is expected to bind may be determined. This method is especially useful, as it allows direct one-to-one estimation of genome loci proximity. However, due to technical limitations such as low-throughput, low resolution of FISH data, and probe requirements for every analysis, it is not optimal for examining multiple positions simultaneously. As a result, the method is not used when studying the organization of chromosomes at a genome-wide scale. Other microscopy techniques that have been developed to study the chromatin organization aimed at providing details about the genome positioning and activities. Some of these methods are called the super resolution microscopy strategies, as they were developed to provide imaging at a high resolution. Examples are saturated structured illumination microscopy (SSIM), stimulated emission depletion (STED), and ground state depletion (GSD) [20, 21]. The introduction of Stochastic super-resolution microscopy techniques such as Photo-activated localization microscopy (PALM or FPALM) and stochastic optical reconstruction microscopy (STORM) produced a different set of ways for investigating the chromatin organization [22, 23]. Generally, the microscopy techniques for studying the chromatin organization could be categorized as light and electron microscopy-based techniques. The more detailed description of the microscopy-based techniques for studying genome organization is given in the section “Genome Organization by microscopy-based techniques”.

In 2002, Dekker et al. [24] developed 3C, a high-throughput methodology that can be used to generate IFs between nearby genomic loci in a cell population. Since then, a number of 3C variants [25,26,27] such as 4C [28], 5C [29], Hi-C [30], TCC [31], ChIA-PET [32, 33] and, later on, single-cell Hi-C [34], have been developed to study the 3-D organization of the chromosome and genome. The development of 3C techniques has substantially benefited the study of the spatial proximity, interaction, and genome conformation of a number of cells. Today, Hi-C is the most widely used and well-known 3C variant. Using next-generation sequencing strategies such as high-throughput and parallel sequencing, Hi-C enables researchers to profile read-pair interactions on an all-versus-all basis—that is, to profile interactions for all read pairs in an entire genome. It also allows them to detect and compute the number of interactions between fragments within a chromosome—i.e., the intra-chromosome interaction frequency (IF)—or between different chromosomes—i.e., the inter-chromosome interaction frequency. Fragments, alternatively known as bins or genomic loci, are the regions to which a chromosome have been divided into. Each fragment has a defined length or size which is the number of base pair (bp) in it. The size of the fragment is determined by the resolution, e.g. a 1 MB resolution signifies that 1,000,000 bp are contained within each fragment.

The IFs obtained are commonly represented in a two-dimensional matrix, also known as a contact matrix, with rows and columns representing the number of fragments in the chromosome or genome.

The Hi-C technique is especially relevant because the IFs it yields can be used to construct 3-D chromosome and genome structures. These structures, in turn, help explain a series of events such as genome folding, gene regulations, the connection between regulatory elements and the higher-order structural features in the nucleus of a cell [1, 2, 14, 35, 36].

Within the past decade, a number of computational methods and algorithms have been proposed for the construction of chromosome and genome 3-D structures from Hi-C data. Most of these methods adopt different strategies for 3-D structure prediction, have different technical requirements for algorithms, and use different noise reduction techniques to analyze Hi-C data. In this review, we categorize these methods based on how they model IF from Hi-C data, highlight a common approach to method evaluation and validation, and finally point to the future direction and challenges of chromosome and genome 3-D structure prediction.

Description of the Hi-C Experiment and Chromosomal Contact Map

Using next generation sequencing technology, the emergence of the Hi-C technique, an extension of 3C, has enabled the identification of the chromosome conformation at a genome wide scale [26, 27, 30, 37, 38]. Compared to other variant of the 3C technique, the Hi-C technique is the first method [30, 38] to capture chromosome conformation on a “all versus all” basis —that is, it can profile interactions for all read pairs in an entire genome. Hi-C protocol begins by using formaldehyde to crosslink the cells, which results in the covalent linking of the chromosomal loci through their protein-DNA interactions. The cross-linked chromatin segment is then cut out with a restriction enzyme, and the segment restriction ends are marked by filing in with biotin-labeled nucleotides [25, 30]. Next, the resulting blunt-end segments are ligated randomly under appropriate condition for ligation events between the cross-linked DNA segments. DNA is purified and sheared, and a biotin pull-down is performed to ensured that only the biotinylated junctions are selected for further high throughput pair-end sequencing and computational analysis. After the sequencing of the pair-reads, the generated output usually in .fastq format is mapped to a reference genome, filtered, and used to create a contact map [39]. Notable tools that support the mapping of the sequenced pair reads to generate contact map are GenomeFlow [40], Juicer [41], HiC-Pro [42], Hi-Cpipe [43], and HiCUP [44].

Interaction frequency, sometimes referred to as contact frequency, is a measure of the number of interactions between a pair of chromosomal or genomic regions in the Hi-C data [45,46,47,48]. The combined contact counts for all pairwise regions or loci may be represented as a symmetric matrix to form an IF matrix of all interacting fragments. The IF matrix is sometimes referred to as a contact matrix or contact map [30, 47]. A chromosome contact matrix is a n-by-n matrix representing the interaction of loci or chromosomal regions as captured in the Hi-C experiment [27, 30, 31, 49]. The rows and columns of the matrix correspond to the index of the equal-sized regions which partition the chromosome. The length of one equal-sized region (e.g., 1 Mb base pair) is referred to as the resolution [30]. Each entry in the matrix represents a count of read pairs that connect two corresponding chromosome regions in a Hi-C experiment [30]. Alternatively, the contacts can be represented in a 3-column sparse matrix [49], where columns 1 and 2 refer to the genomic location or the fragment number of the interacting loci and column 3 represents the IF between them.

Polymer Model

Polymer models are based on the underlying idea that interactions between molecular subunits such as monomers result in large molecular structures known as polymers. This approach was adopted from polymer physics, a branch of statistical physics [50,51,52]. Polymers produced by living organisms are referred to as biopolymers. Two well-known examples of biopolymers are DNA and proteins, with nucleotides and amino acids as their monomers, respectively. Polymerization involves the combination of small molecules through chemical bonding to form a network at equilibrium called a polymer. Various authors have adopted two states of the polymer to model the architecture of chromosomal regions in a cell: the equilibrium globule [53, 54] and the fractal globule [37, 55, 56]. A characteristic feature of the equilibrium globule model is that it is highly knotted [30]. Mirny [37] has pointed out that this configuration is disadvantageous, as it restricts genomic processes such as unfolding—an important property for gene activation—or refolding [57]. Alternatively, Barbieri et al. [55] showed that polymer collapse after exposure to a topological constraint can result in the formation of a long-lived, untangled, non-equilibrium configuration state called a crumpled or fractal globule. A fractal globule is knot-free, and it is organized such that it allows for unfolding or refolding processes while in a highly compact state. Hence, the polymer exhibits a “beads-on-a-string” configuration, with beads representing monomers connected by linkers; DNA connections in eukaryotic chromatin are similarly configured. The fractal globule can be illustrated as a dense multicolor ball of yarn, where each color has its own end, but one can pull out threads with a specific color and put them back in, without disturbing the structure of the overall ball at all. This important property makes the fractal globule suitable for organizing chromatin in a cell because this topology facilitates rapid and easy unfolding, refolding [58], and large-scale opening of genome loci loop that affects and explains biological processes, e.g. the connection of distal single-nucleotide polymorphisms (SNPs) with their target genes, gene activation, gene repression, or the cell cycle [59,60,61,62,63].

When studying these two globules, two biophysical properties are considered: the genomic distance between two loci and the probability of contact between them. It is worth noting that genomic distance (s) is measured by FISH and contact probability is obtained from chromosome conformation methods such as Hi-C. The equilibrium and fractal globules yield different estimates for these properties, and therefore also varying predictions on the three-dimensional distance between pairs of loci. Lieberman-Aiden et al. [30] and Mirny [37] reports, through simulation, that equilibrium and fractal globule scaling for three dimensional-distance are s1/2 and s1/3 (s: genomic distance - number of nucleotides between two loci), respectively. Equilibrium and fractal globule scaling for contact probability are s−3/2 and s−1, respectively. As shown in [37], the properties exhibited by the fractal globule model make it more effective at fitting Hi-C data than the equilibrium globule.

Some methods adopt the knowledge about polymer chain for chromosome structure representation by simulating a physically realistic, bead-chain polymer model of the 30-nm chromatin fiber [64, 65]. As a result, when constructing either a chromosome structure for instance, a locus for a chromosome is represented using a conventional beads-and-spring polymer model, where each bead represents a specific genomic location with well-defined initial and final genomic coordinates. Hence, viewing the chromatin fiber as a polymer model implies that conformation energies such as bending, stretching, and excluding energies of chromatin segments needs to be considered and integrated with the IF for 3-D structure reconstruction (Fig. 1a).

Fig. 1

Chromosome and genome 3-D structure representation for models from Hi-C data. The different models used for representing 3-D chromosome and genome structure by various methods using Hi-C data for modeling chromosomes and genome 3-D structure. a polymer model, b spheres, c points

Spheres and Points

An alternative structure representation model adopted by methods is representing the chromosome region or loci as series of connected spheres or interacting points. Methods using this approach presents the 3-D structure in a simplified model, where the spheres [66,67,68] or points [45, 46, 69, 70] are synonymous to a chromosome region or loci of a chromosome (Fig. 1b, c). Using a beads on string configuration, each bead is modeled as a spherical shape with a defined radius, and an excluded volume used to penalize overlaps between two spheres. The defined radius and the sphere volume could consequently be considered as a restraint to be satisfied during the algorithm’s 3-D structure reconstruction process. The Points representation represent the chromatin region simply as a point, with no radius nor volume, to mark the presence or absence of a loci.

Methodologies for Chromosome and Genome 3-D Structure Reconstruction

The methods for chromosome and genome 3-D structure inference are categorized below based on the IF modeling adopted by them. All methods adopt a stepwise approach to achieve the 3-D structure reconstruction, and a summarization of these steps is provided in Fig. 2. In addition, the key properties of these methods are summarized in Table 1.

Fig. 2

Chromosome and genome 3-D structure reconstruction workflow. A summarization of the steps for genome and chromosome 3-D structure taken by the different methods. Starting from the user input in Step 1: The input preparation, usually, Hi-C contact matrix or sometimes with extra parameters requirement. Step 2: One of the three IF modeling approach is used to represent the IF depending on the method’s algorithm. Step 3: Modeling is done using defined sampling algorithms, and Step 4, a consensus average structure or a group of structure is generated depending on the method’s structure class

Table 1 A comparison of the methods for reconstructing 3-D chromosome and genome structure from Hi-C data

Distance-Based Methods

Over the years, a number of approaches have been proposed for chromosome 3-D structure inference from Hi-C contact data. A group of these methods involve a two-step process: (1) IF is converted to distance, ultimately defining the problem of 3-D genome or chromosome structure reconstruction as a problem of converting distances into 3-D coordinates; and (2) non-linear optimization is subsequently applied to the problem in order to find the genomic coordinates that satisfy converted distances. The most notable differences between these proposed methods are: (1) the way in which IF is converted into distance, and (2) the optimization technique used to infer the 3-D structure from loci distance. The aim of a distance-based modeling is to create a map that shows the relative spatial positioning of a number of objects whose inter-point distance is known. Additionally, representing chromosome structure prediction as a distance-based modeling problem is tempting because methods based on distances are simple and clear: there is no ambiguity regarding metric definition and proximity between objects can eventually be derived. In relation to 3-D genome structure prediction, the distance-based approach makes it easier to handle a large spectrum of modeling problems at different Hi-C data resolutions.

The distance-based approach attempts to reproduce the original metric or distance as accurately as possible. The earliest application of the metric multi-dimensional scaling (MDS) [82, 99] to chromosome 3-D conformation construction, known as 5C3D [45], assumed that the relationship of IF to distance between DNA fragments or loci follows an inverse relation; it then used an optimization approach to find the best 3-D conformation through a misfit objective function of the converted distance and the 3-D Euclidean distance between points. While this method was applied to the 5C variant of 3C data, it could be applied to Hi-C datasets as well. Similarly, in their work based on yeast 3-D genome structure reconstruction, Duan et al. [66] designed a metric that estimated the corresponding Euclidean distance from the mean of the curves obtained from two restriction enzyme libraries for each contact frequency. To aid modeling and ensure that intra- and inter-chromosomal features (e.g., centromeres), distance, and properties were satisfied [66, 67], researchers introduced a series of constraints such as minimum and maximum distances between adjacent beads, minimum distances between pair beads to avoid overlapping and clashes, specific positioning of RNA coding regions, telomeres, and centromeres to guide the construction of the 3-D model; this constituted an improvement over the previous method. Duan et al. used IPOPT [71], an open-source software for nonlinear constrained optimization problems, to minimize the objective function; this ensured that the predicted coordinates of two interacting loci, from which the distance between said loci in the 3-D structure is derived, closely matched the expected distance obtained from IF. Tanizawa et al. [67] developed a method similar to [66] to construct the 3-D structure of the fission yeast genome.

Although Lieberman-Aiden et al. [30] showed that IF can be used to determine the spatial distance between interacting loci, certain factors regarding this conversion are still worth considering. As shown by [76, 100,101,102] in their work, the IF-distance correlation might vary from one dataset resolution to another, and from one organism to another. Hence, an efficient method is required for a distance-based approach to generate a more reliable distance estimate from IF data. To solve this problem, Zhang et al. [76] made two novel propositions for the two-step genome structure prediction pipeline. First, they used a modified version of the golden section search method [103] to determine the best scale parameter, conversion factor (α), to convert IF to its approximate distance equivalent: DijFijα; this ensures that an appropriate conversion factor is obtained for each dataset. Secondly, for the 3-D structure prediction from a distance matrix, they presented an algorithm called ChromSDE (Chromosome Semi-Definite Embedding). Unlike earlier methods, ChromSDE relaxed the optimization problem to a semi-definite programming (SDP) problem. The proposed approach to IF-distance conversion defined by Zhang et al. introduced a new convention for defining the IF-distance relationship, followed by a series of distance-based algorithms that were subsequently developed.

According to Yaffe and Tanay [104], raw Hi-C data obtained from 3C experiments may contain numerous systematic biases, such as GC content, length of restriction fragments, and mappability between fragments. Long-range frequencies are typically noisy and unreliable; this represents a substantial drawback for the construction of 3-D chromosome and genome structures. In order to overcome these limitations, a number of methods have been developed to pre-process Hi-C data through normalization [9, 42, 104,105,106,107,108] before using the data for 3-D reconstruction. Alternatively, certain algorithms for 3-D structure construction incorporate bias removal. Peng et al. [77] proposed a normalization approach to reduce experimental sequencing depth bias, which affects the IF yielded by Hi-C data and makes it hard to compare structures from data obtained from different experiments. The method, called AutoChrom3D, provides an automated pipeline for 3-D modeling, enabling structural comparison at various data resolutions. Two linear transformations were used to determine the frequency-distance correlation, and structure was predicted through nonlinear constrained optimization. Shavit et al. [81] designed an MDS-based optimization approach that used FISH distance to guide the conversion of IF to Hi-C loci distances; this approach aimed to reduce noise, improve the data quality, ensure the consistency of data used for 3-D structure construction, and cover key functionality features in the Hi-C and FISH datasets, which will eventually overlap if these features are vital. Zou et al. [47] designed a flexible algorithm capable of handling biases introduced by restriction enzymes during Hi-C data sequencing. Restriction enzymes are known to have various cutting sites across the genome, so combining different Hi-C tracks provides further information about genomic loci for modeling. The tool developed by Zou et al., called HSA, takes advantage of the uniqueness of the contact map obtained from different restriction enzymes in Hi-C experiments; it creates a generalized linear model through an iterative algorithm that combines simulated annealing and Hamiltonian dynamics. By using HSA, Zou et al. discovered that the obtained 3-D structure fits the contact map obtained from different restriction enzymes. Bau et al. [72] performed a log transformation and the Z-score computation to normalize the contact counts. They converted observed interactions between loci to points and spatial restraints, and used the Integrative Modeling Platform (IMP) [73] to produce possible confirmations that satisfies their defined constraints and maximizes their structure to fit the IF data. Each loci was first represented as a point connected by a “string” to create a pairwise interaction in which the length of the string depended on the number of interactions between the loci.

To date, a number of other distance-based methods have been developed. These algorithms create 3-D models by first converting contact frequency to distance [9, 46, 69, 70, 77, 88, 97, 109, 110] and then apply optimization to predict chromosome structure. Usually, these methods perform chromosome 3-D reconstruction by first defining a random 3-D structure; this structure coordinates are then updated by an objective function that is iteratively optimized until a convergence condition is satisfied. Chromosome3D [46], applied a modified version of the distance geometry simulated annealing (DGSA) based method for chromosome and genome 3-D structure reconstruction from Hi-C data. The DGSA method has been popularly used for protein structure construction over the years and implemented in the Crystallography & NMR System (CNS) suite [111, 112]. The Hi-C distances are used as restraints for the defined simulated annealing (SA) optimization pipeline. SA is carried out through multiple steps of temperature change until the defined structure energy is optimally minimized. Because Chromosome3D uses one of the rigorously tested approaches in protein structure to inferring chromosome and genome 3-D structure, it is reliable and robust against noise in Hi-C data.

LorDG [69] introduced a novel method to address inconsistent chromosomal contacts generated from multi-cell Hi-C data. It used a nonlinear Lorentzian function as the objective function—to enforce the satisfaction of consistent restraints, which is resistant against noisy distance restraints. Unlike the square error function that is susceptible to outliers, LorDG aims to maximize the satisfaction of realistically satisfied restraints rather than unsatisfiable noisy ones. The objective function is optimized by the highly scalable adaptive step-size gradient descent method. Its resilience against noisy contacts and scalability make it a suitable method for constructing the structure of the entire genome involving noisy inter-chromosomal contacts. 3DMax [70] defined a maximum likelihood objective function for chromosome 3-D structure inference from Hi-C data. It is based on the simplified assumption that the contact data is normally distributed and that each Hi-C data point is conditionally independent given a structure. A log likelihood objective function for chromosome structure reconstruction was defined in order to determine the structure that maximizes the likelihood function. 3DMax uses a variant of gradient ascent called Adagrad [113] that adapts the learning rate to each objective function parameter automatically to regulate its learning rate. 3DMax is robust against noise and structural variability, and it is computationally fast and memory efficient.

miniMDS [92] and Hierarchical3DGenome [98] are the distance-based algorithms that reconstruct high-resolution 3-D models at the topologically associating domain (TAD) level. Eventually, these TAD models are assembled to form a complete, high-resolution 3-D chromosomal structure. After the assembly of TAD models, Hierarchical3DGenome uses the contacts between all regions in a chromosome to further refine the assembled whole chromosome model, which leads to high-resolution (e.g. 5 KB) models of good quality.

The conformational space of a chromosomal structure is large, given that Hi-C data are drawn from a population of cells, each with its own independent and unique 3-D structure. Hence, an ensemble of predicted structures obtained through so-called ensemble-based modeling appears to provide a better representation of chromosomal structure than a single structure obtained through consensus modeling. Unfortunately, like Hi-C data at large, this dataset contains a number of biases: the fact that it is noisy, coupled with other technical factors, makes it extremely difficult to determine the various unique 3-D structures of cells used in Hi-C experiments. Due to the drawbacks involved in using multi-cell Hi-C data, studying single-cell Hi-C data has become increasingly relevant [34]. In particular, it does not require designing an algorithm to satisfy the variability of each cell used in the Hi-C experiment. As expected, single-cell Hi-C datasets are sparser than multi-cell Hi-C datasets. Hence, conventional distance- and restraint-based methods are not suitable for 3-D structure reconstruction based on these data. Carstens et al. [90] extended Rieping et al.’s [114] Bayesian probabilistic framework to statistically infer ensembles of 3-D chromosome structures from single-cell Hi-C data using MCMC sampling. They combined single-cell Hi-C contact information with FISH data and a coarse grained model of the chromatin fiber. Lesne et al. [79] formulated a two-step algorithm known as “shortest-path reconstruction in 3-D” (ShRec3D), which combines the shortest-path distance between two points from graph theoretic methods with MDS to achieve chromosome reconstruction. This method is designed for both multi-cell and single-cell Hi-C data. In the case of single-cell Hi-C data, instead of distances between two points, binary numbers signify the presence or absence of interaction. ShRec3D+ [96] extended Lesne et al.’s algorithm by using a golden-section algorithm (an approach similar to Zhang et al. [76]) with an adaptable distance conversion factor for different Hi-C chromosome datasets. Wang et al. [64] proposed a method that combined knowledge of the conformational energy model of a chromatin structure and a Bayesian inference approach. They represented the chromosome structure as a polymer model with a conformational energy, and integrated the IF data as input for an expectation maximization based algorithm under a Bayesian like framework. They took advantage of the prior information about the conformation energy to construct a Bayesian inference of the chromatin structure. An approach proposed by Paulsen et al. [84] employed manifold-based optimization (MBO), which is basically the application of optimization techniques to the manifold of positive semi-definite matrices of fixed rank [115]. Paulsen et al. reported that MBO is capable of generating a consensus 3-D chromosome structure consistent with the original contact map.

Another approach for solving the distance-based problem is called non-metric multidimensional scaling (NMMDS), which assumes that only distance ranks are known; distances themselves are not provided. The method aims to yield a map of these ranks [116, 117]. Using this approach, Ben-Elazar et al. [118] developed a method for structure prediction based on the hypothesis that a pair locus A with a higher IF is closer in 3-D space than any other locus pair B with a lower frequency. Varoquaux et al. [78] also proposed an optimization method to solve the NMDS problem by minimizing the Shepard-Kruskal scaling cost function [119].

Contact Based Methods

Certain methods do not convert IF but use it directly for modeling. These methods are regarded as contact-based methods [15, 80, 83, 91, 93]. MOGEN [49, 80] used contact directly and designed an optimization-based approach that relied mostly on Hi-C intra- and inter-chromosomal contact data to build an ensemble of 3-D conformations for genome and chromosome structures. The contact-based optimization is carried out by the adaptive step-size gradient descent/ascent method that is highly scalable and therefore is well suited for large-scale genome structure modeling. MOGEN does not require two contacted regions to satisfy a specific distance as the distance-based approach does. Instead it only tries to make the distance between the two contact regions below a threshold (i.e. in contact). MOGEN is capable of producing ensemble models that are highly consistent with each other. MOGEN is also robust against noise in the data, particularly the noise in inter-chromosomal contacts, and therefore it is able to build 3-D structures of large genomes such as the entire human genome. Gen3D [83] used a series of meta-heuristic algorithms (e.g. genetic algorithms and simulated annealing) to infer 3-D structure from IF. Zhu et al. [93] proposed a manifold-based framework called GEM, which first uses IF to create an interaction network representing the spatial organization of the loci from Hi-C data. Zhu et al.’s aim was to use a manifold learning algorithm to uncover the low-dimensional (3-D) geometry embedded in a high-dimensional (Hi-C) space, while satisfying certain defined conformation energy requirements. An improvement over this method integrates Hi-C data with FISH data for 3-D structure inference [94]. To ensure the modeling of realistic structures consistent with cellular organization, Paulsen et al. [91] introduced Chrom3D, a genome-modeling algorithm that combines Hi-C and Lamina-associated domain (LAD) information from ChIP-seq data to generate an ensemble of 3-D genome structures in which loci and TAD positioning and interaction requirements are satisfied.

On the other hand, certain methods convert contact frequency into defined spatial restraints. As is the case with distance-based approaches, these restraints are satisfied through an optimization method. In their seminal study, Kalhor et al. [68] developed a 3C variant known as tethered conformation capture (TCC), aimed at increasing the signal-to-noise ratios in conformation capture experiments. This is relevant because it allows for a more accurate representation of IF, especially for genome structure analysis, where low inter-chromosomal interactions are recorded using existing approaches. Using TCC data, researchers proposed a novel modeling approach whereby a variety of genome structures were generated. This approach, called population-based modeling, produces a population of structures representative of genomic configuration and consistent with contact probability. Serra et al. [85] followed certain constraints in order to transform IF into spatial restraints; for instance, consecutive and non-consecutive loci were treated differently. As in the case of Bau et al. [72], these restraints were satisfied by using the IMP.

Probability Based Methods

Methods in this category define a probabilistic measure for contact frequency, hence their name. Using a probabilistic approach to model 3-D structures has a number of advantages; key among them is that such an approach allows uncertainties in experimental Hi-C data to be easily considered through probabilistic representation. In addition, statistical calculations of specific structural properties or noise sources can be carried out. Due to the fact that Hi-C data are drawn from cell populations, IF can be considered as an average; most probability-based methods assume that an ensemble of structures underlies a contact map. In addition, they consider the problem of 3-D structure inference as either a Bayesian inference problem or a maximum likelihood problem. However, some probabilistic modeling may be more time consuming than other methods.

Rousseau et al. [48] developed the first method in this category, called MCMC5C. They defined a probabilistic model of IF and used a Markov chain Monte Carlo (MCMC) sampling to generate an ensemble of structures. MCMC5C through a Gaussian model based on Hi-C data, whose variance was estimated using an improvised approach. A MCMC sampling-based algorithm was selected over alternatives methods because of its inherent ability to estimate the distribution of various structural properties. As previously mentioned, raw Hi-C data contain a number of systematic biases such as GC content, restriction enzyme cutting frequency, and sequence uniqueness [104]. These factors all need to be considered when designing a 3-D genome reconstruction method. To overcome these limitations, Hu et al. [75] proposed two Bayesian models for 3-D genome structure reconstruction from Hi-C data. Their methods combined bias removal with 3-D genome structure construction. They corrected known biases and used a Poisson model to fit contact data, an improvement over MCMC5C when it came to estimating the Gaussian variance. Varoquaux et al. [78] also defined a probabilistic model of IF. Similar to the model defined by Hu et al., it defined the structure inference problem as a maximum likelihood problem and used an optimization method to solve it.

A typical drawback of high-resolution Hi-C data is the sparsity of long-range contacts on the contact matrix and the high proportion of zero-contact counts between loci in the matrix. Hence, certain existing methods might be incapable of modeling at a higher resolution. Park and Lin [87] proposed an algorithm that is robust to resolution specification and corrects known systematic biases. They modeled the contact count using a Poisson distribution and addressed excess zero problems in high resolution datasets. They suggested that these problems could be solved by adjusting the Poisson distribution adopted for modeling.

Nagano et al. and Stevens et al. [34, 120, 121] applied a simulated annealing technique to sample single-cell datasets, while sometimes using contacts as distance restraints at different data resolutions. A novel study by Tjong et al. [86] has proposed a population-based modelling approach called PGS. Different from the ensemble-based approach—where a variety of structures with different variabilities are generated to simulate the heterogeneity of cells in the Hi-C experiment—the population of genome structures generated by PGS is consistent with the normalized contact probability matrix. Tjong et al. have formulated a probabilistic framework that uses an EM algorithm with constraint assignment at the E step and optimization of the structure population through simulated annealing and conjugate gradient descent at the M step. This method takes advantage of other external experimental data, such as lamina information for improved modeling. Rosenthal et al. [95] proposed an approach to recover missing contacts in single-cell Hi-C contact maps by filling missing parts with structures obtained from the corresponding cell populations, while imposing certain penalties on the generated structures.

Correcting Biases in Hi-C Data by Data Normalization

As is the case for most sequencing experiments, raw Hi-C data contain several systematic biases that could potentially affect the 3-D genome reconstruction. An inexhaustive list of these systematic biases include GC content, distance between restriction sites, restriction enzyme cutting frequency, sequence uniqueness, and experimental artifacts [104]. In a Hi-C experiment protocol, a minimum of 25 million cells was used to produce a Hi-C library [27, 30, 38, 69] with the goal of analyzing the contact frequencies between genomic sites in a cell population. One of the reasons for using a population of cells in Hi-C experiments is more sequence reads can be produced from a population of cells than a single cell.

The number of paired-end reads linking two genomic regions is interpreted as the interaction frequency between two genomic regions. This implies that a higher interaction frequency on a contact map means that a higher read count was observed, and that the two regions are spatially close to each other. However, many of these systematic biases affect the observed Hi-C read counts for two interacting regions (or fragments) on a contact matrix [106]. Hence, when these biases are left unhandled, the 3-D model construction is predicated on inaccurate information and consequently may be adversely affected. Additionally, if the effect of duplication, deletion, inversion and ploidy is significant in the pair reads, this could cause a direct effect on the number of paired-end reads linking two genomic regions which will alter the derived contact map. Because the Hi-C contact data is used for 3-D genome modeling, the level of correctness of the Hi-C data largely determines the accuracy of the generated model.

To overcome these limitations, most 3-D reconstruction methods apply normalization methods that focus on removing biases introduced by experimental procedures and by intrinsic properties of the genome to preprocess the data [9, 42, 104,105,106,107,108]. With the application of a normalization and pre-processing technique before 3-D genome reconstruction, the noise and systematic biases introduced by external factors, such as DNA shearing, and cutting, during the Hi-C experiment makes the Hi-C data more suitable for chromosome/genome 3-D structure reconstruction. Alternatively, some probability-based reconstruction methods handle the noise and biases differently by taking the biases into consideration in their algorithm design [75].

A common problem observed in some Hi-C data is the omission of the contact frequency of some genomic positions in the contact matrix. When this occurs, the reconstructed 3-D model from this data varies across the different tools due to difference in the way the methods represent omissions in their 3-D model. Generally, this leaves some doubt about which 3-D model is better when this occurs.

Validation and Evaluation

According to the literature on chromosome and genome 3-D construction methods, algorithms are most often validated by a simulated dataset to assess their reconstruction ability, the consistency with the Hi-C data, known genome and chromosome structural features [49], or Fluorescence in situ hybridization (FISH) data. In the simulation case, most methods use a 3-D polymer model meant to serve as a gold standard model with which to compare the final 3-D reconstructed structure. A set of chromosomal contact data is then simulated from this structure, and a certain degree of Gaussian noise is often added to the data as well. The noise is usually added to assess the methods’ responsiveness and accuracy to noisy data. Eventually, the algorithms’ ability to reconstruct the true model is tested. A commonly used synthetic dataset is the one generated by Trussart et al. [122]. Trussart et al. created a series of simulated Hi-C contact matrices in which genomic architectures are pre-defined, and the noise level and structural variability (SV) are both simulated.

FISH provides a powerful tool for identifying the location of a DNA sequence. It is used to study the 3-D organization of chromosomes and genomes and determine the proximity of a gene relative to other genes through the use of fluorescent probes [123]. It has been determined to be much more accurate, simple, and reliable than all other molecular profiling techniques [124]. Hence, it is often used to determine the distance between loci in a genome and for single-cell analysis of gene and loci positioning [125,126,127,128]. However, its major limitations are low throughput and resolution at higher scales, such as the entire genome or an ensemble of cells. Nonetheless, FISH data can be used to validate the distance between loci in a reconstructed 3-D structure at a lower scale. Given that the FISH method is considered reliable, it is useful in the study of chromosomal and genomic 3-D spatial organization when loci in the structure being evaluated are physically proximal.

Once the structure construction is complete, a method is often needed to assess its accuracy. The most common approach to structure evaluation is to calculate the Pearson correlation coefficient (PCC), the Spearman correlation coefficient (SCC), or the root mean square error (RMSE) of the distance representation of the Hi-C data and the Euclidean distance of the 3-D chromosomal structure. Since these metrics are obtained for distance, they are sometimes referred to as the distance Pearson correlation coefficient (dPCC), the distance Spearman Correlation Coefficient (dSCC), and the distance root mean square error (dRMSE). The value of dSCC and dPCC is in the range of − 1 to + 1, with higher values being preferable. In the case of dRMSE, on the other hand, a lower value is preferred. The latter may vary between 0—which signifies no difference between distances—and a large upper limit dependent on the number of fragments in the structure being compared when they are completely different. The dRMSE is also an appropriate metric to assess the similarity between 3-D structures. In order to do so, a linear transformation that includes translation, orthogonal rotation, and rescaling is performed on one of the structures, so that they are at the same 3-D-coordinate scale as in [49].

Let the pairwise distance between Hi-C data IF be represented by the vector {Di, …, Dn} and the Euclidean distance between loci in a 3-D chromosome model be represented as {EDi, …, EDn}, where n is the number of loci pairwise distances. The dSCC, dPCC, and dRMSE can be computed as shown below:

  1. (1)

    The dPCC is defined as:

    $$ \mathrm{dPCC}=\frac{\sum_{i=1}^n\left({D}_i-\overset{\acute{\mkern6mu}}{D}\right)\left({ED}_i-\overset{\acute{\mkern6mu}}{ED}\right)}{\sqrt{\sum_{i=1}^n{\left({D}_i-\overset{\acute{\mkern6mu}}{D}\right)}^2{\sum}_{i=1}^n{\left({ED}_i-\overset{\acute{\mkern6mu}}{ED}\right)}^2}} $$


  • Di and EDi are single distance samples indexed with i,

  • n is the number of loci pairwise distances,

  • \( \overset{\acute{\mkern6mu}}{D} \) and \( \overset{\acute{\mkern6mu}}{ED} \) represent sample means. \( \overset{\acute{\mkern6mu}}{D}=\frac{1}{n}{\sum}_{i=1}^n{D}_i \), \( \overset{\acute{\mkern6mu}}{ED} \) = \( \frac{1}{n}{\sum}_{i=1}^n{ED}_i \) .

  1. (2)

    The dSCC is defined as:

$$ \mathrm{dSCC}=\frac{\sum_{i=1}^n\left({A}_i-\overset{\acute{\mkern6mu}}{A}\right)\left({B}_i-\overset{\acute{\mkern6mu}}{B}\right)}{\sqrt{\sum_{i=1}^n{\left({A}_i-\overset{\acute{\mkern6mu}}{A}\right)}^2{\sum}_{i=1}^n{\left({B}_i-\overset{\acute{\mkern6mu}}{B}\right)}^2}} $$

dSCC is calculated by converting distance variable Di and EDi into ranked variables Ai and Bi i, and then, computing the dPCC between the ranked variables. Hence, the pairwise distances Di and EDi are converted into ranked variables Ai and Bi respectively,


  • Ai and Bi are the ranks of two distances, Di and EDi respectively.

  • \( \overset{\acute{\mkern6mu}}{A} \) and \( \overset{\acute{\mkern6mu}}{B} \) represent sample means of rank. \( \overset{\acute{\mkern6mu}}{A}=\frac{1}{n}{\sum}_{i=1}^n{A}_i \), \( \overset{\acute{\mkern6mu}}{B} \) = \( \frac{1}{n}{\sum}_{i=1}^n{B}_i \) .

  1. (3)

    The dRMSE is defined as:

$$ \mathrm{dRMSE}=\sqrt{\frac{1}{n}\sum {\left({D}_{ij}-{ED}_{ij}\right)}^2} $$
  • where Dij and EDij represent the pairwise distance between loci i and j of the Hi-C IF data and 3-D structure Euclidean distance

  • n is the number of loci pairwise distances.

Microscopy-Based Techniques for Studying Genome Organization

Although this review highlights the methods for genome structure reconstruction from Hi-C data, it is noteworthy to examine the complementary imaging methods used for studying the genome organization before and after the emergence of high-throughput sequencing techniques. For many years, the structure of the genome has been studied through various microscopy techniques [23, 129,130,131,132,133,134] which can be broadly divided into electron and light microscopy.

The light microscope alternatively referred to as the optical microscope is a well-known research tool that uses visible light to detect small objects. Over the years, light microscopy has greatly enhanced the study of the events and the structural details in the cell nucleus. However, the light microscopy techniques have a well-known limitation for being unable to overcome the diffraction barrier. As a solution to this, several strategies have been proposed to bypass the diffraction barrier of light microscopy and increase resolution. These strategies are called the super resolution microscopy strategies. They include saturated structured illumination microscopy (SSIM), stimulated emission depletion (STED), and ground state depletion (GSD) [20, 21]. The introduction of Stochastic super-resolution microscopy techniques such as Photo-activated localization microscopy (PALM or FPALM) and stochastic optical reconstruction microscopy (STORM) ushered in a new wave of discovery about the genome organization [22, 132, 135]. These techniques allow obtaining images at a higher resolution because they are not limited by the diffraction barrier in optical microscopy. These methods use florescent probes for imaging in multiple colors and support the selection of many fluorescent molecules at a very high resolution to build point by point images that display the relationship between points [135]. The STORM and PALM techniques elevated the visualization of the genome structure to an incredibly high resolution. Ricci et al. [136] used the STORM technique to visualize the chromatin fiber structure of different cells at a nanoscale resolution, single cell level, which revealed nucleosome groups along the chromatin fiber which they called “nucleosome clutches”

One type of light microscopy technique, fluorescence microscope, uses fluorescence and phosphorescence to study the properties of and visualize an object or cellular component of a cell. The fluorescence microscopy technique uses a light intensity that is significantly higher than other light microscopy techniques [137,138,139]. Fluorescence microscopy technique is effective at visualizing fluorescent dyers stains [140,141,142] as well as autofluorescence cellular structures i.e. biological structures which naturally emit fluorescent light [139]. The technique is also used when studying the expression and the localization of proteins using fluorescent antibodies in a biochemical strategy called immunofluorescence. The fluorescence dyers stains are used to determine cellular structure and identify specific targets of interest within a cell. A major limitation of the fluorescence microscopy technique is photobleaching. Photobleaching causes the fading of the dye or a fluorophore molecule making it lose its fluoresce properties, hence rendering the protein molecules or object invisible. Fluorescence recovery after photobleaching (FRAP) [143, 144], and Fluorescence loss in photobleaching (FLIP) [145] analysis are fluorescence microscopy technique used to examine diffusion and molecular movement respectively in a cell. FLAP, FRET and FLIM are also advanced fluorescence microscopy techniques that are used in biological and biomedical research [146].

For some time, the 3-D genome organization was largely discovered through the fluorescent in situ hybridization (FISH) technique. The FISH [2, 18, 19] technique uses a florescence probe to detect specific DNA (or RNA) sequences or selected genome loci in single cell nuclei by light microscopy. Today, there are different types of FISH, each with their specialized function e.g. DNA-FISH, RNA-FISH, cryo-FISH e.t.c. [147]. These variants are more prolific than FISH because of their accuracy, and reliability. FISH techniques allow the conceptualization of the arrangements of genetic materials in the cell nucleus. The FISH technique has revealed that the chromosomes occupy discrete territories in the cell nucleus, referred to as chromosome territories (CT) [2, 148], CT intermingle significantly in the nucleus of human cells [149], the influence of gene density and transcription on chromosome organization in the nucleus [150, 151], and the genome organization in the nucleus based on the partitioning of the chromosomes regions according to the gene distribution [152, 153]. The findings have increased the understanding about the genome architecture and behavior in the nucleus of the cell. However, the FISH technique can only be used to examine predetermined regions in cells. To resolve this, a fully automated FISH-based imaging pipeline called High-throughput imaging position mapping (HIPMap) was developed to perform high-precision, high-throughput, automated fluorescent in situ hybridization imaging of the spatial location of genome regions at large scale [154].

Electron microscopes uses a beam of high energy electrons to examine objects and obtain information about that object or a specimen. This provides the information about the surface characteristics, composition of the elements that makes up the object, the particles within the objects, and the arrangement of the atoms within the object. It was developed due to inability of the light microscope to examine the information about structure of smaller objects. The development of the electron microscope improved the resolution so that tiny objects e.g. atom can be observed under this microscope. To examine objects only observable at the higher resolution e.g. the examination of a cell nucleus, the Electron microscopic techniques such as the, Transmission Electron Microscope (TEM) where instead of using light to illuminate the specimen, a high energy electric beam is used. The scanning electron microscope (SEM), reflection electron microscope (REM), scanning transmission electron microscope (STEM), and the cryo-electron microscopy (Cryo-EM) are other forms electron microscopes techniques each with the unique method for how the structure and composition information is gained from the object [155,156,157,158,159]. Cryo-EM especially has produced very useful insights by enabling the determination of atomic resolution level macromolecular structures [160,161,162]. In protein structure research, Cryo-EM has been used to capture protein structure in its native state. Some methods have been developed to complement the microscopy techniques. Ou et al. [163] combined electron microscopy with a labeling method to reveal the 3-D organization across multiple scales in the cell nucleus They developed a method called ChromEMT that reveals the 3-D packing of the DNA in cells, and through their method reveled information about the DNA folding as it relates to the genome compaction in the cell nucleus.

For many years, the FISH and the microscopy-based techniques have given scientists insight about the spatial organization and architectural arrangement in the nucleus while providing explanation for nuclear positioning in the cell nucleus. Some of these findings include: the discovery of chromosome territories [2], the organization of gene clusters and their influence on transcription in the nucleus [51, 52], the segmentation of chromatin in the cell nucleus, for example, the active euchromatin and inactive heterochromatin occupy separate environments in the nucleus [164], and the existence of unique compartments that influences functional interaction [165,166,167,168].

These methods provide valuable information regarding the genomes organization that can be used as base information when constructing models with Hi-C data. For example, it is common practice to use the results of FISH experiments as validation for chromatin conformation models generated by Hi-C experiments. This can be done by verifying that the spatial distances observed between multiple FISH probes are consistent with the predicted distance between the corresponding genome bins found in the Hi-C conformation model.

Despite the advancements in the FISH and the microscopy approaches, they are limited to studying a region of genome, and do not provide a universal and comprehensive view of the 3-D genome architecture [169] of the whole genome. The need to study the genomic organization at a genome wide scale led to the development of the chromosome conformation capturing techniques. However, it is worth noting that the chromosome conformation capturing techniques and the imaging techniques of probing genome/chromosome structures are complementary and the latter can experimentally validate the former.

Summary and Future Insights

Our review of the methods for reconstructing the 3-D structure of the chromosome and genome has revealed that these methods can be largely categorized into three groups (distance-based methods, contact-based methods and probability-based methods) according to how IF is modeled. For each category, we have discussed their potential strength and weakness in reconstructing 3-D chromosome and genome structures. Although we have primarily grouped methods based on IF modeling, there are other ways they could be categorized. For instance, their classification could be based on the type of structure they generate [72, 78]. Methods that generate a single representative structure for a Hi-C dataset are consensus-based methods [66, 67]. Those that generate a variety of structures to represent the heterogeneity of Hi-C data are ensemble-based methods [45, 48]. Finally, population-based methods [68, 86, 89] generate a population of structures that, as a whole, is statistically consistent with the Hi-C data.

Despite the improvement in 3-D structure modeling approaches, the lack of a real structure with which to contrast these models remains a challenge. In particular, it is currently difficult to confirm the true modeling capability of 3-D genome methods. Although the introduction of 3-D-FISH data and Hi-C data for joint modeling has received some attention recently [94], there is no sufficient 3-D-FISH data to guide most modeling on Hi-C data and to thoroughly validate the quality of computational models. The development of more advanced genome/chromosome imaging techniques will further improve the validation of 3-D genome models. In addition, other high-throughput sequencing data such as functional genomics and epigenomics data can be used to validate the biological validity of 3-D genome/chromosome models by exploring their correlation with 3-D genomes.

Another challenge is to reconstruct high-resolution 3-D models of large genomes from Hi-C data,

which are needed for studying detailed interactions between genes and regulatory elements, due to enormous time complexity and data sparsity associated with high-resolution modeling. Only a few methods [98] was designed to build high-resolution (e.g. 5 KB) models.

Finally, it is important to make 3-D genome modeling methods easy for biomedical scientists to use in their research. To this end, a few tools have been designed to visualize 3-D genome models [88, 89, 170,171,172,173,174]. Recently, GenomeFlow [40] provides a comprehensive graphical environment for users to process Hi-C data, generate chromosomal contact maps, build 3-D models, and apply 3-D models to integrate various omics data. More efforts of making 3-D genome modeling accessible to general users are still needed.



Chromosome conformation capture




Crystallography & NMR System


Distance Pearson correlation coefficient


Distance root mean square error


Distance Spearman’s correlation coefficient


Fluorescence in situ hybridization


The name of a method an extension of 3C method that is capable of identifying read pair interactions on an “all-versus-all” basis—that is, it can profile interactions for all read pairs in an entire genome


Interaction frequency


Integrative Modeling Platform


Lamina-associated domain


Multidimensional scaling


Non-metric multidimensional scaling


Topologically associated domains


Tethered conformation capture


  1. 1.

    Misteli T. Beyond the sequence: cellular organization of genome function. Cell. 2007;128(4):787–800.

    CAS  PubMed  Article  Google Scholar 

  2. 2.

    Cremer T, Cremer C. Chromosome territories, nuclear architecture and gene regulation in mammalian cells. Nat Rev Genet. 2001;2(4):292.

    CAS  PubMed  Article  Google Scholar 

  3. 3.

    Branco MR, Pombo A. Chromosome organization: new facts, new models. Trends Cell Biol. 2007;17(3):127–34.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  4. 4.

    Hakim O, Misteli T. SnapShot: chromosome conformation capture. Cell. 2012;148(5):1068–e1.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Osório J. Chromosome biology: moving a TAD closer to unravelling chromosome architecture. Nat Rev Mol Cell Biol. 2015;16(12):701.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  6. 6.

    Dekker J, Mirny L. The 3D genome as moderator of chromosomal communication. Cell. 2016;164(6):1110–21.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    Shen Y, Yue F, McCleary DF, Ye Z, Edsall L, Kuan S, Wagner U, Dixon J, Lee L, Lobanenkov VV, Ren B. A map of the cis-regulatory sequences in the mouse genome. Nature. 2012;488(7409):116.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  8. 8.

    Makova KD, Hardison RC. The effects of chromatin organization on variation in mutation rates in the genome. Nat Rev Genet. 2015;16(4):213.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  9. 9.

    Cournac A, Marie-Nelly H, Marbouty M, Koszul R, Mozziconacci J. Normalization of a chromosomal contact map. BMC Genomics. 2012;13(1):436.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Wasserman WW, Sandelin A. Applied bioinformatics for the identification of regulatory elements. Nat Rev Genet. 2004;5(4):276.

    CAS  PubMed  Article  Google Scholar 

  11. 11.

    Taberlay PC, Achinger-Kawecka J, Lun AT, Buske FA, Sabir K, Gould CM, Zotenko E, Bert SA, Giles KA, Bauer DC, Smyth GK. Three-dimensional disorganization of the cancer genome occurs coincident with long-range genetic and epigenetic alterations. Genome Res. 2016;26(6):719–31.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Dekker J. Gene regulation in the third dimension. Science. 2008;319(5871):1793–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Dekker J, Marti-Renom MA, Mirny LA. Exploring the three-dimensional organization of genomes: interpreting chromatin interaction data. Nat Rev Genet. 2013;14(6):390.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    de Laat W, Grosveld F. Spatial organization of gene expression: the active chromatin hub. Chromosom Res. 2003;11(5):447–59.

    Article  Google Scholar 

  15. 15.

    Gorkin DU, Leung D, Ren B. The 3D genome in transcriptional regulation and pluripotency. Cell stem cell. 2014;14(6):762–75.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  16. 16.

    Woodcock CL, Dimitrov S. Higher-order structure of chromatin and chromosomes. Curr Opin Genet Dev. 2001;11(2):130–5.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  17. 17.

    Chromatin WA. San Diego: Structure and Function. San Diego, CA: Academic Press; 1998.

    Google Scholar 

  18. 18.

    Langer-Safer PR, Levine M, Ward DC. Immunological method for mapping genes on Drosophila polytene chromosomes. Proc Natl Acad Sci. 1982;79(14):4381–5.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  19. 19.

    Amann R, Fuchs BM. Single-cell identification in microbial communities by improved fluorescence in situ hybridization techniques. Nat Rev Microbiol. 2008;6(5):339.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  20. 20.

    Westphal V, Rizzoli SO, Lauterbach MA, Kamin D, Jahn R, Hell SW. Video-rate far-field optical nanoscopy dissects synaptic vesicle movement. Science. 2008;320(5873):246–9.

    CAS  PubMed  Article  Google Scholar 

  21. 21.

    Hell SW, Wichmann J. Breaking the diffraction resolution limit by stimulated emission: stimulated-emission-depletion fluorescence microscopy. Opt Lett. 1994;19(11):780–2.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Betzig E, Patterson GH, Sougrat R, Lindwasser OW, Olenych S, Bonifacino JS, Davidson MW, Lippincott-Schwartz J, Hess HF. Imaging intracellular fluorescent proteins at nanometer resolution. Science. 2006;313(5793):1642–5.

    CAS  PubMed  Article  Google Scholar 

  23. 23.

    Huang B, Babcock H, Zhuang X. Breaking the diffraction barrier: super-resolution imaging of cells. Cell. 2010;143(7):1047–58.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Dekker J, Rippe K, Dekker M, Kleckner N. Capturing chromosome conformation. Science. 2002;295(5558):1306–11.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    de Wit E, de Laat W. A decade of 3C technologies: insights into nuclear organization. Genes Dev. 2012;26(1):11–24.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  26. 26.

    Han J, Zhang Z, Wang K. 3C and 3C-based techniques: the powerful tools for spatial genome organization deciphering. Mol Cytogenet. 2018;11(1):21.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  27. 27.

    Schmitt AD, Hu M, Ren B. Genome-wide mapping and analysis of chromosome architecture. Nat Rev Mol Cell Biol. 2016;17(12):743.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    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(11):1348.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Dostie J, Richmond TA, Arnaout RA, Selzer RR, Lee WL, Honan TA, Rubio ED, Krumm A, Lamb J, Nusbaum C, Green RD. Chromosome conformation capture carbon copy (5C): a massively parallel solution for mapping interactions between genomic elements. Genome Res. 2006;16(10):1299–309.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. 30.

    Lieberman-Aiden E, Van Berkum NL, Williams L, Imakaev M, Ragoczy T, Telling A, Amit I, Lajoie BR, Sabo PJ, Dorschner MO, Sandstrom R. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009;326(5950):289–93.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    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(1):90.

    CAS  Article  Google Scholar 

  32. 32.

    Fullwood MJ, Liu MH, Pan YF, Liu J, Xu H, Mohamed YB, Orlov YL, Velkov S, Ho A, Mei PH, Chew EG. An oestrogen-receptor-α-bound human chromatin interactome. Nature. 2009;462(7269):58.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Li G, Fullwood MJ, Xu H, Mulawadi FH, Velkov S, Vega V, Ariyaratne PN, Mohamed YB, Ooi HS, Tennakoon C, Wei CL. ChIA-PET tool for comprehensive chromatin interaction analysis with paired-end tag sequencing. Genome Biol. 2010;11(2):R22.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  34. 34.

    Nagano T, Lubling Y, Stevens TJ, Schoenfelder S, Yaffe E, Dean W, Laue ED, Tanay A, Fraser P. Single-cell Hi-C reveals cell-to-cell variability in chromosome structure. Nature. 2013;502(7469):59.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  35. 35.

    Ron G, Globerson Y, Moran D, Kaplan T. Promoter-enhancer interactions identified from Hi-C data using probabilistic models and hierarchical topological domains. Nat Commun. 2017;8(1):2237.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  36. 36.

    Fraser P, Bickmore W. Nuclear organization of the genome and the potential for gene regulation. Nature. 2007;447(7143):413.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  37. 37.

    Mirny LA. The fractal globule as a model of chromatin architecture in the cell. Chromosom Res. 2011;19(1):37–51.

    CAS  Article  Google Scholar 

  38. 38.

    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;6(39):e1869.

    Google Scholar 

  39. 39.

    Ay F, Noble WS. Analysis methods for studying the 3D architecture of the genome. Genome Biol. 2015;16(1):183.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  40. 40.

    Trieu T, Oluwadare O, Wopata J, Cheng J. GenomeFlow: a comprehensive graphical tool for modeling and analyzing 3D genome structure. Bioinformatics. 2018;

    Article  Google Scholar 

  41. 41.

    Durand NC, Shamim MS, Machol I, Rao SS, Huntley MH, Lander ES, Aiden EL. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst. 2016;3(1):95–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Servant N, Varoquaux N, Lajoie BR, Viara E, Chen CJ, Vert JP, Heard E, Dekker J, Barillot E. HiC-Pro: an optimized and flexible pipeline for Hi-C data processing. Genome Biol. 2015;16(1):259.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  43. 43.

    Castellano G, Le Dily F, Pulido AH, Beato M, Roma G. Hi-Cpipe: a pipeline for high-throughput chromosome capture. bioRxiv. 2015:020636.

  44. 44.

    Wingett S, Ewels P, Furlan-Magaril M, Nagano T, Schoenfelder S, Fraser P, Andrews S. HiCUP: pipeline for mapping and processing Hi-C data. F1000Research. 2015:4.

  45. 45.

    Fraser J, Rousseau M, Shenker S, Ferraiuolo MA, Hayashizaki Y, Blanchette M, Dostie J. Chromatin conformation signatures of cellular differentiation. Genome Biol. 2009;10(4):R37.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  46. 46.

    Adhikari B, Trieu T, Cheng J. Chromosome3D: reconstructing three-dimensional chromosomal structures from Hi-C interaction frequency data using distance geometry simulated annealing. BMC genomics. 2016;17(1):886.

    PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Zou C, Zhang Y, Ouyang Z. HSA: integrating multi-track Hi-C data for genome-scale reconstruction of 3D chromatin structure. Genome Biol. 2016;17(1):40.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  48. 48.

    Rousseau M, Fraser J, Ferraiuolo MA, Dostie J, Blanchette M. Three-dimensional modeling of chromatin structure from interaction frequency data using Markov chain Monte Carlo sampling. BMC Bioinformatics. 2011;12(1):414.

    PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Trieu T, Cheng J. Large-scale reconstruction of 3D structures of human chromosomes from chromosomal contact data. Nucleic Acids Res. 2014;42(7):e52.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  50. 50.

    Flory PJ. Principles of Polymer Chemistry. Ithaca: Cornell University Press; 1953.

    Google Scholar 

  51. 51.

    Gennes PG d. Scaling Concepts in Polymer Physics. Ithaca: Cornell University Press; 1979.

    Google Scholar 

  52. 52.

    Doi M, Edwards SF. The Theory of Polymer Dynamic. Oxford: Clarendon; 1986.

    Google Scholar 

  53. 53.

    Mateos-Langerak J, Bohn M, de Leeuw W, Giromus O, Manders EM, Verschure PJ, Indemans MH, Gierman HJ, Heermann DW, Van Driel R, Goetze S. Spatially confined folding of chromatin in the interphase nucleus. Proceedings of the National Academy of Sciences. 2009:pnas-0809501106.

  54. 54.

    Münkel C, Langowski J. Chromosome structure predicted by a polymer model. Phys Rev E. 1998;57(5):5888.

    Article  Google Scholar 

  55. 55.

    Barbieri M, Chotalia M, Fraser J, Lavitas LM, Dostie J, Pombo A, Nicodemi M. A model of the large-scale organization of chromatin. Biochem Soc Trans. 2013;41:508–12.

    CAS  PubMed  Article  Google Scholar 

  56. 56.

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

    CAS  Article  Google Scholar 

  57. 57.

    Bölinger D, Sułkowska JI, Hsu HP, Mirny LA, Kardar M, Onuchic JN, Virnau P. A Stevedore's protein knot. PLoS Comput Biol. 2010;6(4):e1000731.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  58. 58.

    Van Holde KE. Chromatin: Springer series in molecular biology. New York: Springer-Verlag; 1988.

  59. 59.

    Woodcock CL, Ghosh RP. Chromatin higher-order structure and dynamics. Cold Spring Harb Perspect Biol. 2010;2(5):a000596.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  60. 60.

    Sewitz SA, Fahmi Z, Lipkow K. Higher order assembly: folding the chromosome. Curr Opin Struct Biol. 2017;42:162–8.

    CAS  PubMed  Article  Google Scholar 

  61. 61.

    Dina C, Meyre D, Gallina S, Durand E, Körner A, Jacobson P, Carlsson LM, Kiess W, Vatin V, Lecoeur C, Delplanque J. Variation in FTO contributes to childhood obesity and severe adult obesity. Nat Genet. 2007;39(6):724.

    CAS  PubMed  Article  Google Scholar 

  62. 62.

    Scuteri A, Sanna S, Chen WM, Uda M, Albai G, Strait J, Najjar S, Nagaraja R, Orrú M, Usala G, Dei M. Genome-wide association scan shows genetic variants in the FTO gene are associated with obesity-related traits. PLoS Genet. 2007;3(7):e115.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  63. 63.

    Norton HK, Phillips-Cremins JE. Crossed wires: 3D genome misfolding in human disease. J Cell Biol. 2017;216(11):3441–52.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  64. 64.

    Wang S, Xu J, Zeng J. Inferential modeling of 3D chromatin structure. Nucleic acids research. 2015;43(8):e54.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  65. 65.

    Hua N, Tjong H, Shin H, Gong K, Zhou XJ, Alber F. Producing genome structure populations with the dynamic and automated PGS software. Nat Protoc. 2018;13(5):915.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  66. 66.

    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(7296):363.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  67. 67.

    Tanizawa H, Iwasaki O, Tanaka A, Capizzi JR, Wickramasinghe P, Lee M, Fu Z, Noma KI. Mapping of long-range associations throughout the fission yeast genome reveals global genome organization linked to transcriptional regulation. Nucleic Acids Res. 2010;38(22):8164–77.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  68. 68.

    Kalhor R, Tjong H, Jayathilaka N, Alber F, Chen L. Solid-phase chromosome conformation capture for structural characterization of genome architectures. Nat Biotechnol. 2012;30(1):90.

    CAS  Article  Google Scholar 

  69. 69.

    Trieu T, Cheng J. 3D genome structure modeling by Lorentzian objective function. Nucleic Acids Res. 2016;45(3):1049–58.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  70. 70.

    Oluwadare O, Zhang Y, Cheng J. A maximum likelihood algorithm for reconstructing 3D structures of human chromosomes from chromosomal contact data. BMC Genomics. 2018;19(1):161.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  71. 71.

    Wächter A, Biegler LT. On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming. Math Program. 2006;106(1):25–57.

    Article  Google Scholar 

  72. 72.

    Baù D, Sanyal A, Lajoie BR, Capriotti E, Byron M, Lawrence JB, Dekker J, Marti-Renom MA. The three-dimensional folding of the α-globin gene domain reveals formation of chromatin globules. Nat Struct Mol Biol. 2011;18(1):107.

    PubMed  Article  CAS  Google Scholar 

  73. 73.

    Alber F, Dokudovskaya S, Veenhoff LM, Zhang W, Kipper J, Devos D, Suprapto A, Karni-Schmidt O, Williams R, Chait BT, Rout MP. Determining the architectures of macromolecular assemblies. Nature. 2007;450(7170):683.

    CAS  PubMed  Article  Google Scholar 

  74. 74.

    Meluzzi D, Arya G. Recovering ensembles of chromatin conformations from contact probabilities. Nucleic Acids Res. 2012;41(1):63–75.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  75. 75.

    Hu M, Deng K, Qin Z, Dixon J, Selvaraj S, Fang J, Ren B, Liu JS. Bayesian inference of spatial organizations of chromosomes. PLoS Comput Biol. 2013;9(1):e1002893.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  76. 76.

    Zhang Z, Li G, Toh KC, Sung WK. Inference of spatial organizations of chromosomes using semi-definite embedding approach and Hi-C data. In: Annual international conference on research in computational molecular biology. Berlin, Heidelberg: Springer; 2013. p. 317–32.

    Chapter  Google Scholar 

  77. 77.

    Peng C, Fu LY, Dong PF, Deng ZL, Li JX, Wang XT, Zhang HY. The sequencing bias relaxed characteristics of Hi-C derived data and implications for chromatin 3D modeling. Nucleic Acids Res. 2013;41(19):e183.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  78. 78.

    Varoquaux N, Ay F, Noble WS, Vert JP. A statistical approach for inferring the 3D structure of the genome. Bioinformatics. 2014;30(12):i26–33.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  79. 79.

    Lesne A, Riposo J, Roger P, Cournac A, Mozziconacci J. 3D genome reconstruction from chromosomal contacts. Nat Methods. 2014;11(11):1141.

    CAS  PubMed  Article  Google Scholar 

  80. 80.

    Trieu T, Cheng J. MOGEN: a tool for reconstructing 3D models of genomes from chromosomal conformation capturing data. Bioinformatics. 2015;32(9):1286–92.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  81. 81.

    Shavit Y, Hamey FK, Lio P. FisHiCal: an R package for iterative FISH-based calibration of Hi-C data. Bioinformatics. 2014;30(21):3120–2.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  82. 82.

    de Leeuw J. Applications of convex analysis to multidimensional scaling. In: van Cutsem B, et al., editors. Recent advantages in Statistics. Amsterdam: North Holland Publishing Company; 1977.

    Google Scholar 

  83. 83.

    Nowotny J, Ahmed S, Xu L, Oluwadare O, Chen H, Hensley N, Trieu T, Cao R, Cheng J. Iterative reconstruction of three-dimensional models of human chromosomes from chromosomal contact data. BMC Bioinformatics. 2015;16(1):338.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  84. 84.

    Paulsen J, Gramstad O, Collas P. Manifold based optimization for single-cell 3D genome reconstruction. PLoS Comput Biol. 2015;11(8):e1004396.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  85. 85.

    Serra F, Baù D, Goodstadt M, Castillo D, Filion GJ, Marti-Renom MA. Automatic analysis and 3D-modelling of Hi-C data using TADbit reveals structural features of the fly chromatin colors. PLoS Comput Biol. 2017;13(7):e1005665.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  86. 86.

    Tjong H, Li W, Kalhor R, Dai C, Hao S, Gong K, Zhou Y, Li H, Zhou XJ, Le Gros MA, Larabell CA. Population-based 3D genome structure analysis reveals driving forces in spatial genome organization. Proc Natl Acad Sci. 2016;113(12):E1663–72.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  87. 87.

    Park J, Lin S. Impact of data resolution on three-dimensional structure inference methods. BMC Bioinformatics. 2016;17(1):70.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  88. 88.

    Szalaj P, Michalski PJ, Wróblewski P, Tang Z, Kadlof M, Mazzocco G, Ruan Y, Plewczynski D. 3D-GNOME: an integrated web service for structural modeling of the 3D genome. Nucleic Acids Res. 2016;44(W1):W288–93.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  89. 89.

    Szałaj P, Tang Z, Michalski P, Pietal MJ, Luo OJ, Sadowski M, Li X, Radew K, Ruan Y, Plewczynski D. An integrated 3-dimensional genome modeling engine for data-driven simulation of spatial genome organization. Genome Res. 2016;

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  90. 90.

    Carstens S, Nilges M, Habeck M. Inferential structure determination of chromosomes from single-cell Hi-C data. PLoS Comput Biol. 2016;12(12):e1005292.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  91. 91.

    Paulsen J, Sekelja M, Oldenburg AR, Barateau A, Briand N, Delbarre E, Shah A, Sørensen AL, Vigouroux C, Buendia B, Collas P. Chrom3D: three-dimensional genome modeling from Hi-C and nuclear lamin-genome contacts. Genome Biol. 2017;18(1):21.

    PubMed  PubMed Central  Article  Google Scholar 

  92. 92.

    Rieber L, Mahony S. miniMDS: 3D structural inference from high-resolution Hi-C data. Bioinformatics. 2017;33(14):i261–6.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  93. 93.

    Zhu G, Deng W, Hu H, Ma R, Zhang S, Yang J, Peng J, Kaplan T, Zeng J. Reconstructing spatial organizations of chromosomes through manifold learning. Nucleic Acids Res. 2018;46(8):e50.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  94. 94.

    Abbas A, He X, Zhou B, Zhu G, Ma Z, Gao JT, Zhang MQ, Zeng J. Integrating Hi-C and FISH data for modeling 3D organizations of chromosomes. bioRxiv. 2018;1:318493.

    Google Scholar 

  95. 95.

    Rosenthal M, Bryner D, Huffer F, Evans S, Srivastava A, Neretti N. Bayesian Estimation of 3D Chromosomal Structure from Single Cell Hi-C Data. BioRxiv. 2018;1:316265.

    Google Scholar 

  96. 96.

    Li J, Zhang W, Li X. 3D genome reconstruction with ShRec3D+ and Hi-C data. IEEE/ACM Trans Comput Biol Bioinform. 2018;1;15(2):460–8.

    PubMed  Article  Google Scholar 

  97. 97.

    Hua KJ, Ma BG. EVR: Reconstruction of Bacterial Chromosome 3D Structure Using Error-Vector Resultant Algorithm. bioRxiv. 2018;1:401513.

    Google Scholar 

  98. 98.

    Trieu T, Oluwadare O, Cheng J. Hierarchical Reconstruction of High-Resolution 3D Models of Large Chromosomes. Scientific reports. 2019;9(1):4971.

    PubMed  PubMed Central  Article  Google Scholar 

  99. 99.

    Borg I, Groenen P. Modern multidimensional scaling: theory and applications. J Educ Meas. 2003;40(3):277–80.

    Article  Google Scholar 

  100. 100.

    Ay F, Bunnik EM, Varoquaux N, Bol SM, Prudhomme J, Vert JP, Noble WS, Le Roch KG. Three-dimensional modeling of the P. falciparum genome during the erythrocytic cycle reveals a strong connection between genome architecture and gene expression. Genome Res. 2014;24:974.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  101. 101.

    Le TB, Imakaev MV, Mirny LA, Laub MT. High-resolution mapping of the spatial organization of a bacterial chromosome. Science. 2013;342(6159):731–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  102. 102.

    Fudenberg G, Mirny LA. Higher-order chromatin structure: bridging physics and biology. Curr Opin Genet Dev. 2012;22(2):115–24.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  103. 103.

    Kiefer J. Sequential minimax search for a maximum. Proc Am Math Soc. 1953;4(3):502–6.

    Article  Google Scholar 

  104. 104.

    Yaffe E, Tanay A. Probabilistic modeling of Hi-C contact maps eliminates systematic biases to characterize global chromosomal architecture. Nat Genet. 2011;43(11):1059.

    CAS  PubMed  Article  Google Scholar 

  105. 105.

    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(10):999.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  106. 106.

    Hu M, Deng K, Selvaraj S, Qin Z, Ren B, Liu JS. HiCNorm: removing biases in Hi-C data via Poisson regression. Bioinformatics. 2012;28(23):3131–3.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  107. 107.

    Servant N, Varoquaux N, Heard E, Barillot E, Vert JP. Effective normalization for copy number variation in Hi-C data. BMC Bioinformatics. 2018;19(1):313.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  108. 108.

    Stansfield JC, Cresswell KG, Vladimirov VI, Dozmorov MG. HiCcompare: an R-package for joint normalization and comparison of HI-C datasets. BMC Bioinformatics. 2018;19(1):279.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  109. 109.

    Serra F, Di Stefano M, Spill YG, Cuartero Y, Goodstadt M, Baù D, Marti-Renom MA. Restraint-based three-dimensional modeling of genomes and genomic domains. FEBS Lett. 2015;589(20):2987–95.

    CAS  PubMed  Article  Google Scholar 

  110. 110.

    Baù D, Marti-Renom MA. Genome structure determination via 3C-based data integration by the integrative modeling platform. Methods. 2012;58(3):300–6.

    PubMed  Article  CAS  Google Scholar 

  111. 111.

    Brunger AT. Version 1.2 of the Crystallography and NMR system. Nat Protoc. 2007;2(11):2728.

    CAS  PubMed  Article  Google Scholar 

  112. 112.

    Brünger AT, Adams PD, Clore GM, DeLano WL, Gros P, Grosse-Kunstleve RW, Jiang JS, Kuszewski J, Nilges M, Pannu NS, Read RJ. Crystallography & NMR system: a new software suite for macromolecular structure determination. Acta Crystallogr Sect D. 1998;54(5):905–21.

    Article  Google Scholar 

  113. 113.

    Duchi J, Hazan E, Singer Y. Adaptive subgradient methods for online learning and stochastic optimization. J Mach Learn Res. 2011;12(Jul):2121–59.

    Google Scholar 

  114. 114.

    Rieping W, Habeck M, Nilges M. Inferential structure determination. Science. 2005;309(5732):303–6.

    CAS  PubMed  Article  Google Scholar 

  115. 115.

    Mishra B, Meyer G, Sepulchre R. Low-rank optimization for distance matrix completion. In: 50th IEEE Conference on Decision and Control and European Control Conference 2011 Dec 12: IEEE; 2011. p. 4455–60.

  116. 116.

    Kruskal JB. Multidimensional scaling by optimizing goodness of fit to a nonmetric hypothesis. Psychometrika. 1964;29(1):1–27.

    Article  Google Scholar 

  117. 117.

    Shepard RN. The analysis of proximities: multidimensional scaling with an unknown distance function. I Psychometrika. 1962;27(2):125–40.

    Article  Google Scholar 

  118. 118.

    Ben-Elazar S, Yakhini Z, Yanai I. Spatial localization of co-regulated genes exceeds genomic gene clustering in the Saccharomyces cerevisiae genome. Nucleic Acids Res. 2013;41(4):2191–201.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  119. 119.

    Agarwal S, Wills J, Cayton L, Lanckriet G, Kriegman D, Belongie S. Generalized non-metric multidimensional scaling. In: Artificial Intelligence and Statistics; 2007. p. 11–8.

    Google Scholar 

  120. 120.

    Stevens TJ, Lando D, Basu S, Atkinson LP, Cao Y, Lee SF, Leeb M, Wohlfahrt KJ, Boucher W, O’Shaughnessy-Kirwan A, Cramard J. 3D structures of individual mammalian genomes studied by single-cell Hi-C. Nature. 2017;544(7648):59.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  121. 121.

    Nagano T, Lubling Y, Várnai C, Dudley C, Leung W, Baran Y, Cohen NM, Wingett S, Fraser P, Tanay A. Cell-cycle dynamics of chromosomal organization at single-cell resolution. Nature. 2017;547(7661):61.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  122. 122.

    Trussart M, Serra F, Baù D, Junier I, Serrano L, Marti-Renom MA. Assessing the limits of restraint-based 3D modeling of genomes and genomic domains. Nucleic Acids Res. 2015;43(7):3465–77.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  123. 123.

    Osborne CS, Chakalova L, Brown KE, Carter D, Horton A, Debrand E, Goyenechea B, Mitchell JA, Lopes S, Reik W, Fraser P. Active genes dynamically colocalize to shared sites of ongoing transcription. Nat Genet. 2004;36(10):1065.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  124. 124.

    Gozzetti A, Le Beau MM. Fluorescence in situ hybridization: uses and limitations. In Seminars in hematology 2000 Oct 1 (Vol. 37, No. 4, pp. 320–33). WB Saunders.

  125. 125.

    Ferrai C, de Castro IJ, Lavitas L, Chotalia M, Pombo A. Gene positioning. Cold Spring Harb Perspect Biol. 2010;2:a000588.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  126. 126.

    Holwerda S, De Laat W. Chromatin loops, gene positioning, and gene expression. Front Genet. 2012;3:217.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  127. 127.

    Geyer PK, Vitalini MW, Wallrath LL. Nuclear organization: taking a position on gene expression. Curr Opin Cell Biol. 2011;23(3):354–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  128. 128.

    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(6):1239–49.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  129. 129.

    Van Steensel B, Dekker J. Genomics tools for unraveling chromosome architecture. Nat Biotechnol. 2010;28(10):1089.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  130. 130.

    Hell SW. Microscopy and its focal switch. Nat Methods. 2009;6(1):24.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  131. 131.

    Gaietta G, Deerinck TJ, Adams SR, Bouwer J, Tour O, Laird DW, Sosinsky GE, Tsien RY, Ellisman MH. Multicolor and electron microscopic imaging of connexin trafficking. Science. 2002;296(5567):503–7.

    CAS  PubMed  Article  Google Scholar 

  132. 132.

    Rust MJ, Bates M, Zhuang X. Sub-diffraction-limit imaging by stochastic optical reconstruction microscopy (STORM). Nat Methods. 2006;3(10):793.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  133. 133.

    Tam J, Merino D. Stochastic optical reconstruction microscopy (STORM) in comparison with stimulated emission depletion (STED) and other imaging methods. J Neurochem. 2015;135(4):643–58.

    CAS  PubMed  Article  Google Scholar 

  134. 134.

    Daban JR. Electron microscopy and atomic force microscopy studies of chromatin and metaphase chromosome structure. Micron. 2011;42(8):733–50.

    CAS  PubMed  Article  Google Scholar 

  135. 135.

    Hess ST, Girirajan TP, Mason MD. Ultra-high resolution imaging by fluorescence photoactivation localization microscopy. Biophys J. 2006;91(11):4258–72.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  136. 136.

    Ricci MA, Manzo C, García-Parajo MF, Lakadamyali M, Cosma MP. Chromatin fibers are formed by heterogeneous groups of nucleosomes in vivo. Cell. 2015;160(6):1145–58.

    CAS  PubMed  Article  Google Scholar 

  137. 137.

    Ploem JS, Tanke HJ. Introduction to fluorescence microscopy; 1987.

    Google Scholar 

  138. 138.

    Ghiran IC. Introduction to fluorescence microscopy. In: Light microscopy. Totowa: Humana Press; 2011. p. 93–136.

    Chapter  Google Scholar 

  139. 139.

    Lindon JC, Tranter GE, Koppenaal D. Encyclopedia of spectroscopy and spectrometry. London: Academic Press; 2016.

    Google Scholar 

  140. 140.

    Haines AM, Tobe SS, Kobus HJ, Linacre A. Properties of nucleic acid staining dyes used in gel electrophoresis. Electrophoresis. 2015;36(6):941–4.

    CAS  PubMed  Article  Google Scholar 

  141. 141.

    Singer VL, Lawlor TE, Yue S. Comparison of SYBR® Green I nucleic acid gel stain mutagenicity and ethidium bromide mutagenicity in the salmonella/mammalian microsome reverse mutation assay (Ames test). Mutat Res Genet Toxicol Environ Mutagen. 1999;439(1):37–47.

    CAS  Article  Google Scholar 

  142. 142.

    Suzuki T, Fujikura K, Higashiyama T, Takata K. DNA staining for fluorescence and laser confocal microscopy. J Histochem Cytochem. 1997;45(1):49–53.

    CAS  PubMed  Article  Google Scholar 

  143. 143.

    Axelrod D, Koppel DE, Schlessinger J, Elson E, Webb WW. Mobility measurement by analysis of fluorescence photobleaching recovery kinetics. Biophys J. 1976;16(9):1055–69.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  144. 144.

    Sprague BL, Pego RL, Stavreva DA, McNally JG. Analysis of binding reactions by fluorescence recovery after photobleaching. Biophys J. 2004;86(6):3473–95.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  145. 145.

    Wüstner D, Solanko LM, Lund FW, Sage D, Schroll HJ, Lomholt MA. Quantitative fluorescence loss in photobleaching for analysis of protein transport and aggregation. BMC Bioinformatics. 2012;13(1):296.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  146. 146.

    Ishikawa-Ankerhold HC, Ankerhold R, Drummen GP. Advanced fluorescence microscopy techniques—Frap, Flip, Flap, Fret and flim. Molecules. 2012;17(4):4047–132.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  147. 147.

    Ratan ZA, Zaman SB, Mehta V, Haidere MF, Runa NJ, Akter N. Application of fluorescence in situ hybridization (FISH) technique for the detection of genetic aberration in medical science. Cureus. 2017;9(6):e1325.

    PubMed  PubMed Central  Google Scholar 

  148. 148.

    Cremer T, Cremer C, Schneider T, Baumann H, Hens L, Kirsch-Volders M. Analysis of chromosome positions in the interphase nucleus of Chinese hamster cells by laser-UV-microirradiation experiments. Hum Genet. 1982;62(3):201–9.

    CAS  PubMed  Article  Google Scholar 

  149. 149.

    Branco MR, Pombo A. Intermingling of chromosome territories in interphase suggests role in translocations and transcription-dependent associations. PLoS Biol. 2006;4(5):e138.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  150. 150.

    Mahy NL, Perry PE, Bickmore WA. Gene density and transcription influence the localization of chromatin outside of chromosome territories detectable by FISH. J Cell Biol. 2002;159(5):753–63.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  151. 151.

    Chambeyron S, Bickmore WA. Chromatin decondensation and nuclear reorganization of the HoxB locus upon induction of transcription. Genes Dev. 2004;18(10):1119–30.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  152. 152.

    Shopland LS, Lynch CR, Peterson KA, Thornton K, Kepper N, von Hase J, Stein S, Vincent S, Molloy KR, Kreth G, Cremer C. Folding and organization of a contiguous chromosome region according to the gene distribution pattern in primary genomic sequence. J Cell Biol. 2006;174(1):27–38.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  153. 153.

    Brown JM, Green J, das Neves RP, Wallace HA, Smith AJ, Hughes J, Gray N, Taylor S, Wood WG, Higgs DR, Iborra FJ. Association between active genes occurs at nuclear speckles and is modulated by chromatin environment. J Cell Biol. 2008;182(6):1083–97.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  154. 154.

    Shachar S, Voss TC, Pegoraro G, Sciascia N, Misteli T. Identification of gene positioning factors using high-throughput imaging mapping. Cell. 2015;162(4):911–23.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  155. 155.

    Batson PE, Dellby N, Krivanek OL. Sub-ångstrom resolution using aberration corrected electron optics. Nature. 2002;418(6898):617.

    CAS  PubMed  Article  Google Scholar 

  156. 156.

    Erni R, Rossell MD, Kisielowski C, Dahmen U. Atomic-resolution imaging with a sub-50-pm electron probe. Phys Rev Lett. 2009;102(9):096101.

    PubMed  Article  CAS  Google Scholar 

  157. 157.

    Crewe AV, Isaacson M, Johnson D. A simple scanning electron microscope. Rev Sci Instrum. 1969;40(2):241–6.

    Article  Google Scholar 

  158. 158.

    Scherzer O. The theoretical resolution limit of the electron microscope. J Appl Phys. 1949;20(1):20–9.

    CAS  Article  Google Scholar 

  159. 159.

    Haider M, Uhlemann S, Schwan E, Rose H, Kabius B, Urban K. Electron microscopy image enhanced. Nature. 1998;392(6678):768.

    CAS  Article  Google Scholar 

  160. 160.

    Callaway E. The revolution will not be crystallized: a new method sweeps through structural biology. Nature News. 2015;525(7568):172.

    CAS  Article  Google Scholar 

  161. 161.

    Glaeser RM. How good can cryo-EM become? Nat Methods. 2015;13(1):28.

    Article  CAS  Google Scholar 

  162. 162.

    Iacovache I, De Carlo S, Cirauqui N, Dal Peraro M, Van Der Goot FG, Zuber B. Cryo-EM structure of aerolysin variants reveals a novel protein fold and the pore-formation process. Nat Commun. 2016;7:12062.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  163. 163.

    Ou HD, Phan S, Deerinck TJ, Thor A, Ellisman MH, O’shea CC. ChromEMT: Visualizing 3D chromatin structure and compaction in interphase and mitotic cells. Science. 2017;357(6349):eaag0025.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  164. 164.

    Bouwman BA, de Laat W. Architectural hallmarks of the pluripotent genome. FEBS Lett. 2015;589(20):2905–13.

    CAS  PubMed  Article  Google Scholar 

  165. 165.

    Felsenfeld G, Groudine M. Controlling the double helix. Nature. 2003;421(6921):448.

    PubMed  Article  CAS  Google Scholar 

  166. 166.

    Chubb JR, Boyle S, Perry P, Bickmore WA. Chromatin motion is constrained by association with nuclear compartments in human cells. Curr Biol. 2002;12(6):439–45.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  167. 167.

    Walter J, Schermelleh L, Cremer M, Tashiro S, Cremer T. Chromosome order in HeLa cells changes during mitosis and early G1, but is stably maintained during subsequent interphase stages. J Cell Biol. 2003;160(5):685–97.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  168. 168.

    Ramani V, Shendure J, Duan Z. Understanding spatial genome organization: methods and insights. Genomics Proteomics Bioinformatics. 2016;14(1):7–20.

    PubMed  PubMed Central  Article  Google Scholar 

  169. 169.

    Bonev B, Cavalli G. Organization and function of the 3D genome. Nature Reviews Genetics. 2016;17(11):661.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  170. 170.

    Nowotny J, Wells A, Oluwadare O, Xu L, Cao R, Trieu T, He C, Cheng J. GMOL: an interactive tool for 3D genome structure visualization. Scientific Reports. 2016;6:20802.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  171. 171.

    Djekidel MN, Wang M, Zhang MQ, Gao J. HiC-3DViewer: a new tool to visualize Hi-C data in 3D space. Quantitative Biology. 2017;5(2):183–90.

    CAS  Article  Google Scholar 

  172. 172.

    Li R, Liu Y, Li T, Li C. 3Disease Browser: a web server for integrating 3D genome and disease-associated chromosome rearrangement data. Scientific Reports. 2016;6:34651.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  173. 173.

    Asbury TM, Mitman M, Tang J, Zheng WJ. Genome3D: a viewer-model framework for integrating and visualizing multi-scale epigenomic information within a three-dimensional genome. BMC Bioinformatics. 2010;11(1):444.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  174. 174.

    Tang B, Li F, Li J, Zhao W, Zhang Z. Delta: a new web-based 3D genome visualization and analysis platform. Bioinformatics. 2017;34(8):1409–10.

    Article  CAS  Google Scholar 

Download references


Not applicable.


This work was supported by the National Science Foundation (NSF) CAREER award (grant no: DBI1149224) to JC.

Availability of Data and Materials

Not applicable.

Author information




OO and JC designed the manuscript outlines. OO drafted the manuscript, and MH and JC revised it. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jianlin Cheng.

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 (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Oluwadare, O., Highsmith, M. & Cheng, J. An Overview of Methods for Reconstructing 3-D Chromosome and Genome Structures from Hi-C Data. Biol Proced Online 21, 7 (2019).

Download citation


  • Hi-C
  • 3-D chromosome and genome structure
  • Chromosome conformation capture
  • 3-D genome
  • Contact-based modeling
  • Distance-based modeling
  • Optimization