Drosophila CTCF tandemly aligns with other insulator proteins at the borders of H3K27me3 domains

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

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

Evolution of 3D Chromatin Folding

Studies examining the evolution of genomes have focused mainly on sequence conservation. However, the inner working of a cell implies tightly regulated crosstalk between complex gene networks controlled by small dispersed regulatory elements of physically contacting DNA regions. How these different levels of chromatin organization crosstalk in different species underpins the potential for genome evolutionary plasticity. We review the evolution of chromatin organization across the Animal Tree of Life. We introduce general aspects of the mode and tempo of genome evolution to later explore the multiple layers of genome organization. We argue that both genome and chromosome size modulate patterns of chromatin folding and that chromatin interactions facilitate the formation of lineage-specific chromosomal reorganizations, especially in germ cells. Overall, analyzing the mechanistic forces involved in the maintenance of chromatin structure and function of the germ line is critical for understanding genome evolution, maintenance, and inheritance.

Multiscale 3D genome rewiring during PTF1A-mediated somatic cell reprogramming into neural stem cells

The genome is intricately folded into chromatin compartments, topologically associating domains (TADs) and loops unique to each cell type. How this higher-order genome organization regulates cell fate transition remains elusive. Here we show how a single non-neural progenitor transcription factor, PTF1A, reorchestrates the 3D genome during fibroblast transdifferentiation into neural stem cells (NSCs). Multiomics analyses integrating Hi-C data, PTF1A and CTCF DNA-binding profiles, H3K27ac modification, and gene expression, demonstrate that PTF1A binds to subTAD boundaries subsequently associated with elevated CTCF binding and enhanced boundary insulation, and reorganizes chromatin loops, leading to gene expression changes that drive transdifferentiation into NSCs. Moreover, PTF1A activates enhancers and super-enhancers near low-insulation boundaries and modulates H3K27ac deposition, promoting cell fate transitions. Together, our data implicate an involvement of 3D genome in transcriptional and cell fate alterations, and highlight an essential role for PTF1A in gene expression control and multiscale 3D genome remodeling during cell reprogramming.

Members of an array of zinc-finger proteins specify distinct Hox chromatin boundaries

Partitioning of repressive from actively transcribed chromatin in mammalian cells fosters cell-type-specific gene expression patterns. While this partitioning is reconstructed during differentiation, the chromatin occupancy of the key insulator, CCCTC-binding factor (CTCF), is unchanged at the developmentally important Hox clusters. Thus, dynamic changes in chromatin boundaries must entail other activities. Given its requirement for chromatin loop formation, we examined cohesin-based chromatin occupancy without known insulators, CTCF and Myc-associated zinc-finger protein (MAZ), and identified a family of zinc-finger proteins (ZNFs), some of which exhibit tissue-specific expression. Two such ZNFs foster chromatin boundaries at the Hox clusters that are distinct from each other and from MAZ. PATZ1 was critical to the thoracolumbar boundary in differentiating motor neurons and mouse skeleton, while ZNF263 contributed to cervicothoracic boundaries. We propose that these insulating activities act with cohesin, alone or combinatorially, with or without CTCF, to implement precise positional identity and cell fate during development.

Members of an array of zinc finger proteins specify distinct Hox chromatin boundaries

Partitioning of repressive from actively transcribed chromatin in mammalian cells fosters cell-type specific gene expression patterns. While this partitioning is reconstructed during differentiation, the chromatin occupancy of the key insulator, CTCF, is unchanged at the developmentally important Hox clusters. Thus, dynamic changes in chromatin boundaries must entail other activities. Given its requirement for chromatin loop formation, we examined cohesin-based chromatin occupancy without known insulators, CTCF and MAZ, and identified a family of zinc finger proteins (ZNFs), some of which exhibit tissue-specific expression. Two such ZNFs foster chromatin boundaries at the Hox clusters that are distinct from each other and from MAZ. PATZ1 was critical to the thoracolumbar boundary in differentiating motor neurons and mouse skeleton, while ZNF263 contributed to cervicothoracic boundaries. We propose that these insulating activities act with cohesin, alone or combinatorially, with or without CTCF, to implement precise positional identity and cell fate during development.

Chromatin insulator mechanisms ensure accurate gene expression by controlling overall 3D genome organization

Chromatin insulators are DNA-protein complexes that promote specificity of enhancer-promoter interactions and maintain distinct transcriptional states through control of 3D genome organization. In this review, we highlight recent work visualizing how mammalian CCCTC-binding factor acts as a boundary to dynamic DNA loop extrusion mediated by cohesin. We also discuss new studies in both mammals and Drosophila that elucidate biological redundancy of chromatin insulator function and interplay with transcription with respect to topologically associating domain formation. Finally, we present novel concepts in spatiotemporal regulation of chromatin insulator function during differentiation and development and possible consequences of disrupted insulator activity on cellular proliferation.

Metabolic regulation of CTCF expression and chromatin association dictates starvation response in mice and flies

