Increased Adaptive Variation Despite Reduced Overall Genetic Diversity in a Rapidly Adapting Invader

Invasive species often evolve rapidly following introduction despite genetic bottlenecks that may result from small numbers of founders; however, some invasions may not fit this “genetic paradox”. The invasive cane toad (Rhinella marina) displays high phenotypic variation across its introduced Australian range. Here, we used three genome-wide datasets to characterize their population structure and genetic diversity. We found that toads form three genetic clusters: 1) native range toads, 2) toads from the source population in Hawaii and long-established areas near introduction sites in Australia, and 3) toads from more recently established northern Australian sites. Although we find an overall reduction in genetic diversity following introduction, we do not see this reduction in loci putatively under selection, suggesting that genetic diversity may have been maintained at ecologically relevant traits, or that mutation rates were high enough to maintain adaptive potential. Nonetheless, toads encounter novel environmental challenges in Australia, and the transition between genetic clusters occurs at a point along the invasion transect where temperature rises and rainfall decreases. We identify environmentally associated loci known to be involved in resistance to heat and dehydration. This study highlights that natural selection occurs rapidly and plays a vital role in shaping the structure of invasive populations.

[1]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[2]  Jean Thioulouse,et al.  Interactive Multivariate Data Analysis in R with the ade4 and ade4TkGUI Packages , 2007 .

[3]  J. Murray,et al.  The Evolutionary Process , 1985 .

[4]  W. Sherwin,et al.  Invasive species can't cover their tracks: using microsatellites to assist management of starling (Sturnus vulgaris) populations in Western Australia , 2009, Molecular ecology.

[5]  R. Alford,et al.  Mixed population genomics support for the central marginal hypothesis across the invasive range of the cane toad (Rhinella marina) in Australia , 2016, Molecular ecology.

[6]  M. Blaxter,et al.  RADSeq: next-generation population genetics. , 2010, Briefings in functional genomics.

[7]  R. Shine,et al.  Characterisation of Major Histocompatibility Complex Class I in the Australian Cane Toad, Rhinella marina , 2014, PloS one.

[8]  Angel Amores,et al.  Stacks: an analysis tool set for population genomics , 2013, Molecular ecology.

[9]  J. Kiang,et al.  Heat shock protein 70 kDa: molecular biology, biochemistry, and physiology. , 1998, Pharmacology & therapeutics.

[10]  J. Webb,et al.  Adrenocortical stress responses influence an invasive vertebrate's fitness in an extreme environment , 2013, Proceedings of the Royal Society B: Biological Sciences.

[11]  M. Gerstein,et al.  RNA-Seq: a revolutionary tool for transcriptomics , 2009, Nature Reviews Genetics.

[12]  E. Grayhack,et al.  Synonymous Codons: Choose Wisely for Expression. , 2017, Trends in genetics : TIG.

[13]  S. Sultan,et al.  A resurrection study reveals rapid adaptive evolution within populations of an invasive plant , 2012, Evolutionary applications.

[14]  M. Stephens,et al.  fastSTRUCTURE: Variational Inference of Population Structure in Large SNP Data Sets , 2014, Genetics.

[15]  Rory S Telemeco,et al.  Temperature has species-specific effects on corticosterone in alligator lizards. , 2014, General and comparative endocrinology.

[16]  M. Hellberg,et al.  Host‐targeted RAD‐Seq reveals genetic changes in the coral Oculina patagonica associated with range expansion along the Spanish Mediterranean coast , 2018, Molecular ecology.

[17]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[18]  P. Unmack,et al.  dartr: An r package to facilitate analysis of SNP data generated from reduced representation genome sequencing , 2018, Molecular ecology resources.

[19]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[20]  B. Melbourne,et al.  Rapid adaptive evolution in novel environments acts as an architect of population range expansion , 2017, Proceedings of the National Academy of Sciences.

[21]  Cathy H. Wu,et al.  UniProt: the Universal Protein knowledgebase , 2004, Nucleic Acids Res..

[22]  Thibaut Jombart,et al.  adegenet 1.3-1: new tools for the analysis of genome-wide SNP data , 2011, Bioinform..

