HiCUP: pipeline for mapping and processing Hi-C data.

HiCUP is a pipeline for processing sequence data generated by Hi-C and Capture Hi-C (CHi-C) experiments, which are techniques used to investigate three-dimensional genomic organisation. The pipeline maps data to a specified reference genome and removes artefacts that would otherwise hinder subsequent analysis. HiCUP also produces an easy-to-interpret yet detailed quality control (QC) report that assists in refining experimental protocols for future studies. The software is freely available and has already been used for processing Hi-C and CHi-C data in several recently published peer-reviewed studies

Global reorganization of the nuclear landscape in senescent cells.

Cellular senescence has been implicated in tumor suppression, development, and aging and is accompanied by large-scale chromatin rearrangements, forming senescence-associated heterochromatic foci (SAHF). However, how the chromatin is reorganized during SAHF formation is poorly understood. Furthermore, heterochromatin formation in senescence appears to contrast with loss of heterochromatin in Hutchinson-Gilford progeria. We mapped architectural changes in genome organization in cellular senescence using Hi-C. Unexpectedly, we find a dramatic sequence- and lamin-dependent loss of local interactions in heterochromatin. This change in local connectivity resolves the paradox of opposing chromatin changes in senescence and progeria. In addition, we observe a senescence-specific spatial clustering of heterochromatic regions, suggesting a unique second step required for SAHF formation. Comparison of embryonic stem cells (ESCs), somatic cells, and senescent cells shows a unidirectional loss in local chromatin connectivity, suggesting that senescence is an endpoint of the continuous nuclear remodelling process during differentiation

Long-Range Enhancer Interactions Are Prevalent in Mouse Embryonic Stem Cells and Are Reorganized upon Pluripotent State Transition.

Transcriptional enhancers, including super-enhancers (SEs), form physical interactions with promoters to regulate cell-type-specific gene expression. SEs are characterized by high transcription factor occupancy and large domains of active chromatin, and they are commonly assigned to target promoters using computational predictions. How promoter-SE interactions change upon cell state transitions, and whether transcription factors maintain SE interactions, have not been reported. Here, we used promoter-capture Hi-C to identify promoters that interact with SEs in mouse embryonic stem cells (ESCs). We found that SEs form complex, spatial networks in which individual SEs contact multiple promoters, and a rewiring of promoter-SE interactions occurs between pluripotent states. We also show that long-range promoter-SE interactions are more prevalent in ESCs than in epiblast stem cells (EpiSCs) or Nanog-deficient ESCs. We conclude that SEs form cell-type-specific interaction networks that are partly dependent on core transcription factors, thereby providing insights into the gene regulatory organization of pluripotent cells.P.J.R.-G. is supported by the Wellcome Trust (WT093736), Biotechnology and Biological Sciences Research Council (BB/M022285/1 and BB/P013406/1), and the European Commission Network of Excellence EpiGeneSys (HEALTH-F4-2010-257082). This work was also supported by the following grants to P.F.: Medical Research Council (MR/L007150/1, MC_UP_1302/1, MC_UP_1302/3, MC_UP_1302/5), and Biotechnology and Biological Sciences Research Council (BB/J004480/1)

The global and promoter-centric 3D genome organization temporally resolved during a circadian cycle

Funder: FP7 Ideas: European Research Council; doi: http://dx.doi.org/10.13039/100011199; Grant(s): 259743Abstract: Background: Circadian gene expression is essential for organisms to adjust their physiology and anticipate daily changes in the environment. The molecular mechanisms controlling circadian gene transcription are still under investigation. In particular, how chromatin conformation at different genomic scales and regulatory elements impact rhythmic gene expression has been poorly characterized. Results: Here we measure changes in the spatial chromatin conformation in mouse liver using genome-wide and promoter-capture Hi-C alongside daily oscillations in gene transcription. We find topologically associating domains harboring circadian genes that switch assignments between the transcriptionally active and inactive compartment at different hours of the day, while their boundaries stably maintain their structure over time. To study chromatin contacts of promoters at high resolution over time, we apply promoter capture Hi-C. We find circadian gene promoters displayed a maximal number of chromatin contacts at the time of their peak transcriptional output. Furthermore, circadian genes, as well as contacted and transcribed regulatory elements, reach maximal expression at the same timepoints. Anchor sites of circadian gene promoter loops are enriched in DNA binding sites for liver nuclear receptors and other transcription factors, some exclusively present in either rhythmic or stable contacts. Finally, by comparing the interaction profiles between core clock and output circadian genes, we show that core clock interactomes are more dynamic compared to output circadian genes. Conclusion: Our results identify chromatin conformation dynamics at different scales that parallel oscillatory gene expression and characterize the repertoire of regulatory elements that control circadian gene transcription through rhythmic or stable chromatin configurations

