Assessing fine-scale pondscape connectivity with amphibian eyes: An integrative approach using genomic and capture-mark-recapture data.

In the face of habitat loss, preserving functional connectivity is essential to maintain genetic diversity and the demographic dynamics required for the viability of biotic communities. This requires knowledge of the dispersal behaviour of target species, which can be modelled as kernels, or probability density functions of dispersal distances at increasing geographic distances. We present an integrative approach to investigate the relationships between genetic connectivity and demographic parameters in organisms with low vagility focusing on five syntopic pond-breeding amphibians. We genotyped 1056 individuals of two anuran and three urodele species (1732-3913 SNPs per species) from populations located in a landscape comprising 64 ponds to characterize fine-scale genetic structure in a comparative framework, and combined these genetic data with information obtained in a previous 2-year capture-mark-recapture (CMR) study. Specifically, we contrasted graphs reconstructed from genomic data with connectivity graphs based on dispersal kernels and demographic information obtained from CMR data from previous studies, and assessed the effects of population size, population density, geographical distances, inverse movement probabilities and the presence of habitat patches potentially functioning as stepping stones on genetic differentiation. Our results show a significant effect of local population sizes on patterns of genetic differentiation at small spatial scales. In addition, movement records and cluster-derived kernels provide robust inferences on most likely dispersal paths that are consistent with genomic inferences on genetic connectivity. The integration of genetic and CMR data holds great potential for understanding genetic connectivity at spatial scales relevant to individual organisms, with applications for the implementation of management actions at the landscape level.

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.

Mass of genes rather than master genes underlie the genomic architecture of amphibian speciation.

The genetic architecture of speciation, i.e., how intrinsic genomic incompatibilities promote reproductive isolation (RI) between diverging lineages, is one of the best-kept secrets of evolution. To directly assess whether incompatibilities arise in a limited set of large-effect speciation genes, or in a multitude of loci, we examined the geographic and genomic landscapes of introgression across the hybrid zones of 41 pairs of frog and toad lineages in the Western Palearctic region. As the divergence between lineages increases, phylogeographic transitions progressively become narrower, and larger parts of the genome resist introgression. This suggests that anuran speciation proceeds through a gradual accumulation of multiple barrier loci scattered across the genome, which ultimately deplete hybrid fitness by intrinsic postzygotic isolation, with behavioral isolation being achieved only at later stages. Moreover, these loci were disproportionately sex linked in one group (Hyla) but not in others (Rana and Bufotes), implying that large X-effects are not necessarily a rule of speciation with undifferentiated sex chromosomes. The highly polygenic nature of RI and the lack of hemizygous X/Z chromosomes could explain why the speciation clock ticks slower in amphibians compared to other vertebrates. The clock-like dynamics of speciation combined with the analytical focus on hybrid zones offer perspectives for more standardized practices of species delimitation.

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 a comprehensive phylogeographic framework for the A. obstetricans complex, as well as a 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.

Coevolution between MHC Class I and Antigen-Processing Genes in Salamanders.

Proteins encoded by antigen-processing genes (APGs) provide major histocompatibility complex (MHC) class I (MHC-I) with antigenic peptides. In mammals, polymorphic multigenic MHC-I family is served by monomorphic APGs, whereas in certain nonmammalian species both MHC-I and APGs are polymorphic and coevolve within stable haplotypes. Coevolution was suggested as an ancestral gnathostome feature, presumably enabling only a single highly expressed classical MHC-I gene. In this view coevolution, while optimizing some aspects of adaptive immunity, would also limit its flexibility by preventing the expansion of classical MHC-I into a multigene family. However, some nonmammalian taxa, such as salamanders, have multiple highly expressed MHC-I genes, suggesting either that coevolution is relaxed or that it does not prevent the establishment of multigene MHC-I. To distinguish between these two alternatives, we use salamanders (30 species from 16 genera representing six families) to test, within a comparative framework, a major prediction of the coevolution hypothesis: the positive correlation between MHC-I and APG diversity. We found that MHC-I diversity explained both within-individual and species-wide diversity of two APGs, TAP1 and TAP2, supporting their coevolution with MHC-I, whereas no consistent effect was detected for the other three APGs (PSMB8, PSMB9, and TAPBP). Our results imply that although coevolution occurs in salamanders, it does not preclude the expansion of the MHC-I gene family. Contrary to the previous suggestions, nonmammalian vertebrates thus may be able to accommodate diverse selection pressures with flexibility granted by rapid expansion or contraction of the MHC-I family, while retaining the benefits of coevolution between MHC-I and TAPs.

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.

