cis-Decoder discovers constellations of conserved DNA sequences shared among tissue-specific enhancers

A Notch-dependent transcriptional mechanism controls expression of temporal patterning factors in Drosophila medulla

Temporal patterning is an important mechanism for generating a great diversity of neuron subtypes from a seemingly homogenous progenitor pool in both vertebrates and invertebrates. Drosophila neuroblasts are temporally patterned by sequentially expressed Temporal Transcription Factors (TTFs). These TTFs are proposed to form a transcriptional cascade based on mutant phenotypes, although direct transcriptional regulation between TTFs has not been verified in most cases. Furthermore, it is not known how the temporal transitions are coupled with the generation of the appropriate number of neurons at each stage. We use neuroblasts of the Drosophila optic lobe medulla to address these questions and show that the expression of TTFs Sloppy-paired 1/2 (Slp1/2) is directly regulated at the transcriptional level by two other TTFs and the cell-cycle dependent Notch signaling through two cis-regulatory elements. We also show that supplying constitutively active Notch can rescue the delayed transition into the Slp stage in cell cycle arrested neuroblasts. Our findings reveal a novel Notch-pathway dependent mechanism through which the cell cycle progression regulates the timing of a temporal transition within a TTF transcriptional cascade.

Annotating the Insect Regulatory Genome

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

Ultraconserved Non-coding DNA Within Diptera and Hymenoptera

This study has taken advantage of the availability of the assembled genomic sequence of flies, mosquitos, ants and bees to explore the presence of ultraconserved sequence elements in these phylogenetic groups. We compared non-coding sequences found within and flanking Drosophila developmental genes to homologous sequences in Ceratitis capitata and Musca domestica. Many of the conserved sequence blocks (CSBs) that constitute Drosophila cis-regulatory DNA, recognized by EvoPrinter alignment protocols, are also conserved in Ceratitis and Musca. Also conserved is the position but not necessarily the orientation of many of these ultraconserved CSBs (uCSBs) with respect to flanking genes. Using the mosquito EvoPrint algorithm, we have also identified uCSBs shared among distantly related mosquito species. Side by side comparison of bee and ant EvoPrints of selected developmental genes identify uCSBs shared between these two Hymenoptera, as well as less conserved CSBs in either one or the other taxon but not in both. Analysis of uCSBs in these dipterans and Hymenoptera will lead to a greater understanding of their evolutionary origin and function of their conserved non-coding sequences and aid in discovery of core elements of enhancers.

This study applies the phylogenetic footprinting program EvoPrinter to detection of ultraconserved non-coding sequence elements in Diptera, including flies and mosquitos, and Hymenoptera, including ants and bees. EvoPrinter outputs an interspecies comparison as a single sequence in terms of the input reference sequence. Ultraconserved sequences flanking known developmental genes were detected in Ceratitis and Musca when compared with Drosophila species, in Aedes and Culex when compared with Anopheles, and between ants and bees. Our methods are useful in detecting and understanding the core evolutionarily hardened sequences required for gene regulation.

REDfly: the transcriptional regulatory element database for Drosophila

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

CORL Expression in the Drosophila Central Nervous System Is Regulated by Stage Specific Interactions of Intertwined Activators and Repressors

CORL proteins (SKOR in mice and Fussel in humans) are a subfamily of central nervous system (CNS) specific proteins related to Sno/Ski oncogenes. Their developmental and homeostatic roles are largely unknown. We previously showed that Drosophila CORL (dCORL; fussel in Flybase) functions between the Activin receptor Baboon and Ecdysone Receptor-B1 (EcR-B1) activation in mushroom body neurons of third instar larval brains. To better understand dCORL regulation and function we generated a series of reporter genes. We examined the embryonic and larval CNS and found that dCORL is regulated by stage specific interactions between intertwined activators and repressors spanning numerous reporters. The reporter AH.lacZ, which contains sequences 7-11kb upstream of dCORL exon1, reflects dCORL brain expression at all stages. Surprisingly, AH.lacZ was not detected in EcR-B1 expressing mushroom body neurons. In larvae AH.lacZ is coexpressed with Elav and the transcription factor Drifter in dILP2 insulin producing cells of the pars intercerebralis. The presence of dCORL in insulin producing cells suggests that dCORL functions non-autonomously in the regulation of EcR-B1 mushroom body activation via the modulation of insulin signaling. Overall, the high level of sequence conservation seen in all CORL/SKOR/Fussel family members and their common CNS specificity suggest that similarly complex regulation and a potential function in insulin signaling are associated with SKOR/Fussel proteins in mammals.

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.

