Skip to main content

Sex-differential reproduction success and selection on floral traits in gynodioecious Salvia pratensis



Gynodioecy, a sexual system with hermaphrodite and female individuals in a population, raises the question how the two sexual morphs are maintained. Salvia pratensis is a gynodioecious species featured by its modified stamens that act as a lever mechanism in pollination. Given sexual dimorphism in floral traits of the species, it is predictable that two sexual morphs differ in their interplay with pollinators and thus in their fitness. In this study, we investigated sex-specific reproduction success and floral adaptation in a population of S. pratensis.


We found that two sexual morphs in S. pratensis distinctly differed in their floral proportions. Female flowers fitted better to the pollinators than hermaphrodites in terms of touching the stigmas when being probed, and hence were more efficient in pollen deposition. Floral traits overall underwent stronger selection in the population, with stigma position and corolla length subject to disruptive selection mediated by different body-sized bumble bees; some selections on floral traits were significantly different in the strength, even opposite in the direction between two morphs. Flower production tended to be under correlational selection with floral structural traits, implying that a large plant with many flowers did not show an advantage in fitness unless its flower construction mechanically matched the pollinators well.


In conclusion, the pollinator-mediated selection likely played an important role in the evolution and maintenance of sexual dimorphism in the gynodioecious S. pratensis; and sex-divergent mechanical interaction with pollinators served as a critical mechanism by which female individuals were maintained in the population with a female advantage in pollen deposition efficiency (i.e. receiving pollen).


Gynodioecy is a dimorphic sexual system, in which hermaphroditic and female individuals coexist within a population [1]. It is rare, but widely distributed in angiosperm [2]. It occurs in at least 81 angiosperm families, but in far less than 1% of angiosperm species [2,3,4]. Hermaphrodites can transmit their genes via both seeds and pollen, whereas females do it only via seeds [1]. Because of this fitness disadvantage of females, it has been an intriguing topic in evolutionary biology to understand the evolution of gynodioecy, particularly the mechanism of maintaining females in gynodioecious population. Theoretically, it is believed that a female advantage, i.e. a higher female fitness in females over hermaphrodites, is necessary to compensate for the loss of male function in female individuals [5].

Female advantage can be, first, achieved by higher seed set due to more resources which are otherwise reallocated to pollen production [6]. Second, the female can gain the advantage of fitness by avoiding inbreeding depression [7, 8], which is sometimes considered as the main process responsible for gynodioecy [9]. Theoretically, female advantage can also occur by sex-differential interactions with biotic factors (reviewed by [5]), such as sex-biased seed predation [10, 11] and nectar robbing [12]. However, it has not been reported that the sex-specific interplay with pollinators can contribute to female advantage in fitness. On the contrary, the disadvantages in females have been documented much more frequently, which are mainly caused by low pollinator attraction [13,14,15] and/or lack of pollen as a reward [16,17,18,19].

Pollinator-mediated selection has been widely believed to play a key role in floral divergence not only among populations (or species) [20], but also within a population [21, 22]. Gender dimorphism, i.e. between-sex floral divergence beyond sexual organs, universally exists in gynodioecious species [18, 23]. Considering close interactions between flowers and pollinators in the process of pollination, it is predictable that both sexual morphs in gynodioecious population differ in their interplay with pollinators, and thereby subject to different selection pressure on floral traits. Hermaphrodite flowers are usually larger and more attractive than the females due to their greater allocation to attraction (e.g. petal size, nectar or pollen) [24, 25]. Thus, it is likely that the hermaphrodite evolves toward functional male provided it benefits more in male function than female function from increased pollinator visitation [26, 27]. Apart from sexual difference in pollination attraction, the sex-differential mechanical fitting to pollinators could be another source that generates difference in reproductive success between sexes. Therefore, it will be informative to investigate sex-specific adaption of floral traits by pollination ecology and phenotypic selection, and can provide insights into the evolutionary processes and maintenance of sexual dimorphism in gynodioecy.

Salvia pratensis is a perennial gynodioecious species mainly distributed in Europe. The species is self-compatible with a mixed mating system [28], and primarily pollinated by bumble bees and/or honey bees [29, 30], (Fig. 1). As with most species in the genus, two stamens of each hermaphrodite flower are modified to lever-like structures functioning as a lever mechanism in pollination, with the lower theca of each stamen reduced and the upper one fertile [31], (Fig. 1b). In female flowers, the thecae or even entire upper arms of the stamens are reduced, resulting in a nonfunctional lever mechanism. Regarding the lever mechanism, it has been argued that the interplay between flowers and pollinators is vulnerable against variation in each of the interactive parts, and that floral traits are under strong selection, as already documented in S. digitaloides [32], (Fig. 1c). This means, minute changes in the proportions of the pollinators or floral structures may have significant consequences for pollination success and hence the fitness [29].

Fig. 1

Floral traits and dorsal pollination in Salvia pratensis. Flower morphometrics (a), dorsal pollinations for hermaphrodites (b) and females (c). cl, corolla length; tl, corolla tube length; ste, style exsertion; sth, style height; st, stigma; ufa, upper fertile anther; ula, upper lever arm; rua, reduced upper anther