Evaluating surrogates of genetic diversity for conservation planning.

Protected-area systems should conserve intraspecific genetic diversity. Because genetic data require resources to obtain, several approaches have been proposed for generating plans for protected-area systems (prioritizations) when genetic data are not available. Yet such surrogate-based approaches remain poorly tested. We evaluated the effectiveness of potential surrogate-based approaches based on microsatellite genetic data collected across the Iberian Peninsula for 7 amphibian and 3 reptilian species. Long-term environmental suitability did not effectively represent sites containing high genetic diversity (allelic richness). Prioritizations based on long-term environmental suitability had similar performance to random prioritizations. Geographic distances and resistance distances based on contemporary environmental suitability were not always effective surrogates for identification of combinations of sites that contain individuals with different genetic compositions. Our results demonstrate that population genetic data based on commonly used neutral markers can inform prioritizations, and we could not find an adequate substitute. Conservation planners need to weigh the potential benefits of genetic data against their acquisition costs.

Evaluación de los Sustitutos de la Diversidad Genética para la Planeación de la Conservación Resumen Los sistemas de áreas protegidas deberían conservar la diversidad genética intraespecífica. Ya que para obtener datos genéticos se requieren recursos, se han propuesto distintas estrategias para generar los planes para los sistemas de áreas protegidas (priorizaciones) cuando los datos genéticos no están disponibles. A pesar de lo anterior, dichas estrategias basadas en sustitutos han sido poco evaluadas. Evaluamos la efectividad del potencial de las estrategias basadas en sustitutos cuya base son los datos genéticos de microsatélites obtenidos en toda la Península Ibérica y correspondientes a siete especies de anfibios y a tres de reptiles. La idoneidad ambiental a largo plazo no representó efectivamente los sitios que contienen una diversidad genética alta (riqueza de alelos). Las priorizaciones basadas en la idoneidad ambiental a largo plazo tuvieron un desempeño similar a las priorizaciones aleatorias. Las distancias geográficas y las distancias de resistencia basadas en la idoneidad ambiental contemporánea no siempre fueron sustitutos efectivos para la identificación de las combinaciones de sitios que contienen individuos con composiciones genéticas diferentes. Nuestros resultados demuestran que los datos genéticos de una población basados en marcadores neutrales de uso común pueden informar a las priorizaciones y que no pudimos encontrar un sustituto adecuado. Los planificadores de la conservación necesitan sopesar los beneficios potenciales de los datos genéticos contra sus costos de adquisición.

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.

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt.

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south-western and south-eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south-western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south-eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150-km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2-27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population-level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.

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.

Telomere attrition with age in a wild amphibian population.

Telomere shortening with age has been documented in many organisms, but few studies have reported telomere length measurements in amphibians, and no information is available for growth after metamorphosis, nor in wild populations. We provide both cross-sectional and longitudinal evidence of net telomere attrition with age in a wild amphibian population of natterjack toads (Epidalea calamita). Based on age-estimation by skeletochronology and qPCR telomere length measurements in the framework of an individual-based monitoring programme, we confirmed telomere attrition in recaptured males. Our results support that toads experience telomere attrition throughout their ontogeny, and that most attrition occurs during the first 1-2 years. We did not find associations between telomere length and inbreeding or body condition. Our results on telomere length dynamics under natural conditions confirm telomere shortening with age in amphibians and provide quantification of wide telomere length variation within and among age-classes in a wild breeding population.

Comparative landscape genetics of pond-breeding amphibians in Mediterranean temporal wetlands: The positive role of structural heterogeneity in promoting gene flow.

