Reconstructing evolution at the community level: A case study on Mediterranean amphibians

Spatial genetic structure in European marbled newts revealed with target enrichment by sequence capture

European marbled newts come in two species that have abutting ranges. The northern species, Triturus marmoratus, is found in France and the northern part of the Iberian Peninsula, whereas the southern species, T. pygmaeus, is found in the southwestern corner of the Iberian Peninsula. We study the intraspecific genetic differentiation of the group because morphological data show geographical variation and because the Iberian Peninsula is a recognized center of speciation and intraspecific genetic diversity for all kinds of organisms, amphibians included. We use target enrichment by sequence capture to generate c. 7 k nuclear DNA markers. We observe limited genetic exchange between the species, which confirms their distinctiveness. Both species show substantial genetic structuring that is only in part mirrored by morphological variation. Genetically differentiated groups are found in the south (T. marmoratus) and west (T. pygmaeus) of the species ranges. Our observations highlight the position of the Iberian Peninsula as a hotspot for genetic differentiation.

Piecing the barcoding puzzle of Palearctic water frogs (Pelophylax) sheds light on amphibian biogeography and global invasions

Palearctic water frogs (genus Pelophylax) are an outstanding model in ecology and evolution, being widespread, speciose, either threatened or threatening to other species through biological invasions, and capable of siring hybrid offspring that escape the rules of sexual reproduction. Despite half a century of genetic research and hundreds of publications, the diversity, systematics and biogeography of Pelophylax still remain highly confusing, in no small part due to a lack of correspondence between studies. To provide a comprehensive overview, we gathered >13,000 sequences of barcoding genes from >1700 native and introduced localities and built multigene mitochondrial (~17 kb) and nuclear (~10 kb) phylogenies. We mapped all currently recognized taxa and their phylogeographic lineages (>40) to get a grasp on taxonomic issues, cyto-nuclear discordances, the genetic makeup of hybridogenetic hybrids, and the origins of introduced populations. Competing hypotheses for the molecular calibration were evaluated through plausibility tests, implementing a new approach relying on predictions from the anuran speciation continuum. Based on our timetree, we propose a new biogeographic paradigm for the Palearctic since the Paleogene, notably by attributing a prominent role to the dynamics of the Paratethys, a vast paleo-sea that extended over most of Europe. Furthermore, our results show that distinct marsh frog lineages from Eastern Europe, the Balkans, the Near East, and Central Asia (P. ridibundus ssp.) are naturally capable of inducing hybridogenesis with pool frogs (P. lessonae). We identified 14 alien lineages (mostly of P. ridibundus) over ~20 areas of invasions, especially in Western Europe, with genetic signatures disproportionally pointing to the Balkans and Anatolia as the regions of origins, in line with exporting records of the frog leg industry and the stocks of pet sellers. Pelophylax thus emerges as one of the most invasive amphibians worldwide, and deserves much higher conservation concern than currently given by the authorities fighting biological invasions.

Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes)

The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide the most comprehensive phylogeographic framework for the A. obstetricans complex to date, as well as the first fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.

Performance of intron 7 of the β-fibrinogen gene for phylogenetic analysis: An example using gladiator frogs, Boana Gray, 1825 (Anura, Hylidae, Cophomantinae)

Boana, the third largest genus of Hylinae, has cryptic morphological species. The potential applicability of b-fibrinogen intron 7 – FGBI7 is explored to propose a robust phylogeny of Boana. The phylogenetic potential of FGBI7 was evaluated using maximum parsimony, MrBayes, and maximum likelihood analysis. Comparison of polymorphic sites and topologies obtained with concatenated analysis of FGBI7 and other nuclear genes (CXCR4, CXCR4, RHO, SIAH1, TYR, and 28S) allowed evaluation of the phylogenetic signal of FGBI7. Mean evolutionary rates were calculated using the sequences of the mitochondrial genes ND1 and CYTB available for Boana in GenBank. Dating of Boana and some of its groups was performed using the RelTime method with secondary calibration. FGBI7 analysis revealed high values at informative sites for parsimony. The absolute values of the mean evolutionary rate were higher for mitochondrial genes than for FGBI7. Dating of congruent Boana groups for ND1, CYTB, and FGBI7 revealed closer values between mitochondrial genes and slightly different values from those of FGBI7. Divergence times of basal groups tended to be overestimated when mtDNA was used and were more accurate when nDNA was used. Although there is evidence of phylogenetic potential arising from concatenation of specific genes, FGBI7 provides well-resolved independent gene trees. These results lead to a paradigm for linking data in phylogenomics that focuses on the uniqueness of species histories and ignores the multiplicities of individual gene histories.