In the present study, we presume that the gynodioecious S. pratensis may serve as an ideal model for the study on the evolution of sexual dimorphism. Given that floral traits are different between sexes, we predict that two sexual morphs exhibit different interactions with pollinators, and thus different pollination efficiency and reproduction success. Consequently, the floral traits will be subject to sex-differential selection mediated by the pollinators. Specifically, we address the following three questions. 1) Do floral traits exhibit sexual dimorphism in gynodioecious S. pratensis? 2) How does selection act on floral traits and differ between sexual morphs? 3) What is the mechanism of maintaining females in the gynodioecious population?


Sexual dimorphism and pollinator assemblage in the gynodioecious S. pratensis

Floral traits significantly differed between two sexual morphs in the population of S. pratensis (Table 1). The female had on average a shorter corolla and tube, a smaller platform (i.e. flower mouth), and a stigma closer to the platform than hermaphrodite flowers. Style exsertion, flower production and stalk diameter did not significantly differ between two morphs.

Table 1 Morphometric data of floral traits in gynodioecious Salvia pratensis

Pollinator assemblage in the population consisted of both queens and workers of B. terrestris, workers of B. sylvarum and queens of B. lapidarius. Among them, B. terrestris was the dominant species, with a frequency of more than 75.9% (Fig. 2). Workers of B. terrestris had a mean body length of 16.07 ± 0.63 mm, thorax thickness of 5.38 ± 0.22 mm and tongue length of 4.14 ± 0.21 mm (n = 4). The queens of B. terrestris and B. lapidaries accounted for 17.4% of all bumblebees. They were distinctly bigger than the other pollinators, and the body length, thorax thickness and tongue length of B. terrestris queens were 27.2 mm, 7.82 mm and 5.15 mm, respectively.

Fig. 2

Pollinator assemblage in gynodioecious Salvia pratensis

Sex-differential pollination and reproduction success in S. pratensis

The queens of bumble bees were distinctly more successful in touching stigma when foraging than the workers given their relatively bigger body proportions. The successful visits, i.e. the probes touching the stigmas, accounted for 86.3% of total probes (n = 87). Two sexual morphs significantly differed in the percentage of successful visits by the dominant pollinators (i.e. B. terrestris) (Fig. 3). For the female flowers, the percentage of touching-stigma visits was over 65% in total visits, whereas the value was less than 30% for the hermaphrodite flowers.

Fig. 3

Difference in pollination efficiency (by Bombus terrestris) between two sexual morphs in gynodioecious Salvia pratensis

As showed in Fig. 4, the hermaphrodite individuals in the population had a fruit set of 39.5% ± 2.9% (mean ± se, n = 42), and seed set per flower of 0.76 ± 0.07 (mean ± se, n = 42). Female individuals had a higher fruit set of 45.3% ± 3.5% (mean ± se, n = 36), and higher seed set of 0.87 ± 0.08 (mean ± se, n = 36) than the hermaphrodites; but each of them did not significantly differ between two morphs.

Fig. 4

Fruit set (a) and seed (nutlet) set per flower (b) of two sexual morphs in gynodioecious Salvia pratensis

Sex-differential selection on floral traits in S. pratensis

The complete model for selection gradient analysis explained about 86.1% of total variation in female fitness among hermaphrodite group (F = 5.88, p = 0.0001), and 85.4% of total variation among female group. The ANCOVA model for selection comparison between two morphs was also highly significant, explaining about 87.8% of total variation (F = 2.88, p = 0.004).

Among five floral traits determined, only style exsertion was subject to significant directional selection among female plants in the population. Its gradient was 0.52 ± 0.20 (p = 0.034), but not significantly different with that among hermaphrodite plants. As for nonlinear selection, positive quadratic selection was consistently detected on corolla length among each sexual group. The selection gradient was 2.17 ± 0.96 (p = 0.05) among female individuals, significantly higher than 0.58 ± 0.23 (p = 0.022) among hermaphrodites (Table 2, Fig. 5a). Significant quadratic selection was also detected on style exsertion among hermaphrodite individuals (p = 0.006), but the selection gradient did not differ significantly with the female’s. Strong quadratic selection acted consistently on style height in both sexual groups. The selection gradient for the female was 4.43 ± 1.13 (p = 0.004), significantly higher than 0.92 ± 0.22 (p < 0.001) for the hermaphrodite (Table 2, Fig. 5b).

Table 2 Selection gradients for floral traits and comparison between sexual morphs in Salvia pratensis
Fig. 5

Standardized phenotypic selection gradients that were significantly different between two sexual morphs in gynodioecious Salvia pratensis (a-i). Selection gradients were illustrated with the added-variable plots, in which the residuals from a complete model of relative fitness on all traits except the focal trait are plotted against the residuals from a regression model of the focal trait on the other traits. The relationships between relative fitness and trait values were illustrated for hermaphrodite morph with open symbols and dashed line, and for female morph with filled symbols and solid line. The asterisks indicate significant difference between morphs: “*” for p < 0.05, “**” p < 0.01, and “***” p < 0.001