Mechanistic differences in the transcriptional interpretation of local and long-range Shh morphogen signaling

Morphogens orchestrate tissue patterning in a concentration-dependent fashion during vertebrate embryogenesis, yet little is known of how positional information provided by such signals is translated into discrete transcriptional outputs. Here we have identified and characterized cis-regulatory modules (CRMs) of genes operating downstream of graded Shh signaling and bifunctional Gli proteins in neural patterning. Unexpectedly, we find that Gli activators have a noninstructive role in long-range patterning and cooperate with SoxB1 proteins to facilitate a largely concentration-independent mode of gene activation. Instead, the opposing Gli-repressor gradient is interpreted at transcriptional levels, and, together with CRM-specific repressive input of homeodomain proteins, comprises a repressive network that translates graded Shh signaling into regional gene expression patterns. Moreover, local and long-range interpretation of Shh signaling differs with respect to CRM context sensitivity and Gli-activator dependence, and we propose that these differences provide insight into how morphogen function may have mechanistically evolved from an initially binary inductive event.

Identification of cis-Regulatory Elements in the dmyc Gene of Drosophila Melanogaster

Myc is a crucial regulator of growth and proliferation during animal development. Many signals and transcription factors lead to changes in the expression levels of Drosophila myc, yet no clear model exists to explain the complexity of its regulation at the level of transcription. In this study we used Drosophila genetic tools to track the dmyc cis-regulatory elements. Bioinformatics analyses identified conserved sequence blocks in the noncoding regions of the dmyc gene. Investigation of lacZ reporter activity driven by upstream, downstream, and intronic sequences of the dmyc gene in embryonic, larval imaginal discs, larval brain, and adult ovaries, revealed that it is likely to be transcribed from multiple transcription initiation units including a far upstream regulatory region, a TATA box containing proximal complex and a TATA-less downstream promoter element in conjunction with an initiator within the intron 2 region. Our data provide evidence for a modular organization of dmyc regulatory sequences; these modules will most likely be required to generate the tissue-specific patterns of dmyc transcripts. The far upstream region is active in late embryogenesis, while activity of other cis elements is evident during embryogenesis, in specific larval imaginal tissues and during oogenesis. These data provide a framework for further investigation of the transcriptional regulatory mechanisms of dmyc.

Regulation of a duplicated locus: Drosophila sloppy paired is replete with functionally overlapping enhancers

In order to investigate regulation and redundancy within the sloppy paired (slp) locus, we analyzed 30 kilobases of DNA encompassing the tandem, coordinately regulated slp1 and slp2 transcription units. We found a remarkable array of stripe enhancers with overlapping activities surrounding the slp1 transcription unit, and, unexpectedly, glial cell enhancers surrounding slp2. The slp stripe regulatory region generates 7 stripes at blastoderm, and later 14 stripes that persist throughout embryogenesis. Phylogenetic analysis among drosophilids suggests that the multiplicity of stripe enhancers did not evolve through recent duplication. Most of the direct integration among cis-regulatory modules appears to be simply additive, with one notable exception. Despite the apparent redundancy among stripe enhancers, transgenic rescue suggests that most are required for full function, to maintain wingless expression and parasegment boundaries throughout embryogenesis. Transgenic rescue also reveals indirect positive autoregulation by the 7 early stripes, without which alternate stripes within the 14-stripe pattern are lost, leading to embryos with a pair-rule phenotype.

Use of a Drosophila genome-wide conserved sequence database to identify functionally related cis-regulatory enhancers

