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

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.

Mask exhibits trxG-like behavior and associates with H3K27ac marked chromatin

The Trithorax group (trxG) proteins counteract the repressive effect of Polycomb group (PcG) complexes and maintain transcriptional memory of active states of key developmental genes. Although chromatin structure and modifications appear to play a fundamental role in this process, it is not clear how trxG prevents PcG-silencing and heritably maintains an active gene expression state. Here, we report a hitherto unknown role of Drosophila Multiple ankyrin repeats single KH domain (Mask), which emerged as one of the candidate trxG genes in our reverse genetic screen. The genome-wide binding profile of Mask correlates with known trxG binding sites across the Drosophila genome. In particular, the association of Mask at chromatin overlaps with CBP and H3K27ac, which are known hallmarks of actively transcribed genes by trxG. Importantly, Mask predominantly associates with actively transcribed genes in Drosophila. Depletion of Mask not only results in the downregulation of trxG targets but also correlates with diminished levels of H3K27ac. The fact that Mask positively regulates H3K27ac levels in flies was also found to be conserved in human cells. Strong suppression of Pc mutant phenotype by mutation in mask provides physiological relevance that Mask contributes to the anti-silencing effect of trxG, maintaining expression of key developmental genes. Since Mask is a downstream effector of multiple cell signaling pathways, we propose that Mask may connect cell signaling with chromatin mediated epigenetic cell memory governed by trxG.

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.

Interlocking of co-opted developmental gene networks in Drosophila and the evolution of pre-adaptive novelty

The re-use of genes in new organs forms the base of many evolutionary novelties. A well-characterised case is the recruitment of the posterior spiracle gene network to the Drosophila male genitalia. Here we find that this network has also been co-opted to the testis mesoderm where is required for sperm liberation, providing an example of sequentially repeated developmental co-options. Associated to this co-option event, an evolutionary expression novelty appeared, the activation of the posterior segment determinant Engrailed to the anterior A8 segment controlled by common testis and spiracle regulatory elements. Enhancer deletion shows that A8 anterior Engrailed activation is not required for spiracle development but only necessary in the testis. Our study presents an example of pre-adaptive developmental novelty: the activation of the Engrailed transcription factor in the anterior compartment of the A8 segment where, despite having no specific function, opens the possibility of this developmental factor acquiring one. We propose that recently co-opted networks become interlocked, so that any change to the network because of its function in one organ, will be mirrored by other organs even if it provides no selective advantage to them.

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.

Functional Divergence of the Tribolium castaneum engrailed and invected Paralogs

Engrailed (en) and invected (inv) encode paralogous transcription factors found as a closely linked tandem duplication within holometabolous insects. Drosophila en mutants segment normally, then fail to maintain their segments. Loss of Drosophila inv is viable, while loss of both genes results in asegmental larvae. Surprisingly, the knockdown of Oncopeltus inv can result in the loss or fusion of the entire abdomen and en knockdowns in Tribolium show variable degrees of segmental loss. The consequence of losing or knocking down both paralogs on embryogenesis has not been studied beyond Drosophila. To further investigate the relative functions of each paralog and the mechanism behind the segmental loss, Tribolium double and single knockdowns of en and inv were analyzed. The most common cuticular phenotype of the double knockdowns was small, limbless, and open dorsally, with all but a single, segmentally iterated row of bristles. Less severe knockdowns had fused segments and reduced appendages. The Tribolium paralogs appear to act synergistically: the knockdown of either Tribolium gene alone was typically less severe, with all limbs present, whereas the most extreme single knockdowns mimic the most severe double knockdown phenotype. Morphological abnormalities unique to either single gene knockdown were not found. inv expression was not affected in the Tribolium en knockdowns, but hh expression was unexpectedly increased midway through development. Thus, while the segmental expression of en/inv is broadly conserved within insects, the functions of en and inv are evolving independently in different lineages.

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.

Divergent expression of aristaless1 and aristaless2 during embryonic appendage and pupal wing development in butterflies

BACKGROUND

