A position-effect assay for boundaries of higher order chromosomal domains

Enhancer selectivity in space and time: from enhancer-promoter interactions to promoter activation

The primary regulators of metazoan gene expression are enhancers, originally functionally defined as DNA sequences that can activate transcription at promoters in an orientation-independent and distance-independent manner. Despite being crucial for gene regulation in animals, what mechanisms underlie enhancer selectivity for promoters, and more fundamentally, how enhancers interact with promoters and activate transcription, remain poorly understood. In this Review, we first discuss current models of enhancer-promoter interactions in space and time and how enhancers affect transcription activation. Next, we discuss different mechanisms that mediate enhancer selectivity, including repression, biochemical compatibility and regulation of 3D genome structure. Through 3D polymer simulations, we illustrate how the ability of 3D genome folding mechanisms to mediate enhancer selectivity strongly varies for different enhancer-promoter interaction mechanisms. Finally, we discuss how recent technical advances may provide new insights into mechanisms of enhancer-promoter interactions and how technical biases in methods such as Hi-C and Micro-C and imaging techniques may affect their interpretation.

Developmental and Housekeeping Genes: Two Types of Genetic Organization in the Drosophila Genome

We developed a procedure for locating genes on Drosophila melanogaster polytene chromosomes and described three types of chromosome structures (gray bands, black bands, and interbands), which differed markedly in morphological and genetic properties. This was reached through the use of our original methods of molecular and genetic analysis, electron microscopy, and bioinformatics data processing. Analysis of the genome-wide distribution of these properties led us to a bioinformatics model of the Drosophila genome organization, in which the genome was divided into two groups of genes. One was constituted by 65, in which the genome was divided into two groups, 62 genes that are expressed in most cell types during life cycle and perform basic cellular functions (the so-called "housekeeping genes"). The other one was made up of 3162 genes that are expressed only at particular stages of development ("developmental genes"). These two groups of genes are so different that we may state that the genome has two types of genetic organization. Different are the timings of their expression, chromatin packaging levels, the composition of activating and deactivating proteins, the sizes of these genes, the lengths of their introns, the organization of the promoter regions of the genes, the locations of origin recognition complexes (ORCs), and DNA replication timings.

Аpplication of massive parallel reporter analysis in biotechnology and medicine

The development and functioning of an organism relies on tissue-specific gene programs. Genome regulatory elements play a key role in the regulation of such programs, and disruptions in their function can lead to the development of various pathologies, including cancers, malformations and autoimmune diseases. The emergence of high-throughput genomic studies has led to massively parallel reporter analysis (MPRA) methods, which allow the functional verification and identification of regulatory elements on a genome-wide scale. Initially MPRA was used as a tool to investigate fundamental aspects of epigenetics, but the approach also has great potential for clinical and practical biotechnology. Currently, MPRA is used for validation of clinically significant mutations, identification of tissue-specific regulatory elements, search for the most promising loci for transgene integration, and is an indispensable tool for creating highly efficient expression systems, the range of application of which extends from approaches for protein development and design of next-generation therapeutic antibody superproducers to gene therapy. In this review, the main principles and areas of practical application of high-throughput reporter assays will be discussed.

The spatial organization of transcriptional control

In animals, the sequences for controlling gene expression do not concentrate just at the transcription start site of genes, but are frequently thousands to millions of base pairs distal to it. The interaction of these sequences with one another and their transcription start sites is regulated by factors that shape the three-dimensional (3D) organization of the genome within the nucleus. Over the past decade, indirect tools exploiting high-throughput DNA sequencing have helped to map this 3D organization, have identified multiple key regulators of its structure and, in the process, have substantially reshaped our view of how 3D genome architecture regulates transcription. Now, new tools for high-throughput super-resolution imaging of chromatin have directly visualized the 3D chromatin organization, settling some debates left unresolved by earlier indirect methods, challenging some earlier models of regulatory specificity and creating hypotheses about the role of chromatin structure in transcriptional regulation.

No apparent role for the Wari insulator in transcriptional regulation of the endogenous white gene of Drosophila melanogaster

Chromatin insulators have been proposed to play an important role in chromosome organization and local regulatory interactions. In Drosophila , one of these insulators is known as Wari. It is located immediately downstream of the 3' end of the white transcription unit. Wari has been proposed to interact with the white promoter region, thereby facilitating recycling of the RNA polymerase machinery. We have tested this model by deleting the Wari insulator at the endogenous white locus and could not detect a significant effect on eye pigmentation.

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

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

Establishment and inheritance of epigenetic transcriptional memory

