Fifteen shades of green: The evolution of Bufotes toads revisited.

The radiation of Palearctic green toads (Bufotes) holds great potential to evaluate the role of hybridization in phylogeography at multiple stages along the speciation continuum. With fifteen species representing three ploidy levels, this model system is particularly attractive to examine the causes and consequences of allopolyploidization, a prevalent yet enigmatic pathway towards hybrid speciation. Despite substantial efforts, the evolutionary history of this species complex remains largely blurred by the lack of consistency among the corresponding literature. To get a fresh, comprehensive view on Bufotes phylogeography, here we combined genome-wide multilocus analyses (RAD-seq) with an extensive compilation of mitochondrial, genome size, niche modelling, distribution and phenotypic (bioacoustics, morphometrics, toxin composition) datasets, representing hundreds of populations throughout Eurasia. We provide a fully resolved nuclear phylogeny for Bufotes and highlight exceptional cyto-nuclear discordances characteristic of complete mtDNA replacement (in 20% of species), mitochondrial surfing during post-glacial expansions, and the formation of homoploid hybrid populations. Moreover, we traced the origin of several allopolyploids down to species level, showing that all were exclusively fathered by the West Himalayan B. latastii but mothered by several diploid forms inhabiting Central Asian lowlands, an asymmetry consistent with hypotheses on mate choice and Dobzhansky-Muller incompatibilities. Their intermediate call phenotypes potentially allowed for rapid reproductive isolation, while toxin compositions converged towards the ecologically-closest parent. Across the radiation, we pinpoint a stepwise progression of reproductive isolation through time, with a threshold below which hybridizability is irrespective of divergence (<6My), above which species barely admix and eventually evolve different mating calls (6-10My), or can successfully cross-breed through allopolyploidization (>15My). Finally, we clarified the taxonomy of Bufotes (including genetic analyses of type series) and formally described two new species, B. cypriensis sp. nov. (endemic to Cyprus) and B. perrini sp. nov. (endemic to Central Asia). Embracing the genomic age, our framework marks the advent of a new exciting era for evolutionary research in these iconic amphibians.

Radically different phylogeographies and patterns of genetic variation in two European brown frogs, genus Rana.

We reconstruct range-wide phylogeographies of two widespread and largely co-occurring Western Palearctic frogs, Rana temporaria and R. dalmatina. Based on tissue or saliva samples of over 1000 individuals, we compare a variety of genetic marker systems, including mitochondrial DNA, single-copy protein-coding nuclear genes, microsatellite loci, and single nucleotide polymorphisms (SNPs) of transcriptomes of both species. The two focal species differ radically in their phylogeographic structure, with R. temporaria being strongly variable among and within populations, and R. dalmatina homogeneous across Europe with a single strongly differentiated population in southern Italy. These differences were observed across the various markers studied, including microsatellites and SNP density, but especially in protein-coding nuclear genes where R. dalmatina had extremely low heterozygosity values across its range, including potential refugial areas. On the contrary, R. temporaria had comparably high range-wide values, including many areas of probable postglacial colonization. A phylogeny of R. temporaria based on various concatenated mtDNA genes revealed that two haplotype clades endemic to Iberia form a paraphyletic group at the base of the cladogram, and all other haplotypes form a monophyletic group, in agreement with an Iberian origin of the species. Demographic analysis suggests that R. temporaria and R. dalmatina have genealogies of roughly the same time to coalescence (TMRCA ~3.5 mya for both species), but R. temporaria might have been characterized by larger ancestral and current effective population sizes than R. dalmatina. The high genetic variation in R. temporaria can therefore be explained by its early range expansion out of Iberia, with subsequent cycles of differentiation in cryptic glacial refugial areas followed by admixture, while the range expansion of R. dalmatina into central Europe is a probably more recent event.

Fossorial but widespread: the phylogeography of the common spadefoot toad (Pelobates fuscus), and the role of the Po Valley as a major source of genetic variability.

Pelobates fuscus is a fossorial amphibian that inhabits much of the European plain areas. To unveil traces of expansion and contraction events of the species' range, we sequenced 702 bp of the mitochondrial cytochrome b gene. To infer the population history we applied phylogeographical methods, such as nested clade phylogeographical analysis (NCPA), and used summary statistics to analyse population structure under a neutral model of evolution. Populations were assigned to different drainage systems and we tested hypotheses of explicit refugial models using information from analysis of molecular variance, nucleotide diversity, effective population size estimation, NCPA, mismatch distribution and Bayesian dating. Coalescent simulations were used as post hoc tests for plausibility of derived or a priori assumed biogeographical hypotheses. Our combination of all approaches enabled the reconstruction of the colonization history and phylogeography of P. fuscus and confirmed a previous assumption of the existence of two major genetic lineages within P. fuscus. Using the Afro-European vicariance of Pelobates cultripes and Pelobates varaldii and applying Bayesian dating we estimated the divergence of these phylogeographical lineages to the Pliocene. We suggest the existence of three different glacial refugia: (i) the area between the Caspian and Black Seas as the origin for the expansion of the 'eastern lineage'; (ii) the Danube system as a centre of diversity for part of the 'western lineage'; (iii) the Po Valley, the largest centre of genetic variability. This fits the hypothesis that climatic fluctuation was a key event for differentiation processes in P. fuscus.

Correlations of geographic distribution and temperature of embryonic development with the nuclear DNA content in the Salamandridae (Urodela, Amphibia).

We used flow cytometry to measure the nuclear DNA content in erythrocytes of 27 salamandrid species. Across these species, diploid genome size varied more than 2 fold (51.3-104.4 pg). According to genome size and geographic distribution, 3 groups of newt species were recognized: West Palearctics with smaller amounts of nuclear DNA; Nearctic, with intermediate values; and East Asiatic, with higher genome sizes. Viviparous West Palearctic salamanders differed from most of the oviparous West Palearctic newts in possessing larger genome sizes. The nuclear DNA content strongly correlates with species range limits. At the same temperature, embryos of salamandrid species with larger genome sizes have a markedly longer developmental time than those with smaller genomes. We present an analysis of the relationships between the amount of nuclear DNA and water temperature at the breeding sites.