Gene duplication events are critical for the evolution of new gene functions. Aristaless is a major regulator of distinct developmental processes. It is most known for its role during appendage development across animals. However, more recently other distinct biological functions have been described for this gene and its duplicates. Butterflies and moths have two copies of aristaless, aristaless1 (al1) and aristaless2 (al2), as a result of a gene duplication event. Previous work in Heliconius has shown that both copies appear to have novel functions related to wing color patterning. Here we expand our knowledge of the expression profiles associated with both ancestral and novel functions of Al1 across embryogenesis and wing pigmentation. Furthermore, we characterize Al2 expression, providing a comparative framework between gene copies within the same species, allowing us to understand the origin of new functions following gene duplication.

RESULTS

Our work shows that the expression of both Al1 and Al2 is associated with the ancestral function of sensory appendage (leg, mouth, spines, and eyes) development in embryos. Interestingly, Al1 exhibits higher expression earlier in embryogenesis while the highest levels of Al2 expression are shifted to later stages of embryonic development. Furthermore, Al1 localization appears extranuclear while Al2 co-localizes tightly with nuclei earlier, and then also expands outside the nucleus later in development. Cellular expression of Al1 and Al2 in pupal wings is broadly consistent with patterns observed during embryogenesis. We also describe, for the first time, how Al1 localization appears to correlate with zones of anterior/posterior elongation of the body during embryonic growth, showcasing a possible new function related to Aristaless' previously described role in appendage extension.

CONCLUSIONS

Overall, our data suggest that while both gene copies play a role in embryogenesis and wing pigmentation, the duplicates have diverged temporally and mechanistically across those functions. Our study helps clarify principles behind sub-functionalization and gene expression evolution associated with developmental functions following gene duplication events.

The roles of growth regulation and appendage patterning genes in the morphogenesis of treehopper pronota

Treehoppers of the insect family Membracidae have evolved enlarged and elaborate pronotal structures, which is hypothesized to involve co-opted expression of genes that are shared with the wings. Here, we investigate the similarity between the pronotum and wings in relation to growth. Our study reveals that the ontogenetic allometry of the pronotum is similar to that of wings in Membracidae, but not the outgroup. Using transcriptomics, we identify genes related to translation and protein synthesis, which are mutually upregulated. These genes are implicated in the eIF2, eIF4/p70S6K and mTOR pathways, and have known roles in regulating cell growth and proliferation. We find that species-specific differential growth patterning of the pronotum begins as early as the third instar, which suggests that expression of appendage patterning genes occurs long before the metamorphic molt. We propose that a network related to growth and size determination is the more likely mechanism shared with wings. However, regulators upstream of the shared genes in pronotum and wings need to be elucidated to substantiate whether co-option has occurred. Finally, we believe it will be helpful to distinguish the mechanisms leading to pronotal size from those regulating pronotal shape as we make sense of this spectacular evolutionary innovation.

Transcriptional coupling of distant regulatory genes in living embryos

The prevailing view of metazoan gene regulation is that individual genes are independently regulated by their own dedicated sets of transcriptional enhancers. Past studies have reported long-range gene-gene associations^1-3, but their functional importance in regulating transcription remains unclear. Here we used quantitative single-cell live imaging methods to provide a demonstration of co-dependent transcriptional dynamics of genes separated by large genomic distances in living Drosophila embryos. We find extensive physical and functional associations of distant paralogous genes, including co-regulation by shared enhancers and co-transcriptional initiation over distances of nearly 250 kilobases. Regulatory interconnectivity depends on promoter-proximal tethering elements, and perturbations in these elements uncouple transcription and alter the bursting dynamics of distant genes, suggesting a role of genome topology in the formation and stability of co-transcriptional hubs. Transcriptional coupling is detected throughout the fly genome and encompasses a broad spectrum of conserved developmental processes, suggesting a general strategy for long-range integration of gene activity.

A shared ancient enhancer element differentially regulates the bric-a-brac tandem gene duplicates in the developing Drosophila leg