Positive correlational selection was detected on style height and platform, and negative correlational selection on style height and corolla length among either sex-morph individuals; the selection gradient for each combination significantly differed between two sexual morphs (Table 2, Fig. 5c, e). Style exsertion and corolla length was subject to positive correlational selection only among the hermaphrodites, with a selection gradient of 0.77 ± 0.21 (p = 0.001), which was not significantly different from that among the females. Negative correlational selection was detected on the combination of style exsertion and platform among hermaphrodite individuals. Marginally significant correlational selection was detected on the combination of style height and exsertion among female group. Negative correlational selection was consistently detected on flower number (i.e. flower production) and corolla length among either sex-morph group. The selection gradient was − 0.45 ± 0.22 (p = 0.0006) among the hermaphrodites, significantly higher than − 1.06 ± 0.44 (p = 0.042) among the females (Table 2, Fig. 5f). Positive correlational selection was consistently detected on flower production and style height among either sexual group; however, the added-variable plot reflected distinct selection only among female individuals, and there existed significant difference between two morphs (Table 2, Fig. 5i). The combination of flower production and platform was favored only among female group, subject to positive correlational selection. Flower production and style exsertion underwent opposite correlational selection among different sex-morph group (Fig. 5h). The selection gradient was − 0.44 ± 0.14 among the hermaphrodites, whereas it was 0.63 ± 0.23 among the females (Table 2).


Evolutionary divergence in floral traits of gynodioecious S. pratensis

Gynodioecious species usually exhibit gender divergence in their floral traits. Often, the size of flowers is larger in hermaphrodites than in females [18, 23]. In present study, it was also found that floral traits (e.g. corolla, tube length and stigma height) were significantly larger in hermaphrodites than in females of S. pratensis. There have been two nonexclusive explanations for the evolution of dimorphism in flower size of gynodioecious species. Firstly, the species with few ovules usually have relatively small female flowers [18], because they either rarely suffer from pollen limitation and thus need not to invest too much in flower size for pollinator attraction, or have less floral structures to protect [24, 33]. Secondly, small female flowers likely contribute to female advantage by saving resources for seed production (i.e. resource compensation) [34, 35], while the larger hermaphroditic flowers enhance pollinator attraction and thereby male function [36, 37].

Apart from floral size, flower production was also usually sexually dimorphic to different extent in dimorphic species [38, 39]. However, no evidence indicated sexual dimorphism of flower production in the S. pratensis population. Besides the influence of sampling on similarity of plant size in the study, highly intersexual genetic correlation could be one of the main reasons for the absence of sexual dimorphism in flower number [40], which constrained its divergence between two sexual morphs. Further, it was also reasonable that two sexual morphs had relatively consistent flower production in the gynodioecious population. Because flower number of an inflorescence usually affects floral display and thus pollinator attraction, the similar production of flowers in both morphs may assure pollinators foraging all individuals in a sex-unbiased way [19]. Otherwise, if either sexual morph has a much smaller floral display than the other, the resulting decreased visitation by pollinators will affect pollen export for the hermaphrodites, i.e. male fitness, and reduce the opportunity and efficiency of pollination for the females.

Sex-specific reproduction success in gynodioecious S. pratensis

Theoretically, female advantage is a prerequisite for the maintenance of females in a gynodioecious population to compensate for the lost male function [1, 2]. In present study, the female flowers had relatively smaller floral structures (particularly, a lower stigma position) such that they matched distinctly better with the pollinators than did the hermaphrodite flowers. In consequence, the female individuals had an advantage over the hermaphrodites in pollen deposition onto stigmas, i.e. receiving pollen for the realization of female fitness. However, both fruit and seed sets in the females were only slightly higher than those in the hermaphrodites. The reason was probably that the female was prone to pollen limitation due to the lack of reproductive assurance through self-pollination [19, 41]. In the population studied, each sexual morph of S. pratensis had a large floral display, and bumble bees tended to probe flowers in sequence within an inflorescence (also see [42]). Therefore, it was likely that the females underwent local pollen depletion (i.e. carryover) [43] when pollinators constantly probed flowers within or among female plants [44, 45]. As such, the female advantage in pollen deposition might be counteracted by pollen limitation, provided no adequate pollen import.

Although theoretically required for the gynodioecy, female advantage was not always found in natural populations of gynodioecious species [5], such as Beta vulgaris [46], Daphne laureola [47] and Raphanus sativus [48]. For the absence of female advantage, a convincing explanation is that these populations are being in a phase of evolutionary dynamics, which is featured by a high frequency of females. As discussed above, the females are prone to suffer from pollen limitation in such case [49] and hence have a declined and even undetectable female advantage due to overall drop in pollen availability [5]. Nevertheless, such a slight female advantage observed in the S. pratensis population does not still violate the least requirement for the maintenance of nuclear-cytoplasmic gynodioecy [50, 51], although no detailed information has been known on the sex determination of the species.

Selection on floral traits in gynodioecious S. pratensis

Sex-specific selection has been considered as one of the primary evolutionary processes responsible for the evolution and maintenance of sexual dimorphism, interplaying with the breakdown of intersexual genetic correlations that constrain between-sex divergence of traits [40, 52]. In present study, although the sample size was relatively small, significant selections were still detected on floral traits in the gynodioecious S. pratensis; and the strength and pattern of selections for some traits exhibited distinct differences between both sexual morphs.

