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

Insulators play important roles in controlling gene activity and maintaining regulatory independence between neighbouring genes. In this article, we show that the enhancer-blocking activity of the insulator present within the LTR retrotransposon Idefix can be abolished if two copies of the region containing the insulator--specifically, the long terminal repeat (LTR)--are fused to the retrotransposon's 5' untranslated region (5' UTR). The presence of this combination of two [LTR-5' UTR] modules is a prerequisite for the loss of enhancer-blocking activity. We further show that the 5' UTR causes flanking genomic sequences to be displaced to the nuclear periphery, which is not observed when two insulators are present by themselves. This study thus provides a functional link between insulators and independent genomic modules, which may cooperate to allow the specific regulation of defined genomic loci via nuclear repositioning. It further illustrates the complexity of genomic regulation within a chromatic environment with multiple functional elements.

Drosophila insulators were the first DNA elements found to regulate gene expression by delimiting chromatin contacts. We still do not know how many of them exist and what impact they have on the Drosophila genome folding. Contrary to vertebrates, there is no evidence that fly insulators block cohesin-mediated chromatin loop extrusion. Therefore, their mechanism of action remains uncertain. To bridge these gaps, we mapped chromatin contacts in Drosophila cells lacking the key insulator proteins CTCF and Cp190. With this approach, we found hundreds of insulator elements. Their study indicates that Drosophila insulators play a minor role in the overall genome folding but affect chromatin contacts locally at many loci. Our observations argue that Cp190 promotes cobinding of other insulator proteins and that the model, where Drosophila insulators block chromatin contacts by forming loops, needs revision. Our insulator catalog provides an important resource to study mechanisms of genome folding.

Hox genes impart segment identity to body structures along the anterior-posterior axis and are crucial for the proper development of all organisms. Multiple regulatory elements, best defined in Drosophila melanogaster, ensure that Hox expression patterns follow the spatial and temporal colinearity reflected in their tight genomic organization. However, the precise mechanisms that regulate colinear patterns of Hox gene expression remain unclear, especially in higher vertebrates where it is not fully determined how the distinct activation domains of the tightly clustered Hox genes are defined independently of each other. Here, we report the identification of a large number of novel cis-elements at mammalian Hox clusters that can help in regulating their precise expression pattern.

Hypoxia inducible factor (HIF) up-regulates the transcription of a few hundred genes required for the adaptation to hypoxia. This restricted set of targets is in sharp contrast with the widespread distribution of the HIF binding motif throughout the genome. Here, we investigated the transcriptional response of GYS1 and RUVBL2 genes to hypoxia to understand the mechanisms that restrict HIF activity toward specific genes. GYS1 and RUVBL2 genes are encoded by opposite DNA strands and separated by a short intergenic region (~1 kb) that contains a functional hypoxia response element equidistant to both genes. However, hypoxia induced the expression of GYS1 gene only. Analysis of the transcriptional response of chimeric constructs derived from the intergenic region revealed an inhibitory sequence whose deletion allowed RUVBL2 induction by HIF. Enhancer blocking assays, performed in cell culture and transgenic zebrafish, confirmed the existence of an insulator element within this inhibitory region that could explain the differential regulation of GYS1 and RUVBL2 by hypoxia. Hence, in this model, the selective response to HIF is achieved with the aid of insulator elements. This is the first report suggesting a role for insulators in the regulation of differential gene expression in response to environmental signals.

In Drosophila, defects in asymmetric cell division often result in the formation of stem-cell-derived tumours. Here, we show that very similar terminal brain tumour phenotypes arise through a fundamentally different mechanism. We demonstrate that brain tumours in l(3)mbt mutants originate from overproliferation of neuroepithelial cells in the optic lobes caused by derepression of target genes in the Salvador-Warts-Hippo (SWH) pathway. We use ChIP-sequencing to identify L(3)mbt binding sites and show that L(3)mbt binds to chromatin insulator elements. Mutating l(3)mbt or inhibiting expression of the insulator protein gene mod(mdg4) results in upregulation of SWH pathway reporters. As l(3)mbt tumours are rescued by mutations in bantam or yorkie or by overexpression of Expanded, the deregulation of SWH pathway target genes is an essential step in brain tumour formation. Therefore, very different primary defects result in the formation of brain tumours, which behave quite similarly in their advanced stages.

