Classifying leukemia types with chromatin conformation data

A conserved HOTAIRM1-HOXA1 regulatory axis contributes early to neuronal differentiation

HOTAIRM1 is unlike most long non-coding RNAs in that its sequence is highly conserved across mammals. Such evolutionary conservation points to it having a role in key cellular processes. We previously reported that HOTAIRM1 is required to curb premature activation of downstream HOXA genes in a cell model recapitulating their sequential induction during development. We found that it regulates 3' HOXA gene expression by a mechanism involving epigenetic and three-dimensional chromatin changes. Here we show that HOTAIRM1 participates in proper progression through the early stages of neuronal differentiation. We found that it can associate with the HOXA1 transcription factor and contributes to its downstream transcriptional program. Particularly, HOTAIRM1 affects the NANOG/POU5F1/SOX2 core pluripotency network maintaining an undifferentiated cell state. HOXA1 depletion similarly perturbed expression of these pluripotent factors, suggesting that HOTAIRM1 is a modulator of this transcription factor pathway. Also, given that binding of HOTAIRM1 to HOXA1 was observed in different cell types and species, our results point to this ribonucleoprotein complex as an integral part of a conserved HOTAIRM1-HOXA1 regulatory axis modulating the transition from a pluripotent to a differentiated neuronal state.

The retinoic acid receptor co-factor NRIP1 is uniquely upregulated and represents a therapeutic target in acute myeloid leukemia with chromosome 3q rearrangements

Aberrant expression of Ecotropic Viral Integration Site 1 (EVI1) is a hallmark of acute myeloid leukemia (AML) with inv(3) or t(3;3), which is a disease subtype with especially poor outcome. In studying transcriptomes from AML patients with chromosome 3q rearrangements, we identified a significant upregulation of the Nuclear Receptor Interacting Protein 1 (NRIP1) as well as its adjacent non-coding RNA LOC101927745. Utilizing transcriptomic and epigenomic data from over 900 primary samples from patients as well as genetic and transcriptional engineering approaches, we have identified several mechanisms that can lead to upregulation of NRIP1 in AML. We hypothesize that the LOC101927745 transcription start site harbors a context-dependent enhancer that is bound by EVI1, causing upregulation of NRIP1 in AML with chromosome 3 abnormalities. Furthermore, we showed that NRIP1 knockdown negatively affects the proliferation and survival of 3qrearranged AML cells and increases their sensitivity to all-trans retinoic acid, suggesting that NRIP1 is relevant for the pathogenesis of inv(3)/t(3;3) AML and could serve as a novel therapeutic target in myeloid malignancies with 3q abnormalities.

Chromatin Conformation in Development and Disease

Chromatin domains and loops are important elements of chromatin structure and dynamics, but much remains to be learned about their exact biological role and nature. Topological associated domains and functional loops are key to gene expression and hold the answer to many questions regarding developmental decisions and diseases. Here, we discuss new findings, which have linked chromatin conformation with development, differentiation and diseases and hypothesized on various models while integrating all recent findings on how chromatin architecture affects gene expression during development, evolution and disease.

Contribution of CTCF binding to transcriptional activity at the HOXA locus in NPM1-mutant AML cells

Transcriptional regulation of the HOXA genes is thought to involve CTCF-mediated chromatin loops and the opposing actions of the COMPASS and Polycomb epigenetic complexes. We investigated the role of these mechanisms at the HOXA cluster in AML cells with the common NPM1c mutation, which express both HOXA and HOXB genes. CTCF binding at the HOXA locus is conserved across primary AML samples, regardless of HOXA gene expression, and defines a continuous chromatin domain marked by COMPASS-associated histone H3 trimethylation in NPM1-mutant primary AML samples. Profiling of the three-dimensional chromatin architecture in primary AML samples with the NPM1c mutation identified chromatin loops between the HOXA cluster and loci in the SNX10 and SKAP2 genes, and an intergenic region located 1.4 Mbp upstream of the HOXA locus. Deletion of CTCF binding sites in the NPM1-mutant OCI-AML3 AML cell line reduced multiple long-range interactions, but resulted in CTCF-independent loops with sequences in SKAP2 that were marked by enhancer-associated histone modifications in primary AML samples. HOXA gene expression was maintained in CTCF binding site mutants, indicating that transcriptional activity at the HOXA locus in NPM1-mutant AML cells may be sustained through persistent interactions with SKAP2 enhancers, or by intrinsic factors within the HOXA gene cluster.