Gene duplications and transcriptional enhancer emergence/modifications are thought having greatly contributed to phenotypic innovations during animal evolution. Nevertheless, little is known about how enhancers evolve after gene duplication and how regulatory information is rewired between duplicated genes. The Drosophila melanogaster bric-a-brac (bab) complex, comprising the tandem paralogous genes bab1 and bab2, provides a paradigm to address these issues. We previously characterized an intergenic enhancer (named LAE) regulating bab2 expression in the developing legs. We show here that bab2 regulators binding directly the LAE also govern bab1 expression in tarsal cells. LAE excision by CRISPR/Cas9-mediated genome editing reveals that this enhancer appears involved but not strictly required for bab1 and bab2 co-expression in leg tissues. Instead, the LAE enhancer is critical for paralog-specific bab2 expression along the proximo-distal leg axis. Chromatin features and phenotypic rescue experiments indicate that LAE functions partly redundantly with leg-specific regulatory information overlapping the bab1 transcription unit. Phylogenomics analyses indicate that (i) the bab complex originates from duplication of an ancestral singleton gene early on within the Cyclorrhapha dipteran sublineage, and (ii) LAE sequences have been evolutionarily-fixed early on within the Brachycera suborder thus predating the gene duplication event. This work provides new insights on enhancers, particularly about their emergence, maintenance and functional diversification during evolution.

The Hippo pathway coactivator Yorkie can reprogram cell fates and create compartment-boundary-like interactions at clone margins

During development, tissue-specific patterns of gene expression are established by transcription factors and then stably maintained via epigenetic mechanisms. Cancer cells often express genes that are inappropriate for that tissue or developmental stage. Here, we show that high activity levels of Yki, the Hippo pathway coactivator that causes overgrowth in Drosophila imaginal discs, can also disrupt cell fates by altering expression of selector genes like engrailed (en) and Ultrabithorax (Ubx). Posterior clones expressing activated Yki can down-regulate en and express an anterior selector gene, cubitus interruptus (ci). The microRNA bantam and the chromatin regulator Taranis both function downstream of Yki in promoting ci expression. The boundary between Yki-expressing posterior clones and surrounding wild-type cells acquires properties reminiscent of the anteroposterior compartment boundary; Hedgehog signaling pathway activation results in production of Dpp. Thus, at least in principle, heterotypic interactions between Yki-expressing cells and their neighbors could activate boundary-specific signaling mechanisms.

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.

Mimicry diversification in Papilio dardanus via a genomic inversion in the regulatory region of engrailed-invected

Polymorphic Batesian mimics exhibit multiple protective morphs that each mimic a different noxious model. Here, we study the genomic transitions leading to the evolution of different mimetic wing patterns in the polymorphic Mocker Swallowtail Papilio dardanus. We generated a draft genome (231 Mb over 30 chromosomes) and re-sequenced individuals of three morphs. Genome-wide single nucleotide polymorphism (SNP) analysis revealed elevated linkage disequilibrium and divergence between morphs in the regulatory region of engrailed, a developmental gene previously implicated in the mimicry switch. The diverged region exhibits a discrete chromosomal inversion (of 40 kb) relative to the ancestral orientation that is associated with the cenea morph, but not with the bottom-recessive hippocoonides morph or with non-mimetic allopatric populations. The functional role of this inversion in the expression of the novel phenotype is currently unknown, but by preventing recombination, it allows the stable inheritance of divergent alleles enabling geographic spread and local coexistence of multiple adaptive morphs.

Zygotic pioneer factor activity of Odd-paired/Zic is necessary for late function of the Drosophila segmentation network

Because chromatin determines whether information encoded in DNA is accessible to transcription factors, dynamic chromatin states in development may constrain how gene regulatory networks impart embryonic pattern. To determine the interplay between chromatin states and regulatory network function, we performed ATAC-seq on Drosophila embryos during the establishment of the segmentation network, comparing wild-type and mutant embryos in which all graded maternal patterning inputs are eliminated. While during the period between zygotic genome activation and gastrulation many regions maintain stable accessibility, cis-regulatory modules (CRMs) within the network undergo extensive patterning-dependent changes in accessibility. A component of the network, Odd-paired (opa), is necessary for pioneering accessibility of late segmentation network CRMs. opa-driven changes in accessibility are accompanied by equivalent changes in gene expression. Interfering with the timing of opa activity impacts the proper patterning of expression. These results indicate that dynamic systems for chromatin regulation directly impact the reading of embryonic patterning information.

Disparate expression specificities coded by a shared Hox-C enhancer