Comparative landscape genetics studies can provide key information to implement cost-effective conservation measures favouring a broad set of taxa. These studies are scarce, particularly in Mediterranean areas, which include diverse but threatened biological communities. Here, we focus on Mediterranean wetlands in central Iberia and perform a multi-level, comparative study of two endemic pond-breeding amphibians, a salamander (Pleurodeles waltl) and a toad (Pelobates cultripes). We genotyped 411 salamanders from 20 populations and 306 toads from 16 populations at 18 and 16 microsatellite loci, respectively, and identified major factors associated with population connectivity through the analysis of three sets of variables potentially affecting gene flow at increasingly finer levels of spatial resolution. Topographic, land use/cover, and remotely sensed vegetation/moisture indices were used to derive optimized resistance surfaces for the two species. We found contrasting patterns of genetic structure, with stronger, finer scale genetic differentiation in Pleurodeles waltl, and notable differences in the role of fine-scale patterns of heterogeneity in vegetation cover and water content in shaping patterns of regional genetic structure in the two species. Overall, our results suggest a positive role of structural heterogeneity in population connectivity in pond-breeding amphibians, with habitat patches of Mediterranean scrubland and open oak woodlands ("dehesas") facilitating gene flow. Our study highlights the usefulness of remotely sensed continuous variables of land cover, vegetation and water content (e.g., NDVI, NDMI) in conservation-oriented studies aimed at identifying major drivers of population connectivity.

Hybrid zone formation and contrasting outcomes of secondary contact over transects in common toads.

Much progress in speciation research stems from documenting patterns of morphological and genetic variation in hybrid zones. Contrasting patterns of marker introgression in different sections of the contact can provide valuable insights on the relative importance of various evolutionary mechanisms maintaining species differences in the face of hybridization and gene flow and on hybrid zone temporal and spatial dynamics. We studied species interactions in the common toads Bufo bufo and B. spinosus in France and northwestern Italy using morphological and molecular data from the mitochondrial and nuclear genomes in an extensive survey, including two independent transects west and east of the Alps. At both, we found sharp, coincident and concordant nuclear genetic transitions. However, morphological clines were wider or absent and mtDNA introgression was asymmetric. We discuss alternative, nonexclusive hypotheses about evolutionary processes generating these patterns, including drift, selection, long-distance dispersal and spatial shifts in hybrid zone location and structure. The distribution of intraspecific mtDNA lineages supports a scenario in which B. bufo held a local refugium during the last glacial maximum. Present-day genetic profiles are best explained by an advance of B. spinosus from a nearby Iberian refugium, largely superseding the local B. bufo population, followed by an advance of B. bufo from the Balkans, with prongs north and south of the Alps, driving B. spinosus southwards. A pendulum moving hybrid zone, first northwards and then southwards, explains the wide areas of introgression at either side of the current position of the contact zones.

Integrative inference of population history in the Ibero-Maghrebian endemic Pleurodeles waltl (Salamandridae).

Inference of population histories from the molecular signatures of past demographic processes is challenging, but recent methodological advances in species distribution models and their integration in time-calibrated phylogeographic studies allow detailed reconstruction of complex biogeographic scenarios. We apply an integrative approach to infer the evolutionary history of the Iberian ribbed newt (Pleurodeles waltl), an Ibero-Maghrebian endemic with populations north and south of the Strait of Gibraltar. We analyzed an extensive multilocus dataset (mitochondrial and nuclear DNA sequences and ten polymorphic microsatellite loci) and found a deep east-west phylogeographic break in Iberian populations dating back to the Plio-Pleistocene. This break is inferred to result from vicariance associated with the formation of the Guadalquivir river basin. In contrast with previous studies, North African populations showed exclusive mtDNA haplotypes, and formed a monophyletic clade within the Eastern Iberian lineage in the mtDNA genealogy. On the other hand, microsatellites failed to recover Moroccan populations as a differentiated genetic cluster. This is interpreted to result from post-divergence gene flow based on the results of IMA2 and Migrate analyses. Thus, Moroccan populations would have originated after overseas dispersal from the Iberian Peninsula in the Pleistocene, with subsequent gene flow in more recent times, implying at least two trans-marine dispersal events. We modeled the distribution of the species and of each lineage, and projected these models back in time to infer climatically favourable areas during the mid-Holocene, the last glacial maximum (LGM) and the last interglacial (LIG), to reconstruct more recent population dynamics. We found minor differences in climatic favourability across lineages, suggesting intraspecific niche conservatism. Genetic diversity was significantly correlated with the intersection of environmental favourability in the LIG and LGM, indicating that populations of P. waltl are genetically more diverse in regions that have remained environmentally favourable through the last glacial cycle, particularly southern Iberia and northern Morocco. This study provides novel insights into the relative roles of geology and climate on the biogeography of a biodiversity hotspot.

Effects of Sample Size and Full Sibs on Genetic Diversity Characterization: A Case Study of Three Syntopic Iberian Pond-Breeding Amphibians.