[23]  Pablo Cingolani,et al.  © 2012 Landes Bioscience. Do not distribute. , 2022 .

[24]  R. Colautti,et al.  Rapid Adaptation to Climate Facilitates Range Expansion of an Invasive Plant , 2013, Science.

[25]  Frank Wannemaker,et al.  Genetics And Conservation Of Rare Plants , 2016 .

[26]  M. Whitlock,et al.  Reliable Detection of Loci Responsible for Local Adaptation: Inference of a Null Model through Trimming the Distribution of FST* , 2015, The American Naturalist.

[27]  Richard Shine,et al.  Conservation biology: Restoration of an inbred adder population , 1999, Nature.

[28]  Daniel Simberloff,et al.  What are invasive species , 2009 .

[29]  A. Estoup,et al.  Is There a Genetic Paradox of Biological Invasion , 2016 .

[30]  M. Fahim,et al.  Dehydration and food deprivation exacerbate mouse cerebral microvascular responses to local hyperthermia , 1995 .

[31]  F. Valladares,et al.  Invasive species can handle higher leaf temperature under water stress than Mediterranean natives , 2011 .

[32]  Gregory P. Brown,et al.  Locomotor performance of cane toads differs between native-range and invasive populations , 2017, Royal Society Open Science.

[33]  K. Whitney,et al.  Rapid evolution in introduced species, ‘invasive traits’ and recipient communities: challenges for predicting invasive potential , 2008 .

[34]  Cascade J. B. Sorte,et al.  Temperature Tolerance and Stress Proteins as Mechanisms of Invasive Species Success , 2011, PloS one.

[35]  D. H. Reed,et al.  Correlation between Fitness and Genetic Diversity , 2003 .

[36]  Mark R. Christie,et al.  Rapid genetic adaptation to a novel environment despite a genome‐wide reduction in genetic diversity , 2018, Molecular ecology.

[37]  J. Malmkvist,et al.  Effects of thermal environment on hypothalamic-pituitary-adrenal axis hormones, oxytocin, and behavioral activity in periparturient sows. , 2009, Journal of animal science.

[38]  R. Kassen,et al.  Adaptive synonymous mutations in an experimentally evolved Pseudomonas fluorescens population , 2014, Nature Communications.

[39]  K. Holsinger,et al.  Genetic and evolutionary consequences of small population size in plants : implications for conservation , 1991 .

[40]  Gregory P. Brown,et al.  Invader immunology: invasion history alters immune system function in cane toads (Rhinella marina) in tropical Australia. , 2015, Ecology letters.

[41]  R. Tingley,et al.  Hydric balance and locomotor performance of an anuran (Rhinella marina) invading the Australian arid zone , 2012 .

[42]  Gregory P. Brown,et al.  It is lonely at the front: contrasting evolutionary trajectories in male and female invaders , 2016, Royal Society Open Science.

[43]  Gregory P. Brown,et al.  Sexual and geographical divergence in head widths of invasive cane toads, Rhinella marina (Anura: Bufonidae), is driven by both rapid evolution and plasticity , 2018 .

[44]  J. Goudet HIERFSTAT , a package for R to compute and test hierarchical F -statistics , 2005 .

[45]  T. White,et al.  Adaptive evolution during an ongoing range expansion: the invasive bank vole (Myodes glareolus) in Ireland , 2013, Molecular ecology.

[46]  H. MacIsaac,et al.  Is invasion success explained by the enemy release hypothesis , 2004 .

[47]  J. L. Parra,et al.  Very high resolution interpolated climate surfaces for global land areas , 2005 .

[48]  F. Rousset,et al.  Absence of evidence for isolation by distance in an expanding cane toad (Bufo marinus) population: an individual‐based analysis of microsatellite genotypes , 2000, Molecular ecology.

[49]  J. Plotkin,et al.  Synonymous but not the same: the causes and consequences of codon bias , 2011, Nature Reviews Genetics.

[50]  Gregory P. Brown,et al.  Is the behavioural divergence between range-core and range-edge populations of cane toads (Rhinella marina) due to evolutionary change or developmental plasticity? , 2017, Royal Society Open Science.

