Genomic, phenotypic and environmental correlates of speciation in the midwife toads ( Alytes )

Speciation, i.e., the formation of new species, implies that diverging populations evolve genetic and phenotypic factors that promote reproductive isolation (RI), but the adaptive vs. neutral origin of these factors and their relative contributions across the speciation continuum remain elusive. Here we test which of genomic, bioacoustic, morphological and environmental differentiation best predicts RI across the midwife toads (genus Alytes ), a diversification of Mediterranean amphibians. Hybrid zone analyses in the A. obstetricans complex support that without strong geographic barriers to dispersal, the extent of introgression (which should reflect the strength of RI) covaries with genomic divergence irrespective of other factors. Phenotypic divergence become important later along the continuum, namely between non-admixing species attributed to distinct subgenera. Our results suggest that by putatively causing intrinsic incompatibilities in hybrids, the genetic mutations accumulating randomly between allopatric populations act as the initial trigger of RI, while substantial ecological and behavioral differentiation is a long-term consequence of species divergence that ultimately promotes sympatry. Whereas speciation is usually claimed to be primarily adaptive, our study corroborates recent findings that new species may also be a neutral outcome of gradual phylogeographic divergence, which has practical implications for species delimitation.