For certain inducible genes, the rate and molecular mechanism of transcriptional activation depends on the prior experiences of the cell. This phenomenon, called epigenetic transcriptional memory, accelerates reactivation and requires both changes in chromatin structure and recruitment of poised RNA Polymerase II (RNAPII) to the promoter. Forms of epigenetic transcriptional memory have been identified in S. cerevisiae, D. melanogaster, C. elegans, and mammals. A well-characterized model of memory is found in budding yeast where memory of inositol starvation involves a positive feedback loop between gene-and condition-specific transcription factors, which mediate an interaction with the nuclear pore complex and a characteristic histone modification: histone H3 lysine 4 dimethylation (H3K4me2). This histone modification permits recruitment of a memory-specific pre-initiation complex, poising RNAPII at the promoter. During memory, H3K4me2 is essential for recruitment of RNAPII and faster reactivation, but RNAPII is not required for H3K4me2. Unlike the RNAPII-dependent H3K4me2 associated with active transcription, RNAPII-independent H3K4me2 requires Nup100, SET3C, the Leo1 subunit of the Paf1 complex and can be inherited through multiple cell cycles upon disrupting the interaction with the Nuclear Pore Complex. The H3K4 methyltransferase (COMPASS) physically interacts with the potential reader (SET3C), suggesting a molecular mechanism for the spreading and re-incorporation of H3K4me2 following DNA replication. Thus, epigenetic transcriptional memory is a conserved adaptation that utilizes a heritable chromatin state, allowing cells and organisms to alter their gene expression programs in response to recent experiences over intermediate time scales.

Single-nucleus RNA-sequencing in pre-cellularization Drosophila melanogaster embryos

Our current understanding of the regulation of gene expression in the early Drosophila melanogaster embryo comes from observations of a few genes at a time, as with in situ hybridizations, or observation of gene expression levels without regards to patterning, as with RNA-sequencing. Single-nucleus RNA-sequencing however, has the potential to provide new insights into the regulation of gene expression for many genes at once while simultaneously retaining information regarding the position of each nucleus prior to dissociation based on patterned gene expression. In order to establish the use of single-nucleus RNA sequencing in Drosophila embryos prior to cellularization, here we look at gene expression in control and insulator protein, dCTCF, maternal null embryos during zygotic genome activation at nuclear cycle 14. We find that early embryonic nuclei can be grouped into distinct clusters according to gene expression. From both virtual and published in situ hybridizations, we also find that these clusters correspond to spatial regions of the embryo. Lastly, we provide a resource of candidate differentially expressed genes that might show local changes in gene expression between control and maternal dCTCF null nuclei with no detectable differential expression in bulk. These results highlight the potential for single-nucleus RNA-sequencing to reveal new insights into the regulation of gene expression in the early Drosophila melanogaster embryo.

TFIIIC-based chromatin insulators through eukaryotic evolution

Eukaryotic chromosomes are divided into domains with distinct structural and functional properties, such as differing levels of chromatin compaction and gene transcription. Domains of relatively compact chromatin and minimal transcription are termed heterochromatic, whereas euchromatin is more open and actively transcribed. Insulators separate these domains and maintain their distinct features. Disruption of insulators can cause diseases such as cancer. Many insulators contain tRNA genes (tDNAs), examples of which have been shown to block the spread of activating or silencing activities. This characteristic of specific tDNAs is conserved through evolution, such that human tDNAs can serve as barriers to the spread of silencing in fission yeast. Here we demonstrate that tDNAs from the methylotrophic fungus Pichia pastoris can function effectively as insulators in distantly-related budding yeast. Key to the function of tDNAs as insulators is TFIIIC, a transcription factor that is also required for their expression. TFIIIC binds additional loci besides tDNAs, some of which have insulator activity. Although the mechanistic basis of TFIIIC-based insulation has been studied extensively in yeast, it is largely uncharacterized in metazoa. Utilising publicly-available genome-wide ChIP-seq data, we consider the extent to which mechanisms conserved from yeast to man may suffice to allow efficient insulation by TFIIIC in the more challenging chromatin environments of metazoa and suggest features that may have been acquired during evolution to cope with new challenges. We demonstrate the widespread presence at human tDNAs of USF1, a transcription factor with well-established barrier activity in vertebrates. We predict that tDNA-based insulators in higher organisms have evolved through incorporation of modules, such as binding sites for factors like USF1 and CTCF that are absent from yeasts, thereby strengthening function and providing opportunities for regulation between cell types.

Pre-existing chromatin accessibility of switchable repressive compartment delineates cell plasticity

Cell plasticity endows differentiated cells with competence to be reprogrammed to other lineages. Although extrinsic factors driving cell-identity conversion have been extensively characterized, it remains elusive which intrinsic epigenetic attributes, including high-order chromatin organization, delineate cell plasticity. By analysing the transcription-factor-induced transdifferentiation from fibroblasts to hepatocytes, we uncovered contiguous compartment-switchable regions (CSRs) as a unique chromatin unit. Specifically, compartment B-to-A CSRs, enriched with hepatic genes, possessed a mosaic status of inactive chromatin and pre-existing and continuous accessibility in fibroblasts. Pre-existing accessibility enhanced the binding of inducible factor Foxa3, which triggered epigenetic activation and chromatin interaction as well as hepatic gene expression. Notably, these changes were restrained within B-to-A CSR boundaries that were defined by CTCF occupancy. Moreover, such chromatin organization and mosaic status were detectable in different cell types and involved in multiple reprogramming processes, suggesting an intrinsic chromatin attribute in understanding cell plasticity.

