Combinatorial Gene Regulatory Functions Underlie Ultraconserved Elements in Drosophila

Ultraconserved elements (UCEs) are discrete genomic elements conserved across large evolutionary distances. Although UCEs have been linked to multiple facets of mammalian gene regulation their extreme evolutionary conservation remains largely unexplained. Here, we apply a computational approach to investigate this question in Drosophila, exploring the molecular functions of more than 1,500 UCEs shared across the genomes of 12 Drosophila species. Our data indicate that Drosophila UCEs are hubs for gene regulatory functions and suggest that UCE sequence invariance originates from their combinatorial roles in gene control. We also note that the gene regulatory roles of intronic and intergenic UCEs (iUCEs) are distinct from those found in exonic UCEs (eUCEs). In iUCEs, transcription factor (TF) and epigenetic factor binding data strongly support iUCE roles in transcriptional and epigenetic regulation. In contrast, analyses of eUCEs indicate that they are two orders of magnitude more likely than the expected to simultaneously include protein-coding sequence, TF-binding sites, splice sites, and RNA editing sites but have reduced roles in transcriptional or epigenetic regulation. Furthermore, we use a Drosophila cell culture system and transgenic Drosophila embryos to validate the notion of UCE combinatorial regulatory roles using an eUCE within the Hox gene Ultrabithorax and show that its protein-coding region also contains alternative splicing regulatory information. Taken together our experiments indicate that UCEs emerge as a result of combinatorial gene regulatory roles and highlight common features in mammalian and insect UCEs implying that similar processes might underlie ultraconservation in diverse animal taxa.

[1]  D. Haussler,et al.  Ultraconserved Elements in the Human Genome , 2004, Science.

[2]  G. Technau,et al.  The RNA-binding protein ELAV regulates Hox RNA processing, expression and function within the Drosophila nervous system , 2014, Development.

[3]  Laurent Excoffier,et al.  Conserved noncoding sequences are selectively constrained and not mutation cold spots , 2006, Nature Genetics.

[4]  Jürg Müller,et al.  Histone trimethylation and the maintenance of transcriptional ON and OFF states by trxG and PcG proteins. , 2006, Genes & development.

[5]  Alexander van Oudenaarden,et al.  Highly expressed loci are vulnerable to misleading ChIP localization of multiple unrelated proteins , 2013, Proceedings of the National Academy of Sciences.

[6]  A. Kern,et al.  Highly Constrained Intergenic Drosophila Ultraconserved Elements Are Candidate ncRNAs , 2015, Genome biology and evolution.

[7]  M. Ringnér,et al.  An integrative genomics screen uncovers ncRNA T-UCR functions in neuroblastoma tumours , 2010, Oncogene.

[8]  Cole Trapnell,et al.  Ultrafast and memory-efficient alignment of short DNA sequences to the human genome , 2009, Genome Biology.

[9]  M. Gerstein,et al.  Unlocking the secrets of the genome , 2009, Nature.

[10]  Tyson A. Clark,et al.  Ultraconserved elements are associated with homeostatic control of splicing regulators by alternative splicing and nonsense-mediated decay. , 2007, Genes & development.

[11]  Gautier Koscielny,et al.  VectorBase: improvements to a bioinformatics resource for invertebrate vector genomics , 2011, Nucleic Acids Res..

[12]  Michael Pheasant,et al.  Comparison of Ultra-Conserved Elements in Drosophilids and Vertebrates , 2013, PloS one.

[13]  Thomas D. Schmittgen,et al.  Ultraconserved regions encoding ncRNAs are altered in human leukemias and carcinomas. , 2007, Cancer cell.

[14]  Alan M. Moses,et al.  In vivo enhancer analysis of human conserved non-coding sequences , 2006, Nature.

[15]  Aaron R. Quinlan,et al.  Bioinformatics Applications Note Genome Analysis Bedtools: a Flexible Suite of Utilities for Comparing Genomic Features , 2022 .

[16]  S. Brenner,et al.  Unproductive splicing of SR genes associated with highly conserved and ultraconserved DNA elements , 2007, Nature.

[17]  M. Nóbrega,et al.  Comparative genomics at the vertebrate extremes , 2004, Nature Reviews Genetics.

[18]  Robert C. Edgar,et al.  MUSCLE: multiple sequence alignment with high accuracy and high throughput. , 2004, Nucleic acids research.

[19]  Manolis Kellis,et al.  Locating protein-coding sequences under selection for additional, overlapping functions in 29 mammalian genomes. , 2011, Genome research.

[20]  Loqmane Seridi,et al.  The evolution of ultraconserved elements with different phylogenetic origins , 2012, BMC Evolutionary Biology.

