No abstract

In organisms reproducing sexually, speciation occurs when increasing divergence results in pre-or post-zygotic reproductive isolation between lineages. Studies focusing on reproductive isolation origin in early stages of speciation are common and many rely on genomic scans to infer introgression providing limited information on the genomic architecture of reproductive isolation long-term maintenance. This study analyses a natural hybrid zone between two species in a late stage of speciation. We used ddRADseq genotyping in the contact between Podarcis bocagei and P. carbonelli to examine admixture extent, analyse hybrid zone stability and assess genome-wide variation in selection against introgression. We confirmed strong but incomplete reproductive isolation in a bimodal hybrid zone. New findings revealed population genetic structure within P. carbonelli in the contact zone; geographical and genomic clines analysis suggested strong selection against gene flow, but a relatively small proportion of the loci can introgress, mostly within the narrow contact zone. However, geographical clines revealed that a few introgressed loci show signs of potential positive selection, particularly into P. bocagei. Geographical clines also detected a signal of hybrid zone movement towards P. bocagei distribution. Genomic cline analysis revealed heterogeneous patterns of introgression among loci within the syntopy zone, but the majority maintain a strong association with the genomic background of origin. However, incongruences between both cline approaches were found, potentially driven by confounding effects on genomic clines. Last, an important role of the Z chromosome in reproductive isolation is suggested. Importantly, overall patterns of restricted introgression seem to result from numerous strong intrinsic barriers across the genome.

In organisms reproducing sexually, speciation occurs when increasing divergence results in pre- or post-zygotic reproductive isolation between lineages. Studies focusing on reproductive isolation origin in early stages of speciation are common. Many rely on indirect measures of introgression providing limited information on the genomic architecture of reproductive isolation maintenance in the long term. This study focuses on direct measures of introgression between two species in a late stage of speciation. We used ddRADseq genotyping in a natural hybrid zone between Podarcis bocagei and P. carbonelli to examine admixture extent, estimate effective selection, analyse hybrid zone stability, and assess variation in selection against introgression across the genome. Hybridization was confirmed the narrow and bimodal hybrid zone demonstrating the existence of strong mechanisms of reproductive isolation. Simulations suggested that simple premating barriers were not enough to explain the observed distribution of admixture classes, pointing out the role of post-mating isolation. A geographic cline approach confirmed strong reproductive isolation and high effective selection preventing extensive introgression outside of the contact zone; and detected a signal of hybrid zone movement towards P. bocagei distribution. Genomic cline analysis revealed heterogeneous patterns of introgression among loci within the syntopy zone, but most of the loci do not introgress more or less than the genomic average maintaining a strong association with the genomic background of origin. However, genomic clines can be driven by confounding effects resulting in incongruences between both cline approaches. Importantly, overall patterns of introgression seem to result from strong intrinsic barriers across the genome, without clear “islands of differentiation”. Last, an important role of the Z chromosome in reproductive isolation is suggested.