Morphological diversification of Mediterranean anurans: the roles of evolutionary history and climate

Investigation of the ecological and evolutionary mechanisms governing the origin and diversification of species requires integrative approaches that often have to accommodate strong discordance among datasets. A common source of conflict is the combination of morphological and molecular characters with different evolutionary rates. Resolution of these discordances is crucial to assess the relative roles of different processes in generating and maintaining biodiversity. Anuran amphibians provide many examples of morphologically similar, genetically divergent lineages, posing questions about the relative roles of phylogeny and ecological factors in phenotypic evolution. We focused on three circum-Mediterranean anuran genera (Hyla, Alytes and Discoglossus), characterizing morphological and environmental disparity and comparing diversity patterns across biological levels of organization. Using a comparative phylogenetic framework, we tested how shared ancestry and climatic factors come together to shape phenotypic diversity. We found higher morphological differentiation within Hyla and Alytes than in Discoglossus. Body size and limb morphology contributed most to inter- and intraspecific morphological variation in Hyla and Alytes, but there was no strong phylogenetic signal, indicating that shared ancestry does not predict patterns of phenotypic divergence. In contrast, we uncovered a significant association between morphology and climatic descriptors, supporting the hypothesis that morphological disparity between species results from adaptive evolution.

Independent evolutionary transitions to pueriparity across multiple timescales in the viviparous genus Salamandra

The ability to bear live offspring, viviparity, has evolved multiple times across the tree of life and is a remarkable adaptation with profound life-history and ecological implications. Within amphibians the ancestral reproductive mode is oviparity followed by a larval life stage, but viviparity has evolved independently in all three amphibian orders. Two types of viviparous reproduction can be distinguished in amphibians; larviparity and pueriparity. Larviparous amphibians deliver larvae into nearby ponds and streams, while pueriparous amphibians deliver fully developed juveniles and thus do not require waterbodies for reproduction. Among amphibians, the salamander genus Salamandra is remarkable as it exhibits both inter- and intraspecific variation in the occurrence of larviparity and pueriparity. While the evolutionary relationships among Salamandra lineages have been the focus of several recent studies, our understanding of how often and when transitions between modes occurred is still incomplete. Furthermore, in species with intraspecific variation, the reproductive mode of a given population can only be confirmed by direct observation of births and thus the prevalence of pueriparous populations is also incompletely documented. We used sequence capture to obtain 1,326 loci from 94 individuals from across the geographic range of the genus, focusing on potential reproductive mode transition zones. We also report additional direct observations of pueriparous births for 20 new locations and multiple lineages. We identify at least five independent transitions from the ancestral mode of larviparity to pueriparity among and within species, occurring at different evolutionary timescales ranging from the Pliocene to the Holocene. Four of these transitions occurred within species. Based on a distinct set of markers and analyses, we also confirm previous findings of introgression between species and the need for taxonomic revisions in the genus. We discuss the implications of our findings with respect to the evolution of this complex trait, and the potential of using five independent convergent transitions for further studies on the ecological context in which pueriparity evolves and the genetic architecture of this specialized reproductive mode.

Karyological Diversification in the Genus Lyciasalamandra (Urodela: Salamandridae)

Simple Summary