Coordinated temporal control of gene expression is essential for physiological homeostasis, especially during metabolic transitions. However, the interplay between chromatin architectural proteins and metabolism in regulating transcription is less understood. Here, we demonstrate a conserved bidirectional interplay between CTCF (CCCTC-binding factor) expression/function and metabolic inputs during feed-fast cycles. Our results indicate that its loci-specific functional diversity is associated with physiological plasticity in mouse hepatocytes. CTCF differential expression and long non-coding RNA-Jpx mediated changes in chromatin occupancy, unraveled its paradoxical yet tuneable functions, which are governed by metabolic inputs. We illustrate the key role of CTCF in controlling temporal cascade of transcriptional response, with effects on hepatic mitochondrial energetics and lipidome. Underscoring the evolutionary conservation of CTCF-dependent metabolic homeostasis, CTCF knockdown in flies abrogated starvation resistance. In summary, we demonstrate the interplay between CTCF and metabolic inputs that highlights the coupled plasticity of physiological responses and chromatin function.

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

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

Topological screen identifies hundreds of Cp190- and CTCF-dependent Drosophila chromatin insulator 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.

Low-affinity CTCF binding drives transcriptional regulation whereas high-affinity binding encompasses architectural functions

CTCF is a DNA-binding protein which plays critical roles in chromatin structure organization and transcriptional regulation; however, little is known about the functional determinants of different CTCF-binding sites (CBS). Using a conditional mouse model, we have identified one set of CBSs that are lost upon CTCF depletion (lost CBSs) and another set that persists (retained CBSs). Retained CBSs are more similar to the consensus CTCF-binding sequence and usually span tandem CTCF peaks. Lost CBSs are enriched at enhancers and promoters and associate with active chromatin marks and higher transcriptional activity. In contrast, retained CBSs are enriched at TAD and loop boundaries. Integration of ChIP-seq and RNA-seq data has revealed that retained CBSs are located at the boundaries between distinct chromatin states, acting as chromatin barriers. Our results provide evidence that transient, lost CBSs are involved in transcriptional regulation, whereas retained CBSs are critical for establishing higher-order chromatin architecture.

Nickel-induced alterations to chromatin structure and function

Nickel (Ni), a heavy metal is prevalent in the atmosphere due to both natural and anthropogenic activities. Ni is a carcinogen implicated in the development of lung and nasal cancers in humans. Furthermore, Ni exposure is associated with a number of chronic lung diseases in humans including asthma, chronic bronchitis, emphysema, pulmonary fibrosis, pulmonary edema and chronic obstructive pulmonary disease (COPD). While Ni compounds are weak mutagens, a number of studies have demonstrated the potential of Ni to alter the epigenome, suggesting epigenomic dysregulation as an important underlying cause for its pathogenicity. In the eukaryotic nucleus, the DNA is organized in a three-dimensional (3D) space through assembly of higher order chromatin structures. Such an organization is critically important for transcription and other biological activities. Accumulating evidence suggests that by negatively affecting various cellular regulatory processes, Ni could potentially affect chromatin organization. In this review, we discuss the role of Ni in altering the chromatin architecture, which potentially plays a major role in Ni pathogenicity.

Chromatin accessibility-based characterisation of brain gene regulatory networks in three distinct honey bee polyphenisms

The honey bee genome has the capacity to produce three phenotypically distinct organisms (two diploid female castes: queen and worker, and a haploid male drone). Previous studies have implicated metabolic flux acting via epigenetic regulation in directing nutrition-driven phenotypic plasticity in the honey bee. However, the cis-acting DNA regulatory elements that establish tissue and polyphenism -specific epigenomes and gene expression programmes, remain unclear. Using a high resolution multiomic approach including assay for transposase-accessible chromatin by sequencing (ATAC-seq), RNA-seq and ChIP-seq, we produce the first genome-wide maps of the regulatory landscape across all three adult honey bee phenotypes identifying > 5000 regulatory regions in queen, 7500 in worker and 6500 in drone, with the vast majority of these sites located within intronic regions. These regions are defined by positive enrichment of H3K27ac and depletion of H3K4me3 and show a positive correlation with gene expression. Using ATAC-seq footprinting we determine queen, worker and drone -specific transcription factor occupancy and uncover novel phenotype-specific regulatory networks identifying two key nuclear receptors that have previously been implicated in caste-determination and adult behavioural maturation in honey bees; ecdysone receptor and ultraspiracle. Collectively, this study provides novel insights into key gene regulatory networks that are associated with these distinct polyphenisms in the honey bee.

Lint‐O cooperates with L(3)mbt in target gene suppression to maintain homeostasis in fly ovary and brain

