Different enhancer classes in Drosophila bind distinct architectural proteins and mediate unique chromatin interactions and 3D architecture

New Drosophila promoter-associated architectural protein Mzfp1 interacts with CP190 and is required for housekeeping gene expression and insulator activity

In Drosophila, a group of zinc finger architectural proteins recruits the CP190 protein to the chromatin, an interaction that is essential for the functional activity of promoters and insulators. In this study, we describe a new architectural C2H2 protein called Madf and Zinc-Finger Protein 1 (Mzfp1) that interacts with CP190. Mzfp1 has an unusual structure that includes six C2H2 domains organized in a C-terminal cluster and two tandem MADF domains. Mzfp1 predominantly binds to housekeeping gene promoters located in both euchromatin and heterochromatin genome regions. In vivo mutagenesis studies showed that Mzfp1 is an essential protein, and both MADF domains and the CP190 interaction region are required for its functional activity. The C2H2 cluster is sufficient for the specific binding of Mzfp1 to regulatory elements, while the second MADF domain is required for Mzfp1 recruitment to heterochromatin. Mzfp1 binds to the proximal part of the Fub boundary that separates regulatory domains of the Ubx and abd-A genes in the Bithorax complex. Mzfp1 participates in Fub functions in cooperation with the architectural proteins Pita and Su(Hw). Thus, Mzfp1 is a new architectural C2H2 protein involved in the organization of active promoters and insulators in Drosophila.

Drosophila TET acts with PRC1 to activate gene expression independently of its catalytic activity

Enzymes of the ten-eleven translocation (TET) family play a key role in the regulation of gene expression by oxidizing 5-methylcytosine (5mC), a prominent epigenetic mark in many species. Yet, TET proteins also have less characterized noncanonical modes of action, notably in Drosophila, whose genome is devoid of 5mC. Here, we show that Drosophila TET activates the expression of genes required for larval central nervous system (CNS) development mainly in a catalytic-independent manner. Genome-wide profiling shows that TET is recruited to enhancer and promoter regions bound by Polycomb group complex (PcG) proteins. We found that TET interacts and colocalizes on chromatin preferentially with Polycomb repressor complex 1 (PRC1) rather than PRC2. Furthermore, PRC1 but not PRC2 is required for the activation of TET target genes. Last, our results suggest that TET and PRC1 binding to activated genes is interdependent. These data highlight the importance of TET noncatalytic function and the role of PRC1 for gene activation in the Drosophila larval CNS.

Mechanistic drivers of chromatin organization into compartments

The human genome is not just a simple string of DNA, it is a complex and dynamic entity intricately folded within the cell's nucleus. This three-dimensional organization of chromatin, the combination of DNA and proteins in the nucleus, is crucial for many biological processes and has been prominently studied for its intricate relationship to gene expression. Indeed, the transcriptional machinery does not operate in isolation but interacts intimately with the folded chromatin structure. Techniques for chromatin conformation capture, including genome-wide sequencing approaches, have revealed key organizational features of chromatin, such as the formation of loops by CCCTC-binding factor (CTCF) and the division of loci into chromatin compartments. While much of the recent research and reviews have focused on CTCF loops, we discuss several new revelations that have emerged concerning chromatin compartments, with a particular focus on what is known about mechanistic drivers of compartmentalization. These insights challenge the traditional views of chromatin organization and reveal the complexity behind the formation and maintenance of chromatin compartments.

The zinc-finger protein Z4 cooperates with condensin II to regulate somatic chromosome pairing and 3D chromatin organization

Chromosome pairing constitutes an important level of genome organization, yet the mechanisms that regulate pairing in somatic cells and the impact on 3D chromatin organization are still poorly understood. Here, we address these questions in Drosophila, an organism with robust somatic pairing. In Drosophila, pairing preferentially occurs at loci consisting of numerous architectural protein binding sites (APBSs), suggesting a role of architectural proteins (APs) in pairing regulation. Amongst these, the anti-pairing function of the condensin II subunit CAP-H2 is well established. However, the factors that regulate CAP-H2 localization and action at APBSs remain largely unknown. Here, we identify two factors that control CAP-H2 occupancy at APBSs and, therefore, regulate pairing. We show that Z4, interacts with CAP-H2 and is required for its localization at APBSs. We also show that hyperosmotic cellular stress induces fast and reversible unpairing in a Z4/CAP-H2 dependent manner. Moreover, by combining the opposite effects of Z4 depletion and osmostress, we show that pairing correlates with the strength of intrachromosomal 3D interactions, such as active (A) compartment interactions, intragenic gene-loops, and polycomb (Pc)-mediated chromatin loops. Altogether, our results reveal new players in CAP-H2-mediated pairing regulation and the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions.

HiCrayon reveals distinct layers of multi-state 3D chromatin organization

The co-visualization of chromatin conformation with 1D 'omics data is key to the multi-omics driven data analysis of 3D genome organization. Chromatin contact maps are often shown as 2D heatmaps and visually compared to 1D genomic data by simple juxtaposition. While common, this strategy is imprecise, placing the onus on the reader to align features with each other. To remedy this, we developed HiCrayon, an interactive tool that facilitates the integration of 3D chromatin organization maps and 1D datasets. This visualization method integrates data from genomic assays directly into the chromatin contact map by coloring interactions according to 1D signal. HiCrayon is implemented using R shiny and python to create a graphical user interface (GUI) application, available in both web or containerized format to promote accessibility. HiCrayon is implemented in R, and includes a graphical user interface (GUI), as well as a slimmed-down web-based version that lets users quickly produce publication-ready images. We demonstrate the utility of HiCrayon in visualizing the effectiveness of compartment calling and the relationship between ChIP-seq and various features of chromatin organization. We also demonstrate the improved visualization of other 3D genomic phenomena, such as differences between loops associated with CTCF/cohesin vs. those associated with H3K27ac. We then demonstrate HiCrayon's visualization of organizational changes that occur during differentiation and use HiCrayon to detect compartment patterns that cannot be assigned to either A or B compartments, revealing a distinct 3rd chromatin compartment. Overall, we demonstrate the utility of co-visualizing 2D chromatin conformation with 1D genomic signals within the same matrix to reveal fundamental aspects of genome organization. Local version: https://github.com/JRowleyLab/HiCrayon Web version: https://jrowleylab.com/HiCrayon.

Crosstalk between chromatin and Shavenbaby defines transcriptional output along the Drosophila intestinal stem cell lineage

The transcription factor Shavenbaby (Svb), the only member of the OvoL family in Drosophila, controls the fate of various epithelial embryonic cells and adult stem cells. Post-translational modification of Svb produces two protein isoforms, Svb-ACT and Svb-REP, which promote adult intestinal stem cell renewal or differentiation, respectively. To define Svb mode of action, we used engineered cell lines and develop an unbiased method to identify Svb target genes across different contexts. Within a given cell type, Svb-ACT and Svb-REP antagonistically regulate the expression of a set of target genes, binding specific enhancers whose accessibility is constrained by chromatin landscape. Reciprocally, Svb-REP can influence local chromatin marks of active enhancers to help repressing target genes. Along the intestinal lineage, the set of Svb target genes progressively changes, together with chromatin accessibility. We propose that Svb-ACT-to-REP transition promotes enterocyte differentiation of intestinal stem cells through direct gene regulation and chromatin remodeling.

