Assessing models of speciation under different biogeographic scenarios; an empirical study using multi-locus and RNA-seq analyses

Abstract

Evolutionary biology often seeks to decipher the drivers of speciation, and much debate persists over the relative importance of isolation and gene flow in the formation of new species. Genetic studies of closely related species can assess if gene flow was present during speciation, because signatures of past introgression often persist in the genome. We test hypotheses on which mechanisms of speciation drove diversity among three distinct lineages of desert tortoise in the genus Gopherus. These lineages offer a powerful system to study speciation, because different biogeographic patterns (physical vs. ecological segregation) are observed at opposing ends of their distributions. We use 82 samples collected from 38 sites, representing the entire species' distribution and generate sequence data for mtDNA and four nuclear loci. A multilocus phylogenetic analysis in *BEAST estimates the species tree. RNA-seq data yield 20,126 synonymous variants from 7665 contigs from two individuals of each of the three lineages. Analyses of these data using the demographic inference package ∂a∂i serve to test the null hypothesis of no gene flow during divergence. The best-fit demographic model for the three taxa is concordant with the *BEAST species tree, and the ∂a∂i analysis does not indicate gene flow among any of the three lineages during their divergence. These analyses suggest that divergence among the lineages occurred in the absence of gene flow and in this scenario the genetic signature of ecological isolation (parapatric model) cannot be differentiated from geographic isolation (allopatric model).