Flower production determines the capacity of plants for mating opportunities through male and female functions, and thus, tends to be under positive directional selection [53, 54]. However, in present study, significantly directional selection was not detected on flower production among either sex-morph individuals. That is, female fitness did not positively respond to the increase in flower production. Partially, the sampling of plants with similar size in the study might have artificially reduced variation in flower number in both sexual groups, and therefore potentially hided evidence for selection on flower number. Besides, severe pollen limitation might be another cause, due to high frequency of females in the gynodioecious population. In such case, the positive correlation between female fitness and flower production would become weak and even disappeared (see [42]). This argument was further supported by the findings of correlational selections on flower number and other floral traits as followed.

In the population, flower production (i.e. flower number) tended to be under correlational selection with floral structural traits, like platform, and style height and exsertion determining stigma position. This implied that the contribution of flower production to female fitness was affected by a given construction of flowers in the plant. Among hermaphrodite individuals, negative correlational selection acted on flower number and style exsertion, indicating that a large plant with many flowers would has a high female fitness given with a relatively short style exsertion. As hermaphrodite flowers were significantly larger in floral structures than female flowers, such as longer corolla and larger platform (i.e. relatively larger flower mouth), thereby, bumble bees could land and enter the flowers in a horizontal manner when probing. In this case, a relatively shorter style exsertion could facilitate bumble bees probing the flowers, having the stigma positioned right above bumble bee’s back and not at flower entrance as a barrier. In contrast, female flowers were distinctly small in structures, particularly in platform; thereby, bumble bees landed the flowers in a vertical manner and hung on the low lip when probing. In this case, it was not surprising that selection favored positive correlations between flower number with platform size, style exsertion or stigma high among female plants, as the findings showed. That is, a large platform and/or relatively far stigma from platform were beneficial to the fitness realization of a large female plant. It was reasonable because such a flower construction could not only provide enough space for bumble bee landing by a horizontal way, but also make the stigma easier to touch. Taken together, given a large inflorescence, many more flowers produced could not be transformed into the advantage in female fitness unless the flower construction mechanically fitted pollinators well.

Finally, our findings indicated that stigma height tended to be under disruptive selection among either hermaphrodite or female individuals, and the selection strength was much stronger in the female group. It was implied that the plants with either high- or low-positioned stigma were equally favored in the population. Meanwhile, flower size (i.e. corolla length) was also subject to disruptive selection among either sex-morph individuals. Considering two types of bumble bees with different body size in the population, it could be reasonably interpreted that the two types of flowers stood equal chance of pollination success due to their divergent adaption to both sizes of bumble bees. Besides, significantly correlational selection for platform and style height also reflected divergent adaptation of the trait pair to different pollinators. The flowers with a large platform and a high stigma could fit well to big bumble bees, whereas the flowers with the two smaller traits fit better to small bumble bees. As such, the disruptive selection mediated by pollinators likely played an important role in the evolution and maintenance of sexual dimorphism in the gynodioecious S. pratensis.


Gynodioecy has long drawn biologist’s attention to the question how females are maintained in the population. Theoretically, a female advantage is necessary for the maintenance of gynodioecy. In the study, we found that floral traits significantly differed between two sexual morphs in the gynodioecious S. pratensis. Sexual divergence in flower size conferred female individuals an advantage in pollen deposition (i.e. receiving pollen) over the hermaphrodites, because relatively smaller flowers of the female’s fit better to the pollinators. However, the females just gained a slightly higher fitness than the hermaphrodites due to their intrinsic disadvantage in pollen availability during pollination; that is, pollen limitation could be one of the main reasons for the weak difference in female fitness between two sexual morphs. Therefore, it was predictable that the female advantage in fitness varied with population dynamics (e.g. sex ratio) and pollinator’s activity. Floral traits overall underwent strong selection in the gynodioecious population, with flower size and stigma position subject to disruptive selection. Flower production tended to be under correlational selection with floral structural traits, implying that many more flowers in a large plant could not be transformed into the advantage in fitness unless the flower construction mechanically matched pollinators well (i.e. efficient in pollination). In conclusion, the pollinator-mediated selection likely played an important role in the evolution and maintenance of sexual dimorphism in the gynodioecious S. pratensis, and the sex-divergent mechanical interaction with pollinator served as a mechanism by which female individuals, with an advantage in pollen deposition efficiency, were maintained in the gynodioecious population.

Materials and methods

Study species and site

The study was conducted during the main flowering season of Salvia pratensis in May 2016 at Mainz, Germany. The studied population was located close to a small shrubbery in the midst of farmland at Mainz-Hechtsheim (49°56′23″N, 8°15′04″E; 197 m a.s.l.), with an area of 45 × 3.5 m2. The whole area was densely covered by individuals of S. pratensis (over 500 individuals). The female and hermaphrodite individuals were almost evenly distributed in the population, and their relative proportion (i.e. sex ratio) was close to 50%.

Experimental approach

In the mid of the flowering season, we randomly labeled 80 individuals (44 hermaphrodites and 36 females) in the gynodioecious population. These labeled plants were about the same size and flowered synchronously. For morphometric measurements, three or four completely opened flowers were randomly selected in each labeled inflorescence as the representatives of the plant. Five floral structural traits were measured in each flower (see Fig. 1): corolla length (cl), corolla tube length (tl), style height (sth) defined as the distance between the lower lip and the stigma, style exsertion (se), and platform size determined by subtracting “cl” from “tl”. The diameter of the main axis under the first branch was measured as a proxy of plant size.

