Dual functionality of cis-regulatory elements as developmental enhancers and Polycomb response elements

Abundant repressor binding sites in human enhancers are associated with the fine-tuning of gene regulation

The regulation of gene expression relies on the coordinated action of transcription factors (TFs) at enhancers, including both activator and repressor TFs. We employed deep learning (DL) to dissect HepG2 enhancers into positive (PAR), negative (NAR), and neutral activity regions. Sharpr-MPRA and STARR-seq highlight the dichotomy impact of NARs and PARs on modulating and catalyzing the activity of enhancers, respectively. Approximately 22% of HepG2 enhancers, termed "repressive impact enhancers" (RIEs), are predominantly populated by NARs and transcriptional repression motifs. Genes flanking RIEs exhibit a stage-specific decline in expression during late development, suggesting RIEs' role in trimming enhancer activities. About 16.7% of human NARs emerge from neutral rhesus macaque DNA. This gain of repressor binding sites in RIEs is associated with a 30% decrease in the average expression of flanking genes in humans compared to rhesus macaque. Our work reveals modulated enhancer activity and adaptable gene regulation through the evolutionary dynamics of TF binding sites.

Comprehensive Predictions of Mef2-Mediated Chromatin Loops, Which May Inhibit Ubx Binding by Blocking Low-Affinity Binding Sites

Gene regulation depends on the interaction between chromatin-associated factors, such as transcription factors (TFs), which promote chromatin loops to ensure tight contact between enhancer and promoter regions. So far, positive interactions that lead to gene activation have been the main focus of research, but regulations related to blocking or inhibiting factor binding are also essential to maintaining a defined cellular status. To understand these interactions in greater detail, I investigated the possibility of the muscle differentiation factor Mef2 to prevent early Hox factor binding, leading to the proper timing of regulatory processes and the activation of differentiation events. My investigations relied on a collection of publicly available genome-wide binding data sets of Mef2 and Ubx (as the Hox factor), Capture-C interactions, and ATAC-seq analysis in Mef2 mutant cells. The analysis indicated that Mef2 can form possible chromatin loops to Ubx-bound regions. These regions contain low-affinity Ubx binding sites, and the chromatin architecture is independent of Mef2's function. High levels of Ubx may disrupt the loops and allow specific Ubx bindings to regulate defined targets. In summary, my investigations highlight that the use of many publicly available data sets enables computational approaches to make robust predictions and, for the first time, suggest a molecular function of Mef2 as a preventer of Hox binding, indicating that it may act as a timer for muscle differentiation.

The contribution of silencer variants to human diseases

BACKGROUND

Although disease-causal genetic variants have been found within silencer sequences, we still lack a comprehensive analysis of the association of silencers with diseases. Here, we profiled GWAS variants in 2.8 million candidate silencers across 97 human samples derived from a diverse panel of tissues and developmental time points, using deep learning models.

RESULTS

We show that candidate silencers exhibit strong enrichment in disease-associated variants, and several diseases display a much stronger association with silencer variants than enhancer variants. Close to 52% of candidate silencers cluster, forming silencer-rich loci, and, in the loci of Parkinson's-disease-hallmark genes TRIM31 and MAL, the associated SNPs densely populate clustered candidate silencers rather than enhancers displaying an overall twofold enrichment in silencers versus enhancers. The disruption of apoptosis in neuronal cells is associated with both schizophrenia and bipolar disorder and can largely be attributed to variants within candidate silencers. Our model permits a mechanistic explanation of causative SNP effects by identifying altered binding of tissue-specific repressors and activators, validated with a 70% of directional concordance using SNP-SELEX. Narrowing the focus of the analysis to individual silencer variants, experimental data confirms the role of the rs62055708 SNP in Parkinson's disease, rs2535629 in schizophrenia, and rs6207121 in type 1 diabetes.

CONCLUSIONS

In summary, our results indicate that advances in deep learning models for the discovery of disease-causal variants within candidate silencers effectively "double" the number of functionally characterized GWAS variants. This provides a basis for explaining mechanisms of action and designing novel diagnostics and therapeutics.

Hidden secrets of the cancer genome: unlocking the impact of non-coding mutations in gene regulatory elements

Discoveries in the field of genomics have revealed that non-coding genomic regions are not merely "junk DNA", but rather comprise critical elements involved in gene expression. These gene regulatory elements (GREs) include enhancers, insulators, silencers, and gene promoters. Notably, new evidence shows how mutations within these regions substantially influence gene expression programs, especially in the context of cancer. Advances in high-throughput sequencing technologies have accelerated the identification of somatic and germline single nucleotide mutations in non-coding genomic regions. This review provides an overview of somatic and germline non-coding single nucleotide alterations affecting transcription factor binding sites in GREs, specifically involved in cancer biology. It also summarizes the technologies available for exploring GREs and the challenges associated with studying and characterizing non-coding single nucleotide mutations. Understanding the role of GRE alterations in cancer is essential for improving diagnostic and prognostic capabilities in the precision medicine era, leading to enhanced patient-centered clinical outcomes.

Polycomb protein binding and looping in the ON transcriptional state

Polycomb group (PcG) proteins mediate epigenetic silencing of important developmental genes by modifying histones and compacting chromatin through two major protein complexes, PRC1 and PRC2. These complexes are recruited to DNA by CpG islands (CGIs) in mammals and Polycomb response elements (PREs) in Drosophila. When PcG target genes are turned OFF, PcG proteins bind to PREs or CGIs, and PREs serve as anchors that loop together and stabilize gene silencing. Here, we address which PcG proteins bind to PREs and whether PREs mediate looping when their targets are in the ON transcriptional state. While the binding of most PcG proteins decreases at PREs in the ON state, one PRC1 component, Ph, remains bound. Further, PREs can loop to each other and with presumptive enhancers in the ON state and, like CGIs, may act as tethering elements between promoters and enhancers. Overall, our data suggest that PREs are important looping elements for developmental loci in both the ON and OFF states.

Functional categorization of gene regulatory variants that cause Mendelian conditions

