Enhancer-promoter communication at the Drosophila engrailed locus

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

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

The Blueprint of Logical Decisions in a NF-κB Signaling System

Nearly identical cells can exhibit substantially different responses to the same stimulus that causes phenotype diversity. Such interplay between phenotype diversity and the architecture of regulatory circuits is crucial since it determines the state of a biological cell. Here, we theoretically analyze how the circuit blueprints of NF-κB in cellular environments are formed and their role in determining the cells' metabolic state. The NF-κB is a collective name for a developmental conserved family of five different transcription factors that can form homodimers or heterodimers and often promote DNA looping to reprogram the inflammatory gene response. The NF-κB controls many biological functions, including cellular differentiation, proliferation, migration, and survival. Our model shows that nuclear localization of NF-κB differentially promotes logic operations such as AND, NAND, NOR, and OR in its regulatory network. Through the quantitative thermodynamic model of transcriptional regulation and systematic variation of promoter-enhancer interaction modes, we can account for the origin of various logic gates as formed in the NF-κB system. We further show that the interconversion or switching of logic gates yielded under systematic variations of the stimuli activity and DNA looping parameters. Such computation occurs in regulatory and signaling pathways in individual cells at a molecular scale, which one can exploit to design a biomolecular computer.

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.

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 Landscape of the DNA Transposons in the Genome of the Horezu_LaPeri Strain of Drosophila melanogaster

Natural transposons (NTs) represent mobile DNA sequences found in both prokaryotic and eukaryotic genomes. Drosophila melanogaster (the fruit fly) is a eukaryotic model organism with NTs standing for about 20% of its genome and has contributed significantly to the understanding of various aspects of transposon biology. Our study describes an accurate approach designed to map class II transposons (DNA transposons) in the genome of the Horezu_LaPeri fruit fly strain, consecutive to Oxford Nanopore Technology sequencing. A whole genome bioinformatics analysis was conducted using Genome ARTIST_v2, LoRTE and RepeatMasker tools to identify DNA transposons insertions. Then, a gene ontology enrichment analysis was performed in order to evaluate the potential adaptive role of some DNA transposons insertions. Herein, we describe DNA transposon insertions specific for the Horezu_LaPeri genome and a predictive functional analysis of some insertional alleles. The PCR validation of P-element insertions specific for this fruit fly strain, along with a putative consensus sequence for the KP element, is also reported. Overall, the genome of the Horezu_LaPeri strain contains several insertions of DNA transposons associated with genes known to be involved in adaptive processes. For some of these genes, insertional alleles obtained via mobilization of the artificial transposons were previously reported. This is a very alluring aspect, as it suggests that insertional mutagenesis experiments conducting adaptive predictions for laboratory strains may be confirmed by mirroring insertions which are expected to be found at least in some natural fruit fly strains.

Altering enhancer-promoter linear distance impacts promoter competition in cis and in trans

In Drosophila, pairing of maternal and paternal homologs can permit trans-interactions between enhancers on one homolog and promoters on another, an example of a phenomenon called transvection. When chromosomes are paired, promoters in cis and in trans to an enhancer can compete for the enhancer's activity, but the parameters that govern this competition are as yet poorly understood. To assess how the linear spacing between an enhancer and promoter can influence promoter competition in Drosophila, we employed transgenic constructs wherein the eye-specific enhancer GMR is placed at varying distances from a heterologous hsp70 promoter driving a fluorescent reporter. While GMR activates the reporter to a high degree when the enhancer and promoter are spaced by a few hundred base pairs, activation is strongly attenuated when the enhancer is moved 3 kilobases away. By examining transcription of endogenous genes near the point of transgene insertion, we show that linear spacing of 3 kb between GMR and the hsp70 promoter results in elevated transcription of neighboring promoters, suggesting a loss of specificity between the enhancer and its intended transgenic target promoter. Furthermore, increasing spacing between GMR and hsp70 by just 100 bp can enhance transvection, resulting in increased activation of a promoter on a paired homolog at the expense of a promoter in cis to the enhancer. Finally, cis-/trans-promoter competition assays in which one promoter carries mutations to key core promoter elements show that GMR will skew its activity toward a wild type promoter, suggesting that an enhancer is in a balanced competition between its potential target promoters in cis and in trans.

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.

Subunits of the PBAP Chromatin Remodeler Are Capable of Mediating Enhancer-Driven Transcription in Drosophila

The chromatin remodeler SWI/SNF is an important participant in gene activation, functioning predominantly by opening the chromatin structure on promoters and enhancers. Here, we describe its novel mode of action in which SWI/SNF factors mediate the targeted action of an enhancer. We studied the functions of two signature subunits of PBAP subfamily, BAP170 and SAYP, in Drosophila. These subunits were stably tethered to a transgene reporter carrying the hsp70 core promoter. The tethered subunits mediate transcription of the reporter in a pattern that is generated by enhancers close to the insertion site in multiple loci throughout the genome. Both tethered SAYP and BAP170 recruit the whole PBAP complex to the reporter promoter. However, we found that BAP170-dependent transcription is more resistant to the depletion of other PBAP subunits, suggesting that BAP170 may play a more critical role in establishing enhancer-dependent transcription.

Defining the Boundaries of Polycomb Domains in Drosophila

Polycomb group (PcG) proteins are an important group of transcriptional repressors that act by modifying chromatin. PcG target genes are covered by the repressive chromatin mark H3K27me3. Polycomb repressive complex 2 (PRC2) is a multiprotein complex that is responsible for generating H3K27me3. In Drosophila, PRC2 is recruited by Polycomb Response Elements (PREs) and then trimethylates flanking nucleosomes, spreading the H3K27me3 mark over large regions of the genome, the "Polycomb domains." What defines the boundary of a Polycomb domain? There is experimental evidence that insulators, PolII, and active transcription can all form the boundaries of Polycomb domains. Here we divide the boundaries of larval Polycomb domains into six different categories. In one category, genes are transcribed toward the Polycomb domain, where active transcription is thought to stop the spreading of H3K27me3. In agreement with this, we show that introducing a transcriptional terminator into such a transcription unit causes an extension of the Polycomb domain. Additional data suggest that active transcription of a boundary gene may restrict the range of enhancer activity of a Polycomb-regulated gene.

A theoretical model of Polycomb/Trithorax action unites stable epigenetic memory and dynamic regulation