Characterization of GJB2 cis-regulatory elements in the DFNB1 locus

Although most disease-causing variants are within coding region of genes, it is now well established that cis-acting regulatory sequences, depending on 3D-chromatin organization, are required for temporal and spatial control of gene expression. Disruptions of such regulatory elements and/or chromatin conformation are likely to play a critical role in human genetic disease. Hence, recurrent monoallelic cases, who present the most common hereditary type of nonsyndromic hearing loss (i.e., DFNB1), carry only one identified pathogenic allele. This strongly suggests the presence of uncharacterized distal cis-acting elements in the missing allele. Here within, we study the spatial organization of a large DFNB1 locus encompassing the gap junction protein beta 2 (GJB2) gene, the most frequently mutated gene in this inherited hearing loss phenotype, with the chromosome conformation capture carbon copy technology (5C). By combining this approach with functional activity reporter assays and mapping of CCCTC-binding factor (CTCF) along the DFNB1 locus, we identify a novel set of cooperating GJB2 cis-acting elements and suggest a DFNB1 three-dimensional looping regulation model.

Rich Chromatin Structure Prediction from Hi-C Data

Recent studies involving the 3-dimensional conformation of chromatin have revealed the important role it has to play in different processes within the cell. These studies have also led to the discovery of densely interacting segments of the chromosome, called topologically associating domains. The accurate identification of these domains from Hi-C interaction data is an interesting and important computational problem for which numerous methods have been proposed. Unfortunately, most existing algorithms designed to identify these domains assume that they are non-overlapping whereas there is substantial evidence to believe a nested structure exists. We present a methodology to predict hierarchical chromatin domains using chromatin conformation capture data. Our method predicts domains at different resolutions, calculated using intrinsic properties of the chromatin data, and effectively clusters these to construct the hierarchy. At each individual level, the domains are non-overlapping in such a way that the intra-domain interaction frequencies are maximized. We show that our predicted structure is highly enriched for actively transcribing housekeeping genes and various chromatin markers, including CTCF, around the domain boundaries. We also show that large-scale domains, at multiple resolutions within our hierarchy, are conserved across cell types and species. We also provide comparisons against existing tools for extracting hierarchical domains. Our software, Matryoshka, is written in C++11 and licensed under GPL v3; it is available at https://github.com/COMBINE-lab/matryoshka.

Modeling three-dimensional genomic organization in evolution and pathogenesis

The regulation of gene expression is mediated via the complex three-dimensional (3D) conformation of the genetic material and its interactions with various intracellular factors. Various experimental and computational approaches have been developed in recent years for understating the relation between the 3D conformation of the genome and the phenotypes of cells in normal condition and diseases. In this review, we will discuss novel approaches for analyzing and modeling the 3D genomic conformation, focusing on deciphering disease-causing mutations that affect gene expression. We conclude that as this is a very challenging mission, an important direction should involve the comparative analysis of various 3D models from various organisms or cells.

Novel long-range regulatory mechanisms controlling PKD2 gene expression

Background

Cis-regulatory elements control gene expression over large distances through the formation of chromatin loops, which allow contact between enhancers and gene promoters. Alterations in cis-acting regulatory systems could be linked to human genetic diseases. Here, we analyse the spatial organization of a large region spanning the polycystic kidney disease 2 (PKD2) gene, one of the genes responsible of autosomal dominant polycystic kidney disease (ADPKD).

Results

By using chromosome conformation capture carbon copy (5C) technology in primary human renal cyst epithelial cells, we identify novel contacts of the PKD2 promoter with chromatin regions, which display characteristics of regulatory elements. In parallel, by using functional analysis with a reporter assay, we demonstrate that three DNAse I hypersensitive sites regions are involved in the regulation of PKD2 gene expression.