Much of our current understanding of rare human diseases is driven by coding genetic variants. However, non-coding genetic variants play a pivotal role in numerous rare human diseases, resulting in diverse functional impacts ranging from altered gene regulation, splicing, and/or transcript stability. With the increasing use of genome sequencing in clinical practice, it is paramount to have a clear framework for understanding how non-coding genetic variants cause disease. To this end, we have synthesized the literature on hundreds of non-coding genetic variants that cause rare Mendelian conditions via the disruption of gene regulatory patterns and propose a functional classification system. Specifically, we have adapted the functional classification framework used for coding variants (i.e., loss-of-function, gain-of-function, and dominant-negative) to account for features unique to non-coding gene regulatory variants. We identify that non-coding gene regulatory variants can be split into three distinct categories by functional impact: (1) non-modular loss-of-expression (LOE) variants; (2) modular loss-of-expression (mLOE) variants; and (3) gain-of-ectopic-expression (GOE) variants. Whereas LOE variants have a direct corollary with coding loss-of-function variants, mLOE and GOE variants represent disease mechanisms that are largely unique to non-coding variants. These functional classifications aim to provide a unified terminology for categorizing the functional impact of non-coding variants that disrupt gene regulatory patterns in Mendelian conditions.

Chromatin gene-gene loops support the cross-regulation of genes with related function

Chromatin loops between gene pairs have been observed in diverse contexts in both flies and vertebrates. Combining high-resolution Capture-C, DNA fluorescence in situ hybridization, and genetic perturbations, we dissect the functional role of three loops between genes with related function during Drosophila embryogenesis. By mutating the loop anchor (but not the gene) or the gene (but not loop anchor), we disentangle loop formation and gene expression and show that the 3D proximity of paralogous gene loci supports their co-regulation. Breaking the loop leads to either an attenuation or enhancement of expression and perturbs their relative levels of expression and cross-regulation. Although many loops appear constitutive across embryogenesis, their function can change in different developmental contexts. Taken together, our results indicate that chromatin gene-gene loops act as architectural scaffolds that can be used in different ways in different contexts to fine-tune the coordinated expression of genes with related functions and sustain their cross-regulation.

Polycomb protein binding and looping mediated by Polycomb Response Elements in the ON transcriptional state

Polycomb group proteins (PcG) mediate epigenetic silencing of important developmental genes and other targets. In Drosophila, canonical PcG-target genes contain Polycomb Response Elements (PREs) that recruit PcG protein complexes including PRC2 that trimethylates H3K27 forming large H3K27me3 domains. In the OFF transcriptional state, PREs loop with each other and this looping strengthens silencing. Here we address the question of what PcG proteins bind to PREs when canonical PcG target genes are expressed, and whether PREs loop when these genes are ON. Our data show that the answer to this question is PRE-specific but general conclusions can be made. First, within a PcG-target gene, some regulatory DNA can remain covered with H3K27me3 and PcG proteins remain bound to PREs in these regions. Second, when PREs are within H3K27ac domains, PcG-binding decreases, however, this depends on the protein and PRE. The DNA binding protein GAF, and the PcG protein Ph remain at PREs even when other PcG proteins are greatly depleted. In the ON state, PREs can still loop with each other, but also form loops with presumptive enhancers. These data support the model that, in addition to their role in PcG silencing, PREs can act as "promoter-tethering elements" mediating interactions between promoter proximal PREs and distant enhancers.

The activity of engrailed imaginal disc enhancers is modulated epigenetically by chromatin and autoregulation

engrailed (en) encodes a homeodomain transcription factor crucial for the proper development of Drosophila embryos and adults. Like many developmental transcription factors, en expression is regulated by many enhancers, some of overlapping function, that drive expression in spatially and temporally restricted patterns. The en embryonic enhancers are located in discrete DNA fragments that can function correctly in small reporter transgenes. In contrast, the en imaginal disc enhancers (IDEs) do not function correctly in small reporter transgenes. En is expressed in the posterior compartment of wing imaginal discs; in contrast, small IDE-reporter transgenes are expressed mainly in the anterior compartment. We found that En binds to the IDEs and suggest that it may directly repress IDE function and modulate En expression levels. We identified two en IDEs, O and S. Deletion of either of these IDEs from a 79kb HA-en rescue transgene (HAen79) caused a loss-of-function en phenotype when the HAen79 transgene was the sole source of En. In contrast, flies with a deletion of the same IDEs from an endogenous en gene had no phenotype, suggesting a resiliency not seen in the HAen79 rescue transgene. Inserting a gypsy insulator in HAen79 between en regulatory DNA and flanking sequences strengthened the activity of HAen79, giving better function in both the ON and OFF transcriptional states. Altogether our data suggest that the en IDEs stimulate expression in the entire imaginal disc, and that the ON/OFF state is set by epigenetic memory set by the embryonic enhancers. This epigenetic regulation is similar to that of the Ultrabithorax IDEs and we suggest that the activity of late-acting enhancers in other genes may be similarly regulated.

Context-dependent role of Pho binding sites in Polycomb complex recruitment in Drosophila

Polycomb group (PcG) proteins maintain the silenced state of key developmental genes, but how these proteins are recruited to specific regions of the genome is still not completely understood. In Drosophila, PcG proteins are recruited to Polycomb response elements (PREs) comprised of a flexible array of sites for sequence-specific DNA binding proteins, "PcG recruiters," including Pho, Spps, Cg, and GAF. Pho is thought to play a central role in PcG recruitment. Early data showed that mutation of Pho binding sites in PREs in transgenes abrogated the ability of those PREs to repress gene expression. In contrast, genome-wide experiments in pho mutants or by Pho knockdown showed that PcG proteins can bind to PREs in the absence of Pho. Here, we directly addressed the importance of Pho binding sites in 2 engrailed (en) PREs at the endogenous locus and in transgenes. Our results show that Pho binding sites are required for PRE activity in transgenes with a single PRE. In a transgene, 2 PREs together lead to stronger, more stable repression and confer some resistance to the loss of Pho binding sites. Making the same mutation in Pho binding sites has little effect on PcG-protein binding at the endogenous en gene. Overall, our data support the model that Pho is important for PcG binding but emphasize how multiple PREs and chromatin environment increase the ability of PREs to function in the absence of Pho. This supports the view that multiple mechanisms contribute to PcG recruitment in Drosophila.

Polycomb Recruiters Inside and Outside of the Repressed Domains

The establishment and stable inheritance of individual patterns of gene expression in different cell types are required for the development of multicellular organisms. The important epigenetic regulators are the Polycomb group (PcG) and Trithorax group (TrxG) proteins, which control the silenced and active states of genes, respectively. In Drosophila, the PcG/TrxG group proteins are recruited to the DNA regulatory sequences termed the Polycomb response elements (PREs). The PREs are composed of the binding sites for different DNA-binding proteins, the so-called PcG recruiters. Currently, the role of the PcG recruiters in the targeting of the PcG proteins to PREs is well documented. However, there are examples where the PcG recruiters are also implicated in the active transcription and in the TrxG function. In addition, there is increasing evidence that the genome-wide PcG recruiters interact with the chromatin outside of the PREs and overlap with the proteins of differing regulatory classes. Recent studies of the interactomes of the PcG recruiters significantly expanded our understanding that they have numerous interactors besides the PcG proteins and that their functions extend beyond the regulation of the PRE repressive activity. Here, we summarize current data about the functions of the PcG recruiters.

