The repeated transition from outcrossing to selfing is a key topic in evolutionary biology. However, the molecular basis of such shifts has been rarely examined due to lack of knowledge of the genes controlling these transitions. A classic example of mating system transition is the repeated shift from heterostyly to homostyly. Occurring in 28 angiosperm families, heterostyly is characterized by the reciprocal position of male and female sexual organs in two (or three) distinct, usually self-incompatible floral morphs. Conversely, homostyly is characterized by a single, self-compatible floral morph with reduced separation of male and female organs, facilitating selfing. Here, we investigate the origins of homostyly in Primula vulgaris and its microevolutionary consequences by integrating surveys of the frequency of homostyles in natural populations, DNA sequence analyses of the gene controlling the position of female sexual organs (CYPᵀ), and microsatellite genotyping of both progeny arrays and natural populations characterized by varying frequencies of homostyles. As expected, we found that homostyles displace short-styled individuals, but long-style morphs are maintained at low frequencies within populations. We also demonstrated that homostyles repeatedly evolved from short-styled individuals in association with different types of loss-of-function mutations in CYPᵀ. Additionally, homostyly triggers a shift to selfing, promoting increased inbreeding within and genetic differentiation among populations. Our results elucidate the causes and consequences of repeated transitions to homostyly within species, and the putative mechanisms precluding its fixation in P. vulgaris. This study represents a benchmark for future analyses of losses of heterostyly in other angiosperms.

Summary The complex nature of species boundaries has been a central topic in evolutionary biology ever since Darwin. Despite numerous separate studies on reproductive isolation and hybridization, their relationship remains underinvestigated. Are the strengths and asymmetries of reproductive barriers reflected in the extent and directionalities of interspecific genetic exchange? We combined field, experimental, and molecular data to quantify strengths and asymmetries of sympatric reproductive barriers and hybridization between florally heteromorphic primroses. We also assessed whether generalist pollinators discriminate between different floral cues and contribute to reproductive isolation, a long‐debated topic. Sympatric reproductive isolation is high but incomplete, and most phenotypic intermediates are genetic F1 hybrids, whereas backcrosses are rare, revealing low interspecific gene flow. Species integrity rests on multiple barriers, but ethological isolation is among the strongest, demonstrating that even generalist pollinators crucially contribute to the maintenance of species boundaries. Furthermore, reproductive barriers are weaker for Primula veris and short‐styled plants, results corroborated by molecular data. Thus, in florally heteromorphic systems, both species‐ and morph‐dependent asymmetries affect permeability of species boundaries. Our study illustrates how the interactions between complex floral syndromes and pollinators shape species boundaries in unique, previously undescribed ways.

Summary Biodiversity hotspots, such as the Caucasus mountains, provide unprecedented opportunities for understanding the evolutionary processes that shape species diversity and richness. Therefore, we investigated the evolution of Primula sect. Primula, a clade with a high degree of endemism in the Caucasus. We performed phylogenetic and network analyses of whole‐genome resequencing data from the entire nuclear genome, the entire chloroplast genome, and the entire heterostyly supergene. The different characteristics of the genomic partitions and the resulting phylogenetic incongruences enabled us to disentangle evolutionary histories resulting from tokogenetic vs cladogenetic processes. We provide the first phylogeny inferred from the heterostyly supergene that includes all species of Primula sect. Primula. Our results identified recurrent admixture at deep nodes between lineages in the Caucasus as the cause of non‐monophyly in Primula. Biogeographic analyses support the ‘out‐of‐the‐Caucasus’ hypothesis, emphasizing the importance of this hotspot as a cradle for biodiversity. Our findings provide novel insights into causal processes of phylogenetic discordance, demonstrating that genome‐wide analyses from partitions with contrasting genetic characteristics and broad geographic sampling are crucial for disentangling the diversification of species‐rich clades in biodiversity hotspots.

Monomorphic enantiostyly should be resistant to invasion of the dimorphic form because the prevalent ecological conditions favour the maintenance of geitonogamous individuals that are able to take advantage of ecological heterogeneity and generalized pollination limitation.