Conclusions

Finally, through alignment of CCCTC-binding factor (CTCF) sites, we suggest that these novel enhancer elements are brought to the PKD2 promoter by chromatin looping via the recruitment of CTCF.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4892-6) contains supplementary material, which is available to authorized users.

Storage, visualization, and navigation of 3D genomics data

The field of 3D genomics grew at increasing rates in the last decade. The volume and complexity of 2D and 3D data produced is progressively outpacing the capacities of the technology previously used for distributing genome sequences. The emergence of new technologies provides also novel opportunities for the development of innovative approaches. In this paper, we review the state-of-the-art computing technology, as well as the solutions adopted by the platforms currently available.

Alterations of specific chromatin conformation affect ATRA-induced leukemia cell differentiation

Chromatin conformation plays a key role in regulating gene expression and controlling cell differentiation. However, the whole-genome chromatin conformation changes that occur during leukemia cell differentiation are poorly understood. Here, we characterized the changes in chromatin conformation, histone states, chromatin accessibility, and gene expression using an all-trans retinoic acid (ATRA)-induced HL-60 cell differentiation model. The results showed that the boundaries of topological associated domains (TADs) were stable during differentiation; however, the chromatin conformations within several specific TADs were obviously changed. By combining H3K4me3, H3K27ac, and Hi-C signals, we annotated the differential gene-regulatory chromatin interactions upon ATRA induction. The gains and losses of the gene-regulatory chromatin interactions are significantly correlated with gene expression and chromatin accessibility. Finally, we found that the loss of GATA2 expression and DNA binding are crucial for the differentiation process, and changes in the chromatin structure around the GATA2 regulate its expression upon ATRA induction. This study provided both statistical insights and experimental details regarding the relationship between chromatin conformation changes and transcription regulation during leukemia cell differentiation, and the results suggested that the chromatin conformation is a new type of potential drug target for cancer therapy.

Novel insights into chromosomal conformations in cancer

Exploring gene function is critical for understanding the complexity of life. DNA sequences and the three-dimensional organization of chromatin (chromosomal interactions) are considered enigmatic factors underlying gene function, and interactions between two distant fragments can regulate transactivation activity via mediator proteins. Thus, a series of chromosome conformation capture techniques have been developed, including chromosome conformation capture (3C), circular chromosome conformation capture (4C), chromosome conformation capture carbon copy (5C), and high-resolution chromosome conformation capture (Hi-C). The application of these techniques has expanded to various fields, but cancer remains one of the major topics. Interactions mediated by proteins or long noncoding RNAs (lncRNAs) are typically found using 4C-sequencing and chromatin interaction analysis by paired-end tag sequencing (ChIA-PET). Currently, Hi-C is used to identify chromatin loops between cancer risk-associated single-nucleotide polymorphisms (SNPs) found by genome-wide association studies (GWAS) and their target genes. Chromosomal conformations are responsible for altered gene regulation through several typical mechanisms and contribute to the biological behavior and malignancy of different tumors, particularly prostate cancer, breast cancer and hematologic neoplasms. Moreover, different subtypes may exhibit different 3D-chromosomal conformations. Thus, C-tech can be used to help diagnose cancer subtypes and alleviate cancer progression by destroying specific chromosomal conformations. Here, we review the fundamentals and improvements in chromosome conformation capture techniques and their clinical applications in cancer to provide insight for future research.

Dynamic chromatin accessibility modeled by Markov process of randomly-moving molecules in the 3D genome