The activity of engrailed imaginal disc enhancers is modulated epigenetically by chromatin and autoregulation

engrailed (en) encodes a homeodomain transcription factor crucial for the proper development of Drosophila embryos and adults. Like many developmental transcription factors, en expression is regulated by many enhancers, some of overlapping function, that drive expression in spatially and temporally restricted patterns. The en embryonic enhancers are located in discrete DNA fragments that can function correctly in small reporter transgenes. In contrast, the en imaginal disc enhancers (IDEs) do not function correctly in small reporter transgenes. En is expressed in the posterior compartment of wing imaginal disks; small IDE-reporter transgenes are expressed in the anterior compartment, the opposite of what is expected. Our data show that the En protein binds to en IDEs, and we suggest that En directly represses IDE function. We identified two en IDEs, 'O' and 'S'. Deletion of either of these IDEs from a 79kb HA-en rescue transgene (HAen79) caused a loss-of-function en phenotype when the HAen79 transgene was the sole source of En. In contrast, flies with a deletion of the same IDEs from the endogenous en gene had no phenotype, suggesting a resiliency not seen in the HAen79 rescue transgene. Inserting a gypsy insulator in HAen79 between en regulatory DNA and flanking sequences strengthened the activity of HAen79, giving better function in both the ON and OFF transcriptional states. Altogether our data show that the en IDEs stimulate expression in the entire imaginal disc, and that the ON/OFF state is set by epigenetic regulators. Further, the endogenous locus imparts a stability to en function not seen even in a large transgene, reflecting the importance of both positive and negative epigenetic influences that act over relatively large distances in chromatin.

Gene silencing dynamics are modulated by transiently active regulatory elements

Transcriptional enhancers have been extensively characterized, but cis-regulatory elements involved in acute gene repression have received less attention. Transcription factor GATA1 promotes erythroid differentiation by activating and repressing distinct gene sets. Here, we study the mechanism by which GATA1 silences the proliferative gene Kit during murine erythroid cell maturation and define stages from initial loss of activation to heterochromatinization. We find that GATA1 inactivates a potent upstream enhancer but concomitantly creates a discrete intronic regulatory region marked by H3K27ac, short noncoding RNAs, and de novo chromatin looping. This enhancer-like element forms transiently and serves to delay Kit silencing. The element is ultimately erased via the FOG1/NuRD deacetylase complex, as revealed by the study of a disease-associated GATA1 variant. Hence, regulatory sites can be self-limiting by dynamic co-factor usage. Genome-wide analyses across cell types and species uncover transiently active elements at numerous genes during repression, suggesting that modulation of silencing kinetics is widespread.

Distinct roles for canonical and variant histone H3 lysine-36 in Polycomb silencing

Polycomb complexes regulate cell type-specific gene expression programs through heritable silencing of target genes. Trimethylation of histone H3 lysine 27 (H3K27me3) is essential for this process. Perturbation of H3K36 is thought to interfere with H3K27me3. We show that mutants of Drosophila replication-dependent (H3.2K36R) or replication-independent (H3.3K36R) histone H3 genes generally maintain Polycomb silencing and reach later stages of development. In contrast, combined (H3.3K36RH3.2K36R) mutants display widespread Hox gene misexpression and fail to develop past the first larval stage. Chromatin profiling revealed that the H3.2K36R mutation disrupts H3K27me3 levels broadly throughout silenced domains, whereas these regions are mostly unaffected in H3.3K36R animals. Analysis of H3.3 distributions showed that this histone is enriched at presumptive Polycomb response elements located outside of silenced domains but relatively depleted from those inside. We conclude that H3.2 and H3.3 K36 residues collaborate to repress Hox genes using different mechanisms.

Mechanisms of Interaction between Enhancers and Promoters in Three Drosophila Model Systems

In higher eukaryotes, the regulation of developmental gene expression is determined by enhancers, which are often located at a large distance from the promoters they regulate. Therefore, the architecture of chromosomes and the mechanisms that determine the functional interaction between enhancers and promoters are of decisive importance in the development of organisms. Mammals and the model animal Drosophila have homologous key architectural proteins and similar mechanisms in the organization of chromosome architecture. This review describes the current progress in understanding the mechanisms of the formation and regulation of long-range interactions between enhancers and promoters at three well-studied key regulatory loci in Drosophila.

Hox function and specificity – A tissue centric view

Since their discovery, the Hox genes, with their incredible power to reprogram the identity of complete body regions, a phenomenon called homeosis, have captured the fascination of many biologists. Recent research has provided new insights into the function of Hox proteins in different germ layers and the mechanisms they employ to control tissue morphogenesis. We focus in this review on the ectoderm and mesoderm to highlight new findings and discuss them with regards to established concepts of Hox target gene regulation. Furthermore, we highlight the molecular mechanisms involved the transcriptional repression of specific groups of Hox target genes, and summarize the role of Hox mediated gene silencing in tissue development. Finally, we reflect on recent findings identifying a large number of tissue-specific Hox interactor partners, which open up new avenues and directions towards a better understanding of Hox function and specificity in different tissues.

Single-cell multiomics reveals the complexity of TGFβ signalling to chromatin in iPSC-derived kidney organoids

TGFβ1 plays a regulatory role in the determination of renal cell fate and the progression of renal fibrosis. Here we show an association between SMAD3 and the histone methyltransferase, EZH2, during cell differentiation; ChIP-seq revealed that SMAD3 and EZH2 co-occupy the genome in iPSCs and in iPSC-derived nephron progenitors. Through integration of single cell gene expression and epigenome profiling, we identified de novo ACTA2^+ve/POSTN^+ve myofibroblasts in kidney organoids treated with TGFβ1, characterised by increased SMAD3-dependent cis chromatin accessibility and gene expression associated with fibroblast activation. We have identified fibrosis-associated regulons characterised by enrichment of SMAD3, AP1, the ETS family of transcription factors, and NUAK1, CREB3L1, and RARG, corresponding to enriched motifs at accessible loci identified by scATACseq. Treatment with the EZH2 specific inhibitor GSK343, blocked SMAD3-dependent cis co-accessibility and inhibited myofibroblast activation. This mechanism, through which TGFβ signals directly to chromatin, represents a critical determinant of fibrotic, differentiated states.