Polycomb and Trithorax group proteins maintain stable epigenetic memory of gene expression states for some genes, but many targets show highly dynamic regulation. Here we combine experiment and theory to examine the mechanistic basis of these different modes of regulation. We present a mathematical model comprising a Polycomb/Trithorax response element (PRE/TRE) coupled to a promoter and including Drosophila developmental timing. The model accurately recapitulates published studies of PRE/TRE mediated epigenetic memory of both silencing and activation. With minimal parameter changes, the same model can also recapitulate experimental data for a different PRE/TRE that allows dynamic regulation of its target gene. The model predicts that both cell cycle length and PRE/TRE identity are critical for determining whether the system gives stable memory or dynamic regulation. Our work provides a simple unifying framework for a rich repertoire of PRE/TRE functions, and thus provides insights into  genome-wide Polycomb/Trithorax regulation.

Polycomb (PcG) and Trithorax (TrxG) group regulate several hundred target genes with important roles in development and disease. Here the authors combine experiment and theory to provide evidence that the Polycomb/Trithorax system has the potential for a rich repertoire of regulatory modes beyond simple epigenetic memory.

Transcriptional Enhancers in Drosophila

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

Bicoid-Dependent Activation of the Target Gene hunchback Requires a Two-Motif Sequence Code in a Specific Basal Promoter

In complex genetic loci, individual enhancers interact most often with specific basal promoters. Here we investigate the activation of the Bicoid target gene hunchback (hb), which contains two basal promoters (P1 and P2). Early in embryogenesis, P1 is silent, while P2 is strongly activated. In vivo deletion of P2 does not cause activation of P1, suggesting that P2 contains intrinsic sequence motifs required for activation. We show that a two-motif code (a Zelda binding site plus TATA) is required and sufficient for P2 activation. Zelda sites are present in the promoters of many embryonically expressed genes, but the combination of Zelda plus TATA does not seem to be a general code for early activation or Bicoid-specific activation per se. Because Zelda sites are also found in Bicoid-dependent enhancers, we propose that simultaneous binding to both enhancers and promoters independently synchronizes chromatin accessibility and facilitates correct enhancer-promoter interactions.

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.

Identification of evolutionarily conserved downstream core promoter elements required for the transcriptional regulation of Fushi tarazu target genes

The regulation of transcription initiation is critical for developmental and cellular processes. RNA polymerase II (Pol II) is recruited by the basal transcription machinery to the core promoter where Pol II initiates transcription. The core promoter encompasses the region from -40 to +40 bp relative to the +1 transcription start site (TSS). Core promoters may contain one or more core promoter motifs that confer specific properties to the core promoter, such as the TATA box, initiator (Inr) and motifs that are located downstream of the TSS, namely, motif 10 element (MTE), the downstream core promoter element (DPE) and the Bridge, a bipartite core promoter element. We had previously shown that Caudal, an enhancer-binding homeodomain transcription factor and a key regulator of the Hox gene network, is a DPE-specific activator. Interestingly, pair-rule proteins have been implicated in enhancer-promoter communication at the engrailed locus. Fushi tarazu (Ftz) is an enhancer-binding homeodomain transcription factor encoded by the ftz pair-rule gene. Ftz works in concert with its co-factor, Ftz-F1, to activate transcription. Here, we examined whether Ftz and Ftz-F1 activate transcription with a preference for a specific core promoter motif. Our analysis revealed that similarly to Caudal, Ftz and Ftz-F1 activate the promoter containing a TATA box mutation to significantly higher levels than the promoter containing a DPE mutation, thus demonstrating a preference for the DPE motif. We further discovered that Ftz target genes are enriched for a combination of functional downstream core promoter elements that are conserved among Drosophila species. Thus, the unique combination (Inr, Bridge and DPE) of functional downstream core promoter elements within Ftz target genes highlights the complexity of transcriptional regulation via the core promoter in the transcription of different developmental gene regulatory networks.

Structure and function of an ectopic Polycomb chromatin domain

Polycomb group proteins (PcGs) drive target gene repression and form large chromatin domains. In Drosophila, DNA elements known as Polycomb group response elements (PREs) recruit PcGs to the DNA. We have shown that, within the invected-engrailed (inv-en) Polycomb domain, strong, constitutive PREs are dispensable for Polycomb domain structure and function. We suggest that the endogenous chromosomal location imparts stability to this Polycomb domain. To test this possibility, a 79-kb en transgene was inserted into other chromosomal locations. This transgene is functional and forms a Polycomb domain. The spreading of the H3K27me3 repressive mark, characteristic of PcG domains, varies depending on the chromatin context of the transgene. Unlike at the endogenous locus, deletion of the strong, constitutive PREs from the transgene leads to both loss- and gain-of function phenotypes, demonstrating the important role of these regulatory elements. Our data show that chromatin context plays an important role in Polycomb domain structure and function.

Programmable DNA looping using engineered bivalent dCas9 complexes

DNA looping is a ubiquitous and critical feature of gene regulation. Although DNA looping can be efficiently detected, tools to readily manipulate DNA looping are limited. Here we develop CRISPR-based DNA looping reagents for creation of programmable DNA loops. Cleavage-defective Cas9 proteins of different specificity are linked by heterodimerization or translational fusion to create bivalent complexes able to link two separate DNA regions. After model-directed optimization, the reagents are validated using a quantitative DNA looping assay in E. coli. Looping efficiency is ~15% for a 4.7 kb loop, but is significantly improved by loop multiplexing with additional guides. Bivalent dCas9 complexes are also used to activate endogenous norVW genes by rewiring chromosomal DNA to bring distal enhancer elements to the gene promoters. Such reagents should allow manipulation of DNA looping in a variety of cell types, aiding understanding of endogenous loops and enabling creation of new regulatory connections.

DNA loops are a ubiquitious feature of gene regulation across the kingdoms of life. Here the authors design a Cas9-based dimerization system for inducing DNA loops in E. coli, allowing activation and rewiring of gene expression.

The free-energy cost of interaction between DNA loops

From the viewpoint of thermodynamics, the formation of DNA loops and the interaction between them, which are all non-equilibrium processes, result in the change of free energy, affecting gene expression and further cell-to-cell variability as observed experimentally. However, how these processes dissipate free energy remains largely unclear. Here, by analyzing a mechanic model that maps three fundamental topologies of two interacting DNA loops into a 4-state model of gene transcription, we first show that a longer DNA loop needs more mean free energy consumption. Then, independent of the type of interacting two DNA loops (nested, side-by-side or alternating), the promotion between them always consumes less mean free energy whereas the suppression dissipates more mean free energy. More interestingly, we find that in contrast to the mechanism of direct looping between promoter and enhancer, the facilitated-tracking mechanism dissipates less mean free energy but enhances the mean mRNA expression, justifying the facilitated-tracking hypothesis, a long-standing debate in biology. Based on minimal energy principle, we thus speculate that organisms would utilize the mechanisms of loop-loop promotion and facilitated tracking to survive in complex environments. Our studies provide insights into the understanding of gene expression regulation mechanism from the view of energy consumption.

