Architectural proteins Pita, Zw5,and ZIPIC contain homodimerization domain and support specific long-range interactions in Drosophila

The homie insulator has sub-elements with different insulating and long-range pairing properties

Chromatin insulators are major determinants of chromosome architecture. Specific architectures induced by insulators profoundly influence nuclear processes, including how enhancers and promoters interact over long distances and between homologous chromosomes. Insulators can pair with copies of themselves in trans to facilitate homolog pairing. They can also pair with other insulators, sometimes with great specificity, inducing long-range chromosomal loops. Contrary to their canonical function of enhancer blocking, these loops can bring distant enhancers and promoters together to activate gene expression, while at the same time blocking other interactions in cis. The details of these effects depend on the choice of pairing partner, and on the orientation specificity of pairing, implicating the 3D architecture as a major functional determinant. Here, we dissect the homie insulator from the Drosophila even skipped (eve) locus, to understand its substructure. We test pairing function based on homie-carrying transgenes interacting with endogenous eve. The assay is sensitive to both pairing strength and orientation. Using this assay, we found that a Su(Hw) binding site in homie is required for efficient long-range interaction, although some activity remains without it. This binding site also contributes to the canonical insulator activities of enhancer blocking and barrier function. Based on this and other results from our functional dissection, each of the canonical insulator activities, chromosomal loop formation, enhancer blocking, and barrier activity, are partially separable. Our results show the complexity inherent in insulator functions, which can be provided by an array of different proteins with both shared and distinct properties.

The homie insulator has sub-elements with different insulating and long-range pairing properties

Chromatin insulators are major determinants of chromosome architecture. Specific architectures induced by insulators profoundly influence nuclear processes, including how enhancers and promoters interact over long distances and between homologous chromosomes. Insulators can pair with copies of themselves in trans to facilitate homolog pairing. They can also pair with other insulators, sometimes with great specificity, inducing long-range chromosomal loops. Contrary to their canonical function of enhancer blocking, these loops can bring distant enhancers and promoters together to activate gene expression, while at the same time blocking other interactions in cis. The details of these effects depend on the choice of pairing partner, and on the orientation specificity of pairing, implicating the 3-dimensional architecture as a major functional determinant. Here we dissect the homie insulator from the Drosophila even skipped (eve) locus, to understand its substructure. We test pairing function based on homie-carrying transgenes interacting with endogenous eve. The assay is sensitive to both pairing strength and orientation. Using this assay, we found that a Su(Hw) binding site in homie is required for efficient long-range interaction, although some activity remains without it. This binding site also contributes to the canonical insulator activities of enhancer blocking and barrier function. Based on this and other results from our functional dissection, each of the canonical insulator activities, chromosomal loop formation, enhancer blocking, and barrier activity, are partially separable. Our results show the complexity inherent in insulator functions, which can be provided by an array of different proteins with both shared and distinct properties.

Chromosome structure in Drosophila is determined by boundary pairing not loop extrusion

Two different models have been proposed to explain how the endpoints of chromatin looped domains ('TADs') in eukaryotic chromosomes are determined. In the first, a cohesin complex extrudes a loop until it encounters a boundary element roadblock, generating a stem-loop. In this model, boundaries are functionally autonomous: they have an intrinsic ability to halt the movement of incoming cohesin complexes that is independent of the properties of neighboring boundaries. In the second, loops are generated by boundary:boundary pairing. In this model, boundaries are functionally non-autonomous, and their ability to form a loop depends upon how well they match with their neighbors. Moreover, unlike the loop-extrusion model, pairing interactions can generate both stem-loops and circle-loops. We have used a combination of MicroC to analyze how TADs are organized, and experimental manipulations of the even skipped TAD boundary, homie, to test the predictions of the 'loop-extrusion' and the 'boundary-pairing' models. Our findings are incompatible with the loop-extrusion model, and instead suggest that the endpoints of TADs in flies are determined by a mechanism in which boundary elements physically pair with their partners, either head-to-head or head-to-tail, with varying degrees of specificity. Although our experiments do not address how partners find each other, the mechanism is unlikely to require loop extrusion.

Su(Hw) interacts with Combgap to establish long-range chromatin contacts

BACKGROUND

Insulator-binding proteins (IBPs) play a critical role in genome architecture by forming and maintaining contact domains. While the involvement of several IBPs in organising chromatin architecture in Drosophila has been described, the specific contribution of the Suppressor of Hairy wings (Su(Hw)) insulator-binding protein to genome topology remains unclear.

RESULTS

In this study, we provide evidence for the existence of long-range interactions between chromatin bound Su(Hw) and Combgap, which was first characterised as Polycomb response elements binding protein. Loss of Su(Hw) binding to chromatin results in the disappearance of Su(Hw)-Combgap long-range interactions and in a decrease in spatial self-interactions among a subset of Su(Hw)-bound genome sites. Our findings suggest that Su(Hw)-Combgap long-range interactions are associated with active chromatin rather than Polycomb-directed repression. Furthermore, we observe that the majority of transcription start sites that are down-regulated upon loss of Su(Hw) binding to chromatin are located within 2 kb of Combgap peaks and exhibit Su(Hw)-dependent changes in Combgap and transcriptional regulators' binding.