Methods of massive parallel reporter assays for investigation of enhancers

The correct deployment of genetic programs for development and differentiation relies on finely coordinated regulation of specific gene sets. Genomic regulatory elements play an exceptional role in this process. There are few types of gene regulatory elements, including promoters, enhancers, insulators and silencers. Alterations of gene regulatory elements may cause various pathologies, including cancer, congenital disorders and autoimmune diseases. The development of high-throughput genomic assays has made it possible to significantly accelerate the accumulation of information about the characteristic epigenetic properties of regulatory elements. In combination with high-throughput studies focused on the genome-wide distribution of epigenetic marks, regulatory proteins and the spatial structure of chromatin, this significantly expands the understanding of the principles of epigenetic regulation of genes and allows potential regulatory elements to be searched for in silico. However, common experimental approaches used to study the local characteristics of chromatin have a number of technical limitations that may reduce the reliability of computational identification of genomic regulatory sequences. Taking into account the variability of the functions of epigenetic determinants and complex multicomponent regulation of genomic elements activity, their functional verification is often required. A plethora of methods have been developed to study the functional role of regulatory elements on the genome scale. Common experimental approaches for in silico identification of regulatory elements and their inherent technical limitations will be described. The present review is focused on original high-throughput methods of enhancer activity reporter analysis that are currently used to validate predicted regulatory elements and to perform de novo searches. The methods described allow assessing the functional role of the nucleotide sequence of a regulatory element, to determine its exact boundaries and to assess the influence of the local state of chromatin on the activity of enhancers and gene expression. These approaches have contributed substantially to the understanding of the fundamental principles of gene regulation.

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

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

Expression of the foraging gene in adult Drosophila melanogaster

The foraging gene in Drosophila melanogaster, which encodes a cGMP-dependent protein kinase, is a highly conserved, complex gene with multiple pleiotropic behavioral and physiological functions in both the larval and adult fly. Adult foraging expression is less well characterized than in the larva. We characterized foraging expression in the brain, gastric system, and reproductive systems using a T2A-Gal4 gene-trap allele. In the brain, foraging expression appears to be restricted to multiple sub-types of glia. This glial-specific cellular localization of foraging was supported by single-cell transcriptomic atlases of the adult brain. foraging is extensively expressed in most cell types in the gastric and reproductive systems. We then mapped multiple cis-regulatory elements responsible for parts of the observed expression patterns by a nested cloned promoter-Gal4 analysis. The mapped cis-regulatory elements were consistently modular when comparing the larval and adult expression patterns. These new data using the T2A-Gal4 gene-trap and cloned foraging promoter fusion GAL4's are discussed with respect to previous work using an anti-FOR antibody, which we show here to be non-specific. Future studies of foraging's function will consider roles for glial subtypes and peripheral tissues (gastric and reproductive systems) in foraging's pleiotropic behavioral and physiological effects.

3D genome organisation in Drosophila

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

An insulator blocks access to enhancers by an illegitimate promoter, preventing repression by transcriptional interference

Several distinct activities and functions have been described for chromatin insulators, which separate genes along chromosomes into functional units. Here, we describe a novel mechanism of functional separation whereby an insulator prevents gene repression. When the homie insulator is deleted from the end of a Drosophila even skipped (eve) locus, a flanking P-element promoter is activated in a partial eve pattern, causing expression driven by enhancers in the 3’ region to be repressed. The mechanism involves transcriptional read-through from the flanking promoter. This conclusion is based on the following. Read-through driven by a heterologous enhancer is sufficient to repress, even when homie is in place. Furthermore, when the flanking promoter is turned around, repression is minimal. Transcriptional read-through that does not produce anti-sense RNA can still repress expression, ruling out RNAi as the mechanism in this case. Thus, transcriptional interference, caused by enhancer capture and read-through when the insulator is removed, represses eve promoter-driven expression. We also show that enhancer-promoter specificity and processivity of transcription can have decisive effects on the consequences of insulator removal. First, a core heat shock 70 promoter that is not activated well by eve enhancers did not cause read-through sufficient to repress the eve promoter. Second, these transcripts are less processive than those initiated at the P-promoter, measured by how far they extend through the eve locus, and so are less disruptive. These results highlight the importance of considering transcriptional read-through when assessing the effects of insulators on gene expression.