Can a single regulatory sequence be shared by two genes undergoing functional divergence? Here we describe a single promiscuous enhancer within the Drosophila Antennapedia Complex, EO053, that directs aspects of the expression of two adjacent genes, pb (a Hox2 ortholog) and zen2 (a divergent Hox3 paralog), with disparate spatial and temporal expression patterns. We were unable to separate the pb-like and zen2-like specificities within EO053, and we identify sequences affecting both expression patterns. Importantly, genomic deletion experiments demonstrate that EO053 cooperates with additional pb- and zen2-specific enhancers to regulate the mRNA expression of both genes. We examine sequence conservation of EO053 within the Schizophora, and show that patterns of synteny between the Hox2 and Hox3 orthologs in Arthropods are consistent with a shared regulatory relationship extending prior to the Hox3/zen divergence. Thus, EO053 represents an example of two genes having evolved disparate outputs while utilizing this shared regulatory region.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Molecular mechanisms underlying simplification of venation patterns in holometabolous insects

How mechanisms of pattern formation evolve has remained a central research theme in the field of evolutionary and developmental biology. The mechanism of wing vein differentiation in Drosophila is a classic text-book example of pattern formation using a system of positional-information, yet very little is known about how species with a different number of veins pattern their wings, and how insect venation patterns evolved. Here, we examine the expression pattern of genes previously implicated in vein differentiation in Drosophila in two butterfly species with more complex venation Bicyclus anynana and Pieris canidia. We also test the function of some of these genes in B. anynana. We identify both conserved as well as new domains of decapentaplegic, engrailed, invected, spalt, optix, wingless, armadillo, blistered, and rhomboid gene expression in butterflies, and propose how the simplified venation in Drosophila might have evolved via loss of decapentaplegic, spalt and optix gene expression domains, silencing of vein inducing programs at Spalt-expression boundaries, and changes in gene expression of vein maintenance genes.

Two pair-rule responsive enhancers regulate wingless transcription in the Drosophila blastoderm embryo

BACKGROUND

While many developmentally relevant enhancers act in a modular fashion, there is growing evidence for nonadditive interactions between distinct cis-regulatory enhancers. We investigated if nonautonomous enhancer interactions underlie transcription regulation of the Drosophila segment polarity gene, wingless.

RESULTS

We identified two wg enhancers active at the blastoderm stage: wg 3613u, located from -3.6 to -1.3 kb upstream of the wg transcription start site (TSS) and 3046d, located in intron two of the wg gene, from 3.0 to 4.6 kb downstream of the TSS. Genetic experiments confirm that Even Skipped (Eve), Fushi-tarazu (Ftz), Runt, Odd-paired (Opa), Odd-skipped (Odd), and Paired (Prd) contribute to spatially regulated wg expression. Interestingly, there are enhancer specific differences in response to the gain or loss of function of pair-rule gene activity. Although each element recapitulates aspects of wg expression, a composite reporter containing both enhancers more faithfully recapitulates wg regulation than would be predicted from the sum of their individual responses.

CONCLUSION

These results suggest that the regulation of wg by pair-rule genes involves nonadditive interactions between distinct cis-regulatory enhancers.

Gene activation by dCas9-CBP and the SAM system differ in target preference

Gene overexpression through the targeting of transcription activation domains to regulatory DNA via catalytically defective Cas9 (dCas9) represents a powerful approach to investigate gene function as well as the mechanisms of gene control. To date, the most efficient dCas9-based activator is the Synergistic Activation Mediator (SAM) system whereby transcription activation domains are directly fused to dCas9 as well as tethered through MS2 loops engineered into the gRNA. Here, we show that dCas9 fused to the catalytic domain of the histone acetyltransferase CBP is a more potent activator than the SAM system at some loci, but less efficient at other locations in Drosophila cells. Our results suggest that different rate-limiting steps in the transcription cycle are affected by dCas9-CBP and the SAM system, and that comparing these activators may be useful for mechanistic studies of transcription as well as for increasing the number of hits in genome-wide overexpression screens.

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.

Diverse Cis-Regulatory Mechanisms Contribute to Expression Evolution of Tandem Gene Duplicates

