Abstract: The key process in speciation concerns the formation and maintenance of reproductive isolating barriers between diverging lineages. Although species boundaries are frequently investigated between two species across many taxa, reproductive isolating barriers among multiple species (>2) that would represent the most common phenomenon in nature, remain to be clarified. Here, we use double digest restriction-site associated DNA (ddRAD) sequencing to examine patterns of hybridization at a sympatric site where three Ligularia species grow together and verify whether those patterns contribute to the maintenance of boundaries among species. The results based on the RAD SNP datasets indicated hybridization Ligularia cyathiceps × L. duciformis and L. duciformis × L. yunnanensis were both restricted to F1s plus a few first-generation backcrosses and no gene introgression were identified, giving rise to strong reproductive isolation among hybridizing species. Moreover, hybrid swarm simulation, using HYBRIDLAB, indicated the RAD SNP datasets had sufficient discriminatory power for accurate hybrid detection. We conclude that parental species show strong reproductive isolation and they still maintain species boundaries, which may be the key mechanism to maintain species diversity of Ligularia in the eastern Qinghai-Tibetan Plateau and adjacent areas. Moreover, this study highlights the effectiveness of RAD sequencing in hybridization studies.