[51]  D. Heath,et al.  Standing genetic diversity and selection at functional gene loci are associated with differential invasion success in two non‐native fish species , 2018, Molecular ecology.

[52]  O. Gaggiotti,et al.  A Genome-Scan Method to Identify Selected Loci Appropriate for Both Dominant and Codominant Markers: A Bayesian Perspective , 2008, Genetics.

[53]  R. Shine,et al.  A genetic perspective on rapid evolution in cane toads (Rhinella marina) , 2015, Molecular ecology.

[54]  S. Copley,et al.  Synonymous mutations make dramatic contributions to fitness when growth is limited by a weak-link enzyme , 2018, PLoS genetics.

[55]  N. Grünwald,et al.  VcfR : a package to manipulate and visualize VCF data in , 2016 .

[56]  É. Frichot,et al.  LEA: An R package for landscape and ecological association studies , 2015 .

[57]  C. Moritz,et al.  Phylogeography of Bufo marinus from its natural and introduced ranges , 1998, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[58]  Gregory P. Brown,et al.  The thermal dependency of locomotor performance evolves rapidly within an invasive species , 2018, Ecology and evolution.

[59]  P. Legendre,et al.  vegan : Community Ecology Package. R package version 1.8-5 , 2007 .

[60]  Donald Blomqvist,et al.  Trapped in the extinction vortex? Strong genetic effects in a declining vertebrate population , 2010, BMC Evolutionary Biology.

[61]  A. Baer The genetic perspective , 1977 .

[62]  A. Toland,et al.  The Impact of 3′UTR Variants on Differential Expression of Candidate Cancer Susceptibility Genes , 2013, PloS one.

[63]  G. Somero,et al.  Following the heart: temperature and salinity effects on heart rate in native and invasive species of blue mussels (genus Mytilus) , 2006, Journal of Experimental Biology.

[64]  Thomas F. Cross,et al.  diveRsity: An R package for the estimation and exploration of population genetics parameters and their associated errors , 2013 .

[65]  R. Frankham Genetics and extinction , 2005 .

[66]  S. Bonatto,et al.  Multiple introductions and gene flow in subtropical South American populations of the fireweed, Senecio madagascariensis(Asteraceae) , 2016, Genetics and molecular biology.

[67]  Fiona J. Thomson,et al.  High genetic diversity is not essential for successful introduction , 2013, Ecology and evolution.

[68]  Laurent Excoffier,et al.  PGDSpider: an automated data conversion tool for connecting population genetics and genomics programs , 2012, Bioinform..

[69]  Gregory P. Brown,et al.  An evolutionary process that assembles phenotypes through space rather than through time , 2011, Proceedings of the National Academy of Sciences.

[70]  F. Allendorf,et al.  Introduction: Population Biology, Evolution, and Control of Invasive Species , 2003 .

[71]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[72]  Niklaus J Grünwald,et al.  vcfr: a package to manipulate and visualize variant call format data in R , 2017, Molecular ecology resources.

[73]  B. Schmid,et al.  Latitudinal population differentiation in two species of Solidago (Asteraceae) introduced into Europe. , 1998, American journal of botany.

[74]  J. L. Bouzat,et al.  Tracking the long-term decline and recovery of an isolated population , 1998, Science.

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

[76]  C. Goodnight,et al.  GENETIC VARIATION AND EVOLUTION OF SECONDARY COMPOUNDS IN NATIVE AND INTRODUCED POPULATIONS OF THE INVASIVE PLANT MELALEUCA QUINQUENERVIA , 2012, Evolution; international journal of organic evolution.

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

[78]  R. Francis,et al.  pophelper: an R package and web app to analyse and visualize population structure , 2017, Molecular ecology resources.

[79]  J. Cornuet,et al.  Inferring population history from microsatellite and enzyme data in serially introduced cane toads, Bufo marinus. , 2001, Genetics.

[80]  R. Greger Physiology of renal sodium transport. , 2000, The American journal of the medical sciences.

[81]  M. Feder,et al.  Evolution of heat‐shock protein expression underlying adaptive responses to environmental stress , 2018, Molecular ecology.

[82]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.