sexual systems in angiosperms range from hermaphroditism ( monomorphic populations of plants with bisexual flowers ) to dioecy ( dimorphic populations of male and female individuals ) and include almost all imaginable combinations and gradations ( table 1 ; [ 1 , 2 ] and references therein ) . such remarkable diversity of sexual systems has perplexed naturalists and evolutionary biologists for a long time [ 36 ] . the evolution of dioecy from a hermaphroditic ancestor has been particularly difficult to understand because the invasion and maintenance of unisexual mutants in a population of hermaphrodites require that the loss of fitness resulting from the loss of one sexual function be compensated by increased fitness gains through the remaining sexual function of the unisexual mutant [ 7 , 8 ] .
this requirement seems very restrictive , and therefore considerable effort has been devoted towards understanding the conditions under which dioecy can evolve [ 5 , 818 ] . in contrast , the evolution of sexually dimorphic traits following the evolution of dioecy ( successful establishment of only two reciprocal unisexual morphs in a population ) has received less attention . consequently , our current understanding of the evolution of sex - related traits ultimately leading to morphological or physiological differences between unisexual morphs ( i.e. , sexual dimorphism or secondary sexual traits ) is still limited , despite recent advances and excellent syntheses on the topic [ 5 , 6 , 1921 ] .
this paper focuses on one set of traits subject to becoming sexually dimorphic upon the evolution of dioecy : those traits that provide plants with defence against herbivores .  
sex - biased herbivory may be one of the selective pressures conducive to the evolution of dioecy , and it can also be a consequence of sex - specific selection on patterns of resource allocation in dioecious species . considering only the gynodioecy pathway of the evolution of dioecy , we can think of three possible scenarios regarding the role of herbivory in each of the two steps involved in this pathway ( figure 1 ) .
the first step in the gynodioecy pathway to dioecy is the successful establishment of females ( male - sterile mutants ) in a population of hermaphrodites , thus resulting in a gynodioecious population .
as mentioned above , this step requires that females compensate for the fitness loss incurred with the loss of the male function .
the reallocation of resources freed from the male function towards defence may contribute towards fitness compensation if increased defence results in greater fitness for the females .
increased defence may result in lower herbivore damage on females than on hermaphrodites ( figure 1 , path b ) .
defence may be achieved through resistance : traits that reduce the rate of herbivore attack such as low nutritional content of tissues ( particularly , n content ) , secondary metabolites , trichomes , cutin , waxy cuticles , lignin , and volatiles that attract natural enemies of herbivores ; and also through tolerance : traits that mitigate the negative effects of damage on fitness , including higher or lower growth rates , mobilization of stored resources , and activation of apical meristems .
if females reallocate resources to tolerance traits , they could be the morph with greater herbivore damage ( figure 1 , path c ) .
the second step in the gynodioecy pathway to dioecy is the successful establishment of male individuals ( female - sterile mutants ) in a gynodioecious population followed by the loss of the hermaphroditic morph , thus resulting in a dioecious population . upon the evolution of two unisexual morphs
, defensive traits may evolve differently in each sex and eventually become sex linked [ 8486 ] .
the particular way in which defensive traits diverge between sexes will depend on the costs and benefits derived from the specific pattern of resource allocation to growth , reproduction , and defence in each sex .
currently , it is thought that females generally evolve greater resistance than males ( see the following ; figure 1 , path b ) .  
therefore , we will not delve into morph - biased herbivory in gynodioecious species , which would be the topic of a different essay .
however , we do recognize that sex - biased herbivory  indicative of sexual dimorphism in resistance against herbivores  is likely related to morph - biased herbivory in the ancestral gynodioecious population from which it evolved ( figure 1 , path b - b ) .
thus , male - biased herbivory may have its origin in a gynodioecious population where hermaphrodites ( functionally , the male morph ) bear greater levels of herbivory than females .  
the above view for the origin of greater resistance against herbivores in females is based directly on the principle of allocation : resources freed from the male function are used for the female function , growth , and defence .
in contrast to this view , the finding of male - biased herbivory in dioecious populations has been explained on the basis of sex - specific selection of resistance traits , where the main difference between sexes that drives the sex - specific selection is the cost of reproduction . in this alternative view
, female individuals of dioecious species are expected to have lower herbivory levels than males because the higher cost of reproduction of females confers a selective advantage to females with traits that reduce herbivore attack .
the logic of this argument is as follows : since females invest more in reproduction than males , they are left with a smaller pool of resources for growth and therefore must grow more slowly than males [ 20 , 87 ] .
according to the resource availability hypothesis , the fitness cost of losing tissue to herbivores is greater for plants that grow more slowly , thus favouring the evolution of increased defence against herbivores in slow - growing plants . since females tend to grow more slowly , they should be better defended against herbivores than males .
we must note that the argument for greater defence levels in females is usually understood in terms of resistance , but it could also be interpreted in terms of tolerance , in which case the predictions of sex - biased herbivory would be the opposite ( figure 1 , path c ) , as developed below .
the first review of the empirical studies on the topic of sex - biased herbivory concluded that 
males are more likely than females to be preferentially used by herbivores  and suggested that male - biased herbivory was widespread among dioecious species .
the authors , however , recognized that sex - biased herbivory was by no means a unanimous finding across all the dioecious species examined to that date , and that the relative susceptibility of each sex to herbivory could be influenced , among other factors , by fluctuations in ecological tradeoffs between functions ( rather than evolutionary changes in patterns of allocation ) , such as phenological changes in resource allocation to reproduction and growth [ 75 , 89 ] .
therefore , the life stage at which damage measurements are taken can determine whether a study concludes that herbivory is sex biased or not .  
in addition , gren et al . cautioned against publication bias , whereby studies that found differences between genders could be more likely to be published than those that did not ; and taxonomic bias , an overabundance of studies from certain genera or families . in this instance , the taxonomic bias is correlated with an ecological bias for studies of temperate species , despite dioecy being more prevalent in tropical ecosystems .
. called for future studies that ( 1 ) examine the causes of differential palatability between males and females , ( 2 ) measure the fitness consequences of natural levels of herbivory in both sexes , and ( 3 ) determine whether herbivore pressure can actually cause adaptive changes in tissue palatability .
in addition , they urged for broadening the taxonomical scope of the studies . in spite of such encouragement , there still is a paucity of studies that address these issues .
more recently , sex - biased herbivory in dioecious species was tested by means of a meta - analysis of 33 studies encompassing 30 species , 19 of which were previously included in gren et al .
however , they did not emphasize other shortcomings of the dataset and concluded that male - biased herbivory in dioecious species is a rule .
here we propose alternative evolutionary scenarios that could result in female - biased herbivory or lack of intersexual differences in herbivory levels .
we invite the reader to reconsider the evidence for male - biased herbivory in dioecious plants and recommend a standard protocol for evolutionary - ecological studies of sex - biased herbivory in dioecious species that addresses the shortcomings listed below .
we contend that taking for granted the generality of male - biased herbivory in dioecious species is hampering our progress in this field .
  while male - biased herbivory has been explained as a consequence of sex - specific selection of resistance determined by the cost of reproduction of each sex , few of the reviewed studies ( table 2 ) are actually placed within such evolutionary context .
