A Functional Insulator Screen Identifies NURF and dREAM Components to Be Required for Enhancer-Blocking

Chromatin insulators of higher eukaryotes functionally divide the genome into active and inactive domains. Furthermore, insulators regulate enhancer/promoter communication, which is evident from the Drosophila bithorax locus in which a multitude of regulatory elements control segment specific gene activity. Centrosomal protein 190 (CP190) is targeted to insulators by CTCF or other insulator DNA-binding factors. Chromatin analyses revealed that insulators are characterized by open and nucleosome depleted regions. Here, we wanted to identify chromatin modification and remodelling factors required for an enhancer blocking function. We used the well-studied Fab-8 insulator of the bithorax locus to apply a genome-wide RNAi screen for factors that contribute to the enhancer blocking function of CTCF and CP190. Among 78 genes required for optimal Fab-8 mediated enhancer blocking, all four components of the NURF complex as well as several subunits of the dREAM complex were most evident. Mass spectrometric analyses of CTCF or CP190 bound proteins as well as immune precipitation confirmed NURF and dREAM binding. Both co-localise with most CP190 binding sites in the genome and chromatin immune precipitation showed that CP190 recruits NURF and dREAM. Nucleosome occupancy and histone H3 binding analyses revealed that CP190 mediated NURF binding results in nucleosomal depletion at CP190 binding sites. Thus, we conclude that CP190 binding to CTCF or to other DNA binding insulator factors mediates recruitment of NURF and dREAM. Furthermore, the enhancer blocking function of insulators is associated with nucleosomal depletion and requires NURF and dREAM.     

Adaptive Evolution and the Birth of CTCF Binding Sites in the Drosophila Genome

Changes in the physical interaction between cis-regulatory DNA sequences and proteins drive the evolution of gene expression. However, it has proven difficult to accurately quantify evolutionary rates of such binding change or to estimate the relative effects of selection and drift in shaping the binding evolution. Here we examine the genome-wide binding of CTCF in four species of Drosophila separated by between ∼2.5 and 25 million years. CTCF is a highly conserved protein known to be associated with insulator sequences in the genomes of human and Drosophila. Although the binding preference for CTCF is highly conserved, we find that CTCF binding itself is highly evolutionarily dynamic and has adaptively evolved. Between species, binding divergence increased linearly with evolutionary distance, and CTCF binding profiles are diverging rapidly at the rate of 2.22% per million years (Myr). At least 89 new CTCF binding sites have originated in the Drosophila melanogaster genome since the most recent common ancestor with Drosophila simulans. Comparing these data to genome sequence data from 37 different strains of Drosophila melanogaster, we detected signatures of selection in both newly gained and evolutionarily conserved binding sites. Newly evolved CTCF binding sites show a significantly stronger signature for positive selection than older sites. Comparative gene expression profiling revealed that expression divergence of genes adjacent to CTCF binding site is significantly associated with the gain and loss of CTCF binding. Further, the birth of new genes is associated with the birth of new CTCF binding sites. Our data indicate that binding of Drosophila CTCF protein has evolved under natural selection, and CTCF binding evolution has shaped both the evolution of gene expression and genome evolution during the birth of new genes.     

NURF301 contributes to gypsy chromatin insulator-mediated nuclear organization

Chromatin insulators are DNA-protein complexes that can prevent the spread of repressive chromatin and block communication between enhancers and promoters to regulate gene expression. In Drosophila, the gypsy chromatin insulator complex consists of three core proteins: CP190, Su(Hw), and Mod(mdg4)67.2. These factors concentrate at nuclear foci termed insulator bodies, and changes in insulator body localization have been observed in mutants defective for insulator function. Here, we identified NURF301/E(bx), a nucleosome remodeling factor, as a novel regulator of gypsy insulator body localization through a high-throughput RNAi imaging screen. NURF301 promotes gypsy-dependent insulator barrier activity and physically interacts with gypsy insulator proteins. Using ChIP-seq, we found that NURF301 co-localizes with insulator proteins genome-wide, and NURF301 promotes chromatin association of Su(Hw) and CP190 at gypsy insulator binding sites. These effects correlate with NURF301-dependent nucleosome repositioning. At the same time, CP190 and Su(Hw) both facilitate recruitment of NURF301 to chromatin. Finally, Oligopaint FISH combined with immunofluorescence revealed that NURF301 promotes 3D contact between insulator bodies and gypsy insulator DNA binding sites, and NURF301 is required for proper nuclear positioning of gypsy binding sites. Our data provide new insights into how a nucleosome remodeling factor and insulator proteins cooperatively contribute to nuclear organization.     

Functional interaction between the Fab-7 and Fab-8 boundaries and the upstream promoter region in the Drosophila Abd-B gene