Chromatin three-dimensional (3D) structure plays critical roles in gene expression regulation by influencing locus interactions and accessibility of chromatin regions. Here we propose a Markov process model to derive a chromosomal equilibrium distribution of randomly-moving molecules as a functional consequence of spatially organized genome 3D structures. The model calculates steady-state distributions (SSD) from Hi-C data as quantitative measures of each chromatin region's dynamic accessibility for transcription factors and histone modification enzymes. Different from other Hi-C derived features such as compartment A/B and interaction hubs, or traditional methods measuring chromatin accessibility such as DNase-seq and FAIRE-seq, SSD considers both chromatin–chromatin and protein–chromatin interactions. Through our model, we find that SSD could capture the chromosomal equilibrium distributions of activation histone modifications and transcription factors. Compared with compartment A/B, SSD has higher correlations with the binding of these histone modifications and transcription factors. In addition, we find that genes located in high SSD regions tend to be expressed at higher level. Furthermore, we track the change of genome organization during stem cell differentiation, and propose a two-stage model to explain the dynamic change of SSD and gene expression during differentiation, where chromatin organization genes first gain chromatin accessibility and are expressed before lineage-specific genes do. We conclude that SSD is a novel and better measure of dynamic chromatin activity and accessibility.

MLL-AF9 and MLL-AF4 oncofusion proteins bind a distinct enhancer repertoire and target the RUNX1 program in 11q23 acute myeloid leukemia

In 11q23 leukemias, the N-terminal part of the mixed lineage leukemia (MLL) gene is fused to >60 different partner genes. In order to define a core set of MLL rearranged targets, we investigated the genome-wide binding of the MLL-AF9 and MLL-AF4 fusion proteins and associated epigenetic signatures in acute myeloid leukemia (AML) cell lines THP-1 and MV4-11. We uncovered both common as well as specific MLL-AF9 and MLL-AF4 target genes, which were all marked by H3K79me2, H3K27ac and H3K4me3. Apart from promoter binding, we also identified MLL-AF9 and MLL-AF4 binding at specific subsets of non-overlapping active distal regulatory elements. Despite this differential enhancer binding, MLL-AF9 and MLL-AF4 still direct a common gene program, which represents part of the RUNX1 gene program and constitutes of CD34⁺ and monocyte-specific genes. Comparing these data sets identified several zinc finger transcription factors (TFs) as potential MLL-AF9 co-regulators. Together, these results suggest that MLL fusions collaborate with specific subsets of TFs to deregulate the RUNX1 gene program in 11q23 AMLs.

Emerging Roles for Epigenetic Programming in the Control of Inflammatory Signaling Integration in Heath and Disease

Macrophages and dendritic cells initiate the innate immune response to infection and injury and contribute to inflammatory signaling to maintain the homeostasis of various tissues, which includes resident macrophages for the elimination of invading microorganisms and tissue damage. Inappropriate inflammatory signaling can lead to persistent inflammation and further develop into autoimmune and inflammation-associated diseases. Inflammatory signaling pathways have been well characterized, but how these signaling pathways are converted into sustained and diverse patterns of expression of cytokines, chemokines, and other genes in response to environmental challenges is unclear. Emerging evidence suggests the important role of epigenetic mechanisms in finely tuning the outcome of the host innate immune response. An understanding of epigenetic regulation of innate immune cell identity and function will enable the identification of the mechanism between gene-specific host defenses and inflammatory disease and will also allow for exploration of the program of innate immune memory in health and disease. This information could be used to develop therapeutic agents to enhance the host response, preventing chronic inflammation through preserving tissues and signaling integrity.

Disruption of the 3D cancer genome blueprint

Recent advances in chromosome conformation capture technologies are improving the current appreciation of how 3D genome architecture affects its function in different cell types and disease. Long-range chromatin interactions are organized into topologically associated domains, which are known to play a role in constraining gene expression patterns. However, in cancer cells there are alterations in the 3D genome structure, which impacts on gene regulation. Disruption of topologically associated domains architecture can result in alterations in chromatin interactions that bring new regulatory elements and genes together, leading to altered expression of oncogenes and tumor suppressor genes. Here, we discuss the impact of genetic and epigenetic changes in cancer and how this affects the spatial organization of chromatin. Understanding how disruptions to the 3D architecture contribute to the cancer genome will provide novel insights into the principles of epigenetic gene regulation in cancer and mechanisms responsible for cancer associated mutations and rearrangements.

Feasibility of the AML profiler (Skyline™ Array) for patient risk stratification in a multicentre trial: a preliminary comparison with the conventional approach

