The Primula edelbergii S-locus is an example of a jumping supergene.

Research on supergenes, non-recombining genomic regions housing tightly linked genes that control complex phenotypes, has recently gained prominence in genomics. Heterostyly, a floral heteromorphism promoting outcrossing in several angiosperm families, is controlled by the S-locus supergene. The S-locus has been studied primarily in closely related Primula species and, more recently, in other groups that independently evolved heterostyly. However, it remains unknown whether genetic architecture and composition of the S-locus are maintained among species that share a common origin of heterostyly and subsequently diverged across larger time scales. To address this research gap, we present a chromosome-scale genome assembly of Primula edelbergii, a species that shares the same origin of heterostyly with Primula veris (whose S-locus has been characterized) but diverged from it 18 million years ago. Comparative genomic analyses between these two species allowed us to show, for the first time, that the S-locus can 'jump' (i.e. translocate) between chromosomes maintaining its function in controlling heterostyly. Additionally, we found that four S-locus genes were conserved but reshuffled within the supergene, seemingly without affecting their expression, thus we could not detect changes explaining the lack of self-incompatibility in P. edelbergii. Furthermore, we confirmed that the S-locus is not undergoing genetic degeneration. Finally, we investigated P. edelbergii evolutionary history within Ericales in terms of whole genome duplications and transposable element accumulation. In summary, our work provides a valuable resource for comparative analyses aimed at investigating the genetics of heterostyly and the pivotal role of supergenes in shaping the evolution of complex phenotypes.

Survival in nunatak and peripheral glacial refugia of three alpine plant species is partly predicted by altitudinal segregation.

Mountain biota survived the Quaternary cold stages most probably in peripheral refugia and/or ice-free peaks within ice-sheets (nunataks). While survival in peripheral refugia has been broadly demonstrated, evidence for nunatak refugia is still scarce. We generated RADseq data from three mountain plant species occurring at different elevations in the southeastern European Alps to investigate the role of different glacial refugia during the Last Glacial Maximum (LGM). We tested the following hypotheses. (i) The deep Piave Valley forms the deepest genetic split in the species distributed across it, delimiting two peripheral refugia. (ii) The montane to alpine species Campanula morettiana and Primula tyrolensis survived the LGM in peripheral refugia, while high-alpine to subnival Saxifraga facchinii likely survived in several nunatak refugia. (iii) The lower elevation species suffered a strong population decline during the LGM. By contrast, the higher elevation species shows long-term stability of population sizes due to survival on permanently ice-free peaks and small population sizes at present. We found peripheral refugia on both sides of the Piave Valley, which acted as a major genetic barrier. Demographic modelling confirmed nunatak survival not only for S. facchinii but also for montane to alpine C. morettiana. Altitudinal segregation influenced the species' demographic fluctuations, with the lower elevation species showing a significant population increase at the end of the LGM, and the higher elevation species either showing decrease towards the present or stable population sizes with a short bottleneck. Our results highlight the role of nunatak survival and species ecology in the demographic history of mountain species.

Development of polymorphic microsatellite markers for distylous-homostylous Primula secundiflora (Primulaceae) using HiSeq sequencing.

Primula secundiflora is an insect-pollinated, perennial herb belonging to the section Proliferae (Primulaceae) that exhibits considerable variation in its mating system, with predominantly outcrossing populations comprising long-styled and short-styled floral morphs and selfing populations comprising only homostyles. To facilitate future investigations of the population genetics and mating patterns of this species, we developed 25 microsatellite markers from P. secundiflora using next-generation sequencing and measured polymorphism and genetic diversity in a sample of 30 individuals from three natural populations. The markers displayed high polymorphism, with the number of observed alleles per locus ranging from three to 16 (mean = 8.36). The observed and expected heterozygosities ranged from 0.100 to 1.000 and 0.145 to 0.843, respectively. Twenty-one of the loci were also successfully amplified in P. denticulata. These microsatellite markers should provide powerful tools for investigating patterns of population genetic diversity and the evolutionary relationships between distyly and homostyly in P. secundiflora.