Several distinct activities and functions have been described for chromatin insulators, which are regulatory DNA elements that separate genes along chromosomes into functional units. Here, we describe how insulators can prevent repression of one gene by preventing inappropriate transcription of another gene, without blocking read-through of transcription per se. When the insulator homie is deleted from the end of a transgenic eve locus, a flanking transposable element promoter is activated by eve enhancers, causing repression of the eve promoter. The mechanism involves transcriptional read-through from the flanking promoter, which disrupts normal eve enhancer-promoter activities. When the flanking promoter is turned around, repression of eve is minimal. Thus, transcriptional interference, caused by enhancer capture and read-through when the insulator is removed, represses the eve promoter. These results show a novel role for transcriptional read-through in the effects of insulators on gene expression.

Insulator foci distance correlates with cellular and nuclear morphology in early Drosophila embryos

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.

The relationship between genome structure and function

Precise patterns of gene expression in metazoans are controlled by three classes of regulatory elements: promoters, enhancers and boundary elements. During differentiation and development, these elements form specific interactions in dynamic higher-order chromatin structures. However, the relationship between genome structure and its function in gene regulation is not completely understood. Here we review recent progress in this field and discuss whether genome structure plays an instructive role in regulating gene expression or is a reflection of the activity of the regulatory elements of the genome.

Mutations in the insulator protein Suppressor of Hairy wing induce genome instability

Insulator proteins orchestrate the three-dimensional organization of the genome. Insulators function by facilitating communications between regulatory sequences and gene promoters, allowing accurate gene transcription regulation during embryo development and cell differentiation. However, the role of insulator proteins beyond genome organization and transcription regulation remains unclear. Suppressor of Hairy wing [Su(Hw)] is a Drosophila insulator protein that plays an important function in female oogenesis. Here we find that su(Hw) has an unsuspected role in genome stability during cell differentiation. We show that su(Hw) mutant developing egg chambers have poorly formed microtubule organization centers (MTOCs) in the germarium and display mislocalization of the anterior/posterior axis specification factor gurken in later oogenesis stages. Additionally, eggshells from partially rescued su(Hw) mutant female germline exhibit dorsoventral patterning defects. These phenotypes are very similar to phenotypes found in the important class of spindle mutants or in piRNA pathway mutants in Drosophila, in which defects generally result from the failure of germ cells to repair DNA damage. Similarities between mutations in su(Hw) and spindle and piRNA mutants are further supported by an excess of DNA damage in nurse cells, and because Gurken localization defects are partially rescued by mutations in the ATR (mei-41) and Chk1 (grapes) DNA damage response genes. Finally, we also show that su(Hw) mutants produce an elevated number of chromosome breaks in dividing neuroblasts from larval brains. Together, these findings suggest that Su(Hw) is necessary for the maintenance of genome integrity during Drosophila development, in both germline and dividing somatic cells.

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.

Transcriptional Enhancers in Drosophila

Key discoveries in Drosophila have shaped our understanding of cellular "enhancers." With a special focus on the fly, this chapter surveys properties of these adaptable cis-regulatory elements, whose actions are critical for the complex spatial/temporal transcriptional regulation of gene expression in metazoa. The powerful combination of genetics, molecular biology, and genomics available in Drosophila has provided an arena in which the developmental role of enhancers can be explored. Enhancers are characterized by diverse low- or high-throughput assays, which are challenging to interpret, as not all of these methods of identifying enhancers produce concordant results. As a model metazoan, the fly offers important advantages to comprehensive analysis of the central functions that enhancers play in gene expression, and their critical role in mediating the production of phenotypes from genotype and environmental inputs. A major challenge moving forward will be obtaining a quantitative understanding of how these cis-regulatory elements operate in development and disease.

Chromatin-modifying elements for recombinant protein production in mammalian cell systems

Mammalian cells are the preferred choice system for the production of complex molecules, such as recombinant therapeutic proteins. Although the technology for increasing the yield of proteins has improved rapidly, the process of selecting, identifying as well as maintaining high-yield cell clones is still troublesome, time-consuming and usually uncertain. Optimization of expression vectors is one of the most effective methods for enhancing protein expression levels. Several commonly used chromatin-modifying elements, including the matrix attachment region, ubiquitous chromatin opening elements, insulators, stabilizing anti-repressor elements can be used to increase the expression level and stability of recombinant proteins. In this review, these chromatin-modifying elements used for the expression vector optimization in mammalian cells are summarized, and future strategies for the utilization of expression cassettes are also discussed.

Insights into HP1a-Chromatin Interactions

Understanding the packaging of DNA into chromatin has become a crucial aspect in the study of gene regulatory mechanisms. Heterochromatin establishment and maintenance dynamics have emerged as some of the main features involved in genome stability, cellular development, and diseases. The most extensively studied heterochromatin protein is HP1a. This protein has two main domains, namely the chromoshadow and the chromodomain, separated by a hinge region. Over the years, several works have taken on the task of identifying HP1a partners using different strategies. In this review, we focus on describing these interactions and the possible complexes and subcomplexes associated with this critical protein. Characterization of these complexes will help us to clearly understand the implications of the interactions of HP1a in heterochromatin maintenance, heterochromatin dynamics, and heterochromatin's direct relationship to gene regulation and chromatin organization.