REDfly: An Integrated Knowledgebase for Insect Regulatory Genomics

We provide here an updated description of the REDfly (Regulatory Element Database for Fly) database of transcriptional regulatory elements, a unique resource that provides regulatory annotation for the genome of Drosophila and other insects. The genomic sequences regulating insect gene expression-transcriptional cis-regulatory modules (CRMs, e.g., "enhancers") and transcription factor binding sites (TFBSs)-are not currently curated by any other major database resources. However, knowledge of such sequences is important, as CRMs play critical roles with respect to disease as well as normal development, phenotypic variation, and evolution. Characterized CRMs also provide useful tools for both basic and applied research, including developing methods for insect control. REDfly, which is the most detailed existing platform for metazoan regulatory-element annotation, includes over 40,000 experimentally verified CRMs and TFBSs along with their DNA sequences, their associated genes, and the expression patterns they direct. Here, we briefly describe REDfly's contents and data model, with an emphasis on the new features implemented since 2020. We then provide an illustrated walk-through of several common REDfly search use cases.

Comparative interactome analysis of the PRE DNA-binding factors: purification of the Combgap-, Zeste-, Psq-, and Adf1-associated proteins

The Polycomb group (PcG) and Trithorax group (TrxG) proteins are key epigenetic regulators controlling the silenced and active states of genes in multicellular organisms, respectively. In Drosophila, PcG/TrxG proteins are recruited to the chromatin via binding to specific DNA sequences termed polycomb response elements (PREs). While precise mechanisms of the PcG/TrxG protein recruitment remain unknown, the important role is suggested to belong to sequence-specific DNA-binding factors. At the same time, it was demonstrated that the PRE DNA-binding proteins are not exclusively localized to PREs but can bind other DNA regulatory elements, including enhancers, promoters, and boundaries. To gain an insight into the PRE DNA-binding protein regulatory network, here, using ChIP-seq and immuno-affinity purification coupled to the high-throughput mass spectrometry, we searched for differences in abundance of the Combgap, Zeste, Psq, and Adf1 PRE DNA-binding proteins. While there were no conspicuous differences in co-localization of these proteins with other functional transcription factors, we show that Combgap and Zeste are more tightly associated with the Polycomb repressive complex 1 (PRC1), while Psq interacts strongly with the TrxG proteins, including the BAP SWI/SNF complex. The Adf1 interactome contained Mediator subunits as the top interactors. In addition, Combgap efficiently interacted with AGO2, NELF, and TFIID. Combgap, Psq, and Adf1 have architectural proteins in their networks. We further investigated the existence of direct interactions between different PRE DNA-binding proteins and demonstrated that Combgap-Adf1, Psq-Dsp1, and Pho-Spps can interact in the yeast two-hybrid assay. Overall, our data suggest that Combgap, Psq, Zeste, and Adf1 are associated with the protein complexes implicated in different regulatory activities and indicate their potential multifunctional role in the regulation of transcription.

Regulating specificity in enhancer-promoter communication

Enhancers are cis-regulatory elements that can activate transcription remotely to regulate a specific pattern of a gene's expression. Genes typically have many enhancers that are often intermingled in the loci of other genes. To regulate expression, enhancers must therefore activate their correct promoter while ignoring others that may be in closer linear proximity. In this review, we discuss mechanisms by which enhancers engage with promoters, including recent findings on the role of cohesin and the Mediator complex, and how this specificity in enhancer-promoter communication is encoded. Genetic dissection of model loci, in addition to more recent findings using genome-wide approaches, highlight the core promoter sequence, its accessibility, cofactor-promoter preference, in addition to the surrounding genomic context, as key components.

Enhancer-silencer transitions in the human genome

Dual-function regulatory elements (REs), acting as enhancers in some cellular contexts and as silencers in others, have been reported to facilitate the precise gene regulatory response to developmental signals in Drosophila melanogaster. However, with few isolated examples detected, dual-function REs in mammals have yet to be systematically studied. We herein investigated this class of REs in the human genome and profiled their activity across multiple cell types. Focusing on enhancer–silencer transitions specific to the development of T cells, we built an accurate deep learning classifier of REs and identified about 12,000 silencers active in primary peripheral blood T cells that act as enhancers in embryonic stem cells. Compared with regular silencers, these dual-function REs are evolving under stronger purifying selection and are enriched for mutations associated with disease phenotypes and altered gene expression. In addition, they are enriched in the loci of transcriptional regulators, such as transcription factors (TFs) and chromatin remodeling genes. Dual-function REs consist of two intertwined but largely distinct sets of binding sites bound by either activating or repressing TFs, depending on the type of RE function in a given cell line. This indicates the recruitment of different TFs for different regulatory modes and a complex DNA sequence composition of these REs with dual activating and repressive encoding. With an estimated >6% of cell type–specific human silencers acting as dual-function REs, this overlooked class of REs requires a specific investigation on how their inherent functional plasticity might be a contributing factor to human diseases.

An Evolutionary Perspective on Hox Binding Site Preferences in Two Different Tissues

Transcription factor (TF) networks define the precise development of multicellular organisms. While many studies focused on TFs expressed in specific cell types to elucidate their contribution to cell specification and differentiation, it is less understood how broadly expressed TFs perform their precise functions in the different cellular contexts. To uncover differences that could explain tissue-specific functions of such TFs, we analyzed here genomic chromatin interactions of the broadly expressed Drosophila Hox TF Ultrabithorax (Ubx) in the mesodermal and neuronal tissues using bioinformatics. Our investigations showed that Ubx preferentially interacts with multiple yet tissue-specific chromatin sites in putative regulatory regions of genes in both tissues. Importantly, we found the classical Hox/Ubx DNA binding motif to be enriched only among the neuronal Ubx chromatin interactions, whereas a novel Ubx-like motif with rather low predicted Hox affinities was identified among the regions bound by Ubx in the mesoderm. Finally, our analysis revealed that tissues-specific Ubx chromatin sites are also different with regards to the distribution of active and repressive histone marks. Based on our data, we propose that the tissue-related differences in Ubx binding behavior could be a result of the emergence of the mesoderm as a new germ layer in triploblastic animals, which might have required the Hox TFs to relax their binding specificity.

Distinct requirements for Pho, Sfmbt, and Ino80 for cell survival in Drosophila