[21]  Jin Billy Li,et al.  RADAR: a rigorously annotated database of A-to-I RNA editing , 2013, Nucleic Acids Res..

[22]  R. Maeda,et al.  An optimized transgenesis system for Drosophila using germ-line-specific φC31 integrases , 2007, Proceedings of the National Academy of Sciences.

[23]  A. J. Lopez,et al.  Generation of alternative Ultrabithorax isoforms and stepwise removal of a large intron by resplicing at exon-exon junctions. , 1998, Molecular cell.

[24]  J. Shendure,et al.  Systematic Dissection of Coding Exons at Single Nucleotide Resolution Supports an Additional Role in Cell-Specific Transcriptional Regulation , 2014, PLoS genetics.

[25]  V. Pirrotta,et al.  Extensive conservation of sequences and chromatin structures in the bxd polycomb response element among Drosophilid species. , 2002, International Journal of Developmental Biology.

[26]  Israel Steinfeld,et al.  BMC Bioinformatics BioMed Central , 2008 .

[27]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[28]  Timothy B Sackton,et al.  Genome of the house fly, Musca domestica L., a global vector of diseases with adaptations to a septic environment , 2014, Genome Biology.

[29]  B. Cohen,et al.  The strength of selection on ultraconserved elements in the human genome. , 2007, American journal of human genetics.

[30]  M. Esteller,et al.  Regulation of pri-miRNA processing by a long noncoding RNA transcribed from an ultraconserved region. , 2014, Molecular cell.

[31]  C. Shyu,et al.  Long identical multispecies elements in plant and animal genomes , 2012, Proceedings of the National Academy of Sciences.

[32]  A. J. Lopez,et al.  Evolutionary conservation of the structure and expression of alternatively spliced Ultrabithorax isoforms from Drosophila. , 1994, Genetics.

[33]  A. Kern,et al.  Computational analysis and characterization of UCE-like elements (ULEs) in plant genomes , 2012, Genome research.

[34]  Timothy L. Bailey,et al.  Gene expression Advance Access publication May 4, 2011 DREME: motif discovery in transcription factor ChIP-seq data , 2011 .

[35]  Devanand S. Manoli,et al.  Male-specific fruitless specifies the neural substrates of Drosophila courtship behaviour , 2005, Nature.

[36]  M. Nei,et al.  Ultraconserved coding regions outside the homeobox of mammalian Hox genes , 2008, BMC Evolutionary Biology.

[37]  Axel Visel,et al.  Deletion of Ultraconserved Elements Yields Viable Mice , 2007, PLoS biology.

[38]  G. Azzam,et al.  Developmental RNA processing of 3′UTRs in Hox mRNAs as a context-dependent mechanism modulating visibility to microRNAs , 2010, Development.

[39]  Geoffrey H. Siwo,et al.  Genome Sequence of the Tsetse Fly (Glossina morsitans): Vector of African Trypanosomiasis , 2014, Science.

[40]  Ziheng Yang,et al.  PAML: a program package for phylogenetic analysis by maximum likelihood , 1997, Comput. Appl. Biosci..

[41]  Brian S. Clark,et al.  The Evf-2 noncoding RNA is transcribed from the Dlx-5/6 ultraconserved region and functions as a Dlx-2 transcriptional coactivator. , 2006, Genes & development.

[42]  N. Perrimon,et al.  Targeted gene expression as a means of altering cell fates and generating dominant phenotypes. , 1993, Development.

[43]  Tom H. Pringle,et al.  The human genome browser at UCSC. , 2002, Genome research.

[44]  J. Tena,et al.  A functional survey of the enhancer activity of conserved non-coding sequences from vertebrate Iroquois cluster gene deserts. , 2005, Genome research.

[45]  D. Haussler,et al.  Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. , 2005, Genome research.

[46]  L. Pease,et al.  Gene splicing and mutagenesis by PCR-driven overlap extension , 2007, Nature Protocols.

[47]  Boris Lenhard,et al.  The mystery of extreme non-coding conservation , 2013, Philosophical Transactions of the Royal Society B: Biological Sciences.

[48]  I. Makunin,et al.  Underreplicated Regions in Drosophila melanogaster Are Enriched with Fast-Evolving Genes and Highly Conserved Noncoding Sequences , 2014, Genome biology and evolution.

[49]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[50]  J. Mattick,et al.  Evidence for Control of Splicing by Alternative RNA Secondary Structures in Dipteran Homothorax Pre-mRNA , 2006, RNA biology.