An integrated primary care service to reduce cardiovascular disease risk in people with severe mental illness: Primrose-A - thematic analysis of its acceptability, feasibility, and implementation.

BACKGROUND: Cardiovascular disease among patients with severe mental illness in England is a major preventable contributor to premature mortality. To address this, a nurse and peer-coach delivered service (Primrose-A) was implemented in three London general practices from 2019 (implementation continued during COVID-19). This study aimed to conduct interviews with patient and staff to determine the acceptability of, and experiences with, Primrose-A. METHODS: Semi-structured audio-recorded interviews with eight patients who had received Primrose-A, and 3 nurses, 1 GP, and 1 peer-coach who had delivered Primrose-A in three London-based GP surgeries were conducted. Reflexive thematic analysis was used to identify themes from the transcribed interviews. FINDINGS: Overall, Primrose-A was viewed positively by patients and staff, with participants describing success in improving patients' mental health, isolation, motivation, and physical health. Therapeutic relationships between staff and patients, and long regular appointments were important facilitators of patient engagement and acceptance of the intervention. Several barriers to the implementation of Primrose-A were identified, including training, administrative and communication issues, burden of time and resources, and COVID-19. CONCLUSIONS: Intervention acceptability could be enhanced by providing longer-term continuity of care paired with more peer-coaching sessions to build positive relationships and facilitate sustained health behaviour change. Future implementation of Primrose-A or similar interventions should consider: (1) training sufficiency (covering physical and mental health, including addiction), (2) adequate staffing to deliver the intervention, (3) facilitation of clear communication pathways between staff, and (4) supporting administrative processes.

Parallel evolution of morphological and genomic selfing syndromes accompany the breakdown of heterostyly.

Evolutionary transitions from outcrossing to selfing in flowering plants have convergent morphological and genomic signatures and can involve parallel evolution within related lineages. Adaptive evolution of morphological traits is often assumed to evolve faster than nonadaptive features of the genomic selfing syndrome. We investigated phenotypic and genomic changes associated with transitions from distyly to homostyly in the Primula oreodoxa complex. We determined whether the transition to selfing occurred more than once and investigated stages in the evolution of morphological and genomic selfing syndromes using 22 floral traits and both nuclear and plastid genomic data from 25 populations. Two independent transitions were detected representing an earlier and a more recently derived selfing lineage. The older lineage exhibited classic features of the morphological and genomic selfing syndrome. Although features of both selfing syndromes were less developed in the younger selfing lineage, they exhibited parallel development with the older selfing lineage. This finding contrasts with the prediction that some genomic changes should lag behind adaptive changes to morphological traits. Our findings highlight the value of comparative studies on the timing and extent of transitions from outcrossing to selfing between related lineages for investigating the tempo of morphological and molecular evolution.

PfPIN5 promotes style elongation by regulating cell length in Primula forbesii Franch.

BACKGROUND AND AIMS: Style dimorphism is one of the polymorphic characteristics of flowers in heterostylous plants, which have two types of flowers: the pin morph, with long styles and shorter anthers, and the thrum morph, with short styles and longer anthers. The formation of dimorphic styles has received attention in the plant world. Previous studies showed that CYP734A50 in Primula determined style length and limited style elongation and that the brassinosteroid metabolic pathway was involved in regulation of style length. However, it is unknown whether there are other factors affecting the style length of Primula. METHODS: Differentially expressed genes highly expressed in pin morph styles were screened based on Primula forbesii transcriptome data. Virus-induced gene silencing was used to silence these genes, and the style length and anatomical changes were observed 20 days after injection. KEY RESULTS: PfPIN5 was highly expressed in pin morph styles. When PfPIN5 was silenced, the style length was shortened in pin and long-homostyle plants by shortening the length of style cells. Moreover, silencing CYP734A50 in thrum morph plants increased the expression level of PfPIN5 significantly, and the style length increased. The results indicated that PfPIN5, an auxin efflux transporter gene, contributed to regulation of style elongation in P. forbesii. CONCLUSIONS: The results implied that the auxin pathway might also be involved in the formation of styles of P. forbesii, providing a new pathway for elucidating the molecular mechanism of style elongation in P. forbesii.