The Drosophila proteins Pleiohomeotic (Pho) and its paralog Pho-like (Phol) are the homologs of the mammalian transcription factor YY1. Pho and Phol are subunits of the Polycomb group protein complex PhoRC and they are also stably associated with the INO80 nucleosome remodeling complex. Drosophila lacking both Pho and Phol arrest development as larvae with small misshaped imaginal discs. The basis of this phenotype is poorly understood. We find that in pho phol mutant animals cells retain the capacity to proliferate but show a high incidence of apoptotic cell death that results in tissue hypoplasia. Clonal analyses establish that cells stringently require Pho and Phol to survive. In contrast, the PhoRC subunit Sfmbt and the ATP-dependent nucleosome remodeling factor Ino80 are not essential for cell viability. Pho and Phol, therefore, execute their critical role for cell survival through mechanisms that do not involve Sfmbt function or INO80 nucleosome remodeling.

Annotating the Insect Regulatory Genome

An ever-growing number of insect genomes is being sequenced across the evolutionary spectrum. Comprehensive annotation of not only genes but also regulatory regions is critical for reaping the full benefits of this sequencing. Driven by developments in sequencing technologies and in both empirical and computational discovery strategies, the past few decades have witnessed dramatic progress in our ability to identify cis-regulatory modules (CRMs), sequences such as enhancers that play a major role in regulating transcription. Nevertheless, providing a timely and comprehensive regulatory annotation of newly sequenced insect genomes is an ongoing challenge. We review here the methods being used to identify CRMs in both model and non-model insect species, and focus on two tools that we have developed, REDfly and SCRMshaw. These resources can be paired together in a powerful combination to facilitate insect regulatory annotation over a broad range of species, with an accuracy equal to or better than that of other state-of-the-art methods.

Enhancer grammar in development, evolution, and disease: dependencies and interplay

Each language has standard books describing that language's grammatical rules. Biologists have searched for similar, albeit more complex, principles relating enhancer sequence to gene expression. Here, we review the literature on enhancer grammar. We introduce dependency grammar, a model where enhancers encode information based on dependencies between enhancer features shaped by mechanistic, evolutionary, and biological constraints. Classifying enhancers based on the types of dependencies may identify unifying principles relating enhancer sequence to gene expression. Such rules would allow us to read the instructions for development within genomes and pinpoint causal enhancer variants underlying disease and evolutionary changes.

H3K27me3-rich genomic regions can function as silencers to repress gene expression via chromatin interactions

The mechanisms underlying gene repression and silencers are poorly understood. Here we investigate the hypothesis that H3K27me3-rich regions of the genome, defined from clusters of H3K27me3 peaks, may be used to identify silencers that can regulate gene expression via proximity or looping. We find that H3K27me3-rich regions are associated with chromatin interactions and interact preferentially with each other. H3K27me3-rich regions component removal at interaction anchors by CRISPR leads to upregulation of interacting target genes, altered H3K27me3 and H3K27ac levels at interacting regions, and altered chromatin interactions. Chromatin interactions did not change at regions with high H3K27me3, but regions with low H3K27me3 and high H3K27ac levels showed changes in chromatin interactions. Cells with H3K27me3-rich regions knockout also show changes in phenotype associated with cell identity, and altered xenograft tumor growth. Finally, we observe that H3K27me3-rich regions-associated genes and long-range chromatin interactions are susceptible to H3K27me3 depletion. Our results characterize H3K27me3-rich regions and their mechanisms of functioning via looping.

Emerging Roles of Single-Cell Multi-Omics in Studying Developmental Temporal Patterning

The complexity of brain structure and function is rooted in the precise spatial and temporal regulation of selective developmental events. During neurogenesis, both vertebrates and invertebrates generate a wide variety of specialized cell types through the expansion and specification of a restricted set of neuronal progenitors. Temporal patterning of neural progenitors rests on fine regulation between cell-intrinsic and cell-extrinsic mechanisms. The rapid emergence of high-throughput single-cell technologies combined with elaborate computational analysis has started to provide us with unprecedented biological insights related to temporal patterning in the developing central nervous system (CNS). Here, we present an overview of recent advances in Drosophila and vertebrates, focusing both on cell-intrinsic mechanisms and environmental influences. We then describe the various multi-omics approaches that have strongly contributed to our current understanding and discuss perspectives on the various -omics approaches that hold great potential for the future of temporal patterning research.

Polycomb and Trithorax Group Proteins: The Long Road from Mutations in Drosophila to Use in Medicine

Polycomb group (PcG) and Trithorax group (TrxG) proteins are evolutionarily conserved factors responsible for the repression and activation of the transcription of multiple genes in Drosophila and mammals. Disruption of the PcG/TrxG expression is associated with many pathological conditions, including cancer, which makes them suitable targets for diagnosis and therapy in medicine. In this review, we focus on the major PcG and TrxG complexes, the mechanisms of PcG/TrxG action, and their recruitment to chromatin. We discuss the alterations associated with the dysfunction of a number of factors of these groups in oncology and the current strategies used to develop drugs based on small-molecule inhibitors.

Ecdysone-Induced 3D Chromatin Reorganization Involves Active Enhancers Bound by Pipsqueak and Polycomb

Evidence suggests that Polycomb (Pc) is present at chromatin loop anchors in Drosophila. Pc is recruited to DNA through interactions with the GAGA binding factors GAF and Pipsqueak (Psq). Using HiChIP in Drosophila cells, we find that the psq gene, which has diverse roles in development and tumorigenesis, encodes distinct isoforms with unanticipated roles in genome 3D architecture. The BR-C, ttk, and bab domain (BTB)-containing Psq isoform (Psq^L) colocalizes genome-wide with known architectural proteins. Conversely, Psq lacking the BTB domain (Psq^S) is consistently found at Pc loop anchors and at active enhancers, including those that respond to the hormone ecdysone. After stimulation by this hormone, chromatin 3D organization is altered to connect promoters and ecdysone-responsive enhancers bound by Psq^S. Our findings link Psq variants lacking the BTB domain to Pc-bound active enhancers, thus shedding light into their molecular function in chromatin changes underlying the response to hormone stimulus.

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.

Polycomb Repressive Complexes in Hox Gene Regulation: Silencing and Beyond: The Functional Dynamics of Polycomb Repressive Complexes in Hox Gene Regulation

The coordinated expression of the Hox gene family encoding transcription factors is critical for proper embryonic development and patterning. Major efforts have thus been dedicated to understanding mechanisms controlling Hox expression. In addition to the temporal and spatial sequential activation of Hox genes, proper embryonic development requires that Hox genes get differentially silenced in a cell-type specific manner as development proceeds. Factors contributing to Hox silencing include the polycomb repressive complexes (PRCs), which control gene expression through epigenetic modifications. This review focuses on PRC-dependent regulation of the Hox genes and is aimed at integrating the growing complexity of PRC functional properties in the context of Hox regulation. In particular, mechanisms underlying PRC binding dynamics as well as a series of studies that have revealed the impact of PRC on the 3D organization of the genome is discussed, which has a significant role on Hox regulation during development.