Loss‐of‐function mutations in Drosophila lethal(3)malignant brain tumor [l(3)mbt] cause ectopic expression of germline genes and brain tumors. Loss of L(3)mbt function in ovarian somatic cells (OSCs) aberrantly activates germ‐specific piRNA amplification and leads to infertility. However, the underlying mechanism remains unclear. Here, ChIP‐seq for L(3)mbt in cultured OSCs and RNA‐seq before and after L(3)mbt depletion shows that L(3)mbt genomic binding is not necessarily linked to gene regulation and that L(3)mbt controls piRNA pathway genes in multiple ways. Lack of known L(3)mbt co‐repressors, such as Lint‐1, has little effect on the levels of piRNA amplifiers. Identification of L(3)mbt interactors in OSCs and subsequent analysis reveals CG2662 as a novel co‐regulator of L(3)mbt, termed “L(3)mbt interactor in OSCs” (Lint‐O). Most of the L(3)mbt‐bound piRNA amplifier genes are also bound by Lint‐O in a similar fashion. Loss of Lint‐O impacts the levels of piRNA amplifiers, similar to the lack of L(3)mbt. The lint‐O‐deficient flies exhibit female sterility and tumorous brains. Thus, L(3)mbt and its novel co‐suppressor Lint‐O cooperate in suppressing target genes to maintain homeostasis in the ovary and brain.

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.

Implications of Dosage Deficiencies in CTCF and Cohesin on Genome Organization, Gene Expression, and Human Neurodevelopment

Properly organizing DNA within the nucleus is critical to ensure normal downstream nuclear functions. CTCF and cohesin act as major architectural proteins, working in concert to generate thousands of high-intensity chromatin loops. Due to their central role in loop formation, a massive research effort has been dedicated to investigating the mechanism by which CTCF and cohesin create these loops. Recent results lead to questioning the direct impact of CTCF loops on gene expression. Additionally, results of controlled depletion experiments in cell lines has indicated that genome architecture may be somewhat resistant to incomplete deficiencies in CTCF or cohesin. However, heterozygous human genetic deficiencies in CTCF and cohesin have illustrated the importance of their dosage in genome architecture, cellular processes, animal behavior, and disease phenotypes. Thus, the importance of considering CTCF or cohesin levels is especially made clear by these heterozygous germline variants that characterize genetic syndromes, which are increasingly recognized in clinical practice. Defined primarily by developmental delay and intellectual disability, the phenotypes of CTCF and cohesin deficiency illustrate the importance of architectural proteins particularly in neurodevelopment. We discuss the distinct roles of CTCF and cohesin in forming chromatin loops, highlight the major role that dosage of each protein plays in the amplitude of observed effects on gene expression, and contrast these results to heterozygous mutation phenotypes in murine models and clinical patients. Insights highlighted by this comparison have implications for future research into these newly emerging genetic syndromes.

CRISPR and biochemical screens identify MAZ as a cofactor in CTCF-mediated insulation at Hox clusters

CCCTC-binding factor (CTCF) is critical to three-dimensional genome organization. Upon differentiation, CTCF insulates active and repressed genes within Hox gene clusters. We conducted a genome-wide CRISPR knockout (KO) screen to identify genes required for CTCF-boundary activity at the HoxA cluster, complemented by biochemical approaches. Among the candidates, we identified Myc-associated zinc-finger protein (MAZ) as a cofactor in CTCF insulation. MAZ colocalizes with CTCF at chromatin borders and, similar to CTCF, interacts with the cohesin subunit RAD21. MAZ KO disrupts gene expression and local contacts within topologically associating domains. Similar to CTCF motif deletions, MAZ motif deletions lead to derepression of posterior Hox genes immediately after CTCF boundaries upon differentiation, giving rise to homeotic transformations in mouse. Thus, MAZ is a factor contributing to appropriate insulation, gene expression and genomic architecture during development.

Three classes of epigenomic regulators converge to hyperactivate the essential maternal gene deadhead within a heterochromatin mini-domain

The formation of a diploid zygote is a highly complex cellular process that is entirely controlled by maternal gene products stored in the egg cytoplasm. This highly specialized transcriptional program is tightly controlled at the chromatin level in the female germline. As an extreme case in point, the massive and specific ovarian expression of the essential thioredoxin Deadhead (DHD) is critically regulated in Drosophila by the histone demethylase Lid and its partner, the histone deacetylase complex Sin3A/Rpd3, via yet unknown mechanisms. Here, we identified Snr1 and Mod(mdg4) as essential for dhd expression and investigated how these epigenomic effectors act with Lid and Sin3A to hyperactivate dhd. Using Cut&Run chromatin profiling with a dedicated data analysis procedure, we found that dhd is intriguingly embedded in an H3K27me3/H3K9me3-enriched mini-domain flanked by DNA regulatory elements, including a dhd promoter-proximal element essential for its expression. Surprisingly, Lid, Sin3a, Snr1 and Mod(mdg4) impact H3K27me3 and this regulatory element in distinct manners. However, we show that these effectors activate dhd independently of H3K27me3/H3K9me3, and that dhd remains silent in the absence of these marks. Together, our study demonstrates an atypical and critical role for chromatin regulators Lid, Sin3A, Snr1 and Mod(mdg4) to trigger tissue-specific hyperactivation within a unique heterochromatin mini-domain.

