Three-dimensional Genome Organization Research Articles

Noncoding genetic variation in GATA3 increases acute lymphoblastic leukemia risk through local and global changes in chromatin conformation.

Inherited noncoding genetic variants confer significant disease susceptibility to childhood acute lymphoblastic leukemia (ALL) but the molecular processes linking germline polymorphisms with somatic lesions in this cancer are poorly understood. Through targeted sequencing in 5,008 patients, we identified a key regulatory germline variant in GATA3 associated with Philadelphia chromosome-like ALL (Ph-like ALL). Using CRISPR-Cas9 editing and samples from patients with Ph-like ALL, we showed that this variant activated a strong enhancer that upregulated GATA3 transcription. This, in turn, reshaped global chromatin accessibility and three-dimensional genome organization, including regions proximal to the ALL oncogene CRLF2. Finally, we showed that GATA3 directly regulated CRLF2 and potentiated the JAK-STAT oncogenic effects during leukemogenesis. Taken together, we provide evidence for a distinct mechanism by which a germline noncoding variant contributes to oncogene activation, epigenetic regulation and three-dimensional genome reprogramming.

CRISPR and biochemical screens identify MAZ as a cofactor in CTCF-mediated insulation at Hox clusters

CCCTC-binding factor (CTCF) is critical to three-dimensional genome organization. Upon differentiation, CTCF insulates active and repressed genes within Hox gene clusters. We conducted a genome-wide CRISPR knockout (KO) screen to identify genes required for CTCF-boundary activity at the HoxA cluster, complemented by biochemical approaches. Among the candidates, we identified Myc-associated zinc-finger protein (MAZ) as a cofactor in CTCF insulation. MAZ colocalizes with CTCF at chromatin borders and, similar to CTCF, interacts with the cohesin subunit RAD21. MAZ KO disrupts gene expression and local contacts within topologically associating domains. Similar to CTCF motif deletions, MAZ motif deletions lead to derepression of posterior Hox genes immediately after CTCF boundaries upon differentiation, giving rise to homeotic transformations in mouse. Thus, MAZ is a factor contributing to appropriate insulation, gene expression and genomic architecture during development.

Near-atomistic modeling reconciles difference between irregular and regular chromatin

Three-dimensional genome organization lays the foundation for biological processes involving gene expression and epigenetic regulation. Nevertheless, it is unclear how chromatin is structured at the fundamental level. There is a heated debate over the existence of chromatin fibril structure and its regulation by multiple chromatin factors. Controversy also remains on the sol-gel properties of chromatin subject to different environmental conditions. Here, building upon our recently developed near-atomistic chromatin model, we combined exhaustive molecular dynamics simulations with deep-learning techniques to study the detailed folding landscape of a string of nucleosomes.

Contribution of 3D genome topological domains to genetic risk of cancers: a genome-wide computational study

BackgroundGenome-wide association studies have identified statistical associations between various diseases, including cancers, and a large number of single-nucleotide polymorphisms (SNPs). However, they provide no direct explanation of the mechanisms underlying the association. Based on the recent discovery that changes in three-dimensional genome organization may have functional consequences on gene regulation favoring diseases, we investigated systematically the genome-wide distribution of disease-associated SNPs with respect to a specific feature of 3D genome organization: topologically associating domains (TADs) and their borders.ResultsFor each of 449 diseases, we tested whether the associated SNPs are present in TAD borders more often than observed by chance, where chance (i.e., the null model in statistical terms) corresponds to the same number of pointwise loci drawn at random either in the entire genome, or in the entire set of disease-associated SNPs listed in the GWAS catalog. Our analysis shows that a fraction of diseases displays such a preferential localization of their risk loci. Moreover, cancers are relatively more frequent among these diseases, and this predominance is generally enhanced when considering only intergenic SNPs. The structure of SNP-based diseasome networks confirms that localization of risk loci in TAD borders differs between cancers and non-cancer diseases. Furthermore, different TAD border enrichments are observed in embryonic stem cells and differentiated cells, consistent with changes in topological domains along embryogenesis and delineating their contribution to disease risk.ConclusionsOur results suggest that, for certain diseases, part of the genetic risk lies in a local genetic variation affecting the genome partitioning in topologically insulated domains. Investigating this possible contribution to genetic risk is particularly relevant in cancers. This study thus opens a way of interpreting genome-wide association studies, by distinguishing two types of disease-associated SNPs: one with an effect on an individual gene, the other acting in interplay with 3D genome organization.