Pairs of duplicated genes generally display a combination of conserved expression patterns inherited from their unduplicated ancestor and newly acquired domains. However, how the cis-regulatory architecture of duplicated loci evolves to produce these expression patterns is poorly understood. We have directly examined the gene-regulatory evolution of two tandem duplicates, the Drosophila Ly6 genes CG9336 and CG9338, which arose at the base of the drosophilids between 40 and 60 Ma. Comparing the expression patterns of the two paralogs in four Drosophila species with that of the unduplicated ortholog in the tephritid Ceratitis capitata, we show that they diverged from each other as well as from the unduplicated ortholog. Moreover, the expression divergence appears to have occurred close to the duplication event and also more recently in a lineage-specific manner. The comparison of the tissue-specific cis-regulatory modules (CRMs) controlling the paralog expression in the four Drosophila species indicates that diverse cis-regulatory mechanisms, including the novel tissue-specific enhancers, differential inactivation, and enhancer sharing, contributed to the expression evolution. Our analysis also reveals a surprisingly variable cis-regulatory architecture, in which the CRMs driving conserved expression domains change in number, location, and specificity. Altogether, this study provides a detailed historical account that uncovers a highly dynamic picture of how the paralog expression patterns and their underlying cis-regulatory landscape evolve. We argue that our findings will encourage studying cis-regulatory evolution at the whole-locus level to understand how interactions between enhancers and other regulatory levels shape the evolution of gene expression.

Chromatin profiling of Drosophila CNS subpopulations identifies active transcriptional enhancers

One of the key issues in studying transcriptional regulation during development is how to employ genome-wide assays that reveals sites of open chromatin and transcription factor binding to efficiently identify biologically relevant genes and enhancers. Analysis of Drosophila CNS midline cell development provides a useful system for studying transcriptional regulation at the genomic level due to a large, well-characterized set of midline-expressed genes and in vivo validated enhancers. In this study, FAIRE-seq on FACS-purified midline cells was performed and the midline FAIRE data were compared with whole-embryo FAIRE data. We find that regions of the genome with a strong midline FAIRE peak and weak whole-embryo FAIRE peak overlap with known midline enhancers and provide a useful predictive tool for enhancer identification. In a complementary analysis, we compared a large dataset of fragments that drive midline expression in vivo with the FAIRE data. Midline enhancer fragments with a midline FAIRE peak tend to be near midline-expressed genes, whereas midline enhancers without a midline FAIRE peak were often distant from midline-expressed genes and unlikely to drive midline transcription in vivo.

Atrophin controls developmental signaling pathways via interactions with Trithorax-like

Mutations in human Atrophin1, a transcriptional corepressor, cause dentatorubral-pallidoluysian atrophy, a neurodegenerative disease. Drosophila Atrophin (Atro) mutants display many phenotypes, including neurodegeneration, segmentation, patterning and planar polarity defects. Despite Atro's critical role in development and disease, relatively little is known about Atro's binding partners and downstream targets. We present the first genomic analysis of Atro using ChIP-seq against endogenous Atro. ChIP-seq identified 1300 potential direct targets of Atro including engrailed, and components of the Dpp and Notch signaling pathways. We show that Atro regulates Dpp and Notch signaling in larval imaginal discs, at least partially via regulation of thickveins and fringe. In addition, bioinformatics analyses, sequential ChIP and coimmunoprecipitation experiments reveal that Atro interacts with the Drosophila GAGA Factor, Trithorax-like (Trl), and they bind to the same loci simultaneously. Phenotypic analyses of Trl and Atro clones suggest that Atro is required to modulate the transcription activation by Trl in larval imaginal discs. Taken together, these data indicate that Atro is a major Trl cofactor that functions to moderate developmental gene transcription.

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.

Co-option of an Ancestral Hox-Regulated Network Underlies a Recently Evolved Morphological Novelty

The evolutionary origins of complex morphological structures such as the vertebrate eye or insect wing remain one of the greatest mysteries of biology. Recent comparative studies of gene expression imply that new structures are not built from scratch, but rather form by co-opting preexisting gene networks. A key prediction of this model is that upstream factors within the network will activate their preexisting targets (i.e., enhancers) to form novel anatomies. Here, we show how a recently derived morphological novelty present in the genitalia of D. melanogaster employs an ancestral Hox-regulated network deployed in the embryo to generate the larval posterior spiracle. We demonstrate how transcriptional enhancers and constituent transcription factor binding sites are used in both ancestral and novel contexts. These results illustrate network co-option at the level of individual connections between regulatory genes and highlight how morphological novelty may originate through the co-option of networks controlling seemingly unrelated structures.