Insulator-based loops mediate the spreading of H3K27me3 over distant micro-domains repressing euchromatin genes

Background

Chromosomes are subdivided spatially to delimit long-range interactions into topologically associating domains (TADs). TADs are often flanked by chromatin insulators and transcription units that may participate in such demarcation. Remarkably, single-cell Drosophila TAD units correspond to dynamic heterochromatin nano-compartments that can self-assemble. The influence of insulators on such dynamic compartmentalization remains unclear. Moreover, to what extent heterochromatin domains are fully compartmentalized away from active genes remains unclear from Drosophila to human.

Results

Here, we identify H3K27me3 micro-domains genome-wide in Drosophila, which are attributed to the three-dimensional spreading of heterochromatin marks into euchromatin. Whereas depletion of insulator proteins increases H3K27me3 spreading locally, across heterochromatin borders, it concomitantly decreases H3K27me3 levels at distant micro-domains discrete sites. Quantifying long-range interactions suggests that random interactions between heterochromatin TADs and neighbor euchromatin cannot predict the presence of micro-domains, arguing against the hypothesis that they reflect defects in self-folding or in insulating repressive TADs. Rather, micro-domains are predicted by specific long-range interactions with the TAD borders bound by insulator proteins and co-factors required for looping. Accordingly, H3K27me3 spreading to distant sites is impaired by insulator mutants that compromise recruitment of looping co-factors. Both depletions and insulator mutants significantly reduce H3K27me3 micro-domains, deregulating the flanking genes.

Conclusions

Our data highlight a new regulatory mode of H3K27me3 by insulator-based long-range interactions controlling distant euchromatic genes.

Overlapping but Distinct Sequences Play Roles in the Insulator and Promoter Activities of the Drosophila BEAF-Dependent scs' Insulator

Chromatin domain insulators are thought to help partition the genome into genetic units called topologically associating domains (TADs). In Drosophila, TADs are often separated by inter-TAD regions containing active housekeeping genes and associated insulator binding proteins. This raises the question of whether insulator binding proteins are involved primarily in chromosomal TAD architecture or gene activation, or if these two activities are linked. The Boundary Element-Associated Factor of 32 kDa (BEAF-32, or BEAF for short) is usually found in inter-TADs. BEAF was discovered based on binding to the scs' insulator, and is important for the insulator activity of scs' and other BEAF binding sites. There are divergent promoters in scs' with a BEAF binding site by each. Here, we dissect the scs' insulator to identify DNA sequences important for insulator and promoter activity, focusing on the half of scs' with a high affinity BEAF binding site. We find that the BEAF binding site is important for both insulator and promoter activity, as is another sequence we refer to as LS4. Aside from that, different sequences play roles in insulator and promoter activity. So while there is overlap and BEAF is important for both, insulator and promoter activity can be separated.

Promoter-Proximal Chromatin Domain Insulator Protein BEAF Mediates Local and Long-Range Communication with a Transcription Factor and Directly Activates a Housekeeping Promoter in Drosophila

BEAF (Boundary Element-Associated Factor) was originally identified as a Drosophila melanogaster chromatin domain insulator-binding protein, suggesting a role in gene regulation through chromatin organization and dynamics. Genome-wide mapping found that BEAF usually binds near transcription start sites, often of housekeeping genes, suggesting a role in promoter function. This would be a nontraditional role for an insulator-binding protein. To gain insight into molecular mechanisms of BEAF function, we identified interacting proteins using yeast two-hybrid assays. Here, we focus on the transcription factor Serendipity δ (Sry-δ). Interactions were confirmed in pull-down experiments using bacterially expressed proteins, by bimolecular fluorescence complementation, and in a genetic assay in transgenic flies. Sry-δ interacted with promoter-proximal BEAF both when bound to DNA adjacent to BEAF or > 2-kb upstream to activate a reporter gene in transient transfection experiments. The interaction between BEAF and Sry-δ was detected using both a minimal developmental promoter (y) and a housekeeping promoter (RpS12), while BEAF alone strongly activated the housekeeping promoter. These two functions for BEAF implicate it in playing a direct role in gene regulation at hundreds of BEAF-associated promoters.

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.

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 Role of Insulators in Transgene Transvection in Drosophila

