Evolved pesticide tolerance influences susceptibility to parasites in amphibians

Abstract Because ecosystems throughout the globe are contaminated with pesticides, there is a need to understand how natural populations cope with pesticides and the implications for ecological interactions. From an evolutionary perspective, there is evidence that pesticide tolerance can be achieved via two mechanisms: selection for constitutive tolerance over multiple generations or by inducing tolerance within a single generation via phenotypic plasticity. While both mechanisms can allow organisms to persist in contaminated environments, they might result in different performance trade‐offs including population susceptibility to parasites. We have identified 15 wood frog populations that exist along a gradient from close to agriculture and high, constitutive pesticide tolerance to far from agriculture and inducible pesticide tolerance. Using these populations, we investigated the relationship between evolutionary responses to the common insecticide carbaryl and host susceptibility to the trematode Echinoparyphium lineage 3 and ranavirus using laboratory exposure assays. For Echinoparyphium, we discovered that wood frog populations living closer to agriculture with high, constitutive tolerance experienced lower loads than populations living far from agriculture with inducible pesticide tolerance. For ranavirus, we found no relationship between the mechanism of evolved pesticide tolerance and survival, but populations living closer to agriculture with high, constitutive tolerance experienced higher viral loads than populations far from agriculture with inducible tolerance. Land use and mechanisms of evolved pesticide tolerance were associated with susceptibility to parasites, but the direction of the relationship is dependent on the type of parasite, underscoring the complexity between land use and disease outcomes. Collectively, our results demonstrate that evolved pesticide tolerance can indirectly influence host–parasite interactions and underscores the importance of including evolutionary processes in ecotoxicological studies.