Germline stem cells (GSCs) are the progenitor cells of the germline for the lifetime of an animal. In Drosophila, these cells reside in a cellular niche that is required for both their maintenance (self-renewal) and differentiation (asymmetric division resulting in a daughter cell that differs from the GSC). The stem cell-daughter cell transition is tightly regulated by a number of processes, including an array of proteins required for genome stability. The germline stem-cell maintenance factor Stonewall (Stwl) associates with heterochromatin, but its molecular function is poorly understood. We performed RNA-Seq on stwl mutant ovaries and found significant derepression of many transposon families but not heterochromatic genes. We also discovered inappropriate expression of multiple classes of genes. Most prominent are testis-enriched genes, including the male germline sex-determination switch Phf7, the differentiation factor bgcn, and a large testis-specific gene cluster on chromosome 2, all of which are upregulated or ectopically expressed in stwl mutant ovaries. Surprisingly, we also found that RNAi knockdown of stwl in somatic S2 cells results in ectopic expression of these testis genes. Using parallel ChIP-Seq and RNA-Seq experiments in S2 cells, we discovered that Stwl localizes upstream of transcription start sites and at heterochromatic sequences including repetitive sequences associated with telomeres. Stwl is also enriched at bgcn, suggesting that it directly regulates this essential differentiation factor. Finally, we identify Stwl binding motifs that are shared with known insulator binding proteins. We propose that Stwl affects gene regulation, including repression of male transcripts in the female germline, by binding insulators and establishing chromatin boundaries.

Lineage reporters of human embryonic stem cell (hESC) lines are useful for differentiation studies and drug screening. Previously, we created reporter lines driven by an elongation factor 1 alpha (EF1α) promoter at a chromosome 13q32.3 locus in the hESC line WA09 and an abnormal hESC line BG01V in a site-specific manner. Expression of reporters in these lines was maintained in long-term culture at undifferentiated state. However, when these cells were differentiated into specific lineages, reduction in reporter expression was observed, indicating transgene silencing. To develop an efficient and reliable genetic engineering strategy in hESCs, we used chromatin insulator elements to flank single-copy transgenes and integrated the combined expression constructs via PhiC31/R4 integrase-mediated recombination technology to the chromosome 13 locus precisely. Two copies of cHS4 double-insulator sequences were placed adjacent to both 5' and 3' of the promoter reporter constructs. The green fluorescent protein (GFP) gene was driven by EF1α or CMV early enhancer/chicken β actin (CAG) promoter. In the engineered hESC lines, for both insulated CAG-GFP and EF1α-GFP, constitutive expression at the chromosome 13 locus was maintained during prolonged culture and in directed differentiation assays toward diverse types of neurons, pancreatic endoderm, and mesodermal progeny. In particular, described here is the first normal hESC fluorescent reporter line that robustly expresses GFP in both the undifferentiated state and throughout dopaminergic lineage differentiation. The dual strategy of utilizing insulator sequences and integration at the constitutive chromosome 13 locus ensures appropriate transgene expression. This is a valuable tool for lineage development study, gain- and loss-of-function experiments, and human disease modeling using hESCs.