Deciphering the Role of 3D Genome Organization in Breast Cancer Susceptibility.

Cancer risk by environmental exposure is modulated by an individual’s genetics and age at exposure. This age-specific period of susceptibility is referred to as the “Window of Susceptibility” (WOS). Rats have a similar WOS for developing breast cancer. A previous study in rat identified an age-specific long-range regulatory interaction for the cancer gene, Pappa, that is associated with breast cancer susceptibility. However, the global role of three-dimensional genome organization and downstream gene expression programs in the WOS is not known. Therefore, we generated Hi-C and RNA-seq data in rat mammary epithelial cells within and outside the WOS. To systematically identify higher-order changes in 3D genome organization, we developed NE-MVNMF that combines network enhancement followed by multitask non-negative matrix factorization. We examined three-dimensional genome organization dynamics at the level of individual loops as well as higher-order domains. Differential chromatin interactions tend to be associated with differentially up-regulated genes with the WOS and recapitulate several human SNP-gene interactions associated with breast cancer susceptibility. Our approach identified genomic blocks of regions with greater overall differences in contact count between the two time points when the cluster assignments change and identified genes and pathways implicated in early carcinogenesis and cancer treatment. Our results suggest that WOS-specific changes in 3D genome organization are linked to transcriptional changes that may influence susceptibility to breast cancer.

Assessing Specific Networks of Chromatin Interactions with HiChIP.

The introduction of chromosome conformation capture (3C)-based technologies coupled with next-generation sequencing have significantly advanced our understanding of how the genetic material is organized within the eukaryotic nucleus. Three-dimensional (3D) genomic organization occurs at hierarchical levels, ranging from chromosome territories and subnuclear compartments to smaller self-associated domains and fine-scale chromatin interactions. The latter can be further categorized into different subtypes, such as structural or regulatory, based either on their presumed functionality and/or the factors that mediate their formation. Various enrichment strategies coupled with 3C-based technologies have been developed to prospectively isolate and quantify chromatin interactions around regions occupied by specific proteins or marks of interest. These approaches not only enable high-resolution characterization of the selected chromatin contacts at a cost-effective manner, but also offer important biological insights into their organizational principles and regulatory function. In this chapter, we will focus on the recently developed HiChIP technology with an emphasis on the discovery of putative active enhancers and promoter interactions in cell types of interest. We will describe the specific steps for designing, performing and analyzing successful HiChIP experiments as well as important limitations and considerations.

Three-Dimensional Genome Organization in Breast and Gynecological Cancers: How Chromatin Folding Influences Tumorigenic Transcriptional Programs.

A growing body of research on the transcriptome and cancer genome has demonstrated that many gynecological tumor-specific gene mutations are located in cis-regulatory elements. Through chromosomal looping, cis-regulatory elements interact which each other to control gene expression by bringing distant regulatory elements, such as enhancers and insulators, into close proximity with promoters. It is well known that chromatin connections may be disrupted in cancer cells, promoting transcriptional dysregulation and the expression of abnormal tumor suppressor genes and oncogenes. In this review, we examine the roles of alterations in 3D chromatin interactions. This includes changes in CTCF protein function, cancer-risk single nucleotide polymorphisms, viral integration, and hormonal response as part of the mechanisms that lead to the acquisition of enhancers or super-enhancers. The translocation of existing enhancers, as well as enhancer loss or acquisition of insulator elements that interact with gene promoters, is also revised. Remarkably, similar processes that modify 3D chromatin contacts in gene promoters may also influence the expression of non-coding RNAs, such as long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), which have emerged as key regulators of gene expression in a variety of cancers, including gynecological malignancies.

Three-dimensional telomere profiles in papillary thyroid cancer variants: A pilot study.