j. lovett - doust and l. lovett - doust   were the first to argue , citing charnov , that an evolutionary divergence between sexes in resource allocation patterns could result in sex - biased resistance against herbivory .
based their expectation of male - biased herbivory on sexual selection : it may be more advantageous for males to invest less in resistance and more in reproduction when exposed to a cyclic herbivore compared to a noncyclic herbivore because males will lose fitness only once every four years , and their reproductive output during years of little or no herbivory will more than compensate for the fitness lost to herbivores on the heavy herbivory years .
the above explanation would hold only to the extent that there are no carry - over effects from one year to another .  
more recently , mccall   cited bierzychudek and eckhart   and delph   in support of the claim that the reproductive output of females is more limited by resources than that of males .
however , while there may be ample evidence that reproductive allocation is generally greater for the female function , it is possible that , upon the separation of sexes , physiological mechanisms involved in resource acquisition and allocation evolve in such way as to minimize the differences in reproductive effort between sexes .
for example , in ilex glabra , sexes do not differ in total reproductive biomass produced in a growing season because the greater unitary investment in pistillate flowers and fruit development is negated by the sevenfold greater flower production in males .
given the importance of the tenet of greater female reproductive allocation for the expectation of male - biased herbivory , all studies of sex - biased herbivory should test for intersexual differences in reproductive allocation or provide a reference to an empirical study that demonstrates such differences for the species in question . in the measurement of reproductive allocation , particular attention
must be paid to obtaining reliable estimates of male reproductive output ( pollen production ) , which presents its own logistical difficulties .
in addition , resource expenditure in pollinator attraction needs to be considered , as this is another expenditure related to reproduction that may differ between sexes .  
in essence , the presumed chain of evolutionary events that lead to male - biased herbivory in dioecious plants stems from the reallocation of those resources freed upon the loss of a sexual function in unisexual mutants towards defence . in most studies
however , defence may also occur through tolerance [ 22 , 93 ] . in comparison to resistance traits ,
tolerance traits have been more elusive . the capacity to store and mobilize carbohydrates , the presence of meristems and the capacity to activate them in response to damage
growth rate has also been proposed to influence tolerance , although it is controversial whether high or low growth rates favour tolerance [ 9498 ] .
a recent model shows that plants with low growth rates are more tolerant to herbivore damage .
this model also shows that plants that can change their growth rate positively in response to damage will tolerate damage better than those with a different response in growth rate .
the activation of meristems and mobilization of resources in response to the loss of tissue are two well - documented responses [ 100 , 101 ] that could contribute to increased plant growth rates in response to damage .
thus , according to this model , if females grow more slowly than males because of their greater allocation to reproduction , then females should be more tolerant to herbivory than males .
whether this prediction necessarily implies that females should be less resistant than males is not clear at this point since there is increasing evidence that these two modes of defence are not necessarily mutually exclusive [ 99 , 103 , 104 ] .  
there is one other possibility that has not been emphasized enough in the proposed models of the evolution of defence in dioecious species : while one possible consequence of a greater allocation of resources to reproduction in females is reduced allocation to growth , it is also possible that the main reduction in allocation is to defence . in this case , there would be no detectable detriment to growth .
consequently , female plants would suffer more damage ( if they are less resistant ; figure 1 , path c ) or greater fitness losses ( if they are less tolerant ) compared to males .
greater damage in females could lead to fewer , more spaced reproductive events or greater interannual variability in reproductive output either directly through a decrease in the availability of resources for reproduction or indirectly through a decrease in pollinator visitation rates due to a lack of resources needed for floral display or nectar production [ 75 , 105108 ] .
fewer resources available for reproduction could also pose a selective pressure to become choosier about their mates , which may lead to increased fruit or seed abortion [ 75 , 109 , 110 ] .   
evolutionary changes in the rate of resource acquisition in female individuals may occur through increased photosynthetic rates , canopy area , rates of mineral nutrient uptake , as well as greater branching of roots , and enhancement of mycorrhizal associations .
a greater rate of resource acquisition in females would decrease the relative differences in costs of reproduction between sexes : the sex with the greater resource demand for reproduction would have an increased capacity to garner resources . under such scenario , the life - time cost of reproduction at the individual level would be equal between sexes , thus eliminating the source of inequalities in the patterns of allocation between males and females . in summary
, nothing dictates that there is only one evolutionary pathway regarding changes in the patterns of resource allocation among reproduction , growth , and defence following the evolution of unisexuality ( figure 1 ) .  
alternatively , a stage in which female individuals have heavier damage levels because of resource limitation for resistance may be a transient evolutionary stage prior to the invasion of mutants whose greater defence levels are attained at the expense of growth . in this case , we should observe female - biased herbivory in younger dioecious lineages and male - biased herbivory in those lineages in which there has been enough time for selection to reshape the patterns of resource allocation to reproduction , defence , and growth
. we should be able to test this by means of relative dating of lineages with male- or female - biased herbivory .
similarly , as long as there has not been selection on prereproductive growth rates following the evolution of unisexuality , we should not see differences in growth rates or other physiological vegetative traits between males and females before their first reproductive event .
it is difficult to test this prediction without reliable morphological or genetic markers that allow juveniles to be sexed so that their performance can be compared on the basis of sex .
some sex - linked markers may be , effectively , sex - related traits expressed before the onset of reproduction . whether the presence of these markers implies the existence of sex chromosomes is still an area in need of further investigation [ 84 , 85 , 111113 ] .    in short , without fitness gain curves for each sex , it is difficult to predict accurately which sex should evolve greater resistance against herbivores and whether we should expect or not male - biased herbivory in dioecious species .
in fact , we need fitness surfaces in order to include the effects of reductions in leaf area caused by herbivory . moreover
, the fitness surfaces should account for the short- and long - term responses of plants in terms of changes in photosynthetic rates , reallocation of resources to shoot or root , activation of meristems , and delays in phenology or shortening of life span brought about by herbivore damage [ 95 , 99 ] .   
it has not escaped our attention that the evolution of defence in gynodioecious species can be approached from a similar perspective to the one presented above for dioecious species .
it is important to consider that some of gynodioecious species may be in evolutionary transition to dioecy while others are not [ 112 , 114 ] .
another important difference with dioecious species is that , in gynodioecious plants , the morph that performs the male function  the hermaphrodites  may have a greater cost of reproduction because of the expenditure of resources on two sexual functions .
does this mean that hermaphrodites would be less resistant ( figure 1 , path c ) or grow more slowly and evolve greater resistance ( figure 1 , path b ) ? clearly , making predictions with respect to gender dimorphism in defensive traits for bisexual conditions along the gradation from monoecy to dioecy is not straight forward .
the collection of studies cited in the reviews of herbivory in dioecious species [ 72 , 73 ] has , as a group , important shortcomings that weaken the conclusion of male - biased herbivory as a generality in dioecious species .
the main shortcomings are ( 1 ) the taxonomic bias of the sample of species studied ; and ( 2 ) failure to test for or provide references of empirical studies of intersexual differences in ( a ) resistance traits  the purported cause of the intersexual differences in herbivory levels ; ( b ) growth rates  the purported cause of the intersexual differences in resistance to herbivore attack ; and ( c ) reproductive effort  the purported cause of the aforementioned intersexual differences in growth rates .
these deficiencies had been pointed out earlier [ 72 , 89 ] , but judging by statements included in the introduction or discussion of many of the papers published after the 1999 review , such caveats have not been considered to their full extent , and many authors take for granted either the generality of male - biased herbivory in dioecious species or its expectation without reference to any theoretical context .  
cornelissen and stiling 's meta - analysis of sex - biased herbivory includes 30 species , 28 of which are angiosperms . focusing only on angiosperms , 13 of the 28 species were not considered previously in gren et al .
these 30 species represent a total of 20 genera , 18 families , and 10 orders .
nine of those species belong to the same genus : salix . adding to those populus tremula
, the species in the salicaceae represent one - third of all species considered for the meta - analysis .
such distribution contrasts greatly with the taxonomic distribution of dioecy in 14,620 species , 959 genera , 157 families , and 36 orders . of the four dioecious genera ( consisting of solely dioecious species ) with most species ( 400 ) ,
pandanus , diospyros , and litsea , with 700 , 500 + , and 400 species , respectively (; s. renner , university of munich , unpublished data ) have not been studied for sex - biased herbivory yet .
clearly , we need to direct our research efforts to the most understudied orders and families if we want to arrive at generalizations regarding the biology of dioecious species , and particularly the influence of herbivores in their ecology and evolution .
in addition to this taxonomic bias , a critical reexamination of that list of species casts serious doubt on the conclusion that male - biased herbivory is a rule in dioecious species : only 13 of those species were reported invariably to have male - biased herbivory .
this list includes three salix , two freycinetia , and two species for which evidence of male - biased herbivory has not been published : hippophae rhamnoides and rumex acetosa .
( in fact , male - biased herbivory in myrica gale  not included in these 13 species  is also anecdotal . )
greater herbivore damage on females is reported for four species , while the rest show either no intersexual differences ( 16 ) , or variation in the result , depending on different factors ( species of herbivore , kind of herbivore , tissue damaged , time of year , phenological or ontogenetic stage , etc . ) .
moreover , it is possible that the results for those 13 species would show variation with either population or site , had these factors been studied . perhaps the most serious problem with several studies of herbivory and dioecy has been the failure to make the connection between sex - biased herbivore damage and intersexual differences in growth rate , precisely because the latter is the purported cause of the former . of the 30 species of angiosperms in table 2 , either growth rate or a surrogate variable for growth ( e.g. , shoot length ) was measured only in 21 species .
males grew faster in six species , females in two , no difference between sexes was detected in six species , and three species showed variable results .
it must be noted that the same number of species shows no difference between genders as those that show greater growth rate in males .
considering solely the 13 species that invariably showed male - biased herbivory , only two show a greater growth rate in males : acer negundo and hippophae rhamnoides .
however , as the evidence for male - biased herbivory in h. rhamnoides is anecdotal , we are left with only one species for which growth rate was measured in the same study as herbivore damage : acer negundo .  
only 12 of the 30 species listed were assessed for intersexual differences in reproductive allocation in terms of reproductive effort ( the proportion of biomass or other currency devoted to reproductive structures relative to the total biomass or expenditure in the selected currency of an individual ) .
reproductive effort was greater in females of 10 species and in males for the other two species .
in some species , reproductive effort was measured during flowering , but allocation to fruit production was not considered ( e.g. , silene dioica ) . in those cases ,
we are left with an incomplete picture of reproductive allocation , and we can only join the authors in speculating whether species of the same genus have similar patterns of reproductive allocation .
the only species that have been assessed for foliar damage , growth rate , and reproductive allocation in the same study are c. tepejilote , b. halimifolia , i. glabra , n. psychotrioides , and r. alpinus ( table 2 ) .
these studies clearly made the chain of causal connections from sex bias in reproductive allocation all the way to sex bias in some resistance traits ( except for c. tepejilote and i. glabra ) , and , as a consequence of the latter , sex bias in levels of damage .
the study of r. alpinus went even further , comparing these attributes between pre- and postreproductive plants , and thus emphasizing that the root of the differences in growth rates , resistance traits , and leaf damage is in the patterns of reproductive resource expenditures [ 63 , 116 ] . in some species , reproductive allocation , growth rate , and/or resistance
however , even with these studies , the number of species for which we have a more complete picture of the causal links amongst these attributes remains low : nine more species ( c. alternans , b. dracunculifolia .
a. canescens , r. acetosella , s. caprea , s. cinerea , s. lasiolepis , s. sericeae , and h. rhamnoides [ ? ] ; table 2 ) now have published data for damage and growth rate , bringing the number of species in this situation to 15 .
two more species , for a total of three , now have data on damage , growth rate , and reproductive allocation ( c. alternans , s. dioica , and n. psychotrioides )
. two more species now have data on reproductive allocation , growth rate , and resistance apart from herbivore damage , for a total of four species with all four variables measured ( a. triphyllum , b. halimifolia , l. benzoin , and r. alpinus ) . in summary
, the majority of studies on the topic of sex - biased herbivory have neglected the purported causal connections between bias in reproductive allocation , differential growth rate , resistance , and herbivore damage .
also , some authors seemed to confuse theoretical expectations with empirical evidence of greater female reproductive allocation : while lloyd and webb   argue convincingly for the expectation of greater reproductive effort in females , they provided empirical evidence only for rumex acetosella , citing putwain and harper 1972 . therefore , lloyd and webb 's excellent paper can not be cited as solid empirical evidence of greater reproductive effort in females .
lastly , anecdotal evidence should be taken with great caution and always flagged as such until data are published ( see entries marked  ?  in table 2 ) .   using the search terms herbiv * and dioec * for entries between january 1998 and may 2012 on the web of science
, we found nine studies encompassing 14 species that were not included in either of the previous reviews of the topic .
of these , only the study on the three species of chamaedorea palms measured reproductive allocation , growth rate , resistance , and herbivore damage ( table 3 ; n.b .
: one of these species had been studied before : c. alternans = c. tepejilote ) .
only one other study measured damage and reproductive allocation ( sclerocarya birrea , table 3 ) .
similarly , growth rate was assessed in only one other species ( salix arctica ) .
the taxonomic breadth of the studies of herbivory in dioecious species increased only by one family ( in an unplaced order of the euasterids i ) .
the general lack of consistency in the level of detail and the variables that have been measured in all these studies could be addressed if researchers interested in this topic followed a minimally standardized protocol .
  new studies must clearly allude to the theoretical framework from which the prediction of sex - biased herbivory levels ( resistance ) stems  resource allocation theory , in particular , sex allocation .