Efficient chromosomal-scale DNA looping in Escherichia coli using multiple DNA-looping elements

Genes are frequently regulated by interactions between proteins that bind to the DNA near the gene and proteins that bind to DNA sites located far away, with the intervening DNA looped out. But it is not understood how efficient looping can occur when the sites are very far apart. We develop a simple theoretical framework that relates looping efficiency to the energetic cost and benefit of looping, allowing prediction of the efficiency of single or multiple nested loops at different distances. Measurements of absolute loop efficiencies for Lac repressor and λ CI using gene expression reporters in Escherichia coli cells show that, as predicted by the model, long-range DNA looping between a pair of sites can be strongly enhanced by the use of nested DNA loops or by the use of additional protein-binding sequences. A combination of these approaches was able to generate efficient DNA looping at a 200 kb distance.

Enhancer Sharing Promotes Neighborhoods of Transcriptional Regulation Across Eukaryotes

Enhancers physically interact with transcriptional promoters, looping over distances that can span multiple regulatory elements. Given that enhancer-promoter (EP) interactions generally occur via common protein complexes, it is unclear whether EP pairing is predominantly deterministic or proximity guided. Here, we present cross-organismic evidence suggesting that most EP pairs are compatible, largely determined by physical proximity rather than specific interactions. By reanalyzing transcriptome datasets, we find that the transcription of gene neighbors is correlated over distances that scale with genome size. We experimentally show that nonspecific EP interactions can explain such correlation, and that EP distance acts as a scaling factor for the transcriptional influence of an enhancer. We propose that enhancer sharing is commonplace among eukaryotes, and that EP distance is an important layer of information in gene regulation.

The Shh Topological Domain Facilitates the Action of Remote Enhancers by Reducing the Effects of Genomic Distances

Gene expression often requires interaction between promoters and distant enhancers, which occur within the context of highly organized topologically associating domains (TADs). Using a series of engineered chromosomal rearrangements at the Shh locus, we carried out an extensive fine-scale characterization of the factors that govern the long-range regulatory interactions controlling Shh expression. We show that Shh enhancers act pervasively, yet not uniformly, throughout the TAD. Importantly, changing intra-TAD distances had no impact on Shh expression. In contrast, inversions disrupting the TAD altered global folding of the region and prevented regulatory contacts in a distance-dependent manner. Our data indicate that the Shh TAD promotes distance-independent contacts between distant regions that would otherwise interact only sporadically, enabling functional communication between them. In large genomes where genomic distances per se can limit regulatory interactions, this function of TADs could be as essential for gene expression as the formation of insulated neighborhoods.

Formation of a Polycomb-Domain in the Absence of Strong Polycomb Response Elements

Polycomb group response elements (PREs) in Drosophila are DNA-elements that recruit Polycomb proteins (PcG) to chromatin and regulate gene expression. PREs are easily recognizable in the Drosophila genome as strong peaks of PcG-protein binding over discrete DNA fragments; many small but statistically significant PcG peaks are also observed in PcG domains. Surprisingly, in vivo deletion of the four characterized strong PREs from the PcG regulated invected-engrailed (inv-en) gene complex did not disrupt the formation of the H3K27me3 domain and did not affect inv-en expression in embryos or larvae suggesting the presence of redundant PcG recruitment mechanism. Further, the 3D-structure of the inv-en domain was only minimally altered by the deletion of the strong PREs. A reporter construct containing a 7.5kb en fragment that contains three weak peaks but no large PcG peaks forms an H3K27me3 domain and is PcG-regulated. Our data suggests a model for the recruitment of PcG-complexes to Drosophila genes via interactions with multiple, weak PREs spread throughout an H3K27me3 domain.

Functional Dissection of the Blocking and Bypass Activities of the Fab-8 Boundary in the Drosophila Bithorax Complex

Functionally autonomous regulatory domains direct the parasegment-specific expression of the Drosophila Bithorax complex (BX-C) homeotic genes. Autonomy is conferred by boundary/insulator elements that separate each regulatory domain from its neighbors. For six of the nine parasegment (PS) regulatory domains in the complex, at least one boundary is located between the domain and its target homeotic gene. Consequently, BX-C boundaries must not only block adventitious interactions between neighboring regulatory domains, but also be permissive (bypass) for regulatory interactions between the domains and their gene targets. To elucidate how the BX-C boundaries combine these two contradictory activities, we have used a boundary replacement strategy. We show that a 337 bp fragment spanning the Fab-8 boundary nuclease hypersensitive site and lacking all but 83 bp of the 625 bp Fab-8 PTS (promoter targeting sequence) fully rescues a Fab-7 deletion. It blocks crosstalk between the iab-6 and iab-7 regulatory domains, and has bypass activity that enables the two downstream domains, iab-5 and iab-6, to regulate Abdominal-B (Abd-B) transcription in spite of two intervening boundary elements. Fab-8 has two dCTCF sites and we show that they are necessary both for blocking and bypass activity. However, CTCF sites on their own are not sufficient for bypass. While multimerized dCTCF (or Su(Hw)) sites have blocking activity, they fail to support bypass. Moreover, this bypass defect is not rescued by the full length PTS. Finally, we show that orientation is critical for the proper functioning the Fab-8 replacement. Though the inverted Fab-8 boundary still blocks crosstalk, it disrupts the topology of the Abd-B regulatory domains and does not support bypass. Importantly, altering the orientation of the Fab-8 dCTCF sites is not sufficient to disrupt bypass, indicating that orientation dependence is conferred by other factors.

Boundary elements in the Bithorax complex have two seemingly contradictory activities. They must block crosstalk between neighboring regulatory domains, but at the same time be permissive (insulator bypass) for regulatory interactions between the domains and the BX-C homeotic genes. We have used a replacement strategy to investigate how they carry out these two functions. We show that a 337 bp fragment spanning the Fab-8 boundary nuclease hypersensitive site is sufficient to fully rescue a Fab-7 boundary deletion. It blocks crosstalk and supports bypass. As has been observed in transgene assays, blocking activity requires the Fab-8 dCTCF sites, while full bypass activity requires the dCTCF sites plus a small part of PTS. In transgene assays, bypass activity typically depends on the orientation of the two insulators relative to each other. A similar orientation dependence is observed for the Fab-8 replacement in BX-C. When the orientation of the Fab-8 boundary is reversed, bypass activity is lost, while blocking is unaffected. Interestingly, unlike what has been observed in mammals, reversing the orientation of only the Fab-8 dCTCF sites does not affect boundary function. This finding indicates that other Fab-8 factors must play a critical role in determining orientation. Taken together, our findings argue that carrying out the paradoxical functions of the BX-C boundaries does not require any unusual or special properties; rather BX-C boundaries utilize generic blocking and insulator bypass activities that are appropriately adapted to their regulatory context. Thus making them a good model for studying the functional properties of boundaries/insulators in their native setting.

