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.

Mitochondrial phylogeny shows multiple independent ecological transitions and northern dispersion despite of Pleistocene glaciations in meadow and steppe vipers (Vipera ursinii and Vipera renardi).

The phylogeny and historical demography of small Eurasian vipers of the Vipera ursinii and V. renardi complexes were studied using mitochondrial DNA sequences analysed with Bayesian inference, Maximum Likelihood and Maximum Parsimony approaches, and mismatch distributions. Diversification in the group resulted from an initial dispersion in the later Pliocene - Pleistocene in two directions: north-westwards via the Balkans (V. ursinii complex) and north-eastwards from Asia Minor via the Caucasus (V. renardi complex). An independent, comparatively recent transition occurred from montane habitats to lowland grasslands in different mitochondrial lineages during the Late Pleistocene, when representatives of the both complexes had reached lowland steppes to the north. Effective population size showed clear signs of rapid growth in eastern V. renardi, triggered by colonization of vast lowland steppes, but in western V. ursinii complex grew during the Last Glaciation and experienced stabilization in Holocene. Expansion and population growth in lowland lineages of V. renardi was not strongly affected by Pleistocene climatic oscillations, when cold, dry conditions could have favoured species living in open grasslands. The high diversity of closely related haplotypes in the Caucasus and Tien-Shan could have resulted from repetitive expansion-constriction-isolation events in montane regions during Pleistocene climate fluctuations. The mitochondrial phylogeny pattern conflicts with the current taxonomy.