Multicellular development depends on a tight control of gene expression in each cell type. This relies on the coordinated activities of nuclear proteins that interact with DNA or its histone scaffold to promote or restrict gene transcription. For example, we previously showed that the histone modifying enzymes Lid and Sin3A/Rpd3 are required in Drosophila ovaries for the massive expression of deadhead (dhd), a gene encoding for a thioredoxin that is essential for fertility. In this paper, we have further identified two additional dhd regulators, Snr1 and Mod(mdg4) and dissected the mechanism behind hyperactivation of this gene. Using the epigenomic profiling method Cut&Run with a dedicated data analysis approach, we unexpectedly found that dhd is embedded in an unusual chromatin mini-domain featuring repressive histone modifications H3K27me3 and H3K9me3 and flanked by two regulatory elements. However, we further showed that Lid, Sin3A, Snr1 and Mod(mdg4) behave like obligatory activators of dhd independently of this mini-domain. Our study unveils how multiple broad-acting epigenomic effectors operate in non-canonical manners to ensure a critical and specialized gene activation event. These findings challenge our knowledge on these regulatory mechanisms and their roles in development and pathology.

An explainable artificial intelligence approach for decoding the enhancer histone modifications code and identification of novel enhancers in Drosophila

BACKGROUND

Enhancers are non-coding regions of the genome that control the activity of target genes. Recent efforts to identify active enhancers experimentally and in silico have proven effective. While these tools can predict the locations of enhancers with a high degree of accuracy, the mechanisms underpinning the activity of enhancers are often unclear.

RESULTS

Using machine learning (ML) and a rule-based explainable artificial intelligence (XAI) model, we demonstrate that we can predict the location of known enhancers in Drosophila with a high degree of accuracy. Most importantly, we use the rules of the XAI model to provide insight into the underlying combinatorial histone modifications code of enhancers. In addition, we identified a large set of putative enhancers that display the same epigenetic signature as enhancers identified experimentally. These putative enhancers are enriched in nascent transcription, divergent transcription and have 3D contacts with promoters of transcribed genes. However, they display only intermediary enrichment of mediator and cohesin complexes compared to previously characterised active enhancers. We also found that 10-15% of the predicted enhancers display similar characteristics to super enhancers observed in other species.

CONCLUSIONS

Here, we applied an explainable AI model to predict enhancers with high accuracy. Most importantly, we identified that different combinations of epigenetic marks characterise different groups of enhancers. Finally, we discovered a large set of putative enhancers which display similar characteristics with previously characterised active enhancers.

CGGBP1-dependent CTCF-binding sites restrict ectopic transcription

Binding sites of the chromatin regulator protein CTCF function as important landmarks in the human genome. The recently characterized CTCF-binding sites at LINE-1 repeats depend on another repeat-regulatory protein CGGBP1. These CGGBP1-dependent CTCF-binding sites serve as potential barrier elements for epigenetic marks such as H3K9me3. Such CTCF-binding sites are associated with asymmetric H3K9me3 levels as well as RNA levels in their flanks. The functions of these CGGBP1-dependent CTCF-binding sites remain unknown. By performing targeted studies on candidate CGGBP1-dependent CTCF-binding sites cloned in an SV40 promoter-enhancer episomal system we show that these regions act as inhibitors of ectopic transcription from the SV40 promoter. CGGBP1-dependent CTCF-binding sites that recapitulate their genomic function of loss of CTCF binding upon CGGBP1 depletion and H3K9me3 asymmetry in immediate flanks are also the ones that show the strongest inhibition of ectopic transcription. By performing a series of strand-specific reverse transcription PCRs we demonstrate that this ectopic transcription results in the synthesis of RNA from the SV40 promoter in a direction opposite to the downstream reporter gene in a strand-specific manner. The unleashing of the bidirectionality of the SV40 promoter activity and a breach of the transcription barrier seems to depend on depletion of CGGBP1 and loss of CTCF binding proximal to the SV40 promoter. RNA-sequencing reveals that CGGBP1-regulated CTCF-binding sites act as barriers to transcription at multiple locations genome-wide. These findings suggest a role of CGGBP1-dependent binding sites in restricting ectopic transcription.

Evolution and genomic organization of the insect sHSP gene cluster and coordinate regulation in phenotypic plasticity

Background

Conserved syntenic gene complexes are rare in Arthropods and likely only retained due to functional constraint. Numerous sHSPs have been identified in the genomes of insects, some of which are located clustered in close proximity. Previous phylogenetic analyses of these clustered sHSP have been limited to a small number of holometabolous insect species and have not determined the pattern of evolution of the clustered sHSP genes (sHSP-C) in insect or Arthropod lineages.

Results

Using eight genomes from representative insect orders and three non-insect arthropod genomes we have identified that a syntenic cluster of sHSPs (sHSP-C) is a hallmark of most Arthropod genomes. Using 11 genomes from Hymenopteran species our phylogenetic analyses have refined the evolution of the sHSP-C in Hymenoptera and found that the sHSP-C is order-specific with evidence of birth-and-death evolution in the hymenopteran lineage. Finally we have shown that the honeybee sHSP-C is co-ordinately expressed and is marked by genomic features, including H3K27me3 histone marks consistent with coordinate regulation, during honeybee ovary activation.

Conclusions

The syntenic sHSP-C is present in most insect genomes, and its conserved coordinate expression and regulation implies that it is an integral genomic component of environmental response in arthropods.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01885-8.