The genomes of Darwin's primroses reveal chromosome-scale adaptive introgression and differential permeability of species boundaries.

Introgression is an important source of genetic variation that can determine species adaptation to environmental conditions. Yet, definitive evidence of the genomic and adaptive implications of introgression in nature remains scarce. The widespread hybrid zones of Darwin's primroses (Primula elatior, Primula veris, and Primula vulgaris) provide a unique natural laboratory for studying introgression in flowering plants and the varying permeability of species boundaries. Through analysis of 650 genomes, we provide evidence of an introgressed genomic region likely to confer adaptive advantage in conditions of soil toxicity. We also document unequivocal evidence of chloroplast introgression, an important precursor to species-wide chloroplast capture. Finally, we provide the first evidence that the S-locus supergene, which controls heterostyly in primroses, does not introgress in this clade. Our results contribute novel insights into the adaptive role of introgression and demonstrate the importance of extensive genomic and geographical sampling for illuminating the complex nature of species boundaries.

Haplotype-resolved genome assembly provides insights into the evolution of S-locus supergene in distylous Nymphoides indica.

Distyly has evolved independently in numerous animal-pollinated angiosperm lineages. Understanding of its molecular basis has been restricted to a few species, primarily Primula. Here, we investigate the genetic architecture of the single diallelic locus (S-locus) supergene, a linkage group of functionally associated genes, and explore how it may have evolved in distylous Nymphoides indica, a lineage of flowering plants not previously investigated. We assembled haplotype-resolved genomes, used read-coverage-based genome-wide association study (rb-GWAS) to locate the S-locus supergene, co-expression network analysis to explore gene networks underpinning the development of distyly, and comparative genomic analyses to investigate the origins of the S-locus supergene. We identified three linked candidate S-locus genes - NinBAS1, NinKHZ2, and NinS1 - that were only evident in the short-styled morph and were hemizygous. Co-expression network analysis suggested that brassinosteroids contribute to dimorphic sex organs in the short-styled morph. Comparative genomic analyses indicated that the S-locus supergene likely evolved via stepwise duplications and has been affected by transposable element activities. Our study provides novel insight into the structure, regulation, and evolution of the supergene governing distyly in N. indica. It also provides high-quality genomic resources for future research on the molecular mechanisms underlying the striking evolutionary convergence in form and function across heterostylous taxa.

Transcriptome Analyses Reveal the Aroma Terpeniods Biosynthesis Pathways of Primula forbesii Franch. and the Functional Characterization of the PfDXS2 Gene.

Primula forbesii Franch. is a unique biennial herb with a strong floral fragrance, making it an excellent material for studying the aroma characteristics of the genus Primula. The floral scent is an important ornamental trait that facilitates fertilization. However, the molecular mechanism regulating the floral scent in Primula is unknown. In order to better understand the biological mechanisms of floral scents in this species, this study used RNA sequencing analysis to discuss the first transcriptome sequence of four flowering stages of P. forbesii, which generated 12 P. forbesii cDNA libraries with 79.64 Gb of clean data that formed 51,849 unigenes. Moreover, 53.26% of the unigenes were annotated using public databases. P. forbesii contained 44 candidate genes covering all known enzymatic steps for the biosynthesis of volatile terpenes, the major contributor to the flower's scent. Finally, 1-deoxy-d-xylulose 5-phosphate synthase gene of P. forbesii (PfDXS2, MK370094), the first key enzyme gene in the 2-c-methyl-d-erythritol 4-phosphate (MEP) pathway of terpenoids, was cloned and functionally verified using virus-induced gene silencing (VIGs). The results showed that PfDXS2-silencing significantly reduced the relative concentrations of main volatile terpenes. This report is the first to present molecular data related to aroma metabolites biosynthesis pathways and the functional characterization of any P. forbesii gene. The data on RNA sequencing provide comprehensive information for further analysis of other plants of the genus Primula.

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

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.