Besides the two main histologic types of papillary thyroid carcinoma (PTC), the classical PTC (CL-PTC) and the follicular variant PTC (FV-PTC), several other variants are described. The encapsulated FV-PTC variant was recently reclassified as non-invasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP) due to its similarities to benign lesions. Specific molecular signatures, however, are still unavailable. It is well known that improper DNA repair of dysfunctional telomeres may cause telomere-related genome instability. The mechanisms involved in the damaged telomere repair processing may lead to detrimental outcomes, altering the three-dimensional (3D) nuclear telomere and genome organization in cancer cells. This pilot study aimed to evaluate whether specific 3D nuclear telomere architecture might characterize NIFTP, potentially distinguishing it from other PTC histologic variants. Our findings demonstrate that 3D telomere profiles of CL-PTC and FV-PTC were different from NIFTP and that NIFTP more closely resembles follicular thyroid adenoma (FTA). NIFTP has longer telomeres than CL-PTC and FV-PTC samples, and the telomere length of NIFTP overlaps with that of the FTA histotype. In contrast, there was no association between BRAF expression and telomere length in all tested samples. These preliminary findings reinforce the view that NIFTP is closer to non-malignant thyroid nodules and confirm that PTC features short telomeres.

Reorganization of the 3D Genome Pinpoints Noncoding Drivers of Primary Prostate Tumors.

Prostate cancer is a heterogeneous disease whose progression is linked to genome instability. However, the impact of this instability on the noncoding genome and its three-dimensional organization to aid progression is unclear. Using primary benign and tumor tissue, we find a high concordance in higher-order three-dimensional genome organization. This concordance argues for constraints to the topology of prostate tumor genomes. Nonetheless, we identified changes in focal chromatin interactions, typical of loops bridging noncoding cis-regulatory elements, and showed how structural variants can induce these changes to guide cis-regulatory element hijacking. Such events resulted in opposing differential expression of genes found at antipodes of rearrangements. Collectively, these results argue that changes to focal chromatin interactions, as opposed to higher-order genome organization, allow for aberrant gene regulation and are repeatedly mediated by structural variants in primary prostate cancer. SIGNIFICANCE: This work showcases how the noncoding genome can be hijacked by focal insults to its three-dimensional organization that contribute to prostate cancer oncogenesis.

Defective chromatin architectures in embryonic stem cells derived from somatic cell nuclear transfer impair their differentiation potentials

Nuclear transfer embryonic stem cells (ntESCs) hold enormous promise for individual-specific regenerative medicine. However, the chromatin states of ntESCs remain poorly characterized. In this study, we employed ATAC-seq and Hi-C techniques to explore the chromatin accessibility and three-dimensional (3D) genome organization of ntESCs. The results show that the chromatin accessibility and genome structures of somatic cells are re-arranged to ESC-like states overall in ntESCs, including compartments, topologically associating domains (TADs) and chromatin loops. However, compared to fertilized ESCs (fESCs), ntESCs show some abnormal openness and structures that have not been reprogrammed completely, which impair the differentiation potential of ntESCs. The histone modification H3K9me3 may be involved in abnormal structures in ntESCs, including incorrect compartment switches and incomplete TAD rebuilding. Moreover, ntESCs and iPSCs show high similarity in 3D genome structures, while a few differences are detected due to different somatic cell origins and reprogramming mechanisms. Through systematic analyses, our study provides a global view of chromatin accessibility and 3D genome organization in ntESCs, which can further facilitate the understanding of the similarities and differences between ntESCs and fESCs.

Combinatorial Genetic Uncovers DOCK1/RAC2 As Specific Targets for the Treatment of NPM1;Cohesin Mutated AML