CONCLUSIONS

This study demonstrates that Su(Hw) insulator binding protein can form long-range interactions with Combgap, Polycomb response elements binding protein, and that these interactions are associated with active chromatin factors rather than with Polycomb dependent repression.

The N-terminal dimerization domains of human and Drosophila CTCF have similar functionality

BACKGROUND

CTCF is highly likely to be the ancestor of proteins that contain large clusters of C2H2 zinc finger domains, and its conservation is observed across most bilaterian organisms. In mammals, CTCF is the primary architectural protein involved in organizing chromosome topology and mediating enhancer-promoter interactions over long distances. In Drosophila, CTCF (dCTCF) cooperates with other architectural proteins to establish long-range interactions and chromatin boundaries. CTCFs of various organisms contain an unstructured N-terminal dimerization domain (DD) and clusters comprising eleven zinc-finger domains of the C2H2 type. The Drosophila (dCTCF) and human (hCTCF) CTCFs share sequence homology in only five C2H2 domains that specifically bind to a conserved 15 bp motif.

RESULTS

Previously, we demonstrated that CTCFs from different organisms carry unstructured N-terminal dimerization domains (DDs) that lack sequence homology. Here we used the CTCFattP(mCh) platform to introduce desired changes in the Drosophila CTCF gene and generated a series of transgenic lines expressing dCTCF with different variants of the N-terminal domain. Our findings revealed that the functionality of dCTCF is significantly affected by the deletion of the N-terminal DD. Additionally, we observed a strong impact on the binding of the dCTCF mutant to chromatin upon deletion of the DD. However, chromatin binding was restored in transgenic flies expressing a chimeric CTCF protein with the DD of hCTCF. Although the chimeric protein exhibited lower expression levels than those of the dCTCF variants, it efficiently bound to chromatin similarly to the wild type (wt) protein.

CONCLUSIONS

Our findings suggest that one of the evolutionarily conserved functions of the unstructured N-terminal dimerization domain is to recruit dCTCF to its genomic sites in vivo.

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.

lncRNA read-through regulates the BX-C insulator Fub-1

Though long non-coding RNAs (lncRNAs) represent a substantial fraction of the Pol II transcripts in multicellular animals, only a few have known functions. Here we report that the blocking activity of the Bithorax complex (BX-C) Fub-1 boundary is segmentally regulated by its own lncRNA. The Fub-1 boundary is located between the Ultrabithorax (Ubx) gene and the bxd/pbx regulatory domain, which is responsible for regulating Ubx expression in parasegment PS6/segment A1. Fub-1 consists of two hypersensitive sites, HS1 and HS2. HS1 is an insulator while HS2 functions primarily as an lncRNA promoter. To activate Ubx expression in PS6/A1, enhancers in the bxd/pbx domain must be able to bypass Fub-1 blocking activity. We show that the expression of the Fub-1 lncRNAs in PS6/A1 from the HS2 promoter inactivates Fub-1 insulating activity. Inactivation is due to read-through as the HS2 promoter must be directed toward HS1 to disrupt blocking.

Boundary Bypass Activity in the Abdominal-B Region of the Drosophila Bithorax Complex is Position Dependent and Regulated

Expression of Abdominal-B ( Abd-B ) in abdominal segments A5 - A8 is controlled by four regulatory domains, iab-5 - iab-8 . Each domain has an initiator element (which sets the activity state), elements that maintain this state and tissue-specific enhancers. To ensure their functional autonomy, each domain is bracketed by boundary elements ( Mcp , Fab-7 , Fab-7 and Fab-8 ). In addition to blocking crosstalk between adjacent regulatory domains, the Fab boundaries must also have bypass activity so the relevant regulatory domains can "jump over" intervening boundaries and activate the Abd-B promoter. In the studies reported here we have investigated the parameters governing bypass activity. We find that the bypass elements in the Fab-7 and Fab-8 boundaries must be located in the regulatory domain that is responsible for driving Abd-B expression. We suggest that bypass activity may also be subject to regulation.

Summary Statement

Boundaries separating Abd-B regulatory domains block crosstalk between domains and mediate their interactions with Abd-B . The latter function is location but not orientation dependent.

CTCF, BEAF-32, and CP190 are not required for the establishment of TADs in early Drosophila embryos but have locus-specific roles

The boundaries of topologically associating domains (TADs) are delimited by insulators and/or active promoters; however, how they are initially established during embryogenesis remains unclear. Here, we examined this during the first hours of Drosophila embryogenesis. DNA-FISH confirms that intra-TAD pairwise proximity is established during zygotic genome activation (ZGA) but with extensive cell-to-cell heterogeneity. Most newly formed boundaries are occupied by combinations of CTCF, BEAF-32, and/or CP190. Depleting each insulator individually from chromatin revealed that TADs can still establish, although with lower insulation, with a subset of boundaries (~10%) being more dependent on specific insulators. Some weakened boundaries have aberrant gene expression due to unconstrained enhancer activity. However, the majority of misexpressed genes have no obvious direct relationship to changes in domain-boundary insulation. Deletion of an active promoter (thereby blocking transcription) at one boundary had a greater impact than deleting the insulator-bound region itself. This suggests that cross-talk between insulators and active promoters and/or transcription might reinforce domain boundary insulation during embryogenesis.

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.