Decoding transcriptional enhancers: Evolving from annotation to functional interpretation

Deciphering the intricate molecular processes that orchestrate the spatial and temporal regulation of genes has become an increasingly major focus of biological research. The differential expression of genes by diverse cell types with a common genome is a hallmark of complex cellular functions, as well as the basis for multicellular life. Importantly, a more coherent understanding of gene regulation is critical for defining developmental processes, evolutionary principles and disease etiologies. Here we present our current understanding of gene regulation by focusing on the role of enhancer elements in these complex processes. Although functional genomic methods have provided considerable advances to our understanding of gene regulation, these assays, which are usually performed on a genome-wide scale, typically provide correlative observations that lack functional interpretation. Recent innovations in genome editing technologies have placed gene regulatory studies at an exciting crossroads, as systematic, functional evaluation of enhancers and other transcriptional regulatory elements can now be performed in a coordinated, high-throughput manner across the entire genome. This review provides insights on transcriptional enhancer function, their role in development and disease, and catalogues experimental tools commonly used to study these elements. Additionally, we discuss the crucial role of novel techniques in deciphering the complex gene regulatory landscape and how these studies will shape future research.

Effect of Interaction between Chromatin Loops on Cell-to-Cell Variability in Gene Expression

According to recent experimental evidence, the interaction between chromatin loops, which can be characterized by three factors-connection pattern, distance between regulatory elements, and communication form, play an important role in determining the level of cell-to-cell variability in gene expression. These quantitative experiments call for a corresponding modeling effect that addresses the question of how changes in these factors affect variability at the expression level in a systematic rather than case-by-case fashion. Here we make such an effort, based on a mechanic model that maps three fundamental patterns for two interacting DNA loops into a 4-state model of stochastic transcription. We first show that in contrast to side-by-side loops, nested loops enhance mRNA expression and reduce expression noise whereas alternating loops have just opposite effects. Then, we compare effects of facilitated tracking and direct looping on gene expression. We find that the former performs better than the latter in controlling mean expression and in tuning expression noise, but this control or tuning is distance-dependent, remarkable for moderate loop lengths, and there is a limit loop length such that the difference in effect between two communication forms almost disappears. Our analysis and results justify the facilitated chromatin-looping hypothesis.

Establishment of a Developmental Compartment Requires Interactions between Three Synergistic Cis-regulatory Modules

The subdivision of cell populations in compartments is a key event during animal development. In Drosophila, the gene apterous (ap) divides the wing imaginal disc in dorsal vs ventral cell lineages and is required for wing formation. ap function as a dorsal selector gene has been extensively studied. However, the regulation of its expression during wing development is poorly understood. In this study, we analyzed ap transcriptional regulation at the endogenous locus and identified three cis-regulatory modules (CRMs) essential for wing development. Only when the three CRMs are combined, robust ap expression is obtained. In addition, we genetically and molecularly analyzed the trans-factors that regulate these CRMs. Our results propose a three-step mechanism for the cell lineage compartment expression of ap that includes initial activation, positive autoregulation and Trithorax-mediated maintenance through separable CRMs.

The separation of cell populations into distinct functional units is essential for both vertebrate and invertebrate animal development. A classical paradigm for this phenomenon is the establishment of developmental compartments during Drosophila wing development. These compartments depend on the restricted expression of two selector genes, engrailed in the posterior compartment and apterous (ap) in the dorsal compartment. Yet, despite the central role these genes and their restricted expression patterns play in Drosophila development, we still do not understand how these patterns are established or maintained. Here, by dissecting the regulatory sequences required for ap expression, we solve this problem for this critical selector gene. We used a combination of experimental approaches to identify and functionally characterize the cis-regulatory modules (CRMs) that regulate ap expression during Drosophila wing development. For these analyses we implement a novel technique allowing us to study the function of these CRMs in vivo, at the native ap locus. We found three ap CRMs crucial for wing development: the Early (apE) and the D/V (apDV) enhancers and the ap PRE (apP). Only when all three regulatory elements are combined is a uniform and complete ap expression domain generated. In summary, our results indicate that ap is regulated in time and space by a three-step mechanism that generates a lineage compartment by integrating input from separate CRMs for the initiation, refinement and maintenance of its expression.

Identifying transcriptional cis-regulatory modules in animal genomes

Gene expression is regulated through the activity of transcription factors (TFs) and chromatin-modifying proteins acting on specific DNA sequences, referred to as cis-regulatory elements. These include promoters, located at the transcription initiation sites of genes, and a variety of distal cis-regulatory modules (CRMs), the most common of which are transcriptional enhancers. Because regulated gene expression is fundamental to cell differentiation and acquisition of new cell fates, identifying, characterizing, and understanding the mechanisms of action of CRMs is critical for understanding development. CRM discovery has historically been challenging, as CRMs can be located far from the genes they regulate, have few readily identifiable sequence characteristics, and for many years were not amenable to high-throughput discovery methods. However, the recent availability of complete genome sequences and the development of next-generation sequencing methods have led to an explosion of both computational and empirical methods for CRM discovery in model and nonmodel organisms alike. Experimentally, CRMs can be identified through chromatin immunoprecipitation directed against TFs or histone post-translational modifications, identification of nucleosome-depleted 'open' chromatin regions, or sequencing-based high-throughput functional screening. Computational methods include comparative genomics, clustering of known or predicted TF-binding sites, and supervised machine-learning approaches trained on known CRMs. All of these methods have proven effective for CRM discovery, but each has its own considerations and limitations, and each is subject to a greater or lesser number of false-positive identifications. Experimental confirmation of predictions is essential, although shortcomings in current methods suggest that additional means of validation need to be developed. For further resources related to this article, please visit the WIREs website.

CONFLICT OF INTEREST

The authors have declared no conflicts of interest for this article.

Quantitation of interactions between two DNA loops demonstrates loop domain insulation in E. coli cells