Forecasting histone methylation by Polycomb complexes with minute-scale precision

Animals use the Polycomb system to epigenetically repress developmental genes. The repression requires trimethylation of lysine 27 of histone H3 (H3K27me3) by Polycomb Repressive Complex 2 (PRC2), but the dynamics of this process is poorly understood. To bridge the gap, we developed a computational model that forecasts H3K27 methylation in Drosophila with high temporal resolution and spatial accuracy of contemporary experimental techniques. Using this model, we show that pools of methylated H3K27 in dividing cells are defined by the effective concentration of PRC2 and the replication frequency. We find that the allosteric stimulation by preexisting H3K27me3 makes PRC2 better in methylating developmental genes as opposed to indiscriminate methylation throughout the genome. Applied to Drosophila development, our model argues that, in this organism, the intergenerationally inherited H3K27me3 does not "survive" rapid cycles of embryonic chromatin replication and is unlikely to transmit the memory of epigenetic repression to the offspring. Our model is adaptable to other organisms, including mice and humans.

Genome-wide expression QTL mapping reveals the highly dynamic regulatory landscape of a major wheat pathogen

BACKGROUND

In agricultural ecosystems, outbreaks of diseases are frequent and pose a significant threat to food security. A successful pathogen undergoes a complex and well-timed sequence of regulatory changes to avoid detection by the host immune system; hence, well-tuned gene regulation is essential for survival. However, the extent to which the regulatory polymorphisms in a pathogen population provide an adaptive advantage is poorly understood.

RESULTS

We used Zymoseptoria tritici, one of the most important pathogens of wheat, to generate a genome-wide map of regulatory polymorphism governing gene expression. We investigated genome-wide transcription levels of 146 strains grown under nutrient starvation and performed expression quantitative trait loci (eQTL) mapping. We identified cis-eQTLs for 65.3% of all genes and the majority of all eQTL loci are within 2kb upstream and downstream of the transcription start site (TSS). We also show that polymorphism in different gene elements contributes disproportionally to gene expression variation. Investigating regulatory polymorphism in gene categories, we found an enrichment of regulatory variants for genes predicted to be important for fungal pathogenesis but with comparatively low effect size, suggesting a separate layer of gene regulation involving epigenetics. We also show that previously reported trait-associated SNPs in pathogen populations are frequently cis-regulatory variants of neighboring genes with implications for the trait architecture.

CONCLUSIONS

Overall, our study provides extensive evidence that single populations segregate large-scale regulatory variation and are likely to fuel rapid adaptation to resistant hosts and environmental change.

The N-Terminal Part of Drosophila CP190 Is a Platform for Interaction with Multiple Architectural Proteins

CP190 is a co-factor in many Drosophila architectural proteins, being involved in the formation of active promoters and insulators. CP190 contains the N-terminal BTB/POZ (Broad-Complex, Tramtrack and Bric a brac/POxvirus and Zinc finger) domain and adjacent conserved regions involved in protein interactions. Here, we examined the functional roles of these domains of CP190 in vivo. The best-characterized architectural proteins with insulator functions, Pita, Su(Hw), and dCTCF, interacted predominantly with the BTB domain of CP190. Due to the difficulty of mutating the BTB domain, we obtained a transgenic line expressing a chimeric CP190 with the BTB domain of the human protein Kaiso. Another group of architectural proteins, M1BP, Opbp, and ZIPIC, interacted with one or both of the highly conserved regions in the N-terminal part of CP190. Transgenic lines of D. melanogaster expressing CP190 mutants with a deletion of each of these domains were obtained. The results showed that these mutant proteins only partially compensated for the functions of CP190, weakly binding to selective chromatin sites. Further analysis confirmed the essential role of these domains in recruitment to regulatory regions associated with architectural proteins. We also found that the N-terminal of CP190 was sufficient for recruiting Z4 and Chromator proteins and successfully achieving chromatin opening. Taken together, our results and the results of previous studies showed that the N-terminal region of CP190 is a platform for simultaneous interaction with various DNA-binding architectural proteins and transcription complexes.

3D chromatin interactions involving Drosophila insulators are infrequent but preferential and arise before TADs and transcription

In mammals, insulators contribute to the regulation of loop extrusion to organize chromatin into topologically associating domains. In Drosophila the role of insulators in 3D genome organization is, however, under current debate. Here, we addressed this question by combining bioinformatics analysis and multiplexed chromatin imaging. We describe a class of Drosophila insulators enriched at regions forming preferential chromatin interactions genome-wide. Notably, most of these 3D interactions do not involve TAD borders. Multiplexed imaging shows that these interactions occur infrequently, and only rarely involve multiple genomic regions coalescing together in space in single cells. Finally, we show that non-border preferential 3D interactions enriched in this class of insulators are present before TADs and transcription during Drosophila development. Our results are inconsistent with insulators forming stable hubs in single cells, and instead suggest that they fine-tune existing 3D chromatin interactions, providing an additional regulatory layer for transcriptional regulation.

The CLAMP GA-binding transcription factor regulates heat stress-induced transcriptional repression by associating with 3D loop anchors

In order to survive when exposed to heat stress (HS), organisms activate stress response genes and repress constitutive gene expression to prevent the accumulation of potentially toxic RNA and protein products. Although many studies have elucidated the mechanisms that drive HS-induced activation of stress response genes across species, little is known about repression mechanisms or how genes are targeted for activation versus repression context-specifically. The mechanisms of heat stress-regulated activation have been well-studied in Drosophila, in which the GA-binding transcription factor GAF is important for activating genes upon heat stress. Here, we show that a functionally distinct GA-binding transcription factor (TF) protein, CLAMP (Chromatin-linked adaptor for MSL complex proteins), is essential for repressing constitutive genes upon heat stress but not activation of the canonical heat stress pathway. HS induces loss of CLAMP-associated 3D chromatin loop anchors associated with different combinations of GA-binding TFs prior to HS if a gene becomes repressed versus activated. Overall, we demonstrate that CLAMP promotes repression of constitutive genes upon HS, and repression and activation are associated with the loss of CLAMP-associated 3D chromatin loops bound by different combinations of GA-binding TFs.

Underlying causes for prevalent false positives and false negatives in STARR-seq data

Self-transcribing active regulatory region sequencing (STARR-seq) and its variants have been widely used to characterize enhancers. However, it has been reported that up to 87% of STARR-seq peaks are located in repressive chromatin and are not functional in the tested cells. While some of the STARR-seq peaks in repressive chromatin might be active in other cell/tissue types, some others might be false positives. Meanwhile, many active enhancers may not be identified by the current STARR-seq methods. Although methods have been proposed to mitigate systematic errors caused by the use of plasmid vectors, the artifacts due to the intrinsic limitations of current STARR-seq methods are still prevalent and the underlying causes are not fully understood. Based on predicted cis-regulatory modules (CRMs) and non-CRMs in the human genome as well as predicted active CRMs and non-active CRMs in a few human cell lines/tissues with STARR-seq data available, we reveal prevalent false positives and false negatives in STARR-seq peaks generated by major variants of STARR-seq methods and possible underlying causes. Our results will help design strategies to improve STARR-seq methods and interpret the results.