Dimerization Activity of a Disordered N-Terminal Domain from Drosophila CLAMP Protein

In Drosophila melanogaster, CLAMP is an essential zinc-finger transcription factor that is involved in chromosome architecture and functions as an adaptor for the dosage compensation complex. Most of the known Drosophila architectural proteins have structural N-terminal homodimerization domains that facilitate distance interactions. Because CLAMP performs architectural functions, we tested its N-terminal region for the presence of a homodimerization domain. We used a yeast two-hybrid assay and biochemical studies to demonstrate that the adjacent N-terminal region between 46 and 86 amino acids is capable of forming homodimers. This region is conserved in CLAMP orthologs from most insects, except Hymenopterans. Biophysical techniques, including nuclear magnetic resonance (NMR) and small-angle X-ray scattering (SAXS), suggested that this domain lacks secondary structure and has features of intrinsically disordered regions despite the fact that the protein structure prediction algorithms suggested the presence of beta-sheets. The dimerization domain is essential for CLAMP functions in vivo because its deletion results in lethality. Thus, CLAMP is the second architectural protein after CTCF that contains an unstructured N-terminal dimerization domain.

Quantitative-enhancer-FACS-seq (QeFS) reveals epistatic interactions among motifs within transcriptional enhancers in developing Drosophila tissue

Understanding the contributions of transcription factor DNA binding sites to transcriptional enhancers is a significant challenge. We developed Quantitative enhancer-FACS-Seq for highly parallel quantification of enhancer activities from a genomically integrated reporter in Drosophila melanogaster embryos. We investigate the contributions of the DNA binding motifs of four poorly characterized TFs to the activities of twelve embryonic mesodermal enhancers. We measure quantitative changes in enhancer activity and discover a range of epistatic interactions among the motifs, both synergistic and alleviating. We find that understanding the regulatory consequences of TF binding motifs requires that they be investigated in combination across enhancer contexts.

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.

Study of the N-Terminal Domain Homodimerization in Human Proteins with Zinc Finger Clusters

CTCF belongs to a large family of transcription factors with clusters of C2H2-type zinc finger domains (C2H2 proteins) and is a main architectural protein in mammals. Human CTCF has a homodimerizing unstructured domain at the N-terminus which is involved in long-distance interactions. To test the presence of similar N-terminal domains in other human C2H2 proteins, a yeast two-hybrid system was used. In total, the ability of unstructured N-terminal domains to homodimerize was investigated for six human C2H2 proteins with an expression profile similar to CTCF. The data indicate the lack of the homodimerization ability of these domains. On the other hand, three C2H2 proteins containing the structured domain DUF3669 at the N-terminus demonstrated homo- and heterodimerization activity.

Mechanism and functional role of the interaction between CP190 and the architectural protein Pita in Drosophila melanogaster

Background

Pita is required for Drosophila development and binds specifically to a long motif in active promoters and insulators. Pita belongs to the Drosophila family of zinc-finger architectural proteins, which also includes Su(Hw) and the conserved among higher eukaryotes CTCF. The architectural proteins maintain the active state of regulatory elements and the long-distance interactions between them. In particular, Pita is involved in the formation of several boundaries between regulatory domains that controlled the expression of three hox genes in the Bithorax complex (BX-C). The CP190 protein is recruited to chromatin through interaction with the architectural proteins.

Results

Using in vitro pull-down analysis, we precisely mapped two unstructured regions of Pita that interact with the BTB domain of CP190. Then we constructed transgenic lines expressing the Pita protein of the wild-type and mutant variants lacking CP190-interacting regions. We have demonstrated that CP190-interacting region of the Pita can maintain nucleosome-free open chromatin and is critical for Pita-mediated enhancer blocking activity in BX-C. At the same time, interaction with CP190 is not required for the in vivo function of the mutant Pita protein, which binds to the same regions of the genome as the wild-type protein. Unexpectedly, we found that CP190 was still associated with the most of genome regions bound by the mutant Pita protein, which suggested that other architectural proteins were continuing to recruit CP190 to these regions.

Conclusions

The results directly demonstrate role of CP190 in insulation and support a model in which the regulatory elements are composed of combinations of binding sites that interact with several architectural proteins with similar functions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00391-x.

Structural basis of diversity and homodimerization specificity of zinc-finger-associated domains in Drosophila