Background: Phylogenetic footprinting has revealed that cis-regulatory enhancers consist of conserved DNA sequence clusters (CSCs). Currently, there is no systematic approach for enhancer discovery and analysis that takes full-advantage of the sequence information within enhancer CSCs. Results: We have generated a Drosophila genome-wide database of conserved DNA consisting of >100,000 CSCs derived from EvoPrints spanning over 90% of the genome. cis-Decoder database search and alignment algorithms enable the discovery of functionally related enhancers. The program first identifies conserved repeat elements within an input enhancer and then searches the database for CSCs that score highly against the input CSC. Scoring is based on shared repeats as well as uniquely shared matches, and includes measures of the balance of shared elements, a diagnostic that has proven to be useful in predicting cis-regulatory function. To demonstrate the utility of these tools, a temporally-restricted CNS neuroblast enhancer was used to identify other functionally related enhancers and analyze their structural organization. Conclusions:cis-Decoder reveals that co-regulating enhancers consist of combinations of overlapping shared sequence elements, providing insights into the mode of integration of multiple regulating transcription factors. The database and accompanying algorithms should prove useful in the discovery and analysis of enhancers involved in any developmental process. Developmental Dynamics 241:169–189, 2012. © 2011 Wiley Periodicals, Inc.

Functional analysis of conserved sequences within a temporally restricted neural precursor cell enhancer

Many of the key regulators of Drosophila CNS neural identity are expressed in defined temporal orders during neuroblast (NB) lineage development. To begin to understand the structural and functional complexity of enhancers that regulate ordered NB gene expression programs, we have undertaken the mutational analysis of the temporally restricted nerfin-1 NB enhancer. Our previous studies have localized the enhancer to a region just proximal to the nerfin-1 transcription start site. Analysis of this enhancer, using the phylogenetic footprint program EvoPrinter, reveals the presence of multiple sequence blocks that are conserved among drosophilids. cis-Decoder alignments of these conserved sequence blocks (CSBs) has identified shorter elements that are conserved in other Drosophila NB enhancers. Mutagenesis of the enhancer reveals that although each CSB is required for wild-type expression, neither position nor orientation of the CSBs within the enhancer is crucial for enhancer function; removal of less-conserved or non-conserved sequences flanking CSB clusters also does not significantly alter enhancer activity. While all three conserved E-box transcription factor (TF) binding sites (CAGCTG) are required for full function, adding an additional site at different locations within non-conserved sequences interferes with enhancer activity. Of particular note, none of the mutations resulted in ectopic reporter expression outside of the early NB expression window, suggesting that the temporally restricted pattern is defined by transcriptional activators and not by direct DNA binding repressors. Our work also points to an unexpectedly large number of TFs required for optimal enhancer function - mutant TF analysis has identified at least four that are required for full enhancer regulation.

Common motifs shared by conserved enhancers of Drosophila midline glial genes

Coding sequences are usually the most highly conserved sectors of DNA, but genomic regions controlling the expression pattern of certain genes can also be conserved across diverse species. In this study, we identify five enhancers capable of activating transcription in the midline glia of Drosophila melanogaster and each contains sequences conserved across at least 11 Drosophila species. In addition, the conserved sequences contain reiterated motifs for binding sites of the known midline transcriptional activators, Single-minded, Tango, Dichaete, and Pointed. To understand the molecular basis for the highly conserved genomic subregions within enhancers of the midline genes, we tested the ability of various motifs to affect midline expression, both individually and in combination, within synthetic reporter constructs. Multiple copies of the binding site for the midline regulators Single-minded and Tango can drive expression in midline cells; however, small changes to the sequences flanking this transcription factor binding site can inactivate expression in midline cells and activate expression in tracheal cells instead. For the midline genes described in this study, the highly conserved sequences appear to juxtapose positive and negative regulatory factors in a configuration that activates genes specifically in the midline glia, while maintaining them inactive in other tissues, including midline neurons and tracheal cells.

Expression and regulation of vitamin D3 upregulated protein 1 (VDUP1) is conserved in mammalian and insect brain

Originally characterized as a cell-cycle inhibitor induced by vitamin D(3), the tumor suppressor vitamin-D(3) upregulated protein 1 (VDUP1) has increasingly been shown to play major physiological roles in cell differentiation and glucose metabolism. Here we show evolutionarily conserved expression patterns of VDUP1 in Drosophila and rat nervous systems, including subcellular localization--cytoplasmic enrichment in neurons and nuclear expression in glia. These anatomical correlates suggested conservation of VDUP1 regulation, which was investigated both functionally and through promoter studies. Characterization of orthologous vdup1 cis-regulatory regions identified evolutionarily conserved sequence blocks (CSBs) with similarities to neural enhancers, including basic helix-loop-helix (bHLH) transcription factor Neurogenin/Math/atonal and Mash/achaete-scute family members. E-boxes (CANNTG), the binding sites for bHLH proteins, were associated with these CSBs as well, including E-boxes known to mediate glucose-dependent upregulation of VDUP1 in nonneuronal cells. Hyperglycemia-induced upregulation of VDUP1 was observed in brain tumor cells and in the Drosophila nervous system, which resulted in developmental arrest. Taken together, these data demonstrate evolutionary conservation of VDUP1 regulation and function, and suggest an expanding role for VDUP1 in nervous system development.