Each tagged plant was collected about 20 days after flower withering when nutlets were matured but still enclosed in the calyx. Flower number was determined by counting the number of pedicels which remained on the inflorescence after flower withering. Seed (nutlet) number per inflorescence was counted as an estimate of female fitness component for each plant. All traits were measured by using digital calipers with ±0.01 mm of error.

During the experimental period, we chose three sunny days (each time from 10:00 to 14:00) to identify the range and relative frequency of different pollinator species, observing their behavior of pollination. The frequency of different pollinator species was recorded by observing bees presented in a defined area (2 × 2 m2) or along a transverse sector (50 m2) in the population. For the mechanical match between pollinators and flowers, we stochastically targeted a plant (or a flower) being probed (totally for 400 flowers), followed by recording its sex morph and observing whether the forager could touch its stigma. Three to five individuals of each pollinator species were collected for species identification and for morphological measurements including body length, thorax thickness (i.e. body thickness in the part of thorax), thorax width and tongue length (naturally extended length).

Data analysis

We employed linear model in R (3.5.0) to determine difference in floral traits between sexual morphs in the gynodioecious S. pratensis. We calculated fruit set as fruit number per flower and seed set as seed (nutlet) number per flower per inflorescence. Generalized linear regression (GLM) was used to test the differences in fruit set and seed set (quasibinominal family, logit link function) between both morphs. We measured pollination efficiency of pollinators by calculating relative frequency of touching-stigma visits in total visits, and determined its difference between morphs with Pearson’s chi-squared test.

Multivariate regression analysis was used to estimate the strength and pattern of selection on the floral traits [55,56,57]. We employed the most complete regression model for estimates of different selection gradients (see [58]). The complete model includes linear and quadratic terms of each trait, and the products of pairs of all traits included. The partial regression coefficient of each trait serves as the main indicator of directional selection (β), and the double value of coefficient for a trait squared as a measure of nonlinear selection on the trait (γ, i.e. nonlinear selection gradient); and a significant coefficient of the product of a pair of traits as a measure of selection on combination of the two traits (i.e. correlational selection gradient) [58, 59]. We performed separate analyses for each sexual morph. Before the analyses, the seed number of each individual was standardized to the relative seed number (i.e. individual seed number / morph mean) as estimate of the female fitness component. Each trait value of an individual plant was the mean of 3 or 4 flowers. The trait values for each morph were standardized (i.e. mean = 0 and standard deviation = 1) by subtracting morph mean from the individual trait value and then divided by the standard deviation. Because of significant correlations of corolla length with tube length and corolla platform (each correlation coefficient > 0.86), we used relative platform (r-platform) rather than platform by dividing platform by corolla length in selection analysis, as a measurement of flower mouth opening (i.e. 1- tube/corolla). That is, we introduced flower number, style height, style exsertion, corolla and r-platform without corolla tube in the model for selection analysis. In addition, considering the effect of resource availability on plant’s size and thus on both the size of floral traits and fecundity, we introduced stalk diameter as a variable of plant size in the model to eliminate maternal effect on selection for traits [53, 60].

After establishing the full model, we further obtained the least adequate model by step regression. The significant terms in the full model were almost consistent in the number and sign with those in the least adequate model except that the significance was enhanced. Therefore, to maintain completeness and consistence of the models for two morphs, we presented the results from the full model in the study. Finally, we performed comparison of selection on floral traits between two sexual morphs, by establishing ANCOVA models with sexual morph as a factor. For all the models, we checked the normality of error distribution by Shapiro Test (W = 0.96, p = 0.19 for female group; W = 0.98, p = 0.55 for hermaphrodite group), the variance constancy by Non-constant Variance Score Test (Chisquare = 0.01, p = 0.92 for female group; Chisquare = 0.13, p = 0.71 for hermaphrodite group). We also assessed the assumptions of linear model using the Global Test with the Package gvlma in R, and each of assumptions was acceptable (p > 0.05). Kappa values for correlation matrix of the traits included in the models were 8.50 for female group, and 11.3 for hermaphrodite group, so there was no evidence for a problem of multi-collinearity. Finally, we illustrated selection gradients that were significantly different between two sexual morphs with the Added-Variable Plots [61]. Except for the standardization of variables, no transformations were done. For all data analysis, R version 3.5.0 [62] was used.

Availability of data and materials

Data generated or analyzed during this study are either included in this published article or are available from the corresponding author on reasonable request.



Analysis of covariance


Analysis of variance


corolla length


Generalized linear regression


platform size


reduced upper anther


style exsertion


style height, i.e. the distance between the lower lip and the stigma


corolla tube length


upper fertile anther