The Lycian salamanders of the genus Lyciasalamandra are characterized by a debated taxonomy and phylogenetic relationships. They have been the subject of various molecular and phylogenetic analyses, but their chromosomal diversity is completely unknown. We here present a comparative cytogenetic analysis on five out of the seven described species and seven subspecies of Lyciasalamandra, providing the first karyological assessment on the genus and comparing them to closely related representatives of the genus Salamandra. We analyzed the occurrence and distribution of different conserved (chromosome number and morphology) and highly variable karyological features. We found an impressive diversity in the configuration of nucleolus organizing regions (NORs), which alternatively occur either as heteromorphic or homomorphic loci on distinct regions of different chromosome pairs. We highlight that the observed peculiar taxon-specific pattern of chromosome markers supports the taxonomic validity of the different studied evolutionary lineages and is consistent with a scenario of synchronous evolution in the Lycian salamanders.

Abstract

We performed the first cytogenetic analysis on five out of the seven species of the genus Lyciasalamandra, including seven subspecies, and representatives of its sister genus Salamandra. All the studied species have a similar karyotype of 2n = 24, mostly composed of biarmed elements. C-bands were observed on all chromosomes, at centromeric, telomeric and interstitial position. We found a peculiar taxon-specific NOR configuration, including either heteromorphic and homomorphic NORs on distinct regions of different chromosomes. Lyciasalamandra a.antalyana and L. helverseni showed two homomorphic NORs (pairs 8 and 2, respectively), while heteromorphic NORs were found in L. billae (pairs 6, 12), L. flavimembris (pairs 2, 12), L. l. luschani (pairs 2, 12), L. l. basoglui (pairs 6, 12), L. l. finikensis (pairs 2, 6) and S. lanzai (pairs 8, 10). Homomorphic NORs with an additional supernumerary site were shown by S. s. salamandra (pairs 2, 8) and S. s. gigliolii (pairs 2, 10). This unexpected highly variable NOR configuration is probably derived from multiple independent NOR translocations and paracentric inversions and correlated to lineage divergence in Lyciasalamandra. These results support the taxonomic validity of the studied taxa and are consistent with a hypothesized scenario of synchronous evolution in the genus.

Landscape resistance constrains hybridization across contact zones in a reproductively and morphologically polymorphic salamander

Explicitly accounting for phenotypic differentiation together with environmental heterogeneity is crucial to understand the evolutionary dynamics in hybrid zones. Species showing intra-specific variation in phenotypic traits that meet across environmentally heterogeneous regions constitute excellent natural settings to study the role of phenotypic differentiation and environmental factors in shaping the spatial extent and patterns of admixture in hybrid zones. We studied three environmentally distinct contact zones where morphologically and reproductively divergent subspecies of Salamandra salamandra co-occur: the pueriparous S. s. bernardezi that is mostly parapatric to its three larviparous subspecies neighbours. We used a landscape genetics framework to: (i) characterise the spatial location and extent of each contact zone; (ii) assess patterns of introgression and hybridization between subspecies pairs; and (iii) examine the role of environmental heterogeneity in the evolutionary dynamics of hybrid zones. We found high levels of introgression between parity modes, and between distinct phenotypes, thus demonstrating the evolution to pueriparity alone or morphological differentiation do not lead to reproductive isolation between these highly divergent S. salamandra morphotypes. However, we detected substantial variation in patterns of hybridization across contact zones, being lower in the contact zone located on a topographically complex area. We highlight the importance of accounting for spatial environmental heterogeneity when studying evolutionary dynamics of hybrid zones.

Slow and steady wins the race: contrasted phylogeographic signatures in two Alpine amphibians

A deeper phylogeographic structure is expected for slow-dispersing habitat specialists compared to widespread adaptable species, especially in topographically complex regions. We tested this classic assumption by comparing the genomic (RAD-sequencing) phylogeographies of two amphibians inhabiting the Swiss Alps: the mobile, cosmopolitan common frog (Rana temporaria) against the stationary, mountain endemic Alpine salamander (Salamandra atra). Our results ran opposite of predictions: the frog displayed significantly higher genetic divergences and lower within-population variation compared to the salamander. This implies a prominent role for their distinctive glacial histories in shaping intraspecific diversity and structure: diversification and recolonization from several circum-Alpine micro-refugia for the frog versus a single refugium for the salamander, potentially combined with better population connectivity and stability. These striking differences emphasize the great variability of phylogeographic responses to the Quaternary glaciations, hence the complexity to predict general patterns of genetic diversity at the regional scale, and the forces that underlie them.