Multi-Scale Organization of the Drosophila melanogaster Genome

Interphase chromatin, despite its appearance, is a highly organized framework of loops and bends. Chromosomes are folded into topologically associating domains, or TADs, and each chromosome and its homolog occupy a distinct territory within the nucleus. In Drosophila, genome organization is exceptional because homologous chromosome pairing is in both germline and somatic tissues, which promote interhomolog interactions such as transvection that can affect gene expression in trans. In this review, we focus on what is known about genome organization in Drosophila and discuss it from TADs to territory. We start by examining intrachromosomal organization at the sub-chromosome level into TADs, followed by a comprehensive analysis of the known proteins that play a key role in TAD formation and boundary establishment. We then zoom out to examine interhomolog interactions such as pairing and transvection that are abundant in Drosophila but rare in other model systems. Finally, we discuss chromosome territories that form within the nucleus, resulting in a complete picture of the multi-scale organization of the Drosophila genome.

The Dm-Myb Oncoprotein Contributes to Insulator Function and Stabilizes Repressive H3K27me3 PcG Domains

Drosophila Myb (Dm-Myb) encodes a protein that plays a key role in regulation of mitotic phase genes. Here, we further refine its role in the context of a developing tissue as a potentiator of gene expression required for proper RNA polymerase II (RNA Pol II) function and efficient H3K4 methylation at promoters. In contrast to its role in gene activation, Myb is also required for repression of many genes, although no specific mechanism for this role has been proposed. We now reveal a critical role for Myb in contributing to insulator function, in part by promoting binding of insulator proteins BEAF-32 and CP190 and stabilizing H3K27me3 Polycomb-group (PcG) domains. In the absence of Myb, H3K27me3 is markedly reduced throughout the genome, leading to H3K4me3 spreading and gene derepression. Finally, Myb is enriched at boundaries that demarcate chromatin environments, including chromatin loop anchors. These results reveal functions of Myb that extend beyond transcriptional regulation.

Myb has been considered a transcriptional activator of primarily M phase genes. Here, Santana et al. show that Myb also contributes to insulator function, in part by promoting binding of insulator factors, and is required to stabilize repressive domains in the genome.

CTCF loss has limited effects on global genome architecture in Drosophila despite critical regulatory functions

Vertebrate genomes are partitioned into contact domains defined by enhanced internal contact frequency and formed by two principal mechanisms: compartmentalization of transcriptionally active and inactive domains, and stalling of chromosomal loop-extruding cohesin by CTCF bound at domain boundaries. While Drosophila has widespread contact domains and CTCF, it is currently unclear whether CTCF-dependent domains exist in flies. We genetically ablate CTCF in Drosophila and examine impacts on genome folding and transcriptional regulation in the central nervous system. We find that CTCF is required to form a small fraction of all domain boundaries, while critically controlling expression patterns of certain genes and supporting nervous system function. We also find that CTCF recruits the pervasive boundary-associated factor Cp190 to CTCF-occupied boundaries and co-regulates a subset of genes near boundaries together with Cp190. These results highlight a profound difference in CTCF-requirement for genome folding in flies and vertebrates, in which a large fraction of boundaries are CTCF-dependent and suggest that CTCF has played mutable roles in genome architecture and direct gene expression control during metazoan evolution.

SetDB1 and Su(var)3-9 play non-overlapping roles in somatic cell chromosomes of Drosophila melanogaster

We explored functional roles of two H3K9-specific histone methyltransferases of Drosophila melanogaster, SetDB1 (also known as Eggless) and Su(var)3-9. Using the DamID approach, we generated the binding profile for SetDB1 in Drosophila salivary gland chromosomes, and matched it to the profile of Su(var)3-9. Unlike Su(var)3-9, SetDB1 turned out to be an euchromatic protein that is absent from repeated DNA compartments, and is largely restricted to transcription start sites (TSSs) and 5' untranslated regions (5'UTRs) of ubiquitously expressed genes. Significant SetDB1 association is also observed at binding sites for the insulator protein CP190. SetDB1 and H3K9 di- and tri-methylated (me2 and me3)-enriched sites tend to display poor overlap. At the same time, SetDB1 has a clear connection with the distribution of H3K27me3 mark. SetDB1 binds outside the domains possessing this modification, and about half of the borders of H3K27me3 domains are decorated by SetDB1 together with actively transcribed genes. On the basis of poor correlation between the distribution of SetDB1 and H3K9 methylation marks, we speculate that, in somatic cells, SetDB1 may contribute to the methylation of a broader set of chromosomal proteins than just H3K9. In addition, SetDB1 can be expected to play a role in the establishment of chromatin functional domains.

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.

The genetic factors of bilaterian evolution