In arthropods, zinc finger-associated domains (ZADs) are found at the N-termini of many DNA-binding proteins with tandem arrays of Cys2-His2 zinc fingers (ZAD-C2H2 proteins). ZAD-C2H2 proteins undergo fast evolutionary lineage-specific expansion and functional diversification. Here, we show that all ZADs from Drosophila melanogaster form homodimers, but only certain ZADs with high homology can also heterodimerize. CG2712, for example, is unable to heterodimerize with its paralog, the previously characterized insulator protein Zw5, with which it shares 46% homology. We obtained a crystal structure of CG2712 protein's ZAD domain that, in spite of a low sequence homology, has similar spatial organization with the only known ZAD structure (from Grauzone protein). Steric clashes prevented the formation of heterodimers between Grauzone and CG2712 ZADs. Using detailed structural analysis, site-directed mutagenesis, and molecular dynamics simulations, we demonstrated that rapid evolutionary acquisition of interaction specificity was mediated by the more energy-favorable formation of homodimers in comparison to heterodimers, and that this specificity was achieved by multiple amino acid substitutions resulting in the formation or breaking of stabilizing interactions. We speculate that specific homodimerization of ZAD-C2H2 proteins is important for their architectural role in genome organization.

Mechanisms of Enhancer-Promoter Interactions in Higher Eukaryotes

In higher eukaryotes, enhancers determine the activation of developmental gene transcription in specific cell types and stages of embryogenesis. Enhancers transform the signals produced by various transcription factors within a given cell, activating the transcription of the targeted genes. Often, developmental genes can be associated with dozens of enhancers, some of which are located at large distances from the promoters that they regulate. Currently, the mechanisms underlying specific distance interactions between enhancers and promoters remain poorly understood. This review briefly describes the properties of enhancers and discusses the mechanisms of distance interactions and potential proteins involved in this process.

CTCF As an Example of DNA-Binding Transcription Factors Containing Clusters of C2H2-Type Zinc Fingers

In mammals, most of the boundaries of topologically associating domains and all well-studied insulators are rich in binding sites for the CTCF protein. According to existing experimental data, CTCF is a key factor in the organization of the architecture of mammalian chromosomes. A characteristic feature of the CTCF is that the central part of the protein contains a cluster consisting of eleven domains of C2H2-type zinc fingers, five of which specifically bind to a long DNA sequence conserved in most animals. The class of transcription factors that carry a cluster of C2H2-type zinc fingers consisting of five or more domains (C2H2 proteins) is widely represented in all groups of animals. The functions of most C2H2 proteins still remain unknown. This review presents data on the structure and possible functions of these proteins, using the example of the vertebrate CTCF protein and several well- characterized C2H2 proteins in Drosophila and mammals.

Innovation of heterochromatin functions drives rapid evolution of essential ZAD-ZNF genes in Drosophila

Contrary to dogma, evolutionarily young and dynamic genes can encode essential functions. We find that evolutionarily dynamic ZAD-ZNF genes, which encode the most abundant class of insect transcription factors, are more likely to encode essential functions in Drosophila melanogaster than ancient, conserved ZAD-ZNF genes. We focus on the Nicknack ZAD-ZNF gene, which is evolutionarily young, poorly retained in Drosophila species, and evolves under strong positive selection. Yet we find that it is necessary for larval development in D. melanogaster. We show that Nicknack encodes a heterochromatin-localizing protein like its paralog Oddjob, also an evolutionarily dynamic yet essential ZAD-ZNF gene. We find that the divergent D. simulans Nicknack protein can still localize to D. melanogaster heterochromatin and rescue viability of female but not male Nicknack-null D. melanogaster. Our findings suggest that innovation for rapidly changing heterochromatin functions might generally explain the essentiality of many evolutionarily dynamic ZAD-ZNF genes in insects.

The Functions and Mechanisms of Action of Insulators in the Genomes of Higher Eukaryotes

The mechanisms underlying long-range interactions between chromatin regions and the principles of chromosomal architecture formation are currently under extensive scrutiny. A special class of regulatory elements known as insulators is believed to be involved in the regulation of specific long-range interactions between enhancers and promoters. This review focuses on the insulators of Drosophila and mammals, and it also briefly characterizes the proteins responsible for their functional activity. It was initially believed that the main properties of insulators are blocking of enhancers and the formation of independent transcription domains. We present experimental data proving that the chromatin loops formed by insulators play only an auxiliary role in enhancer blocking. The review also discusses the mechanisms involved in the formation of topologically associating domains and their role in the formation of the chromosomal architecture and regulation of gene transcription.

Adaptation of gene loci to heterochromatin in the course of Drosophila evolution is associated with insulator proteins

Pericentromeric heterochromatin is generally composed of repetitive DNA forming a transcriptionally repressive environment. Dozens of genes were embedded into pericentromeric heterochromatin during evolution of Drosophilidae lineage while retaining activity. However, factors that contribute to insusceptibility of gene loci to transcriptional silencing remain unknown. Here, we find that the promoter region of genes that can be embedded in both euchromatin and heterochromatin exhibits a conserved structure throughout the Drosophila phylogeny and carries motifs for binding of certain chromatin remodeling factors, including insulator proteins. Using ChIP-seq data, we demonstrate that evolutionary gene relocation between euchromatin and pericentric heterochromatin occurred with preservation of sites of insulation of BEAF-32 in evolutionarily distant species, i.e. D. melanogaster and D. virilis. Moreover, promoters of virtually all protein-coding genes located in heterochromatin in D. melanogaster are enriched with insulator proteins BEAF-32, GAF and dCTCF. Applying RNA-seq of a BEAF-32 mutant, we show that the impairment of BEAF-32 function has a complex effect on gene expression in D. melanogaster, affecting even those genes that lack BEAF-32 association in their promoters. We propose that conserved intrinsic properties of genes, such as sites of insulation near the promoter regions, may contribute to adaptation of genes to the heterochromatic environment and, hence, facilitate the evolutionary relocation of genes loci between euchromatin and heterochromatin.

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.

