Patterns of Genomic Divergence and Introgression in Two <i>Primulina</i> Hybrid Zones

Hybrid zones have long been promoted as natural laboratories for understanding the mechanisms of speciation. Multiple or replicated hybrid zones are particularly informative, as they allow for assessing the consistency of genomic divergence and introgression across different environmental contexts and demographic histories, thereby improving our understanding of the factors that drive or hinder speciation on a broader scale. Here, using whole-genome resequencing data, we compare the patterns of genomic divergence and introgression in two Primulina hybrid zones. We found that genomic divergence in both hybrid zones is largely shaped by neutral processes, with only a few genomic regions showing signatures of balancing or lineage-specific selection. Genomic cline analyses identified numerous SNPs that showed significantly steeper clines and biased centres than the genome-wide expectation in both hybrid zones, consistent with the existence of reproductive barriers. Within regions of restricted gene flow, we identified 21 genes shared between the two hybrid zones. Annotation of gene function revealed that several genes are involved in reproductive processes. In addition, many zone-specific outlier loci were linked to genes associated with pollen and flower development, suggesting that these barriers may contribute to reproductive isolation under localised ecological conditions. Overall, these findings suggest that while certain reproductive barriers remain consistent across independent hybrid zones, others may be contingent on local environmental contexts. Our results demonstrate that both general and zone-specific mechanisms contribute to reproductive isolation in Primulina, providing empirical evidence that some genomic barriers recur across independent hybrid zones while others arise through localised adaptation.