The Cambrian explosion was a unique animal radiation ~540 million years ago that produced the full range of body plans across bilaterians. The genetic mechanisms underlying these events are unknown, leaving a fundamental question in evolutionary biology unanswered. Using large-scale comparative genomics and advanced orthology evaluation techniques, we identified 157 bilaterian-specific genes. They include the entire Nodal pathway, a key regulator of mesoderm development and left-right axis specification; components for nervous system development, including a suite of G-protein-coupled receptors that control physiology and behaviour, the Robo-Slit midline repulsion system, and the neurotrophin signalling system; a high number of zinc finger transcription factors; and novel factors that previously escaped attention. Contradicting the current view, our study reveals that genes with bilaterian origin are robustly associated with key features in extant bilaterians, suggesting a causal relationship.

A multi-layered structure of the interphase chromocenter revealed by proximity-based biotinylation

During interphase centromeres often coalesce into a small number of chromocenters, which can be visualized as distinct, DAPI dense nuclear domains. Intact chromocenters play a major role in maintaining genome stability as they stabilize the transcriptionally silent state of repetitive DNA while ensuring centromere function. Despite its biological importance, relatively little is known about the molecular composition of the chromocenter or the processes that mediate chromocenter formation and maintenance. To provide a deeper molecular insight into the composition of the chromocenter and to demonstrate the usefulness of proximity-based biotinylation as a tool to investigate those questions, we performed super resolution microscopy and proximity-based biotinylation experiments of three distinct proteins associated with the chromocenter in Drosophila. Our work revealed an intricate internal architecture of the chromocenter suggesting a complex multilayered structure of this intranuclear domain.

Exploring the changing landscape of cell-to-cell variation after CTCF knockdown via single cell RNA-seq

BACKGROUND

CCCTC-Binding Factor (CTCF), also known as 11-zinc finger protein, participates in many cellular processes, including insulator activity, transcriptional regulation and organization of chromatin architecture. Based on single cell flow cytometry and single cell RNA-FISH analyses, our previous study showed that deletion of CTCF binding site led to a significantly increase of cellular variation of its target gene. However, the effect of CTCF on genome-wide landscape of cell-to-cell variation remains unclear.

RESULTS

We knocked down CTCF in EL4 cells using shRNA, and conducted single cell RNA-seq on both wild type (WT) cells and CTCF-Knockdown (CTCF-KD) cells using Fluidigm C1 system. Principal component analysis of single cell RNA-seq data showed that WT and CTCF-KD cells concentrated in two different clusters on PC1, indicating that gene expression profiles of WT and CTCF-KD cells were systematically different. Interestingly, GO terms including regulation of transcription, DNA binding, zinc finger and transcription factor binding were significantly enriched in CTCF-KD-specific highly variable genes, implying tissue-specific genes such as transcription factors were highly sensitive to CTCF level. The dysregulation of transcription factors potentially explains why knockdown of CTCF leads to systematic change of gene expression. In contrast, housekeeping genes such as rRNA processing, DNA repair and tRNA processing were significantly enriched in WT-specific highly variable genes, potentially due to a higher cellular variation of cell activity in WT cells compared to CTCF-KD cells. We further found that cellular variation-increased genes were significantly enriched in down-regulated genes, indicating CTCF knockdown simultaneously reduced the expression levels and increased the expression noise of its regulated genes.

CONCLUSIONS

To our knowledge, this is the first attempt to explore genome-wide landscape of cellular variation after CTCF knockdown. Our study not only advances our understanding of CTCF function in maintaining gene expression and reducing expression noise, but also provides a framework for examining gene function.

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

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

The Role of Insulation in Patterning Gene Expression

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

ATP, Mg2+, Nuclear Phase Separation, and Genome Accessibility

Misregulation of the processes controlling eukaryotic gene expression can result in disease. Gene expression is influenced by the surrounding chromatin; hence the nuclear environment is also of vital importance. Recently, understanding of chromatin hierarchical folding has increased together with the discovery of membrane-less organelles which are distinct, dynamic liquid droplets that merge and expand within the nucleus. These 'sieve'-like regions may compartmentalize and separate functionally distinct regions of chromatin. This article aims to discuss recent studies on nuclear phase within the context of poly(ADP-ribose), ATP, and Mg²⁺ levels, and we propose a combinatorial complex role for these molecules in phase separation and genome regulation. We also discuss the implications of this process for gene regulation and discuss possible strategies to test this.

The same domain of Su(Hw) is required for enhancer blocking and direct promoter repression

Suppressor of Hairy-wing [Su(Hw)] is a DNA-binding architectural protein that participates in the organization of insulators and repression of promoters in Drosophila. This protein contains acidic regions at both ends and a central cluster of 12 zinc finger domains, some of which are involved in the specific recognition of the binding site. One of the well-described in vivo function of Su(Hw) is the repression of transcription of neuronal genes in oocytes. Here, we have found that the same Su(Hw) C-terminal region (aa 720–892) is required for insulation as well as for promoter repression. The best characterized partners of Su(Hw), CP190 and Mod(mdg4)-67.2, are not involved in the repression of neuronal genes. Taken together, these results suggest that an unknown protein or protein complex binds to the C-terminal region of Su(Hw) and is responsible for the direct repression activity of Su(Hw).

Spatial Genome Organization: From Development to Disease

