BACKGROUND
evolutionary biologists have long studied visible polymorphisms as they are excellent model systems to examine micro evolutionary processes  <cit> .
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 arises when individuals of a rare morph experiences a higher fitness than those of a more common type  <cit> .
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  <cit> .
classic examples of nfds include coexistence of different colour morphs, to gain access to mates  <cit> , to challenge predators  <cit>  and to lower sexual conflict intensity  <cit> .
although the idea of nfds has been appreciated for decades  <cit> , 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  <cit> .
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  <cit>  and theoretical  <cit>  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 .
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  <cit> .
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  <cit> , diving beetles  <cit> , african bat bugs  <cit> , damselflies  <cit> .
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  <cit> .
mate searching males are considered to face fewer cognitive challenges when confronted with only one, rather than multiple female phenotypes coexisting within populations  <cit> , some of which may appear like males .
therefore on a phenotypic level, females may experience diversifying sexual selection to avoid sexual harassment  <cit> .
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  <cit> , increased male sexual interest has been observed towards the most common female morph in a given population  <cit> .
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  <cit> .
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 , because male–female interactions occur more frequently under these conditions  <cit> .
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  <cit> , which forms the first aim of the current study.
in addition, sex-limited polymorphisms are frequently considered as alternative reproductive tactics , i.e.
a discontinuous set of selected traits to maximize reproductive output in two or more alternative ways  <cit> .
although repeatedly studied in males, recent observations in several species of owl  <cit> , lizard  <cit>  and insect  <cit>  collectively indicate that female morphs may also represent arts.
for example in the polymorphic lizard uta stansburiana <cit> , combined density-dependent and negative frequency-dependent interactions among conspecifics determine the relative success of orange  and yellow  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  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  <cit>  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  <cit> .
testing the second hypothesis, we explore whether female morphs exhibit arts in which resources are allocated differently into qualitative  or quantitative  fecundity traits  <cit>  and meanwhile account for potential spatial and temporal variation.
methods
model system
female polymorphism has been observed in more than a hundred damselfly and dragonfly species  <cit> .
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  <cit> .
large geographic variation in frequencies and densities of males and female morphs has been described in several species  <cit> , 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 .
nehalennia irene is not an endangered nor a protected species  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  <cit>  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  <cit> , gynomorph females have distinctive yellowish lateral thorax sides and a less conspicuous abdominal melanin pattern; for colour figures see  <cit> , for pictures see  <cit> .
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  <cit> .
this tandem formation can last several hours .
if receptive, the female cooperates by bending her abdomen towards the male’s secondary genitalia  to form a ‘copulation wheel’  <cit> .
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  <cit> .
in n. irene, females lay eggs in floating pieces of dead plant material while the female is tandem-guarded by her last successful male  <cit> .
several clutches of eggs are laid throughout a female’s lifetime  <cit> .
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  <dig> to >90% throughout the species’ distribution range over canada  <cit> .
for our current aims we selected six study populations that differed significantly in social conditions .
frequency and density estimates were obtained in a manner similar to that described in van gossum et al.
<cit> .
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 , operational sex ratio  and male density  were calculated and collectively quantify the social environment.
calculating these parameters including or excluding immature individuals gave very similar results, see also  <cit> .
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  <cit> .
each of the six populations were monitored during the reproductive season over three consecutive years .
estimates of andromorph frequency , operational sex ratio  and male density  per population and province , averaged across three study years .
an extended table with detailed information on sample sizes per population and per year is given in .
sample collection was carried out on several days throughout the reproductive season  always between 9 am and 3 pm.
this is the period before the majority of females start to oviposit .
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  <dig> adult andromorphs and  <dig> gynomorphs in each population for each investigated year.
measurements for relative body mass  were performed in all three years  and the three other fecundity estimates were investigated in two successive years .
for more detailed sample sizes per population and per year, .
mating status at the moment of capture was noted, i.e.
being single or mating .
all individuals were stored for further measurements in 95% ethanol immediately after capture.
fecundity estimates
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  <dig> balance .
next, a digital picture was taken  of the right hind wing.
using imagej  <dig> x  <cit> , wing length was measured from the second antenodal cross vein to the stigma; see  <cit>  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 .
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  <cit> .
hence, rbm is here treated as a coarse measure of overall fecundity.
dissection of specimens was performed under a leica mz  <dig>  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  at 60°c for 12 h and weighted on a sartorius se2f balance .
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  <cit> .
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 .
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  <cit> .
a limited measurement error was observed: body mass , wing length , egg mass , and egg number .
statistical analyses
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 .
this allowed us to simultaneously test the predictive value of these parameters, along with female morph and their morph-specific interactions  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  or quantitative  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  <dig>  .
RESULTS
with regards to our first hypothesis, both egg number  and clutch mass  varied with population morph frequency and in opposite directions for both female morphs , which provides support for nfds on female fecundity.
specifically, andromorph females store  <dig> % more eggs and have  <dig> % higher clutch mass when rare, compared with gynomorphs .
this fecundity advantage for andromorphs gradually switched with rising andromorph frequencies towards a reverse situation with higher fecundity for gynomorphs when rare .
a morph-specific and frequency-dependent effect was not found for egg mass and relative body mass .
however, all four investigated fecundity estimates significantly decreased with operational sex ratio , which had similar effects on andromorphs and gynomorphs .
male density on the other hand, had no effect on both overall fecundity measures, neither as a main effect , nor as an interaction with female morph , see table  <dig>  as body size may vary with latitude reviewed by  <cit>  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  <dig>
explanatory variables include female morph frequencies , male density  and operational sex ratio .
all four analyses are controlled for variation among populations, years and mating status .
the asterisks  indicate tested interactions between two main effects.
numerator degrees of freedom is  <dig> in all cases, df in the table refers to the denominator degrees of freedom.
note that egg mass, egg number and relative clutch mass were measured in two successive years and body mass in three subsequent years.
significant results are indicated in bold.
finally, a negative correlation was observed between egg mass and egg number .
this reproductive trade-off was observed within all populations, except in quebec city .
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 .
however, the resource allocation towards these traits did not differ between female morphs .
finally, when controlling for spatial and temporal variation, female morphs did not differ in egg number  or egg mass .
discussion
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  as an index of male harassment rate.
indeed, obtrusive males may reduce female fitness in a variety of ways  <cit>  ranging from physical damage to inhibiting foraging success, resulting in suboptimal fecundity; e.g.
bees  <cit> , damselflies  <cit> .
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  <cit> .
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 .
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  <cit> , diving beetles  <cit>  and other damselfly species  <cit>  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  <cit> , damselflies  <cit> .
morph frequencies vary little over the years in n. irene  compared to the very large spatial variation .
in fact, such large spatial variation has now been reported repeatedly in polymorphic damselflies and often resembles a geographical cline .
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  <cit>  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  <cit> .
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  <cit> .
previous work further indicated that female phenotypic appearance varies with the population morph frequency in this system  <cit> .
indeed, andromorphs differ overall in body size and shape from gynomorphs  <cit> .
however when common, andromorphs resemble the smaller conspecific male more closely than gynomorphs, which is consistent with theory explaining imperfect mimicry  <cit> .
in addition to direct effects of male harassment , 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 <cit> .
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  or quantity .
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  <cit> .
odonate larvae are generalist predators that interact aggressively towards conspecifics or heterospecifics, leading to high rates of cannibalism  <cit> .
hatching from heavier eggs, resulting in larger larvae may thus be more advantageous in this  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 , a general pattern which is also confirmed in damselflies  <cit> .
taken together, investment in heavier eggs may depend on several ecological variables, including abiotic conditions  <cit>  and degree of the competitive environment  <cit> , 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  <cit>  and lizards  <cit> .
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 <cit> .
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  and morph-specific behavioural strategies to cope with male harassment  <cit> .
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  <cit> .
CONCLUSIONS
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.
<cit>  argues, studies on female arts are largely understudied and should deserve more attention in future research.
abbreviations
afreq: andromorph frequency; arts: alternative reproductive tactics; mdens: male density; nfds: negative frequency-dependent selection; osr: operational sex ratio; rbm: relative body mass.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
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.
supplementary material
additional file 1
detailed overview of sample dates , key estimates of the social environment and female morph fecundity.
click here for file
additional file 2
quantity-quality trade-off within each study population.
click here for file
acknowledgements
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  to tns, the institute for the promotion of innovation through science and technology in flanders , the fund for scientific research-flanders  and a research grant from the university of antwerp .

 evolutionary biologists have long studied visible polymorphisms as they are excellent model systems to examine micro evolutionary processes. 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 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. 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, 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. 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 and theoretical 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. 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. 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. 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. 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. 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. 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, because male–female interactions occur more frequently under these conditions. 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, 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 and insect 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 and yellow 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 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 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. testing the second hypothesis, we explore whether female morphs exhibit arts in which resources are allocated differently into qualitative or quantitative fecundity traits and meanwhile account for potential spatial and temporal variation. methods model system 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, 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. nehalennia irene is not an endangered nor a protected species 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, 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. if receptive, the female cooperates by bending her abdomen towards the male’s secondary genitalia 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 to >90% throughout the species’ distribution range over canada. for our current aims we selected six study populations that differed significantly in social conditions. 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, operational sex ratio and male density 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. each of the six populations were monitored during the reproductive season over three consecutive years. estimates of andromorph frequency, operational sex ratio and male density per population and province, averaged across three study years. an extended table with detailed information on sample sizes per population and per year is given in. sample collection was carried out on several days throughout the reproductive season always between 9 am and 3 pm. this is the period before the majority of females start to oviposit. 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 adult andromorphs and gynomorphs in each population for each investigated year. measurements for relative body mass were performed in all three years and the three other fecundity estimates were investigated in two successive years. for more detailed sample sizes per population and per year,. mating status at the moment of capture was noted, i.e. being single or mating. all individuals were stored for further measurements in 95% ethanol immediately after capture. fecundity estimates 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 balance. next, a digital picture was taken of the right hind wing. using imagej x, 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. 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 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 at 60°c for 12 h and weighted on a sartorius se2f balance. 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. 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. 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, wing length, egg mass, and egg number. statistical analyses 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. this allowed us to simultaneously test the predictive value of these parameters, along with female morph and their morph-specific interactions 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 or quantitative 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.  with regards to our first hypothesis, both egg number and clutch mass varied with population morph frequency and in opposite directions for both female morphs, which provides support for nfds on female fecundity. specifically, andromorph females store % more eggs and have % higher clutch mass when rare, compared with gynomorphs. this fecundity advantage for andromorphs gradually switched with rising andromorph frequencies towards a reverse situation with higher fecundity for gynomorphs when rare. a morph-specific and frequency-dependent effect was not found for egg mass and relative body mass. however, all four investigated fecundity estimates significantly decreased with operational sex ratio, which had similar effects on andromorphs and gynomorphs. male density on the other hand, had no effect on both overall fecundity measures, neither as a main effect, nor as an interaction with female morph, see table  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  explanatory variables include female morph frequencies, male density and operational sex ratio. all four analyses are controlled for variation among populations, years and mating status. the asterisks indicate tested interactions between two main effects. numerator degrees of freedom is in all cases, df in the table refers to the denominator degrees of freedom. note that egg mass, egg number and relative clutch mass were measured in two successive years and body mass in three subsequent years. significant results are indicated in bold. finally, a negative correlation was observed between egg mass and egg number. this reproductive trade-off was observed within all populations, except in quebec city. 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. however, the resource allocation towards these traits did not differ between female morphs. finally, when controlling for spatial and temporal variation, female morphs did not differ in egg number or egg mass. discussion 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 as an index of male harassment rate. indeed, obtrusive males may reduce female fitness in a variety of ways 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. 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 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 compared to the very large spatial variation. in fact, such large spatial variation has now been reported repeatedly in polymorphic damselflies and often resembles a geographical cline. 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. 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. 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. however when common, andromorphs resemble the smaller conspecific male more closely than gynomorphs, which is consistent with theory explaining imperfect mimicry. in addition to direct effects of male harassment, 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 or quantity. 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. hatching from heavier eggs, resulting in larger larvae may thus be more advantageous in this 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, a general pattern which is also confirmed in damselflies. 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. 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 and morph-specific behavioural strategies to cope with male harassment. 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. abbreviations afreq: andromorph frequency; arts: alternative reproductive tactics; mdens: male density; nfds: negative frequency-dependent selection; osr: operational sex ratio; rbm: relative body mass. competing interests the authors declare that they have no competing interests. authors’ contributions 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. supplementary material additional file 1 detailed overview of sample dates, key estimates of the social environment and female morph fecundity. click here for file additional file 2 quantity-quality trade-off within each study population. click here for file acknowledgements 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 to tns, the institute for the promotion of innovation through science and technology in flanders, the fund for scientific research-flanders and a research grant from the university of antwerp.
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. the coexistence of multiple male and/or female morphs is widely believed to be maintained through negative frequency-dependent selection imposed by social interactions. 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. our results further indicate that all investigated fecundity traits are negatively affected by relative male density, which confirms male harassment as selective agent. alternative reproductive tacticscolour polymorphismfitnessmale harassmentodonataquantity-quality trade-offsexual conflict 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. sex-limited polymorphisms have long intrigued evolutionary biologists and have been the subject of long-standing debates. 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. second, we evaluate whether the two female morphs differ in resource allocation among fecundity traits, indicating alternative tactics to maximize reproductive output. this morph-specific fecundity advantage gradually switched with the population morph frequency.
 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, which confirms male harassment as selective agent. finally, we show a clear trade-off between qualitative and quantitative 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. alternative reproductive tacticscolour polymorphismfitnessmale harassmentodonataquantity-quality trade-offsexual conflict 