Non-deep physiological dormancy and germination characteristics of Primula florindae (Primulaceae), a rare alpine plant in the Hengduan Mountains of southwest China.

Timing of seed germination is directly related to the survival probability of seedlings. For alpine plants, autumn-dispersal seeds should not germinate immediately because the cold temperature is not conducive to the survival of seedlings. Seed dormancy is a characteristic of the seed that prevents it from germinating after dispersal. Primula florindae is an alpine perennial forb endemic to eastern Tibet, SW China. We hypothesized that primary dormancy and environmental factors prevent seeds of P. florindae to germinate in autumn and allow them to germinate at the first opportunity in spring. We determined how GA3, light, temperature, dry after-ripening (DAR) and cold-wet stratification (CS) treatments affect seed germination by conducting a series of laboratory experiments. Firstly, the effects of gibberellic acid (GA3; 0, 20, and 200 mg L-1) on germination of freshly shed seeds at alternating temperatures (15/5 and 25/15 °C) were immediately investigated to characterize seed with a physiological dormancy component. Then, the fresh seeds treated with 0, 3, and 6 months of after-ripening (DAR) and cold-wet stratification (CS) were incubated at seven constant (1, 5, 10, 15, 20, 25, and 30 °C) and two alternating temperatures (5/1, 15/5, and 25/15 °C) at light and dark conditions. Fresh seeds were dormant, which only germinated well (>60%) at 20, 25, and 25/15 °C in light but not at ≤15 °C and to higher percentages in light than in dark. GA3 increased germination percentage of fresh seeds, and DAR or CS treatments increased final germination percentage, germination rate (speed), and widened the temperature range for germination from high to low. Moreover, CS treatments reduced the light requirement for germination. Thus, after dormancy release, seeds germinated over a wide range of constant and alternating temperatures, regardless of light conditions. Our results demonstrated that P. florindae seeds have type 2 non-deep physiological dormancy. Timing of germination should be restricted to early spring, ensuring a sufficient length of the growing season for seedling recruitment. These dormancy/germination characteristics prevent seeds from germinating in autumn when temperatures are low but allow them to germinate after snowmelt in spring.

Bisbibenzyls from Serbian Primula veris subsp. Columnae (Ten.) Lȕdi and P. acaulis (L.) L.

Bisbibenzyls are specialized metabolites found exclusively in liverworts, until recently; they represent chemical markers of liverworts. Their occurrence in vascular plants was noticed in 2007, when they were found in Primula veris subsp. macrocalyx from Russia. This report prompted us to chemically analyze the two most common Serbian Primula species, P. veris subsp. columnae and P. acaulis, in order to determine the presence of bisbibenzyls in them. Our study revealed nine structurally distinct bisbibenzyls (1-9), identified based on 1D and 2D NMR, IR, UV and HRESIMS data. Among them were five previously undescribed compounds (2-6). The remaining compounds found and previously described in the literature were: the bisbibenzyls riccardin C (1), isoperrottetin A (7), isoplagiochin E (8) and 11-O-demethylmarchantin I (9), as well as 4-hydroxyphenylmethylketone (10) and 4-hydroxy-3-methoxyphenylmethylketone (11). Riccardin C was the most dominant bisbibenzyl in both species studied. Previously, it was the first bisbibenzyl found in vascular plants (P. veris subsp. macrocalyx). An assessment of the cytotoxic activity of the isolated compounds against A549 lung cancer and healthy MRC5 cell lines was also the subject of our study. Compounds 6 and 9 exhibited significant cytotoxic activity expressed by IC50 values of 12 µM, but the selectivity was not satisfactory.