Eukaryotic gene regulation involves complex patterns of long-range DNA-looping interactions between enhancers and promoters, but how these specific interactions are achieved is poorly understood. Models that posit other DNA loops--that aid or inhibit enhancer-promoter contact--are difficult to test or quantitate rigorously in eukaryotic cells. Here, we use the well-characterized DNA-looping proteins Lac repressor and phage λ CI to measure interactions between pairs of long DNA loops in E. coli cells in the three possible topological arrangements. We find that side-by-side loops do not affect each other. Nested loops assist each other's formation consistent with their distance-shortening effect. In contrast, alternating loops, where one looping element is placed within the other DNA loop, inhibit each other's formation, thus providing clear support for the loop domain model for insulation. Modeling shows that combining loop assistance and loop interference can provide strong specificity in long-range interactions.

Bipartite recognition of DNA by TCF/Pangolin is remarkably flexible and contributes to transcriptional responsiveness and tissue specificity of wingless signaling

The T-cell factor (TCF) family of transcription factors are major mediators of Wnt/β-catenin signaling in metazoans. All TCFs contain a High Mobility Group (HMG) domain that possesses specific DNA binding activity. In addition, many TCFs contain a second DNA binding domain, the C-clamp, which binds to DNA motifs referred to as Helper sites. While HMG and Helper sites are both important for the activation of several Wnt dependent cis-regulatory modules (W-CRMs), the rules of what constitutes a functional HMG-Helper site pair are unknown. In this report, we employed a combination of in vitro binding, reporter gene analysis and bioinformatics to address this question, using the Drosophila family member TCF/Pangolin (TCF/Pan) as a model. We found that while there were constraints for the orientation and spacing of HMG-Helper pairs, the presence of a Helper site near a HMG site in any orientation increased binding and transcriptional response, with some orientations displaying tissue-specific patterns. We found that altering an HMG-Helper site pair from a sub-optimal to optimal orientation/spacing dramatically increased the responsiveness of a W-CRM in several fly tissues. In addition, we used the knowledge gained to bioinformatically identify two novel W-CRMs, one that was activated by Wnt/β-catenin signaling in the prothoracic gland, a tissue not previously connected to this pathway. In sum, this work extends the importance of Helper sites in fly W-CRMs and suggests that the type of HMG-Helper pair is a major factor in setting the threshold for Wnt activation and tissue-responsiveness.

Regulation of gene expression is controlled in large part by proteins known as transcription factors, which bind to specific DNA sequences in the genome. The DNA binding domains of transcription factors recognize short stretches (5–11 base pairs) of DNA with considerable sequence degeneracy. This means that a single DNA binding domain, on its own, cannot find its targets in the vast excess of genomic sequence. We are studying this question using TCF/Pangolin, a Drosophila transcription factor that mediates Wnt/β-catenin signaling, an important developmental cell-cell communication pathway. TCF/Pangolin contains two DNA binding domains that bind to a pair of DNA motifs known as HMG and Helper sites. We used a combination of biochemistry, genetics and bioinformatics to elucidate the spacing and orientation constraints of HMG-Helper site pairs. We found that HMG-Helper site spacing/orientation influenced the sensitivity of a target to Wnt signaling, as well as its tissue-responsiveness. We used this information to improve our ability to search the Drosophila genome for Wnt targets, one of which was activated by the pathway in the fly ring gland, the major endocrine organ in insects. Our work is relevant to related mammalian TCF family members, which are implicated in development, stem cell biology and the progression of cancer.

Co-regulation of invected and engrailed by a complex array of regulatory sequences in Drosophila

invected (inv) and engrailed (en) form a gene complex that extends about 115 kb. These two genes encode highly related homeodomain proteins that are co-regulated in a complex manner throughout development. Our dissection of inv/en regulatory DNA shows that most enhancers are spread throughout a 62 kb region. We used two types of constructs to analyze the function of this DNA: P-element based reporter constructs with small pieces of DNA fused to the en promoter driving lacZ expression and large constructs with HA-tagged en and inv inserted in the genome with the phiC31 system. In addition, we generated deletions of inv and en DNA in situ and assayed their effects on inv/en expression. Our results support and extend our knowledge of inv/en regulation. First, inv and en share regulatory DNA, most of which is flanking the en transcription unit. In support of this, a 79-kb HA-en transgene can rescue inv en double mutants to viable, fertile adults. In contrast, an 84-kb HA-inv transgene lacks most of the enhancers for inv/en expression. Second, there are multiple enhancers for inv/en stripes in embryos; some of these may be redundant but others play discrete roles at different stages of embryonic development. Finally, no small reporter construct gave expression in the posterior compartment of imaginal discs, a hallmark of inv/en expression. Robust expression of HA-en in the posterior compartment of imaginal discs is evident from the 79-kb HA-en transgene, while a 45-kb HA-en transgene gives weaker, variable imaginal disc expression. We suggest that the activity of the imaginal disc enhancer(s) is dependent on the chromatin structure of the inv/en domain.

In search of the determinants of enhancer-promoter interaction specificity

Although it was originally believed that enhancers activate only the nearest promoter, recent global analyses enabled by high-throughput technology suggest that the network of enhancer-promoter interactions is far more complex. The mechanisms that determine the specificity of enhancer-promoter interactions are still poorly understood, but they are thought to include biochemical compatibility, constraints imposed by the three-dimensional architecture of chromosomes, insulator elements, and possibly the effects of local chromatin composition. In this review, we assess the current insights into these determinants, and highlight the functional genomic approaches that will lead the way towards better mechanistic understanding.

In search of the determinants of enhancer-promoter interaction specificity

Although it was originally believed that enhancers activate only the nearest promoter, recent global analyses enabled by high-throughput technology suggest that the network of enhancer-promoter interactions is far more complex. The mechanisms that determine the specificity of enhancer-promoter interactions are still poorly understood, but they are thought to include biochemical compatibility, constraints imposed by the three-dimensional architecture of chromosomes, insulator elements, and possibly the effects of local chromatin composition. In this review, we assess the current insights into these determinants, and highlight the functional genomic approaches that will lead the way towards better mechanistic understanding.

Regulation of gene expression in the genomic context

Metazoan life is dependent on the proper temporal and spatial control of gene expression within the many cells-essentially all with the identical genome-that make up the organism. While much is understood about how individual gene regulatory elements function, many questions remain about how they interact to maintain correct regulation globally throughout the genome. In this review we summarize the basic features and functions of the crucial regulatory elements promoters, enhancers, and insulators and discuss some of the ways in which proper interactions between these elements is realized. We focus in particular on the role of core promoter sequences and propose explanations for some of the contradictory results seen in experiments aimed at understanding insulator function. We suggest that gene regulation depends on local genomic context and argue that more holistic in vivo investigations that take into account multiple local features will be necessary to understand how genome-wide gene regulation is maintained.