A Comprehensive Drosophila melanogaster Transcription Factor Interactome

Combinatorial interactions among transcription factors (TFs) play essential roles in generating gene expression specificity and diversity in metazoans. Using yeast 2-hybrid (Y2H) assays on nearly all sequence-specific Drosophila TFs, we identified 1,983 protein-protein interactions (PPIs), more than doubling the number of currently known PPIs among Drosophila TFs. For quality assessment, we validated a subset of our interactions using MITOMI and bimolecular fluorescence complementation assays. We combined our interactome with prior PPI data to generate an integrated Drosophila TF-TF binary interaction network. Our analysis of ChIP-seq data, integrating PPI and gene expression information, uncovered different modes by which interacting TFs are recruited to DNA. We further demonstrate the utility of our Drosophila interactome in shedding light on human TF-TF interactions. This study reveals how TFs interact to bind regulatory elements in vivo and serves as a resource of Drosophila TF-TF binary PPIs for understanding tissue-specific gene regulation.

Combinatorial regulation by transcription factors (TFs) is one mechanism for achieving condition and tissue-specific gene regulation. Shokri et al. mapped TF-TF interactions between most Drosophila TFs, reporting a comprehensive TF-TF network integrated with previously known interactions. They used this network to discern distinct TF-DNA binding modes.

The insulator functions of the Drosophila polydactyl C2H2 zinc finger protein CTCF: Necessity versus sufficiency

In mammals, a C2H2 zinc finger (C2H2) protein, CTCF, acts as the master regulator of chromosomal architecture and of the expression of Hox gene clusters. Like mammalian CTCF, the Drosophila homolog, dCTCF, localizes to boundaries in the bithorax complex (BX-C). Here, we have determined the minimal requirements for the assembly of a functional boundary by dCTCF and two other C2H2 zinc finger proteins, Pita and Su(Hw). Although binding sites for these proteins are essential for the insulator activity of BX-C boundaries, these binding sites alone are insufficient to create a functional boundary. dCTCF cannot effectively bind to a single recognition sequence in chromatin or generate a functional insulator without the help of additional proteins. In addition, for boundary elements in BX-C at least four binding sites for dCTCF or the presence of additional DNA binding factors is required to generate a functional insulator.

N-terminal domain of the architectural protein CTCF has similar structural organization and ability to self-association in bilaterian organisms

CTCF is the main architectural protein found in most of the examined bilaterian organisms. The cluster of the C2H2 zinc-finger domains involved in recognition of long DNA-binding motif is only part of the protein that is evolutionarily conserved, while the N-terminal domain (NTD) has different sequences. Here, we performed biophysical characterization of CTCF NTDs from various species representing all major phylogenetic clades of higher metazoans. With the exception of Drosophilides, the N-terminal domains of CTCFs show an unstructured organization and absence of folded regions in vitro. In contrast, NTDs of Drosophila melanogaster and virilis CTCFs contain unstructured folded regions that form tetramers and dimers correspondingly in vitro. Unexpectedly, most NTDs are able to self-associate in the yeast two-hybrid and co-immunoprecipitation assays. These results suggest that NTDs of CTCFs might contribute to the organization of CTCF-mediated long-distance interactions and chromosomal architecture.

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.

Complete reconstitution of bypass and blocking functions in a minimal artificial Fab-7 insulator from Drosophila bithorax complex

Boundaries in the bithorax complex (BX-C) delimit autonomous regulatory domains that drive parasegment-specific expression of the Hox genes Ubx, abd-A, and Abd-B The Fab-7 boundary is located between the iab-6 and iab-7 domains and has two key functions: blocking cross-talk between these domains and at the same time promoting communication (boundary bypass) between iab-6 and the Abd-B promoter. Using a replacement strategy, we found that multimerized binding sites for the architectural proteins Pita, Su(Hw), and dCTCF function as conventional insulators and block cross-talk between the iab-6 and iab-7 domains; however, they lack bypass activity, and iab-6 is unable to regulate Abd-B Here we show that an ∼200-bp sequence of dHS1 from the Fab-7 boundary rescues the bypass defects of these multimerized binding sites. The dHS1 sequence is bound in embryos by a large multiprotein complex, Late Boundary Complex (LBC), that contains the zinc finger proteins CLAMP and GAF. Using deletions and mutations in critical GAGAG motifs, we show that bypass activity correlates with the efficiency of recruitment of LBC components CLAMP and GAF to the artificial boundary. These results indicate that LBC orchestrates long-distance communication between the iab-6 regulatory domain and the Abd-B gene, while the Pita, Su(Hw), and dCTCF proteins function to block local cross-talk between the neighboring regulatory domains iab-6 and iab-7.

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.

Functional dissection of the developmentally restricted BEN domain chromatin boundary factor Insensitive

Background