Accurate characterization of genetic diversity is essential for understanding population demography, predicting future trends and implementing efficient conservation policies. For that purpose, molecular markers are routinely developed for nonmodel species, but key questions regarding sampling design, such as calculation of minimum sample sizes or the effect of relatives in the sample, are often neglected. We used accumulation curves and sibship analyses to explore how these 2 factors affect marker performance in the characterization of genetic diversity. We illustrate this approach with the analysis of an empirical dataset including newly optimized microsatellite sets for 3 Iberian amphibian species: Hyla molleri, Epidalea calamita, and Pelophylax perezi. We studied 17-21 populations per species (total n = 547, 652, and 516 individuals, respectively), including a reference locality in which the effect of sample size was explored using larger samples (77-96 individuals). As expected, FIS and tests for Hardy-Weinberg equilibrium and linkage disequilibrium were affected by the presence of full sibs, and most initially inferred disequilibria were no longer statistically significant when full siblings were removed from the sample. We estimated that to obtain reliable estimates, the minimum sample size (potentially including full sibs) was close to 20 for expected heterozygosity, and between 50 and 80 for allelic richness. Our pilot study based on a reference population provided a rigorous assessment of marker properties and the effects of sample size and presence of full sibs in the sample. These examples illustrate the advantages of this approach to produce robust and reliable results for downstream analyses.

Hybridization during altitudinal range shifts: nuclear introgression leads to extensive cyto-nuclear discordance in the fire salamander.

Ecological models predict that, in the face of climate change, taxa occupying steep altitudinal gradients will shift their distributions, leading to the contraction or extinction of the high-elevation (cold-adapted) taxa. However, hybridization between ecomorphologically divergent taxa commonly occurs in nature and may lead to alternative evolutionary outcomes, such as genetic merger or gene flow at specific genes. We evaluate this hypothesis by studying patterns of divergence and gene flow across three replicate contact zones between high- and low-elevation ecomorphs of the fire salamander (Salamandra salamandra) that have experienced altitudinal range shifts over the current postglacial period. Strong population structure with high genetic divergence in mitochondrial DNA suggests that vicariant evolution has occurred over several glacial-interglacial cycles and that it has led to cryptic differentiation within ecomorphs. In current parapatric boundaries, we do not find evidence for local extinction and replacement upon postglacial expansion. Instead, parapatric taxa recurrently show discordance between mitochondrial and nuclear markers, suggesting nuclear-mediated gene flow across contact zones. Isolation with migration models support this hypothesis by showing significant gene flow across all five parapatric boundaries. Together, our results suggest that, while some genomic regions, such as the mitochondria, may follow morphologic species traits and retreat to isolated mountain tops, other genomic regions, such as nuclear markers, may flow across parapatric boundaries, sometimes leading to a complete genetic merger. We show that despite high ecologic and morphologic divergence over prolonged periods of time, hybridization allows for evolutionary outcomes alternative to extinction and replacement of taxa in response to climate change.

The effects of inference method, population sampling, and gene sampling on species tree inferences: an empirical study in slender salamanders (Plethodontidae: Batrachoseps).

Species tree methods are now widely used to infer the relationships among species from multilocus data sets. Many methods have been developed, which differ in whether gene and species trees are estimated simultaneously or sequentially, and in how gene trees are used to infer the species tree. While these methods perform well on simulated data, less is known about what impacts their performance on empirical data. We used a data set including five nuclear genes and one mitochondrial gene for 22 species of Batrachoseps to compare the effects of method of analysis, within-species sampling and gene sampling on species tree inferences. For this data set, the choice of inference method had the largest effect on the species tree topology. Exclusion of individual loci had large effects in *BEAST and STEM, but not in MP-EST. Different loci carried the greatest leverage in these different methods, showing that the causes of their disproportionate effects differ. Even though substantial information was present in the nuclear loci, the mitochondrial gene dominated the *BEAST species tree. This leverage is inherent to the mtDNA locus and results from its high variation and lower assumed ploidy. This mtDNA leverage may be problematic when mtDNA has undergone introgression, as is likely in this data set. By contrast, the leverage of RAG1 in STEM analyses does not reflect properties inherent to the locus, but rather results from a gene tree that is strongly discordant with all others, and is best explained by introgression between distantly related species. Within-species sampling was also important, especially in *BEAST analyses, as shown by differences in tree topology across 100 subsampled data sets. Despite the sensitivity of the species tree methods to multiple factors, five species groups, the relationships among these, and some relationships within them, are generally consistently resolved for Batrachoseps.