upper lever arm


  1. 1.

    Darwin C. The different forms of flowers on plants of the same species. London: John Murray; 1877.

    Google Scholar 

  2. 2.

    Rivkin LR, Case AL, Caruso CM. Why is gynodioecy a rare but widely distributed sexual system? Lessons from the Lamiaceae. New Phytol. 2016;211(2):688–96.

    PubMed  Article  Google Scholar 

  3. 3.

    Dufay M, Champelovier P, Kafer J, Henry J-P, Mousset S, Marais GAB. An angiosperm-wide analysis of the gynodioecy-dioecy pathway. Ann Bot. 2014;114(3):539–48.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  4. 4.

    Caruso CM, Eisen K, Case AL, Mazer SJ. An angiosperm-wide analysis of the correlates of gynodioecy. Int J Plant Sci. 2016;177(2):115–21.

    Article  Google Scholar 

  5. 5.

    Dufay M, Billard E. How much better are females? The occurrence of female advantage, its proximal causes and its variation within and among gynodioecious species. Ann Bot. 2012;109(3):505–19.

    PubMed  Article  Google Scholar 

  6. 6.

    Ashman TL, Williams CF. Constraints on the evolution of males and sexual dimorphism: field estimates of genetic architecture of reproductive traits in three populations of gynodioecious Fragaria virginiana. Evolution. 2003;57(9):2012–25.

    PubMed  Article  Google Scholar 

  7. 7.

    Mutikainen P, Delph LF. Inbreeding depression in gynodioecious Lobelia siphilitica: among-family differences override between-morph differences. Evolution. 1998;52(6):1572–82.

    PubMed  Google Scholar 

  8. 8.

    Weller SG, Sakai AK. Selfing and resource allocation in Schiedea salicaria (Caryophyllaceae), a gynodioecious species. J Evol Biol. 2005;18(2):301–8.

    CAS  PubMed  Article  Google Scholar 

  9. 9.

    Kubota S, Ohara M. Discovery of male sterile plants and their contrasting occurrence between self-compatible and self-incompatible populations of the hermaphroditic perennial Trillium camschatcense. Plant Species Biol. 2009;24(3):169–78.

    Article  Google Scholar 

  10. 10.

    Collin C, Pennings P, Rueffler C, Widmer A, Shykoff J. Natural enemies and sex: how seed predators and pathogens contribute to sex-differential reproductive success in a gynodioecious plant. Oecologia. 2002;131(1):94–102.

    CAS  PubMed  Article  Google Scholar 

  11. 11.

    Marshall M, Ganders FR. Sex-biased seed predation and the maintenance of females in a gynodioecious plant. Am J Bot. 2001;88(8):1437–43.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    Zhang YW, Yang CF, Zhao JM, Guo YH. Selective nectar robbing in a gynodioecious plant (Glechoma longituba) enhances female advantage. J Evol Biol. 2009;22(3):527–35.

    PubMed  Article  Google Scholar 

  13. 13.

    Delph LF, Lively CM. Pollinator visitation, floral display, and nectar production of the sexual morphs of a gynodioecious shrub. Oikos. 1992;63(2):161–70.

    Article  Google Scholar 

  14. 14.

    Williams CF, Kuchenreuther MA, Drew A. Floral dimorphism, pollination, and self-fertilization in gynodioecious Geranium richardsonii (Geraniaceae). Am J Bot. 2000;87(5):661–9.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Ashman TL, Stanton M. Seasonal variation in pollination dynamics of sexually dimorphic Sidalcea oregana ssp. spicata (Malvaceae). Ecology. 1991;72(3):993–1003.

    Article  Google Scholar 

  16. 16.

    Stevens DP. On the gynodioecious polymorphism in Saxifraga granulata L.(Saxifragaceae). Biol J Linn Soc. 1988;35(1):15–28.

    Article  Google Scholar 

  17. 17.

    Jordano P. Pollination biology of Prunus mahaleb L.: deferred consequences of gender variation for fecundity and seed size. Biol J Linn Soc. 1993;50(1):65–84.

    Article  Google Scholar 

  18. 18.

    Shykoff JA, Kolokotronis SO, Collin CL, Lopez-Villavicencio M. Effects of male sterility on reproductive traits in gynodioecious plants: a meta-analysis. Oecologia. 2003;135(1):1–9.

    PubMed  Article  Google Scholar 

  19. 19.

    Arnan X, Escola A, Rodrigo A, Bosch J. Female reproductive success in gynodioecious Thymus vulgaris: pollen versus nutrient limitation and pollinator foraging behaviour. Bot J Linn Soc. 2014;175(3):395–408.

    Article  Google Scholar 

  20. 20.

    Herrera CM, Castellanos MC, Medrano M. Geographical context of floral evolution: towards an improved research programme in floral diversification. In: Harder LD, Barrett SCH, editors. Ecology and evolution of flowers. New York: Oxford University Press; 2006. p. 278–94.

    Google Scholar 

  21. 21.

    Sletvold N, Agren J. There is more to pollinator-mediated selection than pollen limitation. Evolution. 2014;68(7):1907–18.

    PubMed  Article  Google Scholar 

  22. 22.

    Ferdy JB, Gouyon PH, Moret J, Godelle B. Pollinator behavior and deceptive pollination: learning process and floral evolution. Am Nat. 1998;152(5):696–705.

    CAS  PubMed  Article  Google Scholar 

  23. 23.

    Eckhart VM. Sexual dimorphism in flowers and inflorescences. In: Geber MA, Dawson TE, Delph LF, editors. Gender and sexual dimorphism in flowering plants. Berlin Heidelberg: Springer; 1999. p. 123–48.

    Google Scholar 

  24. 24.

    Delph LF, Galloway LF, Stanton ML. Sexual dimorphism in flower size. Am Nat. 1996;148(2):299–320.

    Article  Google Scholar 

  25. 25.

    Asikainen E, Mutikainen P. Preferences of pollinators and herbivores in Gynodioecious Geranium sylvaticum. Ann Bot. 2005;95(5):879–86.

    PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    Ashman TL. Pollinator selectivity and its implications for the evolution of dioecy and sexual dimorphism. Ecology. 2000;81(9):2577–91.

    Article  Google Scholar 

  27. 27.

    Morgan MT. Attractive structures and the stability of hermaphroditic sex expression in flowering plants. Evolution. 1992;46(4):1199–213.

    PubMed  Article  Google Scholar 

  28. 28.

    Van Treuren R, Bijlsma R, Ouborg NJ, Van Delden W. The effects of population size and plant density on outcrossing rates in locally endangered Salvia pratensis. Evolution. 1993;47(4):1094–104.

    PubMed  Article  Google Scholar 

  29. 29.

    Claßen-Bockhoff R, Speck T, Tweraser E, Wester P, Thimm S, Reith M. The staminal lever mechanism in Salvia L. (Lamiaceae): a key innovation for adaptive radiation? Organ Divers Evol. 2004;4(3):189–205.

    Article  Google Scholar 

  30. 30.

    Reith M, Baumann G, Claßen-Bockhoff R, Speck T. New insights into the functional morphology of the lever mechanism of Salvia pratensis (Lamiaceae). Ann Bot. 2007;100(2):393–400.

    PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Claßen-Bockhoff R, Wester P, Tweraser E. The staminal lever mechanism in Salvia L. (Lamiaceae): a review. Plant Biology (Stuttgart). 2003;5(1):33–41.

    Article  Google Scholar 

  32. 32.

    Zhang B, Li QJ. Phenotypic selection on the staminal lever mechanism in Salvia digitaloides (Labiaceae). Evol Ecol. 2014;28(2):373–86.

    Article  Google Scholar 

  33. 33.

    Thompson JD, Rolland AG, Prugnolle F. Genetic variation for sexual dimorphism in flower size within and between populations of gynodioecious Thymus vulgaris. J Evol Biol. 2002;15(3):362–72.

    Article  Google Scholar 

  34. 34.

    Eckhart VM. Resource compensation and the evolution of gynodioecy in Phacelia linearis (Hydrophyllaceae). Evolution. 1992;46(5):1313–28.

    PubMed  Article  Google Scholar 

  35. 35.

    Eckhart VM. Do hermaphrodites of gynodioecious Phacelia linearis (Hydrophyllaceae) trade off seed production to attract pollinators? Biol J Linn Soc. 1993;50(1):47–63.

    Article  Google Scholar 

  36. 36.

    Conner JK, Rush S. Effects of flower size and number on pollinator visitation to wild radish, Raphanus raphanistrum. Oecologia. 1996;105(4):509–16.

    PubMed  Article  Google Scholar 

  37. 37.

    Campbell DR. Comparing pollen dispersal and gene flow in a natural population. Evolution. 1991;45(8):1965–8.

    PubMed  Article  Google Scholar 

  38. 38.

    Delph LF, Knapczyk FN, Taylor DR. Among-population variation and correlations in sexually dimorphic traits of Silene latifolia. J Evol Biol. 2002;15(6):1011–20.

    Article  Google Scholar 

  39. 39.

    Golonka AM, Sakai AK, Weller SG. Wind pollination, sexual dimorphism, and changes in floral traits of Schiedea (Caryophyllaceae). Am J Bot. 2005;92(9):1492–502.

    PubMed  Article  Google Scholar 

  40. 40.

    Barrett SCH, Hough J. Sexual dimorphism in flowering plants. J Exp Bot. 2013;64(1):67–82.

    CAS  PubMed  Article  Google Scholar 

  41. 41.

    Asikainen E, Mutikainen P. Pollen and resource limitation in a gynodioecious species. Am J Bot. 2005;92(3):487–94.

    PubMed  Article  Google Scholar 

  42. 42.

    Zhang B, Claßen-Bockhoff R, Zhang ZQ, Sun S, Luo YJ, Li QJ. Functional implications of the staminal lever mechanism in Salvia cyclostegia (Lamiaceae). Ann Bot. 2011;107(4):621–8.

    PubMed  PubMed Central  Article  Google Scholar 

  43. 43.

    Robertson AW. The relationship between floral display size, pollen carryover and geitonogamy in Myosotis colensoi (Kirk) Macbride (Boraginaceae). Biol J Linn Soc. 1992;46(4):333–49.

    Article  Google Scholar 

  44. 44.

    McCauley DE. The genetic structure of a gynodioecious plant: nuclear and cytoplasmic genes. Evolution. 1998;52(1):255–60.

    PubMed  Article  Google Scholar 

  45. 45.

    Klaas AL, Olson MS. Spatial distributions of cytoplasmic types and sex expression in Alaskan populations of Silene acaulis. Int J Plant Sci. 2006;167(2):179–89.

    CAS  Article  Google Scholar 

  46. 46.

    De Cauwer I, Arnaud JF, Courseaux A, Dufay M. Sex-specific fitness variation in gynodioecious Beta vulgaris ssp. maritima: do empirical observations fit theoretical predictions? J Evol Biol. 2011;24(11):2456–72.

    PubMed  Article  Google Scholar 

  47. 47.

    Alonso C, Herrera CM. Neither vegetative nor reproductive advantages account for high frequency of male-steriles in southern Spanish gynodioecious Daphne laureola (Thymelaeaceae). Am J Bot. 2001;88(6):1016–24.

    CAS  PubMed  Article  Google Scholar 

  48. 48.

    Murayama K, Yahara T, Terachi T. Variation of female frequency and cytoplasmic male-sterility gene frequency among natural gynodioecious populations of wild radish (Raphanus sativus L.). Mol Ecol. 2004;13(8):2459–64.

    CAS  PubMed  Article  Google Scholar 

  49. 49.

    Delph LF, Touzet P, Bailey MF. Merging theory and mechanism in studies of gynodioecy. Trends Ecol Evol. 2007;22(1):17–24.

    PubMed  Article  Google Scholar 

  50. 50.

    Bailey MF, Delph LF, Lively CM. Modeling gynodioecy: novel scenarios for maintaining polymorphism. Am Nat. 2003;161(5):762–76.

    PubMed  Article  Google Scholar 

  51. 51.

    Gouyon PH, Vichot F, Van Damme JMM. Nuclear-cytoplasmic male sterility: single-point equilibria versus limit cycles. Am Nat. 1991;137(4):498–514.

    Article  Google Scholar 

  52. 52.

    Lande R. Sexual dimorphism, sexual selection, and adaptation in polygenic characters. Evolution. 1980;34(2):292–305.

    PubMed  Article  Google Scholar 

  53. 53.

    Toland Ø. Environment-dependent pollen limitation and selection on floral traits in an alpine species. Ecology. 2001;82(8):2233–44.

    Article  Google Scholar 

  54. 54.

    Benitez-Vieyra S, Glinos E, Medina AM, Cocucci AA. Temporal variation in the selection on floral traits in Cyclopogon elatus (Orchidaceae). Evol Ecol. 2012;26(6):1451–68.

    Article  Google Scholar 

  55. 55.

    Arnold SJ, Wade MJ. On the measurement of natural and sexual selection: theory. Evolution. 1984;38(4):709–19.

    PubMed  Article  Google Scholar 

  56. 56.

    Arnold SJ, Wade MJ. On the measurement of natural and sexual selection: applications. Evolution. 1984;38(4):720–34.

    PubMed  Article  Google Scholar 

  57. 57.

    Lande R, Arnold SJ. The measurement of selection on correlated characters. Evolution. 1983;37(6):1210–26.

    PubMed  Article  Google Scholar 

  58. 58.

    Harder LD, Johnson SD. Darwin's beautiful contrivances: evolutionary and functional evidence for floral adaptation. New Phytol. 2009;183(3):530–45.

    PubMed  Article  Google Scholar 

  59. 59.

    Conner JK, Hartl DL. A primer of ecological genetics. Massachusetts, USA: Sinauer Associates; 2004.

    Google Scholar 

  60. 60.

    Rausher MD. The measurement of selection on quantitative traits: biases due to environmental covariances between traits and fitness. Evolution. 1992;46(3):616–26.

    PubMed  Article  Google Scholar 

  61. 61.

    Cook RD, Weisberg S. Applied regression including computing and graphics. New York: Wiley; 1999.

    Google Scholar 

  62. 62.

    R Development Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018.

    Google Scholar 