| INTRODUCTION
For over 70 years, pesticide use has improved human health and enhanced agricultural yields by reducing populations of disease vectors and pest species (Pimentel, 2005). One major consequence of pesticide use is the evolution of pesticide tolerance in pest species (Georghiou, 1990;Hoy, 1998;Weill et al., 2003). For nearly every major class of agrochemical, the evolution of tolerance has been reported in at least one pest species (Pimentel, 2005). Additionally, the evolution of pesticide tolerance in pest species is estimated to cause >$1.5 billion in crop losses each year (Pimentel, 2005). Despite our increasing knowledge of evolved pesticide tolerance in pest species, few studies have explored whether pesticide tolerance has evolved in nontarget species that are inadvertently exposed to these chemicals (Brausch & Smith, 2009a,b;Jansen, Coors, Stoks, & De Meester, 2011). Pesticide application for pest control is typically localized; however, processes such as aerial drift, runoff, and movement through food webs can distribute pesticides across the landscape potentially shaping the evolutionary courses of natural populations (Gilliom, 2007;Stone, Gilliom, & Ryberg, 2014). The potential for nontarget organisms to evolve tolerance is critical to consider in toxicological assays because it affects estimated toxicity metrics and therefore risk assessments conducted for nontarget species.
There is increasing evidence that nontarget populations can evolve increased tolerance to pesticides (Bendis & Relyea, 2014;Brausch & Smith, 2009a;Cothran, Brown, & Relyea, 2013;Hua et al., 2015;Jansen et al., 2011). Theory suggests that the evolution of pesticide tolerance can be achieved via two mechanisms: selection for constitutive tolerance over multiple generations or by inducing tolerance within a single generation via phenotypic plasticity Hua et al., 2015;Jansen et al., 2011). The evolution of constitutive tolerance is more likely when populations are in close proximity to agriculture (<200 m) and consistently exposed to pesticides (Crispo, 2007;Declerck et al., 2006).
In contrast, inducible tolerance is more likely when populations are far from agriculture, infrequently exposed to pesticides, and there are costs of maintaining pesticide tolerance when pesticides are not present (i.e., phenotypic trade-off; Crispo, 2007;Crispo et al., 2010). While constitutive and inducible pesticide tolerance can improve survival following pesticide exposure, both mechanisms can alter how organisms interact with other community members (i.e., competition, predator-prey, host-parasite interactions; Gotthard & Nylin, 1995;Hanazato, 1995;McCarroll & Hemingway, 2002;Raberg, 2014). Thus, to advance the field of toxicology, there is a need to not only consider how evolutionary responses to pesticides affect estimates of pesticide toxicity, but also how these evolutionary responses alter ecological interactions.
The influence of pesticides on host-parasite interactions has been explored in many systems (Blaustein et al., 2011;Dietrich, Van Gaest, Strickland, & Arkoosh, 2014;Doublet, Labarussias, de Miranda, Moritz, & Paxton, 2015). Several studies suggest that the interaction between pesticides and parasites can contribute to declines in wild populations either directly (e.g., mortality) or indirectly when the presence of the pesticide alters the traits of the parasite or the host (e.g., immunosuppression; Christin et al., 2003;Kiesecker, 2002;Rohr, Raffel, Sessions, & Hudson, 2008). Fundamentally different from these mechanisms is the possibility that populations with evolved constitutive or inducible tolerance to pesticides experience altered host-parasite interactions, even in the absence of the pesticide, due to phenotypic trade-offs.
In contrast, A. gambiae that achieve tolerance via target-site insensitivity (i.e., knock-down resistance; Protopopoff et al., 2013) are more susceptible to Plasmodium (Alout et al., 2016;Ndiath et al., 2014). These studies demonstrate that the different mechanisms of pesticide tolerance can differentially shape disease outcomes in target species. However, no studies have explored these questions in nontarget organisms.
Wood frog populations inhabit pond ecosystems that can encounter pesticides following direct application, runoff, or aerial drift (Smalling et al., 2015). Using 15 populations in northwestern Pennsylvania (Hua et al., 2013(Hua et al., , 2014, we found that populations living close to agriculture have higher baseline tolerance to a common acetylcholineesterase (AChE) inhibiting pesticide (carbaryl) compared to populations living far from agriculture . Using the same populations, we also investigated patterns of plasticity to carbaryl and found that the degree of plasticity to carbaryl varied along an agricultural gradient; populations close to agriculture possess high baseline but no inducible tolerance to carbaryl (i.e., constitutive tolerance), whereas populations far from agriculture possess low baseline but high inducible tolerance to carbaryl ( Figure 1). Additionally, we have demonstrated that these landscape patterns were generalizable; wood frog populations with constitutive and inducible tolerance to carbaryl also had constitutive and inducible cross-tolerance to other AChE inhibitors (Hua et al., 2013(Hua et al., , 2014. Collectively, the spatial patterns of constitutive and inducible tolerance provide a useful framework for exploring the relationship between the mechanisms of pesticide tolerance and parasite susceptibility. Using experimental parasite exposures, here we hypothesized that in an environment without pesticides, populations with constitutive tolerance (i.e., populations close to agriculture with high baseline tolerance to pesticides) should be more susceptible to both parasite species compared to populations with inducible tolerance (i.e., living far from agriculture with low baseline tolerance).

| Animal collection and husbandry
We collected wood frog eggs from 15 populations (Table 1) in northwestern Pennsylvania. The populations were separated by at least 4 km as past studies demonstrate that the genetic neighborhood is generally within ~1 km of the breeding pond (Berven & Grudzien, 1990;Semlitsch, 2000;Semlitsch & Bodie, 1998). On April 8, 2014, we collected 10 egg clutches (Gosner stage 3-7;Gosner, 1960) from each population. On 21 April, the eggs (Gosner stage 18) were shipped F I G U R E 1 The relationship between distance to agriculture, magnitude of baseline tolerance, and magnitude of inducible tolerance to carbaryl. Regression lines are reproduced using data from Hua et al. (2015). Distance to agriculture (m) was calculated as the shortest distance from collection location to an agricultural field. The measure of baseline tolerance was calculated as the percent survival following exposure to a lethal concentration of carbaryl. The measure of inducible tolerance is represented by the hazard regression coefficient (more negative values indicate populations with higher plasticity to pesticides). Exact values are reported In d u c ib le to le ra n c e r = .5 6 p = .0 1 B a s e lin e to le ra n c e r = -.4 7 p = .0 4 Reproduced from Hua et al. 2015 T A B L E 1 The abbreviations, coordinates, distance to agriculture (m), measure of constitutive tolerance (percent survival following exposure to carbaryl), and measure of inducible tolerance (hazard regression coefficient; the more negative values indicate populations with higher pesticide-induced tolerance) for the wood frog populations used in this study. All values are derived from Hua et al. (2015) Pond ID Latitude ( We placed all masses into in 100-L pools filled with 90 L of well water, keeping eggs from the same populations together. By 28 April, all tadpoles reached Gosner stage 25 and were fed rabbit chow ad libitum until the start of the experiments.

| Trematode experiment
To obtain Echinoparyphium for the experiments, we collected adult ramshorn snails (Helisoma trivolvis), the first intermediate host of the trematode, from the PWA. We screened the snails for Echinoparyphium infection by isolating individuals in 50-ml tubes filled with 45 ml of aged, UV-irradiated well water and placed the tube under a light source for one hour to induce the shedding of the free-swimming cercariae (Hua, Buss, Kim, Orlofske, & Hoverman, 2016). We selected five snails that shed the highest density of cercariae and housed them separately in 2-L plastic containers with 1.5 L of aged, UV-irradiated well water at 7°C to slow shedding of cercariae until the start of the experiment. Three days prior to the start of the experiment, snails were acclimated to 25°C and fed rabbit chow ad libitum. We identified cercariae as Echinoparyphium lineage 3 using standard molecular sequencing of the ITS1 gene and implementation of Bayesian phylogenetics methods (Detwiler, Bos, & Minchella, 2010;Hua et al., 2016).
To investigate the susceptibility of tadpoles to trematodes, we selected 20 tadpoles of a similar size from each population and allowed them to acclimate to indoor conditions for 24 hr on 18 May. During this acclimation period, we haphazardly distributed tadpoles into 1-L plastic containers filled with 500 ml of aged, UV-irradiated well water (pH 7.4; 19.1°C; 14:10 light:dark cycle). Each container contained a single tadpole, and we fed all tadpoles TetraMin (Tetra Spectrum Brands) ad libitum. On 19 May, we individually exposed 15 tadpoles from each of the 15 populations to 50 free-swimming trematodes (cercariae) for a total of 225 experimental units. To account for any potential background mortality, we also included a no-trematode control for each population, which was replicated five times for an additional 75 experimental units.
To obtain the cercariae, we placed the previously collected infected snails into 50-ml tubes filled with 45 ml of aged, UV-irradiated well water under a light source to induce shedding. We then used a stereo dissecting scope and glass pipette to isolate 50 cercariae into individual 50-ml Falcon tubes containing 45 ml of aged, UV-irradiated aged well water. Within 4 hr of counting, we added the cercariae to each experimental unit. To ensure that all cercariae were added to each treatment, each tube was rinsed three times using water from the experimental unit. For the control treatment, we repeated this procedure by adding the same volume of water from an uninfected snail.
Tadpoles were checked daily for mortality and fed TetraMin ad libitum. We terminated the experiment 3 days postexposure because past studies demonstrate this is enough time for parasites to successfully encyst and before parasite clearance begins (Hoverman, Hoye, & Johnson, 2013). Tadpoles were euthanized using a MS-222 overdose and then preserved in 10% formalin. We weighed, staged, and measured snout-vent length for each individual (Table S1). The tadpoles were then dissected to quantify parasite load. We then calculated the proportion of trematodes encysted (i.e., number of encysted trematodes out of the 50 cercariae added to each experimental unit). We scanned the entire body of the tadpole for encysted Echinoparyphium but found all metacercariae in the kidneys. We also examined all control individuals to confirm the absence of trematode infection; all individuals were negative.

| Ranavirus experiment
We used a FV3-like isolate in our experiment that was obtained Prior to quantifying viral loads, we weighed, staged, and measured snout-vent length for each individual (Table S2). Then, we dissected each individual to remove the kidneys and liver to quantify viral load.

| Statistical analysis
Distance to agriculture, baseline pesticide tolerance, and the magnitude of inducible tolerance are highly correlated factors that contribute to evolved pesticide tolerance in natural wood frog populations . Therefore, we conducted a factor reduction analysis using Principal Axis Factoring with an orthogonal Varimax rotation (SPSS; factor reduction) to facilitate interpretations of the relationship between the mechanism of pesticide tolerance and susceptibility to parasites. Using the Guttman-Kaiser criterion (i.e., the number of factors are determined by those with eigenvalues over 1.0; Yeomans & Golder, 1982), we reduced our three variables into a single predictor (PC-1). Values for distance to agriculture, baseline pesticide tolerance, and inducible tolerance were obtained from a previous study that tested wood frogs from the same egg masses used in this study ; Table 1). Due to differences in the total number of populations tested in the trematode and ranavirus experiment (15 vs. 14, respectively), we conducted two separate factor reduction analysis and calculated two separate PC-1 for the trematode and ranavirus experiments.
For the trematode experiment, we conducted a regression analysis (SPSS Regression) to assess the relationship between PC-1 for the 15 populations and trematode resistance (average proportion of trematodes encysted). We did not conduct a regression analysis assessing the relationship between PC-1 and tadpole survival because all tadpoles survived the trematode and no-trematode exposures.
For the ranavirus experiment, we conducted four regression analyses (SPSS regression). The first analysis assessed the relationship between PC-1 for the 14 populations and average percent tadpole survival following a 120-hr exposure to ranavirus. The second analysis assessed the relationship between PC-1 and average tadpole time to death (TTD). The third analysis assessed the relationship between percent survival and viral load, using a binomial logistic regression (SPSS Logistic regression). Because an individual's survival is directly related to its viral load, we conducted our logistic regression on individual tadpoles rather than population means. To account for the potential effect of experimental units, we included experimental unit as a covariate.
Finally, we conducted a fourth regression analysis to determine the relationship between PC-1 and ranavirus load. For this analysis, we conducted two separate linear regressions to account for differences in the viral load of individuals that survived versus those that did not survive. For both linear regressions, we conducted our analyses on the average load of tadpoles from each population. We accidentally did not add ranavirus to one of the replicates from the SKN population, so we removed this replicate from the analysis. All variables in the linear regression analyses met assumptions of normality (Shapiro-Wilk; p > .05).  Tables S3 and S4).

| Trematode experiment
We had 100% survival in the no-trematode and trematode exposure treatments. The proportion of trematodes encysted varied across the 15 populations with population averages ranging from 23.8% to 52.8% (Figure 2; trematode loads reported in Fig. S1). We found a significant negative relationship between PC-1 and the proportion of trematodes encysted in tadpoles from each population (r = −.55; p = .03; Figure 2). Populations that live closer to agriculture and have higher baseline tolerance had a lower proportion of trematodes encysted compared to tadpoles that live further from agriculture and have higher inducible tolerance.
Therefore, we did not examine the relationship between survival of tadpoles not exposed to ranavirus and PC-1. For tadpoles exposed to ranavirus, average survival, time to death, and log viral load of each population ranged from 0% to 4.7%, 190-224 hr, and 4.6-5.9 viral copies/ng of DNA, respectively (Fig. S2). We found no relationship between PC-1 and mean tadpole survival among populations exposed to ranavirus (r = .11; p = .70; Fig. S3) nor between PC-1 and mean tadpole TTD among the populations exposed to ranavirus (r = .20; p = .42; Figure 3). However, survival was related to ranavirus load (χ 2 = 249.8; p < .001; df = 2; Fig. S4); tadpoles with higher ranavirus loads were less likely to survive than tadpoles with lower loads (p < .001). The logistic regression model explained 61.5% (Nagelkerke R) of the variance in survival and correctly classified 92.8% of cases. There was no significant effect of experimental unit (p = .328).
Finally, for tadpoles that survived ranavirus exposure, we found a significant positive relationship between PC-1 and ranavirus load (r = .60; p = .03; Figure 3). Populations that live closer to agriculture and have higher baseline tolerance had higher viral loads compared to tadpoles that live further from agriculture and have higher inducible tolerance. This relationship was not influenced by differences in the number of individuals surviving from inducible versus constitutive populations (Fig. S3). For tadpoles that died following ranavirus exposure, we found a marginally significant, positive relationship between PC-1 and ranavirus load (r = .48; p = .08; Figure 3).

| DISCUSSION
Using 15 wood frog populations, we demonstrated that differences in evolved pesticide tolerance associated with proximity to agriculture influenced susceptibility to parasites, but the patterns depended on the type of parasite. When exposed to trematodes, populations closer to agriculture with less inducible tolerance and higher baseline tolerance were less susceptible to trematodes (as indicated by a lower proportion of trematodes encysted) compared to populations that are further from agriculture with more inducible tolerance and lower baseline tolerance. In contrast, for individuals that survived the ranavirus experiment, those from populations closer to agriculture with lower inducible tolerance and higher baseline tolerance had higher viral loads compared to those populations further F I G U R E 2 The relationship between PC-1 and the average proportion of trematodes encysted in wood frog tadpoles from 15 populations (r represents the correlation coefficient). To create PC-1, we reduced our three variables (distance to agriculture, magnitude of baseline tolerance, and magnitude of inducible tolerance) into a single predictor (PC-1). Each point represents the average proportion of trematodes that successfully encysted in tadpoles from a population Populations closer to agriculture with less inducible tolerance and higher baseline tolerance were more resistant to Echinoparyphium (lower proportion of trematodes encysted) compared to populations further from agriculture with more inducible tolerance and lower baseline tolerance. While the mechanisms underlying this relationship are unknown, research on mosquitoes has found that constitutive pesticide tolerance was associated with an alteration of the redox potential within mosquito tissue, which ultimately created a toxic environment for embedded roundworm larvae (McCarroll & Hemingway, 2002).
In a study using the toad Rhinella arenarum, researchers found that exposure to malathion, an AChE inhibitor, also caused a modification in cell redox potential (Anguiano, Caballero de Castro, & Pechen de D'Angelo, 2001). Thus, one hypothesis for our results is that selection for constitutive tolerance-but not inducible tolerance-indirectly influences trematode resistance by initiating an internal environment that is toxic to trematodes. Given that host pathology (e.g., hemorrhaging, edema, mortality) is dose-dependent for Echinoparyphium, a reduction in trematode load could increase host fitness (Huffman & Fried, 2012;Orlofske, Belden, & Hopkins, 2009). However, our study focused on short-term resistance to trematode infection. While individuals from populations further from agriculture with higher inducible tolerance and lower baseline tolerance were less resistant to trematodes, they could display higher levels of parasite clearance or tolerance compared to individuals from populations closer to agriculture with lower inducible tolerance and higher baseline tolerance (Raberg, 2014). Integrating short-term (i.e., resistance) and long-term parameters (i.e., clearance and tolerance) will provide a more complete understanding of how evolved pesticide tolerance influences hosttrematode interactions.
We found that 100% of the wood frog tadpoles that were exposed to ranavirus became infected. Average final survival varied minimally across the populations (0% to 4.7%) suggesting that all populations, regardless of pesticide tolerance mechanism, were highly susceptible to ranavirus. This is consistent with past work demonstrating that wood frogs are one of the most susceptible amphibian species to ranavirus infection (Hoverman et al., 2011). Given the limited variability in survival across the populations, it is not surprising that we did not find a relationship with the mechanism of pesticide tolerance. However, we also measured time to death (TTD), which provides an alternative metric for measuring susceptibility to ranavirus. Although average TTD across the 14 populations varied more than survival (TTD ranging from 190 to 224 hr), we still did not find a relationship between the mechanism of pesticide tolerance and susceptibility to ranavirus. virus (e.g., hylids;Hoverman et al., 2011) may be better suited to address our hypotheses regarding the interaction with evolved pesticide tolerance because there is the potential for more variation in survival rates. However, wood frogs are the only amphibian species that has been examined for patterns in evolved pesticide tolerance.

Amphibian species with intermediate levels of susceptibility to rana
While there was no relationship between the mechanism of pesticide tolerance and survival with ranavirus, we found a significant relationship between the mechanism of pesticide tolerance and ranavirus load. Ranaviral load is associated with increased mortality; thus, higher loads may result in decreased survival in the wild (Wuerthner, Hua, & Hoverman, 2017). For tadpoles that survived the ranavirus exposure, individuals from populations closer to agriculture with lower inducible tolerance and higher baseline tolerance had markedly higher viral loads than those further from agriculture with higher inducible tolerance and lower baseline tolerance. Indeed, we found a 230% increase in viral loads from the least to most susceptible population. In nature, wood frog tadpoles interact with a diversity of other amphibian species (Werner et al., 2014) and wood frogs are often implicated as important species in ranavirus transmission (Harp & Petranka, 2006).
As such, the wide variation in viral loads between constitutive versus inducible populations may have broad implications for understanding transmission dynamics of ranavirus within communities.
We also found a significant positive relationship between tadpole mortality and viral load. Previous studies have reported a similar correlation between viral load and mortality (Wuerthner et al., 2017). Our results suggest that regardless of population, wood frogs experience mortality when log viral loads passed a threshold of approximately 5 viral copies/ng of DNA. A similar threshold was detected when log viral loads approached 4.5 viral copies/ng of DNA in gray tree frogs (Hyla versicolor; Wuerthner et al., 2017). Identifying a threshold value can potentially allow researchers to make predictions about the likelihood or timing of a ranavirus die-off in natural populations.
Our results are consistent with the growing body of work demonstrating population-level variation in the susceptibility of hosts to parasites (Bradley et al., 2015;Breitburg et al., 2015;Koprivnikar, Baker, & Forbes, 2006;Matson, 2006;Pearman & Garner, 2005;Schock, Bollinger, Chinchar, Jancovich, & Collins, 2008;Schock, Bollinger, & Collins, 2010). Further, this study provides the first evidence in a vertebrate species that evolutionary responses to pesticides influence susceptibility to parasites. Pesticide-mediated effects on disease outcomes are commonly documented in the literature, yet the results are often variable with pesticides increasing host tolerance to parasite infection in some cases (Coors, Decaestecker, Jansen, & De Meester, 2008;Kiesecker, 2002; but not in others (Griggs & Belden, 2008;Marcogliese et al., 2009;Schotthoefer et al., 2011). We demonstrate that the mechanism underlying the evolution of pesticide tolerance may be critical for assessing susceptibility to parasites. Moreover, we observed that the direction of the relationship varies depending on the parasite species. Thus, considering the evolutionary history of populations may shed light on the equivocal nature of pesticide-parasite interactions documented in the literature.

| Future considerations
While we found that evolved pesticide tolerance can influence pesticide-parasite interactions, there are several important future research directions. For instance, our study focused on two parasite species and a single pesticide, host, and life stage. Because variation in these factors is known to influence disease outcomes (Bridges, 2000; Hammond, Jones, Stephens, & Relyea, 2012; Jones, Hammond, & Relyea, 2009;Kegley, Hill, Orme, & Choi, 2010), future work expanding upon this research is critical for identifying general trends and developing predictions about the role of evolved pesticide tolerance on host-parasite interactions.
We demonstrated a relationship between evolved pesticide tolerance and susceptibility to parasites, but the mechanisms driving these patterns are beyond the scope of this study. The consequences of evolutionary responses to pesticides on parasite susceptibility are more commonly addressed in studies of vector species. For example, in the mosquito, A. gambiae, the evolution of pesticide tolerance can arise via multiple mechanisms (e.g., production of a thicker cuticle or physiological defenses such as target-site insensitivity; Balabanidou et al., 2016) that can lead to direct and indirect selection for traits that also confer resistance or tolerance to parasites. Therefore, interdisciplinary approaches that integrate discoveries in vector research may greatly facilitate our understanding of how evolutionary responses to pesticides influence susceptibility to parasites in nontarget species.
In this study, we focused only on assessing how variation in host tolerance to pesticides and proximity to agriculture affected disease outcomes. However, variation in the host, parasite, and environment (i.e., points in the disease triangle; Stevens, 1960;Budria & Candolin, 2014) all can influence disease outcomes. Similar to population-level variation in the host populations, parasites can also vary in their tolerance to pesticides (Hua et al., 2016). Additionally, we conducted our study in pesticide-free experimental environments, but if pesticides were present in the environment, the outcomes could be altered.
Pesticides can perturb host-parasite interactions by directly reducing host or parasite abundance (i.e., density-mediated effects). Pesticides in the environment may also alter traits, such as host immune responses, parasite infectivity, host behavior, or parasite behavior (i.e., trait-mediated effects Rohr et al., 2003Rohr et al., , 2004Weis, Smith, Zhou, Santiago-Bass, & Weis, 2001). Thus, future studies that consider evolutionary responses of parasites to pesticides and incorporate various environmental conditions would be valuable in determining the generalizability of the relationship we detected between evolved pesticide tolerance and parasite tolerance.
Our study focused on how evolved tolerance to pesticides affects host susceptibility to parasites. However, the reverse could be true as evolutionary responses to parasites could influence pesticide tolerance. For example, in mosquitoes, studies suggest that hosts that expended resources toward resisting infection by Plasmodium (i.e., detoxification enzymes) had increased vulnerability to insecticides (Saddler, Burda, & Koella, 2015). Past amphibian studies have also suggested scenarios that may lead to evolutionary responses to parasites. For example, high-nutrient agricultural run-off may indirectly facilitate trematode abundance by increasing the amount of algal resources available to the first intermediate hosts (snails) of trematodes (Koprivnikar & Redfern, 2012;Koprivnikar et al., 2006;Rohr, Schotthoefer, et al., 2008). Higher abundances of Echinoparyphium may select for increased resistance or tolerance to Echinoparyphium in amphibian populations near agriculture. Therefore, future studies that consider the possibility of evolved responses to parasites and their influence on disease outcomes would complement our understanding of how evolutionary processes affect responses to chemicals.
This study suggests that considering patterns of constitutive and inducible tolerance is critical in studies of toxicology as the predicted toxicity of a pesticide can depend on the evolutionary history of test populations. Evolutionary responses to pesticides not only effect estimates of pesticide toxicity but can also have indirect consequences by altering ecological interactions between hosts and their parasites.
Therefore, considering evolutionary responses in toxicology is critical to understanding how contaminants affect ecosystems.