Acute myeloid leukemia (AML) is an aggressive malignancy of the blood and bone marrow resulting from the accumulation of multiple serially acquired mutations. The mutational complexity of the disease makes AML difficult to treat, contributing to a low 5-year survival rate. A better understanding of how different mutations interact with one another to influence disease characteristics is critical and may result in the development of targeted therapies to improve patient outcome. The most common recurrent somatic mutation in AML affects the gene NPM1 and occurs in 25-30% of patients. This mutation results in the aberrant cytoplasmic localization of the protein and is termed NPM1cA (Falini et al. 2005, Cancer Genome Atlas Research Network 2013). NPM1cA is considered to be a driver of AML, however Npm1 cA/+ mice only develop disease after a long latency (median 18 months), suggesting that other cooperating mutations are required for AML development (Vassiliou et al. 2011).Approximately 50% of patients with an NPM1cA mutation also harbor a mutation in one of four members of the cohesin complex (STAG2, SMC3, SMC1A, and RAD21). Mutations in cohesin genes are mutually exclusive and result in haploinsufficiency of the complex. As cohesin is known to regulate gene expression by facilitating promoter-enhancer interactions and three-dimensional genome organization, we wished to determine how cohesin mutation influences AML biology and gene expression in the presence of NPM1cA. We utilized the inducible Npm1 cAflox/+ and Smc3 flox/+ mouse models to examine this genetic interaction. We and others have shown that cohesin mutations result in enhanced hematopoietic stem and progenitor cell (HSPC) self-renewal (Mazumdar et al. 2015, Viny et al. 2015, Mullenders et al. 2015, Galeev et al. 2016, Fisher et al. 2017). Consistent with this, Npm1 cA/+;Smc3 Δ/+ HSPCs show enhanced self-renewal in vitro over HSPCs harboring either single mutation. Despite a shared role in self-renewal, Npm1 cA/+;Smc3 Δ/+mice developed AML with similar latency as Npm1 cA/+ mice. Interestingly, however, Npm1 cA/+;Smc3 Δ/+ HSPCs exhibited dysregulation of a unique set of genes compared to cells from either single mutant, suggesting that the Npm1 cA;cohesin mutational combination uniquely alters the transcriptional environment. Further, Npm1 cA/+;Smc3 Δ/+ leukemias had a completely unique spectrum of acquired mutations compared to Npm1 cA/+ leukemias, suggesting that the addition of Smc3 haploinsufficiency to Npm1 cA/+alters AML evolution and that different driver mutations result in the accumulation of very different somatic mutations.Among the uniquely upregulated genes in Npm1 cA/+;Smc3 Δ/+ HSPCs is Dock1, a guanine nucleotide exchange factor (GEF) for Rac1/2. As high Dock1 expression has been associated with low overall and disease-free survival in multiple cohorts of AML patients (Lee et al. 2017, Zhang et al. 2019), we hypothesized that Dock1 would be a novel target for the treatment of Npm1cA; cohesin mut leukemias. Consistent with this hypothesis, We found that knockdown of Dock1 resulted in decreased growth and adhesion and increased apoptosis in an Npm1 cA/+;Smc3 Δ/+leukemic line, but not in an Npm1 cAsingle mutant line. Higher Rac activity was also observed in Npm1 cA/+;Smc3 Δ/+ vs. Npm1 cA/+ leukemic lines. We found that DN Rac2 specifically impacted both the growth and apoptosis of an Npm1 cA/+;Smc3 Δ/+line, suggesting that Dock1 functions primarily through Rac2 to regulate survival. Importantly, the Dock1/Rac pathway is specifically targetable in Npm1 cA/+;Smc3 Δ/+ AMLs in vitro and in vivo, as knockdown of Dock1 resulted in prolonged latency in a Npm1 cA/+;Smc3 Δ/+ transplant model and slowed the growth and enhanced apoptosis of an Npm1 cA/+;Smc3 Δ/+, but not an Npm1 cA/+, leukemic line. Further, small molecule inhibitors of Dock and Rac had similar effects. Our results suggest that Dock1/Rac2 represent unique targets for the treatment of patients harboring the NPM1cA;cohesin mutational combination and provide validity to the concept that combinatorial genetics can uncover novel precision oncology targets. DisclosuresVassiliou: Kymab Ltd: Divested equity in a private or publicly-traded company in the past 24 months; STRM.BIO: Consultancy; Astrazeneca: Consultancy.

Three-dimensional genome organization via triplex-forming RNAs.

An increasing number of long noncoding RNAs (lncRNAs) have been proposed to act as nuclear organization factors during interphase. Direct RNA-DNA interactions can be achieved by the formation of triplex helix structures where a single-stranded RNA molecule hybridizes by complementarity into the major groove of double-stranded DNA. However, whether and how these direct RNA-DNA associations influence genome structure in interphase chromosomes remain poorly understood. Here we theorize that RNA organizes the genome in space via a triplex-forming mechanism. To test this theory, we apply a computational modeling approach of chromosomes that combines restraint-based modeling with polymer physics. Our models suggest that colocalization of triplex hotspots targeted by lncRNAs could contribute to large-scale chromosome compartmentalization cooperating, rather than competing, with architectural transcription factors such as CTCF.

Tcf1 and Lef1 provide constant supervision to mature\xa0CD8+ T cell identity and function by organizing genomic architecture