Concordant morphological and molecular clines in a contact zone of the Common and Spined toad (Bufo bufo and B. spinosus) in the northwest of France.

BACKGROUND: Hybrid zones are regions where individuals of two species meet and produce hybrid progeny, and are often regarded as natural laboratories to understand the process of species formation. Two microevolutionary processes can take place in hybrid zones, with opposing effects on population differentiation. Hybridization tends to produce genetic homogenization, reducing species differences, whereas the presence of mechanisms of reproductive isolation result in barriers to gene flow, maintaining or increasing differences between taxa. RESULTS: Here we study a contact zone between two hybridizing toad species, Bufo bufo and B. spinosus, through a combination of molecular (12 polymorphic microsatellites, four nuclear and two mitochondrial SNP markers) and morphological data in a transect in the northwest of France. The results show largely concordant clines across markers, defining a narrow hybrid zone of ca. 30 km wide. Most hybrids in the centre of the contact zone are classified as F2 or backcrossed individuals, with no individuals assigned to the F1 hybrid class. CONCLUSIONS: We discuss the implications of these results for our understanding of the evolutionary history of these species. We anticipate that the toad contact zone here described will become an important asset in the study of hybrid zone dynamics and evolutionary biology because of its easy access and the abundance of the species involved.

Molecular evidence for cryptic candidate species in Iberian Pelodytes (Anura, Pelodytidae).

Species delineation is a central topic in evolutionary biology, with current efforts focused on developing efficient analytical tools to extract the most information from molecular data and provide objective and repeatable results. In this paper we use a multilocus dataset (mtDNA and two nuclear markers) in a geographically comprehensive population sample across Iberia and Western Europe to delineate candidate species in a morphologically cryptic species group, Parsley frogs (genus Pelodytes). Pelodytes is the sole extant representative of an ancient, historically widely distributed anuran clade that currently includes three species: P. caucasicus in the Caucasus; P. punctatus in Western Europe, from Portugal to North-Western Italy; and P. ibericus in Southern Iberia. Phylogenetic analyses recovered four major well-supported haplotype clades in Western Europe, corresponding to well demarcated geographical subdivisions and exhibiting contrasting demographic histories. Splitting times date back to the Plio-Pleistocene and are very close in time. Species-tree analyses recovered one of these species lineages, corresponding to P. ibericus (lineage B), as the sister taxon to the other three major species lineages, distributed respectively in: western Iberian Peninsula, along the Atlantic coast and part of central Portugal (lineage A); Central and Eastern Spain (lineage C); and North-eastern Spain, France and North-western Italy (lineage D). The latter is in turn subdivided into two sub-clades, one in SE France and NW Italy and the other one from NE Spain to NW France, suggesting the existence of a Mediterranean-Atlantic corridor along the Garonne river. An information theory-based validation approach implemented in SpedeSTEM supports an arrangement of four candidate species, suggesting the need for a taxonomic revision of Western European Pelodytes.

Multilocus assessment of phylogenetic relationships in Alytes (Anura, Alytidae).

With the advent of large multilocus datasets, molecular systematics is experiencing very rapid progress, but important challenges remain regarding data analysis and interpretation. Midwife toads (genus Alytes) exemplify two of the most widespread problems for accurate phylogenetic reconstruction: discerning the causes of discordance between gene trees, and resolving short internodes produced during rapid, successive splitting events. The three species in subgenus Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), the sister group to A. obstetricans, have disjunct and highly restricted geographical ranges, which are thought to result from old vicariant events affecting their common ancestor, but their phylogenetic relationships are still unresolved. In this study we re-address the phylogeny of Alytes with a special focus on the relationships in Baleaphryne with a multilocus dataset including >9000 base pairs of mitochondrial DNA and four nuclear markers (3142bp) in all recognized taxa, including all subspecies of A. obstetricans. Both concatenation and species tree analyses suggest that A. muletensis, endemic to the Balearic island of Mallorca, is the sister taxon to a clade comprising the southeastern Iberian endemic A. dickhilleni and the North African A. maurus. This scenario is consistent with palaeogeological evidence associated with the fragmentation of the Betic-Rifean Massif, followed by the opening of the Strait of Gibraltar. On the other hand, analyses of intraspecific variation in A. obstetricans are inconclusive regarding relationships between major clades and conflict with current subspecific taxonomy.