Genomic Takeover by Transposable Elements in the Strawberry Poison Frog

Abstract We sequenced the genome of the strawberry poison frog, Oophaga pumilio, at a depth of 127.5× using variable insert size libraries. The total genome size is estimated to be 6.76 Gb, of which 4.76 Gb are from high copy number repetitive elements with low differentiation across copies. These repeats encompass DNA transposons, RNA transposons, and LTR retrotransposons, including at least 0.4 and 1.0 Gb of Mariner/Tc1 and Gypsy elements, respectively. Expression data indicate high levels of gypsy and Mariner/Tc1 expression in ova of O. pumilio compared with Xenopus laevis. We further observe phylogenetic evidence for horizontal transfer (HT) of Mariner elements, possibly between fish and frogs. The elements affected by HT are present in high copy number and are highly expressed, suggesting ongoing proliferation after HT. Our results suggest that the large amphibian genome sizes, at least partially, can be explained by a process of repeated invasion of new transposable elements that are not yet suppressed in the germline. We also find changes in the spliceosome that we hypothesize are related to permissiveness of O. pumilio to increases in intron length due to transposon proliferation. Finally, we identify the complement of ion channels in the first genomic sequenced poison frog and discuss its relation to the evolution of autoresistance to toxins sequestered in the skin.

[1]  R. Shine,et al.  Improving amphibian genomic resources: a multitissue reference transcriptome of an iconic invader , 2017, GigaScience.

[2]  R. Mueller piRNAs and Evolutionary Trajectories in Genome Size and Content , 2017, Journal of Molecular Evolution.

[3]  E. Betrán,et al.  Transposable Element Domestication As an Adaptation to Evolutionary Conflicts. , 2017, Trends in genetics : TIG.

[4]  B. Venkatesh,et al.  Voltage-gated sodium channel gene repertoire of lampreys: gene duplications, tissue-specific expression and discovery of a long-lost gene , 2017, Proceedings of the Royal Society B: Biological Sciences.

[5]  R. Harris,et al.  Interacting amino acid replacements allow poison frogs to evolve epibatidine resistance , 2017, Science.

[6]  G. Wang,et al.  Single rat muscle Na+ channel mutation confers batrachotoxin autoresistance found in poison-dart frog Phyllobates terribilis , 2017, Proceedings of the National Academy of Sciences.

[7]  C. Feschotte,et al.  Regulatory activities of transposable elements: from conflicts to benefits , 2016, Nature Reviews Genetics.

[8]  R. McLaughlin,et al.  Genetic conflicts: the usual suspects and beyond , 2017, Journal of Experimental Biology.

[9]  C. Richards‐Zawacki,et al.  Poison frog color morphs express assortative mate preferences in allopatry but not sympatry , 2016, Evolution; international journal of organic evolution.

[10]  R. Murphy,et al.  Contribution of Multiple Inter-Kingdom Horizontal Gene Transfers to Evolution and Adaptation of Amphibian-Killing Chytrid, Batrachochytrium dendrobatidis , 2016, Front. Microbiol..

[11]  Yufeng Wu,et al.  REPdenovo: Inferring De Novo Repeat Motifs from Short Sequence Reads , 2016, PloS one.

[12]  D. Cannatella,et al.  Convergent Substitutions in a Sodium Channel Suggest Multiple Origins of Toxin Resistance in Poison Frogs. , 2016, Molecular biology and evolution.

[13]  Erik Kaestner,et al.  The Origins Of Genome Architecture , 2016 .

[14]  Juan C. Santos,et al.  A Review of Chemical Defense in Poison Frogs (Dendrobatidae): Ecology, Pharmacokinetics, and Autoresistance , 2016 .

[15]  O. Kohany,et al.  Repbase Update, a database of repetitive elements in eukaryotic genomes , 2015, Mobile DNA.

[16]  Jacob S. Vestergaard,et al.  Number of genes controlling a quantitative trait in a hybrid zone of the aposematic frog Ranitomeya imitator , 2015, Proceedings of the Royal Society B: Biological Sciences.

[17]  Jun Wang,et al.  Whole-genome sequence of the Tibetan frog Nanorana parkeri and the comparative evolution of tetrapod genomes , 2015, Proceedings of the National Academy of Sciences.

[18]  Cédric Feschotte,et al.  Ancient Transposable Elements Transformed the Uterine Regulatory Landscape and Transcriptome during the Evolution of Mammalian Pregnancy , 2015, Cell reports.

[19]  R. Mueller,et al.  Low Levels of LTR Retrotransposon Deletion by Ectopic Recombination in the Gigantic Genomes of Salamanders , 2015, Journal of Molecular Evolution.

[20]  Jürgen Gadau,et al.  Transposable element islands facilitate adaptation to novel environments in an invasive species , 2014, Nature Communications.

[21]  P. Davis,et al.  De novo Assembly and Analysis of the Northern Leopard Frog Rana pipiens Transcriptome , 2014, Journal of genomics.