CTCF and R-loops are boundaries of cohesin-mediated DNA looping

Cohesin and CCCTC-binding factor (CTCF) are key regulatory proteins of three-dimensional (3D) genome organization. Cohesin extrudes DNA loops that are anchored by CTCF in a polar orientation. Here, we present direct evidence that CTCF binding polarity controls cohesin-mediated DNA looping. Using single-molecule imaging, we demonstrate that a critical N-terminal motif of CTCF blocks cohesin translocation and DNA looping. The cryo-EM structure of the cohesin-CTCF complex reveals that this CTCF motif ahead of zinc fingers can only reach its binding site on the STAG1 cohesin subunit when the N terminus of CTCF faces cohesin. Remarkably, a C-terminally oriented CTCF accelerates DNA compaction by cohesin. DNA-bound Cas9 and Cas12a ribonucleoproteins are also polar cohesin barriers, indicating that stalling may be intrinsic to cohesin itself. Finally, we show that RNA-DNA hybrids (R-loops) block cohesin-mediated DNA compaction in vitro and are enriched with cohesin subunits in vivo, likely forming TAD boundaries.

Chromatin alternates between A and B compartments at kilobase scale for subgenic organization

Nuclear compartments are prominent features of 3D chromatin organization, but sequencing depth limitations have impeded investigation at ultra fine-scale. CTCF loops are generally studied at a finer scale, but the impact of looping on proximal interactions remains enigmatic. Here, we critically examine nuclear compartments and CTCF loop-proximal interactions using a combination of in situ Hi-C at unparalleled depth, algorithm development, and biophysical modeling. Producing a large Hi-C map with 33 billion contacts in conjunction with an algorithm for performing principal component analysis on sparse, super massive matrices (POSSUMM), we resolve compartments to 500 bp. Our results demonstrate that essentially all active promoters and distal enhancers localize in the A compartment, even when flanking sequences do not. Furthermore, we find that the TSS and TTS of paused genes are often segregated into separate compartments. We then identify diffuse interactions that radiate from CTCF loop anchors, which correlate with strong enhancer-promoter interactions and proximal transcription. We also find that these diffuse interactions depend on CTCF's RNA binding domains. In this work, we demonstrate features of fine-scale chromatin organization consistent with a revised model in which compartments are more precise than commonly thought while CTCF loops are more protracted.

Study of the in Vivo Functional Role of Mutations in the BTB Domain of the CP190 Protein of Drosophila melanogaster

The Drosophila transcription factor СР190 is one of the key proteins that determine the activity of housekeeping gene promoters and insulators. CP190 has an N-terminal BTB domain that allows for dimerization. Many of known Drosophila architectural proteins interact with the hydrophobic peptide-binding groove in the BTB domain, which is presumably a mechanisms for recruiting CP190 to regulatory elements. To study the role of the BTB domain in the interaction with architectural proteins, we obtained transgenic flies expressing CP190 variants with mutations in the peptide-binding groove, which disrupts their interaction with architectural proteins. As a result of the studies, it was found that mutations in the BTB domain do not affect binding of the CP190 protein to polytene chromosomes. Thus, our studies confirm the previously obtained data that CP190 is recruited to regulatory elements by several transcription factors interacting, in addition to BTB, with other CP190 domains.

Mechanisms of Interaction between Enhancers and Promoters in Three Drosophila Model Systems

In higher eukaryotes, the regulation of developmental gene expression is determined by enhancers, which are often located at a large distance from the promoters they regulate. Therefore, the architecture of chromosomes and the mechanisms that determine the functional interaction between enhancers and promoters are of decisive importance in the development of organisms. Mammals and the model animal Drosophila have homologous key architectural proteins and similar mechanisms in the organization of chromosome architecture. This review describes the current progress in understanding the mechanisms of the formation and regulation of long-range interactions between enhancers and promoters at three well-studied key regulatory loci in Drosophila.

A Comparison of Topologically Associating Domain Callers Based on Hi-C Data

Topologically associating domains (TADs) are local chromatin interaction domains, which have been shown to play an important role in gene expression regulation. TADs were originally discovered in the investigation of 3D genome organization based on High-throughput Chromosome Conformation Capture (Hi-C) data. Continuous considerable efforts have been dedicated to developing methods for detecting TADs from Hi-C data. Different computational methods for TADs identification vary in their assumptions and criteria in calling TADs. As a consequence, the TADs called by these methods differ in their similarities and biological features they are enriched in. In this work, we performed a systematic comparison of twenty-six TAD callers. We first compared the TADs and gaps between adjacent TADs across different methods, resolutions, and sequencing depths. We then assessed the quality of TADs and TAD boundaries according to three criteria: the decay of contact frequencies over the genomic distance, enrichment and depletion of regulatory elements around TAD boundaries, and reproducibility of TADs and TAD boundaries in replicate samples. Last, due to the lack of a gold standard of TADs, we also evaluated the performance of the methods on synthetic datasets. We discussed the key principles of TAD callers, and pinpointed current situation in the detection of TADs. We provide a concise, comprehensive, and systematic framework for evaluating the performance of TAD callers, and expect our work will provide useful guidance in choosing suitable approaches for the detection and evaluation of TADs.

Mechanisms of enhancer-promoter communication and chromosomal architecture in mammals and Drosophila

The spatial organization of chromosomes is involved in regulating the majority of intranuclear processes in higher eukaryotes, including gene expression. Drosophila was used as a model to discover many transcription factors whose homologs play a key role in regulation of gene expression in mammals. According to modern views, a cohesin complex mostly determines the architecture of mammalian chromosomes by forming chromatin loops on anchors created by the CTCF DNA-binding architectural protein. The role of the cohesin complex in chromosome architecture is poorly understood in Drosophila, and CTCF is merely one of many Drosophila architectural proteins with a proven potential to organize specific long-range interactions between regulatory elements in the genome. The review compares the mechanisms responsible for long-range interactions and chromosome architecture between mammals and Drosophila.

Unveiling dynamic enhancer–promoter interactions in Drosophila melanogaster

Proper enhancer–promoter interactions are essential to maintaining specific transcriptional patterns and preventing ectopic gene expression. Drosophila is an ideal model organism to study transcriptional regulation due to extensively characterized regulatory regions and the ease of implementing new genetic and molecular techniques for quantitative analysis. The mechanisms of enhancer–promoter interactions have been investigated over a range of length scales. At a DNA level, compositions of both enhancer and promoter sequences affect transcriptional dynamics, including duration, amplitude, and frequency of transcriptional bursting. 3D chromatin topology is also important for proper enhancer–promoter contacts. By working competitively or cooperatively with one another, multiple, simultaneous enhancer–enhancer, enhancer–promoter, and promoter–promoter interactions often occur to maintain appropriate levels of mRNAs. For some long-range enhancer–promoter interactions, extra regulatory elements like insulators and tethering elements are required to promote proper interactions while blocking aberrant ones. This review provides an overview of our current understanding of the mechanism of enhancer–promoter interactions and how perturbations of such interactions affect transcription and subsequent physiological outcomes.

