Integrative taxonomy in a rapid speciation group associated with mating system transition: A case study in the Primula cicutariifolia complex

Abstract

Accurate species delimitation is the key to biodiversity conservation and is fundamental to most branches of biology. However, species delimitation remains challenging in those evolutionary radiations associated with mating system transition from outcrossing to self-fertilization, which have frequently occurred in angiosperms and are usually accompanied by rapid speciation. Here, using the Primula cicutariifolia complex as a case, we integrated molecular, morphological and reproductive isolation evidence to test and verify whether its outcrossing (distylous) and selfing (homostylous) populations have developed into independent evolutionary lineages. Phylogenetic trees based on whole plastomes and SNPs of the nuclear genome both indicated that the distylous and homostylous populations grouped into two different clades. Multispecies coalescent, gene flow and genetic structure analyses all supported such two clades as two different genetic entities. In morphology, as expected changes in selfing syndrome, homostylous populations have significantly fewer umbel layers and smaller flower and leaf sizes compared to distylous populations, and the variation range of some floral traits, such as corolla diameter and umbel layers, show obvious discontinuity. Furthermore, hand-pollinated hybridization between the two clades produced almost no seeds, indicating that well post-pollination reproductive isolation has been established between them. Therefore, the distylous and homostylous populations in this studied complex are two independent evolutionary lineages, and thus these distylous populations should be treated as a distinct species, here named Primula qiandaoensis W. Zhang & J.W. Shao sp. nov.. Our empirical study of the P. cicutariifolia complex highlights the importance of applying multiple lines of evidence, in particular genomic data, to delimit species in pervasive evolutionary plant radiations associated with mating system transition.