[22]  S. Edwards,et al.  Parallel Evolution of Tetrodotoxin Resistance in Three Voltage-Gated Sodium Channel Genes in the Garter Snake Thamnophis sirtalis , 2014, Molecular biology and evolution.

[23]  C. Dulac,et al.  Neural control of maternal and paternal behaviors , 2014, Science.

[24]  P. Stynoski,et al.  Maternally derived chemical defences are an effective deterrent against some predators of poison frog tadpoles (Oophaga pumilio) , 2014, Biology Letters.

[25]  J. Bennetzen,et al.  The contributions of transposable elements to the structure, function, and evolution of plant genomes. , 2014, Annual review of plant biology.

[26]  J. Stynoski,et al.  Evidence of maternal provisioning of alkaloid-based chemical defenses in the strawberry poison frog Oophaga pumilio. , 2014, Ecology.

[27]  A. Gregory Matera,et al.  A day in the life of the spliceosome , 2014, Nature Reviews Molecular Cell Biology.

[28]  Xun Xu,et al.  SOAPdenovo-Trans: de novo transcriptome assembly with short RNA-Seq reads , 2013, Bioinform..

[29]  Josefa González,et al.  The impact of transposable elements in environmental adaptation , 2013, Molecular ecology.

[30]  K. Katoh,et al.  MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability , 2013, Molecular biology and evolution.

[31]  R. Mueller,et al.  Slow DNA Loss in the Gigantic Genomes of Salamanders , 2012, Genome biology and evolution.

[32]  K. Summers,et al.  Mate choice and the genetic basis for colour variation in a polymorphic dart frog: inferences from a wild pedigree , 2012, Molecular ecology.

[33]  Maxim Teslenko,et al.  MrBayes 3.2: Efficient Bayesian Phylogenetic Inference and Model Choice Across a Large Model Space , 2012, Systematic biology.

[34]  M. A. Donnelly,et al.  A review of chemical ecology in poison frogs , 2012, Chemoecology.

[35]  Francis M. Jiggins,et al.  Successive Increases in the Resistance of Drosophila to Viral Infection through a Transposon Insertion Followed by a Duplication , 2011, PLoS genetics.

[36]  N. Friedman,et al.  Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data , 2011, Nature Biotechnology.

[37]  H. Zakon,et al.  Expansion of voltage-dependent Na+ channel gene family in early tetrapods coincided with the emergence of terrestriality and increased brain complexity. , 2011, Molecular biology and evolution.

[38]  Ramón Doallo,et al.  ProtTest 3: fast selection of best-fit models of protein evolution , 2011, Bioinform..

[39]  Qin Ding,et al.  AlienG: An Effective Computational Tool for Phylogenomic Identification of Horizontally Transferred Genes , 2011, Bioinformatics and Computational Biology.

[40]  Cédric Feschotte,et al.  Promiscuous DNA: horizontal transfer of transposable elements and why it matters for eukaryotic evolution. , 2010, Trends in ecology & evolution.

[41]  O. Gascuel,et al.  New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. , 2010, Systematic biology.

[42]  J. Daly,et al.  Roughing it: a mantellid poison frog shows greater alkaloid diversity in some disturbed habitats. , 2010, Journal of natural products.

[43]  Y. Nomura,et al.  Functional Expression of an Arachnid Sodium Channel Reveals Residues Responsible for Tetrodotoxin Resistance in Invertebrate Sodium Channels* , 2009, The Journal of Biological Chemistry.

[44]  D. Hillis,et al.  Toxin-resistant sodium channels: parallel adaptive evolution across a complete gene family. , 2008, Molecular biology and evolution.

[45]  P. Capy,et al.  Revisiting horizontal transfer of transposable elements in Drosophila , 2008, Heredity.

[46]  Ian A. Swinburne,et al.  Intron delays and transcriptional timing during development. , 2008, Developmental cell.

[47]  J. Daly,et al.  Spatial and temporal patterns of alkaloid variation in the poison frog Oophaga pumilio in Costa Rica and Panama over 30 years. , 2007, Toxicon : official journal of the International Society on Toxinology.

[48]  Rodrigo Lopez,et al.  Clustal W and Clustal X version 2.0 , 2007, Bioinform..

[49]  Ziheng Yang PAML 4: phylogenetic analysis by maximum likelihood. , 2007, Molecular biology and evolution.

[50]  Jane Mitchell,et al.  Irreversible Block of Cardiac Mutant Na+ Channels by Batrachotoxin , 2007, Channels.

[51]  Jane Mitchell,et al.  Serine-401 as a batrachotoxin- and local anesthetic-sensing residue in the human cardiac Na+ channel , 2007, Pflügers Archiv - European Journal of Physiology.

[52]  Jane Mitchell,et al.  How Batrachotoxin Modifies the Sodium Channel Permeation Pathway: Computer Modeling and Site-Directed Mutagenesis , 2006, Molecular Pharmacology.