A global high-density chromatin interaction network reveals functional long-range and trans-chromosomal relationships

BACKGROUND

Chromatin contacts are essential for gene-expression regulation; however, obtaining a high-resolution genome-wide chromatin contact map is still prohibitively expensive owing to large genome sizes and the quadratic scale of pairwise data. Chromosome conformation capture (3C)-based methods such as Hi-C have been extensively used to obtain chromatin contacts. However, since the sparsity of these maps increases with an increase in genomic distance between contacts, long-range or trans-chromatin contacts are especially challenging to sample.

RESULTS

Here, we create a high-density reference genome-wide chromatin contact map using a meta-analytic approach. We integrate 3600 human, 6700 mouse, and 500 fly Hi-C experiments to create species-specific meta-Hi-C chromatin contact maps with 304 billion, 193 billion, and 19 billion contacts in respective species. We validate that meta-Hi-C contact maps are uniquely powered to capture functional chromatin contacts in both cis and trans. We find that while individual dataset Hi-C networks are largely unable to predict any long-range coexpression (median 0.54 AUC), meta-Hi-C networks perform comparably in both cis and trans (0.65 AUC vs 0.64 AUC). Similarly, for long-range expression quantitative trait loci (eQTL), meta-Hi-C contacts outperform all individual Hi-C experiments, providing an improvement over the conventionally used linear genomic distance-based association. Assessing between species, we find patterns of chromatin contact conservation in both cis and trans and strong associations with coexpression even in species for which Hi-C data is lacking.

CONCLUSIONS

We have generated an integrated chromatin interaction network which complements a large number of methodological and analytic approaches focused on improved specificity or interpretation. This high-depth "super-experiment" is surprisingly powerful in capturing long-range functional relationships of chromatin interactions, which are now able to predict coexpression, eQTLs, and cross-species relationships. The meta-Hi-C networks are available at https://labshare.cshl.edu/shares/gillislab/resource/HiC/ .

The Drosophila BEAF insulator protein interacts with the polybromo subunit of the PBAP chromatin remodeling complex

The Drosophila Boundary Element-Associated Factor of 32 kDa (BEAF) binds in promoter regions of a few thousand mostly housekeeping genes. BEAF is implicated in both chromatin domain boundary activity and promoter function, although molecular mechanisms remain elusive. Here, we show that BEAF physically interacts with the polybromo subunit (Pbro) of PBAP, a SWI/SNF-class chromatin remodeling complex. BEAF also shows genetic interactions with Pbro and other PBAP subunits. We examine the effect of this interaction on gene expression and chromatin structure using precision run-on sequencing and micrococcal nuclease sequencing after RNAi-mediated knockdown in cultured S2 cells. Our results are consistent with the interaction playing a subtle role in gene activation. Fewer than 5% of BEAF-associated genes were significantly affected after BEAF knockdown. Most were downregulated, accompanied by fill-in of the promoter nucleosome-depleted region and a slight upstream shift of the +1 nucleosome. Pbro knockdown caused downregulation of several hundred genes and showed a correlation with BEAF knockdown but a better correlation with promoter-proximal GAGA factor binding. Micrococcal nuclease sequencing supports that BEAF binds near housekeeping gene promoters while Pbro is more important at regulated genes. Yet there is a similar general but slight reduction of promoter-proximal pausing by RNA polymerase II and increase in nucleosome-depleted region nucleosome occupancy after knockdown of either protein. We discuss the possibility of redundant factors keeping BEAF-associated promoters active and masking the role of interactions between BEAF and the Pbro subunit of PBAP in S2 cells. We identify Facilitates Chromatin Transcription (FACT) and Nucleosome Remodeling Factor (NURF) as candidate redundant factors.

RNA Polymerase II “Pause” Prepares Promoters for Upcoming Transcription during Drosophila Development

According to previous studies, during Drosophila embryogenesis, the recruitment of RNA polymerase II precedes active gene transcription. This work is aimed at exploring whether this mechanism is used during Drosophila metamorphosis. In addition, the composition of the RNA polymerase II “paused” complexes associated with promoters at different developmental stages are described in detail. For this purpose, we performed ChIP-Seq analysis using antibodies for various modifications of RNA polymerase II (total, Pol II CTD Ser5P, and Pol II CTD Ser2P) as well as for subunits of the NELF, DSIF, and PAF complexes and Brd4/Fs(1)h that control transcription elongation. We found that during metamorphosis, similar to mid-embryogenesis, the promoters were bound by RNA polymerase II in the “paused” state, preparing for activation at later stages of development. During mid-embryogenesis, RNA polymerase II in a “pause” state was phosphorylated at Ser5 and Ser2 of Pol II CTD and bound the NELF, DSIF, and PAF complexes, but not Brd4/Fs(1)h. During metamorphosis, the “paused” RNA polymerase II complex included Brd4/Fs(1)h in addition to NELF, DSIF, and PAF. The RNA polymerase II in this complex was phosphorylated at Ser5 of Pol II CTD, but not at Ser2. These results indicate that, during mid-embryogenesis, RNA polymerase II stalls in the “post-pause” state, being phosphorylated at Ser2 of Pol II CTD (after the stage of p-TEFb action). During metamorphosis, the “pause” mechanism is closer to classical promoter-proximal pausing and is characterized by a low level of Pol II CTD Ser2P.

3D chromatin remodeling potentiates transcriptional programs driving cell invasion

The contribution of deregulated chromatin architecture, including topologically associated domains (TADs), to cancer progression remains ambiguous. CCCTC-binding factor (CTCF) is a central regulator of higher-order chromatin structure that undergoes copy number loss in over half of all breast cancers, but the impact of this defect on epigenetic programming and chromatin architecture remains unclear. We find that under physiological conditions, CTCF organizes subTADs to limit the expression of oncogenic pathways, including phosphatidylinositol 3-kinase (PI3K) and cell adhesion networks. Loss of a single CTCF allele potentiates cell invasion through compromised chromatin insulation and a reorganization of chromatin architecture and histone programming that facilitates de novo promoter-enhancer contacts. However, this change in the higher-order chromatin landscape leads to a vulnerability to inhibitors of mTOR. These data support a model whereby subTAD reorganization drives both modification of histones at de novo enhancer-promoter contacts and transcriptional up-regulation of oncogenic transcriptional networks.

Application of the 3C Method to Study the Developmental Genes in Drosophila Larvae

A transition from one developmental stage to another is accompanied by activation of developmental programs and corresponding gene ensembles. Changes in the spatial conformation of the corresponding loci are associated with this activation and can be investigated with the help of the Chromosome Conformation Capture (3C) methodology. Application of 3C to specific developmental stages is a sophisticated task. Here, we describe the use of the 3C method to study the spatial organization of developmental loci in Drosophila larvae. We critically analyzed the existing protocols and offered our own solutions and the optimized protocol to overcome limitations. To demonstrate the efficiency of our procedure, we studied the spatial organization of the developmental locus Dad in 3rd instar Drosophila larvae. Differences in locus conformation were found between embryonic cells and living wild-type larvae. We also observed the establishment of novel regulatory interactions in the presence of an adjacent transgene upon activation of its expression in larvae. Our work fills the gap in the application of the 3C method to Drosophila larvae and provides a useful guide for establishing 3C on an animal model.