Phylogenomic inference of species and subspecies diversity in the Palearctic salamander genus Salamandra

The salamander genus Salamandra is widespread across Europe, North Africa, and the Near East and is renowned for its conspicuous and polymorphic colouration and diversity of reproductive modes. The phylogenetic relationships within the genus, and especially in the highly polymorphic species S. salamandra, have been very challenging to elucidate, leaving its real evolutionary history and classification at species and subspecies levels a topic of debate and contention. However, the distribution of diversity and species delimitation within the genus are critically important for identifying evolutionarily significant units for conservation and management, especially in light of threats posed by the pathogenic chytrid fungus Batrachochytrium salamandrivorans that is causing massive declines of S. salamandra populations in central Europe. Here, we conducted a phylogenomic analysis from across the taxonomic and geographic breadth of the genus Salamandra in its entire range. Bayesian, maximum likelihood and network-based phylogenetic analyses of up to 4905 ddRADseq-loci (294,300 nucleotides of sequence) supported the distinctiveness of all currently recognised species (Salamandra algira, S. atra, S. corsica, S. infraimmaculata, S. lanzai, and S. salamandra), and all five species for which we have multiple exemplars were confirmed as monophyletic. Within S. salamandra, two main clades can be distinguished: one clade with the Apenninic subspecies S. s. gigliolii nested within the Iberian S. s. bernardezi/fastuosa; and a second clade comprising all other Iberian, Central and East European subspecies. Our analyses revealed that some of the currently recognized subspecies of S. salamandra are paraphyletic and may require taxonomic revision, with the Central- and Eastern-European subspecies all being poorly differentiated at the analysed genomic markers. Salamandra s. longirostris - sometimes considered a separate species - was nested within S. salamandra, consistent with its subspecies status. The relationships identified within and between Salamandra species provide valuable context for future systematic and biogeographic studies, and help elucidate critical evolutionary units for conservation and taxonomy.

Diversification in Qinghai-Tibet Plateau: Orchidinae (Orchidaceae) clades exhibiting pre-adaptations play critical role

We explore the origins of the extraordinary plant diversity in the Qinghai-Tibetan Plateau (QTP) using Orchidinae (Orchidaceae) as a model. Our results indicate that six major clades in Orchidinae exhibited substantial variation in the temporal and spatial sequence of diversification. Our time-calibrated phylogenetic model suggests that the species-richness of Orchidinae arose through a combination of in situ diversification, colonisation, and local recruitment. There are multiple origins of species-richness of Orchidinae in the QTP, and pre-adaptations in clades from North Temperate and alpine regions were crucial for in situ diversification. The geographic analysis identified 29 dispersals from Asia, Africa and Europe into the QTP and 15 dispersals out. Most endemic species of Orchidinae evolved within the past six million years.

Unequal sisters - Past and potential future range development of Anatolian and Hyrcanian brown frogs

Phylogeography can reconstruct historical evolutionary processes by comparing historical patterns of gene flow, divergence among species and by using species distribution models (SDM) upon geographic distribution. We investigate the phylogeographic patterns of Anatolian brown frogs including R. macrocnemis and R. tavasensis as well as the Hyrcanian brown frog, R. pseudodalmatina, using a fragment of the mitochondrial 16S rRNA gene for 145 specimens across the entire range of these frogs. We calculate parameters of molecular diversity, such as the number of variable sites (S), the number of haplotypes (h), haplotype diversity (Hd) and nucleotide diversity (π). We generated a haplotype network and used three methods (Neutrality tests, mismatch distributions and Bayesian skyline plots) to reconstruct the demographic histories of R. macrocnemis and R. pseudodalmatina. Finally, we used SDMs to predict the habitat suitability for three periods: The Present Day, the Last Glacial Maximum (LGM) and the future until 2070 for R. macrocnemis and R. pseudodalmatina. Our phylogenetic analyses support a late Miocene origin of Anatolian and Hyrcanian lineages. Hyrcanian brown frogs were enclosed in lowlands of the southern coast of the Caspian Sea after the uplift of the Elburz range and the Armenian plateau. The formation of a salinity belt from the north Aegean corridor (the south western Turkey) to northward during the Late Tortonian led to the subdivision of ancestor of the Anatolian lineage into today isolated western and eastern populations. The salinity belt had a considerable impact on the divergence of R. tavasensis from R. macrocnemis. Combined historical demographic analyses and SDMs revealed a rapid expansion occurring during the Pleistocene in R. macrocnemis and R. pseudodalmatina. Currently, suitable habitat for R. macrocnemis has declined compared to the LGM, and the species is predicted to do even worse under future climatic conditions. In contrast, R. pseudodalmatina found suitable habitat from the LGM to present within its restricted distribution area; it is predicted to do fine even under future climate.

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Subdivided Pleistocene glacial refugia, best known as "refugia within refugia", provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and "ghost" mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD-sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply-diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation-by-distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.