Boundaries in the Drosophila bithorax complex delimit autonomous regulatory domains that activate the parasegment (PS)-specific expression of homeotic genes. The Fab-7 boundary separates the iab-6 and iab-7 regulatory domains that control Abd-B expression in PS11 and PS12. This boundary is composed of multiple functionally redundant elements and has two key activities: it blocks crosstalk between iab-6 and iab-7 and facilitates boundary bypass.

Results

Here, we have used a structure–function approach to elucidate the biochemical properties and the in vivo activities of a conserved BEN domain protein, Insensitive, that is associated with Fab-7. Our biochemical studies indicate that in addition to the C-terminal BEN DNA-binding domain, Insv has two domains that mediate multimerization: one is a coiled-coil domain in the N-terminus, and the other is next to the BEN domain. These multimerization domains enable Insv to bind simultaneously to two canonical 8-bp recognition motifs, as well as to a ~ 100-bp non-canonical recognition sequence. They also mediate the assembly of higher-order multimers in the presence of DNA. Transgenic proteins lacking the N-terminal coiled-coil domain are compromised for boundary function in vivo. We also show that Insv interacts directly with CP190, a protein previously implicated in the boundary functions of several DNA-binding proteins, including Su(Hw) and dCTCF. While CP190 interaction is required for Insv binding to a subset of sites on polytene chromosomes, it has only a minor role in the boundary activity of Insv in the context of Fab-7.

Conclusions

The subdivision of eukaryotic chromosomes into discrete topological domains depends upon the pairing of boundary elements. In flies, pairing interactions are specific and typically orientation dependent. They occur in cis between neighboring heterologous boundaries, and in trans between homologous boundaries. One potential mechanism for ensuring pairing-interaction specificity is the use of sequence-specific DNA-binding proteins that can bind simultaneously with two or more recognition sequences. Our studies indicate that Insv can assemble into a multivalent DNA-binding complex and that the N-terminal Insv multimerization domain is critical for boundary function.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0249-2) contains supplementary material, which is available to authorized users.

4C-seq characterization of Drosophila BEAF binding regions provides evidence for highly variable long-distance interactions between active chromatin

Chromatin organization is crucial for nuclear functions such as gene regulation, DNA replication and DNA repair. Insulator binding proteins, such as the Drosophila Boundary Element-Associated Factor (BEAF), are involved in chromatin organization. To further understand the role of BEAF, we detected cis- and trans-interaction partners of four BEAF binding regions (viewpoints) using 4C (circular chromosome conformation capture) and analyzed their association with different genomic features. Previous genome-wide mapping found that BEAF usually binds near transcription start sites, often of housekeeping genes, so our viewpoints were selected to reflect this. Our 4C data show the interaction partners of our viewpoints are highly variable and generally enriched for active chromatin marks. The most consistent association was with housekeeping genes, a feature in common with our viewpoints. Fluorescence in situ hybridization indicated that the long-distance interactions occur even in the absence of BEAF. These data are most consistent with a model in which BEAF is redundant with other factors found at active promoters. Our results point to principles of long-distance interactions made by active chromatin, supporting a previously proposed model in which condensed chromatin is sticky and associates into topologically associating domains (TADs) separated by active chromatin. We propose that the highly variable long-distance interactions we detect are driven by redundant factors that open chromatin to promote transcription, combined with active chromatin filling spaces between TADs while packing of TADs relative to each other varies from cell to cell.

The BEN Domain Protein Insensitive Binds to the Fab-7 Chromatin Boundary To Establish Proper Segmental Identity in Drosophila

Boundaries (insulators) in the Drosophila bithorax complex (BX-C) delimit autonomous regulatory domains that orchestrate the parasegment (PS)-specific expression of the BX-C homeotic genes. The Fab-7 boundary separates the iab-6 and iab-7 regulatory domains, which control Abd-B expression in PS11 and PS12, respectively. This boundary is composed of multiple functionally redundant elements and has two key functions: it blocks cross talk between iab-6 and iab-7 and facilitates boundary bypass. Here, we show that two BEN domain protein complexes, Insensitive and Elba, bind to multiple sequences located in the Fab-7 nuclease hypersensitive regions. Two of these sequences are recognized by both Insv and Elba and correspond to a CCAATTGG palindrome. Elba also binds to a related CCAATAAG sequence, while Insv does not. However, the third Insv recognition sequences is ∼100 bp in length and contains the CCAATAAG sequence at one end. Both Insv and Elba are assembled into large complexes (∼420 and ∼265-290 kDa, respectively) in nuclear extracts. Using a sensitized genetic background, we show that the Insv protein is required for Fab-7 boundary function and that PS11 identity is not properly established in insv mutants. This is the first demonstration that a BEN domain protein is important for the functioning of an endogenous fly boundary.

Multiple interactions are involved in a highly specific association of the Mod(mdg4)-67.2 isoform with the Su(Hw) sites in Drosophila

