Hybridization outcomes reflect context-dependent reproductive isolation in two damselfly hybrid zones

Pair of species with multiple hybrid zones provide powerful systems to examine how reproductive isolation and reinforcement may depend on ecological and demographic contexts. We conducted a comparative analysis of two hybrid zones between Ischnura elegans and Ischnura graellsii in North-central (NC) and North-western (NW) Spain to evaluate how differences in reproductive barriers shape patterns of hybridization and admixture. First, we quantified five prezygotic barriers in the NC hybrid zone--where no barrier estimates previously existed--and compared them with published data from the NW hybrid zone. Second, we conducted genomic analyses using 5,702 SNPs to assess genetic structure, differentiation, and diversity between species; and 381 SNPs with species-specific alleles to characterize hybrid classes and introgression across both regions. Reproductive isolation was strongly context dependent. Mechanical reinforcement acted in the same cross direction in both hybrid zones, but the NC hybrid zone also showed stronger gametic (fertility) isolation. Consequently, total prezygotic isolation was higher in the NC hybrid zone. These differences were also observed genomically: the NC hybrid zone exhibited limited and largely unidirectional introgression, whereas the NW hybrid zone showed more extensive admixture, bidirectional gene flow, and reduced interspecific differentiation, consistent with a longer or more sustained history of hybridization. Our results demonstrate that hybridization outcomes can differ markedly between hybrid zones involving the same pair of species, underscoring the importance of evaluating spatial variation in reproductive barriers and introgression to assess the presence, asymmetry, and repeatability of reinforcement, as well as the microevolutionary processes shaping hybrid zone structure.