A role for insulator elements in the regulation of gene expression response to hypoxia

Hypoxia inducible factor (HIF) up-regulates the transcription of a few hundred genes required for the adaptation to hypoxia. This restricted set of targets is in sharp contrast with the widespread distribution of the HIF binding motif throughout the genome. Here, we investigated the transcriptional response of GYS1 and RUVBL2 genes to hypoxia to understand the mechanisms that restrict HIF activity toward specific genes. GYS1 and RUVBL2 genes are encoded by opposite DNA strands and separated by a short intergenic region (~1 kb) that contains a functional hypoxia response element equidistant to both genes. However, hypoxia induced the expression of GYS1 gene only. Analysis of the transcriptional response of chimeric constructs derived from the intergenic region revealed an inhibitory sequence whose deletion allowed RUVBL2 induction by HIF. Enhancer blocking assays, performed in cell culture and transgenic zebrafish, confirmed the existence of an insulator element within this inhibitory region that could explain the differential regulation of GYS1 and RUVBL2 by hypoxia. Hence, in this model, the selective response to HIF is achieved with the aid of insulator elements. This is the first report suggesting a role for insulators in the regulation of differential gene expression in response to environmental signals

An insulator embedded in the chicken α-globin locus regulates chromatin domain configuration and differential gene expression

Genome organization into transcriptionally active domains denotes one of the first levels of gene expression regulation. Although the chromatin domain concept is generally accepted, only little is known on how domain organization impacts the regulation of differential gene expression. Insulators might hold answers to address this issue as they delimit and organize chromatin domains. We have previously identified a CTCF-dependent insulator with enhancer-blocking activity embedded in the 5′ non-coding region of the chicken α-globin domain. Here, we demonstrate that this element, called the αEHS-1.4 insulator, protects a transgene against chromosomal position effects in stably transfected cell lines and transgenic mice. We found that this insulator can create a regulated chromatin environment that coincides with the onset of adult α-globin gene expression. Furthermore, such activity is in part dependent on the in vivo regulated occupancy of CTCF at the αEHS-1.4 element. Insulator function is also regulated by CTCF poly(ADP-ribosyl)ation. Our results suggest that the αEHS-1.4 insulator contributes in organizing the chromatin structure of the α-globin gene domain and prevents activation of adult α-globin gene expression at the erythroblast stage via CTCF

Gain of DNA methylation is enhanced in the absence of CTCF at the human retinoblastoma gene promoter

<p>Abstract</p> <p>Background</p> <p>Long-term gene silencing throughout cell division is generally achieved by DNA methylation and other epigenetic processes. Aberrant DNA methylation is now widely recognized to be associated with cancer and other human diseases. Here we addressed the contribution of the multifunctional nuclear factor CTCF to the epigenetic regulation of the human <it>retinoblastoma </it>(<it>Rb</it>) gene promoter in different tumoral cell lines.</p> <p>Methods</p> <p>To assess the DNA methylation status of the <it>Rb </it>promoter, genomic DNA from stably transfected human erythroleukemic K562 cells expressing a <it>GFP </it>reporter transgene was transformed with sodium bisulfite, and then PCR-amplified with modified primers and sequenced. Single- and multi-copy integrants with the CTCF binding site mutated were isolated and characterized by Southern blotting. Silenced transgenes were reactivated using 5-aza-2'-deoxycytidine and Trichostatin-A, and their expression was monitored by fluorescent cytometry. <it>Rb </it>gene expression and protein abundance were assessed by RT-PCR and Western blotting in three different glioma cell lines, and DNA methylation of the promoter region was determined by sodium bisulfite sequencing, together with CTCF dissociation and methyl-CpG-binding protein incorporation by chromatin immunoprecipitation assays.</p> <p>Results</p> <p>We found that the inability of CTCF to bind to the <it>Rb </it>promoter causes a dramatic loss of gene expression and a progressive gain of DNA methylation.</p> <p>Conclusions</p> <p>This study indicates that CTCF plays an important role in maintaining the <it>Rb </it>promoter in an optimal chromatin configuration. The absence of CTCF induces a rapid epigenetic silencing through a progressive gain of DNA methylation. Consequently, CTCF can now be seen as one of the epigenetic components that allows the proper configuration of tumor suppressor gene promoters. Its aberrant dissociation can then predispose key genes in cancer cells to acquire DNA methylation and epigenetic silencing.</p