MafK Mediates Chromatin Remodeling to Silence IRF8 Expression in Non-immune Cells in a Cell Type-SpecificManner

The regulation of gene expression is a result of a complex interplay between chromatin remodeling, transcription factors, and signaling molecules. Cell differentiation is accompanied by chromatin remodeling of specific loci to permanently silence genes that are not essential for the differentiated cell activity. The molecular cues that recruit the chromatin remodeling machinery are not well characterized. IRF8 is an immune-cell specific transcription factor and its expression is augmented by interferon-γ. Therefore, it serves as a model gene to elucidate the molecular mechanisms governing its silencing in non-immune cells. Ahigh-throughput shRNA library screen in IRF8 expression-restrictive cells enabled the identification of MafK as modulator of IRF8 silencing, affecting chromatin architecture. ChIP-Seq analysis revealed three MafK binding regions (-25 kb, -20 kb, and IRF8 6th intron) within the IRF8 locus. These MafK binding sites are sufficient to repress a reporter gene when cloned in genome-integrated lentiviral reporter constructs in only expression-restrictive cells. Conversely, plasmid-based constructs do not demonstrate such repressive effect. These results highlight the role of these MafK binding sites in mediating repressed chromatin assembly. Finally, a more thorough genomic analysis was performed, using CRISPR-Cas9 to delete MafK-int6 binding region in IRF8 expression-restrictive cells. Deleted clones exhibited an accessible chromatin conformation within the IRF8 locus that was accompanied by a significant increase in basal expression of IRF8 that was further induced by interferon-γ. Taken together, we identified and characterized several MafK binding elements within the IRF8 locus that mediate repressive chromatin conformation resulting in the silencing of IRF8 expression in a celltype-specific manner.

Dynamic Competition of Polycomb and Trithorax in Transcriptional Programming

Predicting regulatory potential from primary DNA sequences or transcription factor binding patterns is not possible. However, the annotation of the genome by chromatin proteins, histone modifications, and differential compaction is largely sufficient to reveal the locations of genes and their differential activity states. The Polycomb Group (PcG) and Trithorax Group (TrxG) proteins are the central players in this cell type-specific chromatin organization. PcG function was originally viewed as being solely repressive and irreversible, as observed at the homeotic loci in flies and mammals. However, it is now clear that modular and reversible PcG function is essential at most developmental genes. Focusing mainly on recent advances, we review evidence for how PcG and TrxG patterns change dynamically during cell type transitions. The ability to implement cell type-specific transcriptional programming with exquisite fidelity is essential for normal development.

REDfly: the transcriptional regulatory element database for Drosophila

The REDfly database provides a comprehensive curation of experimentally-validated Drosophila transcriptional cis-regulatory elements and includes information on DNA sequence, experimental evidence, patterns of regulated gene expression, and more. Now in its thirteenth year, REDfly has grown to over 23 000 records of tested reporter gene constructs and 2200 tested transcription factor binding sites. Recent developments include the start of curation of predicted cis-regulatory modules in addition to experimentally-verified ones, improved search and filtering, and increased interaction with the authors of curated papers. An expanded data model that will capture information on temporal aspects of gene regulation, regulation in response to environmental and other non-developmental cues, sexually dimorphic gene regulation, and non-endogenous (ectopic) aspects of reporter gene expression is under development and expected to be in place within the coming year. REDfly is freely accessible at http://redfly.ccr.buffalo.edu, and news about database updates and new features can be followed on Twitter at @REDfly_database.

RNA-DNA strand exchange by the Drosophila Polycomb complex PRC2

Polycomb Group (PcG) proteins form memory of transient transcriptional repression that is necessary for development. In Drosophila, DNA elements termed Polycomb Response Elements (PREs) recruit PcG proteins. How PcG activities are targeted to PREs to maintain repressed states only in appropriate developmental contexts has been difficult to elucidate. PcG complexes modify chromatin, but also interact with both RNA and DNA, and RNA is implicated in PcG targeting and function. Here we show that R-loops form at many PREs in Drosophila embryos, and correlate with repressive states. In vitro, both PRC1 and PRC2 can recognize R-loops and open DNA bubbles. Unexpectedly, we find that PRC2 drives formation of RNA-DNA hybrids, the key component of R-loops, from RNA and dsDNA. Our results identify R-loop formation as a feature of Drosophila PREs that can be recognized by PcG complexes, and RNA-DNA strand exchange as a PRC2 activity that could contribute to R-loop formation.

Silencers, Enhancers, and the Multifunctional Regulatory Genome

Negative regulation of gene expression by transcriptional silencers has been difficult to study due to limited defined examples. A new study by Gisselbrecht et al. has dramatically increased the number of identified silencers and reveals that they are bifunctional regulatory sequences that also act as gene expression-promoting enhancers.

Genomic organization of Polycomb Response Elements and its functional implication in Drosophila and other insects

The epigenetic memory is an essential aspect of multicellular organisms to maintain several cell types and their gene expression pattern. This complex process uses a number of protein factors and specific DNA elements within the developmental cues to achieve this. The protein factors involved in the process are the Polycomb group (PcG) members, and, accordingly, the DNA sequences that interact with these proteins are called Polycomb Response Elements (PREs). Since the PcG proteins are highly conserved among higher eukaryotes, including insects, and function at thousands of sites in the genomes, it is expected that PREs mayalso be present across the genome.However, the studies on PREs in insect species, other thanDrosophila, is currently lacking.We took a bioinformatics approach to develop an inclusive PRE prediction tool, 'PRE Mapper', to address this need. By applying this tool on the Drosophila melanogaster genome, we predicted greater than 20,000 PREs.When comparedwith the available PRE prediction methods, this tool shows far better performance by correctly identifying the in vivo binding sites of PcG proteins, identified by genome-scale ChIP experiments. Further analysis of the predicted PREs shows their cohabitation with chromatin domain boundary elements at several places in the Drosophila genome, possibly defining a composite epigenetic module.We analysed 10 insect genomes in this context and find several conserved features in PREs across the insect species with some variations in their occurrence frequency. These analyses leading to the identification of PREin insect genomes contribute to our understanding of epigenetic mechanisms in these organisms.

Widespread activation of developmental gene expression characterized by PRC1-dependent chromatin looping