Onset of atonal expression in Drosophila retinal progenitors involves redundant and synergistic contributions of Ey/Pax6 and So binding sites within two distant enhancers

Proneural transcription factors drive the generation of specialized neurons during nervous system development, and their dynamic expression pattern is critical to their function. The activation of the proneural gene atonal (ato) in the Drosophila eye disc epithelium represents a critical step in the transition from retinal progenitor cell to developing photoreceptor neuron. We show here that the onset of ato transcription depends on two distant enhancers that function differently in subsets of retinal progenitor cells. A detailed analysis of the crosstalk between these enhancers identifies a critical role for three binding sites for the Retinal Determination factors Eyeless (Ey) and Sine oculis (So). We show how these sites interact to induce ato expression in distinct regions of the eye field and confirm them to be occupied by endogenous Ey and So proteins in vivo. Our study suggests that Ey and So operate differently through the same 3' cis-regulatory sites in distinct populations of retinal progenitors.

Architectural and functional diversity of polycomb group response elements in Drosophila

Polycomb group response elements (PREs) play an essential role in gene regulation by the Polycomb group (PcG) repressor proteins in Drosophila. PREs are required for the recruitment and maintenance of repression by the PcG proteins. PREs are made up of binding sites for multiple DNA-binding proteins, but it is still unclear what combination(s) of binding sites is required for PRE activity. Here we compare the binding sites and activities of two closely linked yet separable PREs of the Drosophila engrailed (en) gene, PRE1 and PRE2. Both PRE1 and PRE2 contain binding sites for multiple PRE-DNA-binding proteins, but the number, arrangement, and spacing of the sites differs between the two PREs. These differences have functional consequences. Both PRE1 and PRE2 mediate pairing-sensitive silencing of mini-white, a functional assay for PcG repression; however, PRE1 requires two binding sites for Pleiohomeotic (Pho), whereas PRE2 requires only one Pho-binding site for this activity. Furthermore, for full pairing-sensitive silencing activity, PRE1 requires an AT-rich region not found in PRE2. These two PREs behave differently in a PRE embryonic and larval reporter construct inserted at an identical location in the genome. Our data illustrate the diversity of architecture and function of PREs.

Cohesin and polycomb proteins functionally interact to control transcription at silenced and active genes

Cohesin is crucial for proper chromosome segregation but also regulates gene transcription and organism development by poorly understood mechanisms. Using genome-wide assays in Drosophila developing wings and cultured cells, we find that cohesin functionally interacts with Polycomb group (PcG) silencing proteins at both silenced and active genes. Cohesin unexpectedly facilitates binding of Polycomb Repressive Complex 1 (PRC1) to many active genes, but their binding is mutually antagonistic at silenced genes. PRC1 depletion decreases phosphorylated RNA polymerase II and mRNA at many active genes but increases them at silenced genes. Depletion of cohesin reduces long-range interactions between Polycomb Response Elements in the invected-engrailed gene complex where it represses transcription. These studies reveal a previously unrecognized role for PRC1 in facilitating productive gene transcription and provide new insights into how cohesin and PRC1 control development.

Cis-regulatory complexity within a large non-coding region in the Drosophila genome

Analysis of cis-regulatory enhancers has revealed that they consist of clustered blocks of highly conserved sequences. Although most characterized enhancers reside near their target genes, a growing number of studies have shown that enhancers located over 50 kb from their minimal promoter(s) are required for appropriate gene expression and many of these ‘long-range’ enhancers are found in genomic regions that are devoid of identified exons. To gain insight into the complexity of Drosophila cis-regulatory sequences within exon-poor regions, we have undertaken an evolutionary analysis of 39 of these regions located throughout the genome. This survey revealed that within these genomic expanses, clusters of conserved sequence blocks (CSBs) are positioned once every 1.1 kb, on average, and that a typical cluster contains multiple (5 to 30 or more) CSBs that have been maintained for at least 190 My of evolutionary divergence. As an initial step toward assessing the cis-regulatory activity of conserved clusters within gene-free genomic expanses, we have tested the in-vivo enhancer activity of 19 consecutive CSB clusters located in the middle of a 115 kb gene-poor region on the 3^rd chromosome. Our studies revealed that each cluster functions independently as a specific spatial/temporal enhancer. In total, the enhancers possess a diversity of regulatory functions, including dynamically activating expression in defined patterns within subsets of cells in discrete regions of the embryo, larvae and/or adult. We also observed that many of the enhancers are multifunctional–that is, they activate expression during multiple developmental stages. By extending these results to the rest of the Drosophila genome, which contains over 70,000 non-coding CSB clusters, we suggest that most function as enhancers.

Enhancer trap infidelity in Drosophila optomotor-blind

Reporter gene activity in enhancer trap lines is often implicitly assumed to mirror quite faithfully the endogenous expression of the "trapped" gene, even though there are numerous examples of enhancer trap infidelity. optomotor-blind (omb) is a 160 kb gene in which 16 independent P-element enhancer trap insertions of three different types have been mapped in a range of more than 60 kb. We have determined the expression pattern of these elements in wing, eye-antennal and leg imaginal discs as well as in the pupal tergites. We noted that one pGawB insertion (omb (P4) ) selectively failed to report parts of the omb pattern even though the missing pattern elements were apparent in all other 15 lines. We ruled out that omb (P4) was defective in the Gal4 promoter region or had inactivated genomic enhancers in the integration process. We propose that the Gal4 reporter gene in pGawB may be sensitive to orientation or promoter proximity effects.

Polycomb group response elements in Drosophila and vertebrates

Polycomb group genes (PcG) encode a group of about 16 proteins that were first identified in Drosophila as repressors of homeotic genes. PcG proteins are present in all metazoans and are best characterized as transcriptional repressors. In Drosophila, these proteins are known as epigenetic regulators because they remember, but do not establish, the patterned expression state of homeotic genes throughout development. PcG proteins, in general, are not DNA binding proteins, but act in protein complexes to repress transcription at specific target genes. How are PcG proteins recruited to the DNA? In Drosophila, there are specific regulatory DNA elements called Polycomb group response elements (PREs) that bring PcG protein complexes to the DNA. Drosophila PREs are made up of binding sites for a complex array of DNA binding proteins. Functional PRE assays in transgenes have shown that PREs act in the context of other regulatory DNA and PRE activity is highly dependent on genomic context. Drosophila PREs tend to regulate genes with a complex array of regulatory DNA in a cell or tissue-specific fashion and it is the interplay between regulatory DNA that dictates PRE function. In mammals, PcG proteins are more diverse and there are multiple ways to recruit PcG complexes, including RNA-mediated recruitment. In this review, we discuss evidence for PREs in vertebrates and explore similarities and differences between Drosophila and vertebrate PREs.