Transvection is the phenomenon where a transcriptional enhancer activates a promoter located on the homologous chromosome. It has been amply documented in Drosophila where homologs are closely paired in most, if not all, somatic nuclei, but it has been known to rarely occur in mammals as well. We have taken advantage of site-directed transgenesis to insert reporter constructs into the same genetic locus in Drosophila and have evaluated their ability to engage in transvection by testing many heterozygous combinations. We find that transvection requires the presence of an insulator element on both homologs. Homotypic trans-interactions between four different insulators can support transvection: the gypsy insulator (GI), Wari, Fab-8 and 1A2; GI and Fab-8 are more effective than Wari or 1A2 We show that, in the presence of insulators, transvection displays the characteristics that have been previously described: it requires homolog pairing, but can happen at any of several loci in the genome; a solitary enhancer confronted with an enhancerless reporter is sufficient to drive transcription; it is weaker than the action of the same enhancer-promoter pair in cis, and it is further suppressed by cis-promoter competition. Though necessary, the presence of homotypic insulators is not sufficient for transvection; their position, number and orientation matters. A single GI adjacent to both enhancer and promoter is the optimal configuration. The identity of enhancers and promoters in the vicinity of a trans-interacting insulator pair is also important, indicative of complex insulator-enhancer-promoter interactions.

Using a phiC31 "Disintegrase" to make new attP sites in the Drosophila genome at locations showing chromosomal position effects

An engineered phiC31 “Disintegrase” able to make an attP site in Drosophila out of an attR-attL pair is described. This was used to generate attP sites at genomic locations where a mini-white (mini-w) transgene was subject to chromosomal position effects (CPE). The first step was random genomic integration of a P-element-based transposon with an insulated mini-w transgene. We then removed the upstream insulator using FLP recombinase to detect CPE. Next mini-w and the downstream insulator were “dis-integrated” leaving behind an attP site. The location is marked by a yellow⁺ transgene that is flanked by loxP sites, so it can also be removed. Using this system, we generated 10 new attP landing platforms. Three of these showing strong activating CPE were selected for further analysis. We show that the attP sites are functional by integrating in plasmids with attB sites. The CPE is recapitulated and can be blocked by insulators. We show that a dimerized 215 bp fragment of the 500 bp BEAF-dependent scs’ insulator containing a high affinity BEAF binding site blocks the CPE, while a monomer of the sequence is less effective. This indicates that two BEAF binding sites make a stronger insulator than a single site. This system could be useful for generating attP sites at prescreened sites for other purposes, such as studying CPE in embryos or other tissues or for use with “trapped” enhancers of interest.

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.

New Insights Into the Long Non-coding RNA SRA: Physiological Functions and Mechanisms of Action

Long non-coding RNAs (lncRNA) are emerging as new genetic/epigenetic regulators that can impact almost all physiological functions. Here, we focus on the long non-coding steroid receptor RNA activator (SRA), including new insights into its effects on gene expression, the cell cycle, and differentiation; how these relate to physiology and disease; and the mechanisms underlying these effects. We discuss how SRA acts as an RNA coactivator in nuclear receptor signaling; its effects on steroidogenesis, adipogenesis, and myocyte differentiation; the impact on breast and prostate cancer tumorigenesis; and, finally, its ability to modulate hepatic steatosis through several signaling pathways. Genome-wide analysis reveals that SRA regulates hundreds of target genes in adipocytes and breast cancer cells and binds to thousands of genomic sites in human pluripotent stem cells. Recent studies indicate that SRA acts as a molecular scaffold and forms networks with numerous coregulators and chromatin-modifying regulators in both activating and repressive complexes. We discuss how modifications to SRA's unique stem-loop secondary structure are important for SRA function, and highlight the various SRA isoforms and mutations that have clinical implications. Finally, we discuss the future directions for better understanding the molecular mechanisms of SRA action and how this might lead to new diagnostic and therapeutic approaches.

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.

Chromatin Architectural Changes during Cellular Senescence and Aging

Chromatin 3D structure is highly dynamic and associated with many biological processes, such as cell cycle progression, cellular differentiation, cell fate reprogramming, cancer development, cellular senescence, and aging. Recently, by using chromosome conformation capture technologies, tremendous findings have been reported about the dynamics of genome architecture, their associated proteins, and the underlying mechanisms involved in regulating chromatin spatial organization and gene expression. Cellular senescence and aging, which involve multiple cellular and molecular functional declines, also undergo significant chromatin structural changes, including alternations of heterochromatin and disruption of higher-order chromatin structure. In this review, we summarize recent findings related to genome architecture, factors regulating chromatin spatial organization, and how they change during cellular senescence and aging.

Identification and characterization of cis-regulatory elements 'insulator and repressor' in PPARD gene

AIM

Identification and functional characterization of cis-regulatory elements in human PPARD gene.

METHODS

We used various bioinformatic tools on the publicly available human genome and Encyclopedia of DNA Elements databases to explore potential cis-regulatory elements in PPARD gene region.

RESULTS

We predicted an insulator and an enhancer element in intron 2 of PPARD gene. Functional characterization using transient transfection, reporter assay and CTCF binding confirmed the insulator status. However, the predicted enhancer element showed repressor/silencer activity. Finally, we observed a potential interaction between these two cis-regulatory elements which is in agreement with 5C-Encyclopedia of DNA Elements data.

CONCLUSION