The best-studied Drosophila insulator complex consists of two BTB-containing proteins, the Mod(mdg4)-67.2 isoform and CP190, which are recruited to the chromatin through interactions with the DNA-binding Su(Hw) protein. It was shown previously that Mod(mdg4)-67.2 is critical for the enhancer-blocking activity of the Su(Hw) insulators and it differs from more than 30 other Mod(mdg4) isoforms by the C-terminal domain required for a specific interaction with Su(Hw) only. The mechanism of the highly specific association between Mod(mdg4)-67.2 and Su(Hw) is not well understood. Therefore, we have performed a detailed analysis of domains involved in the interaction of Mod(mdg4)-67.2 with Su(Hw) and CP190. We found that the N-terminal region of Su(Hw) interacts with the glutamine-rich domain common to all the Mod(mdg4) isoforms. The unique C-terminal part of Mod(mdg4)-67.2 contains the Su(Hw)-interacting domain and the FLYWCH domain that facilitates a specific association between Mod(mdg4)-67.2 and the CP190/Su(Hw) complex. Finally, interaction between the BTB domain of Mod(mdg4)-67.2 and the M domain of CP190 has been demonstrated. By using transgenic lines expressing different protein variants, we have shown that all the newly identified interactions are to a greater or lesser extent redundant, which increases the reliability in the formation of the protein complexes.

DIVERSITY in binding, regulation, and evolution revealed from high-throughput ChIP

Genome-wide in vivo protein-DNA interactions are routinely mapped using high-throughput chromatin immunoprecipitation (ChIP). ChIP-reported regions are typically investigated for enriched sequence-motifs, which are likely to model the DNA-binding specificity of the profiled protein and/or of co-occurring proteins. However, simple enrichment analyses can miss insights into the binding-activity of the protein. Note that ChIP reports regions making direct contact with the protein as well as those binding through intermediaries. For example, consider a ChIP experiment targeting protein X, which binds DNA at its cognate sites, but simultaneously interacts with four other proteins. Each of these proteins also binds to its own specific cognate sites along distant parts of the genome, a scenario consistent with the current view of transcriptional hubs and chromatin loops. Since ChIP will pull down all X-associated regions, the final reported data will be a union of five distinct sets of regions, each containing binding sites of one of the five proteins, respectively. Characterizing all five different motifs and the corresponding sets is important to interpret the ChIP experiment and ultimately, the role of X in regulation. We present diversity which attempts exactly this: it partitions the data so that each partition can be characterized with its own de novo motif. Diversity uses a Bayesian approach to identify the optimal number of motifs and the associated partitions, which together explain the entire dataset. This is in contrast to standard motif finders, which report motifs individually enriched in the data, but do not necessarily explain all reported regions. We show that the different motifs and associated regions identified by diversity give insights into the various complexes that may be forming along the chromatin, something that has so far not been attempted from ChIP data. Webserver at http://diversity.ncl.res.in/; standalone (Mac OS X/Linux) from https://github.com/NarlikarLab/DIVERSITY/releases/tag/v1.0.0.

A high-throughput chromatin immunoprecipitation (ChIP) experiment identifies genomic regions bound by a protein in vivo. Current motif-discovery approaches seek an enriched motif signature in the reported regions, which they can attribute to the protein’s binding preferences. However, Diversity models the fact that since a ChIP experiment pulls down regions participating in all complexes involving the profiled protein, the reported regions are in all likelihood, a collection of different types of protein-DNA contacts. Diversity asks a different question: what sequence component caused a specific region to be reported in a ChIP experiment? The answer, in combination with additional data such as sequence conservation, SNPs, chromatin structure, downstream gene-expression, etc. can yield insights into the diverse regulatory mechanisms at play. The added benefits of a webserver and a standalone parallel version make diversity a practical tool for discovering new biology from ChIP experiments.

Role of Su(Hw) zinc finger 10 and interaction with CP190 and Mod(mdg4) proteins in recruiting the Su(Hw) complex to chromatin sites in Drosophila

Su(Hw) belongs to the class of proteins that organize chromosome architecture and boundaries/insulators between regulatory domains. This protein contains a cluster of 12 zinc finger domains most of which are responsible for binding to three different modules in the consensus site. Su(Hw) forms a complex with CP190 and Mod(mdg4)-67.2 proteins that binds to well-known Drosophila insulators. To understand how Su(Hw) performs its activities and binds to specific sites in chromatin, we have examined the previously described su(Hw)^f mutation that disrupts the 10th zinc finger (ZF10) responsible for Su(Hw) binding to the upstream module. The results have shown that Su(Hw)^f loses the ability to interact with CP190 in the absence of DNA. In contrast, complete deletion of ZF10 does not prevent the interaction between Su(Hw)^Δ10 and CP190. Having studied insulator complex formation in different mutant backgrounds, we conclude that both association with CP190 and Mod(mdg4)-67.2 partners and proper organization of DNA binding site are essential for the efficient recruitment of the Su(Hw) complex to chromatin insulators.

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 .

Opbp is a new architectural/insulator protein required for ribosomal gene expression