Hedgehog-dependent down-regulation of the tumor suppressor, vitamin D3 up-regulated protein 1 (VDUP1), precedes lamina development in Drosophila

The tumor suppressor vitamin D(3) up-regulated protein 1 (VDUP1) is expressed throughout the developing and mature Drosophila nervous system, but its regulatory pathways are not well understood. Within the developing Drosophila visual system, down-regulation of VDUP1 in lamina precursor cells (LPCs) coincided with the arrival of retinal axons into the lamina target field, suggesting VDUP1 regulation by an axonally transmitted signal. Hedgehog (Hh) is a signal well known to coordinate LPC proliferation and differentiation in response to retinal axon innervation, and analysis of orthologous dvdup1 promoters identified an evolutionarily conserved binding site for the Hh-dependent transcription factor cubitus interruptus (Ci). Hh-dependent regulation of VDUP1 in the developing lamina was confirmed in Hh loss-of-function backgrounds where VDUP1 expression was maintained in LPCs, inhibiting both cell proliferation and lamina neurogenesis. This putative coupling of VDUP1 to the Hh signaling pathway may provide novel insights into the mechanisms controlling brain growth and development.

Conserved sequence block clustering and flanking inter-cluster flexibility delineate enhancers that regulate nerfin-1 expression during Drosophila CNS development

We have identified clusters of conserved sequences constituting discrete modular enhancers within the Drosophilanerfin-1 locus. nerfin-1 encodes a Zn-finger transcription factor that directs pioneer interneuron axon guidance. nerfin-1 mRNA is detected in many early delaminating neuroblasts, ganglion mother cells and transiently in nascent neurons. The comparative genomics analysis program EvoPrinter revealed conserved sequence blocks both upstream and downstream of the transcribed region. By using the aligning regions of different drosophilids as the reference DNA, EvoPrinter detects sequence length flexibility between clusters of conserved sequences and thus facilitates differentiation between closely associated modular enhancers. Expression analysis of enhancer-reporter transgenes identified enhancers that drive expression in different regions of the developing embryonic and adult nervous system, including subsets of embryonic CNS neuroblasts, GMCs, neurons and PNS neurons. In summary, EvoPrinter facilitates the discovery and analysis of enhancers that control crucial aspects of nerfin-1 expression.

Sequence conservation and combinatorial complexity of Drosophila neural precursor cell enhancers

Background

The presence of highly conserved sequences within cis-regulatory regions can serve as a valuable starting point for elucidating the basis of enhancer function. This study focuses on regulation of gene expression during the early events of Drosophila neural development. We describe the use of EvoPrinter and cis-Decoder, a suite of interrelated phylogenetic footprinting and alignment programs, to characterize highly conserved sequences that are shared among co-regulating enhancers.

Results

Analysis of in vivo characterized enhancers that drive neural precursor gene expression has revealed that they contain clusters of highly conserved sequence blocks (CSBs) made up of shorter shared sequence elements which are present in different combinations and orientations within the different co-regulating enhancers; these elements contain either known consensus transcription factor binding sites or consist of novel sequences that have not been functionally characterized. The CSBs of co-regulated enhancers share a large number of sequence elements, suggesting that a diverse repertoire of transcription factors may interact in a highly combinatorial fashion to coordinately regulate gene expression. We have used information gained from our comparative analysis to discover an enhancer that directs expression of the nervy gene in neural precursor cells of the CNS and PNS.

Conclusion

The combined use EvoPrinter and cis-Decoder has yielded important insights into the combinatorial appearance of fundamental sequence elements required for neural enhancer function. Each of the 30 enhancers examined conformed to a pattern of highly conserved blocks of sequences containing shared constituent elements. These data establish a basis for further analysis and understanding of neural enhancer function.

Rapid detection and curation of conserved DNA via enhanced-BLAT and EvoPrinterHD analysis