Boundary elements have been found in the regulatory region of the Drosophila melanogaster Abdominal-B (Abd-B) gene, which is subdivided into a series of iab domains. The best-studied Fab-7 and Fab-8 boundaries flank the iab-7 enhancer and isolate it from the four promoters regulating Abd-B expression. Recently binding sites for the Drosophila homolog of the vertebrate insulator protein CTCF (dCTCF) were identified in the Fab-8 boundary and upstream of Abd-B promoter A, with no binding of CTCF to the Fab-7 boundary being detected either in vivo or in vitro. Taking into account the inability of the yeast GAL4 activator to stimulate the white promoter when its binding sites are separated by a 5-kb yellow gene, we have tested the functional interactions between the Fab-7 and Fab-8 boundaries and between these boundaries and the upstream promoter A region containing a dCTCF binding site. It has been found that dCTCF binding sites are essential for pairing between two Fab-8 insulators. However, a strong functional interaction between the Fab-7 and Fab-8 boundaries suggests that additional, as yet unidentified proteins are involved in long-distance interactions between them. We have also shown that Fab-7 and Fab-8 boundaries effectively interact with the upstream region of the Abd-B promoter.     

Genome-Wide Identification of Chromatin Transitional Regions Reveals Diverse Mechanisms Defining the Boundary of Facultative Heterochromatin

Due to the self-propagating nature of the heterochromatic modification H3K27me3, chromatin barrier activities are required to demarcate the boundary and prevent it from encroaching into euchromatic regions. Studies in Drosophila and vertebrate systems have revealed several important chromatin barrier elements and their respective binding factors. However, epigenomic data indicate that the binding of these factors are not exclusive to chromatin boundaries. To gain a comprehensive understanding of facultative heterochromatin boundaries, we developed a two-tiered method to identify the Chromatin Transitional Region (CTR), i.e. the nucleosomal region that shows the greatest transition rate of the H3K27me3 modification as revealed by ChIP-Seq. This approach was applied to identify CTRs in Drosophila S2 cells and human HeLa cells. Although many insulator proteins have been characterized in Drosophila, less than half of the CTRs in S2 cells are associated with known insulator proteins, indicating unknown mechanisms remain to be characterized. Our analysis also revealed that the peak binding of insulator proteins are usually 1–2 nucleosomes away from the CTR. Comparison of CTR-associated insulator protein binding sites vs. those in heterochromatic region revealed that boundary-associated binding sites are distinctively flanked by nucleosome destabilizing sequences, which correlates with significant decreased nucleosome density and increased binding intensities of co-factors. Interestingly, several subgroups of boundaries have enhanced H3.3 incorporation but reduced nucleosome turnover rate. Our genome-wide study reveals that diverse mechanisms are employed to define the boundaries of facultative heterochromatin. In both Drosophila and mammalian systems, only a small fraction of insulator protein binding sites co-localize with H3K27me3 boundaries. However, boundary-associated insulator binding sites are distinctively flanked by nucleosome destabilizing sequences, which correlates with significantly decreased nucleosome density and increased binding of co-factors.     

A possible role of Drosophila CTCF in mitotic bookmarking and maintaining chromatin domains during the cell cycle

Background

The CCCTC-binding factor (CTCF) is a highly conserved insulator protein that plays various roles in many cellular processes. CTCF is one of the main architecture proteins in higher eukaryotes, and in combination with other architecture proteins and regulators, also shapes the three-dimensional organization of a genome. Experiments show CTCF partially remains associated with chromatin during mitosis. However, the role of CTCF in the maintenance and propagation of genome architectures throughout the cell cycle remains elusive.

Results

We performed a comprehensive bioinformatics analysis on public datasets of Drosophila CTCF (dCTCF). We characterized dCTCF-binding sites according to their occupancy status during the cell cycle, and identified three classes: interphase-mitosis-common (IM), interphase-only (IO) and mitosis-only (MO) sites. Integrated function analysis showed dCTCF-binding sites of different classes might be involved in different biological processes, and IM sites were more conserved and more intensely bound. dCTCF-binding sites of the same class preferentially localized closer to each other, and were highly enriched at chromatin syntenic and topologically associating domains boundaries.

Conclusions

Our results revealed different functions of dCTCF during the cell cycle and suggested that dCTCF might contribute to the establishment of the three-dimensional architecture of the Drosophila genome by maintaining local chromatin compartments throughout the whole cell cycle.

Electronic supplementary material

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

A Comprehensive Map of Insulator Elements for the Drosophila Genome

Insulators are DNA sequences that control the interactions among genomic regulatory elements and act as chromatin boundaries. A thorough understanding of their location and function is necessary to address the complexities of metazoan gene regulation. We studied by ChIP–chip the genome-wide binding sites of 6 insulator-associated proteins—dCTCF, CP190, BEAF-32, Su(Hw), Mod(mdg4), and GAF—to obtain the first comprehensive map of insulator elements in Drosophila embryos. We identify over 14,000 putative insulators, including all classically defined insulators. We find two major classes of insulators defined by dCTCF/CP190/BEAF-32 and Su(Hw), respectively. Distributional analyses of insulators revealed that particular sub-classes of insulator elements are excluded between cis-regulatory elements and their target promoters; divide differentially expressed, alternative, and divergent promoters; act as chromatin boundaries; are associated with chromosomal breakpoints among species; and are embedded within active chromatin domains. Together, these results provide a map demarcating the boundaries of gene regulatory units and a framework for understanding insulator function during the development and evolution of Drosophila.     