[53]  J. Inoue,et al.  The mitochondrial genome of spotted green pufferfish Tetraodon nigroviridis (Teleostei: Tetraodontiformes) and divergence time estimation among model organisms in fishes. , 2006, Genes & genetic systems.

[54]  Burkhard Morgenstern,et al.  Gene prediction in eukaryotes with a generalized hidden Markov model that uses hints from external sources , 2006, BMC Bioinformatics.

[55]  B. Birren,et al.  Structure and Architecture of the Maize Genome1[W] , 2005, Plant Physiology.

[56]  J. Daly,et al.  Alkaloids from amphibian skin: a tabulation of over eight-hundred compounds. , 2005, Journal of natural products.

[57]  D. Petrov,et al.  References and Notes Materials and Methods Tables S1 and S2 References and Notes Pesticide Resistance via Transposition-mediated Adaptive Gene Truncation in Drosophila , 2022 .

[58]  J. Jurka,et al.  Repbase Update, a database of eukaryotic repetitive elements , 2005, Cytogenetic and Genome Research.

[59]  Marc A. Schaub,et al.  Subdivision of Large Introns in Drosophila by Recursive Splicing at Nonexonic Elements , 2005, Genetics.

[60]  P. Ruben,et al.  Evolutionary diversification of TTX-resistant sodium channels in a predator–prey interaction , 2005, Nature.

[61]  J. Daly,et al.  Dietary source for skin alkaloids of poison frogs (Dendrobatidae)? , 1994, Journal of Chemical Ecology.

[62]  Pavel A. Pevzner,et al.  De novo identification of repeat families in large genomes , 2005, ISMB.

[63]  Burkhard Morgenstern,et al.  AUGUSTUS: a web server for gene finding in eukaryotes , 2004, Nucleic Acids Res..

[64]  P. Weygoldt Complex brood care and reproductive behaviour in captive poison-arrow frogs, Dendrobates pumilio O. Schmidt , 1980, Behavioral Ecology and Sociobiology.

[65]  A. J. Crawford,et al.  Relative Rates of Nucleotide Substitution in Frogs , 2003, Journal of Molecular Evolution.

[66]  Mario Stanke,et al.  Gene prediction with a hidden Markov model and a new intron submodel , 2003, ECCB.

[67]  D. Ragsdale,et al.  The batrachotoxin receptor on the voltage-gated sodium channel is guarded by the channel activation gate. , 2002, Molecular pharmacology.

[68]  International Human Genome Sequencing Consortium Initial sequencing and analysis of the human genome , 2001, Nature.

[69]  K. Summers,et al.  Visual mate choice in poison frogs , 1999, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[70]  M. Joron,et al.  Diversity in mimicry: paradox or paradigm? , 1998, Trends in ecology & evolution.

[71]  J. Thompson,et al.  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. , 1994, Nucleic acids research.

[72]  F. Gusovsky,et al.  First occurrence of tetrodotoxin in a dendrobatid frog (Colostethus inguinalis), with further reports for the bufonid genus Atelopus. , 1994, Toxicon : official journal of the International Society on Toxinology.

[73]  I. Mattaj,et al.  The determinants for Sm protein binding to Xenopus U1 and U5 snRNAs are complex and non‐identical. , 1993, The EMBO journal.

[74]  S. K. Sessions,et al.  Amphibian cytogenetics and evolution , 1992 .

[75]  C. Langley,et al.  Chromosome rearrangement by ectopic recombination in Drosophila melanogaster: genome structure and evolution. , 1991, Genetics.

[76]  J. Bogart CHAPTER 11 – The Influence of Life History on Karyotypic Evolution in Frogs , 1991 .

[77]  L. Lowcock,et al.  Genome size and metabolic rate in salamanders , 1991 .

[78]  E. Lund,et al.  The transcription of Xenopus laevis embryonic U1 snRNA genes changes when oocytes mature into eggs. , 1987, Genes & development.

[79]  E. Albuquerque,et al.  Levels of batrachotoxin and lack of sensitivity to its action in poison-dart frogs (Phyllobates). , 1980, Science.

[80]  E. Albuquerque,et al.  The pharmacology of batrachotoxin. V. A comparative study of membrane properties and the effect of batrachotoxin on sartorius muscles of the frogs Phyllobates aurotaenia and Rana pipiens. , 1973, The Journal of pharmacology and experimental therapeutics.

[81]  J. Daly,et al.  The structure of batrachotoxin, a steroidal alkaloid from the Colombian arrow poison frog, Phyllobates aurotaenia, and partial synthesis of batrachotoxin and its analogs and homologs. , 1969, Journal of the American Chemical Society.

[82]  J. Daly,et al.  Toxicity of Panamanian Poison Frogs (Dendrobates): Some Biological and Chemical Aspects , 1967, Science.