ZipHiC: a novel Bayesian framework to identify enriched interactions and experimental biases in Hi-C data

Motivation

Several computational and statistical methods have been developed to analyze data generated through the 3C-based methods, especially the Hi-C. Most of the existing methods do not account for dependency in Hi-C data.

Results

Here, we present ZipHiC, a novel statistical method to explore Hi-C data focusing on the detection of enriched contacts. ZipHiC implements a Bayesian method based on a hidden Markov random field (HMRF) model and the Approximate Bayesian Computation (ABC) to detect interactions in two-dimensional space based on a Hi-C contact frequency matrix. ZipHiC uses data on the sources of biases related to the contact frequency matrix, allows borrowing information from neighbours using the Potts model and improves computation speed using the ABC model. In addition to outperforming existing tools on both simulated and real data, our model also provides insights into different sources of biases that affects Hi-C data. We show that some datasets display higher biases from DNA accessibility or Transposable Elements content. Furthermore, our analysis in Drosophila melanogaster showed that approximately half of the detected significant interactions connect promoters with other parts of the genome indicating a functional biological role. Finally, we found that the micro-C datasets display higher biases from DNA accessibility compared to a similar Hi-C experiment, but this can be corrected by ZipHiC.

Availability and implementation

The R scripts are available at https://github.com/igosungithub/HMRFHiC.git.

Supplementary information

Supplementary data are available at Bioinformatics online.

3DGenBench: a web-server to benchmark computational models for 3D Genomics

Modeling 3D genome organisation has been booming in the last years thanks to the availability of experimental datasets of genomic contacts. However, the field is currently missing the standardisation of methods and metrics to compare predictions and experiments. We present 3DGenBench, a web server available at https://inc-cost.eu/benchmarking/, that allows benchmarking computational models of 3D Genomics. The benchmark is performed using a manually curated dataset of 39 capture Hi-C profiles in wild type and genome-edited mouse cells, and five genome-wide Hi-C profiles in human, mouse, and Drosophila cells. 3DGenBench performs two kinds of analysis, each supplied with a specific scoring module that compares predictions of a computational method to experimental data using several metrics. With 3DGenBench, the user obtains model performance scores, allowing an unbiased comparison with other models. 3DGenBench aims to become a reference web server to test new 3D genomics models and is conceived as an evolving platform where new types of analysis will be implemented in the future.

3D chromatin architecture and transcription regulation in cancer

Chromatin has distinct three-dimensional (3D) architectures important in key biological processes, such as cell cycle, replication, differentiation, and transcription regulation. In turn, aberrant 3D structures play a vital role in developing abnormalities and diseases such as cancer. This review discusses key 3D chromatin structures (topologically associating domain, lamina-associated domain, and enhancer-promoter interactions) and corresponding structural protein elements mediating 3D chromatin interactions [CCCTC-binding factor, polycomb group protein, cohesin, and Brother of the Regulator of Imprinted Sites (BORIS) protein] with a highlight of their associations with cancer. We also summarise the recent development of technologies and bioinformatics approaches to study the 3D chromatin interactions in gene expression regulation, including crosslinking and proximity ligation methods in the bulk cell population (ChIA-PET and HiChIP) or single-molecule resolution (ChIA-drop), and methods other than proximity ligation, such as GAM, SPRITE, and super-resolution microscopy techniques.

The role of insulators and transcription in 3D chromatin organization of flies

The DNA in many organisms, including humans, is shown to be organized in topologically associating domains (TADs). In Drosophila, several architectural proteins are enriched at TAD borders, but it is still unclear whether these proteins play a functional role in the formation and maintenance of TADs. Here, we show that depletion of BEAF-32, Cp190, Chro, and Dref leads to changes in TAD organization and chromatin loops. Their depletion predominantly affects TAD borders located in regions moderately enriched in repressive modifications and depleted in active ones, whereas TAD borders located in euchromatin are resilient to these knockdowns. Furthermore, transcriptomic data has revealed hundreds of genes displaying differential expression in these knockdowns and showed that the majority of differentially expressed genes are located within reorganized TADs. Our work identifies a novel and functional role for architectural proteins at TAD borders in Drosophila and a link between TAD reorganization and subsequent changes in gene expression.

Implications of Dosage Deficiencies in CTCF and Cohesin on Genome Organization, Gene Expression, and Human Neurodevelopment

Properly organizing DNA within the nucleus is critical to ensure normal downstream nuclear functions. CTCF and cohesin act as major architectural proteins, working in concert to generate thousands of high-intensity chromatin loops. Due to their central role in loop formation, a massive research effort has been dedicated to investigating the mechanism by which CTCF and cohesin create these loops. Recent results lead to questioning the direct impact of CTCF loops on gene expression. Additionally, results of controlled depletion experiments in cell lines has indicated that genome architecture may be somewhat resistant to incomplete deficiencies in CTCF or cohesin. However, heterozygous human genetic deficiencies in CTCF and cohesin have illustrated the importance of their dosage in genome architecture, cellular processes, animal behavior, and disease phenotypes. Thus, the importance of considering CTCF or cohesin levels is especially made clear by these heterozygous germline variants that characterize genetic syndromes, which are increasingly recognized in clinical practice. Defined primarily by developmental delay and intellectual disability, the phenotypes of CTCF and cohesin deficiency illustrate the importance of architectural proteins particularly in neurodevelopment. We discuss the distinct roles of CTCF and cohesin in forming chromatin loops, highlight the major role that dosage of each protein plays in the amplitude of observed effects on gene expression, and contrast these results to heterozygous mutation phenotypes in murine models and clinical patients. Insights highlighted by this comparison have implications for future research into these newly emerging genetic syndromes.

OUP accepted manuscript

Eukaryotic chromosomes are spatially segregated into topologically associating domains (TADs). Some TADs are attached to the nuclear lamina (NL) through lamina-associated domains (LADs). Here, we identified LADs and TADs at two stages of Drosophila spermatogenesis – in bam^Δ86 mutant testes which is the commonly used model of spermatogonia (SpG) and in larval testes mainly filled with spermatocytes (SpCs). We found that initiation of SpC-specific transcription correlates with promoters’ detachment from the NL and with local spatial insulation of adjacent regions. However, this insulation does not result in the partitioning of inactive TADs into sub-TADs. We also revealed an increased contact frequency between SpC-specific genes in SpCs implying their de novo gathering into transcription factories. In addition, we uncovered the specific X chromosome organization in the male germline. In SpG and SpCs, a single X chromosome is stronger associated with the NL than autosomes. Nevertheless, active chromatin regions in the X chromosome interact with each other more frequently than in autosomes. Moreover, despite the absence of dosage compensation complex in the male germline, randomly inserted SpG-specific reporter is expressed higher in the X chromosome than in autosomes, thus evidencing that non-canonical dosage compensation operates in SpG.