Deoxyribonucleic acid microarrays allow researchers to measure mRNA levels of thousands of genes in a single experiment and could be useful for diagnostic purposes in patients with acute myeloid leukaemia (AML). We assessed the feasibility of the AML profiler (Skyline™ Array) in genetic stratification of patients with de novo AML and compared the results with those obtained using the standard cytogenetic and molecular approach. Diagnostic bone marrow from 31 consecutive de novo AML cases was used to test MLL-PTD, FLT3-ITD and TKD, NPM1 and CEBPAdm mutations. Purified RNA was used to assess RUNX1-RUNX1T1, PML-RARα and CBFβ-MYH11 rearrangements. RNA remnants underwent gene expression profiling analysis using the AML profiler, which detects chromosomal aberrations: t(8;21), t(15;17), inv(16), mutations (CEBPAdm, ABD-NPM1) and BAALC and EVI1 expression. Thirty cases were successfully analysed with both methods. Five cases had FLT3-ITD. In one case, a t(8;21) was correctly detected by both methods. Four cases had inv(16); in one, the RNA quality was unsatisfactory and it was not hybridized, and in the other three, the AML profiler detected the genetic lesion - this being a rare type I translocation in one case. Two cases with acute promyelocytic leukaemia were diagnosed by both methods. Results for NPM1 mutations were concordant in all but two cases (2/11, non-ABD mutations). Analysis of costs and turnaround times showed that the AML profiler was no more expensive than the conventional molecular approach. These results suggest that the AML profiler could be useful in multicentre trials to rapidly identify patients with AML with a good prognosis. Copyright © 2016 John Wiley & Sons, Ltd.

Spatial Organization of Epigenomes

The role of genome architecture in transcription regulation has become the focus of an increasing number of studies over the past decade. Chromatin organization can have a significant impact on gene expression by promoting or restricting the physical proximity between regulatory DNA elements. Given that any change in chromatin state has the potential to alter DNA folding and the proximity between control elements, the spatial organization of chromatin is inherently linked to its molecular composition. In this review, we explore how modulators of chromatin state and organization might keep gene expression in check. We discuss recent findings and present some of the less well-studied aspects of spatial genome organization such as chromatin dynamics and regulation by non-coding RNAs.

The three-dimensional cancer genome

The past decade of cancer research has ushered in a comprehensive understanding of the way that the sequence of the genome can be co-opted during the process of tumorigenesis. However, only recently has the epigenome, and in particular the three-dimensional topology of chromatin, been implicated in cancer progression. Here we review recent findings of how the cancer genome is regulated and dysregulated to effect changes in 3D genome topology. We discuss the impact of the spatial organization of the genome on the frequency of tumorigenic chromosomal translocations and the effects of disruption of the proteins responsible for the establishment of chromatin loops. Alteration of the three-dimensional cancer genome is a rapidly emerging hallmark of multiple cancer subtypes.

C-ing the Genome: A Compendium of Chromosome Conformation Capture Methods to Study Higher-Order Chromatin Organization

Three-dimensional organization of the chromatin has important roles in transcription, replication, DNA repair, and pathologic events such as translocations. There are two fundamental ways to study higher-order chromatin organization: microscopic and molecular approaches. In this review, we briefly introduce the molecular approaches, focusing on chromosome conformation capture or "3C" technology and its derivatives, which can be used to probe chromatin folding at resolutions beyond that provided by microscopy techniques. We further discuss the different types of data generated by the 3C-based methods and how they can be used to answer distinct biological questions.

Analysis of long-range interactions in primary human cells identifies cooperative CFTR regulatory elements

A mechanism by which control DNA elements regulate transcription over large linear genomic distances is by achieving close physical proximity with genes, and looping of the intervening chromatin paths. Alterations of such regulatory 'chromatin looping' systems are likely to play a critical role in human genetic disease at large. Here, we studied the spatial organization of a ≈790 kb locus encompassing the cystic fibrosis transmembrane conductance regulator (CFTR) gene. Dysregulation of CFTR is responsible for cystic fibrosis, which is the most common lethal genetic disorder in Caucasian populations. CFTR is a relatively large gene of 189 kb with a rather complex tissue-specific and temporal expression profile. We used chromatin conformation at the CFTR locus to identify new DNA sequences that regulate its transcription. By comparing 5C chromatin interaction maps of the CFTR locus in expressing and non-expressing human primary cells, we identified several new contact points between the CFTR promoter and its surroundings, in addition to regions featuring previously described regulatory elements. We demonstrate that two of these novel interacting regions cooperatively increase CFTR expression, and suggest that the new enhancer elements located on either side of the gene are brought together through chromatin looping via CTCF.