Many low-penetrance breast cancer susceptibility loci are found to be located in non-protein-coding regions, suggesting their involvement in gene expression regulation. We identified the human/rat-conserved breast cancer susceptibility locus MCS5A/Mcs5a. This locus has been shown to act in a non-mammary cell-autonomous fashion through the immune system. The resistant Mcs5a allele from the Wistar-Kyoto (WKy) rat strain consists of two non-protein-coding genetic elements that must be located on the same chromosome to elicit the phenotype. In this study, we show the presence of a conserved higher order chromatin structure in MCS5A/Mcs5a located in between the synthetically interacting genetic elements. The looped elements are shown to be bound by CTCF and cohesin. We identify the downregulation of Fbxo10 expression in T cells as a strong candidate mechanism through which the interacting genetic elements of the resistant Mcs5a allele modulate mammary carcinoma susceptibility. Finally, we show that the human MCS5A polymorphisms associated with breast cancer risk are located at both sides of the looped structure and functionally interact to downregulate transcriptional activity, similar to rat Mcs5a. We propose a mechanistic model for MCS5a/Mcs5a in which a CTCF-mediated insulator loop encompassing the TOMM5/Tomm5 gene, resides in between and brings into closer physical proximity the synthetically and functionally interacting resistant genetic variants.

The specificity of interactions between genomic regulatory elements and potential target genes is influenced by the binding of insulator proteins such as CTCF, which can act as potent enhancer blockers when interposed between an enhancer and a promoter in a reporter assay. But not all CTCF sites genome-wide function as insulator elements, depending on cellular and genomic context. To dissect the influence of genomic context on enhancer blocker activity, we integrated reporter constructs with promoter-only, promoter and enhancer, and enhancer blocker configurations at hundreds of thousands of genomic sites using the Sleeping Beauty transposase. Deconvolution of reporter activity by genomic position reveals distinct expression patterns subject to genomic context, including a compartment of enhancer blocker reporter integrations with robust expression. The high density of integration sites permits quantitative delineation of characteristic genomic context sensitivity profiles and their decomposition into sensitivity to both local and distant DNase I hypersensitive sites. Furthermore, using a single-cell expression approach to test the effect of integrated reporters for differential expression of nearby endogenous genes reveals that CTCF insulator elements do not completely abrogate reporter effects on endogenous gene expression. Collectively, our results lend new insight into genomic regulatory compartmentalization and its influence on the determinants of promoter-enhancer specificity.

Alteration of regulatory DNA elements or their binding proteins may have drastic consequences for morphological evolution. Chromatin insulators are one example of such proteins and play a fundamental role in organizing gene expression. While a single insulator protein, CTCF (CCCTC-binding factor), is known in vertebrates, Drosophila melanogaster utilizes six additional factors. We studied the evolution of these proteins and show here that-in contrast to the bilaterian-wide distribution of CTCF-all other D. melanogaster insulators are restricted to arthropods. The full set is present exclusively in the genus Drosophila whereas only two insulators, Su(Hw) and CTCF, existed at the base of the arthropod clade and all additional factors have been acquired successively at later stages. Secondary loss of factors in some lineages further led to the presence of different insulator subsets in arthropods. Thus, the evolution of insulator proteins within arthropods is an ongoing and dynamic process that reshapes and supplements the ancient CTCF-based system common to bilaterians. Expansion of insulator systems may therefore be a general strategy to increase an organism's gene regulatory repertoire and its potential for morphological plasticity.

Insulators are DNA sequences thought to be important for the establishment and maintenance of cell-type specific nuclear architecture. In Drosophila there are several classes of insulators that appear to have unique roles in gene expression. The mechanisms involved in determining and regulating the specific roles of these insulator classes are not understood. Here we report that DNA Topoisomerase II modulates the activity of the Su(Hw) insulator. Downregulation of Topo II by RNAi or mutations in the Top2 gene result in disruption of Su(Hw) insulator function. This effect is mediated by the Mod(mdg4)2.2 protein, which is a unique component of the Su(Hw) insulator complex. Co-immunoprecipitation and yeast two-hybrid experiments show that Topo II and Mod(mdg4)2.2 proteins directly interact. In addition, mutations in Top2 cause a slight decrease of Mod(mdg4)2.2 transcript but have a dramatic effect on Mod(mdg4)2.2 protein levels. In the presence of proteasome inhibitors, normal levels of Mod(mdg4)2.2 protein and its binding to polytene chromosomes are restored. Thus, Topo II is required to prevent Mod(mdg4)2.2 degradation and, consequently, to stabilize Su(Hw) insulator-mediated chromatin organization.