We report two functionally validated cis-regulatory elements in PPARD gene which will aid in understanding its regulation and role in metabolic functions.

Role of CTCF in DNA damage response

CCCTC-binding factor (CTCF) is a highly conserved, ubiquitously expressed zinc finger protein. CTCF is a multifunctional protein, associated with a number of vital cellular processes such as transcriptional activation, repression, insulation, imprinting and genome organization. Emerging evidence indicates that CTCF is also involved in DNA damage response. In this review, we focus on the newly identified role of CTCF in facilitating DNA double-strand break repair. Due to the large number of cellular processes in which CTCF is involved, factors that functionally affect CTCF could have serious implications on genomic stability. It is becoming increasingly clear that exposure to environmental toxicants could have adverse effects on CTCF functions. Here we discuss the various ways that environmental toxicants could impact CTCF functions and the potential consequences on DNA damage response.

Convergence of topological domain boundaries, insulators, and polytene interbands revealed by high-resolution mapping of chromatin contacts in the early Drosophila melanogaster embryo

High-throughput assays of three-dimensional interactions of chromosomes have shed considerable light on the structure of animal chromatin. Despite this progress, the precise physical nature of observed structures and the forces that govern their establishment remain poorly understood. Here we present high resolution Hi-C data from early Drosophila embryos. We demonstrate that boundaries between topological domains of various sizes map to DNA elements that resemble classical insulator elements: short genomic regions sensitive to DNase digestion that are strongly bound by known insulator proteins and are frequently located between divergent promoters. Further, we show a striking correspondence between these elements and the locations of mapped polytene interband regions. We believe it is likely this relationship between insulators, topological boundaries, and polytene interbands extends across the genome, and we therefore propose a model in which decompaction of boundary-insulator-interband regions drives the organization of interphase chromosomes by creating stable physical separation between adjacent domains.

Drosophila Dosage Compensation Loci Associate with a Boundary-Forming Insulator Complex

Chromatin entry sites (CES) are 100- to 1,500-bp elements that recruit male-specific lethal (MSL) complexes to the X chromosome to upregulate expression of X-linked genes in male flies. CES contain one or more ∼20-bp GA-rich sequences called MSL recognition elements (MREs) that are critical for dosage compensation. Recent studies indicate that CES also correspond to boundaries of X-chromosomal topologically associated domains (TADs). Here, we show that an ∼1,000-kDa complex called the late boundary complex (LBC), which is required for the functioning of the Bithorax complex boundary Fab-7, interacts specifically with a special class of CES that contain multiple MREs. Mutations in the MRE sequences of three of these CES that disrupt function in vivo abrogate interactions with the LBC. Moreover, reducing the levels of two LBC components compromises MSL recruitment. Finally, we show that several of the CES that are physically linked to each other in vivo are LBC interactors.

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.

Drosophila CP190- and dCTCF-mediated enhancer blocking is augmented by SUMOylation

Background

Chromatin insulators shield promoters and chromatin domains from neighboring enhancers or chromatin regions with opposing activities. Insulator-binding proteins and their cofactors mediate the boundary function. In general, covalent modification of proteins by the small ubiquitin-like modifier (SUMO) is an important mechanism to control the interaction of proteins within complexes.

Results

Here we addressed the impact of dSUMO in respect of insulator function, chromatin binding of insulator factors and formation of insulator speckles in Drosophila. SUMOylation augments the enhancer blocking function of four different insulator sequences and increases the genome-wide binding of the insulator cofactor CP190.

Conclusions

These results indicate that enhanced chromatin binding of SUMOylated CP190 causes fusion of insulator speckles, which may allow for more efficient insulation.

Electronic supplementary material

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

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.

X chromosome inactivation: new players in the initiation of gene silencing

X chromosome inactivation (XCI) is a dosage compensation process that was adopted by female mammals to balance gene dosage between XX females and XY males. XCI starts with the upregulation of the non-coding RNA Xist, after which most X-linked genes are silenced and acquire a repressive chromatin state. Even though the chromatin marks of the inactive X have been fairly well described, the mechanisms responsible for the initiation of XCI remain largely unknown. In this review, we discuss recent developments that revealed unexpected factors playing a role in XCI and that might be of crucial importance to understand the mechanisms responsible for the very first steps of this chromosome-wide gene-silencing event.

Establishment of expression-state boundaries by Rif1 and Taz1 in fission yeast