T cell identity is established during thymic development, but how it is maintained in the periphery remains unknown. Here we show that ablating Tcf1 and Lef1 transcription factors in mature CD8+ T cells aberrantly induces genes from non-T cell lineages. Using high-throughput chromosome-conformation-capture sequencing, we demonstrate that Tcf1/Lef1 are important for maintaining three-dimensional genome organization at multiple scales in CD8+ T cells. Comprehensive network analyses coupled with genome-wide profiling of chromatin accessibility and Tcf1 occupancy show the direct impact of Tcf1/Lef1 on the T cell genome is to promote formation of extensively interconnected hubs through enforcing chromatin interaction and accessibility. The integrative mechanisms utilized by Tcf1/Lef1 underlie activation of T cell identity genes and repression of non-T lineage genes, conferring fine control of various T cell functionalities. These findings suggest that Tcf1/Lef1 control global genome organization and help form intricate chromatin-interacting hubs to facilitate promoter-enhancer/silencer contact, hence providing constant supervision of CD8+ T cell identity and function.

BRD4 orchestrates genome folding to promote neural crest differentiation.

Higher-order chromatin structure regulates gene expression, and mutations in proteins mediating genome folding underlie developmental disorders known as cohesinopathies. However, the relationship between three-dimensional genome organization and embryonic development remains unclear. Here we define a role for bromodomain-containing protein 4 (BRD4) in genome folding, and leverage it to understand the importance of genome folding in neural crest progenitor differentiation. Brd4 deletion in neural crest results in cohesinopathy-like phenotypes. BRD4 interacts with NIPBL, a cohesin agonist, and BRD4 depletion or loss of the BRD4-NIPBL interaction reduces NIPBL occupancy, suggesting that BRD4 stabilizes NIPBL on chromatin. Chromatin interaction mapping and imaging experiments demonstrate that BRD4 depletion results in compromised genome folding and loop extrusion. Finally, mutation of individual BRD4 amino acids that mediate an interaction with NIPBL impedes neural crest differentiation into smooth muscle. Remarkably, loss of WAPL, a cohesin antagonist, rescues attenuated smooth muscle differentiation resulting from BRD4 loss. Collectively, our data reveal that BRD4 choreographs genome folding and illustrates the relevance of balancing cohesin activity for progenitor differentiation.

Profiling of 3D Genome Organization in Nasopharyngeal Cancer Needle Biopsy Patient Samples by a Modified Hi-C Approach.

Nasopharyngeal cancer (NPC), a cancer derived from epithelial cells in the nasopharynx, is a cancer common in China, Southeast Asia, and Africa. The three-dimensional (3D) genome organization of nasopharyngeal cancer is poorly understood. A major challenge in understanding the 3D genome organization of cancer samples is the lack of a method for the characterization of chromatin interactions in solid cancer needle biopsy samples. Here, we developed Biop-C, a modified in situ Hi-C method using solid cancer needle biopsy samples. We applied Biop-C to characterize three nasopharyngeal cancer solid cancer needle biopsy patient samples. We identified topologically associated domains (TADs), chromatin interaction loops, and frequently interacting regions (FIREs) at key oncogenes in nasopharyngeal cancer from the Biop-C heatmaps. We observed that the genomic features are shared at some important oncogenes, but the patients also display extensive heterogeneity at certain genomic loci. On analyzing the super enhancer landscape in nasopharyngeal cancer cell lines, we found that the super enhancers are associated with FIREs and can be linked to distal genes via chromatin loops in NPC. Taken together, our results demonstrate the utility of our Biop-C method in investigating 3D genome organization in solid cancers.

DNA methylation is required to maintain both DNA replication timing precision and 3D genome organization integrity.

DNA replication timing and three-dimensional (3D) genome organization are associated with distinct epigenome patterns across large domains. However, whether alterations in the epigenome, in particular cancer-related DNA hypomethylation, affects higher-order levels of genome architecture is still unclear. Here, using Repli-Seq, single-cell Repli-Seq, and Hi-C, we show that genome-wide methylation loss is associated with both concordant loss of replication timing precision and deregulation of 3D genome organization. Notably, we find distinct disruption in 3D genome compartmentalization, striking gains in cell-to-cell replication timing heterogeneity and loss of allelic replication timing in cancer hypomethylation models, potentially through the gene deregulation of DNA replication and genome organization pathways. Finally, we identify ectopic H3K4me3-H3K9me3 domains from across large hypomethylated domains, where late replication is maintained, which we purport serves to protect against catastrophic genome reorganization and aberrant gene transcription. Our results highlight a potential role for the methylome in the maintenance of 3D genome regulation.

Proteome Analysis of Condensed Barley Mitotic Chromosomes.