Integrating hybrid zone analyses in species delimitation: lessons from two anuran radiations of the Western Mediterranean

Molecular ecologists often rely on phylogenetic evidence for assessing the species-level systematics of newly discovered lineages. Alternatively, the extent of introgression at phylogeographic transitions can provide a more direct test to assign candidate taxa into subspecies or species categories. Here, we compared phylogenetic versus hybrid zone approaches of species delimitation in two groups of frogs from the Western Mediterranean region (Discoglossus and Pelodytes), by using genomic data (ddRAD). In both genera, coalescent analyses recovered almost all nominal taxa as "species". However, the least-diverged pairs D. g. galganoi/jeanneae and P. punctatus/hespericus admix over hundreds of kilometers, suggesting that they have not yet developed strong reproductive isolation and should be treated as conspecifics. In contrast, the comparatively older D. scovazzi/pictus and P. atlanticus/ibericus form narrow contact zones, consistent with species distinctiveness. Due to their complementarity, we recommend taxonomists to combine phylogenomics with hybrid zone analyses to scale the gray zone of speciation, i.e., the evolutionary window separating widely admixing lineages versus nascent reproductively isolated species. The radically different transitions documented here conform to the view that genetic incompatibilities accumulating with divergence generate a weak barrier to gene flow for long periods of time, until their effects multiply and the speciation process then advances rapidly. Given the variability of the gray zone among taxonomic groups, at least from our current abilities to measure it, we recommend to customize divergence thresholds within radiations to categorize lineages for which no direct test of speciation is possible.

Phylogeographic analyses point to long-term survival on the spot in micro-endemic Lycian salamanders

Lycian salamanders (genus Lyciasalamandra) constitute an exceptional case of micro-endemism of an amphibian species on the Asian Minor mainland. These viviparous salamanders are confined to karstic limestone formations along the southern Anatolian coast and some islands. We here study the genetic differentiation within and among 118 populations of all seven Lyciasalamandra species across the entire genus’ distribution. Based on circa 900 base pairs of fragments of the mitochondrial 16SrDNA and ATPase genes, we analysed the spatial haplotype distribution as well as the genetic structure and demographic history of populations. We used 253 geo-referenced populations and CHELSA climate data to infer species distribution models which we projected on climatic conditions of the Last Glacial Maximum (LGM). Within all but one species, distinct phyloclades were identified, which only in parts matched current taxonomy. Most haplotypes (78%) were private to single populations. Sometimes population genetic parameters showed contradicting results, although in several cases they indicated recent population expansion of phyloclades. Climatic suitability of localities currently inhabited by salamanders was significantly lower during the LGM compared to recent climate. All data indicated a strong degree of isolation among Lyciasalamandra populations, even within phyloclades. Given the sometimes high degree of haplotype differentiation between adjacent populations, they must have survived periods of deteriorated climates during the Quaternary on the spot. However, the alternative explanation of male biased dispersal combined with a pronounced female philopatry can only be excluded if independent nuclear data confirm this result.

An amphibian species pushed out of Britain by a moving hybrid zone