Hybrid incompatibility between Drosophila melanogaster and D. simulans is caused by a lethal interaction of the proteins encoded by the Hmr and Lhr genes. In D. melanogaster the loss of HMR results in mitotic defects, an increase in transcription of transposable elements and a deregulation of heterochromatic genes. To better understand the molecular mechanisms that mediate HMR's function, we measured genome-wide localization of HMR in D. melanogaster tissue culture cells by chromatin immunoprecipitation. Interestingly, we find HMR localizing to genomic insulator sites that can be classified into two groups. One group belongs to gypsy insulators and another one borders HP1a bound regions at active genes. The transcription of the latter group genes is strongly affected in larvae and ovaries of Hmr mutant flies. Our data suggest a novel link between HMR and insulator proteins, a finding that implicates a potential role for genome organization in the formation of species.

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.

The binding affinity of DNA-binding proteins such as transcription factors is mainly determined by the base composition of the corresponding binding site on the DNA strand. Most proteins do not bind only a single sequence, but rather a set of sequences, which may be modeled by a sequence motif. Algorithms for de novo motif discovery differ in their promoter models, learning approaches, and other aspects, but typically use the statistically simple position weight matrix model for the motif, which assumes statistical independence among all nucleotides. However, there is no clear justification for that assumption, leading to an ongoing debate about the importance of modeling dependencies between nucleotides within binding sites. In the past, modeling statistical dependencies within binding sites has been hampered by the problem of limited data. With the rise of high-throughput technologies such as ChIP-seq, this situation has now changed, making it possible to make use of statistical dependencies effectively. In this work, we investigate the presence of statistical dependencies in binding sites of the human enhancer-blocking insulator protein CTCF by using the recently developed model class of inhomogeneous parsimonious Markov models, which is capable of modeling complex dependencies while avoiding overfitting. These findings lead to a more detailed characterization of the CTCF binding motif, which is only poorly represented by independent nucleotide frequencies at several positions, predominantly at the 3' end.

The chromosomes of multicellular animals are organized into a series of topologically independent looped domains. This domain organization is critical for the proper utilization and propagation of the genetic information encoded by the chromosome. A special set of architectural elements, called boundaries or insulators, are responsible both for subdividing the chromatin into discrete domains and for determining the topological organization of these domains. Central to the architectural functions of insulators are homologous and heterologous insulator:insulator pairing interactions. The former (pairing between copies of the same insulator) dictates the process of homolog alignment and pairing in trans, while the latter (pairing between different insulators) defines the topology of looped domains in cis. To elucidate the principles governing these architectural functions, we use two insulators, Homie and Nhomie, that flank the Drosophila even skipped locus. We show that homologous insulator interactions in trans, between Homie on one homolog and Homie on the other, or between Nhomie on one homolog and Nhomie on the other, mediate transvection. Critically, these homologous insulator:insulator interactions are orientation-dependent. Consistent with a role in the alignment and pairing of homologs, self-pairing in trans is head-to-head. Head-to-head self-interactions in cis have been reported for other fly insulators, suggesting that this is a general principle of self-pairing. Homie and Nhomie not only pair with themselves, but with each other. Heterologous Homie-Nhomie interactions occur in cis, and we show that they serve to delimit a looped chromosomal domain that contains the even skipped transcription unit and its associated enhancers. The topology of this loop is defined by the heterologous pairing properties of Homie and Nhomie. Instead of being head-to-head, which would generate a circular loop, Homie-Nhomie pairing is head-to-tail. Head-to-tail pairing in cis generates a stem-loop, a configuration much like that observed in classical lampbrush chromosomes. These pairing principles provide a mechanistic underpinning for the observed topologies within and between chromosomes.

