Abstract
Xenobiotic resistance is commonly found in species with short generation times such as bacteria, annual plants, and insects. Nevertheless, the fundamental evolutionary principles that govern the spread of resistance alleles hold true for species with longer generation times. One such example could occur with sea lamprey (Petromyzon marinus), a parasitic invasive species in the Laurentian Great Lakes that decimated native fish populations prior to its control with the pesticide 3-trifluoromethyl-4-nitrophenol (TFM). Since the 1950s, tributaries have been treated annually with TFM, where treatments effectively remove most, but not all, larval sea lamprey. We developed an eco-genetic model of sea lamprey to examine factors affecting the evolution of resistance and found that resistance alleles rapidly rise to fixation after 40–80 years of treatment, despite the species’ relatively long generation time (4–7 years). The absence of natal homing allows resistant individuals to spread quickly throughout the entire system, but also makes the early detection of resistance challenging. High costs of resistance and density independent reproduction can delay, but not prevent, the onset of resistance. These results illustrate that sea lamprey have the potential to evolve resistance to their primary control agent in the near future, highlighting the urgent need for alternative controls.
Similar content being viewed by others
Introduction
The rapid evolution of resistance to xenobiotic compounds has been widely documented in microbes1,2, fungi3, invertebrates4,5, and plants6,7. Resistance has also been documented in vertebrates8,9,10, but has largely been constrained to taxa with high fecundity and short generation times. Nevertheless, the fundamental principles governing the evolution of resistance still apply to species with longer generation times. If the selection pressure imposed by a xenobiotic is strong, wide-spread, and consistently applied year after year, then resistance may still evolve.
The need to control nuisance pest species through various means including xenobiotics continues to rise with globalized trade and climate change introducing species into new habitats. Aquatic species are no exception, with some of the more recent threats including invasive fishes11,12. Invasive species cost billions of dollars to control and they can impact ecosystems in myriad ways13. Thus, there is a need to consider the pace at which a broader range of taxa, including those with longer generation times such as fishes, could evolve resistance to control measures.
In this study, we examine the timelines of pesticide resistance evolution in an aquatic invasive fish species, the sea lamprey (Petromyzon marinus). In the Laurentian Great Lakes, a powerful chemical pesticide has been applied to control invasive sea lamprey for over 60 years14. The pesticide, 3-trifluoromethyl-4-nitrophenol (hereafter TFM), has been applied on an annual basis since the 1950s to Great Lakes tributaries, targeting sea lamprey during their larval stage. Because the pesticide kills most, but not all, larval sea lamprey15,16, some individuals are likely to survive exposure. If this process is repeated over a long enough period, then resistant individuals could increase in frequency – a scenario documented in many systems where pests have been controlled by chemical means17. The evolution of resistance would greatly reduce the effectiveness of TFM, which remains the primary control agent, and could result in large declines to native fish populations due to increases in sea lamprey abundance. These declines in fish abundance would not only have large ecological repercussions, but could also result in billions of dollars in economic damage18,19.
As juveniles, sea lamprey feed on the blood and tissues of diverse host species, including species of commercial importance such as salmon and lake trout. Invasive sea lamprey impact native fish populations by wounding and often killing the host fishes that they parasitize18,20. Following their invasion in the late 1930s, sea lamprey contributed to the catastrophic loss of economically valuable commercial and recreational fisheries in Canada and the United States21. In response to the proliferation of sea lamprey throughout the Great Lakes, there was an immediate and concerted effort to develop efficient means of control. One effort, initiated in the 1950s, involved testing over 6,600 chemical compounds on sea lamprey and other fish species22. The organic compound TFM was found to effectively kill larval sea lamprey and had few detectable effects on other fish species at low concentrations. TFM control was deemed highly successful and sea lamprey abundance was reported to have been reduced by approximately 90%, resulting in recovery of economically important fisheries23. This outcome is one of the few documented cases of an invasive vertebrate species being successfully controlled at an ecosystem scale by a pesticide. Alternative sea lamprey control measures have been implemented on smaller scales, but none have proven as effective as TFM, with the possible exception of existing barriers that prevent upstream migration to spawning habitat24,25,26. However, barriers are not present on most tributaries and there is increased pressure from stakeholders to remove existing barriers in order to restore natural connectivity in waterways27. Thus, there is a continued reliance on TFM to manage invasive sea lamprey populations in the Great Lakes. This reliance on a single control measure can be risky; only a handful of manufacturers produce TFM and thus both prices and supplies of TFM could change. More importantly, reliance on a single chemical control measure has been shown to increase the likelihood of resistance evolution28,29.
There are several reasons why invasive sea lamprey could be evolving resistance to TFM. First, sea lamprey, like many pest species, have high fecundity. A single female lamprey can produce over 100,000 eggs in the Great Lakes30 and substantially more in its native range31. With such high fecundity, it may not take long for novel mutations that confer resistance to appear in sea lamprey offspring. Alternatively, there may be resistant individuals already present within the population (i.e., selection may act on standing genetic variation). Second, the selection pressure imposed by TFM may be substantial. Similar to many plant and insect pests, sea lamprey larvae can be found at high density; sea lamprey larvae congregate in lentic, sandy tributaries and this life-history strategy means that large numbers of individuals can be targeted with a single application of TFM. Estimates of larval mortality per stream treated with TFM range from 80–99%32, representing a substantial selective event. Third, the selection pressure from TFM has been present for nearly 60 years. TFM has been applied to Great Lakes tributaries almost every year since 195814. Thus, the combination of high fecundity, strong selection, and long-term treatment make the evolution of resistance a distinct possibility.
Despite the possibility for the evolution of resistance, there are several reasons why resistance evolution may be prevented or delayed. One reason is that sea lamprey have overlapping generations and they can live for a relatively long time (e.g., up to 19 years in its native range33). This longevity means that sea lamprey generation times are longer and their response to selection slower in comparison with many other pest species. Assuming an average generation time of between 3 and 5 years, invasive sea lamprey may have been exposed to TFM for 12–20 generations, which may not be sufficient time for resistance to evolve (but see34,35,36). Another reason that resistance may not evolve is the high rate of dispersal and subsequent gene flow found in this species. Adult sea lamprey do not return to their natal streams to spawn, instead choosing spawning sites via the detection of pheromones released by larvae37. Streams with more larvae release greater pheromone concentrations and attract more adults to spawn. This mixing of adults during reproduction results in high gene flow and near-panmixia throughout the Great Lakes38. The effects of this lack of natal philopatry on the evolution of resistance remain largely unknown, but the high gene flow between treated and non-treated streams may allow for reservoir populations to limit the spread of resistant genotypes39. Previous empirical work has not found any evidence of TFM resistance in invasive sea lamprey34. However, there were limitations to the empirical analysis and, as we show in this manuscript, detection of resistance in anything but its latest stages requires vast sample sizes.
In this study, we used the sea lamprey as a model system to examine the likelihood and timelines of resistance evolution in a species with a relatively long generation time and that has experienced a potentially potent selective pressure from the persistent application of a chemical pesticide. To predict the evolution of resistance in this system, we built an eco-genetic model40 that includes important aspects of sea lamprey life history, spatial structure, and TFM control. Specifically, we answer five questions: 1. How quickly can resistance evolve under varied strengths of selection? 2. What level of treatment intensity results in the fewest lamprey over the long term? 3. How does the lack of natal philopatry affect the development and spread of resistance? 4. How do costs of resistance affect the persistence and timing of resistance? and 5. How does the relative magnitude and timing of density dependent and density independent process affect resistance? We conclude with a discussion of caveats and model assumptions, but nevertheless suggest that if agencies are as effective at controlling invasive lamprey as they believe, then the evolution of resistance is imminent.
Materials and Methods
To examine the evolution of TFM resistance in sea lamprey, we constructed a forward-time, individual-based model (IBM). The model is a type of eco-genetic model40, which follows the evolution of key traits while also including important ecological feedbacks such as density-dependent reproduction and survival. The model is fully age and stage structured (including larvae, parasitic juveniles, and reproductively mature adults) and is based on life history characteristics of sea lamprey from the Great Lakes (see Supporting Information for detailed life history information; see Table 1 for default parameters). Each simulation was initiated by creating larval and juvenile individuals of varying ages and sexes (50:50 sex ratio). Larvae were split among tributaries and juveniles were assigned to a single parasitic population20,38. All individuals within the model were assigned a growth parameter, k, as a random deviate from a normal distribution with a mean of 0.001 and a variance of 0.0002. This growth parameter remained with the individual for their lifetime and was used in the Von Bertalanffy growth equation to determine size at age41:
where \(L(t)\) equals the body length of an individual at time t, \({L}_{\infty }\) equals the maximum length, here set to 150 mm for larvae and 800 mm for juveniles42, k equals the individual growth parameter described above, t equals the individual’s age in days, and where \({t}_{0}\) was set to −15. Examples of individual growth curves for larvae and juveniles using these parameters are illustrated in Fig. 1. Size-at-age was used to determine when individuals would transition between stages. When a larva was greater than or equal to 120 mm in length, the individual transformed into a parasitic juvenile15,32. Likewise, when a parasitic juvenile was greater than or equal to 450 mm in length, the individual transitioned to reproductive maturity42,43.
The first step in the model was ageing, where every individual was aged by one year (Fig. 1). This step was followed by transformation (i.e., larvae to juveniles) or maturation (i.e., juveniles to adults) as dictated by size-at-age. After transformation, new juveniles were added to the single, panmictic lake-inhabiting population. After maturation, adults were assigned back to larval populations (i.e., tributaries) based on the relative abundances of larvae currently inhabiting each tributary. That is, in accordance with their life history, more adults were assigned back to large larval populations than to small larval populations with the following equation:
where \({N}_{t}\) equals the total number of adults assigned back to tributary t, \({l}_{i}\) equals the total number of larvae in tributary t and \({N}_{adults}\) equals the total number of adults ready to spawn. This step was implemented to mimic the observation that adult sea lamprey, unlike salmon and many other anadromous fishes, do not return to their natal streams (see Supporting Information for detailed life history information). After randomly allocating adults to tributaries based on larval population sizes, males and females were paired and allowed to reproduce. Reproduction was modeled as both density dependent and density independent (though the default was density dependent; Table 1, Fig. 1) and the number of offspring per pair was set as:
where \({K}_{t}\) equals the tributary specific carrying capacity, \({N}_{t}\) was calculated with Eq. (2) and \({n}_{m}\) and \({n}_{a}\) are constants that control the strength of density dependence. Previous empirical work in this system has shown that Great Lakes sea lamprey have highly compensatory population dynamics44, and so our default parameters assumed that egg and larval survival through the first year was also density dependent (Table 1). During reproduction each pair created \({N}_{offspring}\) with each offspring inheriting, in Mendelian fashion, one randomly selected allele from each of their parents to create genotypes that determine resistance status (see details below). Variance in the number of offspring per pair was introduced with subsequent larval mortality. After reproduction, all adults were removed from the model simulating the strict semelparity found in adult sea lamprey.
Modeling larval mortality was the next step. This natural mortality occurred regardless of resistance status. For each tributary, a random deviate from a normal distribution was used to determine the number of surviving larvae. The mean of the normal distribution was equal to the carrying capacity for that particular tributary and the default standard deviation was set to 70 allowing for both density dependent and density independent mortality. Juvenile mortality was procedurally identical to larval mortality except that mortality was applied to the juvenile (parasitic, lake-inhabiting) population (substituting a lake-wide carrying capacity, Table 1). The steps (growth, transformation, maturation, reproduction, and mortality) were repeated for the first 50 years of each simulation to allow the age and stage structure to fully develop. At year 50, TFM was applied after the natural larval mortality step (Fig. 1). A parameter within the model allowed us to set the number of tributaries to treat with TFM each year. Throughout the Great Lakes, larval abundance surveys are conducted annually, and the tributaries are rank ordered; tributaries with highest lamprey abundance are often given priority for treatment14. Similarly, in our model we rank ordered each population by the number of larvae and treated the top ranked tributaries. For example, if we set the number of treated tributaries per year to 5, we would treat the 5 largest larval populations. We assumed that TFM mortality was 0.934, meaning that 90% of susceptible larvae within a tributary would die after exposure to TFM. The combination of both the number of treated tributaries and the TFM mortality rate determined the total number of larval lamprey killed by TFM each year.
Resistance was modeled as a single locus trait as small numbers of loci of large effect have dominated resistance studies across diverse taxa45. Here, resistance was modeled as a dominant trait; individuals heterozygous for the resistance allele were also conferred full resistance though they only passed on the resistance allele to half of their offspring. There were no resistant individuals at the beginning of each simulation. At year 70, a single resistant adult was added to the population of juveniles that had just transitioned to the adult stage. The resistant individual was added 20 years after the start of TFM to allow us to disentangle the effects of treatment and resistance, but the actual year a resistant adult was added had no effect on the results (e.g., resistance developed x years after the introduction of a resistant adult regardless of whether it was added at year 50, 70, or 90). We also added a parameter to examine costs associated with resistance, which are often found in empirical studies of resistance19. Cost of resistance was modeled as a decrease in reproductive success. For example, if the parameter was set to 0.1, then a resistant individual would have a 10% reduction in reproductive success (here, set as the number of offspring) relative to a non-resistant individual.
Parameter exploration
To test the effect of TFM treatment on the development of resistance we first varied the number of treated tributaries holding all other parameters constant (default values; Table 1). We examined the effects of applying TFM each year to 5, 10, 20, and 30% of all tributaries. For each year, we measured the number of parasitic juveniles and the number of resistant larvae. To ensure that the results scaled with the larger absolute population sizes observed in the field, we re-ran these analyses setting larval and juvenile carrying capacities to be an order of magnitude larger than the default vales (Table 1). We next examined the effects of natal homing on the spread of resistance using the default parameters (Table 1) except for adult migration. We first examined a linear stepping stone model of migration46 and varied gene flow of adults among natal tributaries from 1 to 20% per year. This alternative stepping stone scenario was run in addition to the default no-homing migration pattern unique to sea lamprey. We kept track of the number and geographic spread of resistant larvae and, for each level of gene flow, measured the proportion of 100 replicates that had at least one resistant larva at model year 200.
We next examined the costs of resistance by modeling three costs of resistance (0, 0.1, and 0.4 reductions in fitness) and stopping treatment at year 100 (30 years after the introduction of a resistant adult). By stopping treatment, we could directly examine how quickly resistant individuals were removed from the population after the survival advantage associated with resistance was removed. To explore costs of resistance in more detail, we varied the cost of resistance from 0 to 0.5 in intervals of 0.1 and measured the proportion of resistant larvae as a function of sample year. We first assumed a perfect detection probability (i.e., all resistant lamprey were successfully identified), but later relaxed this assumption to examine the relationship between sampling effort (% of total larvae sampled), cost of resistance, and probability of detection. Lastly, we varied the strength of density-dependent reproduction by varying the parameter \({n}_{m}\) in Eq. 3 (see Supporting Information for details). The model was developed in R 3.5.047 and run on a high-performance computing cluster. For every unique combination of parameters, 100 independent simulations were performed.
Results
Treating tributaries with TFM effectively reduces the number of parasitic juveniles found in the system, however, this effectiveness is directly tied to the number of larval populations treated each year. When only 5% of tributaries are treated each year, there is only a small reduction in the number of parasitic juveniles (Fig. 2A). Conversely, when 30% of the tributaries are treated with TFM each year there was an 81% reduction in the total number of parasitic juveniles (Fig. 2G). This result was consistent regardless of the total number of individuals or tributaries included in the model (Figs. S1; S2). We also found that the total number of resistant larvae at year 200 (the last year in our model simulations) varied with the total number of treated tributaries (Fig. 2B,D,F,H). Under the high selection intensity of 30% of tributaries treated, resistance took an average of 30 years to develop. Under more moderate selection intensities, resistance took between 51 and 109 years to develop (Fig. 2D,F).
The proportion of simulations that did not develop resistance was higher when the number of treated tributaries was lower (Fig. 3). In other words, when selection intensity was low, the likelihood of successful reproduction by the single resistant adult introduced at year 70 was reduced (Fig. 3). When resistance does evolve, the total (i.e., cumulative) number of parasites in the system over time is equivalent between heavily-treated systems where resistance often evolves and lightly-treated systems where resistance does not evolve as often (Fig. 3). Combined, these results illustrate that a low selection intensity regime may result in the fewest parasitic lamprey over the long term because: 1. resistance is less likely to evolve (Fig. 3A), and 2. when resistance does evolve it takes a much longer time for the number of resistant larvae to increase in frequency (Fig. 2B). Thus, while a high selection intensity treatment regime will result in fewer parasitic lamprey over the short term, a low selection regime can ultimately result in fewer total lamprey in the system over the long term.
We also found that the sea lamprey life history strategy for identifying spawning habitat by cuing in on larval abundance greatly increases both the speed at which resistant individuals spread throughout the system and the proportion of simulations where resistance developed. For example, 15 years after the introduction of a resistant adult, resistant larvae had spread to 70% of larval tributaries when sea lamprey adult abundance was dictated by larval abundance (Fig. 4C). By contrast, if larvae followed an isolation-by-distance pattern of dispersal, then only 5% of larval populations would be colonized 15 years after the introduction of a resistant adult (Fig. 4B). Furthermore, even if gene flow among neighboring populations was as high as 20% per year (which would be high for most vertebrates), the proportion of simulations where resistance developed remained considerably lower than the actual migration strategy employed by lamprey (48% vs. 82%; Fig. 4D). These results illustrate that the lack of natal philopatry in invasive sea lamprey greatly increase the probability of resistance spreading throughout the entire system.
The lack of natal philopatry also makes the early detection of resistance challenging. Assuming that a resistant larva could be detected in a given sample, an immense sampling and screening effort is required to detect resistance in its early stages (e.g., 10 or 20 years after the introduction of a single resistant adult) (Fig. 5). For example, if 0.5% of all larvae were sampled 20 years after the release of a single resistant adult, the probability of the sample including a single resistant larva would only equal 50%. Given that larval abundances are on the order of tens millions of individuals, a 0.5% sample size would involve screening a minimum of fifty thousand individuals. Furthermore, the probability of detection decreases as the costs of resistance increase (Fig. S3).
We next examined the costs of resistance in this system by penalizing the fecundity and thus the reproductive success of resistant adult lamprey. To examine the consequences on larval abundance, we stopped TFM treatment at year 100. We found that when there was no cost of resistance, the total number of resistant larvae remained relatively constant for the next 100 years (Fig. 6A). From a pragmatic standpoint, this result means that if there is no cost of resistance, then resistant genotypes cannot be easily purged from the system even in the absence of treatment. If there is a moderate cost to resistance (i.e., 10%), then the number of resistant larvae will decline through time in the absence of treatment, but not in a way that is meaningful from a management perspective (i.e., it can take 100 years or more for all resistant larvae to be eliminated from the system) (Fig. 6B). By contrast, if there is a high cost to resistance (40%), then stopping treatment could be an effective management strategy and resistant larvae can be eliminated from the system in as little as 10 years (Fig. S4).
We also examined how costs of resistance influence both the probability and the time it takes for resistance to develop in the system. We found that as the cost of resistance increases, the proportion of resistant larvae decreases (Fig. 6C). We next illustrate a strong tradeoff between the proportion of parasites killed by TFM and the probability of resistance evolving (Fig. 6D). As selection intensity increases, the proportion of the parasite population that is initially reduced increases, but so does the probability of resistance evolving (Fig. 6D). As the cost of resistance increases, the probability of resistance evolving decreases, but only up to a point (see also Fig. S5). Thus, costs of resistance can reduce, but not eliminate, both the likelihood and speed at which resistance develops in this system.
Discussion
Our results indicate that sea lamprey, an invasive fish species successfully treated with a pesticide for over 60 years in the Great Lakes, could evolve resistance between 40 and 80 years after the introduction of a resistant individual. The single largest determinant of whether resistance takes a comparatively short (e.g., 30 years) or longer time (e.g., 80 years) to develop is the annual number of tributaries treated. Across all Great Lakes, sea lamprey have been treated with TFM for an average of 50 years34, which is towards the middle of when our model predicts resistance will evolve. As such, sea lamprey in the Great Lakes could be in the incipient stages of resistance evolving. If alternative control measures can’t be found, the spread of resistance in Great Lakes sea lamprey could threaten the viability of one of the world’s most successful invasive vertebrate control programs and cause substantial declines in economically valuable commercial and recreational fish populations. Compared to most other systems studied to date, resistance in the lamprey-TFM system can be viewed as taking a relatively long time to develop (as measured in years, not generations), but rapidly accelerates after a certain threshold is reached (Fig. 2). The relatively long time it takes for resistance to develop is because: i. sea lamprey can spend up to seven years in the larval stage before transforming into a parasitic adult and ii. the lack of natal philopatry results in resistant lamprey remaining at low abundance while spreading throughout the entire system. These life history characteristics mean that resistant larvae and adults initially increase in abundance very slowly in the system.
This slow initial accumulation of resistant individuals also means that it can be very challenging to detect resistance in its incipient stages. In fact, it can take 40 years or more before the proportion of resistant larvae increases above 50% (Fig. 6C). Ten years after the introduction of a resistant adult, our model predicts that, for a moderate selection intensity, screening 0.5% of the larval population for resistance would only result in a 10% probability of including a single resistant individual in the sample. Because sea lamprey are highly fecund, there can be tens of millions of larvae in the population at any given point in time (GLFC control board, personal communication). Thus, a screening procedure that sampled a minimum of 50,000 larvae would only have around a 10% probability of having at least one resistant larvae 10 years after the introduction of a resistant adult (Fig. 5). Our results also assume perfect detection; we assume that if a resistant individual was included in the sample then it would be correctly identified as such. In reality, type I errors (classifying a non-resistant individual as resistant) and type II errors (failing to detect a resistant individual when it is present in the sample) would jointly affect the probability of successful detection. Technical issues aside, the fact remains that a very large sample of larvae must be screened to detect resistance during its early stages. This detection challenge presents a quandary; how much effort and resources should be devoted to screening for resistance versus developing new control measures?
Additional evidence for the evolution of resistance can be obtained by tracking changes through time in the lethal concentrations of a pesticide required to kill fixed percentages of test organisms. In the Great Lakes, the lethal concentration of TFM required for 50% mortality (LC50) of sea lamprey ranges from 1 to 3.5 mg/L34. As resistant genotypes increase in frequency, the prediction is that the LC50 would increase through time. Taking this approach, we recently reviewed the outcomes of historical TFM toxicity tests that have regularly been conducted on sea lamprey larvae since 195634. Over the 57-year period of toxicity testing, no increase in the TFM concentration at which 50% of the sample lamprey died (the LC50 over 19–24 h;) was detected (Fig. 2 in34). However, these findings should be interpreted cautiously when drawing conclusions about the evidence for resistance. The analyses relied on historical data that were collected for other purposes, not with the intent of specifically investigating the evolution of resistance and, given the small sample sizes, the ability to detect resistance would be limited to all but the latest stages of resistance. Appropriately designed experiments that couple toxicology, gene expression, and populations with varying durations of exposure are needed to rule out the possibility of resistance in sea lamprey.
Including a fitness cost for resistance decreased the time for resistance to be eradicated from the system via natural selection once control with TFM ceased. However, resistant individuals were eliminated from the system within a short enough time frame acceptable for management (e.g., 20 years or less) only if the cost of resistance was very high (e.g., a 40% reduction in lifetime reproductive success). We further found that the onset of resistance is delayed and the total number of simulations that developed resistance decreased with increasing costs of resistance. However, regardless of the costs of resistance, there was a strong trade-off between the proportion of parasites initially controlled by TFM and the probability of resistance evolving (Fig. 6D). The relative roles of density dependence and density independence in resistance evolution has often been ignored (but see48,49). Here, we show that as the relative strength of density dependent reproduction decreases, the proportion of simulations that develop resistance decreases and, for those that do develop resistance, the number of years until resistance develops increases (Fig. S6). Although, it is believed that sea lamprey in the Great Lakes are highly regulated by compensatory dynamics44, robust empirical data regarding density-dependent dynamics during reproduction and shortly thereafter are largely lacking. Thus, both high costs of resistance and reduced density dependent reproduction can delay, but not entirely prevent, the evolution of resistance.
Caveats, assumptions, and management recommendations
Our model results suggest that the onset of resistance is either ongoing or imminent for invasive Great Lakes sea lamprey that have been treated with TFM for over 60 years. There are, however, a few caveats worth considering. First, it is possible that the total percentage of larval habitat treated with TFM each year is overestimated. Adults may spawn in undocumented areas and there could be large areas of larval habitat that go unsurveyed (e.g., small reaches of tributaries, lentic bays). If the percentage of treated larval habitat is smaller than expected, the overall mortality rate imposed on the population by TFM could be lower than expected which would reduce the selection intensity (Figs. 2 and 6). However, the ability of the control program to keep adult sea lamprey numbers at relatively low levels implies a mortality rate high enough to impose some selective pressure. Second, we have assumed that resistance will evolve via genes of large effect, the most common form of resistance documented in the literature50,51. However, it is possible that resistance could evolve more along the lines of a quantitative trait with hundreds of loci of smaller effect. If this were to happen, then the response to selection is expected to be much slower52, and the onset of resistance could be delayed by tens of generations. Of course, resistance can evolve in multiple ways and at any point in time, especially given the lamprey’s high fecundity, and so the evolution of a different class of resistant individuals with a gene of large effect may be more likely to occur before the quantitative trait resistance alleles change in frequency. Third, we don’t consider the effects of prey abundance (or other factors like climate) on the sea lamprey population size. Given that sea lamprey population sizes do not affect the timing or magnitude of resistance (Figs. S1, S2), we do not expect to see much of an effect for prey availability on resistance evolution. Nonetheless, small effective populations may limit the response to selection (but see Wood, et al.53) or there could be unaccounted eco-evolutionary feedbacks such as how energy reserves affect fecundity or timing of reproduction. Lastly, it is possible that a severe founder effect or a unique mode of action for TFM may limit the likelihood of an evolutionary response. We think these possibilities are unlikely. Microsatellite data suggest that genetic diversity in Great Lakes sea lamprey is not substantially lower than lamprey in their native range38. Furthermore, because the mode of action for TFM largely targets ATP synthesis54 there are many potential modes of resistance evolution34.
Considering these results, we have several management recommendations. First, we suggest that alternative control methods with different modes of action and high target specificity be developed as quickly as possible. In the meantime, one option is to reduce the strength of selection by lowering the number of treated streams. This decision would have to be carefully considered against the shorter-term economic and ecological costs of potentially higher numbers of parasites. Through a risk assessment and cost-benefits analysis, it should be possible to identify a target level of streams to treat that maximizes the tradeoff between parasite abundance and resistance evolution (e.g., Fig. 3). Second, it would be prudent to develop a high throughput screening program aimed at detecting resistant larvae. Finally, there are several important research gaps that, if addressed, would enable a more accurate assessment of the timelines of resistance evolution in this system and enable more effective management55; these include quantifying the degree of density-dependent reproduction, characterizing fine-scale genetic spatial structure, measuring the true strength of selection, and studying potential costs of resistance. Appropriate action to impede the onset of this incipient crisis should remain a high priority for the management and conservation of Great Lakes fishes. Such action could not only prevent the loss of billions of dollars from the economy but could also prevent large-scale and undesirable ecological shifts in abundance and community composition.
Data availability
Model code available at https://github.com/ChristieLab.
References
Davies, J. & Davies, D. Origins and evolution of antibiotic resistance. Microbiology and molecular biology reviews 74, 417–433 (2010).
Neu, H. C. The crisis in antibiotic resistance. Science 257, 1064–1073 (1992).
Ma, Z. & Michailides, T. J. Advances in understanding molecular mechanisms of fungicide resistance and molecular detection of resistant genotypes in phytopathogenic fungi. Crop Protection 24, 853–863 (2005).
Roush, R. T. & McKenzie, J. A. Ecological genetics of insecticide and acaricide resistance. Annual review of entomology 32, 361–380 (1987).
Hemingway, J. & Ranson, H. Insecticide resistance in insect vectors of human disease. Annual review of entomology 45, 371–391 (2000).
Jasieniuk, M., Brûlé-Babel, A. L. & Morrison, I. N. The evolution and genetics of herbicide resistance in weeds. Weed science, 176–193 (1996).
Owen, M. D. & Zelaya, I. A. Herbicide-resistant crops and weed resistance to herbicides. Pest management science 61, 301–311 (2005).
Thijssen, H. H. Warfarin-based rodenticides: mode of action and mechanism of resistance. Pesticide science 43, 73–78 (1995).
Kohn, M. H., Pelz, H.-J. & Wayne, R. K. Natural selection mapping of the warfarin-resistance gene. Proceedings of the National Academy of Sciences 97, 7911–7915 (2000).
Twigg, L. E., Martin, G. R. & Lowe, T. J. Evidence of pesticide resistance in medium-sized mammalian pests: a case study with 1080 poison and Australian rabbits. Journal of Applied Ecology 39, 549–560 (2002).
Lohmeyer, A. M. & Garvey, J. E. Placing the North American invasion of Asian carp in a spatially explicit context. Biological Invasions 11, 905–916 (2009).
Mooney, H. A. & Cleland, E. E. The evolutionary impact of invasive species. Proceedings of the National Academy of Sciences 98, 5446–5451 (2001).
Pimentel, D., Zuniga, R. & Morrison, D. Update on the environmental and economic costs associated with alien-invasive species in the United States. Ecological economics 52, 273–288 (2005).
Christie, G. C. & Goddard, C. I. Sea Lamprey International Symposium (SLIS II): Advances in the Integrated Management of Sea Lamprey in the Great Lakes. Journal of Great Lakes Research 29, 1–14, https://doi.org/10.1016/S0380-1330(03)70474-2 (2003).
Purvis, H. A. Variations in growth, age at transformation, and sex ratio of sea lampreys reestablished in chemically treated tributaries of the Upper Great Lakes. (Great Lakes Fishery Commission Ann Arbor, MI, 1979).
Youson, J. H. The biology of metamorphosis in sea lampreys: endocrine, environmental, and physiological cues and events, and their potential application to lamprey control. Journal of Great Lakes Research 29, 26–49 (2003).
Whalon, M., Mota-Sanchez, D. & Hollingworth, R. Analysis of global pesticide resistance in arthropods, pp. 5–31. (CABI International, United Kingdom, 2008).
Guo, Z., Andreou, D. & Britton, J. Sea lamprey Petromyzon marinus biology and management across their native and invasive ranges: promoting conservation by knowledge transfer. Reviews in Fisheries Science & Aquaculture 25, 84–99 (2017).
Coustau, C. & Chevillon, C. Resistance to xenobiotics and parasites: can we count the cost? Trends in Ecology & Evolution 15, 378–383 (2000).
Lawrie, A. The sea lamprey in the Great Lakes. Transactions of the American Fisheries Society 99, 766–775 (1970).
Larson, G. L. et al. The history of sea lamprey control in Lake Ontario and updated estimates of suppression targets. Journal of Great Lakes Research 29, 637–654 (2003).
Applegate, V. C., Howell, J. H., Hall, A. E. & Smith, M. A. Toxicity of 4346 chemicals to larval lamprey and fishes. U.S. Fish Wildl. Serv., Spec. Rep., Fish 207 (1957).
Smith, B. & Tibbles, J. Sea lamprey (Petromyzon marinus) in Lakes Huron, Michigan, and Superior: history of invasion and control, 1936–78. Canadian Journal of Fisheries and Aquatic Sciences 37, 1780–1801 (1980).
Sawyer, A. J. Prospects for integrated pest management of the sea lamprey (Petromyzon marinus). Canadian Journal of Fisheries and Aquatic Sciences 37, 2081–2092 (1980).
Sorensen, P. & Hoye, T. A critical review of the discovery and application of a migratory pheromone in an invasive fish, the sea lamprey Petromyzon marinus L. Journal of Fish Biology 71, 100–114 (2007).
Siefkes, M. J., Steeves, T. B., Sullivan, W. P., Twohey, M. B. & Li, W. Sea lamprey control: past, present, and future. Great Lakes Fisheries Policy and Management. Michigan State University Press, East Lansing, 651–704 (2013).
McLaughlin, R. L. et al. Unintended consequences and trade-offs of fish passage. Fish and Fisheries 14, 580–604 (2013).
Georghiou, G. P. & Taylor, C. E. Factors influencing the evolution of resistance. Pesticide resistance: strategies and tactics for management, 157–169 (1986).
Kogan, M. Integrated pest management: historical perspectives and contemporary developments. Annual review of entomology 43, 243–270 (1998).
Manion, P. J. Fecundity of the sea lamprey (Petromyzon marinus) in Lake Superior. Transactions of the American Fisheries Society 101, 718–720 (1972).
Beamish, F. Biology of the North American anadromous sea lamprey, Petromyzon marinus. Canadian Journal of Fisheries and Aquatic Sciences 37, 1924–1943 (1980).
Johnson, N. S. et al. Survival and metamorphosis of low-density populations of larval sea lampreys (Petromyzon marinus) in streams following lampricide treatment. Journal of Great Lakes Research 40, 155–163 (2014).
Renaud, C. B. Lampreys of the world. An annotated and illustrated catalogue of lamprey species known to date. (Food and Agriculture Organization of the United Nations, 2011).
Dunlop, E. S. et al. Rapid evolution meets invasive species control: the potential for pesticide resistance in sea lamprey. Canadian Journal of Fisheries and Aquatic Sciences 75, 152–168 (2017).
Christie, M. R., Marine, M. L., French, R. A. & Blouin, M. S. Genetic adaptation to captivity can occur in a single generation. Proceedings of the National Academy of Sciences 109, 238–242 (2012).
Whitehead, A., Pilcher, W., Champlin, D. & Nacci, D. Common mechanism underlies repeated evolution of extreme pollution tolerance. Proceedings of the Royal Society of London B: Biological Sciences 279, 427–433 (2012).
Bjerselius, R. et al. Direct behavioral evidence that unique bile acids released by larval sea lamprey (Petromyzon marinus) function as a migratory pheromone. Canadian Journal of Fisheries and Aquatic Sciences 57, 557–569 (2000).
Bryan, M. et al. Patterns of invasion and colonization of the sea lamprey (Petromyzon marinus) in North America as revealed by microsatellite genotypes. Molecular Ecology 14, 3757–3773 (2005).
Spencer, J. L. & Raghu, S. Refuge or reservoir? The potential impacts of the biofuel crop Miscanthus x giganteus on a major pest of maize. PLoS One 4, e8336 (2009).
Dunlop, E. S., Heino, M. & Dieckmann, U. Eco-genetic modeling of contemporary life‐history evolution. Ecological Applications 19, 1815–1834 (2009).
Von Bertalanffy, L. A quantitative theory of organic growth (inquiries on growth laws. II). Human biology 10, 181–213 (1938).
Hansen, M. J. et al. Population ecology of the sea lamprey (Petromyzon marinus) as an invasive species in the Laurentian Great Lakes and an imperiled species in Europe. Reviews in fish biology and fisheries 26, 509–535 (2016).
Manion, P. J. & Hanson, L. H. Spawning behavior and fecundity of lampreys from the upper three Great Lakes. Canadian Journal of Fisheries and Aquatic Sciences 37, 1635–1640 (1980).
Jones, M. L. et al. Compensatory mechanisms in Great Lakes sea lamprey populations: implications for alternative control strategies. Journal of Great Lakes Research 29, 113–129 (2003).
Li, X., Schuler, M. A. & Berenbaum, M. R. Molecular mechanisms of metabolic resistance to synthetic and natural xenobiotics. Annu. Rev. Entomol. 52, 231–253 (2007).
Kimura, M. & Weiss, G. H. The stepping stone model of population structure and the decrease of genetic correlation with distance. Genetics 49, 561 (1964).
Team, R. C. R: A Language and Environment for Statistical Computing. Vienna, Austria, https://www.R-project.org/ (2018).
Churcher, T. S. & Basáñez, M. G. Density dependence and the spread of anthelmintic resistance. Evolution: International Journal of organic evolution 62, 528–537 (2008).
Baskett, M. L. & Waples, R. S. Evaluating alternative strategies for minimizing unintended fitness consequences of cultured individuals on wild populations. Conservation Biology 27, 83–94 (2013).
Taylor, M. & Feyereisen, R. Molecular biology and evolution of resistance of toxicants. Molecular biology and evolution 13, 719–734 (1996).
Feyereisen, R. Molecular biology of insecticide resistance. Toxicology letters 82, 83–90 (1995).
Lynch, M. & Walsh, B. Genetics and analysis of quantitative traits. Vol. 1 (Sinauer Sunderland, MA, 1998).
Wood, J. L., Yates, M. C. & Fraser, D. J. Are heritability and selection related to population size in nature? Meta-analysis and conservation implications. Evolutionary applications 9, 640–657 (2016).
Birceanu, O., McClelland, G. B., Wang, Y. S., Brown, J. C. & Wilkie, M. P. The lampricide 3-trifluoromethyl-4-nitrophenol (TFM) uncouples mitochondrial oxidative phosphorylation in both sea lamprey (Petromyzon marinus) and TFM-tolerant rainbow trout (Oncorhynchus mykiss). Comparative Biochemistry and Physiology Part C: Toxicology & Pharmacology 153, 342–349 (2011).
Fussell, K. M. D. et al. A perspective on needed research, modeling, and management approaches that can enhance Great Lakes fisheries management under changing ecosystem conditions. Journal of Great Lakes Research 42, 743–752, https://doi.org/10.1016/j.jglr.2016.04.007 (2016).
Acknowledgements
We thank A. Martinez, A. Muir, C. Searle, and M. Siefkes for feedback and assistance with this project. This research was facilitated by the Rosen Center for Advanced Computing at Purdue, West Lafayette, IN. This project was funded by support to MRC and MSS from the Great Lakes Fishery Commission (Project ID: 2016_CHR_54053).
Author information
Authors and Affiliations
Contributions
M.R.C., M.S.S. and E.S.D. developed the conceptual framework, guided model development, and wrote the manuscript. M.R.C. wrote the model and analyzed the data.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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 images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Christie, M.R., Sepúlveda, M.S. & Dunlop, E.S. Rapid resistance to pesticide control is predicted to evolve in an invasive fish. Sci Rep 9, 18157 (2019). https://doi.org/10.1038/s41598-019-54260-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-54260-5
This article is cited by
-
The effect of putrescine on space use and activity in sea lamprey (Petromyzon marinus)
Scientific Reports (2022)
-
Reduced thiamine utilization by Seneca Lake lake trout embryos and potential implications to restoration of lake trout in the Great Lakes
Environmental Biology of Fishes (2021)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.