Download references


We thank Laura Leineweber and Joshua Wachlin for their assistance in the experiment, and also thank anonymous Reviewers and the Editors for their valuable comments and suggestions on the manuscript.


This research was funded by the National Natural Science Foundation of China (41461014; 31760703) and the China Scholarship Council, and supported by the Science and Technology Innovation Funds of Gansu Agricultural University-Special funds for discipline construction (GAU-XKJS-2018-001). Funding bodies were not involved in the design of the study; collection, analysis, interpretation of data and in writing the manuscript.

Author information




BZ developed the concept and performed the field work and data analysis. BZ and RCB interpreted the results and drafted the manuscript. Both authors reviewed and approved the final version of this manuscript.

Corresponding authors

Correspondence to Bo Zhang or Regine Claßen-Bockhoff.

Ethics declarations

Ethics approval and consent to participate

Salvia pratensis was not under protection in Germany, and to collect data of the species complied with the national guidelines. As to bees, it was permitted to collect specimens for research reasons in and around Mainz by the Struktur - und Genehmigungsdirektion Süd, Oppenheim.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhang, B., Claßen-Bockhoff, R. Sex-differential reproduction success and selection on floral traits in gynodioecious Salvia pratensis. BMC Plant Biol 19, 375 (2019).

Download citation


  • Evolutionary divergence
  • Female advantage
  • Gynodioecy
  • Pollination
  • Phenotypic selection
  • Reproductive performance
  • Salvia
  • Sexual dimorphism