Every living organism, from bacteria to humans, contains DNA encoding anything from a few hundred genes in intracellular parasites such as Mycoplasma, up to several tens of thousands in many higher organisms. The first observations indicating that the nucleus had some kind of organization were made over a hundred years ago. Understanding of its significance is both limited and aided by the development of techniques, in particular electron microscopy, fluorescence in situ hybridization, DamID and most recently HiC. As our knowledge about genome organization grows, it becomes apparent that the mechanisms are conserved in evolution, even if the individual players may vary. These mechanisms involve DNA binding proteins such as histones, and a number of architectural proteins, some of which are very much conserved, with some others having diversified and multiplied, acquiring specific regulatory functions. In this review we will look at the principles of genome organization in a hierarchical manner, from DNA packaging to higher order genome associations such as TADs, and the significance of radial positioning of genomic loci. We will then elaborate on the dynamics of genome organization during development, and how genome architecture plays an important role in cell fate determination. Finally, we will discuss how misregulation can be a factor in human disease.

IDH1: Linking Metabolism and Epigenetics

Mutations in genes encoding enzymes of the tricarboxylic acid cycle often contribute to cancer development and progression by disrupting cell metabolism and altering the epigenetic landscape. This is exemplified by the isoforms of isocitrate dehydrogenase (IDH1/2), which metabolize isocitrate to α-Ketoglutarate (α-KG). Gain of function mutations in IDH1 or IDH2 result in reduced levels of α-KG as a result of increased formation of D-2-Hydroxyglutarate (2-HG). α-KG is an essential co-factor for certain histone and DNA demethylases, while 2-HG is a competitive inhibitor. These IDH1/2 mutations are thought to result in hypermethylated histones and DNA which in turn alters gene expression and drives cancer progression. While this model seems to be generally accepted in the field, the exact molecular mechanisms still remain elusive. How much of this model has been rigorously demonstrated and what is just being assumed? Are the effects genome-wide or focused on specific loci? This Perspective aims at elucidating the key questions that remain to be addressed, the experimental techniques that could be used to gain further insight into the molecular mechanisms involved and the additional consequences of these mutations beyond DNA and protein methylation.

Selective interactions between diverse STEs organize the ANT-C Hox cluster

The three-dimensional organization of the eukaryotic genome is important for its structure and function. Recent studies indicate that hierarchies of chromatin loops underlie important aspects of both genomic organization and gene regulation. Looping between insulator or boundary elements interferes with enhancer-promoter communications and limits the spread active or repressive organized chromatin. We have used the SF1 insulator in the Drosophila Antennapedia homeotic gene complex (ANT-C) as a model to study the mechanism and regulation of chromatin looping events. We reported previously that SF1 tethers a transient chromatin loop in the early embryo that insulates the Hox gene Sex comb reduce from the neighbor non-Hox gene fushi tarazu for their independent regulation. To further probe the functional range and connectivity of SF1, we used high-resolution chromosomal conformation capture (3C) to search for SF1 looping partners across ANT-C. We report here the identification of three distal SF1 Tether Elements (STEs) located in the labial, Deformed and Antennapedia Hox gene regions, extending the range of SF1 looping network to the entire complex. These novel STEs are bound by four different combinations of insulator proteins and exhibit distinct behaviors in enhancer block, enhancer-bypass and boundary functions. Significantly, the six STEs we identified so far map to all but one of the major boundaries between repressive and active histone domains, underlining the functional relevance of these long-range chromatin loops in organizing the Hox complex. Importantly, SF1 selectively captured with only 5 STEs out of ~20 sites that display similar insulator binding profiles, indicating that presence of insulator proteins alone is not sufficient to determine looping events. These findings suggest that selective interaction among diverse STE insulators organize the Drosophila Hox genes in the 3D nuclear space.

The Ubx Polycomb response element bypasses an unpaired Fab-8 insulator via cis transvection in Drosophila

Chromatin insulators or boundary elements protect genes from regulatory activities from neighboring genes or chromatin domains. In the Drosophila Abdominal-B (Abd-B) locus, the deletion of such elements, such as Frontabdominal-7 (Fab-7) or Fab-8 led to dominant gain of function phenotypes, presumably due to the loss of chromatin barriers. Homologous chromosomes are paired in Drosophila, creating a number of pairing dependent phenomena including transvection, and whether transvection may affect the function of Polycomb response elements (PREs) and thus contribute to the phenotypes are not known. Here, we studied the chromatin barrier activity of Fab-8 and how it is affected by the zygosity of the transgene, and found that Fab-8 is able to block the silencing effect of the Ubx PRE on the DsRed reporter gene in a CTCF binding sites dependent manner. However, the blocking also depends on the zygosity of the transgene in that the barrier activity is present when the transgene is homozygous, but absent when the transgene is heterozygous. To analyze this effect, we performed chromatin immunoprecipitation and quantitative PCR (ChIP-qPCR) experiments on homozygous transgenic embryos, and found that H3K27me3 and H3K9me3 marks are restricted by Fab-8, but they spread beyond Fab-8 into the DsRed gene when the two CTCF binding sites within Fab-8 were mutated. Consistent with this, the mutation reduced H3K4me3 and RNA Pol II binding to the DsRed gene, and consequently, DsRed expression. Importantly, in heterozygous embryos, Fab-8 is unable to prevent the spread of H3K27me3 and H3K9me3 marks from crossing Fab-8 into DsRed, suggesting an insulator bypass. These results suggest that in the Abd-B locus, deletion of the insulator in one copy of the chromosome could lead to the loss of insulator activity on the homologous chromosome, and in other loci where chromosomal deletion created hemizygous regions of the genome, the chromatin barrier could be compromised. This study highlights a role of homologous chromosome pairing in the regulation of gene expression in the Drosophila genome.

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

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

