- Research article
- Open Access
Negative frequency-dependent selection or alternative reproductive tactics: maintenance of female polymorphism in natural populations
BMC Evolutionary Biology volume 13, Article number: 139 (2013)
Sex-limited polymorphisms have long intrigued evolutionary biologists and have been the subject of long-standing debates. The coexistence of multiple male and/or female morphs is widely believed to be maintained through negative frequency-dependent selection imposed by social interactions. However, remarkably few empirical studies have evaluated how social interactions, morph frequencies and fitness parameters relate to one another under natural conditions. Here, we test two hypotheses proposed to explain the maintenance of a female polymorphism in a species with extreme geographical variation in morph frequencies. We first elucidate how fecundity traits of the morphs vary in relation to the frequencies and densities of males and female morphs in multiple sites over multiple years. Second, we evaluate whether the two female morphs differ in resource allocation among fecundity traits, indicating alternative tactics to maximize reproductive output.
We present some of the first empirical evidence collected under natural conditions that egg number and clutch mass was higher in the rarer female morph. This morph-specific fecundity advantage gradually switched with the population morph frequency. Our results further indicate that all investigated fecundity traits are negatively affected by relative male density (i.e. operational sex ratio), which confirms male harassment as selective agent. Finally, we show a clear trade-off between qualitative (egg mass) and quantitative (egg number) fecundity traits. This trade-off, however, is not morph-specific.
Our reported frequency- and density-dependent fecundity patterns are consistent with the hypothesis that the polymorphism is driven by a conflict between sexes over optimal mating rate, with costly male sexual harassment driving negative frequency-dependent selection on morph fecundity.
Evolutionary biologists have long studied visible polymorphisms as they are excellent model systems to examine micro evolutionary processes [1–3]. Polymorphisms with morphs co-existing at relatively stable frequencies appear to be common, but this phenomenon can only persist under a limited range of conditions, one of these being negative frequency-dependent selection (NFDS). NFDS arises when individuals of a rare morph experiences a higher fitness than those of a more common type . Over generations, and in absence of other mechanisms, NFDS should lead to a balanced polymorphism, typically with limited fluctuations along an equilibrium frequency of the involved morphs [5–7]. Classic examples of NFDS include coexistence of different colour morphs, to gain access to mates , to challenge predators  and to lower sexual conflict intensity . Although the idea of NFDS has been appreciated for decades [4, 11], relatively few empirical studies have tested the validity of this concept under natural, unmanipulated field conditions. Especially rare are studies which relate natural geographical variation in morph frequencies, a putative selective agent and fitness parameters of the involved morphs with one another [7, 10, 12–14].
Sex-limited polymorphisms represent excellent model systems to study the nature of diversifying selection and consequently they have been subject to a variety of experimental [15–17] and theoretical [18, 19] studies. Particularly popular are studies on male polymorphisms, whose maintenance tend to be explained by a fitness advantage to the rare morph relative to the common phenotypes in the competition over mates (e.g. sneakers do better when territorials predominate; reviewed by Oliveira et al. ). Over the last few decades, however, it has become clear that polymorphisms restricted to the female sex are more common in nature than previously thought [20, 21]. Yet the underlying mechanisms that maintain phenotypic and genetic variation within females remain unresolved in many cases. Female polymorphisms are often considered to have evolved as a counter adaptation to reduce costs of harassment imposed by mate-searching males; e.g. butterflies , diving beetles , African bat bugs , damselflies . The wider context of this proposed mechanism is sexual conflict over optimal mating rate. Evidently, females need males to fulfil their reproductive needs. However, obtrusive males may reduce female fitness by exceeding the females’ optimal number of matings [25–27]. Mate searching males are considered to face fewer cognitive challenges when confronted with only one, rather than multiple female phenotypes coexisting within populations , some of which may appear like males (i.e. andromorphs [29, 30]). Therefore on a phenotypic level, females may experience diversifying sexual selection to avoid sexual harassment .
If the above arguments hold, then it is highly likely that female polymorphism is maintained by NFDS driven by social interactions between sexes. In much the same way that predators form a search image for the most common cryptic prey type , increased male sexual interest has been observed towards the most common female morph in a given population [19, 32–34]. Recent studies showed an inverse relationship between morph-specific fecundity and morph-specific frequency in the population and suggested male harassment as most likely selective agent [7, 10]. Although quantification of this selective agent has been understudied in past studies, support may come from density-dependent effects on fitness. Indeed, the overall intensity of male harassment may rise with male density, either absolute or relative to female density (i.e. operational sex ratio), because male–female interactions occur more frequently under these conditions [35–37]. Thus a thorough investigation of NFDS in female polymorphic systems entails evaluating the role of female morph frequency together with male densities on fitness related parameters , which forms the first aim of the current study.
In addition, sex-limited polymorphisms are frequently considered as alternative reproductive tactics (ARTs), i.e. a discontinuous set of selected traits to maximize reproductive output in two or more alternative ways . Although repeatedly studied in males, recent observations in several species of owl , lizard [15, 38] and insect [21, 39, 40] collectively indicate that female morphs may also represent ARTs. For example in the polymorphic lizard Uta stansburiana, combined density-dependent and negative frequency-dependent interactions among conspecifics determine the relative success of orange (produce many small eggs, r-strategy) and yellow (fewer but larger eggs, K-strategy) throated females. Female morphs may therefore differently allocate resources towards fitness related traits, potentially resulting in trade-offs among life-history and/or physiological traits. Although the majority of the studies cited above provide a new promising research avenue, they should be treated with caution, since many of them are performed with limited spatial replicates, without temporal replicates and/or with small sample sizes within populations. This makes it difficult to reach firm conclusions as to whether female morphs represent ARTs, especially in spatial and temporal heterogeneous environments.
In this study, we examine morph-specific variation in fecundity (i.e. egg number, egg mass, clutch mass and relative body mass) under natural conditions in multiple years and across six populations, which show extreme variation in morph frequencies. The aim of this study is to evaluate two hypotheses which try to explain maintenance of female polymorphism. Based on the first hypothesis we expect that relative male density, as a proxy for intensity of male harassment [35–37] and meanwhile the suggested selective force, negatively affects overall female fecundity. In doing so, we simultaneously test whether a frequency-dependent fecundity advantage exists for the rare female morph due to a lower positive frequency-dependent male detection rate [7, 10, 34]. Testing the second hypothesis, we explore whether female morphs exhibit ARTs in which resources are allocated differently into qualitative (egg mass) or quantitative (egg number) fecundity traits [15, 38, 39] and meanwhile account for potential spatial and temporal variation.
Female polymorphism has been observed in more than a hundred damselfly and dragonfly species . Phenotypic ratios in laboratory cross experiments are consistent with the hypothesis that this polymorphism is genetically controlled by a single autosomal locus, with a number of alleles equal to the number of female morphs; reviewed by . Large geographic variation in frequencies and densities of males and female morphs has been described in several species [43–45], which allows us to investigate the role of social interactions in maintaining sex-limited polymorphisms. We studied the common North American damselfly, Nehalennia irene, for which male harassment estimates towards female morphs has been shown to vary in a positive frequency-dependent manner [34, Bots J, Iserbyt A, Hammers M, Van Gossum H, Sherratt TN: Frequency-dependent sexual selection in two intra-specific mimicry systems, in preparation]. Nehalennia irene is not an endangered nor a protected species (see COSEWIC, federal government Canada) and therefore our research complies with the Convention on Biological Diversity and the Convention on the Trade in Endangered Species of Wild Fauna and Flora. It is a small non-territorial species, which inhabits marshy or boggy waters  and exhibits a discrete polymorphism restricted to the female sex. Female morphs are easily classified into andromorphs or gynomorphs based on their body coloration. Thus, while andromorph females resemble the conspecific male’s body blue coloration and melanin pattern [47, 48], gynomorph females have distinctive yellowish lateral thorax sides and a less conspicuous abdominal melanin pattern; for colour figures see , for pictures see . The species has one generation per year, with reproduction occurring between early June and mid-August. After locating a potential mate, a male will attempt grasping the individual in the so-called tandem formation, i.e. when the male succeeds in attaching his anal appendages to the individual’s prothorax . This tandem formation can last several hours (AI, personal observation). If receptive, the female cooperates by bending her abdomen towards the male’s secondary genitalia (2nd abdominal segment) to form a ‘copulation wheel’ . Very little additional information on reproductive biology is known for this species, but we expect that similar to related damselfly species, ovarial follicles of adult females can develop within one or a few days into roughly two hundred mature eggs . In N. irene, females lay eggs in floating pieces of dead plant material while the female is tandem-guarded by her last successful male . Several clutches of eggs are laid throughout a female’s lifetime .
Study sites and sample procedures
Our previous work with N. irene indicated large spatial, yet limited temporal variation in female morph frequencies among populations. Specifically, andromorph frequencies range from 0 to >90% throughout the species’ distribution range over Canada [42, 54]. For our current aims we selected six study populations that differed significantly in social conditions (see Table 1). Frequency and density estimates were obtained in a manner similar to that described in Van Gossum et al.. In short, individuals were randomly captured with an insect net while walking slowly through the reproductive area, sweeping eight-shaped figures and recording the time elapsed. All caught males, andromorphs and gynomorphs were counted and marked with a permanent marker prior to release to avoid multiple counts. Andromorph frequency (proportion andromorph females), operational sex ratio (OSR, proportion males relative to females in the reproductive zone) and male density (number of males caught per minute) were calculated and collectively quantify the social environment. Calculating these parameters including or excluding immature individuals gave very similar results, see also . Hence, we here use data based on mature, thus reproductively active, individuals given that the aim of the current study deals with sexually active individuals. OSR and male density can be used as a proxy for overall male harassment in a given population [34, 35, 55, 56]. Each of the six populations were monitored during the reproductive season over three consecutive years (2007–2009).
Sample collection was carried out on several days throughout the reproductive season (mean ± SE: 4 ± 0.5 sample days; see Additional file 1) always between 9 am and 3 pm. This is the period before the majority of females start to oviposit (AI, personal observations). Andromorphs and gynomorphs were collected in an alternating manner to maintain a balanced design and to control for potential diurnal variation in egg number. We aimed to collect 25 adult andromorphs and 25 gynomorphs in each population for each investigated year. Measurements for relative body mass (see below) were performed in all three years (N = 43 ± 3 females per population and year; total N = 772) and the three other fecundity estimates were investigated in two successive years (N = 44 ± 3 females per population and year; total N = 547). For more detailed sample sizes per population and per year, (see Additional file 1). Mating status at the moment of capture was noted, i.e. being single or mating (i.e. involved in tandem or copulation). All individuals were stored for further measurements in 95% ethanol immediately after capture.
First, an individual was placed on a sheet of absorbent paper for exactly two minutes to allow standardised evaporation and absorption of most of the ethanol. Then it was weighted on a Kern & Sohn GmbH 870 balance (accuracy 0.1 mg). Next, a digital picture was taken (Nikon D70/Tamron macro lens 90 mm 1:2.8) of the right hind wing. Using ImageJ 1.38x , wing length was measured from the second antenodal cross vein to the stigma; see  for more details. Residuals of body mass were calculated by regressing body mass against wing length and were used as a measure for relative body mass (RBM). Positive values indicate relative heavy individuals for a given wing size, while negative values indicate relatively light individuals. RBM is not only considered an estimate for body condition, it is also suggested to increase with female fecundity in various insect taxa . Hence, RBM is here treated as a coarse measure of overall fecundity.
Dissection of specimens was performed under a Leica MZ 12.5 stereomicroscope. Abdominal sternites were removed and fifty developed eggs were isolated on a pre-weighted aluminum foil. This high number of eggs was chosen to account for potential variation in weight among eggs of the same clutch. Eggs were then dried in an incubator (Binder – APT.line™) at 60°C for 12 h and weighted on a Sartorius SE2F balance (accuracy 1 μg). Average dry weight of a single egg was calculated and further considered as a measure of egg quality since more nutrients are expected in heavier eggs [59, 60]. The total number of developed eggs was also counted for each specimen, and considered a quantitative measure of fecundity. Multiplying egg number with egg mass gave clutch mass as a comprehensive measure of fecundity (i.e. quality*quantity).
By measuring twice a set of randomly chosen individuals among populations and years, the repeatability of our measures could be evaluated. Repeatability was calculated as the proportion of the variation between individuals to the total variation, i.e. between and within individuals . A limited measurement error was observed: body mass (R = 0.90, N = 106), wing length (R = 0.96, N = 136), egg mass (R = 0.78, N = 46), and egg number (R = 1.00, N = 66).
To evaluate our first hypothesis, we initially tested for statistical dependence in our estimates of the social environment. Using single regressions and adding study population as a repeated measure, andromorph frequency, OSR and male density appeared not significantly related to one another (all P ≥ 0.17 and Spearman R2 ≤ 0.32). This allowed us to simultaneously test the predictive value of these parameters, along with female morph and their morph-specific interactions (see Table 1) in the same mixed ANCOVA models. Four such models were fitted to investigate variation in the fecundity measures, while controlling for spatial and temporal variation, as well as potential differences between mated and single females by adding study population, year and mating status as random variables.
With regards to our second hypothesis, we tested whether both morphs differentially invested in qualitative (egg mass) or quantitative (egg number) reproductive traits. In doing so, we fitted an ANCOVA model with egg number treated as response variable and female morph, egg mass and their interactions as explanatory variables. A significant effect of the interaction would suggest a morph-specific trade-off among both reproductive traits. Meanwhile, we controlled for annual and spatial variation, as well as potential effects of mating status in this analysis by adding year, study population and mating status as random variables to the models. All analyses were performed in SAS 9.2 (SAS Institute Inc, Cary, NC, USA).
With regards to our first hypothesis, both egg number (F1,534 = 4.2; P = 0.04) and clutch mass (F1,515 = 5.6; P = 0.02) varied with population morph frequency and in opposite directions for both female morphs (see Morph*Afreq, Table 2), which provides support for NFDS on female fecundity. Specifically, andromorph females store 8.7% more eggs and have 5.4% higher clutch mass when rare, compared with gynomorphs (Figure 1B,D). This fecundity advantage for andromorphs gradually switched with rising andromorph frequencies towards a reverse situation with higher fecundity for gynomorphs when rare (14.4% and 16.5% difference in egg number and egg mass, respectively). A morph-specific and frequency-dependent effect was not found for egg mass and relative body mass (Table 2; Figure 1A,C). However, all four investigated fecundity estimates significantly decreased with operational sex ratio (P ≤ 0.03; Figure 1E-H), which had similar effects on andromorphs and gynomorphs (see Morph*OSR, Table 2). Male density on the other hand, had no effect on both overall fecundity measures, neither as a main effect (P ≥ 0.40), nor as an interaction with female morph (P ≥ 0.16), see Table 2. As body size may vary with latitude reviewed by  and therefore potentially may influence egg number, we performed separate analyses which also included wing length. However, none of these additional analyses altered the outcome of the analyses presented in Table 2.
Finally, a negative correlation was observed between egg mass and egg number (F1,518 = 10.2, P = 0.002). This reproductive trade-off was observed within all populations, except in Quebec City (see Additional file 2). Interestingly, our results also show a trade-off across populations, in which females from different populations seem to invest more into either egg number or egg mass (Figure 2). However, the resource allocation towards these traits did not differ between female morphs (egg mass*morph: F1,516 = 0.43, P = 0.51; Figure 2). Finally, when controlling for spatial and temporal variation, female morphs did not differ in egg number (F1,537 = 0.08, P = 0.78) or egg mass (F1,518 = 0.99, P = 0.32).
Our study provides some of the first empirical evidence for NFDS on female fecundity in natural conditions and as such provides an important key to understanding the maintenance of intra-sexual phenotypic and genetic variation in this species. Our conclusion was influenced by our observation of a significant inverse relationship between all four of our fecundity measures and relative male density (OSR) as an index of male harassment rate. Indeed, obtrusive males may reduce female fitness in a variety of ways [25, 26, 63] ranging from physical damage to inhibiting foraging success, resulting in suboptimal fecundity; e.g. bees , damselflies . Thus our reported relationships are consistent with the idea that male sexual harassment comes with a fitness costs in females and therefore acts as a major selective force in this study system; see also . Second, we show that the rare female morph has a higher egg load relative to the common one. Clutch mass, as our comprehensive fecundity measure, shows a similar frequency-dependent relationship and is most likely driven by the pattern in egg load. Intriguingly recent work with N. irene, involving exactly the same populations and the same years as the current study, clearly indicated that males prefer to mate with the most common female morph [34, Bots J, Iserbyt A, Hammers M, Van Gossum H, Sherratt TN: Frequency-dependent sexual selection in two intra-specific mimicry systems, in preparation]. All of the above observations collectively support our first hypothesis that positive frequency-dependent male harassment translates into the currently presented negative frequency-dependent patterns in female morph fecundity. Together with the work on polymorphic lizards , diving beetles  and other damselfly species [7, 10] our work emphasises the importance of costly frequency -dependent social interactions as a balancing mechanism to explain intra-specific polymorphisms.
Theoretical and empirical studies indicate that in absence of other mechanisms, NFDS could over generations lead to limited temporal frequency fluctuations along an equilibrium frequency; e.g. fishes , damselflies . Morph frequencies vary little over the years in N. irene (0 - 25%) compared to the very large spatial variation (0 - ~100% , AI unpublished results). In fact, such large spatial variation has now been reported repeatedly in polymorphic damselflies and often resembles a geographical cline (Ischnura elegans, I. senegalensis, N. irene, Megalagrion calliphya). A single equilibrium frequency is thus clearly not reached. This may indicate that besides NFDS operating within populations, additional mechanisms may influence the currently observed population morph frequencies. An obvious suggestion is that divergent selection or gene-by-environment interactions  affect the precise equilibrium, favouring certain morphs under a given set of ecological conditions such as particular densities of con- and heterospecific damselflies, differences in climate, predation rate or parasite load among populations. In addition, historical and present-day stochastic mechanisms have been proposed as well to explain the observed large geographic frequency variation, at least in some parts of a species’ distribution range [66, 67]. Thus, although NFDS on female fecundity appears to be a key balancing mechanism operating within populations, the importance of additional mechanisms are still under debate and should deserve more attention in future research to explain this natural phenomenon thoroughly.
Interestingly, the fecundity advantage for andromorphs is much smaller when rare compared to gynomorphs, indicating an asymmetry in this relationship. This is surprising since it can be expected that andromorphs experience lower harassment rates and associated costs due to their male-like appearance, especially in populations where they are the rare morph [29, 30, 68, 69]. Previous work further indicated that female phenotypic appearance varies with the population morph frequency in this system . Indeed, andromorphs differ overall in body size and shape from gynomorphs [70, 71]. However when common, andromorphs resemble the smaller conspecific male more closely than gynomorphs, which is consistent with theory explaining imperfect mimicry [30, 72]. In addition to direct effects of male harassment (discussed above), it is likely that female fecundity is also shaped by morphological constraints, with smaller or slender females being limited in the number of eggs they can store; see also Ischnura elegans. Thus, although male-like appearance of andromorphs may be beneficial in terms of lower detection rate by harassing males, mimicry may on the other hand come along with a fecundity cost due to allometric associations. The interplay between costs and benefits of mimicry may perhaps explain the asymmetry in our fecundity relationship. Taken all together, present-day patterns in female morph fecundity may relate to direct effects of costly male harassment, perhaps combined with indirect consequences of long-term selection that altered female morph morphology.
Our work also highlights a resource allocation trade-off, in which females invest in either quality (egg mass) or quantity (egg number). This trade-off not only holds within populations, but also differs among populations and may perhaps relate to the inhospitable aquatic environment of the larvae in terms of predation . Odonate larvae are generalist predators that interact aggressively towards conspecifics or heterospecifics, leading to high rates of cannibalism [74–76]. Hatching from heavier eggs, resulting in larger larvae may thus be more advantageous in this (geographically heterogeneous ) competitive environment. Furthermore, we observed increasing egg mass and decreasing egg number towards Northern latitudes. It has been suggested before that body size may increase with latitude due to temperature related physical constraints of growth and development (Bergmann’s rule, reviewed by ), a general pattern which is also confirmed in damselflies [45, 50]. Taken together, investment in heavier eggs may depend on several ecological variables, including abiotic conditions  and degree of the competitive environment , but will most likely come at the expense of producing numerous eggs.
This reproductive trade-off, however, did not differ between morphs, which contrasts with hypothesis two and earlier observations in polymorphic birds  and lizards [15, 38]. These former studies indicated that female morphs may differently allocate resources among life-history traits. The current result also contrasts with recently reported alternative physiological optima in N. irene. Specifically, this latter study showed that andromorphs, compared to gynomorphs, invest more in traits related to immune function and less in flight muscles. In fact, numerous studies with female polymorphic damselflies investigated different life-history components, physiological traits (summarized by ) and morph-specific behavioural strategies to cope with male harassment [33, 80, 81]. Mixed results have been reported previously, but often in favour of the hypothesis that female damselfly morphs represent ARTs in order to escape from excessive male harassment. The combination of traits that differ between the female morphs may be context dependent and/or species specific. Therefore, we suggest that future studies on ARTs focus on an integrated set of behavioural and fitness related traits in different social contexts, and perhaps quantify life-time reproductive success by means of molecular tools .
In conclusion, we provide some of the first empirical evidence collected in natural populations with extreme variation in morph frequencies demonstrating NFDS on female fecundity. As fecundity is a key component of fitness, our results explain an important part of the mechanism maintaining intra-sexual polymorphisms. Frequency-dependent male sexual harassment may well be the driving force of this pattern in our study system, either directly or indirectly affecting egg number and clutch mass of female morphs in a frequency-dependent way. We also show that female morphs do not differently allocate resources into quantitative or qualitative fecundity traits, although this does not exclude the potential for alternative reproductive tactics in this system. As Oliveira et al.  argues, studies on female ARTs are largely understudied and should deserve more attention in future research.
Alternative reproductive tactics
Negative frequency-dependent selection
Operational sex ratio
Relative body mass.
Roulin A: The evolution, maintenance and adaptive function of genetic colour polymorphism in birds. Biol Rev. 2004, 79: 815-848.
Gray SM, McKinnon JS: Linking color polymorphism maintenance and speciation. Trends Ecol Evol. 2007, 22: 71-79.
Forsman A, Ahnesjö J, Caesar S, Karlsson M: A model of ecological and evolutionary consequences of color polymorphism. Ecology. 2008, 89: 34-40.
Ayala FJ, Campbell CA: Frequency-dependent selection. Annu Rev Ecol Syst. 1974, 5: 115-138.
Hori M: Frequency-dependent natural selection in the handedness of scale-eating cichlid fish. Science. 1993, 260: 216-219.
Maynard Smith J: Evolutionary genetics. 1998, New York: Oxford University Press
Takahashi Y, Yoshimura J, Morita S, Watanabe M: Negative frequency-dependent selection in female color polymorphism of a damselfly. Evolution. 2010, 64: 3620-3628.
Oliveira RF, Taborsky M, Brockmann HJ: Alternative reproductive tactics. 2008, Cambridge: Cambridge University Press
Punzalan D, Rodd FH, Hughes KA: Perceptual processes and the maintenance of polymorphism through frequency-dependent predation. Evol Ecol. 2005, 19: 303-320.
Svensson EI, Abbott J, Hardling R: Female polymorphism, frequency dependence, and rapid evolutionary dynamics in natural populations. Am Nat. 2005, 165: 567-576.
Clarke B, Donald PO: Frequency-dependent selection. Heredity. 1962, 19: 201-206.
Sinervo B, Lively CM: The rock-paper-scissors game and the evolution of alternative male strategies. Nature. 1996, 380: 240-243.
Sinervo B, Calsbeek R: The developmental, physiological, neural, and genetical causes and consequences of frequency-dependent selection in the wild. Annu Rev Ecol Syst. 2006, 37: 581-610.
Calsbeek R, Bonvini L, Cox RM: Geographic variation, frequency-dependent selection, and the maintenance of a female-limited polymorphism. Evolution. 2010, 64: 116-125.
Sinervo B, Svensson E, Comendant T: Density cycles and an offspring quantity and quality game driven by natural selection. Nature. 2000, 406: 985-988.
Olendorf R, Rodd FH, Punzalan D, Houde AE, Hurt C, Reznick DN, Hughes KA: Frequency-dependent survival in natural guppy populations. Nature. 2006, 441: 633-636.
Bleay C, Comendant T, Sinervo B: An experimental test of frequency-dependent selection on male mating strategy in the field. Proc Roy Soc B. 2007, 274: 2019-2025.
Gross MR: Alternative reproductive strategies and tactics: diversity within sexes. Trends Ecol Evol. 1996, 11: 92-98.
Härdling R, Bergsten J: Nonrandom mating preserves intrasexual polymorphism and stops population differentiation in sexual conflict. Am Nat. 2006, 167: 401-409.
Kunte K: Female-limited mimetic polymorphism: a review of theories and a critique of sexual selection as balancing selection. Anim Behav. 2009, 78: 1029-1036.
Svensson EI, Abbott JK, Gosden TP, Coreau A: Female polymorphisms, sexual conflict and limits to speciation processes in animals. Evol Ecol. 2009, 23: 93-108.
Cook SE, Vernon JG, Bateson M, Guilford T: Mate choice in the polymorphic African swallowtail butterfly, Papilio dardanus: male-like females may avoid sexual harassment. Anim Behav. 1994, 47: 389-397.
Reinhardt K, Harney E, Naylor R, Gorb S, Siva-Jothy MT: Female-limited polymorphism in the copulatory organ of a traumatically inseminating insect. Am Nat. 2007, 170: 931-935.
Van Gossum H, Sherratt TN, Cordero-Rivera A: The evolution of sex-limited colour polymorphism. Dragonflies: Model Organisms for Ecological and Evolutionary Research. Edited by: Córdoba-Aguilar A. 2008, Oxford: Oxford University Press, 219-229.
Clutton-Brock TH, Parker G a: Sexual coercion in animal societies. Anim Behav. 1995, 49: 1345-1365.
Watson PJ, Arnqvist G, Stallmann RR: Sexual conflict and the energetic costs of mating and mate choice in water striders. Am Nat. 1998, 151: 46-58.
Córdoba-Aguilar A: A female evolutionary response when survival is at risk : male harassment mediates early reallocation of resources to increase egg number and size. Behav Ecol Sociobiol. 2009, 63: 751-763.
Fincke OM: Polymorphic signals of harassed female odonates and the males that learn them support a novel frequency-dependent model. Anim Behav. 2004, 67: 833-845.
Robertson HM: Female dimorphism and mating behaviour in a damselfly, Ischnura ramburi: females mimicking males. Anim Behav. 1985, 33: 805-809.
Sherratt TN: The evolution of female-limited polymorphisms in damselflies: a signal detection model. Ecol Lett. 2001, 4: 22-29.
Gavrilets S, Waxman D: Sympatric speciation by sexual conflict. Proc Natl Acad Sci USA. 2002, 99: 10533-10538.
Van Gossum H, Stoks R, Matthysen E, Valck F, De Bruyn L: Male choice for female colour morphs in Ischnura elegans (Odonata, Coenagrionidae): testing the hypotheses. Anim Behav. 1999, 57: 1229-1232.
Gosden TP, Svensson EI: Density-dependent male mating harassment, female resistance, and male mimicry. Am Nat. 2009, 173: 709-721.
Ting JJ, Bots J, Pérez Jvostov F, Van Gossum H, Sherratt TN: Effects of extreme variation in female morph frequencies on the mating behaviour of male damselflies. Behav Ecol Sociobiol. 2009, 64: 225-236.
Le Galliard J-F, Fitze PS, Ferrière R, Clobert J: Sex ratio bias, male aggression, and population collapse in lizards. Proc Natl Acad Sci USA. 2005, 102: 18231-18236.
Kokko H, Rankin DJ: Lonely hearts or sex in the city? Density-dependent effects in mating systems. Philos Trans Roy Soc B. 2006, 361: 319-334.
Van Gossum H, Sherratt TN: A dynamical model of sexual harassment in damselflies and its implications for female-limited polymorphism. Ecol Model. 2008, 210: 212-220.
Vercken E, Massot M, Sinervo B, Clobert J: Colour variation and alternative reproductive strategies in females of the common lizard Lacerta vivipara. J Evol Biol. 2007, 20: 221-232.
Takahashi Y, Watanabe M: Morph-specific fecundity and egg size in the female-dimorphic damselfly Ischnura senegalensis. Zool Sci. 2010, 27: 325-329.
Iserbyt A, Van Gossum H, Stoks R: Biogeographical survey identifies consistent alternative physiological optima and a minor role for environmental drivers in maintaining a polymorphism. PLoS One. 2012, 7: e32648-
Fincke OM, Jödicke R, Paulson DR, Schultz TD: The evolution and frequency of female color morphs in Holarctic Odonata: why are male-like females typically the minority?. Int J Odonatol. 2005, 8: 183-212.
Sánchez-Guillén R a, Van Gossum H, Cordero Rivera a: Hybridization and the inheritance of female colour polymorphism in two ischnurid damselflies (Odonata: Coenagrionidae). Biol J Linn Soc. 2005, 85: 471-481.
Iserbyt A, Bots J, Ting J, Jvostov FP, Forbes MR, Sherratt TN, Van Gossum H: Multi-annual variation in female morph frequencies of the polymorphic damselfly, Nehalennia irene, at continental and regional scales. Anim Biol. 2009, 59: 313-326.
Gosden TP, Stoks R, Svensson EI: Range limits, large-scale biogeographic variation, and localized evolutionary dynamics in a polymorphic damselfly. Biol J Linn Soc. 2011, 102: 775-785.
Takahashi Y, Morita S, Yoshimura J, Watanabe M: A geographic cline induced by negative frequency-dependent selection A geographic cline induced by negative frequency-dependent selection. BMC Evol Biol. 2011, 11: 256-
Westfall MJ, May ML: Damselflies of North America. 1996, Gainesville: Scientific Publishers
Forbes MRL, Richarson JML, Baker RL: Frequency of female morphs is related to an index of male density in the damselfly, Nehalennia irene (Hagen). Ecoscience. 1995, 2: 28-33.
Van Gossum H, Robb T, Forbes MR, Rasmussen L: Female-limited polymorphism in a widespread damselfly: morph frequencies, male density, and phenotypic similarity of andromorphs to males. Can J Zool. 2008, 86: 1131-1138.
Lam E: Damselflies of the Northeast. A Guide to the Species of Eastern Canada and the Northeastern United States. 2004, New York: Biodiversity Books
Iserbyt A, Bots J, Van Dongen S, Ting JJ, Van Gossum H, Sherratt TN: Frequency-dependent variation in mimetic fidelity in an intraspecific mimicry system. Proc Roy Soc B. 2011, 278: 3116-3122.
Corbet PS: Dragonflies: Behaviour and Ecology of Odonata. 1999, Essex: Harley Books
Takahashi Y, Watanabe M: Female reproductive success is affected by selective male harassment in the damselfly Ischnura senegalensis. Anim Behav. 2010, 79: 211-216.
Dubois B: Damselflies of the North Woods. 2005, Duluth: Kollath-Stensaas Publishing
Van Gossum H, Beirinckx K, Forbes MR, Sherratt TN: Do current hypotheses explain continental and seasonal variation in female morph frequencies of the damselfly, Nehalennia irene?. Biol J Linn Soc. 2007, 90: 501-508.
Bots J, De Bruyn L, Van Dongen S, Smolders R, Van Gossum H: Female polymorphism, condition differences, and variation in male harassment and ambient temperature. Biol J Linn Soc. 2009, 97: 545-554.
Xu M, Fincke OM: Tests of the harassment-reduction function and frequency-dependent maintenance of a female-specific color polymorphism in a damselfly. Behav Ecol Sociobiol. 2011, 65: 1215-1227.
Abramoff MD, Magelhaes PJ, Ram SJ: Image processing with ImageJ. Biophoton Int. 2004, 11: 36-42.
Honěk A: Intraspecific variation in body size and fecundity in insects: a general relationship. Oikos. 1993, 66: 483-492.
Lease HM, Wolf BO: Lipid content of terrestrial arthropods in relation to body size, phylogeny, ontogeny and sex. Physiol Entomol. 2011, 36: 29-38.
Wheeler D: The role of nourishment in oogenesis. Annu Rev Entomol. 1996, 41: 407-431.
Lessells CM, Boag PT: Unrepeatable repeatabilities: a common mistake. Auk. 1987, 104: 116-121.
Chown SL, Gaston KJ: Body size variation in insects: a macroecological perspective. Biol Rev. 2010, 85: 139-169.
Rankin DJ, Dieckmann U, Kokko H: Sexual conflict and the tragedy of the commons. Am Nat. 2011, 177: 780-791.
Rossi BH, Nonacs P, Pitts-Singer TL: Sexual harassment by males reduces female fecundity in the alfalfa leafcutting bee, Megachile rotundata. Anim Behav. 2010, 79: 165-171.
Cooper IA: Ecology of sexual dimorphism and clinal variation of coloration in a damselfly. Am Nat. 2010, 176: 566-572.
Iserbyt A, Bots J, Van Gossum H, Jordaens K: Did historical events shape current geographic variation in morph frequencies of a polymorphic damselfly?. J Zool. 2010, 282: 256-265.
Sánchez-Guillén RA, Hansson B, Wellenreuther M, Svensson EI, Cordero-Rivera A: The influence of stochastic and selective forces in the population divergence of female colour polymorphism in damselflies of the genus Ischnura. Heredity. 2011, 107: 513-522.
Cordero Rivera A, Sánchez-Guillén RA: Male-like females of a damselfly are not preferred by males even if they are the majority morph. Anim Behav. 2007, 74: 247-252.
Van Gossum H, Bots J, Heusden JV, Hammers M, Huyghe K, Morehouse NI: Reflectance spectra and mating patterns support intraspecific mimicry in the colour polymorphic damselfly Ischnura elegans. Evol Ecol. 2011, 25: 139-154.
Abbott JK, Gosden TP: Correlated morphological and colour differences among females of the damselfly Ischnura elegans. Ecol Entomol. 2009, 34: 378-386.
Bots J, Breuker CJ, Van Kerkhove a, Van Dongen S, De Bruyn L, Van Gossum H: Variation in flight morphology in a female polymorphic damselfly: intraspecific, intrasexual, and seasonal differences. Can J Zool. 2009, 87: 86-94.
Sherratt TN: The evolution of imperfect mimicry. Behav Ecol. 2002, 13: 821-826.
Stoks R, Córdoba-Aguilar A: Evolutionary ecology of Odonata: a complex life cycle perspective. Ann Rev Entomol. 2012, 57: 249-265.
Anholt BR: Cannibalism and early instar survival in a larval damselfly. Oecologia. 1994, 99: 60-65.
Padeffke T, Suhling F: Temporal priority and intra-guild predation in temporary waters: an experimental study using Namibian desert dragonflies. Ecol Entomol. 2003, 28: 340-347.
Schenk K, Söndgerath D: Influence of egg size differences within egg clutches on larval parameters in nine libellulid species (Odonata). Ecol Entomol. 2005, 30: 456-463.
McPeek MA: Determination of species composition in the Enallagma damselfly assemblages of permanent lakes. Ecology. 1990, 71: 83-98.
Rollinson N, Hutchings JA: Environmental quality predicts optimal offspring size in the wild. Am Nat. 2013, in press
Bots J: The maintenance of female colour polymorphism in a damselfly - PhD thesis. 2008, Antwerp: Universiteit Antwerpen
Van Gossum H, Stoks R, De Bruyn L: Frequency-dependent male mate harassment and intra-specific variation in its avoidance by females of the damselfly Ischnura elegans. Behav Ecol Sociobiol. 2001, 51: 69-75.
Iserbyt A, Van Gossum H: Unexpected absence of behavioural differences between female damselfly colour morphs. Anim Behav. 2009, 78: 1463-1469.
Thompson DJ, Hassall C, Lowe CD, Watts PC: Field estimates of reproductive success in a model insect: behavioural surrogates are poor predictors of fitness. Ecol Lett. 2011, 14: 905-913.
We would like to express our gratitude to J. Ting, F.P. Jvostov and M. Op De Beeck for their assistance in data collection, to S. Lindgren and the Zwozdesky family for local support in BC and the staff of Queens University Biological Station in ON. Three anonymous referees provided constructive comments on an earlier version of the manuscript. This research was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) to TNS, the Institute for the Promotion of Innovation through Science and Technology in Flanders (IWT-Vlaanderen research grant to JB and AI), the Fund for Scientific Research-Flanders (postdoctoral FWO research grant to HVG and JB) and a research grant from the University of Antwerp (BOF-NOI, project number: 2401 to HVG).
The authors declare that they have no competing interests.
All authors participated in designing the study, performing field work, providing financial support and drafting this manuscript. In addition, AI performed all dissections, statistically analysed the dataset and created all graphs. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: sw ), key estimates of the social environment and female morph fecundity.(DOCX 29 KB)
About this article
Cite this article
Iserbyt, A., Bots, J., Van Gossum, H. et al. Negative frequency-dependent selection or alternative reproductive tactics: maintenance of female polymorphism in natural populations. BMC Evol Biol 13, 139 (2013). https://doi.org/10.1186/1471-2148-13-139
- Alternative reproductive tactics
- Colour polymorphism
- Male harassment
- Quantity-quality trade-off
- Sexual conflict