the claim that male - biased herbivory is expected because it has been reported as a pattern , whether implicit or explicit , lacks heuristic value because it does not address the causes of such pattern .
moreover , a plethora of factors may modify the expected pattern , as shown above .  
clearly , we need to increase the taxonomic breadth of the studies of herbivory in dioecious species .
we could direct our attention to those families with the greatest number of dioecious species or those with the greatest proportion of dioecious species .
the first alternative will miss families with low species richness that may have a high proportion of dioecious species .
the second method will miss families with high species richness but low proportion of dioecious species .
one possible compromise is to focus our studies on the families with the greatest number of dioecious species among those with a large proportion of dioecious species , for instance , 50% or more ( table 4 ) .
so far , we have studied only 2% of the dioecious species in the most studied family ( salicaceae ) .
if we took that as a target , we would have to study about 155 species for the 30 most dioecious species - rich families of angiosperms .
however , by this method we would include only one species per family for many families , thus failing to achieve adequate representation of those families .
in addition , we must consider that the conditions that determine sex - biased herbivory can change with habitat , and therefore some species may need to be studied in several habitats .   
in addition to the taxonomic bias , there is a preponderance of studies of woody plants .
while this is understandable because most dioecious species are woody , we should strive for representation of the herbaceous component .
with increased research on herbaceous dioecious species , we can address the influence of life history traits on the evolution of dioecy and defence .  
lastly , we propose that all studies aimed at assessing whether herbivory levels differ between sexes and whether these differences are a consequence of differential growth rates ( in turn resulting from differential allocation to reproduction ) should conduct , at least , the following measurements and observations : ( 1 ) levels of herbivory , measured as precisely as possible ( preferably for more than one growing season in perennials ) ; ( 2 ) species of herbivores responsible for most of the damage ; ( 3 ) growth rates , measured either as rgr for whole individuals or from increments in branch length or leaf production ; ( 4 ) reproductive allocation , measured both as the number of reproductive structures ( flowers and pollen production for males , flowers , and fruits for females ) , and also as reproductive effort ( the proportion of individual or shoot biomass allocated to reproduction , and when possible n and p allocation to reproductive structures ) ; and finally ( 5 ) the most important resistance characters that could be influencing the levels of herbivory and measure them quantitatively .
in addition , these studies could add an experimental component in which plants are damaged at least at the highest rate seen in the surveys of natural damage , so as to measure tolerance to herbivory as well as resistance [ 75 , 102 ] .
ideally , damage should be performed by placing natural herbivores on the plants because mechanical damage does not necessarily elicit the same physiological responses as herbivore damage .
also , these studies should consider that resistance and tolerance may vary both with ontogeny and with respect to other reproductive phenology because the acquisition and expenditure of resources vary at different stages of development and life history [ 119121 ] .
we must reiterate that other authors have emphasized the need to address several of the points outlined above .
it is our hope that future studies take these recommendations seriously so that we have to assume and speculate less , and we have empirical data to further our understanding of the evolution of defence in dioecious species .
  the study of the evolution of sex - biased herbivory is hampered by the notion that male - biased herbivory in dioecious species is a rule .
