Inferring the degree of incipient speciation in secondary contact zones of closely related lineages of Palearctic green toads (Bufo viridis subgroup)

Reproductive isolation between lineages is expected to accumulate with divergence time, but the time taken to speciate may strongly vary between different groups of organisms. In anuran amphibians, laboratory crosses can still produce viable hybrid offspring >20 My after separation, but the speed of speciation in closely related anuran lineages under natural conditions is poorly studied. Palearctic green toads (Bufo viridis subgroup) offer an excellent system to address this question, comprising several lineages that arose at different times and form secondary contact zones. Using mitochondrial and nuclear markers, we previously demonstrated that in Sicily, B. siculus and B. balearicus developed advanced reproductive isolation after Plio-Pleistocene divergence (2.6 My, 3.3–1.9), with limited historic mtDNA introgression, scarce nuclear admixture, but low, if any, current gene flow. Here, we study genetic interactions between younger lineages of early Pleistocene divergence (1.9 My, 2...