Polycomb repressive complexes 1 and 2 have been historically described as transcriptional repressors, but recent reports suggest that PRC1 might also support activation, although the underlying mechanisms remain elusive. Here, we show that stage-specific PRC1 binding at a subset of active promoters and enhancers during Drosophila development coincides with the formation of three-dimensional (3D) loops, an increase in expression during development and repression in PRC1 mutants. Dissection of the dachshund locus indicates that PRC1-anchored loops are versatile architectural platforms that persist when surrounding genes are transcriptionally active and fine-tune their expression. The analysis of RING1B binding profiles and 3D contacts during neural differentiation in mice suggests that this role is conserved in mammals.

DNA sequence models of genome-wide Drosophila melanogaster Polycomb binding sites improve generalization to independent Polycomb Response Elements

Polycomb Response Elements (PREs) are cis-regulatory DNA elements that maintain gene transcription states through DNA replication and mitosis. PREs have little sequence similarity, but are enriched in a number of sequence motifs. Previous methods for modelling Drosophila melanogaster PRE sequences (PREdictor and EpiPredictor) have used a set of 7 motifs and a training set of 12 PREs and 16-23 non-PREs. Advances in experimental methods for mapping chromatin binding factors and modifications has led to the publication of several genome-wide sets of Polycomb targets. In addition to the seven motifs previously used, PREs are enriched in the GTGT motif, recently associated with the sequence-specific DNA binding protein Combgap. We investigated whether models trained on genome-wide Polycomb sites generalize to independent PREs when trained with control sequences generated by naive PRE models and including the GTGT motif. We also developed a new PRE predictor: SVM-MOCCA. Training PRE predictors with genome-wide experimental data improves generalization to independent data, and SVM-MOCCA predicts the majority of PREs in three independent experimental sets. We present 2908 candidate PREs enriched in sequence and chromatin signatures. 2412 of these are also enriched in H3K4me1, a mark of Trithorax activated chromatin, suggesting that PREs/TREs have a common sequence code.

Transcriptional Silencers in Drosophila Serve a Dual Role as Transcriptional Enhancers in Alternate Cellular Contexts

A major challenge in biology is to understand how complex gene expression patterns are encoded in the genome. While transcriptional enhancers have been studied extensively, few transcriptional silencers have been identified, and they remain poorly understood. Here, we used a novel strategy to screen hundreds of sequences for tissue-specific silencer activity in whole Drosophila embryos. Almost all of the transcriptional silencers that we identified were also active enhancers in other cellular contexts. These elements are bound by more transcription factors than non-silencers. A subset of these silencers forms long-range contacts with promoters. Deletion of a silencer caused derepression of its target gene. Our results challenge the common practice of treating enhancers and silencers as separate classes of regulatory elements and suggest the possibility that thousands or more bifunctional CRMs remain to be discovered in Drosophila and 10⁴-10⁵ in humans.

The multiscale effects of polycomb mechanisms on 3D chromatin folding

Polycomb group (PcG) proteins silence master regulatory genes required to properly confer cell identity during the development of both Drosophila and mammals. They may act through chromatin compaction and higher-order folding of chromatin inside the cell nucleus. During the last decade, analysis on interphase chromosome architecture discovered self-interacting regions named topologically associated domains (TADs). TADs result from the 3D chromatin folding of a succession of transcribed and repressed epigenomic domains and from loop extrusion mediated by cohesin/CTCF in mammals. Polycomb silenced chromatin constitutes one type of repressed epigenomic domains which form compacted nano-compartments inside cell nuclei. Recruitment of canonical PcG proteins on chromatin relies on initial binding to discrete elements and further spreading into large chromatin domains covered with H3K27me3. Some of these discrete elements have a bivalent nature both in mammals and Drosophila and are dynamically regulated during development. Loops can occur between them, suggesting that their interaction plays both functional and structural roles. Formation of large chromatin domains covered by H3K27me3 seems crucial for PcG silencing and PcG proteins might exert their function through compaction of these domains in both mammals and flies, rather than by directly controlling the nucleosomal accessibility of discrete regulatory elements. In addition, PcG chromatin domains interact over long genomic distances, shaping a higher-order chromatin network. Therefore, PcG silencing might rely on multiscale chromatin folding to maintain cell identity during differentiation.

Two-step regulation of trachealess ensures tight coupling of cell fate with morphogenesis in the Drosophila trachea

During organogenesis, inductive signals cause cell differentiation and morphogenesis. However, how these phenomena are coordinated to form functional organs is poorly understood. Here, we show that cell differentiation of the Drosophila trachea is sequentially determined in two steps and that the second step is synchronous with the invagination of the epithelial sheet. The master gene trachealess is dispensable for the initiation of invagination, while it is essential for maintaining the invaginated structure, suggesting that tracheal morphogenesis and differentiation are separately induced. trachealess expression starts in bipotential tracheal/epidermal placode cells. After invagination, its expression is maintained in the invaginated cells but is extinguished in the remaining sheet cells. A trachealess cis-regulatory module that shows both tracheal enhancer activity and silencer activity in the surface epidermal sheet was identified. We propose that the coupling of trachealess expression with the invaginated structure ensures that only invaginated cells canalize robustly into the tracheal fate.

Cells in developing organs have two important decisions to make: where to be and what cell type to become. If cells end up in the wrong places, they can stop an organ from working, so it is vital that one decision depends upon the other. The so-called progenitor cells responsible for forming the trachea, for example, can either become part of a flat sheet or part of a tube. The cells on the sheet need to become epidermal cells, while the cells in the tube need to become tracheal cells. Work on fruit flies found that a gene called 'trachealess' plays an important role in this process. Without it, developing flies cannot make a trachea at all.

At the start of trachea development, some of the cells form thickened structures called placodes. The progenitor cells in the placodes start to divide, and the structures buckle inwards to form pockets. These pockets then lengthen into tubes. The trachealess gene codes for a protein that works as a genetic switch. It turns other genes on or off, helping the progenitor cells inside the pockets to become tracheal cells. But, it is not clear whether trachealess drives the formation of the pockets: the progenitor cells first decide what to be; or whether pocket formation tells the cells to use trachealess: the progenitor cells first decide where to be.

To find out, Kondo and Hayashi imaged developing fly embryos and saw that the trachealess gene does not start pocket formation, but that it is essential to maintain the pockets. Flies without the gene managed to form pockets, but they did not last long. Looking at embryos with defects in other genes involved in pocket formation revealed why. In these flies, some of the progenitor cells using trachealess got left behind when the pockets started to form. But rather than forming pockets of their own (as they might if trachealess were driving pocket formation), they turned their trachealess gene off. Progenitor cells in the fly trachea seem to decide where to be before they decide what cell type to become. This helps to make sure that trachea cells do not form in the wrong places.