Mechanisms of CP190 Interaction with Architectural Proteins in Drosophila Melanogaster

Most of the known Drosophila architectural proteins interact with an important cofactor, CP190, that contains three domains (BTB, M, and D) that are involved in protein–protein interactions. The highly conserved N-terminal CP190 BTB domain forms a stable homodimer that interacts with unstructured regions in the three best-characterized architectural proteins: dCTCF, Su(Hw), and Pita. Here, we identified two new CP190 partners, CG4730 and CG31365, that interact with the BTB domain. The CP190 BTB resembles the previously characterized human BCL6 BTB domain, which uses its hydrophobic groove to specifically associate with unstructured regions of several transcriptional repressors. Using GST pull-down and yeast two-hybrid assays, we demonstrated that mutations in the hydrophobic groove strongly affect the affinity of CP190 BTB for the architectural proteins. In the yeast two-hybrid assay, we found that architectural proteins use various mechanisms to improve the efficiency of interaction with CP190. Pita and Su(Hw) have two unstructured regions that appear to simultaneously interact with hydrophobic grooves in the BTB dimer. In dCTCF and CG31365, two adjacent regions interact simultaneously with the hydrophobic groove of the BTB and the M domain of CP190. Finally, CG4730 interacts with the BTB, M, and D domains of CP190 simultaneously. These results suggest that architectural proteins use different mechanisms to increase the efficiency of interaction with CP190.

An explainable artificial intelligence approach for decoding the enhancer histone modifications code and identification of novel enhancers in Drosophila

BACKGROUND

Enhancers are non-coding regions of the genome that control the activity of target genes. Recent efforts to identify active enhancers experimentally and in silico have proven effective. While these tools can predict the locations of enhancers with a high degree of accuracy, the mechanisms underpinning the activity of enhancers are often unclear.

RESULTS

Using machine learning (ML) and a rule-based explainable artificial intelligence (XAI) model, we demonstrate that we can predict the location of known enhancers in Drosophila with a high degree of accuracy. Most importantly, we use the rules of the XAI model to provide insight into the underlying combinatorial histone modifications code of enhancers. In addition, we identified a large set of putative enhancers that display the same epigenetic signature as enhancers identified experimentally. These putative enhancers are enriched in nascent transcription, divergent transcription and have 3D contacts with promoters of transcribed genes. However, they display only intermediary enrichment of mediator and cohesin complexes compared to previously characterised active enhancers. We also found that 10-15% of the predicted enhancers display similar characteristics to super enhancers observed in other species.

CONCLUSIONS

Here, we applied an explainable AI model to predict enhancers with high accuracy. Most importantly, we identified that different combinations of epigenetic marks characterise different groups of enhancers. Finally, we discovered a large set of putative enhancers which display similar characteristics with previously characterised active enhancers.

Genomic organization of the autonomous regulatory domain of eyeless locus in Drosophila melanogaster

In Drosophila, expression of eyeless (ey) gene is restricted to the developing eyes and central nervous system. However, the flanking genes, myoglianin (myo), and bent (bt) have different temporal and spatial expression patterns as compared to the ey. How distinct regulation of ey is maintained is mostly unknown. Earlier, we have identified a boundary element intervening myo and ey genes (ME boundary) that prevents the crosstalk between the cis-regulatory elements of myo and ey genes. In the present study, we further searched for the cis-elements that define the domain of ey and maintain its expression pattern. We identify another boundary element between ey and bt, the EB boundary. The EB boundary separates the regulatory landscapes of ey and bt genes. The two boundaries, ME and EB, show a long-range interaction as well as interact with the nuclear architecture. This suggests functional autonomy of the ey locus and its insulation from differentially regulated flanking regions. We also identify a new Polycomb Response Element, the ey-PRE, within the ey domain. The expression state of the ey gene, once established during early development is likely to be maintained with the help of ey-PRE. Our study proposes a general regulatory mechanism by which a gene can be maintained in a functionally independent chromatin domain in gene-rich euchromatin.

Investigation of the Basic Steps in the Chromosome Conformation Capture Procedure

A constellation of chromosome conformation capture methods (С-methods) are an important tool for biochemical analysis of the spatial interactions between DNA regions that are separated in the primary sequence. All these methods are based on the long sequence of basic steps of treating cells, nuclei, chromatin, and finally DNA, thus representing a significant technical challenge. Here, we present an in-depth study of the basic steps in the chromatin conformation capture procedure (3С), which was performed using Drosophila Schneider 2 cells as a model. We investigated the steps of cell lysis, nuclei washing, nucleoplasm extraction, chromatin treatment with SDS/Triton X-100, restriction enzyme digestion, chromatin ligation, reversion of cross-links, DNA extraction, treatment of a 3C library with RNases, and purification of the 3C library. Several options were studied, and optimal conditions were found. Our work contributes to the understanding of the 3C basic steps and provides a useful guide to the 3C procedure.

Drosophila architectural protein CTCF is not essential for fly survival and is able to function independently of CP190

CTCF is the most likely ancestor of proteins that contain large clusters of C2H2 zinc finger domains (C2H2) and is conserved among most bilateral organisms. In mammals, CTCF functions as the main architectural protein involved in the organization of topology-associated domains (TADs). In vertebrates and Drosophila, CTCF is involved in the regulation of homeotic genes. Previously, it was found that null mutations in the dCTCF gene died as pharate adults, which failed to eclose from their pupal case, or shortly after hatching of adults. Here, we obtained several new null dCTCF mutations and found that the complete inactivation of dCTCF appears is limited mainly to phenotypic manifestations of the Abd-B gene and fertility of adult flies. Many modifiers that are not associated with an independent phenotypic manifestation can significantly enhance the expressivity of the null dCTCF mutations, indicating that other architectural proteins are able to functionally compensate for dCTCF inactivation in Drosophila. We also mapped the 715-735 aa region of dCTCF as being essential for the interaction with the BTB (Broad-Complex, Tramtrack, and Bric a brac) and microtubule-targeting (M) domains of the CP190 protein, which binds to many architectural proteins. However, the mutational analysis showed that the interaction with CP190 was not important for the functional activity of dCTCF in vivo.

3D genome organisation in Drosophila

Ever since Thomas Hunt Morgan's discovery of the chromosomal basis of inheritance by using Drosophila melanogaster as a model organism, the fruit fly has remained an essential model system in studies of genome biology, including chromatin organisation. Very much as in vertebrates, in Drosophila, the genome is organised in territories, compartments and topologically associating domains (TADs). However, these domains might be formed through a slightly different mechanism than in vertebrates due to the presence of a large and potentially redundant set of insulator proteins and the minor role of dCTCF in TAD boundary formation. Here, we review the different levels of chromatin organisation in Drosophila and discuss mechanisms and factors that might be involved in TAD formation. The dynamics of TADs and enhancer-promoter interactions in the context of transcription are covered in the light of currently conflicting results. Finally, we illustrate the value of polymer modelling approaches to infer the principles governing the three-dimensional organisation of the Drosophila genome.

