Integrative species delimitation uncovers hidden diversity within the Pithecopus hypochondrialis species complex (Hylidae, Phyllomedusinae) and its phylogeography reveals Plio-Pleistocene connectivity among Neotropical savannas

Abreu-Jardim, 2023, Neogene-Quaternary tectonic, eustatic and climatic events shaped the evolution of a South American treefrog, J. Biogeogr., 50, 987, 10.1111/jbi.14578 Akinwande, 2015, Variance Inflation Factor: As a Condition for the Inclusion of Suppressor Variable(s) in Regression Analysis, Open J. Stat., 5, 754, 10.4236/ojs.2015.57075 Aktas, 2020, haplotypes: manipulating DNA sequences and estimating unambiguous haplotype network with statistical parsimony, R Package Version, 1, 2 Altschul, 1990, Basic local alignment search tool, J. Mol. Biol., 215, 403, 10.1016/S0022-2836(05)80360-2 Andrade, 2020, A new cryptic species of Pithecopus (Anura, Phyllomedusidae) in north-eastern Brazil, Eur. J. Taxon., 723, 108, 10.5852/ejt.2020.723.1147 Avila-Pires, 1995, Lizards of Brazilian Amazonia (Reptilia: Squamata), Zool. Verh. Leiden, 299, 1 Azevedo, 2020, On the young savannas in the land of ancient forest, 271 Baêta, 2009, Phyllomedusa itacolomi Caramaschi, Cruz & Feio, 2006, as junior synonym of Phyllomedusa ayeaye (B. Lutz, 1966) (Hylidae, Phyllomedusinae), Zootaxa, 2226, 58, 10.11646/zootaxa.2226.1.5 Bandeira, 2021, Different elevational environments dictate contrasting patterns of niche evolution in Neotropical Pithecopus treefrog species, Biotropica, 53, 1042, 10.1111/btp.12929 Barrio-Amorós, 2009, Distribución y aspectos de la historia natural de las ranas lémur (Hylidae: Phyllomedusinae) en Venezuela, Mem. Fund. La Salle Cienc. Nat., 171, 19 Bergsten, 2012, The effect of geographical scale of sampling on DNA barcoding, Syst. Biol., 61, 851, 10.1093/sysbio/sys037 Bernarde, 2007, Calling period and environments of anuran fauna in municipality of Espigão do Oeste, Rondônia, Southeastern of Amazonian – Brazil (Amphibia: Anura), Biota Neotrop., 7, 88 Bielejec, 2011, SPREAD: spatial phylogenetic reconstruction of evolutionary dynamics, Bioinformatics, 27, 2910, 10.1093/bioinformatics/btr481 Bouckaert, 2017, bModelTest: Bayesian phylogenetic site model averaging and model comparison, BMC Evol. Biol., 17, 42, 10.1186/s12862-017-0890-6 Bouckaert, 2019, BEAST 2.5: an advanced software platform for Bayesian evolutionary analysis, PLoS Comput. Biol., 15, e1006650, 10.1371/journal.pcbi.1006650 Brandão, 2002, A new species of Phyllomedusa Wagler, 1830 (Anura: Hylidae) from Central Brazil, J. Herp., 36, 571, 10.1670/0022-1511(2002)036[0571:ANSOPW]2.0.CO;2 Brown, J.L., Hill, D.J., Dolan, A.M., Carnaval, A.C., Haywood, A.M., 2018. Data Descriptor: PaleoClim, high spatial resolution paleoclimate surfaces for global land areas. Sci. Data 5, 180254 (2018). 10.1038/sdata.2018.254. Brunes, 2014, Species limits, phylogeographic and hybridization patterns in Neotropical leaf frogs (Phyllomedusinae), Zool. Scr., 43, 586, 10.1111/zsc.12079 Bruschi, 2012, Cytogenetic analysis of two species in the Phyllomedusa hypochondrialis group, Hereditas, 149, 34, 10.1111/j.1601-5223.2010.02236.x Bruschi, 2013, Evaluation of the taxonomic status of populations assigned to Phyllomedusa hypochondrialis (Anura, Hylidae, Phyllomedusinae) based on molecular, chromosomal, and morphological approach, BMC Genet., 14, 70, 10.1186/1471-2156-14-70 Bueno, 2017, Effects of Quanternary climatic fluctuations on the distribution of Neotropical savanna tree species, Ecography, 39, 1 Buzatti, 2017, Multiple Pleistocene refugia in the Brazilian cerrado: evidence from phylogeography and climatic niche modelling of two Qualea species (Vochysiaceae), Bot. J. Linn. Soc., 185, 307, 10.1093/botlinnean/box062 Buzatti, 2018, Genetic and historical colonization analyses of an endemic savanna tree, Qualea grandiflora, reveal ancient connections between amazonian savannas and Cerrado core, Front. Plant Sci., 9, 981, 10.3389/fpls.2018.00981 Caramaschi, 2006, Redefinição do grupo de Phyllomedusa hypochondrialis, com redescrição de P. megacephala (Miranda-Ribeiro, 1926), revalidação de P. azurea Cope, 1862 e descrição de uma nova espécie (Amphibia, Anura, Hylidae), Arq. Mus. Nac. (rio De J.), 64, 159 Carnaval, 2008, Historical climate modelling predicts patterns of current biodiversity in the Brazilian Atlantic forest, J. Biogeogr., 35, 1187, 10.1111/j.1365-2699.2007.01870.x Carstens, 2013, How to fail at species delimitation, Mol. Ecol., 22, 4369, 10.1111/mec.12413 Chambers, 2020, The multispecies coalescent over-splits species in the case of geographically widespread taxa, Syst. Biol., 69, 184, 10.1093/sysbio/syz042 Cheng, 2013, Climate change patterns in Amazonia and biodiversity, Nat. Commun., 4, 1411, 10.1038/ncomms2415 Cobos, 2019, kuenm: an R package for detailed development of ecological niche models using Maxent, PeerJ, 7, e6281, 10.7717/peerj.6281 Crawford, 2003, Huge populations and old species of Costa Rican and Panamanian dirt frogs inferred from mitochondrial and nuclear gene sequencies, Mol. Ecol., 12, 2525, 10.1046/j.1365-294X.2003.01910.x Dolan, 2015, Modelling the enigmatic Late Pliocene Glacial Event – Marine Isotope Stage M2, Glob. Planet. Change, 128, 47, 10.1016/j.gloplacha.2015.02.001 Dowsett, 2016, The PRISM4 (mid-Piacenzian) paleoenvironmental reconstruction, Clim. past., 12, 1519, 10.5194/cp-12-1519-2016 Drummond, 2006, Relaxed phylogenetics and dating with confidence, PLoS Biol., 4, e88, 10.1371/journal.pbio.0040088 Drummond, 2012, Bayesian phylogenetics with BEAUti and the BEAST 1.7, Mol. Biol. Evol., 29, 1969, 10.1093/molbev/mss075 Duellman, 1994 Duellman, 2016, Phylogenetics, classification, and biogeography of the treefrogs (Amphibia: Anura: Arboranae), Zootaxa, 4104, 1, 10.11646/zootaxa.4104.1.1 Duellman, W.E., 1970. The Hylid Frogs of Middle America, first ed. University of Kansas, Lawrence, Kansas. 10.5962/bhl.title.2835. Elith, 2011, A statistical explanation of MaxEnt for ecologists, Divers. Distrib., 17, 43, 10.1111/j.1472-4642.2010.00725.x Excoffier, 2010, Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows, Mol. Ecol. Res., 10, 564, 10.1111/j.1755-0998.2010.02847.x Faivovich, J., Haddad, C.F.B., Baêta, D., Jungfer., Álvares, G.F.R, Brandão, R.A., Sheil, C., Barrientos, L.S., Barrio-Amorós, C.L., Cruz, C.A.G., Wheeler, W.C., 2010. The phylogenetic relationships of the charismatic poster frogs, Phyllomedusinae (Anura, Hylidae). Cladistics 26, 227–261. 10.1111/j.1096-0031.2009.00287.x. Flot, 2010, SeqPHASE: a web tool for interconverting phase input/output files and fsata sequence alignments, Mol. Ecol. Resour., 10, 162, 10.1111/j.1755-0998.2009.02732.x Fonte, 2019, Long-Distance Dispersal in Amphibians, Front. Biogeogr., 11, e44577, 10.21425/F5FBG44577 Fordham, 2017, PaleoView: a tool for generating continuous climatic projections spanning the last 21 000 years at regional and global scales, Ecography, 40, 1348, 10.1111/ecog.03031 Fox, 2018 Frost, D.R., 2023. Amphibian Species of the World: an online reference. Version 6.1. Am. Mus. Nat. Hist. https://amphibiansoftheworld.amnh.org/index.php (accessed 9 March 2023). 10.5531/db.vz.0001. Gehara, 2014, High levels of diversity uncovered in a widespread nominal taxon: continental phylogeography of the Neotropical tree frog Dendropsophus minutus, PLoS ONE, 9, e103958, 10.1371/journal.pone.0103958 Gernhard, 2008, The conditioned reconstructed process, J. Theor. Biol., 253, 769, 10.1016/j.jtbi.2008.04.005 Giaretta, 2007, A new Phyllomedusa Wagler (Anura, Hylidae) with reticulated pattern on flanks from Southeastern Brazil, Zootaxa, 1614, 31, 10.11646/zootaxa.1614.1.3 Guarnizo, 2013, Genetic divergence within frog species is greater in topographically more complex regions, J. Zoolog. Syst. Evol. Res., 51, 333, 10.1111/jzs.12027 Guarnizo, 2015, DNA barcoding survey of anurans across the eastern cordillera of the Andes on cryptic diversity, PLoS ONE, 10, e0127312, 10.1371/journal.pone.0127312 Guillot, 2005, GENELAND: a computer package for landscape genetics, Mol. Ecol. Notes, 5, 712, 10.1111/j.1471-8286.2005.01031.x Guillot, 2012, A unifying model for the analysis of phenotypic, genetic, and geographic data, Syst. Biol., 61, 897, 10.1093/sysbio/sys038 Haga, 2017, Unrevealing the leaf frogs Cerrado diversity: a new species of Pithecopus (Anura, Arboranae, Phyllomedusidae) from Mato Grosso state, Brazil. Plos ONE, 12, e0184631, 10.1371/journal.pone.0184631 Heyer, 1990, Frogs of Boraceia. Arq. Zool. (são Paulo), 31, 231 Hill, 2015, The non-analogue nature of Pliocene temperature gradients, Earth Planet. Sci. Lett., 425, 232, 10.1016/j.epsl.2015.05.044 IBGE (Instituto Brasileiro de Geografia e Estatística), 2019. Biomas e sistema costeiro-marinho do Brasil: compatível com a escala 1:250000. IBGE, Coordenação de Recursos Naturais e Estudos Ambientais, Rio de Janeiro, Brazil. Jackson, 2017, Species delimitation with gene flow, Syst. Biol., 66, 799 Jones, 2017, Algorithmic improvements to species delimitation and phylogeny estimation under the multispecies coalescent, J. Math. Biol., 74, 447, 10.1007/s00285-016-1034-0 Kalkauskas, 2021, Sampling bias and model choice in continuous phylogeography: getting lost on random walk, PLoS Comput. Biol., 17, e1008561, 10.1371/journal.pcbi.1008561 Katoh, 2019, MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization, Brief. Bioinform., 20, 1160, 10.1093/bib/bbx108 Katoh, 2013, MAFFT multiple sequence alignment software version 7: improvements in performance and usability, Mol. Biol. Evol., 30, 772, 10.1093/molbev/mst010 Katoh, 2008, Recent developments in the MAFFT multiple sequence alignment program, Brief. Bioinform., 9, 286, 10.1093/bib/bbn013 Kern, 2023, A 1.8 million year history of Amazon vegetation, Quat. Sci. Rev., 299, 10.1016/j.quascirev.2022.107867 Kumar, 2018, MEGA X: molecular evolutionary genetics analysis across computing platforms, Mol. Biol. Evol., 35, 1547, 10.1093/molbev/msy096 Langone, 2016, High genetic diversity but low population structure in the frog Pseudopaludicola falcipes (Hensel, 1867) (Amphibia, Anura) from the Pampas of South America, Mol. Phylogenet. Evol., 95, 137, 10.1016/j.ympev.2015.11.012 Leaché, 2019, The spectre of too many species, Syst. Biol., 68, 168, 10.1093/sysbio/syy051 Leite, 2015, Evolutionary affinities of the ‘Lost World’ mouse suggest a late Pliocene connection between the Guiana and Brazilian shields, J. Biogeogr., 42, 706, 10.1111/jbi.12461 Lemey, 2010, Phylogeography takes a relaxed random walk in continuous space and time, Mol. Biol. Evol., 27, 1877, 10.1093/molbev/msq067 Lescure, 1995, Contribution à l’étude des amphibiens de Guyane française. X. Les Phyllomedusa (Anura, Hylidae), Revue Fr. Aquariol., 22, 35 Lima, 2017, Amphibians found in the Amazonian Savanna of the Rio Curiaú environmental protection area in Amapá, Brazil. Biota Neotrop., 17, e20160252 Lima-Rezende, 2019, The spatial genetic structure of the white-banded tanager (Aves, Passeriformes) in fragmented Neotropical savannas suggests two evolutionarily significant units, Biotropica, 51, 234, 10.1111/btp.12623 Liu, 2013, Selecting thresholds for the prediction of species occurrence with presence-only data, J. Biogeogr., 40, 778, 10.1111/jbi.12058 Lleonart, 2000, Removing allometric effects of body size in morphological analysis, J. Theor. Biol., 205, 85, 10.1006/jtbi.2000.2043 Lohse, K., 2009. Can mtDNA barcodes be used to delimit species? A response to Pons et al. (2006). Syst. Biol. 58, 439–442. 10.1093/sysbio/syp039. Magalhães, 2018, Integrative taxonomy helps to assess the extinction risk of anuran species, J. Nat. Conserv., 45, 1, 10.1016/j.jnc.2018.07.001 Magalhães, 2021, Evidence of introgression in endemic frogs from the campo rupestre contradicts the reduced hybridization hypothesis, Biol. J. Linn. Soc., 133, 561, 10.1093/biolinnean/blaa142 Mantel, 1967, The detection of disease clustering and a generalized regression approach, Cancer Res., 27, 209 Martins, 2021, Whiptail lizard lineage delimitation and population expansion as windows into the history of Amazonian open ecosystems, Syst. Biod., 19, 957, 10.1080/14772000.2021.1953185 Mason, 2020, Coalescent-based species delimitation is sensitive to geographic sampling and isolation by distance, Syst. Biod., 18, 269, 10.1080/14772000.2020.1730475 Nascimento, 2019, Extensive cryptic diversity within the Physalaemus cuvieri-Physalaemus ephippifer species complex (Amphibia, Anura) revealed by cytogenetic, mitochondrial, and genomic markers, Front. Genet., 10, 719, 10.3389/fgene.2019.00719 Nylinder, 2014, On the biogeography of Centipeda: a species-tree diffusion approach, Syst. Biol., 63, 178, 10.1093/sysbio/syt102 Ogilvie, H.A., Bouckaert, R.R., Drummond, A.J., 2017. StarBEAST2 brings faster species tree inference and accurate estimates of substitution rates. Mol. Biol. Evol. 34, 2101–2114. 10.1093/molbev/msx126. Olson, 2001, Terrestrial ecoregions of the world: a new map of life on Earth, Bioscience, 51, 933, 10.1641/0006-3568(2001)051[0933:TEOTWA]2.0.CO;2 Oswald, 2022, Colonization rather than fragmentation explains the geographical distribution and diversification of treefrogs endemic to Brazilian shield sky islands, J. Biogeogr., 49, 682, 10.1111/jbi.14320 Otto-Bliesner, B.L., Marshall, S.J., Overpeck, J.T., Miller, G.H., Hu, A., CAPE Last Interglacial Project Members, 2006. Simulating artic climate warmth and icefield retreat in the last interglaciation. Science 311, 1751–1753. 10.1126/science.1120808. Phillips, 2006, Maximum entropy modeling of species geographic distributions, Ecol. Modell., 190, 231, 10.1016/j.ecolmodel.2005.03.026 Pons, 2006, Sequence-based species delimitation for the DNA taxonomy of undescribed insects, Syst. Biol., 55, 595, 10.1080/10635150600852011 Posada, 2000, GeoDis: a program for the cladistic nested analysis of the geographical distribution of genetic haplotypes, Mol. Ecol., 9, 487, 10.1046/j.1365-294x.2000.00887.x Prado, 2012, Cryptic lineages and Pleistocene population expansion in a Brazilian Cerrado frog, Mol. Ecol., 21, 921, 10.1111/j.1365-294X.2011.05409.x Rambaut, 2018, Posterior summarization in Bayesian phylogenetics using Tracer 1.7, Syst. Biol., 67, 901, 10.1093/sysbio/syy032 Ramos, 2019, Cryptic diversity in Brazilian endemic monkey frogs (Hylidae, Phyllomedusinae, Pithecopus) revealed by multispecies coalescent and integrative approaches, Mol. Phylogenet. Evol., 132, 105, 10.1016/j.ympev.2018.11.022 Reid, 2012, Phylogenetic estimation error can decrease the accuracy of species delimitation: a Bayesian implementation of the general mixed Yule-coalescent model, BMC Evol. Biol., 12, 196, 10.1186/1471-2148-12-196 Rodríguez, 2015, Genetic divergence in tropical anurans: deeper phylogeographic structure in forest specialists and in topographically complex regions, Evol. Ecol., 29, 765, 10.1007/s10682-015-9774-7 Roelants, 2007, Global patterns of diversification in the history of modern amphibians, PNAS, 104, 887, 10.1073/pnas.0608378104 Rosenberg, 2004, DISTRUCT: a program for the graphical display of population structure, Mol. Ecol. Notes, 4, 137, 10.1046/j.1471-8286.2003.00566.x Salzburger, 2011, The performance of phylogenetic algorithms in estimating haplotype genealogies with migration, Mol. Ecol., 20, 1952, 10.1111/j.1365-294X.2011.05066.x Santoro, 2014, Reproductive modes, habitat use, and richness of anurans from Chapada dos Veadeiros, central Brazil, North-West. J. Zool., 10, 365 Santos, 2018, Description of the tadpole of Phyllomedusa azurea from the Brazilian Cerrado, with a description of the internal oral morphology of Phyllomedusa oreades, Herpetologica, 74, 50, 10.1655/HERPETOLOGICA-D-16-00050.1 Schmid, 2014, Chromosome banding in Amphibia. XXXI. The Neotropical anuran families Centrolenidae and Allophrynidae, Cytogenet. Genome Res., 142, 268, 10.1159/000362216 Silva, 2002, Biogeographic patterns and conservation in the South American Cerrado: a tropical savanna hotspot, BioScience, 52, 225, 10.1641/0006-3568(2002)052[0225:BPACIT]2.0.CO;2 Simpson, 1951, The species concept, Evolution, 5, 285, 10.2307/2405675 Smith, 2005, Dispersal and the metapopulation paradigm in amphibian ecology and conservation: are all amphibian populations metapopulations?, Ecography, 28, 110, 10.1111/j.0906-7590.2005.04042.x Solís-Lemus, 2015, Bayesian species delimitation combining multiple genes and traits in a unified framework, Evolution, 69, 492, 10.1111/evo.12582 Stephens, 2001, A new statistical method for haplotype reconstruction from population data, Am. J. Hum. Genet., 68, 978, 10.1086/319501 Talavera, 2013, Factors affecting species delimitations with the GMYC model: insights from a butterfly survey, Methods Ecol. Evol., 4, 1101, 10.1111/2041-210X.12107 Thomassen, 2020 Vacher, 2020, Large-scale DNA-based survey of frogs in Amazonia suggests a vast underestimation of species richness and endemism, J. Biogeogr., 47, 1781, 10.1111/jbi.13847 Van Bocxlaer, 2009, Toad radiation reveals into-India dispersal as a source of endemism in the Western Ghats-Sri Lanka biodiversity hotspot, BMC Evol. Biol., 9, 131, 10.1186/1471-2148-9-131 Vasconcellos, 2019, Isolation by instability: historical climate change shapes population structure and genomic divergence of treefrogs in the Neotropical Cerrado savanna, Mol. Ecol., 28, 1748, 10.1111/mec.15045 Vaz-Silva, W., Maciel, N.M., Nomura, F., Morais, A.R., Batista, V.G., Santos, D.L., Andrade, S.P., Oliveira, A.Â.B., Brandão, R.A., Bastos, R.P., 2020. Guia de identificação das espécies de anfíbios (Anura e Gymnophiona) do estado de Goiás e Distrito Federal, Brasil Central, first ed. Sociedade Brasileira de Zoologia, Curitiba, Paraná. 10.7476/9786587590011. Webb, 1991, Ecography and the Great American Interchange, Paleobiology, 17, 266, 10.1017/S0094837300010605 Wells, K. D., 2007. The ecology and behavior of amphibians, first ed. The University of Chicago Press, Chicago, Illinois. http://dx.doi.org/10.7208/chicago/9780226893334.001.0001. Werneck, 2012, Climatic stability in the Brazilian Cerrado: implications for biogeographical connections of South American savannas, species richness and conservation in a biodiversity hotspot, J. Biogeogr., 39, 1695, 10.1111/j.1365-2699.2012.02715.x Wynn, 2001, Do geographically widespread species of tropical amphibians exist? An estimate of genetic relatedness within the neotropical frog Leptodactylus fuscus (Schneider 1799) (Anura Leptodactylidae), Trop. Zool., 14, 255, 10.1080/03946975.2001.10531157 Yeates, 2011, Integrative taxonomy, or iterative taxonomy?, Syst. Ent., 36, 209, 10.1111/j.1365-3113.2010.00558.x