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. Significance Reproductive isolation is instrumental to the formation of new species (speciation), but it remains largely enigmatic how many incompatibilities are required to prevent hybridization and where they lie across the genome. By studying patterns of admixture in amphibian hybrid zones, we found that reproductive isolation is initiated by numerous small-effect incompatibilities scattered across the genome rather than concentrated in a few important genes. Unlike mammals and birds, in which Y/W degeneracy is a major cause of hybrid dysfunctions, the undifferentiated sex chromosomes of amphibians do not always host more genetic incompatibilities than other chromosomes. These combined results might explain why amphibian speciation is relatively slow, and its clock-like dynamics offer practical perspectives to categorize evolutionary lineages into species or subspecies.

The major climatic oscillations during the Quaternary Ice Age heavily influenced the distribution of species and left their mark on intraspecific genetic diversity. Past range shifts can be reconstructed with the aid of species distribution modeling and phylogeographical analyses. We test the responses of the different members of the genus Triturus (i.e. the marbled and crested newts) as the climate shifted from the previous glacial period (the Last Glacial Maximum, ~21 Ka) to the current interglacial.

Interspecies transfer of mitochondrial (mt) DNA is a common phenomenon in plants, invertebrates and vertebrates, normally linked with hybridization of closely related species in zones of sympatry or parapatry. In central Europe, in an area north of 48°N latitude and between 8° and 22°E longitude, western Palaearctic water frogs show massive unidirectional introgression of mtDNA: 33.7% of 407 Rana ridibunda possessed mtDNA specific for Rana lessonae. By contrast, no R. lessonae with R. ridibunda mtDNA was observed. That R. ridibunda with introgressed mitochondrial genomes were found exclusively within the range of the hybrid Rana esculenta and that most hybrids had lessonae mtDNA (90.4% of 335 individuals investigated) is evidence that R. esculenta serves as a vehicle for transfer of lessonae mtDNA into R. ridibunda. Such introgression has occurred several times independently. The abundance and wide distribution of individuals with introgressed mitochondrial genomes show that R. lessonae mt genomes work successfully in a R. ridibunda chromosomal background despite their high sequence divergence from R. ridibunda mtDNAs (14.2–15.2% in the ND2/ND3 genes). Greater effectiveness of enzymes encoded by R. lessonae mtDNA may be advantageous to individuals of R. ridibunda and probably R. esculenta in the northern parts of their ranges.