Botanical Briefs: Primula obconica Dermatitis.

Primula obconica, a household plant originally found in China that was introduced in Europe in the 1880s, has been reported to cause plant-induced contact dermatitis (CD). The condition more commonly is reported in Europe and less frequently in the United States, where the plant is not commonly included in patch testing protocols. Clinical features of P obconica CD can include facial and hand as well as fingertip dermatitis. The main allergens known to cause these findings are primin and miconidin. Treatment of P obconica CD mainly involves avoiding contact with the plant and applying a topical steroid.

Skewed morph ratios lead to lower genetic diversity of the heterostylous Primula veris in fragmented grasslands.

Populations of heterostylous plant species are ideally composed of equal frequencies of two (distylous) or three (tristylous) morphologically different floral morphs. Intra-morph incompatibility helps to avoid inbreeding and to maintain genetic diversity, supporting plant fitness and long-term viability. Habitat fragmentation can lead to skewed morph ratios and thereby reduce the abundance of compatible mates. This, in turn, can result in a loss of genetic diversity. We tested whether the genetic diversity of heterostylous plants is affected by morph ratio bias using populations of the distylous grassland plant Primula veris in recently fragmented grasslands. We recorded morph frequencies and population sizes in 30 study populations of P. veris on two Estonian islands characterised by different degrees of habitat fragmentation. Examining variation of thousands of single nucleotide polymorphisms (SNPs) and heterostyly-specific genetic markers, we quantified overall and morph-specific genetic diversity and differentiation in these populations. Morph frequencies deviated more in smaller populations. Skewed morph ratios had a negative effect on the genetic diversity of P. veris in more fragmented grasslands. In the populations of better-connected grassland systems, genetic differentiation among S-morphs was higher than among L-morphs. Our study shows that deviations from morph balance are stronger in small populations and have a negative impact on the genetic diversity of the distylous plant P. veris. Together with the direct negative effects of habitat loss and decreased population size on the genetic diversity of plants, morph ratio bias may intensify the process of genetic erosion, thus exacerbating the local extinction of heterostylous species.

Differential physiological responses and tolerance to potentially toxic elements in Primula forbesii Franch.

Environmental pollution caused by potentially toxic elements (PTEs) has become a global problem that endangers environmental sustainability due to industrial, agricultural, and urban pollution. Primula forbesii Franch. (a synonym of Primula filipes G. Watt.) is a biennial flower native to China with excellent stress resistance and ornamental value. In this study, we examined the phenotypic traits, growth indexes, and physiological properties of P. forbesii in response to five representative PTEs (Cd, Ni, Cr(III), Cu, and Zn) under hydroponic culture conditions. High concentrations of Zn and Cr had little effect on the growth and physiological properties of P. forbesii, indicating that the species has strong tolerance to Zn and Cr stress. Alternatively, high concentrations of Cd, Ni, and Cu seriously affected plant growth and development, resulting in leaf chlorosis and even death, and therefore may have a serious negative impact on the growth of P. forbesii. However, activity levels of some antioxidant enzymes and osmotic regulatory substances remained high, indicating that P. forbesii resisted PTE stress by regulating physiological and biochemical metabolism to a certain extent. Furthermore, principal component analysis and membership function were used to comprehensively evaluate P. forbesii resistance to PTEs. These analyses revealed that P. forbesii exhibits distinct sensitivities and physiological responses to different PTEs and suggested that the resistance to five PTEs in decreasing order is Zn > Cr > Cd > Cu > Ni. These results provide a theoretical basis for the future application of P. forbesii in environments with PTE pollution and may expand its practical utilization.

Inheritance of distyly and homostyly in self-incompatible Primula forbesii.