M1BP cooperates with CP190 to activate transcription at TAD borders and promote chromatin insulator activity

Genome organization is driven by forces affecting transcriptional state, but the relationship between transcription and genome architecture remains unclear. Here, we identified the Drosophila transcription factor Motif 1 Binding Protein (M1BP) in physical association with the gypsy chromatin insulator core complex, including the universal insulator protein CP190. M1BP is required for enhancer-blocking and barrier activities of the gypsy insulator as well as its proper nuclear localization. Genome-wide, M1BP specifically colocalizes with CP190 at Motif 1-containing promoters, which are enriched at topologically associating domain (TAD) borders. M1BP facilitates CP190 chromatin binding at many shared sites and vice versa. Both factors promote Motif 1-dependent gene expression and transcription near TAD borders genome-wide. Finally, loss of M1BP reduces chromatin accessibility and increases both inter- and intra-TAD local genome compaction. Our results reveal physical and functional interaction between CP190 and M1BP to activate transcription at TAD borders and mediate chromatin insulator-dependent genome organization.

Transcriptional state plays a role in genome organization, however factors that link these processes are not well known. Here, the authors show Drosophila transcription factor Motif 1-binding protein (M1BP) interacts with the insulator protein CP190 to promote insulator function and activate Motif 1-dependent transcription at topologically associating domain (TAD) borders.

Principles of 3D compartmentalization of the human genome

Chromatin is organized in the nucleus via CTCF loops and compartmental domains. Here, we compare different cell types to identify distinct paradigms of compartmental domain formation in human tissues. We identify and quantify compartmental forces correlated with histone modifications characteristic of transcriptional activity and previously underappreciated roles for distinct compartmental domains correlated with the presence of H3K27me3 and H3K9me3, respectively. We present a computer simulation model capable of predicting compartmental organization based on the biochemical characteristics of independent chromatin features. Using this model, we show that the underlying forces responsible for compartmental domain formation in human cells are conserved and that the diverse compartmentalization patterns seen across cell types are due to differences in chromatin features. We extend these findings to Drosophila to suggest that the same principles are at work beyond humans. These results offer mechanistic insights into the fundamental forces driving the 3D organization of the genome.

Using high-resolution Hi-C data and computer simulations, Nichols and Corces show that compartments arise as a consequence of interactions among proteins that correlate with the presence of H3K27ac, H3K27me3, and H3K9me3, suggesting that human cells contain at least three distinct compartments. The same principles apply to other organisms.

Mechanisms of enhancer action: the known and the unknown

Differential gene expression mechanisms ensure cellular differentiation and plasticity to shape ontogenetic and phylogenetic diversity of cell types. A key regulator of differential gene expression programs are the enhancers, the gene-distal cis-regulatory sequences that govern spatiotemporal and quantitative expression dynamics of target genes. Enhancers are widely believed to physically contact the target promoters to effect transcriptional activation. However, our understanding of the full complement of regulatory proteins and the definitive mechanics of enhancer action is incomplete. Here, we review recent findings to present some emerging concepts on enhancer action and also outline a set of outstanding questions.

The Dm-Myb Oncoprotein Contributes to Insulator Function and Stabilizes Repressive H3K27me3 PcG Domains

Drosophila Myb (Dm-Myb) encodes a protein that plays a key role in regulation of mitotic phase genes. Here, we further refine its role in the context of a developing tissue as a potentiator of gene expression required for proper RNA polymerase II (RNA Pol II) function and efficient H3K4 methylation at promoters. In contrast to its role in gene activation, Myb is also required for repression of many genes, although no specific mechanism for this role has been proposed. We now reveal a critical role for Myb in contributing to insulator function, in part by promoting binding of insulator proteins BEAF-32 and CP190 and stabilizing H3K27me3 Polycomb-group (PcG) domains. In the absence of Myb, H3K27me3 is markedly reduced throughout the genome, leading to H3K4me3 spreading and gene derepression. Finally, Myb is enriched at boundaries that demarcate chromatin environments, including chromatin loop anchors. These results reveal functions of Myb that extend beyond transcriptional regulation.

Myb has been considered a transcriptional activator of primarily M phase genes. Here, Santana et al. show that Myb also contributes to insulator function, in part by promoting binding of insulator factors, and is required to stabilize repressive domains in the genome.

Subunits of the PBAP Chromatin Remodeler Are Capable of Mediating Enhancer-Driven Transcription in Drosophila

The chromatin remodeler SWI/SNF is an important participant in gene activation, functioning predominantly by opening the chromatin structure on promoters and enhancers. Here, we describe its novel mode of action in which SWI/SNF factors mediate the targeted action of an enhancer. We studied the functions of two signature subunits of PBAP subfamily, BAP170 and SAYP, in Drosophila. These subunits were stably tethered to a transgene reporter carrying the hsp70 core promoter. The tethered subunits mediate transcription of the reporter in a pattern that is generated by enhancers close to the insertion site in multiple loci throughout the genome. Both tethered SAYP and BAP170 recruit the whole PBAP complex to the reporter promoter. However, we found that BAP170-dependent transcription is more resistant to the depletion of other PBAP subunits, suggesting that BAP170 may play a more critical role in establishing enhancer-dependent transcription.

CTCF loss has limited effects on global genome architecture in Drosophila despite critical regulatory functions

Vertebrate genomes are partitioned into contact domains defined by enhanced internal contact frequency and formed by two principal mechanisms: compartmentalization of transcriptionally active and inactive domains, and stalling of chromosomal loop-extruding cohesin by CTCF bound at domain boundaries. While Drosophila has widespread contact domains and CTCF, it is currently unclear whether CTCF-dependent domains exist in flies. We genetically ablate CTCF in Drosophila and examine impacts on genome folding and transcriptional regulation in the central nervous system. We find that CTCF is required to form a small fraction of all domain boundaries, while critically controlling expression patterns of certain genes and supporting nervous system function. We also find that CTCF recruits the pervasive boundary-associated factor Cp190 to CTCF-occupied boundaries and co-regulates a subset of genes near boundaries together with Cp190. These results highlight a profound difference in CTCF-requirement for genome folding in flies and vertebrates, in which a large fraction of boundaries are CTCF-dependent and suggest that CTCF has played mutable roles in genome architecture and direct gene expression control during metazoan evolution.

Independence of chromatin conformation and gene regulation during Drosophila dorsoventral patterning

The relationship between chromatin organization and gene regulation remains unclear. While disruption of chromatin domains and domain boundaries can lead to misexpression of developmental genes, acute depletion of regulators of genome organization has a relatively small effect on gene expression. It is therefore uncertain whether gene expression and chromatin state drive chromatin organization or whether changes in chromatin organization facilitate cell-type-specific activation of gene expression. Here, using the dorsoventral patterning of the Drosophila melanogaster embryo as a model system, we provide evidence for the independence of chromatin organization and dorsoventral gene expression. We define tissue-specific enhancers and link them to expression patterns using single-cell RNA-seq. Surprisingly, despite tissue-specific chromatin states and gene expression, chromatin organization is largely maintained across tissues. Our results indicate that tissue-specific chromatin conformation is not necessary for tissue-specific gene expression but rather acts as a scaffold facilitating gene expression when enhancers become active.

