Avoid, attack or do both? Behavioral and physiological adaptations in natural enemies faced with novel hosts
© Vacher et al; licensee BioMed Central Ltd. 2005
Received: 04 February 2005
Accepted: 04 November 2005
Published: 04 November 2005
Confronted with well-defended, novel hosts, should an enemy invest in avoidance of these hosts (behavioral adaptation), neutralization of the defensive innovation (physiological adaptation) or both? Although simultaneous investment in both adaptations may first appear to be redundant, several empirical studies have suggested a reinforcement of physiological resistance to host defenses with additional avoidance behaviors. To explain this paradox, we develop a mathematical model describing the joint evolution of behavioral and physiological adaptations on the part of natural enemies to their host defenses. Our specific goals are (i) to derive the conditions that may favor the simultaneous investment in avoidance and physiological resistance and (ii) to study the factors that govern the relative investment in each adaptation mode.
Our results show that (i) a simultaneous investment may be optimal if the fitness costs of the adaptive traits are accelerating and the probability of encountering defended hosts is low. When (i) holds, we find that (ii) the more that defended hosts are rare and/or spatially aggregated, the more behavioral adaptation is favored.
Despite their interference, physiological resistance to host defensive innovations and avoidance of these same defenses are two strategies in which it may be optimal for an enemy to invest in simultaneously. The relative allocation to each strategy greatly depends on host spatial structure. We discuss the implications of our findings for the management of invasive plant species and the management of pest resistance to new crop protectants or varieties.
In natural antagonistic systems such as host-parasite, plant-herbivore, and predator-prey systems (hereafter called 'host-enemy'), enemies may frequently be confronted with hosts expressing novel defenses. For instance, in agro-ecosystems, herbivorous insects are confronted with novel plant defenses each time a new chemical pesticide or each time a new toxic cultivar is introduced. This situation is frequent in natural areas too, when an herbivorous insect's foraging area is invaded by a novel toxic plant variety or species (e.g. invasion of Rocky mountain meadows by Thlaspi arvense, a crucifer lethal to Pieris larvae ). 'Avoid' and 'attack' are two basic strategies that enemies may evolve to cope with these novel defenses. In other words, enemies may evolve the ability to discriminate between defended and undefended hosts and preferentially avoid defended ones (behavioral adaptation), or develop a direct counter-adaptation allowing the successful attack of defended hosts (physiological adaptation). Both adaptation modes have been extensively studied and reviewed, especially in arthropod systems [2–5].
Although simultaneous investment in both adaptations may first appear to be redundant (since avoiding defended hosts is unnecessary if defenses have been overcome anyway), some empirical studies have suggested a reinforcement of physiological resistance to host defenses with additional avoidance behaviors. The most convincing example is perhaps the study of Pluthero & Threlkeld (1981) . These authors measured the levels of behavioral avoidance and physiological resistance in eight strains of wild-caught Drosophila melanogaster tested for their responses to the insecticide malathion. They showed that the most resistant line had also the highest degree of avoidance, suggesting a reinforcement of physiological resistance with an additional avoidance behavior in a wild population. However they did not find any significant correlation between these two modes of insecticide resistance in the eight strains, which indicates that the mechanisms involved were genetically independent from each other. The independence between the two adaptation modes is also suggested by other studies that have demonstrated the evolution of either physiological resistance or behavioral avoidance, but not both. For instance, in Plutella xylostella, exposure to transgenic plants expressing toxins from Bacillus thuringiensis (Bt) induced physiological resistance to Bt toxins without discrimination between transgenic and non-transgenic plants , whereas exposure to toxic baits altered the behavior of German cockroaches, but not their physiological resistance .
To explain the paradox that apparently redundant adaptations may evolve, a number of authors [8–10] have employed population genetics models. These models assume that physiological and behavioral responses are governed by two independent loci each bearing two co-dominant alleles. Nine strategy sets (i.e. nine genotypes) corresponding to the combination of three physiological adaptation levels (high, mild or null) and three behavioral adaptation levels (high, mild or null) are assumed. The results highlight the impact of population-genetic and population-dynamic factors on behavioral and physiological adaptations. In particular, a mixed strategy may be stable depending on the relative fitnesses of the nine genotypes and the initial allelic frequencies [8, 9] and the mode of population regulation can have a striking impact on the likelihood of behavioral adaptation to evolve . Simulations by Rausher  reveal that under the "hard selection" mode (i.e. regulatory factors act on the population as a whole ) the pure behavioral strategy evolves almost as frequently as the pure physiological strategy, whereas under "soft selection" (i.e. the subpopulations of each host are regulated independently ) the pure behavioral strategy never evolves. A potential explanation for this difference is that contrary to "hard selection", "soft selection" leads to overcrowding on the most suitable host and undercrowding on the other host . Behavioral generalists therefore have a higher fitness than choosy enemies, because the former suffer lower levels of intra-specific competition. A less frequent outcome (c. 10% of the runs) is the evolution of a mixed strategy (i.e. a strategy where one of the loci at least is polymorphic). Under "hard selection", this outcome is favored by the absence of fitness costs to both traits .
The objective of this study is to (i) explore in more detail the conditions favoring the occurrence of a mixed strategy and (ii) study the factors that govern the relative investment into each adaptation mode. Specifically, we study the impact of the shape of the fitness cost functions and the impact of spatial heterogeneity in host defenses. Indeed, while several theoretical models have investigated the impact of the abundance and spatial distribution of suitable hosts on physiological resistance evolution [e.g. [12–17]] or on the evolution of host selection behavior [e.g. [18–21]], the impact of space on the joint evolution of physiological and behavioral adaptations to host defenses has been largely ignored . In our model, physiological resistance and behavioral avoidance are represented as quantitative traits and we adopt an optimality approach to identify the conditions favoring the investments in both traits. Contrary to the population genetics models listed above, this approach permits quantitative predictions of the optimal relative investments in both forms of adaptation. Fitness costs of the adaptations and host spatial structure are assumed to be constant.
Below we show that (i) simultaneous investment may be optimal if the fitness costs of the adaptive traits are accelerating and the probability of encountering defended hosts is low. When (i) holds, we find that (ii) the more that defended hosts are rare and/or spatially aggregated, the more behavioral adaptation is favored.
What conditions favor the simultaneous investment in physiological and behavioral adaptations to host defenses?
Accelerating fitness costs of adaptations
Rare and/or spatially aggregated defended hosts
What factors govern the relative investment in each adaptation mode?
Spatial configuration of hosts
Magnitude of search costs
In agreement with previous studies [8–10], we find that the simultaneous investment in avoidance and physiological resistance can be an optimal strategy despite interference between both adaptations (i.e., investment in one adaptation mode decreases the efficiency of investment in the other). We identify two conditions that must be fulfilled: maintenance costs of both adaptations must be accelerating (Fig. 1) and the probability of encountering defended hosts must be low (Fig. 2). Under all other conditions, pure strategies are favored. This parallels the results of Poitrineau and colleagues on host investment in defenses against multiple enemies in the case where defenses interfere with each other .
Our results show that in the case where a pure strategy is optimal, both the pure physiological resistance strategy and the pure behavioral avoidance strategy may evolve depending on initial investment levels (Fig. 1). This result is in agreement with previous population genetics models that explicitly assume "hard selection" , and could be due to our implicit assumption of "hard selection". Indeed, since there is no explicit function of population regulation in our model and because the fitness functions are not dependent on the number or strategy of local competitors, any regulation must occur on the global level.
Moreover, when the mixed investment is optimal, we found that the absolute and relative investments in each adaptation mode are sensitive to the spatial configuration of hosts (Fig. 2). When defended hosts are abundant and/or randomly distributed (i.e., when the probability of encountering defended hosts is high), it is optimal for the enemy to invest mainly in physiological resistance. We have shown that the low investment in behavioral adaptation is due to the costs of actively searching undefended hosts (Fig. 4). As the probability of encountering defended hosts decreases, the absolute investment in behavioral adaptation increases and goes through a maximum. Thereafter, optimal investments in both adaptation modes decrease. In the extreme case, total investment approaches zero because defended hosts are almost never encountered. Looking at the relative investments in each adaptation mode gives a different picture (Fig. 3). Relative investment in physiological resistance is always greater, but the difference tends to disappear with increasing rarity and/or spatial aggregation of defended hosts.
We describe the consequences of interference between two adaptive traits on the joint evolution of these traits in the first part of the Discussion below. In the second and third parts, we highlight the relevance of these findings to plant-herbivore interactions in natural and managed ecosystems.
Simultaneous investment in avoidance of and resistance to host defenses: a paradox?
Interference between physiological and behavioral adaptations is an emergent property of our model, clearly identifiable from the net fitness benefits of investment in both adaptation modes. Simplifying the enemy fitness W for null fitness costs of adaptation (see Methods, Eq. 1 a for k = e = 0) gives
W = W0+dP(1+T-I)
where T = R+B is the total investment in physiological and behavioral adaptations (respectively R and B) and I = BR is the interference between physiological and behavioral adaptations. Thus, for a given total investment T, the net benefits are maximal when all resources are invested in one mode of adaptation only (i.e., interference I is zero).
However, the optimality of a strategy also depends on the fitness costs associated with the adaptive traits. Consider an enemy investing an intermediate amount of resources in one adaptation mode and facing an increase in the frequency of encounters with defended hosts. When fitness costs are decelerating, investing in this same mode of adaptation is not only the most efficient but also the least costly way to reinforce adaptation to host defenses. Thus, continuing to invest in the prevailing adaptation mode is better than developing another from zero (i.e., pure resistance strategies are favored). In contrast, for accelerating costs, investing in the prevailing mode of adaptation is more efficient but also more costly. Thus, investing in another adaptation mode can be an optimal choice (i.e., mixed resistance strategies can be optimal). Mixed resistance strategies tend to disappear when the frequency of encounters with defended hosts increases, because investing simultaneously in two interfering adaptive traits is increasingly wasteful. In the extreme case (R•1, B•1), one half of the investment is useless because of interference.
Avoidance of aggregated novel hosts: a factor in biological invasions?
Our findings are relevant to adaptation in ecological communities, for example when a habitat is invaded by a plant variety or species that is toxic to a resident herbivore. One commonly observed life history trait in invasive plant species is clonal reproduction . This reproductive mode leads to the spatial aggregation of invaders. Thus, during the initial steps of a biological invasion, clonally invasive plants are rare and spatially aggregated. Our model suggests that natural enemies should invest a low amount of resources into adaptations to these novel hosts, and allocate non-trivial amount of these resources to behavioral adaptation. Consequently, we suggest that selection for avoidance of the toxic compounds produced by rare, clumped invasive plants could be a cause for the so-called "ecological release" experienced by these plants. The enemy release hypothesis states that plant species, on introduction to an exotic region, experience a decrease in regulation by herbivores (in particular, specialist herbivores) and other natural enemies, resulting in a rapid increase in distribution and abundance [24, 25]. Comparisons of the parasitic load and the number of pathogens in native versus introduced regions support this hypothesis [26, 27], as well as comparison of the plant anti-herbivore compounds . One approach to test our hypothesis would be to compare herbivore behavior in native and introduced ranges of invasive plants.
Pest management: how to limit physiological resistance to new crop protectants or varieties
Finally, our results are relevant to certain forms of pest management, where one attempts to conserve the efficiency of a new toxic cultivar or a new chemical pesticide. Although models of pest resistance evolution to chemical pesticides or genetically-engineered toxins have long been acknowledged as a tool for pest management , host preferences have rarely been incorporated in theoretical developments (but see [10, 29, 30]). Our model suggests that using rare and aggregated treated/toxic plants during the first years of commercialization may curtail a pest's investment in physiological resistance, and favor the evolution of avoidance of treated/toxic areas. By reducing the frequency of encounters with the new pesticide/toxin, this initial step of behavioral adaptation might delay the evolution of physiological resistance if the treated/toxic plants are subsequently used more extensively (see also ). Refuges (i.e., non-treated/toxic host plants maintained in close proximity to treated/toxic crops to delay physiological resistance evolution ) would then serve as insect traps. This potential role of refuges has rarely been studied (but see ), since in population genetics models refuges are considered to be a source of susceptible insects. Moreover, it is noteworthy that since the commercialization of insect-resistant GM crops, the optimal spatial distribution of refuges for sustainable pest control has received much attention [11–13], but their optimal temporal distribution has rarely been investigated (but see ). Based on our findings, we suggest that more research should be conducted to define this optimal temporal distribution of refuges, when taking into account the evolution of pest specialization.
The originality of our study is to have linked together physiology, behavior and landscape structure into a general model describing the adaptation of natural enemies to their hosts. Our model predicts that the optimal strategy for a natural enemy when confronted with well-defended, novel hosts subtly depends on the fitness costs of the adaptations to host defenses and the spatial distribution of defended hosts. Interestingly, under certain conditions (i.e., maintenance costs of the adaptations are accelerating and the probability of encountering defended hosts is low), a reinforcement of physiological resistance to host defenses by avoidance of these same defenses may be optimal. In this latter case, investment in physiological resistance is favored when the novel hosts are abundant at regional scales because the active search of undefended hosts is costly. It is also favored when the host type that might be encountered is difficult to predict by the enemy because the host from which the enemy emerges cannot be used as a cue (i.e., the aggregation level of host types is low).
Although they remain to be confirmed by empirical data, our theoretical results could have important implications for the management of invasive plant species and the management of pest resistance to new crop protectants or varieties. A logical next step in the model analysis would be to enable the host level evolve (see ). Indeed, since the host level is currently assumed to be non-evolving, our predictions are only relevant to enemy evolution over short time scales (e.g. evolution of herbivorous insects during the initial steps of invasion of their foraging range by new plants) or to systems in which host levels can be managed (e.g. agricultural systems). Relaxing this assumption may be useful to analyze adaptive patterns in coevolving systems having contrasting spatial structures, such as plant-herbivore interactions in tropical and temperate forests [33–36].
1. attack of an undefended host randomly encountered in the foraging area
2. attack of a defended host randomly encountered in the foraging area
3. attack of an undefended host encountered by active search in the foraging area.
Physiological and behavioral adaptations to host defenses are assumed to be quantitative traits that respectively decrease the deleterious effects of host defenses during a type (2) attack and increase the frequency of type (3) attacks.
The enemy's fitness W is stated as
W = W G *(1-C R )*(1-C B ) (1a)
where W G is a function of R and B representing the average fitness gain resulting from host attack and C R and C B are increasing functions of R and B representing the fitness costs of adaptations. The average fitness gain resulting from host attack W G is stated as
W G = (1-B)(1-P)W0+ (1-B)P(W0-H) +B(W0-S) (1b)
where the three terms respectively reflect fitness gains resulting from attacks of type (1), (2) and (3). Each term is detailed below.
The first term corresponds to the attack of an undefended host encountered by a random move in the foraging area. The fitness gain resulting from this type of attack equals W0, which is the maximal potential fitness gain resulting from a host attack. The probability of this event is (1-B)(1-P), where B is the probability for an enemy to engage in the active search of undefended hosts during the foraging step and P is the probability for an enemy to randomly encounter a defended host in its foraging area (Fig. 5). B corresponds to the investment in behavioral adaptation whereas P depends on the spatial configuration of defended and undefended hosts. This latter is described by the frequency f of defended hosts and their level a of spatial aggregation at the scale of an enemy's foraging range, which is the scale of spatial aggregation the most relevant to our study [37–39]. Host aggregation at the scale of the enemy's foraging range is described by the average frequencies of the three possible host-pair types (i.e. defended-defended, undefended-defended and undefended-undefended), when considering only host-pairs between which distance is inferior to the maximal foraging distance . The higher the frequency of homologous pairs, the greater the spatial aggregation. Making a parallel with Wright's inbreeding coefficient in population genetics , we define the aggregation level a of defended hosts at the scale of the enemy's foraging range as the deficit in heterologous host pairs. Defended-defended, undefended-defended and undefended-undefended host-pair average frequencies, f dd , f ud and f uu respectively, are defined as
f dd = f 2 +af(1-f) (2a)
f ud = 2f(1-f)(1-a) (2b)
f uu = (1-f) 2 +af(1-f). (2c)
The mean probability P that an enemy emerging from an undefended host encounters a defended host after a random move in its foraging area is
P = 1/2 f ud /(f uu + 1/2 f ud ) (3a)
which simplifies to
P = f(1-a). (3b)
The second term of equation (1b) corresponds to the attack of a defended host encountered by a random move in the foraging area. The probability of this event is (1-B)P. The fitness gain W0 is decreased by H, H being the fitness loss due to the deleterious effects of host defenses (Fig. 5). H is assumed to decrease with the level R of physiological enemy resistance and increase with the level d of host defense. We chose the simplest function to describe H . Hence
H = d(1-R). (4)
The third term of equation (1b) corresponds to the attack of an undefended host encountered by active search. The probability of this event is B. The fitness gain W0 is decreased by S, S being the fitness loss due to active searching [43, 44] (Fig. 5). S is assumed to increase with the probability P of encountering a defended host by any random move. We also chose the simplest function to describe S. Hence
S = eP (5)
where e is the search cost coefficient.
Let us now describe the constitutive fitness costs of adaptations, C R and C B (Eq. 1a). Evidence of fitness costs of adaptation to host defenses is scarce  and a fortiori, the shape of the fitness cost functions (i.e., variations in the cost magnitude with the level of investment in the adaptive trait) is largely unknown [22, 45]. However, it is reasonable to assume that physiological and behavioral adaptations to host defenses have constitutive fitness costs: a few studies show that physiological resistance to toxic compounds results from permanent metabolic changes that can reduce fitness [47–49] and obviously, discrimination between defended and undefended hosts involves energy allocation to sensors and neural cells allowing the detection and treatment of signals. Consequently, we chose simple functions to describe the constitutive fitness costs of physiological and behavioral adaptations to host defenses and we assessed the robustness of model predictions to the shape of these functions. The cost functions are taken as
C R = k R x (6a)
C B = k B x (6b)
where k is the cost coefficient and x controls the form of the function. If x > 1 then the cost is accelerating, if x = 1 then it is linear, whereas if x < 1 then it is decelerating.
Finally, note that we assume that the two adaptive traits are independent [8–10, 50]. Partial derivatives of enemy fitness W (Eq. 1a) with respect to the level of physiological adaptation R and the level of behavioral adaptation B give
∂W/∂R = (1-C B )((1-B)dP(xC R +R(1-C R -xC R ))-xC R W0)/R + ePB(1-C B )xC R /R (7a)
∂W/∂B = (1-C R )((1-R)dP(xC B +B(1-C B -xC B ))-xC B W0)/B - ePB(1-C R )(1-(1+x)C B )/B (7b)
Model parameters, their definitions, range of values employed, and notes on their use.
Defense level of defended hosts
Held at 1
Frequency of defended hosts in the environment
Spatial aggregation level of defended hosts
a = 0 when defended and undefended hosts are randomly distributed and a→1 when defended and undefended hosts form two distinct patches.
Probability of encountering a defended host during a random move in the foraging area
Increases with the frequency f of defended hosts and decreases with their aggregation level a P = f(1-a)
Maximal potential fitness gain when attacking an host
Held at 1
Corresponds to the fitness gain when attacking a randomly encountered undefended host
Fitness loss due to host defense
Decreases with the level of physiological resistance R of the enemy and increases with the level of defense d of the host H = d(1-R)
Fitness loss due to active searching of undefended hosts
Increases with the difficulty in finding undefended hosts, i.e., with the probability P of encountering a defended host during a random move S = eP
Search cost coefficient
Physiological adaptation level
Physiological adaptation reduces the fitness loss H when attacking a defended victim.
Physiological adaptation maintenance cost
Increases with the level of physiological adaptation R CR = kRx
Behavioral adaptation level
Behavioral adaptation corresponds to the probability of engaging in the active search of undefended hosts.
Behavioral adaptation maintenance cost
Increases with the level of behavioral adaptation B CB = kBx
Maintenance cost coefficient
Shape coefficient of the maintenance cost functions
Held to 1/2, 1 or 2
Maintenance costs increase with the level of adaptation in an accelerating (x = 2), linear (x = 1) or decelerating (x = 1/2) way
We thank Anthony R. Ives, Laurent Lapchin and Thomas Guillemaud for helpful comments. We acknowledge financial support from the French Ministry of Research and the CNRS (ACI "Impact of Biotechnologies on Agro-Ecosystem Functioning"). Sam P. Brown was supported by a Marie Curie Fellowship. Corinne Vacher is supported by an INRA postdoctoral fellowship.
- Chew FS: Coevolution of pierid butterflies and their cruciferous foodplants. II. The distribution of eggs on potential foodplants. Evolution. 1977, 31: 568-579.View ArticleGoogle Scholar
- Pluthero FG, Threlkeld SFH: Genetic differences in malathion avoidance and resistance in Drosophila. J Econ Entomol. 1981, 74: 736-740.View ArticlePubMedGoogle Scholar
- Pluthero FG, Sing RS, Threlkeld SFH: The behavioural and physiological components of malathion resistance in Drosophila melanogaster. Can J Genet Cytol. 1982, 24: 807-815.View ArticleGoogle Scholar
- Lockwood JA, Sparks TC, Story RN: Evolution of insect resistance to insecticides: a reevaluation of the roles of physiology and behavior. Bull Entomol Soc Am. 1984, 41-51.Google Scholar
- Hoy CW, Head GP, Hall FR: Spatial heterogeneity and insect adaptation to toxins. Ann Rev Entomol. 1998, 43: 571-94. 10.1146/annurev.ento.43.1.571.View ArticleGoogle Scholar
- Ramachandran S, Buntin GD, Tabashnik BE, Raymer PL, Adang MJ, Pulliam DA, Stewart CN: Survival, development, and oviposition of resistant diamondback moth (Lepidoptera: Plutellidae) on trangenic canola producing a Bacillus thuringiensis toxin. J Econ Entomol. 1998, 91: 1239-1244.View ArticleGoogle Scholar
- Negus TF, Ross MH: The response of German cockroaches to toxic baits: strin differences and the effects of selection pressure. Entomol Exp Appl. 1997, 82: 247-253. 10.1023/A:1002978101698.View ArticleGoogle Scholar
- Gould F: Role of behaviour in the evolution of insect adaptation to insecticides and resistant host plants. Bull Entomol Soc Am. 1984, 30: 34-41.Google Scholar
- Castillo-Chavez C, Levin SA, Gould F: Physiological and behavioral adaptation to varying environnements: a mathematical model. Evolution. 1988, 42: 986-994.View ArticleGoogle Scholar
- Rausher MD: The evolution of habitat preference: avoidance and adaptation. Evolution of insect pests: patterns of variation. Edited by: Kim KC, McPheron BA. 1994, 259-283.Google Scholar
- Wallace B: Topics in Population Genetics. 1968, New-York: W. W. NortonGoogle Scholar
- Peck SL, Ellner SP: The effect of economic thresholds and life-history parameters on the evolution of pesticide resistance in a regional setting. Am Nat. 1996, 149: 42-63.Google Scholar
- Lenormand T, Raymond M: Resistance management: the stable zone strategy. Proc R Soc Lond B Biol Sci. 1998, 265: 1985-1990. 10.1098/rspb.1998.0529.View ArticleGoogle Scholar
- Peck SL, Gould F, Ellner SP: Spread of resistance in spatially extended regions of trangenic cotton: implications for management of Heliothis virescens (Lepidoptera: Noctuidae). J Econ Entomol. 1999, 92: 1-16.View ArticleGoogle Scholar
- Caprio MA: Source-sink dynamics between transgenic and non-transgenic habitats and their role in the evolution of resistance. J Econ Entomol. 2001, 94: 698-705.View ArticlePubMedGoogle Scholar
- Vacher C, Bourguet D, Rousset F, Chevillon C, Hochberg ME: Modelling the spatial configuration of refuges for a sustainable control of pests: a case study of Bt cotton. J Evol Biol. 2003, 16: 378-387. 10.1046/j.1420-9101.2003.00553.x.View ArticlePubMedGoogle Scholar
- Cerda H, Wright DJ: Modeling the spatial and temporal location of refugia to manage resistance in Bt transgenic crops. Agric Ecosyst Environ. 2004, 102: 163-174. 10.1016/j.agee.2003.08.004.View ArticleGoogle Scholar
- Diehl SR, Bush GL: The role of habitat preference in adaptation and speciation. Speciation and its consequences. Edited by: Otte D, Endler JA. 1989, Sunderland, MA: Sinauer, 345-355.Google Scholar
- Cunningham JP, West SA: Host selection in phytophageous insects: a new explanation for learning in adults. Oikos. 2001, 95: 537-543. 10.1034/j.1600-0706.2001.950319.x.View ArticleGoogle Scholar
- Onstad DW, Spencer JL, Guse CA, Levine E, Isard SA: Modeling evolution of behavioral resistance by an insect to crop rotation. Entomol Exp Appl. 2001, 100: 195-201. 10.1023/A:1019289030248.View ArticleGoogle Scholar
- West SA, Cunningham JP: A general model for host plant selection in phytophagous insects. J Theor Biol. 2002, 214: 499-513. 10.1006/jtbi.2001.2475.View ArticlePubMedGoogle Scholar
- Poitrineau K, Brown SP, Hochberg ME: Defence against multiple enemies. J Evol Biol. 2003, 16: 1319-1327. 10.1046/j.1420-9101.2003.00585.x.View ArticlePubMedGoogle Scholar
- Sakai AK, Allendorf FW, Holt JS, Lodge DM, Molofsky J, With KA, Baughman S, Cabin RJ, Cohen JE, Ellstrand NC, McCauley DE, O'Neil P, Parker IM, Thompson JN, Weller SG: The population biology of invasive species. Annu Rev Ecol Syst. 2001, 32: 305-332. 10.1146/annurev.ecolsys.32.081501.114037.View ArticleGoogle Scholar
- Keane RM, Crawley MJ: Exotic plant invasions and the enemy release hypothesis. TREE. 2002, 17: 164-169.Google Scholar
- Joshi J, Vrieling K: The enemy release and EICA hypothesis revisited: incorporating fundamental difference between specialist and generalist herbivores. Ecol Lett. 2005, 8: 704-714. 10.1111/j.1461-0248.2005.00769.x.View ArticleGoogle Scholar
- Mitchell CE, Power AG: Release of invasive plants from fungal and viral pathogens. Nature. 2003, 421: 625-627. 10.1038/nature01317.View ArticlePubMedGoogle Scholar
- Torchin ME, Lafferty KD, Dobson AP, McKenzie VJ, Kuris AM: Introduced species and their missing parasites. Nature. 2003, 421: 628-630. 10.1038/nature01346.View ArticlePubMedGoogle Scholar
- Tabashnik B: Computer simulations as a tool for pesticide resistance management. Pesticide Resistance: Strategies and Tactics for Management. 1986, Washington, USA: National Academy Press, 194-206.Google Scholar
- Mallet J, Porter P: Preventing insect adaptation to insect-resistant crops: are seed mixtures or refugia the best strategy?. Proc R Soc Lond B Biol Sci. 1992, 250: 165-169.View ArticleGoogle Scholar
- Alstad DN, Andow DA: Managing the evolution of insect resistance to transgenic plants. Science. 1995, 268: 1894-1896.View ArticlePubMedGoogle Scholar
- Gould F: Sustainability of transgenic insecticidal cultivars. Integrating pest genetics and ecology. Annu Rev Entomol. 1998, 43: 701-26. 10.1146/annurev.ento.43.1.701.View ArticlePubMedGoogle Scholar
- Livingston MJ, Carlson GA, Fackler PL: Managing resistance evolution in two pests to two toxins with refugia. Am J Agric Econ. 2004, 86: 1-13. 10.1111/j.0092-5853.2004.00558.x.View ArticleGoogle Scholar
- Levin DA: Alkaloid-bearing plants: an ecogeographic perspective. Am Nat. 1976, 110: 261-284. 10.1086/283063.View ArticleGoogle Scholar
- Coley PD, Barone JA: Herbivory and plant defenses in tropical forests. Annu Rev Ecol Syst. 1996, 27: 305-335. 10.1146/annurev.ecolsys.27.1.305.View ArticleGoogle Scholar
- Gaston KJ: Global patterns in biodiversity. Nature. 2000, 405: 220-227. 10.1038/35012228.View ArticlePubMedGoogle Scholar
- Condit R, Ashtton PS, Baker P, Bunyavejchewin S, Gunatilleke S, Gunatilleke N, Hubbel SP, Foster RB, Itoh A, LaFrankie JV, Lee HS, Losos E, Manokaran N, Sukumar R, Yamakura T: Spatial patterns in the distribution of tropical tree species. Science. 2000, 288: 1414-1418. 10.1126/science.288.5470.1414.View ArticlePubMedGoogle Scholar
- Wiens JA: Spatial scaling in ecology. Funct Ecol. 1989, 3: 385-397.View ArticleGoogle Scholar
- Lima SL, Zollner PA: Towards a behavioral ecology of ecological landscapes. TREE. 1996, 11: 131-135.PubMedGoogle Scholar
- Kawata M, Agawa H: Perceptual scales of spatial heterogeneity of periphyton for freshwater snails. Ecol Lett. 1999, 2: 210-214. 10.1046/j.1461-0248.1999.00071.x.View ArticleGoogle Scholar
- Law R, Dieckmann U: Moment approximations of individual-based models. The geometry of ecological interactions: simplifying spatial complexity. Edited by: Dieckmann U, Law R, Metz JAJ. 2000, Cambridge: Cambridge University PressGoogle Scholar
- Hartl DL, Clark AG: Principles of Population Genetics. 1989, Sunderland, MA: Sinauer Associates, 2Google Scholar
- Kawecki TJ: Red queen meets Santa Rosalia: arms races and the evolution of host specialization in organisms with parasitic lifestyles. Am Nat. 1998, 152: 635-651. 10.1086/286195.View ArticlePubMedGoogle Scholar
- Lucas JR: Time constraints and diet choice: different predictions from different constraints. Am Nat. 1985, 126: 680-705. 10.1086/284447.View ArticleGoogle Scholar
- Crawley MJ, Krebs JR: Foraging theory. Natural enemies: the population biology of predators, parasites and diseases. Edited by: Crawley MJ. 1992, Oxford: Blackwell Scientific Publications, 90-114.View ArticleGoogle Scholar
- Coustau C, Chevillon C, Ffrench-Constant R: Resistance to xenobiotics and parasites: can we count the cost?. TREE. 2000, 15: 378-383.PubMedGoogle Scholar
- Boots M, Haraguchi Y: The evolution of costly resistance in host-parasite systems. Am Nat. 1999, 153: 359-370. 10.1086/303181.View ArticleGoogle Scholar
- McKenzie JA, Purvis A: Chromosomal localization of fitness modifiers of diazon resistance genotypes of Lucilia cuprina. Heredity. 1984, 53: 625-634.View ArticleGoogle Scholar
- Groeters FR, Tabashnik BE, Finson N, Johnson MW: Fitness costs of resistance to Bacillus thuringiensis in the diamondback moth (Plutella xylostella). Evolution. 1994, 48: 197-201.View ArticleGoogle Scholar
- Raymond M, Chevillon C, Guillemaud T, Lenormand T, Pasteur N: An overview of the evolution of overproduced esterases in the mosquito Culex pipiens. Philos Trans R Soc Lond B Biol Sci. 1998, 353: 1707-1711. 10.1098/rstb.1998.0322.PubMed CentralView ArticlePubMedGoogle Scholar
- Hochberg ME: Hide or Fight? The competitive evolution of concealment and encapsulation in parasitoid-host associations. Oikos. 1997, 80: 342-352.View ArticleGoogle Scholar
- Murray JD: Continuous models for interacting populations. Mathematical Biology, second, corrected edition. Edited by: Murray JD. 1993, Berlin/Heidelberg: Springer-VerlagGoogle Scholar
- Wolfram S: Mathematica Version 184.108.40.206. 1999, Champaign, Ill: Wolfram ResearchGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.