Coordinated control of dCTCF and gypsy chromatin insulators in Drosophila

CTCF plays a central role in vertebrate insulators and forms part of the Fab-8 insulator in Drosophila. dCTCF is present at hundreds of sites in the Drosophila genome, where it is located at the boundaries between bands and interbands in polytene chromosomes. dCTCF colocalizes with CP190, which is required for proper binding of dCTCF to chromatin, but not with the other gypsy insulator proteins Su(Hw) or Mod(mdg4)2.2. Mutations in the CP190 gene affect Fab-8 insulator activity, suggesting that CP190 is an essential component of both gypsy and dCTCF insulators. dCTCF is present at specific nuclear locations, forming large insulator bodies that overlap with those formed by Su(Hw), Mod(mdg4)2.2, and CP190. The results suggest that Su(Hw) and dCTCF may be the DNA-binding components of two different subsets of insulators that share CP190 and cooperate in the formation of insulator bodies to regulate the organization of the chromatin fiber in the nucleus.     

Chromatin insulator and the promoter targeting sequence modulate the timing of long-range enhancer-promoter interactions in the Drosophila embryo

The homeotic genes are essential to the patterning of the anterior-posterior axis along the developing Drosophila embryo. The expression timing and levels of these genes are crucial for the correct specification of segmental identity. The Abdominal-B (Abd-B) gene is first detected in the most posterior abdominal segments at high levels and gradually appears in progressively anterior abdominal segments in lower amounts. Regulatory mutations affecting this expression pattern produce homeotic transformations in the abdomen. The promoter targeting sequences (PTS) from Abd-B locus overcome the enhancer blocking effect of insulators and facilitate long-range enhancer-promoter interactions in transgenic flies (1, 2). In this study, we found that transgene activation by the IAB5 enhancer can be delayed by inserting a 9.5 kb 3′ Abd-B regulatory region containing the Frontabdominal-8 (Fab-8) insulator and the PTS element. We found that the delay is caused by the PTS and an insulator, and it is not specific to the enhancer or the promoter tested. Based on these findings, we hypothesize that the delay of remote enhancers is responsible for the Abd-B expression pattern, which is at least in part due to the regulatory activities of the PTS elements and chromatin boundaries.     

Functional Requirements for Fab-7 Boundary Activity in the Bithorax Complex

Chromatin boundaries are architectural elements that determine the three-dimensional folding of the chromatin fiber and organize the chromosome into independent units of genetic activity. The Fab-7 boundary from the Drosophila bithorax complex (BX-C) is required for the parasegment-specific expression of the Abd-B gene. We have used a replacement strategy to identify sequences that are necessary and sufficient for Fab-7 boundary function in the BX-C. Fab-7 boundary activity is known to depend on factors that are stage specific, and we describe a novel ∼700-kDa complex, the late boundary complex (LBC), that binds to Fab-7 sequences that have insulator functions in late embryos and adults. We show that the LBC is enriched in nuclear extracts from late, but not early, embryos and that it contains three insulator proteins, GAF, Mod(mdg4), and E(y)2. Its DNA binding properties are unusual in that it requires a minimal sequence of >65 bp; however, other than a GAGA motif, the three Fab-7 LBC recognition elements display few sequence similarities. Finally, we show that mutations which abrogate LBC binding in vitro inactivate the Fab-7 boundary in the BX-C.     

Multiple Promoter Targeting Sequences exist in Abdominal-B to regulate long-range gene activation

In complex genomes, insulators set up chromatin domain boundaries and protect promoters from inappropriate activation by enhancers from neighboring genes. The Drosophila Abdominal-B locus uses insulator elements to organize its large regulatory region into several body segment-specific chromatin domains. This organization leads to a problem in enhancer-promoter communication, that is, how do distal enhancers activate the Abd-B promoter when there are several insulators in between? This issue is partially resolved by the Promoter Targeting Sequence, which can overcome the enhancer blocking effect of an insulator. In this study, we describe a new Promoter Targeting Sequence, PTS-6, from the Abd-B 3' regulatory region. PTS-6, comprised of approximately 200 bp, was found to bypass both homologous Abdominal-B insulators, such as Fab-7 and Fab-8, and a heterologous insulator, suHw. Most importantly, it also overcomes a combination of two insulators such as Fab-7/Fab-8. Thus, PTS-6 could, in principle, target remote enhancers that are separated from the Abd-B promoter by multiple insulators. In addition, PTS-6 selectively targets the distal enhancer to only one transgenic promoter, and it strongly facilitates Abd-B enhancers. These results suggest that promoter targeting is necessary for long-range enhancer-promoter communication in Abd-B, and PTS elements could be a common occurrence in large, complex genetic loci.