Classical theory states that hybrid zones will be stable in troughs of low population density where dispersal is hampered. Yet, evidence for moving hybrid zones is mounting. One possible reason that moving zones have been underappreciated is that they may drive themselves into oblivion and with just the superseding species remaining, morphological and genetic signals of past species replacement may be difficult to appreciate. Using genetic data (32 diagnostic single nucleotide polymorphisms) from a clinal hybrid zone of the common toad (Bufo bufo) and the spined toad (Bufo spinosus) in France for comparison, alleles of the latter species were documented in common toads in the south of Great Britain, at frequencies in excess of 10%. Because long distance dispersal across the Channel is unlikely, the conclusion reached was that the continental toad hybrid zone which previously extended into Britain, moved southwards and extirpated B. spinosus. Species distribution models for the mid-Holocene and the present support that climate has locally changed in favour of B. bufo. The system bears resemblance with the demise of Homo neanderthalensis and the rise of Homo sapiens and provides an example that some paleoanthropologists demanded in support of a hominin "leaky replacement" scenario. The toad example is informative just because surviving pure B. spinosus and an extant slowly moving interspecific hybrid zone are available for comparison.

Call a spade a spade: taxonomy and distribution of Pelobates, with description of a new Balkan endemic

The genomic era contributes to update the taxonomy of many debated terrestrial vertebrates. In an accompanying work, we provided a comprehensive molecular assessment of spadefoot toads (Pelobates) using genomic data. Our results call for taxonomic updates in this group. First, nuclear phylogenomics confirmed the species-level divergence between the Iberian P.cultripes and its Moroccan relative P.varaldii. Second, we inferred that P.fuscus and P.vespertinus, considered subspecies until recently, feature partial reproductive isolation and thus deserve a specific level. Third, we evidenced cryptic speciation and diversification among deeply diverged lineages collectively known as Pelobatessyriacus. Populations from the Near East correspond to the Eastern spadefoot toad P.syriacus sensu stricto, which is represented by two subspecies, one in the Levant (P.s.syriacus) and the other in the rest of the range (P.s.boettgeri). Populations from southeastern Europe correspond to the Balkan spadefoot toad, P.balcanicus. Based on genetic evidence, this species is also polytypic: the nominal P.b.balcanicus inhabits the Balkan Peninsula; a new subspecies P.b.chloeaessp. nov. appears endemic to the Peloponnese. In this paper, we provide an updated overview of the taxonomy and distribution of all extant Pelobates taxa and describe P.b.chloeaessp. nov.

Phylogeography of a cryptic speciation continuum in Eurasian spadefoot toads (Pelobates)

Cryptic phylogeographic diversifications provide unique models to examine the role of phylogenetic divergence on the evolution of reproductive isolation, without extrinsic factors such as ecological and behavioural differentiation. Yet, to date very few comparative studies have been attempted within such radiations. Here, we characterize a new speciation continuum in a group of widespread Eurasian amphibians, the Pelobates spadefoot toads, by conducting multilocus (restriction site associated DNA sequencing and mitochondrial DNA) phylogenetic, phylogeographic and hybrid zone analyses. Within the P. syriacus complex, we discovered species-level cryptic divergences (>5 million years ago [My]) between populations distributed in the Near-East (hereafter P. syriacus sensu stricto [s.s.]) and southeastern Europe (hereafter P. balcanicus), each featuring deep intraspecific lineages. Altogether, we could scale hybridizability to divergence time along six different stages, spanning from sympatry without gene flow (P. fuscus and P. balcanicus, >10 My), parapatry with highly restricted hybridization (P. balcanicus and P. syriacus s.s., >5 My), narrow hybrid zones (~15 km) consistent with partial reproductive isolation (P. fuscus and P. vespertinus, ~3 My), to extensive admixture between Pleistocene and refugial lineages (≤2 My). This full spectrum empirically supports a gradual build up of reproductive barriers through time, reversible up until a threshold that we estimate at ~3 My. Hence, cryptic phylogeographic lineages may fade away or become reproductively isolated species simply depending on the time they persist in allopatry, and without definite ecomorphological divergence.