Sub-kb Hi-C in D. melanogaster reveals conserved characteristics of TADs between insect and mammalian cells

Topologically associating domains (TADs) are fundamental elements of the eukaryotic genomic structure. However, recent studies suggest that the insulating complexes, CTCF/cohesin, present at TAD borders in mammals are absent from those in Drosophila melanogaster, raising the possibility that border elements are not conserved among metazoans. Using in situ Hi-C with sub-kb resolution, here we show that the D. melanogaster genome is almost completely partitioned into >4000 TADs, nearly sevenfold more than previously identified. The overwhelming majority of these TADs are demarcated by the insulator complexes, BEAF-32/CP190, or BEAF-32/Chromator, indicating that these proteins may play an analogous role in flies as that of CTCF/cohesin in mammals. Moreover, extended regions previously thought to be unstructured are shown to consist of small contiguous TADs, a property also observed in mammals upon re-examination. Altogether, our work demonstrates that fundamental features associated with the higher-order folding of the genome are conserved from insects to mammals.

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

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

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

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

Evolutionarily Conserved Principles Predict 3D Chromatin Organization

Topologically associating domains (TADs), CTCF loop domains, and A/B compartments have been identified as important structural and functional components of 3D chromatin organization, yet the relationship between these features is not well understood. Using high-resolution Hi-C and HiChIP, we show that Drosophila chromatin is organized into domains we term compartmental domains that correspond precisely with A/B compartments at high resolution. We find that transcriptional state is a major predictor of Hi-C contact maps in several eukaryotes tested, including C. elegans and A. thaliana. Architectural proteins insulate compartmental domains by reducing interaction frequencies between neighboring regions in Drosophila, but CTCF loops do not play a distinct role in this organism. In mammals, compartmental domains exist alongside CTCF loop domains to form topological domains. The results suggest that compartmental domains are responsible for domain structure in all eukaryotes, with CTCF playing an important role in domain formation in mammals.

Chromatin Domains: The Unit of Chromosome Organization

How eukaryotic chromosomes fold inside the nucleus is an age-old question that remains unanswered today. Early biochemical and microscopic studies revealed the existence of chromatin domains and loops as a pervasive feature of interphase chromosomes, but the biological implications of such organizational features were obscure. Genome-wide analysis of pair-wise chromatin interactions using chromatin conformation capture (3C)-based techniques has shed new light on the organization of chromosomes in interphase nuclei. Particularly, the finding of cell-type invariant, evolutionarily conserved topologically associating domains (TADs) in a broad spectrum of cell types has provided a new molecular framework for the study of animal development and human diseases. Here, we review recent progress in characterization of such chromatin domains and delineation of mechanisms of their formation in animal cells.

Chromatin binding of Gcn5 in Drosophila is largely mediated by CP190

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

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

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

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.

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

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

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

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

Chromatin boundary elements organize genomic architecture and developmental gene regulation in Drosophila Hox clusters

The three-dimensional (3D) organization of the eukaryotic genome is critical for its proper function. Evidence suggests that extensive chromatin loops form the building blocks of the genomic architecture, separating genes and gene clusters into distinct functional domains. These loops are anchored in part by a special type of DNA elements called chromatin boundary elements (CBEs). CBEs were originally found to insulate neighboring genes by blocking influences of transcriptional enhancers or the spread of silent chromatin. However, recent results show that chromatin loops can also play a positive role in gene regulation by looping out intervening DNA and “delivering” remote enhancers to gene promoters. In addition, studies from human and model organisms indicate that the configuration of chromatin loops, many of which are tethered by CBEs, is dynamically regulated during cell differentiation. In particular, a recent work by Li et al has shown that the SF1 boundary, located in the Drosophila Hox cluster, regulates local genes by tethering different subsets of chromatin loops: One subset enclose a neighboring gene ftz, limiting its access by the surrounding Scr enhancers and restrict the spread of repressive histones during early embryogenesis; and the other loops subdivide the Scr regulatory region into independent domains of enhancer accessibility. The enhancer-blocking activity of these CBE elements varies greatly in strength and tissue distribution. Further, tandem pairing of SF1 and SF2 facilitate the bypass of distal enhancers in transgenic flies, providing a mechanism for endogenous enhancers to circumvent genomic interruptions resulting from chromosomal rearrangement. This study demonstrates how a network of chromatin boundaries, centrally organized by SF1, can remodel the 3D genome to facilitate gene regulation during development.