Proteins play a major role in the three-dimensional organization of nuclear genome and its function. While histones arrange DNA into a nucleosome fiber, other proteins contribute to higher-order chromatin structures in interphase nuclei, and mitotic/meiotic chromosomes. Despite the key role of proteins in maintaining genome integrity and transferring hereditary information to daughter cells and progenies, the knowledge about their function remains fragmentary. This is particularly true for the proteins of condensed chromosomes and, in particular, chromosomes of plants. Here, we purified barley mitotic metaphase chromosomes by a flow cytometric sorting and characterized their proteins. Peptides from tryptic protein digests were fractionated either on a cation exchanger or reversed-phase microgradient system before liquid chromatography coupled to tandem mass spectrometry. Chromosomal proteins comprising almost 900 identifications were classified based on a combination of software prediction, available database localization information, sequence homology, and domain representation. A biological context evaluation indicated the presence of several groups of abundant proteins including histones, topoisomerase 2, POLYMERASE 2, condensin subunits, and many proteins with chromatin-related functions. Proteins involved in processes related to DNA replication, transcription, and repair as well as nucleolar proteins were found. We have experimentally validated the presence of FIBRILLARIN 1, one of the nucleolar proteins, on metaphase chromosomes, suggesting that plant chromosomes are coated with proteins during mitosis, similar to those of human and animals. These results improve significantly the knowledge of plant chromosomal proteins and provide a basis for their functional characterization and comparative phylogenetic analyses.

Single-cell measurement of higher-order 3D genome organization with scSPRITE.

Although three-dimensional (3D) genome organization is central to many aspects of nuclear function, it has been difficult to measure at the single-cell level. To address this, we developed 'single-cell split-pool recognition of interactions by tag extension' (scSPRITE). scSPRITE uses split-and-pool barcoding to tag DNA fragments in the same nucleus and their 3D spatial arrangement. Because scSPRITE measures multiway DNA contacts, it generates higher-resolution maps within an individual cell than can be achieved by proximity ligation. We applied scSPRITE to thousands of mouse embryonic stem cells and detected known genome structures, including chromosome territories, active and inactive compartments, and topologically associating domains (TADs) as well as long-range inter-chromosomal structures organized around various nuclear bodies. We observe that these structures exhibit different levels of heterogeneity across the population, with TADs representing dynamic units of genome organization across cells. We expect that scSPRITE will be a critical tool for studying genome structure within heterogeneous populations.

CTCF chromatin residence time controls three-dimensional genome organization, gene expression and DNA methylation in pluripotent cells.

The 11 zinc finger (ZF) protein CTCF regulates topologically associating domain formation and transcription through selective binding to thousands of genomic sites. Here, we replaced endogenous CTCF in mouse embryonic stem cells with green-fluorescent-protein-tagged wild-type or mutant proteins lacking individual ZFs to identify additional determinants of CTCF positioning and function. While ZF1 and ZF8-ZF11 are not essential for cell survival, ZF8 deletion strikingly increases the DNA binding off-rate of mutant CTCF, resulting in reduced CTCF chromatin residence time. Loss of ZF8 results in widespread weakening of topologically associating domains, aberrant gene expression and increased genome-wide DNA methylation. Thus, important chromatin-templated processes rely on accurate CTCF chromatin residence time, which we propose depends on local sequence and chromatin context as well as global CTCF protein concentration.

Hi\u2013C interaction graph analysis reveals the impact of histone modifications in chromatin shape

Chromosome conformation capture experiments such as Hi–C map the three-dimensional spatial organization of genomes in a genome-wide scale. Even though Hi–C interactions are not biased towards any of the histone modifications, previous analysis has revealed denser interactions around many histone modifications. Nevertheless, simultaneous effects of these modifications in Hi–C interaction graph have not been fully characterized yet, limiting our understanding of genome shape. Here, we propose ChromatinCoverage and its extension TemporalPrizeCoverage methods to decompose Hi–C interaction graph in terms of known histone modifications. Both methods are based on set multicover with pairs, where each Hi–C interaction is tried to be covered by histone modification pairs. We find 4 histone modifications H3K4me1, H3K4me3, H3K9me3, H3K27ac to be significantly predictive of most Hi–C interactions across species, cell types and cell cycles. The proposed methods are quite effective in predicting Hi–C interactions and topologically-associated domains in one species, given it is trained on another species or cell types. Overall, our findings reveal the impact of subset of histone modifications in chromatin shape via Hi–C interaction graph.