Exploring the effects of polymorphisms on cis-regulatory signal transduction response

cis-Regulatory sequences (CRSs) direct cell-specific and inducible gene expression in response to signal transduction networks, and it is becoming apparent that many cases of disease susceptibility and drug response stratification are due to polymorphisms that alter CRS responses in a context-dependent manner. In the current review, we describe successful methods for identifying CRSs and analyzing the effects of allelic variation on their responses to signal transduction. The technologies described build on the successes of ENCODE (ENCyclopedia Of DNA Elements) by exploring the effects of polymorphisms on CRS context dependency. This understanding is essential to uncover the genomic basis of disease susceptibility and will play a major role in delivering on the promise of personalized medicine.

Predicting spatial and temporal gene expression using an integrative model of transcription factor occupancy and chromatin state

Precise patterns of spatial and temporal gene expression are central to metazoan complexity and act as a driving force for embryonic development. While there has been substantial progress in dissecting and predicting cis-regulatory activity, our understanding of how information from multiple enhancer elements converge to regulate a gene's expression remains elusive. This is in large part due to the number of different biological processes involved in mediating regulation as well as limited availability of experimental measurements for many of them. Here, we used a Bayesian approach to model diverse experimental regulatory data, leading to accurate predictions of both spatial and temporal aspects of gene expression. We integrated whole-embryo information on transcription factor recruitment to multiple cis-regulatory modules, insulator binding and histone modification status in the vicinity of individual gene loci, at a genome-wide scale during Drosophila development. The model uses Bayesian networks to represent the relation between transcription factor occupancy and enhancer activity in specific tissues and stages. All parameters are optimized in an Expectation Maximization procedure providing a model capable of predicting tissue- and stage-specific activity of new, previously unassayed genes. Performing the optimization with subsets of input data demonstrated that neither enhancer occupancy nor chromatin state alone can explain all gene expression patterns, but taken together allow for accurate predictions of spatio-temporal activity. Model predictions were validated using the expression patterns of more than 600 genes recently made available by the BDGP consortium, demonstrating an average 15-fold enrichment of genes expressed in the predicted tissue over a naïve model. We further validated the model by experimentally testing the expression of 20 predicted target genes of unknown expression, resulting in an accuracy of 95% for temporal predictions and 50% for spatial. While this is, to our knowledge, the first genome-wide approach to predict tissue-specific gene expression in metazoan development, our results suggest that integrative models of this type will become more prevalent in the future.

Development is a complex process in which a single cell gives rise to a multi-cellular organism comprised of diverse cell types and well-organized tissues. This transformation requires tightly coordinated expression, both spatially and temporally, of hundreds to thousands of genes specific to any given tissue. To orchestrate these patterns, gene expression is regulated at multiple steps, from TF binding to cis-regulatory modules, general transcription factor and RNA polymerase II recruitment to promoters, chromatin remodeling, and three-dimensional looping interactions. Despite this level of complexity, the regulation of gene expression is typically modeled in the context of transcription factor binding and a single enhancer's activity as this is where the majority of experimental data is available. Recent advances in the measurement of chromatin modifications and insulator binding during embryogenesis provide new datasets that can be used for modeling gene expression. Here we use a Bayesian approach to integrate all three levels of information to combine the activity of multiple regulatory elements into a single model of a gene's expression, implementing an expectation maximization strategy to overcome the problem of missing data. Importantly, while the data for histone modifications and insulator binding represents merged signals from all cells in the embryo, the model can extract cell type specific and stage-specific predictions on gene expression for hundreds of genes of unknown expression.

A Hh-driven gene network controls specification, pattern and size of the Drosophila simple eyes

During development, extracellular signaling molecules interact with intracellular gene networks to control the specification, pattern and size of organs. One such signaling molecule is Hedgehog (Hh). Hh is known to act as a morphogen, instructing different fates depending on the distance to its source. However, how Hh, when signaling across a cell field, impacts organ-specific transcriptional networks is still poorly understood. Here, we investigate this issue during the development of the Drosophila ocellar complex. The development of this sensory structure, which is composed of three simple eyes (or ocelli) located at the vertices of a triangular patch of cuticle on the dorsal head, depends on Hh signaling and on the definition of three domains: two areas of eya and so expression--the prospective anterior and posterior ocelli--and the intervening interocellar domain. Our results highlight the role of the homeodomain transcription factor engrailed (en) both as a target and as a transcriptional repressor of hh signaling in the prospective interocellar region. Furthermore, we identify a requirement for the Notch pathway in the establishment of en maintenance in a Hh-independent manner. Therefore, hh signals transiently during the specification of the interocellar domain, with en being required here for hh signaling attenuation. Computational analysis further suggests that this network design confers robustness to signaling noise and constrains phenotypic variation. In summary, using genetics and modeling we have expanded the ocellar gene network to explain how the interaction between the Hh gradient and this gene network results in the generation of stable mutually exclusive gene expression domains. In addition, we discuss some general implications our model may have in some Hh-driven gene networks.

Polycomb group proteins bind an engrailed PRE in both the "ON" and "OFF" transcriptional states of engrailed

Polycomb group (PcG) and trithorax Group (trxG) proteins maintain the “OFF” and “ON” transcriptional states of HOX genes and other targets by modulation of chromatin structure. In Drosophila, PcG proteins are bound to DNA fragments called Polycomb group response elements (PREs). The prevalent model holds that PcG proteins bind PREs only in cells where the target gene is “OFF”. Another model posits that transcription through PREs disrupts associated PcG complexes, contributing to the establishment of the “ON” transcriptional state. We tested these two models at the PcG target gene engrailed. engrailed exists in a gene complex with invected, which together have 4 well-characterized PREs. Our data show that these PREs are not transcribed in embryos or larvae. We also examined whether PcG proteins are bound to an engrailed PRE in cells where engrailed is transcribed. By FLAG-tagging PcG proteins and expressing them specifically where engrailed is “ON” or “OFF”, we determined that components of three major PcG protein complexes are present at an engrailed PRE in both the “ON” and “OFF” transcriptional states in larval tissues. These results show that PcG binding per se does not determine the transcriptional state of engrailed.