A question that still remains is how do the cells know when they are inside a pocket? It is possible that the cells are sensing different mechanical forces or different chemical signals. Further research could help scientists to understand how organs form in living animals, and how they might better recreate that process in the laboratory.

From insects to mammals: regulation of genome architecture in neural development

One of the hallmarks of the metazoan genome is that genes are non-randomly positioned within the cell nucleus; in fact, the entire genome is packaged in a highly organized manner to orchestrate proper gene expression for each cell type. This is an especially daunting task for the development of the brain, which consists of an incredibly diverse population of neural cells. How genome architecture is established, maintained, and regulated to promote diverse cell fates and functions are fascinating questions with important implications in development and disease. The explosion in various biochemical and imaging techniques to analyze chromatin is now making it possible to interrogate the genome at an unprecedented resolution. Here we will focus on current advances in understanding genome architecture and gene regulation in the context of neural development.

Polycomb Assemblies Multitask to Regulate Transcription

The Polycomb system is made of an evolutionary ancient group of proteins, present throughout plants and animals. Known initially from developmental studies with the fly Drosophila melanogaster, they were associated with stable sustainment of gene repression and maintenance of cell identity. Acting as multiprotein assemblies with an ability to modify chromatin, through chemical additions to histones and organization of topological domains, they have been involved subsequently in control of developmental transitions and in cell homeostasis. Recent work has unveiled an association of Polycomb components with transcriptionally active loci and the promotion of gene expression, in clear contrast with conventional recognition as repressors. Focusing on mammalian models, I review here advances concerning roles in transcriptional control. Among new findings highlighted is the regulation of their catalytic properties, recruiting to targets, and activities in chromatin organization and compartmentalization. The need for a more integrated approach to the study of the Polycomb system, given its fundamental complexity and its adaptation to cell context, is discussed.

Dangerous liaisons: interplay between SWI/SNF, NuRD, and Polycomb in chromatin regulation and cancer

In this review, Bracken et al. discuss the functional organization and biochemical activities of remodelers and Polycomb and explore how they work together to control cell differentiation and the maintenance of cell identity. They also discuss how mutations in the genes encoding these various chromatin regulators contribute to oncogenesis by disrupting the chromatin equilibrium.

Changes in chromatin structure mediated by ATP-dependent nucleosome remodelers and histone modifying enzymes are integral to the process of gene regulation. Here, we review the roles of the SWI/SNF (switch/sucrose nonfermenting) and NuRD (nucleosome remodeling and deacetylase) and the Polycomb system in chromatin regulation and cancer. First, we discuss the basic molecular mechanism of nucleosome remodeling, and how this controls gene transcription. Next, we provide an overview of the functional organization and biochemical activities of SWI/SNF, NuRD, and Polycomb complexes. We describe how, in metazoans, the balance of these activities is central to the proper regulation of gene expression and cellular identity during development. Whereas SWI/SNF counteracts Polycomb, NuRD facilitates Polycomb repression on chromatin. Finally, we discuss how disruptions of this regulatory equilibrium contribute to oncogenesis, and how new insights into the biological functions of remodelers and Polycombs are opening avenues for therapeutic interventions on a broad range of cancer types.

Two modes of transvection at the eyes absent gene of Drosophila demonstrate plasticity in transcriptional regulatory interactions in cis and in trans

For many genes, proper gene expression requires coordinated and dynamic interactions between multiple regulatory elements, each of which can either promote or silence transcription. In Drosophila, the complexity of the regulatory landscape is further complicated by the tight physical pairing of homologous chromosomes, which can permit regulatory elements to interact in trans, a phenomenon known as transvection. To better understand how gene expression can be programmed through cis- and trans-regulatory interactions, we analyzed transvection effects for a collection of alleles of the eyes absent (eya) gene. We find that trans-activation of a promoter by the eya eye-specific enhancers is broadly supported in many allelic backgrounds, and that the availability of an enhancer to act in trans can be predicted based on the molecular lesion of an eya allele. Furthermore, by manipulating promoter availability in cis and in trans, we demonstrate that the eye-specific enhancers of eya show plasticity in their promoter preference between two different transcriptional start sites, which depends on promoter competition between the two potential targets. Finally, we show that certain alleles of eya demonstrate pairing-sensitive silencing resulting from trans-interactions between Polycomb Response Elements (PREs), and genetic and genomic data support a general role for PcG proteins in mediating transcriptional silencing at eya. Overall, our data highlight how eya gene regulation relies upon a complex but plastic interplay between multiple enhancers, promoters, and PREs.

Gene regulation requires interactions between regions of DNA known as regulatory elements, which, in combination, determine where and when a gene will be active or silenced. Some genes use just a few regulatory elements, whereas others rely on highly complex interactions between many different elements that are poorly understood. While we typically imagine regulatory elements interacting with one another along the length of a single chromosome, in a curious phenomenon called transvection, elements can communicate between two different chromosomes that are held in close proximity. Here, we use the study of transvection to better understand how different regulatory elements contribute to the expression of eyes absent (eya), a gene required for proper eye development in Drosophila. Our data show that a class of elements that initiate eya gene expression, called promoters, will compete with one another for activation by eya’s enhancers, a second class of regulatory element, with the promoter that is closest to the enhancers being the favored target for activation. Furthermore, our study of transvection uncovers an important role for a silencing element, called a PRE, in opposing eya gene expression. Overall, our study sheds new light on how different elements combine to produce patterned expression of eya.

The Hox transcription factor Ubx stabilizes lineage commitment by suppressing cellular plasticity in Drosophila

During development cells become restricted in their differentiation potential by repressing alternative cell fates, and the Polycomb complex plays a crucial role in this process. However, how alternative fate genes are lineage-specifically silenced is unclear. We studied Ultrabithorax (Ubx), a multi-lineage transcription factor of the Hox class, in two tissue lineages using sorted nuclei and interfered with Ubx in mesodermal cells. We find that depletion of Ubx leads to the de-repression of genes normally expressed in other lineages. Ubx silences expression of alternative fate genes by retaining the Polycomb Group protein Pleiohomeotic at Ubx targeted genomic regions, thereby stabilizing repressive chromatin marks in a lineage-dependent manner. Our study demonstrates that Ubx stabilizes lineage choice by suppressing the multipotency encoded in the genome via its interaction with Pho. This mechanism may explain why the Hox code is maintained throughout the lifecycle, since it could set a block to transdifferentiation in adult cells.