Background

Multi-genome comparative analysis has yielded important insights into the molecular details of gene regulation. We have developed EvoPrinter, a web-accessed genomics tool that provides a single uninterrupted view of conserved sequences as they appear in a species of interest. An EvoPrint reveals with near base-pair resolution those sequences that are essential for gene function.

Results

We describe here EvoPrinterHD, a 2^nd-generation comparative genomics tool that automatically generates from a single input sequence an enhanced view of sequence conservation between evolutionarily distant species. Currently available for 5 nematode, 3 mosquito, 12 Drosophila, 20 vertebrate, 17 Staphylococcus and 20 enteric bacteria genomes, EvoPrinterHD employs a modified BLAT algorithm [enhanced-BLAT (eBLAT)], which detects up to 75% more conserved bases than identified by the BLAT alignments used in the earlier EvoPrinter program. The new program also identifies conserved sequences within rearranged DNA, highlights repetitive DNA, and detects sequencing gaps. EvoPrinterHD currently holds over 112 billion bp of indexed genomes in memory and has the flexibility of selecting a subset of genomes for analysis. An EvoDifferences profile is also generated to portray conserved sequences that are uniquely lost in any one of the orthologs. Finally, EvoPrinterHD incorporates options that allow for (1) re-initiation of the analysis using a different genome's aligning region as the reference DNA to detect species-specific changes in less-conserved regions, (2) rapid extraction and curation of conserved sequences, and (3) for bacteria, identifies unique or uniquely shared sequences present in subsets of genomes.

Conclusion

EvoPrinterHD is a fast, high-resolution comparative genomics tool that automatically generates an uninterrupted species-centric view of sequence conservation and enables the discovery of conserved sequences within rearranged DNA. When combined with cis-Decoder, a program that discovers sequence elements shared among tissue specific enhancers, EvoPrinterHD facilitates the analysis of conserved sequences that are essential for coordinate gene regulation.

Phylogenetic footprinting analysis in the upstream regulatory regions of the Drosophila enhancer of split genes

During Drosophila development Suppressor of Hairless [Su(H)]-dependent Notch activation upregulates transcription of the Enhancer of split-Complex [E(spl)-C] genes. Drosophila melanogaster E(spl) genes share common transcription regulators including binding sites for Su(H), proneural, and E(spl) basic-helix-loop-helix (bHLH) proteins. However, the expression patterns of E(spl) genes during development suggest that additional factors are involved. To better understand regulators responsible for these expression patterns, recently available sequence and annotation data for multiple Drosophila genomes were used to compare the E(spl) upstream regulatory regions from more than nine Drosophila species. The mgamma and mbeta regulatory regions are the most conserved of the bHLH genes. Fine analysis of Su(H) sites showed that high-affinity Su(H) paired sites and the Su(H) paired site plus proneural site (SPS + A) architecture are completely conserved in a subset of Drosophila E(spl) genes. The SPS + A module is also present in the upstream regulatory regions of the more ancient mosquito and honeybee E(spl) bHLH genes. Additional transcription factor binding sites were identified upstream of the E(spl) genes and compared between species of Drosophila. Conserved sites provide new understandings about E(spl) regulation during development. Conserved novel sequences found upstream of multiple E(spl) genes may play a role in the expression of these genes.

REDfly 2.0: an integrated database of cis-regulatory modules and transcription factor binding sites in Drosophila

The identification and study of the cis-regulatory elements that control gene expression are important areas of biological research, but few resources exist to facilitate large-scale bioinformatics studies of cis-regulation in metazoan species. Drosophila melanogaster, with its well-annotated genome, exceptional resources for comparative genomics and long history of experimental studies of transcriptional regulation, represents the ideal system for regulatory bioinformatics. We have merged two existing Drosophila resources, the REDfly database of cis-regulatory modules and the FlyReg database of transcription factor binding sites (TFBSs), into a single integrated database containing extensive annotation of empirically validated cis-regulatory modules and their constituent binding sites. With the enhanced functionality made possible through this integration of TFBS data into REDfly, together with additional improvements to the REDfly infrastructure, we have constructed a one-stop portal for Drosophila cis-regulatory data that will serve as a powerful resource for both computational and experimental studies of transcriptional regulation. REDfly is freely accessible at http://redfly.ccr.buffalo.edu.