Wapl antagonizes cohesin binding and promotes Polycomb-group silencing in Drosophila

Wapl protein regulates binding of the cohesin complex to chromosomes during interphase and helps remove cohesin from chromosomes at mitosis. We isolated a dominant mutation in wapl (wapl(AG)) in a screen for mutations that counteract silencing mediated by an engrailed Polycomb-group response element. wapl(AG) hemizygotes die as pharate adults and have an extra sex combs phenotype characteristic of males with mutations in Polycomb-group (PcG) genes. The wapl gene encodes two proteins, a long form and a short form. wapl(AG) introduces a stop codon at amino acid 271 of the long form and produces a truncated protein. The expression of a transgene encoding the truncated Wapl-AG protein causes an extra-sex-comb phenotype similar to that seen in the wapl(AG) mutant. Mutations in the cohesin-associated genes Nipped-B and pds5 suppress and enhance wapl(AG) phenotypes, respectively. A Pds5-Wapl complex (releasin) removes cohesin from DNA, while Nipped-B loads cohesin. This suggests that Wapl-AG might exert its effects through changes in cohesin binding. Consistent with this model, Wapl-AG was found to increase the stability of cohesin binding to polytene chromosomes. Our data suggest that increasing cohesin stability interferes with PcG silencing at genes that are co-regulated by cohesin and PcG proteins.

Comparing enhancer action in cis and in trans

Studies from diverse systems have shown that distinct interchromosomal interactions are a central component of nuclear organization. In some cases, these interactions allow an enhancer to act in trans, modulating the expression of a gene encoded on a separate chromosome held in close proximity. Despite recent advances in uncovering such phenomena, our understanding of how a regulatory element acts on another chromosome remains incomplete. Here, we describe a transgenic approach to better understand enhancer action in trans in Drosophila melanogaster. Using phiC31-based recombinase-mediated cassette exchange (RMCE), we placed transgenes carrying combinations of the simple enhancer GMR, a minimal promoter, and different fluorescent reporters at equivalent positions on homologous chromosomes so that they would pair via the endogenous somatic pairing machinery of Drosophila. Our data demonstrate that the enhancer GMR is capable of activating a promoter in trans and does so in a variegated pattern, suggesting stochastic interactions between the enhancer and the promoter when they are carried on separate chromosomes. Furthermore, we quantitatively assessed the impact of two concurrent promoter targets in cis and in trans to GMR, demonstrating that each promoter is capable of competing for the enhancer's activity, with the presence of one negatively affecting expression from the other. Finally, the single-cell resolution afforded by our approach allowed us to show that promoters in cis and in trans to GMR can both be activated in the same nucleus, implying that a single enhancer can share its activity between multiple promoter targets carried on separate chromosomes.

Shadow enhancers: frequently asked questions about distributed cis-regulatory information and enhancer redundancy

This paper, in the form of a frequently asked questions page (FAQ), addresses outstanding questions about "shadow enhancers", quasi-redundant cis-regulatory elements, and their proposed roles in transcriptional control. Questions include: What exactly are shadow enhancers? How many genes have shadow/redundant/distributed enhancers? How redundant are these elements? What is the function of distributed enhancers? How modular are enhancers? Is it useful to study a single enhancer in isolation? In addition, a revised definition of "shadow enhancers" is proposed, and possible mechanisms of shadow enhancer function and evolution are discussed.

P-element homing is facilitated by engrailed polycomb-group response elements in Drosophila melanogaster

P-element vectors are commonly used to make transgenic Drosophila and generally insert in the genome in a nonselective manner. However, when specific fragments of regulatory DNA from a few Drosophila genes are incorporated into P-transposons, they cause the vectors to be inserted near the gene from which the DNA fragment was derived. This is called P-element homing. We mapped the minimal DNA fragment that could mediate homing to the engrailed/invected region of the genome. A 1.6 kb fragment of engrailed regulatory DNA that contains two Polycomb-group response elements (PREs) was sufficient for homing. We made flies that contain a 1.5kb deletion of engrailed DNA (en^Δ1.5) in situ, including the PREs and the majority of the fragment that mediates homing. Remarkably, homing still occurs onto the en^{Δ1. 5} chromosome. In addition to homing to en, P[en] inserts near Polycomb group target genes at an increased frequency compared to P[EPgy2], a vector used to generate 18,214 insertions for the Drosophila gene disruption project. We suggest that homing is mediated by interactions between multiple proteins bound to the homing fragment and proteins bound to multiple areas of the engrailed/invected chromatin domain. Chromatin structure may also play a role in homing.

Disruption of the abdominal-B promoter tethering element results in a loss of long-range enhancer-directed Hox gene expression in Drosophila

There are many examples within gene complexes of transcriptional enhancers interacting with only a subset of target promoters. A number of molecular mechanisms including promoter competition, insulators and chromatin looping are thought to play a role in regulating these interactions. At the Drosophila bithorax complex (BX-C), the IAB5 enhancer specifically drives gene expression only from the Abdominal-B (Abd-B) promoter, even though the enhancer and promoter are 55 kb apart and are separated by at least three insulators. In previous studies, we discovered that a 255 bp cis-regulatory module, the promoter tethering element (PTE), located 5′ of the Abd-B transcriptional start site is able to tether IAB5 to the Abd-B promoter in transgenic embryo assays. In this study we examine the functional role of the PTE at the endogenous BX-C using transposon-mediated mutagenesis. Disruption of the PTE by P element insertion results in a loss of enhancer-directed Abd-B expression during embryonic development and a homeotic transformation of abdominal segments. A partial deletion of the PTE and neighboring upstream genomic sequences by imprecise excision of the P element also results in a similar loss of Abd-B expression in embryos. These results demonstrate that the PTE is an essential component of the regulatory network at the BX-C and is required in vivo to mediate specific long-range enhancer-promoter interactions.

Chromosomal organization at the level of gene complexes

Metazoan genomes primarily consist of non-coding DNA in comparison to coding regions. Non-coding fraction of the genome contains cis-regulatory elements, which ensure that the genetic code is read properly at the right time and space during development. Regulatory elements and their target genes define functional landscapes within the genome, and some developmentally important genes evolve by keeping the genes involved in specification of common organs/tissues in clusters and are termed gene complex. The clustering of genes involved in a common function may help in robust spatio-temporal gene expression. Gene complexes are often found to be evolutionarily conserved, and the classic example is the hox complex. The evolutionary constraints seen among gene complexes provide an ideal model system to understand cis and trans-regulation of gene function. This review will discuss the various characteristics of gene regulatory modules found within gene complexes and how they can be characterized.