Developmental enhancers initiate transcription and are fundamental to our understanding of developmental networks, evolution and disease. Despite their importance, the properties governing enhancer-promoter interactions and their dynamics during embryogenesis remain unclear. At the β-globin locus, enhancer-promoter interactions appear dynamic and cell-type specific, whereas at the HoxD locus they are stable and ubiquitous, being present in tissues where the target genes are not expressed. The extent to which preformed enhancer-promoter conformations exist at other, more typical, loci and how transcription is eventually triggered is unclear. Here we generated a high-resolution map of enhancer three-dimensional contacts during Drosophila embryogenesis, covering two developmental stages and tissue contexts, at unprecedented resolution. Although local regulatory interactions are common, long-range interactions are highly prevalent within the compact Drosophila genome. Each enhancer contacts multiple enhancers, and promoters with similar expression, suggesting a role in their co-regulation. Notably, most interactions appear unchanged between tissue context and across development, arising before gene activation, and are frequently associated with paused RNA polymerase. Our results indicate that the general topology governing enhancer contacts is conserved from flies to humans and suggest that transcription initiates from preformed enhancer-promoter loops through release of paused polymerase.

To investigate the mechanism that drives dramatic mistargeting of active chromatin in NUT midline carcinoma (NMC), we have identified protein interactions unique to the BRD4-NUT fusion oncoprotein compared with wild-type BRD4. Using cross-linking, affinity purification, and mass spectrometry, we identified the EP300 acetyltransferase as uniquely associated with BRD4 through the NUT fusion in both NMC and non-NMC cell types. We also discovered ZNF532 associated with BRD4-NUT in NMC patient cells but not detectable in 293T cells. EP300 and ZNF532 are both implicated in feed-forward regulatory loops leading to propagation of the oncogenic chromatin complex in BRD4-NUT patient cells. Adding key functional significance to our biochemical findings, we independently discovered a ZNF532-NUT translocation fusion in a newly diagnosed NMC patient. ChIP sequencing of the major players NUT, ZNF532, BRD4, EP300, and H3K27ac revealed the formation of ZNF532-NUT-associated hyperacetylated megadomains, distinctly localized but otherwise analogous to those found in BRD4-NUT patient cells. Our results support a model in which NMC is dependent on ectopic NUT-mediated interactions between EP300 and components of BRD4 regulatory complexes, leading to a cascade of misregulation.BioTAP-XL | hyperacetylation | ZNF532-NUT | topological domains | BRD4

Developmental gene expression is tightly regulated through enhancer elements, which initiate dynamic spatio–temporal expression, and Polycomb response elements (PREs), which maintain stable gene silencing. These two cis-regulatory functions are thought to operate through distinct dedicated elements. By examining the occupancy of the Drosophila pleiohomeotic repressive complex (PhoRC) during embryogenesis, we revealed extensive co-occupancy at developmental enhancers. Using an established in vivo assay for PRE activity, we demonstrated that a subset of characterized developmental enhancers can function as PREs, silencing transcription in a Polycomb-dependent manner. Conversely, some classic Drosophila PREs can function as developmental enhancers in vivo, activating spatio–temporal expression. This study therefore uncovers elements with dual function: activating transcription in some cells (enhancers) while stably maintaining transcriptional silencing in others (PREs). Given that enhancers initiate spatio–temporal gene expression, reuse of the same elements by the Polycomb group (PcG) system may help fine-tune gene expression and ensure the timely maintenance of cell identities.

Motivation: Circularized Chromosome Conformation Capture (4C) is a powerful technique for studying the spatial interactions of a specific genomic region called the ‘viewpoint’ with the rest of the genome, both in a single condition or comparing different experimental conditions or cell types. Observed ligation frequencies typically show a strong, regular dependence on genomic distance from the viewpoint, on top of which specific interaction peaks are superimposed. Here, we address the computational task to find these specific peaks and to detect changes between different biological conditions.Results: We model the overall trend of decreasing interaction frequency with genomic distance by fitting a smooth monotonically decreasing function to suitably transformed count data. Based on the fit, z-scores are calculated from the residuals, and high z-scores are interpreted as peaks providing evidence for specific interactions. To compare different conditions, we normalize fragment counts between samples, and call for differential contact frequencies using the statistical method DESeq2 adapted from RNA-Seq analysis.Availability and implementation: A full end-to-end analysis pipeline is implemented in the R package FourCSeq available at www.bioconductor.org.Contact: felix.klein@embl.de or whuber@embl.deSupplementary information: Supplementary data are available at Bioinformatics online.

RUNX1 is known to be an essential transcription factor for generating hematopoietic stem cells (HSC), but much less is known about its role in the downstream process of hematopoietic differentiation. RUNX1 has been shown to be part of a large transcription factor complex, together with LDB1, GATA1, TAL1, and ETO2 (N. Meier et al., Development 133 :4913–4923, 2006) in erythroid cells. We used a tagging strategy to show that RUNX1 interacts with two novel protein partners, LSD1 and MYEF2, in erythroid cells. MYEF2 is bound in undifferentiated cells and is lost upon differentiation, whereas LSD1 is bound in differentiated cells. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) and microarray expression analysis were used to show that RUNX1 binds approximately 9,000 target sites in erythroid cells and is primarily active in the undifferentiated state. Functional analysis shows that a subset of the target genes is suppressed by RUNX1 via the newly identified partner MYEF2. Knockdown of Myef2 expression in developing zebrafish results in a reduced number of HSC.