- Research article
- Open Access
Unravelling the role of host plant expansion in the diversification of a Neotropical butterfly genus
© The Author(s). 2016
- Received: 23 March 2016
- Accepted: 7 June 2016
- Published: 16 June 2016
Understanding the processes underlying diversification is a central question in evolutionary biology. For butterflies, access to new host plants provides opportunities for adaptive speciation. On the one hand, locally abundant host species can generate ecologically significant selection pressure. But a diversity of host plant species within the geographic range of each population and/or species might also eliminate any advantage conferred by specialization. This paper focuses on four Melinaea species, which are oligophagous on the family Solanaceae: M. menophilus, M. satevis, M. marsaeus, and finally, M. mothone. We examined both female preference and larval performance on two host plant species that commonly occur in this butterfly’s native range, Juanulloa parasitica and Trianaea speciosa, to determine whether the different Melinaea species show evidence of local adaptation.
In choice experiments, M. mothone females used both host plants for oviposition, whereas all other species used J. parasitica almost exclusively. In no choice experiment, M. mothone was the only species that readily accepted T. speciosa as a larval host plant. Larval survival was highest on J. parasitica (82.0 % vs. 60.9 %) and development took longer on T. speciosa (14.12 days vs. 13.35 days), except for M. mothone, which did equally well on both host plants. For all species, average pupal weight was highest on J. parasitica (450.66 mg vs. 420.01 mg), although this difference was least apparent in M. mothone.
We did not find that coexisting species of Melinaea partition host plant resources as expected if speciation is primarily driven by host plant divergence. Although M. mothone shows evidence of local adaptation to a novel host plant, T. speciosa, which co-occurs, it does not preferentially lay more eggs on or perform better on this host plant than on host plants used by other Melinaea species and not present in its distributional range. It is likely that diversification in this genus is driven by co-occurring Müllerian mimics and the resulting predation pressure, although this is also likely made possible by greater niche diversity as a consequence of plasticity for potential hosts.
- Adaptive divergence
- Aposematic mimicry
- Diet breadth
- Ecological speciation
- Host shift
- Niche width
- Oviposition preference
- Resource use
Understanding the processes underlying diversification is a central question in evolutionary biology. In butterflies, access to new host plants provides opportunities for adaptive speciation, and is thought to be the primary driving mechanism for the diversification of phytophagous insects . Indeed, there exists a positive correlation between host plant diversity and butterfly species diversity in general . Janz et al.  postulated that the diversification of phytophagous insects has been driven by oscillations in host plant ranges, which can lead to increased species distributions and subsequently opportunities for secondary specialization on novel hosts (see also ). Although multiple-host use implies increased plasticity, which should counterbalance local adaptation, Janz & Nylin  suggested that host expansions can also lead to other factors causing divergent selection to increase population fragmentation and for diversification to occur.
As such, locally abundant host plant species might generate ecologically significant selection pressures, and adaptive evolution in response to locally abundant plant species has been shown to occur, even in the face of high gene flow. For example, Eastern tiger swallowtail butterflies (Papilio glaucus) showed greater oviposition preference and larval performance on host plants that were more abundant at their collection sites . However, a diversity of host plant species within the geographic range of each population and/or species can also eliminate any advantages conferred by specialization. For example, Ladner & Altizer  found that although host plant influenced oviposition and larval performance in the Monarch (Danaus plexippus), they could find no evidence for local adaptation to hosts found within each of the butterfly population’s breeding range. Indeed, Chew  noted that species that lay eggs singly and visit many different host species might be slow to evolve consistent preferences.
The genus Melinaea belongs to the Neotropical tribe Ithomiini (Nymphalidae: Danainae), which are aposematic butterflies extensively involved in Müllerian mimicry rings. Butterfly species that are part of a Müllerian mimicry ring possess the same warning signal, and effectively share the cost of educating predators of their unpalatability. In the genus Melinaea, species often consist of multiple subspecies characterised by different wing colour patterns which are associated with distinct mimetic communities. The systematics for this genus remains unclear. Previous studies using mitochondrial and nuclear genes, and rapidly evolving microsatellite markers, show little genetic differentiation among taxa, which suggests recent speciation with incomplete lineage sorting, possibly combined with recent gene flow [8–11]. One possible explanation for this rapid diversification is radiation across different host plants. Records of host plant usage suggest that Melinaea are oligophagous on the family Solanaceae, specifically the widespread tribe Juanulloeae, which might be the reason for their large geographic distribution . However, different populations of different species are likely to encounter only a subset of potential hosts, as no single host plant overlaps the entire geographic range. For example, in Ecuador, M. menophilus has been shown to use both Juanulloa ochracea (K. Willmott, pers. com.) and Juanulloa mexicana  in the lowlands, but Markea sp.  in the montane. Because many environmental variables are correlated with elevation, adaptation and divergent selection have been shown to occur along altitudinal ranges (e.g., ).
Studying host plant choice in the genus Melinaea may provide insights into how adaptation can influence the diversification of mimetic butterflies. For example, detailed studies of heliconiine communities demonstrate that coexisting species of Heliconius partition host plant resources , as expected if speciation is driven by host plant divergence. However, speciation in Heliconius has also been shown to be possible without host shifts , and competitive exclusion on larval hosts may instead be driving the observed patterns of distribution. The Melinaea species in north-eastern Peru (San Martín and Loreto department) are of particular interest for such studies, as multiple species, many consisting of different subspecies, are present and overlap in distribution (see ).
Here we examine whether the different Melinaea species show evidence of local adaptation to host plant species that commonly occur in their native range. Some herbivores are known to feed on plants for which larval performance (i.e., survivorship and development time) is suboptimal, raising the possibility that other extrinsic factors, such as competition or enemy-free space, may contribute to insect fitness and influence the evolution of diet breadth . As such, we look at both female preference and larval performance. We predict that exposure to different host plant species in their respective distribution ranges could select for divergent host use traits, so that butterflies should preferentially lay more eggs on, and larvae should perform better on, species common to their native habitats. We also attempt to determine whether the directionality of larval performance corresponded with female oviposition preference.
This paper focuses on four Melinaea species: M. menophilus (which consists of the two subspecies M. menophilus ssp. n 1 and M. m. hicetas, found in transitional forest and lowland forest respectively), M. mothone, which is found in high altitudinal habitats, M. marsaeus (consisting of the two subspecies M. marsaeus phasiana and M. m. rileyi, found in transitional forest and lowland forest respectively), and M. satevis cydon, a lowland species. In a previous study , we extensively surveyed host plant use in different habitats and found that all the different Melinaea species in this study utilize Juanulloa parasitica (Solanaceae) as a host plant, except for M. mothone, which uses the novel host plant Trianaea speciosa (Solanaceae). Indeed, because J. parasitica did not co-occur with M. mothone, this butterfly species was never collected on this host plant. Similarly, because T. speciosa is only present at higher altitudes, it is present throughout M. mothone’s distribution, but only co-occurs with M. menophilus ssp. n 1 at the very edge of its distribution, and does not co-occur with any of the other Melinaea species in this study. Furthermore, the only species ever collected from T. speciosa is M. mothone, even when M. menophilus was present at low densities. As such, there appears to be an ecological difference related to altitude and host plant use, in addition to genetic differences , between M. mothone and the other Melinaea species of this study. Because the different subspecies of both M. menophilus and M. marsaeus share similar habitat preference, host plant use, and have overlapping distribution , comparisons were only made between the species rather than between the subspecies (i.e., subspecies were pooled).
Gravid wild caught female butterflies were collected in north-eastern Peru (San Martín and Loreto department), in 2012–2013. Collection localities consisted of transitional forest habitats surrounding Tarapoto (Rio Shilcayo basin: 6°27’30”S 76°21’00”W alt 460 m), Carachamera (6°25’85”S 76°15’27”W alt 280 m) and Shapaja (6°36’56”S 76°09’61”W alt 195 m); high altitudinal habitat in the Cordillera Escalera near Tarapoto (the Tunel ridge: 6°27’11”S 76°17’11”W alt 1090 m) and in Moyobamba (6°04’34”S 76°57’27”W alt 1130 m); and lowland forest on Pongo-Baranquita road (6°17’53”S 76°14’38”W alt 200 m) and Shucushyacu (5°57’20”S 75°53’06”W alt 183 m). Because each species has a limited distribution, butterflies of a given species were only collected at a limited number of sites, often in close proximity to each other. Furthermore, previous studies have failed to find any geographical population structure . As such, butterflies from different collection sites were pooled for the analysis. Females were kept in 2 × 2 × 2 m outdoor insectaries under ambient conditions in Tarapoto, San Martín, and were provided with sugar water solution and bee pollen for nourishment.
To test for female oviposition preference across different hosts, we performed two types of trials, as follows: 1) choice experiments, in which a single mated female butterfly was placed in a cage containing both a potted J. parasitica and a T. speciosa, and 2) no-choice experiments, in which females were placed in cages with two T. speciosa plants. For all trials, care was taken so that the plants were of similar size, and that the plants did not touch. No-choice oviposition experiments examined whether female butterflies can accept the novel host plant T. speciosa when no alternatives are present. Strict preference of normal hosts in a choice situation does not preclude the possibility of accepting a novel host when no alternative is present. All trials lasted for four days and 18 females per species were tested for each treatment.
The number of eggs laid on each of the plants was recorded at the same time daily and larvae were collected shortly after hatching. Larvae collected from at least eight different unrelated females used in the preference experiment (so as to reduce any maternal or family effects) per species were used for tests of larval performance (N = 32/species/host plant). Larvae were reared individually in transparent plastic containers in the shade behind a nearby building under ambient conditions. Leaves of either J. parasitica or of T. speciosa were offered ad libitum, and larvae were checked daily for food replacement and cleaning, as well as survival. Duration of larval development was recorded and pupae were weighed the morning following pupation using a small portable scale.
A generalized linear model (GzLM) was done using SPSS (SPSS Inc., Chicago, Illinois) to determine if the number of eggs laid on either J. parasitica or T. speciosa in the choice experiment was affected by either the Melinaea species tested or the host plant, and to test for a possible interaction between these two factors. Because the number of eggs for most of the species included many zero counts, a negative binomial distribution was used. The total number of eggs laid on both T. speciosa in the no-choice experiment was compared between Melinaea species using a non-parametric Kruskal-Wallis test. Larval survival was compared between species for each host plant using a multivariate Cox proportional hazards model. Duration of larval development and pupal weight were compared between species and host plants using a 2-way ANOVA to test for an overall effect of treatment.
Results of a multivariate Cox proportional hazards models for larval mortality of four Melinaea species (M. menophilus, M. mothone, M. marsaeus, and M. satevis) reared on either Juanulloa parasitica or Trianaea speciosa (N = 32 larvae/Melinaea species/host plant)
(95 % confidence interval)
Most death occurred during the early stages of larval development and it occurred more frequently in larvae reared on T. speciosa, except for M. mothone, where it occurred equally in both groups. Mean age at death was of 9.9 days for all species reared on J. parasitica. However, all species other than M. mothone reared on T. speciosa died at a much younger age (6.7 days vs. 9.3 days for larvae of M. mothone reared on T. speciosa).
Determining which traits and corresponding selective pressures initiate divergence is key to understanding the causes of ecological speciation, and diversification in phytophagous insects has often been thought to be driven by shifting and adapting to new host plants . Although the genus Melinaea is characterized by a recent and rapid diversification across much of the Neotropics [8–11], we found that most species in this study co-occur with and utilize the same host plant, J. parasitica, suggesting that diversification and speciation in this genus has mostly occurred without changes in host plants. We did, however, find large differences in performance on the novel host plant T. speciosa between the co-occurring M. mothone and the other species, consistent with local adaptation to host plant species.
But the fact that M. mothone has retained the ability to use and do equally well on both J. parasitica and T. speciosa, despite the fact that the former does not co-occur, and in addition to the absence of oviposition preference for the co-occurring host species, suggests that the use of a novel host plant is more likely due to niche expansion as a consequence of plasticity for potential hosts. We also found a positive relationship between oviposition preference and larval performance for the other Melinaea species, but the reduction in larval performance on the novel host plant Trianaea was not as pronounced as the decline in oviposition. Evidence for asymmetry between oviposition preference and larval performance on novel hosts has been documented for other herbivores, as a result of, for example, plant chemistry, competition or enemy-free space, host abundance, larval conditioning during development [19–21]. It has been postulated that oviposition on novel hosts that are suitable for larval development, be it a result of oviposition mistakes or the result of a labile oviposition strategy, might enable host range expansions [21–23]. Shifts by herbivorous insects are sometimes restricted to related plant species, but also occur on unrelated plants, and can be mediated by chemical similarities of the new host, or can be explained by patterns of parallel cladogenesis, and/or increased ecological opportunities (i.e., hosts that are geographically available) (see  and references therein). Our results suggest that some potential to oviposit, feed, and survive on the related but novel host plant Trianaea is already present in the genus Melinaea, which is consistent with the use of multiple hosts by this group of butterflies throughout their geographic range (see ).
Furthermore, a key prediction of ecological speciation is that traits that prevent gene flow from eroding adaptation are likely those that evolved early, either directly as the result of adaptive divergence, or as the result of selection preventing the formation of maladaptive hybrids [25–27]. In the genus Melinaea, McClure & Elias  observed strong assortative mating among taxa, possibly as a result of homogamy for colour pattern, even in very recently diverged lineages such as the M. marsaeus subspecies, and even in the absence of host plant differentiation. In this genus, it is likely that plasticity in host use has enabled for a large geographical distribution, which in turn likely exposed populations to a different suite of potential Müllerian mimics. For example, Melinaea mothone has distinct elevational preferences and is associated with typically Andean mimicry complexes, including another Melinaea species which shares the same colour pattern, M. isocomma. Because predation pressure differs for different mimetic communities , spatially segregated populations of Melinaea butterflies are likely selected to harbour different colour patterns that coincide with those patterns that are most common within their given mimetic environment [29, 30].
Through the maintenance of a spatial mosaic of mimetic colour patterns, predation on Müllerian mimics constrains geographical distribution and allows for different species or subspecies, even those with similar ecological niches, to exist in different regions . Migrants between populations suffer reduced survival because they are non-mimetic outside their habitat and suffer higher levels of predation attacks, which can directly reduce gene flow between populations by lowering the rate of heterospecific encounters . Switches in mimicry can also lead to pleiotropic changes in mate choice, as assortative mating often coevolves with colour pattern, and reinforcement against maladaptive non-mimetic hybrids [33, 34].
Over time, accumulated differences of other ecological aspects, including but not limited to host plant use, will ultimately accumulate, leading to reproductive isolation and speciation. Nevertheless, the expansion of the potential host plant repertoire, rather than host shifts per se, may have been an important driver of diversification, because more potential host plants means a larger area of distribution and, as such, a larger number of potential niches [3, 4, 35].
Although access to new host plants can provide opportunities for ecological speciation, we did not find that coexisting species of Melinaea partition host plant resources as expected if speciation is primarily driven by host plant divergence. Furthermore, although M. mothone shows evidence of local adaptation to a novel host plant, T. speciosa, which co-occurs, it does not preferentially lay more eggs on or perform better on this host plant than on host plants used by other Melinaea species and not present in its distributional range. Rather, it is likely that diversification in this genus is driven by co-occurring Müllerian mimics and the resulting predation pressure, although this is also likely made possible by greater niche diversity as a consequence of plasticity for potential hosts. To understand the causes and consequences of evolution in ecological traits, more studies are needed of groups in which diversification is recent or ongoing, and for which multiple ecological traits are well-described.
We thank the Peruvian authorities and Dr. Gerardo Lamas (Museo de Historia Natural, Univesidad Mayor de San Marcos) for research permits (236-2012-AG-DGFFS-DGEFFS and 201-2013-MINAGRI-DGFFS/DGEFFS), and Mathieu Chouteau and Ronald Mori-Pezo for their precious help in the field.
This research was funded by le Fonds Québecois de la Recherche sur la Nature et les Technologies as a PDF award to MM, and by a CNRS ATIP grant and ANR grant (SPECREP) awarded to ME.
Availability of data and material
The data set supporting the results of this article are available in the Dryad repository, doi:10.5061/dryad.qt104.
MM and ME conceived the study. MM performed all the experiments and analysed the data. MM and ME drafted the manuscript. Both authors contributed in the discussion, writing and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
No specific approval was needed to conduct this study.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ehrlich PR, Raven PH. Butterflies and plants: a study in coevolution. Evolution. 1964;18:586–608.View ArticleGoogle Scholar
- Janz N, Nylin S, Wahlberg N. Diversity begets diversity: host expansions and the diversification of plant-feeding insects. BMC Evol Biol. 2006;6:4.View ArticlePubMedPubMed CentralGoogle Scholar
- Nylin S, Slove J, Niklas J. Host plant utilization, host range oscillations and diversification in Nymphalid butterflies: a phylogenetic investigation. Evolution. 2014;68:105–24.View ArticlePubMedGoogle Scholar
- Janz N, Nylin S. The oscillation hypothesis of host plant-range and speciation. In: Tilmon KJ, editor. Specialization, Speciation, and Radiation: the Evolutionary Biology of Herbivorous Insects. Berkeley: University of California Press; 2008. p. 203–15.Google Scholar
- Bossart JL, Scriber JM. Maintenance of ecologically significant genetic variation in the tiger swallowtail butterfly through differential selection and gene flow. Evolution. 1995;49:1163–71.View ArticleGoogle Scholar
- Ladner DT, Altizer S. Oviposition preference and larval performance of North American monarch butterflies on four Asclepias species. Entomol Exp Appl. 2005;116:9–20.View ArticleGoogle Scholar
- Chew F. Coevolution of pierid butterflies and their cruciferous foodplants. I. The relative quality of available resources. Oecologia. 1975;29:117–27.View ArticleGoogle Scholar
- Whinnett A, Zimmermann M, Willmott KR, Herrera N, Mallarino R, Simpson F, Joron M, Lamas G, Mallet J. Strikingly variable divergence times inferred across an Amazonian butterfly 'suture zone'. Proc R Soc Lond Ser B Biol Sci. 2005;272:2525–33.View ArticleGoogle Scholar
- Elias M, Hill RI, Willmott KR, Dasmahapatra KK, Brower AVZ, Mallet J, Jiggins CD. Limited performance of DNA barcoding in a diverse community of tropical butterflies. Proc R Soc Lond Ser B Biol Sci. 2007;274:2881–9.View ArticleGoogle Scholar
- Dasmahapatra KK, Lamas G, Simpson F, Mallet J. The anatomy of a "suture zone" in Amazonian butterflies: a coalescent-based test for vicariant geographic divergence and speciation. Mol Ecol. 2010;19:4283–301.View ArticlePubMedGoogle Scholar
- McClure M, Elias M. Ecology, life history, and genetic differentiation in Neotropical Melinaea (Nymphalidae: ithomiini) butterflies from north-eastern Peru. Zool J Linn Soc. 2016. In press. doi:10.1111/zoj.12433.
- Slove J, Janz N. The relationship between diet breadth and geographic range size in the butterfly subfamily Nymphalinae - a study of global scale. PLoS One. 2011;6:e16057.View ArticlePubMedPubMed CentralGoogle Scholar
- Drummond III BA, Brown KS. Ithomiinae (Lepidoptera: Nymphalidae): Summary of known larval food plants. Ann Mo Bot Gard. 1987;74:341–58.View ArticleGoogle Scholar
- Willmott KR, Mallet J. Correlations between adult mimicry and larval host plants in ithomiine butterflies. Proc R Soc Lond Ser B Biol Sci. 2004;271:S266–9.View ArticleGoogle Scholar
- Funk WC, Murphy MA, Hoke KL, Muths E, Amburgey SM, Lemmon EM, Lemmon AR. Elevational speciation in action? Restricted gene flow associated with adaptive divergence across an altitudinal gradient. J Evol Biol. 2016;29:241–52.View ArticlePubMedGoogle Scholar
- Gilbert LE. Biodiversity of a Central American Heliconius community: pattern, process, and problems. In: Price PW, Lewinsohn TM, Fernandes GW, Benson WW, editors. Plant-animal interactions: evolutionary ecology in tropical and temperate regions. New York: Wiley; 1991. p. 403–27.Google Scholar
- Jiggins CD, McMillan WO, Mallet J. Host plant adaptation has not played a role in the recent speciation of Heliconius himera and Heliconius erato. Ecol Entomol. 1997;22:361–5.View ArticleGoogle Scholar
- Vasconcellos-Neto J. Interactions between ithomiine butterflies and Solanaceae: feeding and reproductive strategies. In: Price PW, Lewinsohn TM, Fernandes GW, Benson WW, editors. Plant-animal interactions: evolutionary ecology in tropical and temperate regions. New York: John Wiley and Sons, Inc; 1991. p. 291–313.Google Scholar
- Gratton C, Welter SC. Oviposition preference and larval performance of Liriomyza helianthi (Diptera: Agromyzidae) on normal and novel host plants. Environ Entomol. 1998;27:926–35.View ArticleGoogle Scholar
- Gómez Jiménez MI, Sarmiento CE, Díaz MF, Chautá A, Peraza A, Ramírez A, et al. Oviposition, larval preference, and larval performance in two polyphagous species: does the larva know best? Entomol Exp Appl. 2014;153:24–33.View ArticleGoogle Scholar
- Murphy SM. Inconsistent use of host plants by the Alaskan swallowtail butterfly: adult preference experiments suggest labile oviposition strategy. Ecol Entomol. 2007;32:143–52.View ArticleGoogle Scholar
- Courtney SP, Kibota TT. Mother doesn't know best: selection of hosts by ovipositing insects. In: Bernays EA, editor. Insect-plant interactions, vol. 2. Boca Raton: CRC; 1990. p. 162–88.Google Scholar
- Kibota TT, Courtney SP. Jack of one trade, master of none: host choice by Drosophila magnaquinaria. Oecologia. 1991;86:251–60.View ArticleGoogle Scholar
- Becerra JX, Venable DL. Macroevolution of insect-plant associations: the relevance of host biogeography to host affiliation. Proc Natl Acad Sci U S A. 1999;96:12626–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Matsubayashi KW, Ohshima I, Nosil P. Ecological speciation in phytophagous insects. Entomol Exp Appl. 2010;134:1–27.View ArticleGoogle Scholar
- McBride CS, Singer MC. Field studies reveal strong postmating isolation between ecologically divergent butterfly populations. PLoS Biol. 2010;8:e1000529. doi: 1000510.1001371/journal.pbio.1000529.
- Sandoval CP, Nosil P. Counterating selective regimes and host preference evolution in ecotypes of two species of walking-sticks. Evolution. 2005;59:2405–13.View ArticlePubMedGoogle Scholar
- Chouteau M, Arias M, Joron M. Warning signals are under positive frequency-dependent selection in nature. PNAS. 2016. doi: 10.1073/pnas.1519216113.
- Chazot N, Willmott K, Santacruz Endera PG, Toporov A, Hill RI, Jiggins CD, Elias M. Filtering by elevation and mutualistic mimicry shape the structure of Andean butterfly communities. Am Nat. 2014;183:26–39.View ArticlePubMedGoogle Scholar
- Elias M, Gompert Z, Jiggins C, Willmott K. Mutualistic interactions drive ecological niche convergence in a diverse butterfly community. PLoS Biol. 2008;6:2642–9.View ArticlePubMedGoogle Scholar
- Sherratt TN. Spatial mosaic formation through frequency-dependent selection in Müllerian mimicry complexes. J Theor Biol. 2006;240:165–74.View ArticlePubMedGoogle Scholar
- Nosil P. Reproductive isolation caused by visual predation on migrants between divergent environments. Proc R Soc B. 2004;271:1521–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Jiggins CD, Naisbit RE, Coe RL, Mallet J. Reproductive isolation caused by colour pattern mimicry. Nature. 2001;411:302–5.View ArticlePubMedGoogle Scholar
- Naisbit RE, Jiggins CD, Mallet J. Disruptive sexual selection against hybrids contributes to speciation between Heliconius cydno and H. melpomene. Proc R Soc Lond Ser B Biol Sci. 2001;268:1849–54.View ArticleGoogle Scholar
- Nylin S, Janz N. Ecology and evolution of host plant range: butterflies as a model group. In: Olff H, Brown VK, Drent RH, editors. Herbivores: Between Plants and Predators. Oxford: Blackwell; 1999. p. 31–54.Google Scholar