The Shelterin component Rif1 has emerged as a global regulator of the replication-timing program in all eukaryotes examined to date, possibly by modulating the 3D-organization of the genome. In fission yeast a second Shelterin component, Taz1, might share similar functions. Here, we identified unexpected properties for Rif1 and Taz1 by conducting high-throughput genetic screens designed to identify cis- and trans-acting factors capable of creating heterochromatin-euchromatin boundaries in fission yeast. The preponderance of cis-acting elements identified in the screens originated from genomic loci bound by Taz1 and associated with origins of replication whose firing is repressed by Taz1 and Rif1. Boundary formation and gene silencing by these elements required Taz1 and Rif1 and coincided with altered replication timing in the region. Thus, small chromosomal elements sensitive to Taz1 and Rif1 (STAR) could simultaneously regulate gene expression and DNA replication over a large domain, at the edge of which they established a heterochromatin-euchromatin boundary. Taz1, Rif1, and Rif1-associated protein phosphatases Sds21 and Dis2 were each sufficient to establish a boundary when tethered to DNA. Moreover, efficient boundary formation required the amino-terminal domain of the Mcm4 replicative helicase onto which the antagonistic activities of the replication-promoting Dbf4-dependent kinase and Rif1-recruited phosphatases are believed to converge to control replication origin firing. Altogether these observations provide an insight into a coordinated control of DNA replication and organization of the genome into expression domains.

Boundaries of loop domains (insulators): Determinants of chromosome form and function in multicellular eukaryotes

Chromosomes in multicellular animals are subdivided into a series of looped domains. In addition to being the underlying principle for organizing the chromatin fiber, looping is critical for processes ranging from gene regulation to recombination and repair. The subdivision of chromosomes into looped domains depends upon a special class of architectural elements called boundaries or insulators. These elements are distributed throughout the genome and are ubiquitous building blocks of chromosomes. In this review, we focus on features of boundaries that are critical in determining the topology of the looped domains and their genetic properties. We highlight the properties of fly boundaries that are likely to have an important bearing on the organization of looped domains in vertebrates, and discuss the functional consequences of the observed similarities and differences.

Position effects influence HIV latency reversal

The main obstacle to curing HIV is the presence of latent proviruses in the bodies of infected patients. The partial success of reactivation therapies suggests that the genomic context of integrated proviruses can interfere with treatment. Here we developed a method called Barcoded HIV ensembles (B-HIVE) to map the chromosomal locations of thousands of individual proviruses while tracking their transcriptional activities in an infected cell population. B-HIVE revealed that, in Jurkat cells, the expression of HIV is strongest close to endogenous enhancers. The insertion site also affects the response to latency-reversing agents, because we found that phytohemagglutinin and vorinostat reactivated proviruses inserted at distinct genomic locations. From these results, we propose that combinations of drugs targeting all areas of the genome will be most effective. Overall, our data suggest that the insertion context of HIV is a critical determinant of the viral response to reactivation therapies.

Insulated Neighborhoods: Structural and Functional Units of Mammalian Gene Control

Understanding how transcriptional enhancers control over 20,000 protein-coding genes to maintain cell-type-specific gene expression programs in all human cells is a fundamental challenge in regulatory biology. Recent studies suggest that gene regulatory elements and their target genes generally occur within insulated neighborhoods, which are chromosomal loop structures formed by the interaction of two DNA sites bound by the CTCF protein and occupied by the cohesin complex. Here, we review evidence that insulated neighborhoods provide for specific enhancer-gene interactions, are essential for both normal gene activation and repression, form a chromosome scaffold that is largely preserved throughout development, and are perturbed by genetic and epigenetic factors in disease. Insulated neighborhoods are a powerful paradigm for gene control that provides new insights into development and disease.

Transcriptomic data from panarthropods shed new light on the evolution of insulator binding proteins in insects : Insect insulator proteins

Background

Body plan development in multi-cellular organisms is largely determined by homeotic genes. Expression of homeotic genes, in turn, is partially regulated by insulator binding proteins (IBPs). While only a few enhancer blocking IBPs have been identified in vertebrates, the common fruit fly Drosophila melanogaster harbors at least twelve different enhancer blocking IBPs. We screened recently compiled insect transcriptomes from the 1KITE project and genomic and transcriptomic data from public databases, aiming to trace the origin of IBPs in insects and other arthropods.

Results

Our study shows that the last common ancestor of insects (Hexapoda) already possessed a substantial number of IBPs. Specifically, of the known twelve insect IBPs, at least three (i.e., CP190, Su(Hw), and CTCF) already existed prior to the evolution of insects. Furthermore we found GAF orthologs in early branching insect orders, including Zygentoma (silverfish and firebrats) and Diplura (two-pronged bristletails). Mod(mdg4) is most likely a derived feature of Neoptera, while Pita is likely an evolutionary novelty of holometabolous insects. Zw5 appears to be restricted to schizophoran flies, whereas BEAF-32, ZIPIC and the Elba complex, are probably unique to the genus Drosophila. Selection models indicate that insect IBPs evolved under neutral or purifying selection.

Conclusions

Our results suggest that a substantial number of IBPs either pre-date the evolution of insects or evolved early during insect evolution. This suggests an evolutionary history of insulator binding proteins in insects different to that previously thought. Moreover, our study demonstrates the versatility of the 1KITE transcriptomic data for comparative analyses in insects and other arthropods.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3205-1) contains supplementary material, which is available to authorized users.