An HDAC9-MALAT1-BRG1 complex mediates smooth muscle dysfunction in thoracic aortic aneurysm

Thoracic aortic aneurysm (TAA) has been associated with mutations affecting members of the TGF-β signaling pathway, or components and regulators of the vascular smooth muscle cell (VSMC) actomyosin cytoskeleton. Although both clinical groups present similar phenotypes, the existence of potential common mechanisms of pathogenesis remain obscure. Here we show that mutations affecting TGF-β signaling and VSMC cytoskeleton both lead to the formation of a ternary complex comprising the histone deacetylase HDAC9, the chromatin-remodeling enzyme BRG1, and the long noncoding RNA MALAT1. The HDAC9–MALAT1–BRG1 complex binds chromatin and represses contractile protein gene expression in association with gain of histone H3-lysine 27 trimethylation modifications. Disruption of Malat1 or Hdac9 restores contractile protein expression, improves aortic mural architecture, and inhibits experimental aneurysm growth. Thus, we highlight a shared epigenetic pathway responsible for VSMC dysfunction in both forms of TAA, with potential therapeutic implication for other known HDAC9-associated vascular diseases

Lineage-specific dynamic and pre-established enhancer–promoter contacts cooperate in terminal differentiation

Chromosome conformation is an important feature of metazoan gene regulation; however, enhancer–promoter contact remodeling during cellular differentiation remains poorly understood. To address this, genome-wide promoter capture Hi-C (CHi-C) was performed during epidermal differentiation. Two classes of enhancer–promoter contacts associated with differentiation-induced genes were identified. The first class ('gained') increased in contact strength during differentiation in concert with enhancer acquisition of the H3K27ac activation mark. The second class ('stable') were pre-established in undifferentiated cells, with enhancers constitutively marked by H3K27ac. The stable class was associated with the canonical conformation regulator cohesin, whereas the gained class was not, implying distinct mechanisms of contact formation and regulation. Analysis of stable enhancers identified a new, essential role for a constitutively expressed, lineage-restricted ETS-family transcription factor, EHF, in epidermal differentiation. Furthermore, neither class of contacts was observed in pluripotent cells, suggesting that lineage-specific chromatin structure is established in tissue progenitor cells and is further remodeled in terminal differentiation

5C analysis of the Epidermal Differentiation Complex locus reveals distinct chromatin interaction networks between gene-rich and gene-poor TADs in skin epithelial cells

YesMammalian genomes contain several dozens of large (>0.5 Mbp) lineage-specific gene loci harbouring functionally related genes. However, spatial chromatin folding, organization of the enhancer-promoter networks and their relevance to Topologically Associating Domains (TADs) in these loci remain poorly understood. TADs are principle units of the genome folding and represents the DNA regions within which DNA interacts more frequently and less frequently across the TAD boundary. Here, we used Chromatin Conformation Capture Carbon Copy (5C) technology to characterize spatial chromatin interaction network in the 3.1 Mb Epidermal Differentiation Complex (EDC) locus harbouring 61 functionally related genes that show lineage-specific activation during terminal keratinocyte differentiation in the epidermis. 5C data validated by 3D-FISH demonstrate that the EDC locus is organized into several TADs showing distinct lineage-specific chromatin interaction networks based on their transcription activity and the gene-rich or gene-poor status. Correlation of the 5C results with genome-wide studies for enhancer-specific histone modifications (H3K4me1 and H3K27ac) revealed that the majority of spatial chromatin interactions that involves the gene-rich TADs at the EDC locus in keratinocytes include both intra- and inter-TAD interaction networks, connecting gene promoters and enhancers. Compared to thymocytes in which the EDC locus is mostly transcriptionally inactive, these interactions were found to be keratinocyte-specific. In keratinocytes, the promoter-enhancer anchoring regions in the gene-rich transcriptionally active TADs are enriched for the binding of chromatin architectural proteins CTCF, Rad21 and chromatin remodeler Brg1. In contrast to gene-rich TADs, gene-poor TADs show preferential spatial contacts with each other, do not contain active enhancers and show decreased binding of CTCF, Rad21 and Brg1 in keratinocytes. Thus, spatial interactions between gene promoters and enhancers at the multi-TAD EDC locus in skin epithelial cells are cell type-specific and involve extensive contacts within TADs as well as between different gene-rich TADs, forming the framework for lineage-specific transcription.This study was supported by the grants 5R01AR064580 and 1RO1AR071727 to VAB, TKS and AAS, as well as by the grants from MRC (MR/ M010015/1) and BBSRC (BB/K010050/1) to VAB