Dataset Information

MrTADFinder: A network modularity based approach to identify topologically associating domains in multiple resolutions.

ABSTRACT: Genome-wide proximity ligation based assays such as Hi-C have revealed that eukaryotic genomes are organized into structural units called topologically associating domains (TADs). From a visual examination of the chromosomal contact map, however, it is clear that the organization of the domains is not simple or obvious. Instead, TADs exhibit various length scales and, in many cases, a nested arrangement. Here, by exploiting the resemblance between TADs in a chromosomal contact map and densely connected modules in a network, we formulate TAD identification as a network optimization problem and propose an algorithm, MrTADFinder, to identify TADs from intra-chromosomal contact maps. MrTADFinder is based on the network-science concept of modularity. A key component of it is deriving an appropriate background model for contacts in a random chain, by numerically solving a set of matrix equations. The background model preserves the observed coverage of each genomic bin as well as the distance dependence of the contact frequency for any pair of bins exhibited by the empirical map. Also, by introducing a tunable resolution parameter, MrTADFinder provides a self-consistent approach for identifying TADs at different length scales, hence the acronym "Mr" standing for Multiple Resolutions. We then apply MrTADFinder to various Hi-C datasets. The identified domain boundaries are marked by characteristic signatures in chromatin marks and transcription factors (TF) that are consistent with earlier work. Moreover, by calling TADs at different length scales, we observe that boundary signatures change with resolution, with different chromatin features having different characteristic length scales. Furthermore, we report an enrichment of HOT (high-occupancy target) regions near TAD boundaries and investigate the role of different TFs in determining boundaries at various resolutions. To further explore the interplay between TADs and epigenetic marks, as tumor mutational burden is known to be coupled to chromatin structure, we examine how somatic mutations are distributed across boundaries and find a clear stepwise pattern. Overall, MrTADFinder provides a novel computational framework to explore the multi-scale structures in Hi-C contact maps.

SUBMITTER: Yan KK

PROVIDER: S-EPMC5546724 | biostudies-literature | 2017 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

MrTADFinder: A network modularity based approach to identify topologically associating domains in multiple resolutions.

Yan Koon-Kiu KK Lou Shaoke S Gerstein Mark M

PLoS computational biology 20170724 7

Genome-wide proximity ligation based assays such as Hi-C have revealed that eukaryotic genomes are organized into structural units called topologically associating domains (TADs). From a visual examination of the chromosomal contact map, however, it is clear that the organization of the domains is not simple or obvious. Instead, TADs exhibit various length scales and, in many cases, a nested arrangement. Here, by exploiting the resemblance between TADs in a chromosomal contact map and densely co ...[more]

PMID: 28742097

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:Eukaryotic genomes are structurally organized via the formation of multiple loops that create gene expression regulatory units called topologically associating domains (TADs). Here we revealed the KSHV TAD structure at 500 bp resolution and constructed a 3D KSHV genomic structural model with 2 kb binning. The latent KSHV genome formed very similar genomic architectures in three different naturally infected PEL cell lines and in an experimentally infected epithelial cell line. The majority of the TAD boundaries were occupied by structural maintenance of chromosomes (SMC1) cohesin complex and CCCTC-binding factor (CTCF), and the KSHV transactivator was recruited to those sites during reactivation. Triggering KSHV gene expression decreased prewired genomic loops within the regulatory unit, while contacts extending outside of regulatory borders increased, leading to formation of a larger regulatory unit with a shift from repressive to active compartments (B to A). The 3D genomic structural model proposes that the immediate early promoter region is localized on the periphery of the 3D viral genome during latency, while highly inducible noncoding RNA regions moved toward the inner space of the structure, resembling the configuration of a "bird cage" during reactivation. The compartment-like properties of viral episomal chromatin structure and its reorganization during the transition from latency may help facilitate viral gene transcription. IMPORTANCE The 3D architecture of chromatin allows for efficient arrangement, expression, and replication of genetic material. The genomes of all organisms studied to date have been found to be organized through some form of tiered domain structures. However, the architectural framework of the genomes of large double-stranded DNA viruses such as the herpesvirus family has not been reported. Prior studies with Kaposi's sarcoma-associated herpesvirus (KSHV) have indicated that the viral chromatin shares many biological properties exhibited by the host cell genome, essentially behaving as a mini human chromosome. Thus, we hypothesized that the KSHV genome may be organized in a similar manner. In this report, we describe the domain structure of the latent and lytic KSHV genome at 500 bp resolution and present a 3D genomic structural model for KSHV under each condition. These results add new insights into the complex regulation of the viral life cycle.

Dataset Information

MrTADFinder: A network modularity based approach to identify topologically associating domains in multiple resolutions.

Publications

MrTADFinder: A network modularity based approach to identify topologically associating domains in multiple resolutions.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets