Dataset Information

Identifying topologically associating domains and subdomains by Gaussian Mixture model And Proportion test.

ABSTRACT: The spatial organization of the genome plays a critical role in regulating gene expression. Recent chromatin interaction mapping studies have revealed that topologically associating domains and subdomains are fundamental building blocks of the three-dimensional genome. Identifying such hierarchical structures is a critical step toward understanding the three-dimensional structure-function relationship of the genome. Existing computational algorithms lack statistical assessment of domain predictions and are computationally inefficient for high-resolution Hi-C data. We introduce the Gaussian Mixture model And Proportion test (GMAP) algorithm to address the above-mentioned challenges. Using simulated and experimental Hi-C data, we show that domains identified by GMAP are more consistent with multiple lines of supporting evidence than three state-of-the-art methods. Application of GMAP to normal and cancer cells reveals several unique features of subdomain boundary as compared to domain boundary, including its higher dynamics across cell types and enrichment for somatic mutations in cancer.Spatial organization of the genome plays a crucial role in regulating gene expression. Here the authors introduce GMAP, the Gaussian Mixture model And Proportion test, to identify topologically associating domains and subdomains in Hi-C data.

SUBMITTER: Yu W

PROVIDER: S-EPMC5599511 | biostudies-literature | 2017 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Identifying topologically associating domains and subdomains by Gaussian Mixture model And Proportion test.

Yu Wenbao W He Bing B Tan Kai K

Nature communications 20170914 1

The spatial organization of the genome plays a critical role in regulating gene expression. Recent chromatin interaction mapping studies have revealed that topologically associating domains and subdomains are fundamental building blocks of the three-dimensional genome. Identifying such hierarchical structures is a critical step toward understanding the three-dimensional structure-function relationship of the genome. Existing computational algorithms lack statistical assessment of domain predicti ...[more]

PMID: 28912419

Similar Datasets

Project description:Eukaryotic chromosomes replicate in a temporal order known as the replication-timing program. In mammals, replication timing is cell-type-specific with at least half the genome switching replication timing during development, primarily in units of 400-800 kilobases ('replication domains'), whose positions are preserved in different cell types, conserved between species, and appear to confine long-range effects of chromosome rearrangements. Early and late replication correlate, respectively, with open and closed three-dimensional chromatin compartments identified by high-resolution chromosome conformation capture (Hi-C), and, to a lesser extent, late replication correlates with lamina-associated domains (LADs). Recent Hi-C mapping has unveiled substructure within chromatin compartments called topologically associating domains (TADs) that are largely conserved in their positions between cell types and are similar in size to replication domains. However, TADs can be further sub-stratified into smaller domains, challenging the significance of structures at any particular scale. Moreover, attempts to reconcile TADs and LADs to replication-timing data have not revealed a common, underlying domain structure. Here we localize boundaries of replication domains to the early-replicating border of replication-timing transitions and map their positions in 18 human and 13 mouse cell types. We demonstrate that, collectively, replication domain boundaries share a near one-to-one correlation with TAD boundaries, whereas within a cell type, adjacent TADs that replicate at similar times obscure replication domain boundaries, largely accounting for the previously reported lack of alignment. Moreover, cell-type-specific replication timing of TADs partitions the genome into two large-scale sub-nuclear compartments revealing that replication-timing transitions are indistinguishable from late-replicating regions in chromatin composition and lamina association and accounting for the reduced correlation of replication timing to LADs and heterochromatin. Our results reconcile cell-type-specific sub-nuclear compartmentalization and replication timing with developmentally stable structural domains and offer a unified model for large-scale chromosome structure and function.

Project description:Experiments based on chromosome conformation capture have shown that mammalian genomes are partitioned into topologically associating domains (TADs), within which the chromatin fiber preferentially interacts. TADs may provide three-dimensional scaffolds allowing genes to contact their appropriate distal regulatory DNA sequences (e.g., enhancers) and thus to be properly regulated. Understanding the cell-to-cell and temporal variability of the chromatin fiber within TADs, and what determines them, is thus of great importance to better understand transcriptional regulation. We recently described an equilibrium polymer model that can accurately predict cell-to-cell variation of chromosome conformation within single TADs, from chromosome conformation capture-based data. Here we further analyze the conformational and energetic properties of our model. We show that the chromatin fiber within TADs can easily fluctuate between several conformational states, which are hierarchically organized and are not separated by important free energy barriers, and that this is facilitated by the fact that the chromatin fiber within TADs is close to the onset of the coil-globule transition. We further show that in this dynamic state the properties of the chromatin fiber, and its contact probabilities in particular, are determined in a nontrivial manner not only by site-specific interactions between strongly interacting loci along the fiber, but also by nonlocal correlations between pairs of contacts. Finally, we use live-cell experiments to measure the dynamics of the chromatin fiber in mouse embryonic stem cells, in combination with dynamical simulations, and predict that conformational changes within one TAD are likely to occur on timescales that are much shorter than the duration of one cell cycle. This suggests that genes and their regulatory elements may come together and disassociate several times during a cell cycle. These results have important implications for transcriptional regulation as they support the concept of highly dynamic interactions driven by a complex interplay between site-specific interactions and the intrinsic biophysical properties of the chromatin fiber.

Project description:BACKGROUND:Topologically associating domains (TADs) are genomic regions with varying lengths. The interactions within TADs are more frequent than those between different TADs. TADs or sub-TADs are considered the structural and functional units of the mammalian genomes. Although TADs are important for understanding how genomes function, we have limited knowledge about their 3D structural properties. RESULTS:In this study, we designed and benchmarked three metrics for capturing the three-dimensional and two-dimensional structural signatures of TADs, which can help better understand TADs' structural properties and the relationships between structural properties and genetic and epigenetic features. The first metric for capturing 3D structural properties is radius of gyration, which in this study is used to measure the spatial compactness of TADs. The mass value of each DNA bead in a 3D structure is novelly defined as one or more genetic or epigenetic feature(s). The second metric is folding degree. The last metric is exponent parameter, which is used to capture the 2D structural properties based on TADs' Hi-C contact matrices. In general, we observed significant correlations between the three metrics and the genetic and epigenetic features. We made the same observations when using H3K4me3, transcription start sites, and RNA polymerase II to represent the mass value in the modified radius-of-gyration metric. Moreover, we have found that the TADs in the clusters of depleted chromatin states apparently correspond to smaller exponent parameters and larger radius of gyrations. In addition, a new objective function of multidimensional scaling for modelling chromatin or TADs 3D structures was designed and benchmarked, which can handle the DNA bead-pairs with zero Hi-C contact values. CONCLUSIONS:The web server for reconstructing chromatin 3D structures using multiple different objective functions and the related source code are publicly available at http://dna.cs.miami.edu/3DChrom/.

Dataset Information

Identifying topologically associating domains and subdomains by Gaussian Mixture model And Proportion test.

Publications

Identifying topologically associating domains and subdomains by Gaussian Mixture model And Proportion test.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets