Regulation of chromatin organization and inducible gene expression by a Drosophila insulator

SetDB1 and Su(var)3-9 are essential for late stages of larval development of Drosophila melanogaster

Methylation of H3K9 histone residue is a marker of gene silencing in eukaryotes. Three enzymes responsible for adding this modification - G9a, SetDB1/Egg, and Su(var)3-9 - are known in Drosophila. To understand how simultaneous mutations of SetDB1 and Su(var)3-9 may affect the fly development, appropriate combinations were obtained. Double mutants egg; Su(var)3-9 displayed pronounced embryonic lethality, slower larval growth and died before or during metamorphosis. Analysis of transcription in larval salivary glands and wing imaginal disks indicated that the effect of double mutation is tissue-specific. In salivary gland chromosomes, affected genes display low H3K9me2 enrichment and are rarely bound by SetDB1 or Su(var)3-9. We suppose that each of these enzymes directly or indirectly controls its own set of gene targets in different organs, and double mutation results in an imbalanced developmental program. This also indicates that SetDB1 and Su(var)3-9 may affect transcription via H3K9-independent mechanisms. Unexpectedly, in double and triple mutants, amount of di- and tri-methylated H3K9 is drastically reduced, but not completely absent. We hypothesize that this residual methylation implies the existence of additional H3K9-specific methyltransferase in Drosophila.

Coregulators Reside within Drosophila Ecdysone-Inducible Loci before and after Ecdysone Treatment

Ecdysone signaling in Drosophila remains a popular model for investigating the mechanisms of steroid action in eukaryotes. The ecdysone receptor EcR can effectively bind ecdysone-response elements with or without the presence of a hormone. For years, EcR enhancers were thought to respond to ecdysone via recruiting coactivator complexes, which replace corepressors and stimulate transcription. However, the exact mechanism of transcription activation by ecdysone remains unclear. Here, we present experimental data on 11 various coregulators at ecdysone-responsive loci of Drosophila S2 cells. We describe the regulatory elements where coregulators reside within these loci and assess changes in their binding levels following 20-hydroxyecdysone treatment. In the current study, we detected the presence of some coregulators at the TSSs (active and inactive) and boundaries marked with CP190 rather than enhancers of the ecdysone-responsive loci where EcR binds. We observed minor changes in the coregulators' binding level. Most were present at inducible loci before and after 20-hydroxyecdysone treatment. Our findings suggest that: (1) coregulators can activate a particular TSS operating from some distal region (which could be an enhancer, boundary regulatory region, or inactive TSS); (2) coregulators are not recruited after 20-hydroxyecdysone treatment to the responsive loci; rather, their functional activity changes (shown as an increase in H3K27 acetylation marks generated by CBP/p300/Nejire acetyltransferase). Taken together, our findings imply that the 20-hydroxyecdysone signal enhances the functional activity of coregulators rather than promoting their binding to regulatory regions during the ecdysone response.

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.

Essential role of Cp190 in physical and regulatory boundary formation

Boundaries in animal genomes delimit contact domains with enhanced internal contact frequencies and have debated functions in limiting regulatory cross-talk between domains and guiding enhancers to target promoters. Most mammalian boundaries form by stalling of chromosomal loop-extruding cohesin by CTCF, but most Drosophila boundaries form CTCF independently. However, how CTCF-independent boundaries form and function remains largely unexplored. Here, we assess genome folding and developmental gene expression in fly embryos lacking the ubiquitous boundary-associated factor Cp190. We find that sequence-specific DNA binding proteins such as CTCF and Su(Hw) directly interact with and recruit Cp190 to form most promoter-distal boundaries. Cp190 is essential for early development and prevents regulatory cross-talk between specific gene loci that pattern the embryo. Cp190 was, in contrast, dispensable for long-range enhancer-promoter communication at tested loci. Cp190 is thus currently the major player in fly boundary formation and function, revealing that diverse mechanisms evolved to partition genomes into independent regulatory domains.

Transcriptional coupling of distant regulatory genes in living embryos

The prevailing view of metazoan gene regulation is that individual genes are independently regulated by their own dedicated sets of transcriptional enhancers. Past studies have reported long-range gene-gene associations^1-3, but their functional importance in regulating transcription remains unclear. Here we used quantitative single-cell live imaging methods to provide a demonstration of co-dependent transcriptional dynamics of genes separated by large genomic distances in living Drosophila embryos. We find extensive physical and functional associations of distant paralogous genes, including co-regulation by shared enhancers and co-transcriptional initiation over distances of nearly 250 kilobases. Regulatory interconnectivity depends on promoter-proximal tethering elements, and perturbations in these elements uncouple transcription and alter the bursting dynamics of distant genes, suggesting a role of genome topology in the formation and stability of co-transcriptional hubs. Transcriptional coupling is detected throughout the fly genome and encompasses a broad spectrum of conserved developmental processes, suggesting a general strategy for long-range integration of gene activity.

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.

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.

Proximity-dependent biotin labelling reveals CP190 as an EcR/Usp molecular partner

Proximity-dependent biotin labelling revealed undescribed participants of the ecdysone response in Drosophila. Two labelling enzymes (BioID2 and APEX2) were fused to EcR or Usp to biotin label the surrounding proteins. The EcR/Usp heterodimer was found to collaborate with nuclear pore subunits, chromatin remodelers, and architectural proteins. Many proteins identified through proximity-dependent labelling with EcR/Usp were described previously as functional components of an ecdysone response, corroborating the potency of this labelling method. A link to ecdysone response was confirmed for some newly discovered regulators by immunoprecipitation of prepupal nuclear extract with anti-EcR antibodies and functional experiments in Drosophila S2 cells. A more in-depth study was conducted to clarify the association of EcR/Usp with one of the detected proteins, CP190, a well-described cofactor of Drosophila insulators. CP190 was found to co-immunoprecipitate with the EcR subunit of EcR/Usp in a 20E-independent manner. ChIP-Seq experiments revealed only partial overlapping between CP190 and EcR bound sites in the Drosophila genome and complete absence of CP190 binding at 20E-dependent enhancers. Analysis of Hi-C data demonstrated an existence of remote interactions between 20E-dependent enhancers and CP190 sites which suggests formation of a protein complex between EcR/Usp and CP190 through the space. Our results support the previous concept that CP190 has a role in stabilization of specific chromatin loops for proper activation of transcription of genes regulated by 20E hormone.

In situ dissection of domain boundaries affect genome topology and gene transcription in Drosophila

Chromosomes are organized into high-frequency chromatin interaction domains called topologically associating domains (TADs), which are separated from each other by domain boundaries. The molecular mechanisms responsible for TAD formation are not yet fully understood. In Drosophila, it has been proposed that transcription is fundamental for TAD organization while the participation of genetic sequences bound by architectural proteins (APs) remains controversial. Here, we investigate the contribution of domain boundaries to TAD organization and the regulation of gene expression at the Notch gene locus in Drosophila. We find that deletion of domain boundaries results in TAD fusion and long-range topological defects that are accompanied by loss of APs and RNA Pol II chromatin binding as well as defects in transcription. Together, our results provide compelling evidence of the contribution of discrete genetic sequences bound by APs and RNA Pol II in the partition of the genome into TADs and in the regulation of gene expression in Drosophila.

In Drosophila, transcription is thought to be required for TAD formation, while the role of architectural proteins remains controversial. Here, the authors find that deletion of domain boundaries at the fly Notch locus results in TAD fusion and long-range topological defects, loss of architectural protein and RNA Pol II chromatin binding, and transcription defects.

BEN-solo factors partition active chromatin to ensure proper gene activation in Drosophila

The Drosophila genome encodes three BEN-solo proteins including Insensitive (Insv), Elba1 and Elba2 that possess activities in transcriptional repression and chromatin insulation. A fourth protein-Elba3-bridges Elba1 and Elba2 to form an ELBA complex. Here, we report comprehensive investigation of these proteins in Drosophila embryos. We assess common and distinct binding sites for Insv and ELBA and their genetic interdependencies. While Elba1 and Elba2 binding generally requires the ELBA complex, Elba3 can associate with chromatin independently of Elba1 and Elba2. We further demonstrate that ELBA collaborates with other insulators to regulate developmental patterning. Finally, we find that adjacent gene pairs separated by an ELBA bound sequence become less differentially expressed in ELBA mutants. Transgenic reporters confirm the insulating activity of ELBA- and Insv-bound sites. These findings define ELBA and Insv as general insulator proteins in Drosophila and demonstrate the functional importance of insulators to partition transcription units.

Characterization of Button Loci that Promote Homologous Chromosome Pairing and Cell-Type-Specific Interchromosomal Gene Regulation

Homologous chromosomes colocalize to regulate gene expression in processes including genomic imprinting, X-inactivation, and transvection. In Drosophila, homologous chromosomes pair throughout development, promoting transvection. The "button" model of pairing proposes that specific regions along chromosomes pair with high affinity. Here, we identify buttons interspersed across the fly genome that pair with their homologous sequences, even when relocated to multiple positions in the genome. A majority of transgenes that span a full topologically associating domain (TAD) function as buttons, but not all buttons contain TADs. Additionally, buttons are enriched for insulator protein clusters. Fragments of buttons do not pair, suggesting that combinations of elements within a button are required for pairing. Pairing is necessary but not sufficient for transvection. Additionally, pairing and transvection are stronger in some cell types than in others, suggesting that pairing strength regulates transvection efficiency between cell types. Thus, buttons pair homologous chromosomes to facilitate cell-type-specific interchromosomal gene regulation.

The Role of Insulation in Patterning Gene Expression

Development is orchestrated by regulatory elements that turn genes ON or OFF in precise spatial and temporal patterns. Many safety mechanisms prevent inappropriate action of a regulatory element on the wrong gene promoter. In flies and mammals, dedicated DNA elements (insulators) recruit protein factors (insulator binding proteins, or IBPs) to shield promoters from regulatory elements. In mammals, a single IBP called CCCTC-binding factor (CTCF) is known, whereas genetic and biochemical analyses in Drosophila have identified a larger repertoire of IBPs. How insulators function at the molecular level is not fully understood, but it is currently thought that they fold chromosomes into conformations that affect regulatory element-promoter communication. Here, we review the discovery of insulators and describe their properties. We discuss recent genetic studies in flies and mice to address the question: Is gene insulation important for animal development? Comparing and contrasting observations in these two species reveal that they have different requirements for insulation, but that insulation is a conserved and critical gene regulation strategy.

TAD-free analysis of architectural proteins and insulators

The three-dimensional (3D) organization of the genome is intimately related to numerous key biological functions including gene expression and DNA replication regulations. The mechanisms by which molecular drivers functionally organize the 3D genome, such as topologically associating domains (TADs), remain to be explored. Current approaches consist in assessing the enrichments or influences of proteins at TAD borders. Here, we propose a TAD-free model to directly estimate the blocking effects of architectural proteins, insulators and DNA motifs on long-range contacts, making the model intuitive and biologically meaningful. In addition, the model allows analyzing the whole Hi-C information content (2D information) instead of only focusing on TAD borders (1D information). The model outperforms multiple logistic regression at TAD borders in terms of parameter estimation accuracy and is validated by enhancer-blocking assays. In Drosophila, the results support the insulating role of simple sequence repeats and suggest that the blocking effects depend on the number of repeats. Motif analysis uncovered the roles of the transcriptional factors pannier and tramtrack in blocking long-range contacts. In human, the results suggest that the blocking effects of the well-known architectural proteins CTCF, cohesin and ZNF143 depend on the distance between loci, where each protein may participate at different scales of the 3D chromatin organization.

The zinc-finger protein CLAMP promotes gypsy chromatin insulator function in Drosophila

Chromatin insulators are DNA-protein complexes that establish independent higher-order DNA domains to influence transcription. Insulators are functionally defined by two properties: they can block communication between an enhancer and a promoter, and also act as a barrier between heterochromatin and euchromatin. In Drosophila, the gypsy insulator complex contains three core components; Su(Hw), CP190 and Mod(mdg4)67.2. Here, we identify a novel role for Chromatin-linked adaptor for MSL proteins (CLAMP) in promoting gypsy chromatin insulator function. When clamp is knocked down, gypsy-dependent enhancer-blocking and barrier activities are strongly reduced. CLAMP associates physically with the core gypsy insulator complex, and ChIP-seq analysis reveals extensive overlap, particularly with promoter-bound CP190 on chromatin. Depletion of CLAMP disrupts CP190 binding at a minority of shared sites, whereas depletion of CP190 results in extensive loss of CLAMP chromatin association. Finally, reduction of CLAMP disrupts CP190 localization within the nucleus. Our results support a positive functional relationship between CLAMP and CP190 to promote gypsy chromatin insulator activity.

Tethering of CHROMATOR and dCTCF proteins results in decompaction of condensed bands in the Drosophila melanogaster polytene chromosomes but does not affect their transcription and replication timing

Instulator proteins are central to domain organization and gene regulation in the genome. We used ectopic tethering of CHROMATOR (CHRIZ/CHRO) and dCTCF to pre-defined regions of the genome to dissect the influence of these proteins on local chromatin organization, to analyze their interaction with other key chromatin proteins and to evaluate the effects on transcription and replication. Specifically, using UAS-GAL4DBD system, CHRO and dCTCF were artificially recruited into highly compacted polytene chromosome bands that share the features of silent chromatin type known as intercalary heterochromatin (IH). This led to local chromatin decondensation, formation of novel DHSes and recruitment of several "open chromatin" proteins. CHRO tethering resulted in the recruitment of CP190 and Z4 (PZG), whereas dCTCF tethering attracted CHRO, CP190, and Z4. Importantly, formation of a local stretch of open chromatin did not result in the reactivation of silent marker genes yellow and mini-white immediately adjacent to the targeting region (UAS), nor did RNA polII become recruited into this chromatin. The decompacted region retained late replicated, similarly to the wild-type untargeted region.

High-resolution TADs reveal DNA sequences underlying genome organization in flies

Despite an abundance of new studies about topologically associating domains (TADs), the role of genetic information in TAD formation is still not fully understood. Here we use our software, HiCExplorer (hicexplorer.readthedocs.io) to annotate >2800 high-resolution (570 bp) TAD boundaries in Drosophila melanogaster. We identify eight DNA motifs enriched at boundaries, including a motif bound by the M1BP protein, and two new boundary motifs. In contrast to mammals, the CTCF motif is only enriched on a small fraction of boundaries flanking inactive chromatin while most active boundaries contain the motifs bound by the M1BP or Beaf-32 proteins. We demonstrate that boundaries can be accurately predicted using only the motif sequences at open chromatin sites. We propose that DNA sequence guides the genome architecture by allocation of boundary proteins in the genome. Finally, we present an interactive online database to access and explore the spatial organization of fly, mouse and human genomes, available at http://chorogenome.ie-freiburg.mpg.de .

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

Eukaryotic gene expression is regulated by enhancer–promoter interactions but the molecular mechanisms that govern specificity have remained elusive. Genome-wide studies utilizing STARR-seq identified two enhancer classes in Drosophila that interact with different core promoters: housekeeping enhancers (hkCP) and developmental enhancers (dCP). We hypothesized that the two enhancer classes are occupied by distinct architectural proteins, affecting their enhancer–promoter contacts. By evaluating ChIP-seq occupancy of architectural proteins, typical enhancer-associated proteins, and histone modifications, we determine that both enhancer classes are enriched for RNA Polymerase II, CBP, and architectural proteins but there are also distinctions. hkCP enhancers contain H3K4me3 and exclusively bind Cap-H2, Chromator, DREF and Z4, whereas dCP enhancers contain H3K4me1 and are more enriched for Rad21 and Fs(1)h-L. Additionally, we map the interactions of each enhancer class utilizing a Hi-C dataset with <1 kb resolution. Results suggest that hkCP enhancers are more likely to form multi-TSS interaction networks and be associated with topologically associating domain (TAD) borders, while dCP enhancers are more often bound to one or two TSSs and are enriched at chromatin loop anchors. The data support a model suggesting that the unique architectural protein occupancy within enhancers is one contributor to enhancer–promoter interaction specificity.

Chromatin binding of Gcn5 in Drosophila is largely mediated by CP190

Centrosomal 190 kDa protein (CP190) is a promoter binding factor, mediates long-range interactions in the context of enhancer-promoter contacts and in chromosomal domain formation. All Drosophila insulator proteins bind CP190 suggesting a crucial role in insulator function. CP190 has major effects on chromatin, such as depletion of nucleosomes, high nucleosomal turnover and prevention of heterochromatin expansion. Here, we searched for enzymes, which might be involved in CP190 mediated chromatin changes. Eighty percent of the genomic binding sites of the histone acetyltransferase Gcn5 are colocalizing with CP190 binding. Depletion of CP190 reduces Gcn5 binding to chromatin. Binding dependency was further supported by Gcn5 mediated co-precipitation of CP190. Gcn5 is known to activate transcription by histone acetylation. We used the dCas9 system to target CP190 or Gcn5 to a Polycomb repressed and H3K27me3 marked gene locus. Both, CP190 as well as Gcn5, activate this locus, thus supporting the model that CP190 recruits Gcn5 and thereby activates chromatin.

Polycomb-mediated chromatin loops revealed by a subkilobase-resolution chromatin interaction map

The locations of chromatin loops in Drosophila were determined by Hi-C (chemical cross-linking, restriction digestion, ligation, and high-throughput DNA sequencing). Whereas most loop boundaries or "anchors" are associated with CTCF protein in mammals, loop anchors in Drosophila were found most often in association with the polycomb group (PcG) protein Polycomb (Pc), a subunit of polycomb repressive complex 1 (PRC1). Loops were frequently located within domains of PcG-repressed chromatin. Promoters located at PRC1 loop anchors regulate some of the most important developmental genes and are less likely to be expressed than those not at PRC1 loop anchors. Although DNA looping has most commonly been associated with enhancer-promoter communication, our results indicate that loops are also associated with gene repression.

Tracking Multiple Genomic Elements Using Correlative CRISPR Imaging and Sequential DNA FISH

Live imaging of genome has offered important insights into the dynamics of the genome organization and gene expression. The demand to image simultaneously multiple genomic loci has prompted a flurry of exciting advances in multicolor CRISPR imaging, although color-based multiplexing is limited by the need for spectrally distinct fluorophores. Here we introduce an approach to achieve highly multiplexed live recording via correlative CRISPR imaging and sequential DNA fluorescence in situ hybridization (FISH). This approach first performs one-color live imaging of multiple genomic loci and then uses sequential rounds of DNA FISH to determine the loci identity. We have optimized the FISH protocol so that each round is complete in 1 min, demonstrating the identification of seven genomic elements and the capability to sustain reversible staining and washing for up to 20 rounds. We have also developed a correlation-based algorithm to faithfully register live and FISH images. Our approach keeps the rest of the color palette open to image other cellular phenomena of interest, as demonstrated by our simultaneous live imaging of genomic loci together with a cell cycle reporter. Furthermore, the algorithm to register faithfully between live and fixed imaging is directly transferrable to other systems such as multiplex RNA imaging with RNA-FISH and multiplex protein imaging with antibody-staining.

Chromatin Architecture Emerges during Zygotic Genome Activation Independent of Transcription

Chromatin architecture is fundamental in regulating gene expression. To investigate when spatial genome organization is first established during development, we examined chromatin conformation during Drosophila embryogenesis and observed the emergence of chromatin architecture within a tight time window that coincides with the onset of transcription activation in the zygote. Prior to zygotic genome activation, the genome is mostly unstructured. Early expressed genes serve as nucleation sites for topologically associating domain (TAD) boundaries. Activation of gene expression coincides with the establishment of TADs throughout the genome and co-localization of housekeeping gene clusters, which remain stable in subsequent stages of development. However, the appearance of TAD boundaries is independent of transcription and requires the transcription factor Zelda for locus-specific TAD boundary insulation. These results offer insight into when spatial organization of the genome emerges and identify a key factor that helps trigger this architecture.

Uncovering direct and indirect molecular determinants of chromatin loops using a computational integrative approach

Chromosomal organization in 3D plays a central role in regulating cell-type specific transcriptional and DNA replication timing programs. Yet it remains unclear to what extent the resulting long-range contacts depend on specific molecular drivers. Here we propose a model that comprehensively assesses the influence on contacts of DNA-binding proteins, cis-regulatory elements and DNA consensus motifs. Using real data, we validate a large number of predictions for long-range contacts involving known architectural proteins and DNA motifs. Our model outperforms existing approaches including enrichment test, random forests and correlation, and it uncovers numerous novel long-range contacts in Drosophila and human. The model uncovers the orientation-dependent specificity for long-range contacts between CTCF motifs in Drosophila, highlighting its conserved property in 3D organization of metazoan genomes. Our model further unravels long-range contacts depending on co-factors recruited to DNA indirectly, as illustrated by the influence of cohesin in stabilizing long-range contacts between CTCF sites. It also reveals asymmetric contacts such as enhancer-promoter contacts that highlight opposite influences of the transcription factors EBF1, EGR1 or MEF2C depending on RNA Polymerase II pausing.

Metazoan Nuclear Pores Provide a Scaffold for Poised Genes and Mediate Induced Enhancer-Promoter Contacts

Nuclear pore complex components (Nups) have been implicated in transcriptional regulation, yet what regulatory steps are controlled by metazoan Nups remains unclear. We identified the presence of multiple Nups at promoters, enhancers, and insulators in the Drosophila genome. In line with this binding, we uncovered a functional role for Nup98 in mediating enhancer-promoter looping at ecdysone-inducible genes. These genes were found to be stably associated with nuclear pores before and after activation. Although changing levels of Nup98 disrupted enhancer-promoter contacts, it did not affect ongoing transcription but instead compromised subsequent transcriptional activation or transcriptional memory. In support of the enhancer-looping role, we found Nup98 to gain and retain physical interactions with architectural proteins upon stimulation with ecdysone. Together, our data identify Nups as a class of architectural proteins for enhancers and supports a model in which animal genomes use the nuclear pore as an organizing scaffold for inducible poised genes.

Regions of very low H3K27me3 partition the Drosophila genome into topological domains

It is now well established that eukaryote genomes have a common architectural organization into topologically associated domains (TADs) and evidence is accumulating that this organization plays an important role in gene regulation. However, the mechanisms that partition the genome into TADs and the nature of domain boundaries are still poorly understood. We have investigated boundary regions in the Drosophila genome and find that they can be identified as domains of very low H3K27me3. The genome-wide H3K27me3 profile partitions into two states; very low H3K27me3 identifies Depleted (D) domains that contain housekeeping genes and their regulators such as the histone acetyltransferase-containing NSL complex, whereas domains containing moderate-to-high levels of H3K27me3 (Enriched or E domains) are associated with regulated genes, irrespective of whether they are active or inactive. The D domains correlate with the boundaries of TADs and are enriched in a subset of architectural proteins, particularly Chromator, BEAF-32, and Z4/Putzig. However, rather than being clustered at the borders of these domains, these proteins bind throughout the H3K27me3-depleted regions and are much more strongly associated with the transcription start sites of housekeeping genes than with the H3K27me3 domain boundaries. While we have not demonstrated causality, we suggest that the D domain chromatin state, characterised by very low or absent H3K27me3 and established by housekeeping gene regulators, acts to separate topological domains thereby setting up the domain architecture of the genome.

Combgap Promotes Ovarian Niche Development and Chromatin Association of EcR-Binding Regions in BR-C

The development of niches for tissue-specific stem cells is an important aspect of stem cell biology. Determination of niche size and niche numbers during organogenesis involves precise control of gene expression. How this is achieved in the context of a complex chromatin landscape is largely unknown. Here we show that the nuclear protein Combgap (Cg) supports correct ovarian niche formation in Drosophila by controlling ecdysone-Receptor (EcR)- mediated transcription and long-range chromatin contacts in the broad locus (BR-C). Both cg and BR-C promote ovarian growth and the development of niches for germ line stem cells. BR-C levels were lower when Combgap was either reduced or over-expressed, indicating an intricate regulation of the BR-C locus by Combgap. Polytene chromosome stains showed that Cg co-localizes with EcR, the major regulator of BR-C, at the BR-C locus and that EcR binding to chromatin was sensitive to changes in Cg levels. Proximity ligation assay indicated that the two proteins could reside in the same complex. Finally, chromatin conformation analysis revealed that EcR-bound regions within BR-C, which span ~30 KBs, contacted each other. Significantly, these contacts were stabilized in an ecdysone- and Combgap-dependent manner. Together, these results highlight Combgap as a novel regulator of chromatin structure that promotes transcription of ecdysone target genes and ovarian niche formation.

Germ line stem cells (GSCs) supply either eggs or sperm throughout the life- time of many organisms, including mammals. For their function, GSCs require input from somatic niche cells. Understanding how niches form during development is an important initial step in understanding how stem cell units form, and by extension, how they may regenerate. In this work we describe a new function for the chromatin binding protein Combgap in ovarian niche formation of the model organism Drosophila melanogaster. Combgap is required for the correct expression of another factor, BR-C, in somatic ovarian cells. BR-C is one of the central target genes of the steroid hormone ecdysone, and its expression is controlled by the ecdysone receptor (EcR). Interestingly, EcR-enriched regions within the BR-C locus are engaged in long-range contacts that are stabilized by ecdysone in a Combgap-depended manner. We also found that EcR binding to chromatin depends on WT levels of Combgap. BR-C regulates GSC unit establishment, intestinal stem cells, immune responses, and many other processes. Understanding Combgaps’ function in shaping the BR-C chromatin landscape is a first step towards better appreciation of how this important locus is controlled, and the general machinery coupling gene expression to 3D chromatin structure.

CRISPR/Cas9 genome editing throws descriptive 3-D genome folding studies for a loop

CRISPR/Cas9 genome editing studies have recently shed new light into the causal link between the linear DNA sequence and 3-D chromatin architecture. Here we describe current models for the folding of genomes into a nested hierarchy of higher-order structures and discuss new insights into the organizing principles governing genome folding at each length scale. WIREs Syst Biol Med 2016, 8:286-299. doi: 10.1002/wsbm.1338 For further resources related to this article, please visit the WIREs website.

The insulator protein BEAF-32 is required for Hippo pathway activity in the terminal differentiation of neuronal subtypes

The Hippo pathway is crucial for not only normal growth and apoptosis but also cell fate specification during development. What controls Hippo pathway activity during cell fate specification is incompletely understood. In this article, we identify the insulator protein BEAF-32 as a regulator of Hippo pathway activity in Drosophila photoreceptor differentiation. Though morphologically uniform, the fly eye is composed of two subtypes of R8 photoreceptor neurons defined by expression of light-detecting Rhodopsin proteins. In one R8 subtype, active Hippo signaling induces Rhodopsin 6 (Rh6) and represses Rhodopsin 5 (Rh5), whereas in the other subtype, inactive Hippo signaling induces Rh5 and represses Rh6. The activity state of the Hippo pathway in R8 cells is determined by the expression of warts, a core pathway kinase, which interacts with the growth regulator melted in a double-negative feedback loop. We show that BEAF-32 is required for expression of warts and repression of melted Furthermore, BEAF-32 plays a second role downstream of Warts to induce Rh6 and prevent Rh5 fate. BEAF-32 is dispensable for Warts feedback, indicating that BEAF-32 differentially regulates warts and Rhodopsins. Loss of BEAF-32 does not noticeably impair the functions of the Hippo pathway in eye growth regulation. Our study identifies a context-specific regulator of Hippo pathway activity in post-mitotic neuronal fate, and reveals a developmentally specific role for a broadly expressed insulator protein.

CTCF-dependent co-localization of canonical Smad signaling factors at architectural protein binding sites in D. melanogaster

The transforming growth factor β (TGF-β) and bone morphogenetic protein (BMP) pathways transduce extracellular signals into tissue-specific transcriptional responses. During this process, signaling effector Smad proteins translocate into the nucleus to direct changes in transcription, but how and where they localize to DNA remain important questions. We have mapped Drosophila TGF-β signaling factors Mad, dSmad2, Medea, and Schnurri genome-wide in Kc cells and find that numerous sites for these factors overlap with the architectural protein CTCF. Depletion of CTCF by RNAi results in the disappearance of a subset of Smad sites, suggesting Smad proteins localize to CTCF binding sites in a CTCF-dependent manner. Sensitive Smad binding sites are enriched at low occupancy CTCF peaks within topological domains, rather than at the physical domain boundaries where CTCF may function as an insulator. In response to Decapentaplegic, CTCF binding is not significantly altered, whereas Mad, Medea, and Schnurri are redirected from CTCF to non-CTCF binding sites. These results suggest that CTCF participates in the recruitment of Smad proteins to a subset of genomic sites and in the redistribution of these proteins in response to BMP signaling.

Computational Identification of Genomic Features That Influence 3D Chromatin Domain Formation

Recent advances in long-range Hi-C contact mapping have revealed the importance of the 3D structure of chromosomes in gene expression. A current challenge is to identify the key molecular drivers of this 3D structure. Several genomic features, such as architectural proteins and functional elements, were shown to be enriched at topological domain borders using classical enrichment tests. Here we propose multiple logistic regression to identify those genomic features that positively or negatively influence domain border establishment or maintenance. The model is flexible, and can account for statistical interactions among multiple genomic features. Using both simulated and real data, we show that our model outperforms enrichment test and non-parametric models, such as random forests, for the identification of genomic features that influence domain borders. Using Drosophila Hi-C data at a very high resolution of 1 kb, our model suggests that, among architectural proteins, BEAF-32 and CP190 are the main positive drivers of 3D domain borders. In humans, our model identifies well-known architectural proteins CTCF and cohesin, as well as ZNF143 and Polycomb group proteins as positive drivers of domain borders. The model also reveals the existence of several negative drivers that counteract the presence of domain borders including P300, RXRA, BCL11A and ELK1.

Chromosomal DNA is tightly packed up in 3D such that around 2 meters of this long molecule fits into the microscopic nucleus of every cell. The genome packing is not random, but instead structured in 3D domains that are essential to numerous key processes in the cell, such as for the regulation of gene expression or for the replication of DNA. A current challenge is to identify the key molecular drivers of this higher-order chromosome organization. Here we propose a novel computational integrative approach to identify proteins and DNA elements that positively or negatively influence the establishment or maintenance of 3D domains. Analysis of Drosophila data at very high resolution suggests that among architectural proteins, BEAF-32 and CP190 are the main positive drivers of 3D domains. In humans, our results highlight the roles of CTCF, cohesin, ZNF143 and Polycomb group proteins as positive drivers of 3D domains, in contrast to P300, RXRA, BCL11A and ELK1 that act as negative drivers.

Defining the multivalent functions of CTCF from chromatin state and three-dimensional chromatin interactions

CCCTC-binding factor (CTCF) is a multi-functional protein that is assigned various, even contradictory roles in the genome. High-throughput sequencing-based technologies such as ChIP-seq and Hi-C provided us the opportunity to assess the multivalent functions of CTCF in the human genome. The location of CTCF-binding sites with respect to genomic features provides insights into the possible roles of this protein. Here we present the first genome-wide survey and characterization of three important functions of CTCF: enhancer insulator, chromatin barrier and enhancer linker. We developed a novel computational framework to discover the multivalent functions of CTCF based on chromatin state and three-dimensional chromatin architecture. We applied our method to five human cell lines and identified ∼46 000 non-redundant CTCF sites related to the three functions. Disparate effects of these functions on gene expression were found and distinct genomic features of these CTCF sites were characterized in GM12878 cells. Finally, we investigated the cell-type specificities of CTCF sites related to these functions across five cell types. Our study provides new insights into the multivalent functions of CTCF in the human genome.

Genome-Wide Mapping Targets of the Metazoan Chromatin Remodeling Factor NURF Reveals Nucleosome Remodeling at Enhancers, Core Promoters and Gene Insulators

