Abstract: Natural hybridization is a crucial evolutionary process and a long-standing topic of study in evolutionary biology. Hybrid zones, where two congeneric species interact, can provide insight into the process of natural hybridization, especially with respect to how taxon diversity is maintained. In this study, we used double digest restriction-site associated DNA sequencing technology (ddRAD-seq) to examine genetic structure and estimate introgression in four hybrid zones of Ligularia tongolensis and Ligularia cymbulifera. Our analysis demonstrated that parental species were highly differentiated, whereas pairwise F_ST between parents and their hybrids was low, indicating that sympatric sites can form hybrid zones. As most F₁ hybrid individuals were observed within these zones, our finding also implied the presence of substantial barriers to interbreeding. Furthermore, some individuals that possessed the typical morphology of the parental species belonged to the F₁ generation. Genomic clines analysis revealed that a large fraction of single nucleotide polymorphisms (SNPs) deviated from a model of neutral introgression in the four hybrid zones, and most SNPs exhibited selection favoring the L. cymbulifera genotype. Bidirectional but asymmetric introgression was revealed as evident in the four hybrid zones. Habitat differences between the four hybrid zones may affect isolation barriers between both species. Taken together, these findings suggest that where incomplete reproductive barriers allow natural hybridization, the introgression between species generates rich genetic recombination that contributes to the fast adaptation and diversification of the widespread Ligularia in the Hengduan Mountains Region (HMR).

Key words: ddRAD, hybrid zones, Ligularia, natural hybridization, reproductive isolation