A Role for the Mre11-Rad50-Xrs2 Complex in Gene Expression and Chromosome Organization

Mre11-Rad50-Xrs2 (MRX) is a highly conserved complex with key roles in various aspects of DNA repair. Here, we report a new function for MRX in limiting transcription in budding yeast. We show that MRX interacts physically and colocalizes on chromatin with the transcriptional co-regulator Mediator. MRX restricts transcription of coding and noncoding DNA by a mechanism that does not require the nuclease activity of Mre11. MRX is required to tether transcriptionally active loci to the nuclear pore complex (NPC), and it also promotes large-scale gene-NPC interactions. Moreover, MRX-mediated chromatin anchoring to the NPC contributes to chromosome folding and helps to control gene expression. Together, these findings indicate that MRX has a role in transcription and chromosome organization that is distinct from its known function in DNA repair.

bsAS, an antisense long non-coding RNA, essential for correct wing development through regulation of blistered/DSRF isoform usage

Natural Antisense Transcripts (NATs) are long non-coding RNAs (lncRNAs) that overlap coding genes in the opposite strand. NATs roles have been related to gene regulation through different mechanisms, including post-transcriptional RNA processing. With the aim to identify NATs with potential regulatory function during fly development, we generated RNA-Seq data in Drosophila developing tissues and found bsAS, one of the most highly expressed lncRNAs in the fly wing. bsAS is antisense to bs/DSRF, a gene involved in wing development and neural processes. bsAS plays a crucial role in the tissue specific regulation of the expression of the bs/DSRF isoforms. This regulation is essential for the correct determination of cell fate during Drosophila development, as bsAS knockouts show highly aberrant phenotypes. Regulation of bs isoform usage by bsAS is mediated by specific physical interactions between the promoters of these two genes, which suggests a regulatory mechanism involving the collision of RNA polymerases transcribing in opposite directions. Evolutionary analysis suggests that bsAS NAT emerged simultaneously to the long-short isoform structure of bs, preceding the emergence of wings in insects.

Molecular basis of CTCF binding polarity in genome folding

Current models propose that boundaries of mammalian topologically associating domains (TADs) arise from the ability of the CTCF protein to stop extrusion of chromatin loops by cohesin. While the orientation of CTCF motifs determines which pairs of CTCF sites preferentially stabilize loops, the molecular basis of this polarity remains unclear. By combining ChIP-seq and single molecule live imaging we report that CTCF positions cohesin, but does not control its overall binding dynamics on chromatin. Using an inducible complementation system, we find that CTCF mutants lacking the N-terminus cannot insulate TADs properly. Cohesin remains at CTCF sites in this mutant, albeit with reduced enrichment. Given the orientation of CTCF motifs presents the N-terminus towards cohesin as it translocates from the interior of TADs, these observations explain how the orientation of CTCF binding sites translates into genome folding patterns.

Chromatin Remodelers in the 3D Nuclear Compartment

Chromatin remodeling complexes (CRCs) use ATP hydrolysis to maintain correct expression profiles, chromatin stability, and inherited epigenetic states. More than 20 CRCs have been described to date, which encompass four large families defined by their ATPase subunits. These complexes and their subunits are conserved from yeast to humans through evolution. Their activities depend on their catalytic subunits which through ATP hydrolysis provide the energy necessary to fulfill cellular functions such as gene transcription, DNA repair, and transposon silencing. These activities take place at the first levels of chromatin compaction, and CRCs have been recognized as essential elements of chromatin dynamics. Recent studies have demonstrated an important role for these complexes in the maintenance of higher order chromatin structure. In this review, we present an overview of the organization of the genome within the cell nucleus, the different levels of chromatin compaction, and importance of the architectural proteins, and discuss the role of CRCs and how their functions contribute to the dynamics of the 3D genome organization.

Ecdysone-Induced 3D Chromatin Reorganization Involves Active Enhancers Bound by Pipsqueak and Polycomb

Evidence suggests that Polycomb (Pc) is present at chromatin loop anchors in Drosophila. Pc is recruited to DNA through interactions with the GAGA binding factors GAF and Pipsqueak (Psq). Using HiChIP in Drosophila cells, we find that the psq gene, which has diverse roles in development and tumorigenesis, encodes distinct isoforms with unanticipated roles in genome 3D architecture. The BR-C, ttk, and bab domain (BTB)-containing Psq isoform (Psq^L) colocalizes genome-wide with known architectural proteins. Conversely, Psq lacking the BTB domain (Psq^S) is consistently found at Pc loop anchors and at active enhancers, including those that respond to the hormone ecdysone. After stimulation by this hormone, chromatin 3D organization is altered to connect promoters and ecdysone-responsive enhancers bound by Psq^S. Our findings link Psq variants lacking the BTB domain to Pc-bound active enhancers, thus shedding light into their molecular function in chromatin changes underlying the response to hormone stimulus.

Direct DNA crosslinking with CAP-C uncovers transcription-dependent chromatin organization at high resolution

Determining the spatial organization of chromatin in cells mainly relies on crosslinking-based chromosome conformation capture techniques, but resolution and signal-to-noise ratio of these approaches is limited by interference from DNA-bound proteins. Here we introduce chemical-crosslinking assisted proximity capture (CAP-C), a method that uses multifunctional chemical crosslinkers with defined sizes to capture chromatin contacts. CAP-C generates chromatin contact maps at subkilobase (sub-kb) resolution with low background noise. We applied CAP-C to formaldehyde prefixed mouse embryonic stem cells (mESCs) and investigated loop domains (median size of 200 kb) and nonloop domains (median size of 9 kb). Transcription inhibition caused a greater loss of contacts in nonloop domains than loop domains. We uncovered conserved, transcription-state-dependent chromatin compartmentalization at high resolution that is shared from Drosophila to human, and a transcription-initiation-dependent nuclear subcompartment that brings multiple nonloop domains in close proximity. We also showed that CAP-C could be used to detect native chromatin conformation without formaldehyde prefixing.

Evolution of Regulated Transcription

The genomes of all organisms abound with various cis-regulatory elements, which control gene activity. Transcriptional enhancers are a key group of such elements in eukaryotes and are DNA regions that form physical contacts with gene promoters and precisely orchestrate gene expression programs. Here, we follow gradual evolution of this regulatory system and discuss its features in different organisms. In eubacteria, an enhancer-like element is often a single regulatory element, is usually proximal to the core promoter, and is occupied by one or a few activators. Activation of gene expression in archaea is accompanied by the recruitment of an activator to several enhancer-like sites in the upstream promoter region. In eukaryotes, activation of expression is accompanied by the recruitment of activators to multiple enhancers, which may be distant from the core promoter, and the activators act through coactivators. The role of the general DNA architecture in transcription control increases in evolution. As a whole, it can be seen that enhancers of multicellular eukaryotes evolved from the corresponding prototypic enhancer-like regulatory elements with the gradually increasing genome size of organisms.