Phenotypic Divergence and Genomic Architecture Between Parallel Ecotypes at Two Different Points on the Speciation Continuum in a Marine Snail.

Speciation is rarely observable directly. A way forward is to compare pairs of ecotypes that evolved in parallel in similar contexts but have reached different degrees of reproductive isolation. Such comparisons are possible in the marine snail Littorina saxatilis by contrasting barriers to gene flow between parallel ecotypes in Spain and Sweden. In both countries, divergent ecotypes have evolved to withstand either crab predation or wave action. Here, we explore transects spanning contact zones between the Crab and the Wave ecotypes using low-coverage whole-genome sequencing, morphological and behavioural traits. Despite parallel phenotypic divergence, distinct patterns of differentiation between the ecotypes emerged: a continuous cline in Sweden indicating a weak barrier to gene flow, but two highly genetically and phenotypically divergent, and partly spatially overlapping clusters in Spain suggesting a much stronger barrier to gene flow. The absence of Spanish early-generation hybrids supported strong isolation, but a low level of gene flow is evident from molecular data. In both countries, highly differentiated loci were located in both shared and country-specific chromosomal inversions but were also present in collinear regions. Despite being considered the same species and showing similar levels of phenotypic divergence, the Spanish ecotypes are much closer to full reproductive isolation than the Swedish ones. Barriers to gene flow of very different strengths between ecotypes within the same species might be explained by dissimilarities in the spatial arrangement of habitats, the selection gradients or the ages of the systems.