Chromatin interaction analysis reveals changes in small chromosome and telomere clustering between epithelial and breast cancer cells

Background

Higher-order chromatin structure is often perturbed in cancer and other pathological states. Although several genetic and epigenetic differences have been charted between normal and breast cancer tissues, changes in higher-order chromatin organization during tumorigenesis have not been fully explored. To probe the differences in higher-order chromatin structure between mammary epithelial and breast cancer cells, we performed Hi-C analysis on MCF-10A mammary epithelial and MCF-7 breast cancer cell lines.

Results

Our studies reveal that the small, gene-rich chromosomes chr16 through chr22 in the MCF-7 breast cancer genome display decreased interaction frequency with each other compared to the inter-chromosomal interaction frequency in the MCF-10A epithelial cells. Interestingly, this finding is associated with a higher occurrence of open compartments on chr16–22 in MCF-7 cells. Pathway analysis of the MCF-7 up-regulated genes located in altered compartment regions on chr16–22 reveals pathways related to repression of WNT signaling. There are also differences in intra-chromosomal interactions between the cell lines; telomeric and sub-telomeric regions in the MCF-10A cells display more frequent interactions than are observed in the MCF-7 cells.

Conclusions

We show evidence of an intricate relationship between chromosomal organization and gene expression between epithelial and breast cancer cells. Importantly, this work provides a genome-wide view of higher-order chromatin dynamics and a resource for studying higher-order chromatin interactions in two cell lines commonly used to study the progression of breast cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0768-0) contains supplementary material, which is available to authorized users.

An Overview of Genome Organization and How We Got There: from FISH to Hi-C

In humans, nearly two meters of genomic material must be folded to fit inside each micrometer-scale cell nucleus while remaining accessible for gene transcription, DNA replication, and DNA repair. This fact highlights the need for mechanisms governing genome organization during any activity and to maintain the physical organization of chromosomes at all times. Insight into the functions and three-dimensional structures of genomes comes mostly from the application of visual techniques such as fluorescence in situ hybridization (FISH) and molecular approaches including chromosome conformation capture (3C) technologies. Recent developments in both types of approaches now offer the possibility of exploring the folded state of an entire genome and maybe even the identification of how complex molecular machines govern its shape. In this review, we present key methodologies used to study genome organization and discuss what they reveal about chromosome conformation as it relates to transcription regulation across genomic scales in mammals.

Structural heterogeneity and functional diversity of topologically associating domains in mammalian genomes

Recent chromosome conformation capture (3C) derived techniques have revealed that topologically associating domain (TAD) is a pervasive element in chromatin three-dimensional (3D) organization. However, there is currently no parameter to quantitatively measure the structural characteristics of TADs, thus obscuring our understanding on the structural and functional differences among TADs. Based on our finding that there exist intrinsic chromatin interaction patterns in TADs, we define a theoretical parameter, called aggregation preference (AP), to characterize TAD structures by capturing the interaction aggregation degree. Applying this defined parameter to 11 Hi-C data sets generated by both traditional and in situ Hi-C experimental pipelines, our analyses reveal that heterogeneous structures exist among TADs, and this structural heterogeneity is significantly correlated to DNA sequences, epigenomic signals and gene expressions. Although TADs can be stable in genomic positions across cell lines, structural comparisons show that a considerable number of stable TADs undergo significantly structural rearrangements during cell changes. Moreover, the structural change of TAD is tightly associated with its transcription remodeling. Altogether, the theoretical parameter defined in this work provides a quantitative method to link structural characteristics and biological functions of TADs, and this linkage implies that chromatin interaction pattern has the potential to mark transcription activity in TADs.