NURF is a conserved higher eukaryotic ISWI-containing chromatin remodeling complex that catalyzes ATP-dependent nucleosome sliding. By sliding nucleosomes, NURF is able to alter chromatin dynamics to control transcription and genome organization. Previous biochemical and genetic analysis of the specificity-subunit of Drosophila NURF (Nurf301/Enhancer of Bithorax (E(bx)) has defined NURF as a critical regulator of homeotic, heat-shock and steroid-responsive gene transcription. It has been speculated that NURF controls pathway specific transcription by co-operating with sequence-specific transcription factors to remodel chromatin at dedicated enhancers. However, conclusive in vivo demonstration of this is lacking and precise regulatory elements targeted by NURF are poorly defined. To address this, we have generated a comprehensive map of in vivo NURF activity, using MNase-sequencing to determine at base pair resolution NURF target nucleosomes, and ChIP-sequencing to define sites of NURF recruitment. Our data show that, besides anticipated roles at enhancers, NURF interacts physically and functionally with the TRF2/DREF basal transcription factor to organize nucleosomes downstream of active promoters. Moreover, we detect NURF remodeling and recruitment at distal insulator sites, where NURF functionally interacts with and co-localizes with DREF and insulator proteins including CP190 to establish nucleosome-depleted domains. This insulator function of NURF is most apparent at subclasses of insulators that mark the boundaries of chromatin domains, where multiple insulator proteins co-associate. By visualizing the complete repertoire of in vivo NURF chromatin targets, our data provide new insights into how chromatin remodeling can control genome organization and regulatory interactions.

Active chromatin and transcription play a key role in chromosome partitioning into topologically associating domains

Recent advances enabled by the Hi-C technique have unraveled many principles of chromosomal folding that were subsequently linked to disease and gene regulation. In particular, Hi-C revealed that chromosomes of animals are organized into topologically associating domains (TADs), evolutionary conserved compact chromatin domains that influence gene expression. Mechanisms that underlie partitioning of the genome into TADs remain poorly understood. To explore principles of TAD folding in Drosophila melanogaster, we performed Hi-C and poly(A)⁺ RNA-seq in four cell lines of various origins (S2, Kc167, DmBG3-c2, and OSC). Contrary to previous studies, we find that regions between TADs (i.e., the inter-TADs and TAD boundaries) in Drosophila are only weakly enriched with the insulator protein dCTCF, while another insulator protein Su(Hw) is preferentially present within TADs. However, Drosophila inter-TADs harbor active chromatin and constitutively transcribed (housekeeping) genes. Accordingly, we find that binding of insulator proteins dCTCF and Su(Hw) predicts TAD boundaries much worse than active chromatin marks do. Interestingly, inter-TADs correspond to decompacted inter-bands of polytene chromosomes, whereas TADs mostly correspond to densely packed bands. Collectively, our results suggest that TADs are condensed chromatin domains depleted in active chromatin marks, separated by regions of active chromatin. We propose the mechanism of TAD self-assembly based on the ability of nucleosomes from inactive chromatin to aggregate, and lack of this ability in acetylated nucleosomal arrays. Finally, we test this hypothesis by polymer simulations and find that TAD partitioning may be explained by different modes of inter-nucleosomal interactions for active and inactive chromatin.

Integrated tRNA, transcript, and protein profiles in response to steroid hormone signaling

The accurate and efficient transfer of genetic information into amino acid sequences is carried out through codon-anticodon interactions between mRNA and tRNA, respectively. In this way, tRNAs function at the interface between gene expression and protein synthesis. Whether tRNA levels are dynamically regulated and to what degree tRNA abundance influences the cellular proteome remains largely unexplored. Here we profile tRNA, transcript and protein levels in Drosophila Kc167 cells, a plasmatocyte cell line that, upon treatment with 20-hydroxyecdysone, differentiates into macrophages. We find that high abundance tRNAs associate with codons that are overrepresented in the Kc167 cell proteome, whereas tRNAs that are in low supply associate with codons that are underrepresented. Ecdysone-induced differentiation of Kc167 cells leads to changes in mRNA codon usage in a manner consistent with the developmental progression of the cell. At both early and late time points, ecdysone treatment concomitantly increases the abundance of tRNAThr(CGU), which decodes a differentiation-associated codon that becomes enriched in the macrophage proteome. These results together suggest that tRNA levels may provide a meaningful regulatory mechanism for defining the cellular proteomic landscape.

Chromatin Insulators and Topological Domains: Adding New Dimensions to 3D Genome Architecture

The spatial organization of metazoan genomes has a direct influence on fundamental nuclear processes that include transcription, replication, and DNA repair. It is imperative to understand the mechanisms that shape the 3D organization of the eukaryotic genomes. Chromatin insulators have emerged as one of the central components of the genome organization tool-kit across species. Recent advancements in chromatin conformation capture technologies have provided important insights into the architectural role of insulators in genomic structuring. Insulators are involved in 3D genome organization at multiple spatial scales and are important for dynamic reorganization of chromatin structure during reprogramming and differentiation. In this review, we will discuss the classical view and our renewed understanding of insulators as global genome organizers. We will also discuss the plasticity of chromatin structure and its re-organization during pluripotency and differentiation and in situations of cellular stress.

Probing long-range interactions by extracting free energies from genome-wide chromosome conformation capture data

Background

A variety of DNA binding proteins are involved in regulating and shaping the packing of chromatin. They aid the formation of loops in the DNA that function to isolate different structural domains. A recent experimental technique, Hi-C, provides a method for determining the frequency of such looping between all distant parts of the genome. Given that the binding locations of many chromatin associated proteins have also been measured, it has been possible to make estimates for their influence on the long-range interactions as measured by Hi-C. However, a challenge in this analysis is the predominance of non-specific contacts that mask out the specific interactions of interest.

Results

We show that transforming the Hi-C contact frequencies into free energies gives a natural method for separating out the distance dependent non-specific interactions. In particular we apply Principal Component Analysis (PCA) to the transformed free energy matrix to identify the dominant modes of interaction. PCA identifies systematic effects as well as high frequency spatial noise in the Hi-C data which can be filtered out. Thus it can be used as a data driven approach for normalizing Hi-C data. We assess this PCA based normalization approach, along with several other normalization schemes, by fitting the transformed Hi-C data using a pairwise interaction model that takes as input the known locations of bound chromatin factors. The result of fitting is a set of predictions for the coupling energies between the various chromatin factors and their effect on the energetics of looping. We show that the quality of the fit can be used as a means to determine how much PCA filtering should be applied to the Hi-C data.

Conclusions

We find that the different normalizations of the Hi-C data vary in the quality of fit to the pairwise interaction model. PCA filtering can improve the fit, and the predicted coupling energies lead to biologically meaningful insights for how various chromatin bound factors influence the stability of DNA loops in chromatin.

Electronic supplementary material

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

The insulator protein CTCF regulates Drosophila steroidogenesis

The steroid hormone ecdysone is a central regulator of insect development. In this report we show that CTCF expression in the prothoracic gland is required for full transcriptional activation of the Halloween genes spookier, shadow and noppera-bo, which encode ecdysone biosynthetic enzymes, and for proper timing of ecdysone-responsive gene expression. Loss of CTCF results in delayed and less synchronized larval development that can only be rescued by feeding larvae with both, the steroid hormone 20-hydroxyecdysone and cholesterol. Moreover, CTCF-knockdown in prothoracic gland cells leads to increased lipid accumulation. In conclusion, the insulator protein CTCF is required for Halloween gene expression and cholesterol homeostasis in ecdysone-producing cells controlling steroidogenesis.

Widespread rearrangement of 3D chromatin organization underlies polycomb-mediated stress-induced silencing

Chromosomes of metazoan organisms are partitioned in the interphase nucleus into discrete topologically associating domains (TADs). Borders between TADs are formed in regions containing active genes and clusters of architectural protein binding sites. The transcription of most genes is repressed after temperature stress in Drosophila. Here we show that temperature stress induces relocalization of architectural proteins from TAD borders to inside TADs, and this is accompanied by a dramatic rearrangement in the 3D organization of the nucleus. TAD border strength declines, allowing for an increase in long-distance inter-TAD interactions. Similar but quantitatively weaker effects are observed upon inhibition of transcription or depletion of individual architectural proteins. Heat shock-induced inter-TAD interactions result in increased contacts among enhancers and promoters of silenced genes, which recruit Pc and form Pc bodies in the nucleolus. These results suggest that the TAD organization of metazoan genomes is plastic and can be reconfigured quickly.

Roles of chromatin insulators in the formation of long-range contacts

Chromatin insulators are factors involved in higher-order, genome-wide organization of chromatin, and play key roles in regulating transcriptional programs. In this review, we discuss recent studies on the diverse composition of insulator complexes, and on the mechanism by which they establish long-range DNA interactions. Particularly, we describe new biophysical methods that allow for the study of the composition of large molecular complexes, and for defining the factors potentially required to establish long-range DNA contacts.

A variably occupied CTCF binding site in the ultrabithorax gene in the Drosophila bithorax complex

Although the majority of genomic binding sites for the insulator protein CCCTC-binding factor (CTCF) are constitutively occupied, a subset show variable occupancy. Such variable sites provide an opportunity to assess context-specific CTCF functions in gene regulation. Here, we have identified a variably occupied CTCF site in the Drosophila Ultrabithorax (Ubx) gene. This site is occupied in tissues where Ubx is active (third thoracic leg imaginal disc) but is not bound in tissues where the Ubx gene is repressed (first thoracic leg imaginal disc). Using chromatin conformation capture, we show that this site preferentially interacts with the Ubx promoter region in the active state. The site lies close to Ubx enhancer elements and is also close to the locations of several gypsy transposon insertions that disrupt Ubx expression, leading to the bx mutant phenotype. gypsy insertions carry the Su(Hw)-dependent gypsy insulator and were found to affect both CTCF binding at the variable site and the chromatin topology. This suggests that insertion of the gypsy insulator in this region interferes with CTCF function and supports a model for the normal function of the variable CTCF site as a chromatin loop facilitator, promoting interaction between Ubx enhancers and the Ubx transcription start site.

Expanding the roles of chromatin insulators in nuclear architecture, chromatin organization and genome function

Of the numerous classes of elements involved in modulating eukaryotic chromosome structure and function, chromatin insulators arguably remain the most poorly understood in their contribution to these processes in vivo. Indeed, our view of chromatin insulators has evolved dramatically since their chromatin boundary and enhancer blocking properties were elucidated roughly a quarter of a century ago as a result of recent genome-wide, high-throughput methods better suited to probing the role of these elements in their native genomic contexts. The overall theme that has emerged from these studies is that chromatin insulators function as general facilitators of higher-order chromatin loop structures that exert both physical and functional constraints on the genome. In this review, we summarize the result of recent work that supports this idea as well as a number of other studies linking these elements to a diverse array of nuclear processes, suggesting that chromatin insulators exert master control over genome organization and behavior.

Chromatin insulator factors involved in long-range DNA interactions and their role in the folding of the Drosophila genome

Chromatin insulators are genetic elements implicated in the organization of chromatin and the regulation of transcription. In Drosophila, different insulator types were characterized by their locus-specific composition of insulator proteins and co-factors. Insulators mediate specific long-range DNA contacts required for the three dimensional organization of the interphase nucleus and for transcription regulation, but the mechanisms underlying the formation of these contacts is currently unknown. Here, we investigate the molecular associations between different components of insulator complexes (BEAF32, CP190 and Chromator) by biochemical and biophysical means, and develop a novel single-molecule assay to determine what factors are necessary and essential for the formation of long-range DNA interactions. We show that BEAF32 is able to bind DNA specifically and with high affinity, but not to bridge long-range interactions (LRI). In contrast, we show that CP190 and Chromator are able to mediate LRI between specifically-bound BEAF32 nucleoprotein complexes in vitro. This ability of CP190 and Chromator to establish LRI requires specific contacts between BEAF32 and their C-terminal domains, and dimerization through their N-terminal domains. In particular, the BTB/POZ domains of CP190 form a strict homodimer, and its C-terminal domain interacts with several insulator binding proteins. We propose a general model for insulator function in which BEAF32/dCTCF/Su(HW) provide DNA specificity (first layer proteins) whereas CP190/Chromator are responsible for the physical interactions required for long-range contacts (second layer). This network of organized, multi-layer interactions could explain the different activities of insulators as chromatin barriers, enhancer blockers, and transcriptional regulators, and suggest a general mechanism for how insulators may shape the organization of higher-order chromatin during cell division.

dREAM co-operates with insulator-binding proteins and regulates expression at divergently paired genes

dREAM complexes represent the predominant form of E2F/RBF repressor complexes in Drosophila. dREAM associates with thousands of sites in the fly genome but its mechanism of action is unknown. To understand the genomic context in which dREAM acts we examined the distribution and localization of Drosophila E2F and dREAM proteins. Here we report a striking and unexpected overlap between dE2F2/dREAM sites and binding sites for the insulator-binding proteins CP190 and Beaf-32. Genetic assays show that these components functionally co-operate and chromatin immunoprecipitation experiments on mutant animals demonstrate that dE2F2 is important for association of CP190 with chromatin. dE2F2/dREAM binding sites are enriched at divergently transcribed genes, and the majority of genes upregulated by dE2F2 depletion represent the repressed half of a differentially expressed, divergently transcribed pair of genes. Analysis of mutant animals confirms that dREAM and CP190 are similarly required for transcriptional integrity at these gene pairs and suggest that dREAM functions in concert with CP190 to establish boundaries between repressed/activated genes. Consistent with the idea that dREAM co-operates with insulator-binding proteins, genomic regions bound by dREAM possess enhancer-blocking activity that depends on multiple dREAM components. These findings suggest that dREAM functions in the organization of transcriptional domains.

Insulator function and topological domain border strength scale with architectural protein occupancy

BACKGROUND

Chromosome conformation capture studies suggest that eukaryotic genomes are organized into structures called topologically associating domains. The borders of these domains are highly enriched for architectural proteins with characterized roles in insulator function. However, a majority of architectural protein binding sites localize within topological domains, suggesting sites associated with domain borders represent a functionally different subclass of these regulatory elements. How topologically associating domains are established and what differentiates border-associated from non-border architectural protein binding sites remain unanswered questions.

RESULTS

By mapping the genome-wide target sites for several Drosophila architectural proteins, including previously uncharacterized profiles for TFIIIC and SMC-containing condensin complexes, we uncover an extensive pattern of colocalization in which architectural proteins establish dense clusters at the borders of topological domains. Reporter-based enhancer-blocking insulator activity as well as endogenous domain border strength scale with the occupancy level of architectural protein binding sites, suggesting co-binding by architectural proteins underlies the functional potential of these loci. Analyses in mouse and human stem cells suggest that clustering of architectural proteins is a general feature of genome organization, and conserved architectural protein binding sites may underlie the tissue-invariant nature of topologically associating domains observed in mammals.

CONCLUSIONS

We identify a spectrum of architectural protein occupancy that scales with the topological structure of chromosomes and the regulatory potential of these elements. Whereas high occupancy architectural protein binding sites associate with robust partitioning of topologically associating domains and robust insulator function, low occupancy sites appear reserved for gene-specific regulation within topological domains.

Insulators recruit histone methyltransferase dMes4 to regulate chromatin of flanking genes

Chromosomal domains in Drosophila are marked by the insulator-binding proteins (IBPs) dCTCF/Beaf32 and cofactors that participate in regulating long-range interactions. Chromosomal borders are further enriched in specific histone modifications, yet the role of histone modifiers and nucleosome dynamics in this context remains largely unknown. Here, we show that IBP depletion impairs nucleosome dynamics specifically at the promoters and coding sequence of genes flanked by IBP binding sites. Biochemical purification identifies the H3K36 histone methyltransferase NSD/dMes-4 as a novel IBP cofactor, which specifically co-regulates the chromatin accessibility of hundreds of genes flanked by dCTCF/Beaf32. NSD/dMes-4 presets chromatin before the recruitment of transcriptional activators including DREF that triggers Set2/Hypb-dependent H3K36 trimethylation, nucleosome positioning, and RNA splicing. Our results unveil a model for how IBPs regulate nucleosome dynamics and gene expression through NSD/dMes-4, which may regulate H3K27me3 spreading. Our data uncover how IBPs dynamically regulate chromatin organization depending on distinct cofactors.

The RNA-binding protein Rumpelstiltskin antagonizes gypsy chromatin insulator function in a tissue-specific manner

Chromatin insulators are DNA-protein complexes that are situated throughout the genome that are proposed to contribute to higher-order organization and demarcation into distinct transcriptional domains. Mounting evidence in different species implicates RNA and RNA-binding proteins as regulators of chromatin insulator activities. Here, we identify the Drosophila hnRNP M homolog Rumpelstiltskin (Rump) as an antagonist of gypsy chromatin insulator enhancer-blocking and barrier activities. Despite ubiquitous expression of Rump, decreasing Rump levels leads to improvement of barrier activity only in tissues outside of the central nervous system (CNS). Furthermore, rump mutants restore insulator body localization in an insulator mutant background only in non-CNS tissues. Rump associates physically with core gypsy insulator proteins, and chromatin immunoprecipitation and sequencing analysis of Rump demonstrates extensive colocalization with a subset of insulator sites across the genome. The genome-wide binding profile and tissue specificity of Rump contrast with that of Shep, a recently identified RNA-binding protein that antagonizes gypsy insulator activity primarily in the CNS. Our findings indicate parallel roles for RNA-binding proteins in mediating tissue-specific regulation of chromatin insulator activity.

CTCF: an architectural protein bridging genome topology and function

The eukaryotic genome is organized in the three-dimensional nuclear space in a specific manner that is both a cause and a consequence of its function. This organization is partly established by a special class of architectural proteins, of which CCCTC-binding factor (CTCF) is the best characterized. Although CTCF has been assigned various roles that are often contradictory, new results now help to draw a unifying model to explain the many functions of this protein. CTCF creates boundaries between topologically associating domains in chromosomes and, within these domains, facilitates interactions between transcription regulatory sequences. Thus, CTCF links the architecture of the genome to its function.

Ibf1 and Ibf2 are novel CP190-interacting proteins required for insulator function

Insulators are DNA-protein complexes that play a central role in chromatin organization and regulation of gene expression. In Drosophila different proteins, dCTCF, Su(Hw), and BEAF bind to specific subsets of insulators most of them having in common CP190. It has been shown that there are a number of CP190-binding sites that are not shared with any other known insulator protein, suggesting that other proteins could cooperate with CP190 to regulate insulator activity. Here we report on the identification of two previously uncharacterized proteins as CP190-interacting proteins, that we have named Ibf1 and Ibf2. These proteins localize at insulator bodies and associate with chromatin at CP190-binding sites throughout the genome. We also show that Ibf1 and Ibf2 are DNA-binding proteins that form hetero-oligomers that mediate CP190 binding to chromatin. Moreover, Ibf1 and Ibf2 are necessary for insulator activity in enhancer-blocking assays and Ibf2 null mutation cause a homeotic phenotype. Taken together our data reveal a novel pathway of CP190 recruitment to chromatin that is required for insulator activity.

Chromatin immunoprecipitation indirect peaks highlight long-range interactions of insulator proteins and Pol II pausing

Eukaryotic chromosomes are partitioned into topologically associating domains (TADs) that are demarcated by distinct insulator-binding proteins (IBPs) in Drosophila. Whether IBPs regulate specific long-range contacts and how this may impact gene expression remains unclear. Here we identify "indirect peaks" of multiple IBPs that represent their distant sites of interactions through long-range contacts. Indirect peaks depend on protein-protein interactions among multiple IBPs and their common cofactors, including CP190, as confirmed by high-resolution analyses of long-range contacts. Mutant IBPs unable to interact with CP190 impair long-range contacts as well as the expression of hundreds of distant genes that are specifically flanked by indirect peaks. Regulation of distant genes strongly correlates with RNAPII pausing, highlighting how this key transcriptional stage may trap insulator-based long-range interactions. Our data illustrate how indirect peaks may decipher gene regulatory networks through specific long-range interactions.