A special class of poorly characterized architectural proteins is required for chromatin topology and enhancer–promoter interactions. Here, we identify Opbp as a new Drosophila architectural protein, interacting with CP190 both in vivo and in vitro. Opbp binds to a very restrictive set of genomic regions, through a rare sequence specific motif. These sites are co-bound by CP190 in vivo, and generally located at bidirectional promoters of ribosomal protein genes. We show that Opbp is essential for viability, and loss of opbp function, or destruction of its motif, leads to reduced ribosomal protein gene expression, indicating a functional role in promoter activation. As characteristic of architectural/insulator proteins, the Opbp motif is sufficient for distance-dependent reporter gene activation and enhancer-blocking activity, suggesting an Opbp-mediated enhancer–promoter interaction. Rather than having a constitutive role, Opbp represents a new type of architectural protein with a very restricted, yet essential, function in regulation of housekeeping gene expression.

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.

Gene functioning and storage within a folded genome

In mammals, genomic DNA that is roughly 2 m long is folded to fit the size of the cell nucleus that has a diameter of about 10 μm. The folding of genomic DNA is mediated via assembly of DNA-protein complex, chromatin. In addition to the reduction of genomic DNA linear dimensions, the assembly of chromatin allows to discriminate and to mark active (transcribed) and repressed (non-transcribed) genes. Consequently, epigenetic regulation of gene expression occurs at the level of DNA packaging in chromatin. Taking into account the increasing attention of scientific community toward epigenetic systems of gene regulation, it is very important to understand how DNA folding in chromatin is related to gene activity. For many years the hierarchical model of DNA folding was the most popular. It was assumed that nucleosome fiber (10-nm fiber) is folded into 30-nm fiber and further on into chromatin loops attached to a nuclear/chromosome scaffold. Recent studies have demonstrated that there is much less regularity in chromatin folding within the cell nucleus. The very existence of 30-nm chromatin fibers in living cells was questioned. On the other hand, it was found that chromosomes are partitioned into self-interacting spatial domains that restrict the area of enhancers action. Thus, TADs can be considered as structural-functional domains of the chromosomes. Here we discuss the modern view of DNA packaging within the cell nucleus in relation to the regulation of gene expression. Special attention is paid to the possible mechanisms of the chromatin fiber self-assembly into TADs. We discuss the model postulating that partitioning of the chromosome into TADs is determined by the distribution of active and inactive chromatin segments along the chromosome.

This article was specially invited by the editors and represents work by leading researchers.

Architectural protein Pita cooperates with dCTCF in organization of functional boundaries in Bithorax complex

Boundaries in the Bithorax complex (BX-C) of Drosophila delimit autonomous regulatory domains that drive parasegment-specific expression of homeotic genes. BX-C boundaries have two crucial functions: they must block crosstalk between adjacent regulatory domains and at the same time facilitate boundary bypass. The C2H2 zinc-finger protein Pita binds to several BX-C boundaries, including Fab-7 and Mcp To study Pita functions, we have used a boundary replacement strategy by substituting modified DNAs for the Fab-7 boundary, which is located between the iab-6 and iab-7 regulatory domains. Multimerized Pita sites block iab-6↔iab-7 crosstalk but fail to support iab-6 regulation of Abd-B (bypass). In the case of Fab-7, we used a novel sensitized background to show that the two Pita-binding sites contribute to its boundary function. Although Mcp is from BX-C, it does not function appropriately when substituted for Fab-7: it blocks crosstalk but does not support bypass. Mutation of the Mcp Pita site disrupts blocking activity and also eliminates dCTCF binding. In contrast, mutation of the Mcp dCTCF site does not affect Pita binding, and this mutant boundary retains partial function.

Three-Dimensional Genome Organization and Function in Drosophila

Understanding how the metazoan genome is used during development and cell differentiation is one of the major challenges in the postgenomic era. Early studies in Drosophila suggested that three-dimensional (3D) chromosome organization plays important regulatory roles in this process and recent technological advances started to reveal connections at the molecular level. Here we will consider general features of the architectural organization of the Drosophila genome, providing historical perspective and insights from recent work. We will compare the linear and spatial segmentation of the fly genome and focus on the two key regulators of genome architecture: insulator components and Polycomb group proteins. With its unique set of genetic tools and a compact, well annotated genome, Drosophila is poised to remain a model system of choice for rapid progress in understanding principles of genome organization and to serve as a proving ground for development of 3D genome-engineering techniques.

C2H2 Zinc Finger Proteins: The Largest but Poorly Explored Family of Higher Eukaryotic Transcription Factors

The emergence of whole-genome assays has initiated numerous genome-wide studies of transcription factor localizations at genomic regulatory elements (enhancers, promoters, silencers, and insulators), as well as facilitated the uncovering of some of the key principles of chromosomal organization. However, the proteins involved in the formation and maintenance of the chromosomal architecture and the organization of regulatory domains remain insufficiently studied. This review attempts to collate the available data on the abundant but still poorly understood family of proteins with clusters of the C2H2 zinc finger domains. One of the best known proteins of this family is a well conserved protein known as CTCF, which plays a key role in the establishment of the chromosomal architecture in vertebrates. The distinctive features of C2H2 zinc finger proteins include strong and specific binding to a long and unique DNA recognition target sequence and rapid expansion within various animal taxa during evolution. The reviewed data support a proposed model according to which many of the C2H2 proteins have functions that are similar to those of the CTCF in the organization of the chromatin architecture.