Chromatin insulator elements and associated proteins have been proposed to partition eukaryotic genomes into sets of independently regulated domains. Here we test this hypothesis by quantitative genome-wide analysis of insulator protein binding to Drosophila chromatin. We find distinct combinatorial binding of insulator proteins to different classes of sites and uncover a novel type of insulator element that binds CP190 but not any other known insulator proteins. Functional characterization of different classes of binding sites indicates that only a small fraction act as robust insulators in standard enhancer-blocking assays. We show that insulators restrict the spreading of the H3K27me3 mark but only at a small number of Polycomb target regions and only to prevent repressive histone methylation within adjacent genes that are already transcriptionally inactive. RNAi knockdown of insulator proteins in cultured cells does not lead to major alterations in genome expression. Taken together, these observations argue against the concept of a genome partitioned by specialized boundary elements and suggest that insulators are reserved for specific regulation of selected genes.

The organization of the genome into transcriptionally active and inactive chromatin domains requires well-delineated chromatin boundaries and insulator functions in order to maintain the identity of adjacent genomic loci with antagonistic chromatin marks and functionality. In plants that lack known chromatin insulators, the mechanisms that prevent heterochromatin spreading into euchromatin remain to be identified. Here, we show that DNA Topoisomerase VI participates in a chromatin boundary function that safeguards the expression of genes in euchromatin islands within silenced heterochromatin regions. While some transposable elements are reactivated in mutants of the Topoisomerase VI complex, genes insulated in euchromatin islands within heterochromatic regions of the Arabidopsis thaliana genome are specifically down-regulated. H3K9me2 levels consistently increase at euchromatin island loci and decrease at some transposable element loci. We further show that Topoisomerase VI physically interacts with S-adenosylmethionine synthase methionine adenosyl transferase 3 (MAT3), which is required for H3K9me2. A Topoisomerase VI defect affects MAT3 occupancy on heterochromatic elements and its exclusion from euchromatic islands, thereby providing a possible mechanistic explanation to the essential role of Topoisomerase VI in the delimitation of chromatin domains.

Chromatin insulators, or boundary elements, affect promoter-enhancer interactions and buffer transgenes from position effects. The gypsy insulator of Drosophila is bound by a protein complex with two characterized components, the zinc finger protein Suppressor of Hairy-wing [Su(Hw)] and Mod(mdg4)2.2, which is one of the multiple spliced variants encoded by the modifier of mdg4 [mod(mdg4)] gene. A genetic screen for dominant enhancers of the mod(mdg4) phenotype identified the Centrosomal Protein 190 (CP190) as an essential constituent of the gypsy insulator. The function of the centrosome is not affected in CP190 mutants whereas gypsy insulator activity is impaired. CP190 associates physically with both Su(Hw) and Mod(mdg4)2.2 and colocalizes with both proteins on polytene chromosomes. CP190 does not interact directly with insulator sequences present in the gypsy retrotransposon but binds to a previously characterized endogenous insulator, and it is necessary for the formation of insulator bodies. The results suggest that endogenous gypsy insulators contain binding sites for CP190, which is essential for insulator function, and may or may not contain binding sites for Su(Hw) and Mod(mdg4)2.2.

The three-dimensional (3D) organization of the genome is highly dynamic, changing during development and varying across different tissues and cell types. Recent studies indicate that these changes alter regulatory interactions, leading to changes in gene expression. Despite its importance, the mechanisms that influence genomic organization remain poorly understood. We have previously identified a network of chromatin boundary elements, or insulators, in the Drosophila Antennapedia homeotic complex (ANT-C). These genomic elements interact with one another to tether chromatin loops that could block or promote enhancer-promoter interactions. To understand the function of these insulators, we assessed their interactions by measuring their 3D nuclear distance in developing animal tissues. Our data suggest that the ANT-C Hox complex might be in a folded or looped configuration rather than in a random or extended form. The architecture of the ANT-C complex, as read out by the pair-wise distance between insulators, undergoes a strong compression during late embryogenesis, coinciding with the reduction of cell and nuclear diameters due to continued cell divisions in post-cleavage cells. Our results suggest that genomic architecture and gene regulation may be influenced by cellular morphology and movement during development.