The evolutionary transition from self-incompatible distyly to self-compatible homostyly frequently occurs in heterostylous taxa. Although the inheritance of distyly and homostyly has been deeply studied, our understanding on modifications of the classical simple Mendelian model is still lacking. Primula forbesii, a biennial herb native to southwest China, is a typical distylous species, but after about 20 years of cultivation with open pollination, self-compatible homostyly appeared, providing ideal material for the study of the inheritance of distyly and homostyly. In this study, exogenous homobrassinolide was used to break the heteromorphic incompatibility of P. forbesii. Furthermore, we performed artificial pollination and open-pollination experiments to observe the distribution of floral morphs in progeny produced by different crosses. The viability of seeds from self-pollination was always the lowest among all crosses, and the homozygous S-morph plants (S/S) occurred in artificial pollination experiments but may experience viability selection. The distyly of P. forbesii is governed by a single S-locus, with S-morph dominant hemizygotes (S/-) and L-morph recessive homozygotes (-/-). Homostylous plants have a genotype similar to L-morph plants, and homostyly may be caused by one or more unlinked modifier genes outside the S-locus. Open pollinations confirm that autonomous self-pollination occurs frequently in L-morphs and homostylous plants. This study deepens the understanding of the inheritance of distyly and details a case of homostyly that likely originated from one or more modifier genes.

Whole-genome analyses disentangle reticulate evolution of primroses in a biodiversity hotspot.

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.

Primrose homostyles: A classic case of possible balancing selection revisited.

In this issue of Molecular Ecology, Mora-Carrera et al. (2022) revisit a case of the loss of an outcrossing system in primroses, which has been studied as an example of balancing selection in the wild since the 1940s. Molecular variants in the gene involved in the mutant self-fertile phenotype, which is now known, help towards understanding this textbook example of breakdown of an outcrossing system. However, as often happens, new information also raises further questions.

Soil water and nutrient availability interactively modify pollinator-mediated directional and correlational selection on floral display.

The individual and combined effects of abiotic factors on pollinator-mediated selection on floral traits are not well documented. To examine potential interactive effects of water and nutrient availability on pollinator-mediated selection on three floral display traits of Primula tibetica, we manipulated pollination and nutrient availability in a factorial experiment, conducted at two common garden sites with different soil water content (natural vs addition). We found that both water and nutrient availability affected floral trait expression in P. tibetica and that hand pollination increased seed production most when both nutrient content and water content were high, indicating joint pollen and resource limitation. We documented selection on all floral traits, and pollinators contributed to directional and correlational selection on plant height and number of flowers. Soil water and nutrient availability interactively influenced the strength of both pollinator-mediated directional and correlational selection, with significant selection observed when nutrient or water availability was high, but not when none or both were added. The results suggest that resource limitation constrains the response of P. tibetica to among-individual variation in pollen receipt, that addition of nutrients or water leads to pollinator-mediated selection and that effects of the two abiotic factors are nonadditive.

Genomic evidence supports the genetic convergence of a supergene controlling the distylous floral syndrome.

Heterostyly, a plant sexual polymorphism controlled by the S-locus supergene, has evolved numerous times among angiosperm lineages and represents a classic example of convergent evolution in form and function. Determining whether underlying molecular convergence occurs could provide insights on constraints to floral evolution. Here, we investigated S-locus genes in distylous Gelsemium (Gelsemiaceae) to determine whether there is evidence of molecular convergence with unrelated distylous species. We used several approaches, including anatomical measurements of sex-organ development and transcriptome and whole-genome sequencing, to identify components of the S-locus supergene. We also performed evolutionary analysis with candidate S-locus genes and compared them with those reported in Primula and Turnera. The candidate S-locus supergene of Gelsemium contained four genes, of which three appear to have originated from gene duplication events within Gelsemiaceae. The style-length genes GeCYP in Gelsemium and CYP734A50 in Primula likely arose from duplication of the same gene, CYP734A1. Three out of four S-locus genes in Gelsemium elegans were hemizygous, as previously reported in Primula and Turnera. We provide genomic evidence on the genetic convergence of the supergene underlying distyly among distantly related angiosperm lineages and help to illuminate the genetic architecture involved in the evolution of heterostyly.