[51]  D. Hogness,et al.  Immunochemical dissection of the Ultrabithorax homeoprotein family in Drosophila melanogaster. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[52]  Colin N. Dewey,et al.  Discovery of functional elements in 12 Drosophila genomes using evolutionary signatures , 2007, Nature.

[53]  Robert A. H. White,et al.  Alternative Splicing Modulates Ubx Protein Function in Drosophila melanogaster , 2010, Genetics.

[54]  B. Birren,et al.  Sequencing and comparison of yeast species to identify genes and regulatory elements , 2003, Nature.

[55]  Klaudia Walter,et al.  Highly Conserved Non-Coding Sequences Are Associated with Vertebrate Development , 2004, PLoS biology.

[56]  Boris Lenhard,et al.  Arrays of ultraconserved non-coding regions span the loci of key developmental genes in vertebrate genomes , 2004, BMC Genomics.

[57]  Obi L. Griffith,et al.  ORegAnno: an open-access community-driven resource for regulatory annotation , 2007, Nucleic Acids Res..

[58]  W. Gehring,et al.  Isolation of caudal, a Drosophila homeo box‐containing gene with maternal expression, whose transcripts form a concentration gradient at the pre‐blastoderm stage , 1985, The EMBO journal.

[59]  E. Demerath,et al.  Ultraconserved Elements in the Human Genome: Association and Transmission Analyses of Highly Constrained Single-Nucleotide Polymorphisms , 2012, Genetics.

[60]  A. Visel,et al.  Ultraconservation identifies a small subset of extremely constrained developmental enhancers , 2008, Nature Genetics.

[61]  Alex P. Reynolds,et al.  Exonic Transcription Factor Binding Directs Codon Choice and Affects Protein Evolution , 2013, Science.

[62]  T. Ørntoft,et al.  Identification of expressed and conserved human noncoding RNAs , 2014, RNA.

[63]  C. Nüsslein-Volhard,et al.  Mutations affecting segment number and polarity in Drosophila , 1980, Nature.

[64]  L. Ringrose,et al.  What are memories made of? How Polycomb and Trithorax proteins mediate epigenetic memory , 2014, Nature Reviews Molecular Cell Biology.

[65]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[66]  P. Radice,et al.  SNPs in ultraconserved elements and familial breast cancer risk. , 2009, Carcinogenesis.

[67]  Gill Bejerano,et al.  Ultraconserved elements in insect genomes: a highly conserved intronic sequence implicated in the control of homothorax mRNA splicing. , 2005, Genome research.

[68]  Stephen M. Mount,et al.  The genome sequence of Drosophila melanogaster. , 2000, Science.

[69]  M. Heller,et al.  Systematic Functional Analysis of Bicaudal-D Serine Phosphorylation and Intragenic Suppression of a Female Sterile Allele of BicD , 2009, PloS one.

[70]  Angela N. Brooks,et al.  Identification and experimental validation of splicing regulatory elements in Drosophila melanogaster reveals functionally conserved splicing enhancers in metazoans. , 2011, RNA.

[71]  M. Ekker,et al.  An SNP in an ultraconserved regulatory element affects Dlx5/Dlx6 regulation in the forebrain , 2010, Development.

[72]  Alessandro Vullo,et al.  Ensembl 2015 , 2014, Nucleic Acids Res..

[73]  E. Sánchez-Herrero,et al.  Integration of RNA processing and expression level control modulates the function of the Drosophila Hox gene Ultrabithorax during adult development , 2011, Development.

[74]  Melanie A. Huntley,et al.  Evolution of genes and genomes on the Drosophila phylogeny , 2007, Nature.

[75]  Matthias Mann,et al.  A DNA-centric protein interaction map of ultraconserved elements reveals contribution of transcription factor binding hubs to conservation. , 2013, Cell reports.

[76]  T Marty,et al.  Regulation of Hox target genes by a DNA bound Homothorax/Hox/Extradenticle complex. , 1999, Development.

[77]  Inna Dubchak,et al.  Comparative genome sequencing of Drosophila pseudoobscura: chromosomal, gene, and cis-element evolution. , 2005, Genome research.

[78]  David E. Leib,et al.  Evf2 (Dlx6as) lncRNA regulates ultraconserved enhancer methylation and the differential transcriptional control of adjacent genes , 2013, Development.

[79]  R. Rezsohazy,et al.  An ultraconserved Hox–Pbx responsive element resides in the coding sequence of Hoxa2 and is active in rhombomere 4 , 2008, Nucleic acids research.

[80]  Megan F. Cole,et al.  Control of Developmental Regulators by Polycomb in Human Embryonic Stem Cells , 2006, Cell.