we have presented other possible evolutionary outcomes with regards to sex - biased herbivory in the transition from hermaphroditic populations to dioecious ones .
we have also discussed how these different outcomes can be predicted under different theoretical assumptions .
therefore , future studies of herbivory in dioecious species should be based on a clear theoretical framework . in particular , we urge that all new studies of herbivory in dioecious species include assessments of reproductive allocation , growth rates , and resistance traits deemed to differ between sexes and , therefore , determine sex - biased herbivory .
in addition , tolerance should also be considered as a potentially important defence mode that can vary between sexes . in this manner
the advancement of our knowledge about sex - related defence in plants should help us gain a better understanding of the evolution of sex - related traits in general .
 defensive traits may evolve differently between sexes in dioecious plant species . 
 our current understanding of this process hinges on a partial view of the evolution of resistance traits that may result in male - biased herbivory in dioecious populations . here , we present a critical summary of the current state of the knowledge of herbivory in dioecious species and propose alternative evolutionary scenarios that have been neglected . 
 these scenarios consider the potential evolutionary and functional determinants of sexual dimorphism in patterns of resource allocation to reproduction , growth , and defence . 
 we review the evidence upon which two previous reviews of sex - biased herbivory have concluded that male - biased herbivory is a rule for dioecious species , and we caution readers about a series of shortcomings of many of these studies . 
 lastly , we propose a minimal standard protocol that should be followed in any studies that intend to elucidate the ( co)evolution of interactions between dioecious plants and their herbivores . 