References April, J., R. H. Hanner, A. M. Dion-Cote, and L. Bernatchez. 2013. Glacial cycles as an allopatric speciation pump in north-eastern American freshwater fishes. Mol. Ecol. 22:409–422. Arnold, M. L. 2007. Evolution through genetic exchange. Oxford Univ. Press, London. Avise, J. C., B. W. Bowen, T. Lamb, A. B. Meylan, and E. Bermingham. 1992. Mitochondrial DNA evolution at a turtle's pace: evidence for low genetic variability and reduced microevolutionary rate in the Testudines. Mol. Biol. Evol. 9:457–473. Barton, N. H., and G. M. Hewitt. 1985. Analysis of hybrid zones. Annu. Rev. Ecol. Syst. 16:113–148. Barton, N. H., and G. M. Hewitt. 1989. Adaptation, speciation, and hybrid zones. Nature 341:497–503. Bernardi, G. 2013. Speciation in fishes. Mol. Ecol. 22:5487–5502. Bolger, A. M., M. Lohse, and B. Usadel. 2014. Trimmomatic: a flexible trimmer for illumina sequence data. BioInformatics 30:2114–2120. doi: 10.1093/bioinformatics/btu170 Bouckaert, R. R., and J. Heled. 2014. DensiTree 2: seeing trees through the forest. bioRxiv, http://dx.doi.org/10.1101/012401 Bouckaert, R. R., J. Heled, D. Kühnert, T. Vaughan, C.-H. Wu, D. Xie, et al. 2014. BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Comput. Biol. 10:e1003537. Bramble, D. M., and J. H. Hutchison. 2014. Morphology, taxonomy, and distribution of North American tortoises; an evolutionary perspective. Pp. 1–12 in D. C. Rostal, E. D. McCoy and H. R. Mushinsky, eds. Biology and conservation of North American tortoises. Johns Hopkins Univ. Press, Baltimore, MD. Brito, P. H., and S. V. Edwards. 2009. Multilocus phylogeography and phylogenetics using sequence-based markers. Genetica 135:439–455. Brown, D. E., and C. H. Lowe. 1980. Biotic communities of the southwest. Rocky Mountain Forest and Range Experiment Station, U.S. Forest Service. General Technical Report RM-78. Map. Búrquez, A., A. Martínez-Yrízar, R. S. Felger, and D. Yetman. 1999. Vegetation and habitat diversity at the southern edge of the Sonoran Desert. Pp. 36–67 in R. H. Robichaux, ed. Ecology of Sonoran desert plants and plant communities. University of Arizona Press, Tucson, AZ. Cahais, V., P. Gayral, G. Tsagkogeorga, J. Melo-Ferreira, M. Ballenghien, L. Weinert, et al. 2012. Reference-free transcriptome assembly in non-model animals from next-generation sequencing data. Mol. Ecol. Resour. 12:834–845. Cooke, G. M., E. L. Landguth, and L. B. Beheregaray. 2014. Riverscape genetics identifies replicated ecological divergence across an Amazonian ecotone. Evolution 68:1947–1960. Coyne, J. A., and H. A. Orr. 2004. Speciation. Sinauer Associates, Inc., Sunderland, MA. Cruickshank, T. E., and M. W. Hahn. 2014. Reanalysis suggests that genomic islands of speciation are due to reduced diversity, not reduced gene flow. Mol. Ecol. 23:3133–3157. De Wit, P., M. H. Pespeni, J. T. Ladner, D. J. Barshis, F. Seneca, H. Jaris, et al. 2012. The simple fool's guide to population genomics via RNA-Seq: an introduction to high-throughput sequencing data analysis. Mol. Ecol. Resour. 12:1058–1067. Degnan, J. H., and N. A. Rosenberg. 2009. Gene tree discordance, phylogenetic inference and the multispecies coalescent. Trends Ecol. Evol. 24:332–340. DePristo, M., E. Banks, R. Poplin, K. Garimella, J. Maguire, C. Hartl, et al. 2011. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43:491–498. Doran, A. G., and C. J. Creevey. 2013. Snpdat: easy and rapid annotation of results from de novo snp discovery projects for model and non-model organisms. BMC Bioinformatics 14:45. Drummond, A. J., S. Y. W. Ho, M. J. Phillips, and A. Rambaut. 2006. PLoS Biol. 4:e88. Edwards, T. 2003. Desert tortoise conservation genetics, Unpublished M.S. Thesis, University of Arizona, Tucson, USA. Edwards, T., and K. H. Berry. 2013. Are captive tortoises a reservoir for conservation? An assessment of genealogical affiliation of captive Gopherus agassizii to local, wild populations. Conserv. Genet. 14:649–659. Edwards, T., and S. Harrison. 2014. Population and conservation genetics. Pp. 127–133 in D. C. Rostal, E. D. McCoy and H. R. Mushinsky, eds. Biology and conservation of North American tortoises. Johns Hopkins Univ. Press, Baltimore, MD. Edwards, T., C. R. Schwalbe, D. E. Swann, and C. S. Goldberg. 2004. Implications of anthropogenic landscape change on inter-population movements of the desert tortoise (Gopherus agassizii). Conserv. Genet. 5:485–499. Edwards, T., C. J. Jarchow, C. A. Jones, and K. E. Bonine. 2010. Tracing genetic lineages of captive desert tortoises in Arizona. J. Wildl. Manage. 74:801–807. Edwards, T., M. Vaughn, Ma. C. Meléndez-Torres, A. E. Karl, P. R. Rosen, K. H. Berry, et al. 2012. A biogeographic perspective of speciation among desert tortoises in the genus Gopherus: a preliminary evaluation. pp. 243–247 in G. J. Gottfried, P. F. Ffolliott, B. S. Gebow, L. G. Eskew, eds. Merging science and management in a rapidly changing world: biodiversity and management of the Madrean Archipelago III, 2012 May 1-5, Tucson, AZ. Proceedings RMRS-P-67. U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station, Fort Collins, CO. Edwards, T., E. C. Cox, V. Buzzard, C. Wiese, L. S. Hillard, and R. W. Murphy. 2014. Genetic assessments and parentage analysis of captive Bolson tortoises (Gopherus flavomarginatus) inform their “rewilding” in New Mexico. PLoS ONE 9:e102787. Edwards, T., K. H. Berry, R. D. Inman, T. C. Esque, K. E. Nussear, C. A. Jones, et al. 2015a. Testing taxon tenacity of tortoises: evidence for a geographical-selection gradient at a secondary contact zone. Ecol. Evol. 5:2095–2114. Edwards, T., M. Vaughn, P. C. Rosen, C. Meléndez-Torres, A. E. Karl, M. Culver, et al. 2015b. Shaping species with ephemeral boundaries; the distribution and genetic structure of the desert tortoise (Gopherus morafkai) in the Sonoran Desert. J. Biogeogr., doi: 10.1111/jbi.12664. Endler, J. A. 1977. Geographic variation, speciation, and clines. Princeton Univ. Press, Princeton, NJ. Feder, J. L., and P. Nosil. 2010. The efficacy of divergence hitchhiking in generating genomic islands during ecological speciation. Evolution 64:1729–1747. Felger, R. S., D. F. Austin, T. R. Van Devender, J. Sánchez-Escalante, and M. Costea. 2012. Convolvulaceae of Sonora, Mexico – I. Convolvulus, Cressa, Dichondra, Evolvulus, Ipomoea, Jacquemontia, Merremia, Operculina. J. Bot. Res. Inst. Tex. 6:459–527. Fisher-Reid, M. C., T. N. Engstrom, C. A. Kuczynski, P. R. Stephens, and J. J. Wiens. 2013. Parapatric divergence of sympatric morphs in a salamander: incipient speciation on Long Island? Mol. Ecol. 22:4681–4694. Fujita, M. K., T. N. Engstrom, D. E. Starkey, and H. B. Shaffer. 2004. Turtle phylogeny: insights from a novel nuclear intron. Mol. Phylogenet. Evol. 31:1031–1040. Geraldes, A., N. Ferrand, and M. W. Nachman. 2006. Contrasting patterns of introgression at X-linked loci across the hybrid zone between subspecies of the European rabbit (Oryctolagus cuniculus). Genetics 173:919–933. Germano, D. J., R. B. Bury, T. C. Esque, T. H. Fritts, and P. A. Medica. 1994. Range and habitats of the desert tortoise. Pp. 73–84 in R. B. Bury, D. J. Germano, eds. Biology of North American tortoises, vol. 13. U.S. Department of Interior National Biological Survey, Fish and Wildlife Research, Washington, DC. Grabherr, M. G., B. J. Haas, M. Yassour, J. Z. Levin, D. A. Thompson, I. Amit, et al. 2011. Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat. Biotechnol. 29:644–652. Gutenkunst, R. N., R. D. Hernandez, S. H. Williams, and C. D. Bustamante. 2009. Inferring the joint demographic history of multiple populations from multidimensional SNP frequency data. PLoS Genet. 5:e1000695. Haas, B. J., A. Papanicolaou, M. Yassour, M. Grabherr, P. D. Blood, J. Bowden, et al. 2013. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 8:1494–1512. Hagerty, B. E., K. E. Nussear, T. C. Esque, and C. R. Tracy. 2011. Making molehills out of mountain landscape genetics of the Mojave desert tortoise. Landscape Ecol. 2:267–280. Hearn, J., G. N. Stone, L. Bunnefeld, J. A. Nicholls, N. H. Barton, and K. Lohse. 2014. Likelihood-based inference of population history from low-coverage de novo genome assemblies. Mol. Ecol. 23:198–211. Heled, J., and A. J. Drummond. 2008. Bayesian inference of population size history from multiple loci. BMC Evol. Biol. 8:289. Heled, J., and A. J. Drummond. 2010. Bayesian inference of species trees from multilocus data. Mol. Biol. Evol. 27:570–580. Hewitt, G. M. 1996. Some genetic consequences of ice ages, and their role, in divergence and speciation. Biol. J. Linn. Soc. 58:247–276. Hudson, R. R., and M. Turelli. 2003. Stochasticity overrules the “three-times rule”: genetic drift, genetic draft, and coalescence times for nuclear loci versus mitochondrial DNA. Evolution 57:182–190. Jones, F. C., M. G. Grabherr, Y. F. Chan, P. Russel, E. Mauceli, J. Johnson, et al. 2012. The genomic basis of adaptive evolution in threespine sticklebacks. Nature 484:55–61. Kisel, Y., and T. G. Barraclough. 2010. Speciation has a spatial scale that depends on levels of gene flow. Am. Nat. 175:316–334. Kulathinal, R. J., L. S. Stevison, and M. A. F. Noor. 2009. The genomics of speciation in Drosophila: diversity, divergence, and introgression estimated using low-coverage genome sequencing. PLoS Genet. 5:e1000550. Lamb, T., and C. Lydeard. 1994. A molecular phylogeny of the gopher tortoises, with comments on familial relationships within the Testudinoidea. Mol. Phylogenet. Evol. 3:283–291. Lamb, T., J. C. Avise, and J. W. Gibbons. 1989. Phylogeographic patterns in mitochondrial DNA of the desert tortoise (Xerobates agassizii), and evolutionary relationships among the North American gopher tortoises. Evolution 43:76–87. Leaché, A. D., and J. A. McGuire. 2006. Phylogenetic relationships of horned lizards (Phrynosoma) based on nuclear and mitochondrial data: evidence for a misleading mitochondrial gene tree. Mol. Phylogenet. Evol. 39:628–644. Leaché, A. D., and B. Rannala. 2011. The accuracy of species tree estimation under simulation: a comparison of methods. Syst. Biol. 60:126–137. Li, H., and R. Durbin. 2009. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25:1754–1760. Li, H., B. Handsaker, A. Wysoker, T. Fennell, J. Ruan, N. Homer, et al. (2009) The sequence alignment/map (SAM) format and SAMtools. Bioinformatics, 25, 2078–2079. Librado, P., and J. Rozas. 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25:1451–1452. Lohse, K., R. J. Harrison, and N. H. Barton. 2011. A general method for calculating likelihoods under the coalescent process. Genetics 58:977–987. Magrane, M., and UniProt Consortium. 2011. UniProt knowledgebase: a hub of integrated protein data. Database: J Biol Databases Curation. 2011:bar009. doi: 10.1093/database/bar009. McCormack, J. E., H. Huang, and L. L. Knowles. 2009. Maximum likelihood estimates of species trees: how accuracy of phylogenetic inference depends upon the divergence history and sampling design. Syst. Biol. 58:501–508. McGettigan, P. A. 2013. Transcriptomics in the RNA-seq era. Curr. Opin. Chem. Biol. 17:4–11. McKenna, A., M. Hanna, E. Banks, A. Sivachenko, K. Cibulskis, A. Kernytsky, et al. 2010. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20:1297–1303. McLuckie, A. M., T. Lamb, C. R. Schwalbe, and R. D. McCord. 1999. Genetic and morphometric assessment of an unusual tortoise (Gopherus agassizii) population in the Black Mountains of Arizona. J. Herpetol. 33:36–44. Meitern, R., R. Andreson, and P. Horak. 2014. Profile of whole blood gene expression following immune stimulation in a wild passerine. BMC Genom. 15:533. Melo-Ferreira, J., P. C. Alves, H. Freitas, N. Ferrand, and P. Boursot. 2009. The genomic legacy from the extinct Lepus timidus to the three hare species of Iberia: contrast between mtDNA, sex chromosomes and autosomes. Mol. Ecol. 18:2643–2658. Mendez, F. L., J. C. Watkins, and M. F. Hammer. 2012. A haplotype at STAT2 introgressed from Neanderthals and serves as a candidate of positive selection in Papua New Guinea. Am. J. Hum. Genet. 91:265–274. Morafka, D. J., G. Aguirre, and R. W. Murphy. 1994. Allozyme differentiation among gopher tortoises (Gopherus) - conservation genetics and phylogenetic and taxonomic implications. Can. J. Zool. 72:1665–1671. Murphy, R. W., K. H. Berry, T. Edwards, and A. M. McLuckie. 2007. A genetic assessment of the recovery units for the Mojave population of desert tortoises, Gopherus agassizii. Chelonian Conserv. Biol. 6:229–251. Murphy, R. W., K. H. Berry, T. Edwards, A. E. Leviton, A. Lathrop, and J. D. Riedle. 2011. The dazed and confused identity of Agassiz's land tortoise, Gopherus agassizii (Testudines: Testudinidae) with the description of a new species, and its consequences for conservation. ZooKeys 113:39–71. Nachman, M. W., and B. A. Payseur. 2012. Recombination rate variation and speciation: theoretical predictions and empirical results from rabbits and mice. Philos. Trans. R. Soc. B Biol. Sci. 367:409–421. Nylander, J. A. A. 2004. MrModeltest v2. Program distributed by the author. Evolutionary Biology Centre, Uppsala University, Uppsala, Sweden. Patterson, N., A. Price, and D. Reich. 2006. Population structure and eigenanalysis. PLoS Genet. 2:e190. Payseur, B. A. 2010. Using differential introgression in hybrid zones to identify genomic regions involved in speciation. Mol. Ecol. 10:806–820. Payseur, B. A., J. G. Krenz, and M. W. Nachman. 2004. Differential patterns of introgression across the X chromosome in a hybrid zone between two species of house mice. Evolution 58:2064–2078. Pinho, C., and J. Hey. 2010. Divergence with gene flow: models and data. Annu. Rev. Ecol. Evol. Syst. 41:215–230. Pollard, D. A., V. N. Iyer, A. M. Moses, and M. B. Eisen. 2006. Widespread discordance of gene trees with species tree in Drosophila: evidence for incomplete lineage sorting. PLoS Genet. 2:e173. Rambaut, A. 2006–2012. FigTree; tree figure drawing tool. Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, United Kingdom. Rambaut, A., and A. J. Drummond. 2002–2012. TreeAnnotator; MCMC output analysis. Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, United Kingdom. Rambaut, A., M. Suchard, W. Xie, and A. J. Drummond. 2003–2013. Tracer: MCMC trace analysis tool. Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, United Kingdom. Riddle, B. R., and D. J. Hafner. 2006. A step-wise approach to integrating phylogeographic and phylogenetic biogeographic perspectives on the history of a core North American warm deserts biota. J. Arid Environ. 66:435–461. Robinson, J. D., A. J. Coffman, M. J. Hickerson, and R. N. Gutenkunst. 2014. Sampling strategies for frequency spectrum-based population genomic inference. BMC Evol. Biol. 14:254. Rundle, H. D., and P. Nosil. 2005. Ecological speciation. Ecol. Lett. 8:336–352. Smadja, C. M., and R. K. Butlin. 2011. A framework for comparing processes of speciation in the presence of gene flow. Mol. Ecol. 20:5123–5140. Spinks, P. Q., H. B. Shaffer, J. B. Iverson, and W. P. McCord. 2004. Phylogenetic hypotheses for the turtle family Geoemydidae. Mol. Phylogenet. Evol. 32:164–182. Stephens, M., and P. Donnelly. 2003. A comparison of Bayesian methods for haplotype reconstruction from population genotype data. Am. J. Hum. Genet. 73:1162–1169. Swofford, D. L. 2002. PAUP* Phylogenetic analysis using parsimony (*and other methods). Version 4. Sinauer Associates, Sunderland, MA. Teeter, K. C., B. A. Payseur, L. W. Harris, M. A. Bakewell, L. M. Thibodeau, J. E. O'Brien, et al. 2008. Genome-wide patterns of gene flow across a house mouse hybrid zone. Genome Res. 18:67–76. Thomson, R. C., A. M. Shedlock, S. V. Edwards, and H. B. Shaffer. 2008. Developing markers for multilocus phylogenetics in non-model organisms: a test case with turtles. Mol. Phylogenet. Evol. 49:514–525. Trapnell, C., B. A. Williams, G. Pertea, A. Mortazavi, G. Kwan, M. J. van Baren, et al. 2010. Transcript assembly and quantification by RNA-seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28:511–515. U.S. Fish and Wildlife Service. 1994. Desert tortoise (Mojave population) recovery plan. U.S. Fish and Wildlife Service, Portland, OR. U.S. Fish, Wildlife Service. 2011. Revised recovery plan for the Mojave population of the desert tortoise (Gopherus agassizii). U.S. Fish and Wildlife Service, Pacific Southwest Region, Sacramento, CA. Van Devender, T. R. 2000. The deep history of the Sonoran Desert. Pp. 61–69 in S. J. Phillips and P. W. Comus, eds. A natural history of the Sonoran desert. University of California Press, Berkeley, CA. Via, S. 2009. Natural selection in action during speciation. Proc. Natl Acad. Sci. USA 106:9939–9946. Wang, Y., and J. Hey. 2010. Estimating divergence parameters with small samples from a large number of loci. Genetics 184:363–379. Williamson, S. H., R. Hernandez, A. Fledel-Alon, L. Zhu, R. Nielsen, and C. D. Bustamante. 2005. Simultaneous inference of selection and population growth from patterns of variation in the human genome. Proc. Natl Acad. Sci. USA 102:7882–7887. Yang, Z. 2007. PAML 4: a program package for phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24:1586–1591. Zhang, L. X. 2011. From gene trees to species trees II: species tree inference by minimizing deep coalescence events. IEEE/ACM Trans. Comput. Biol. Bioinform. 8:1685–1691.