Using Agent-Based Modelling to Predict the Role of Wild Refugia in the Evolution of Resistance of Sea Lice to Chemotherapeutants

A modeling study of the impact of treatment policies on the evolution of resistance in sea lice on salmon farms

Salmonid aquaculture is an important source of nutritious food with more than 2 million tonnes of fish produced each year (Food and Agriculture Organisation of the United Nations, 2019). In most salmon producing countries, sea lice represent a major barrier to the sustainability of salmonid aquaculture. This issue is exacerbated by widespread resistance to chemical treatments on both sides of the Atlantic. Regulation for sea lice management mostly involves reporting lice counts and treatment thresholds, which depending on interpretation may encourage preemptive treatments. We have developed a stochastic simulation model of sea lice infestation including the lice life-cycle, genetic resistance to treatment, a wildlife reservoir, salmon growth and stocking practices in the context of infestation, and coordination of treatment between farms. Farms report infestation levels to a central organisation, and may then cooperate or not when coordinated treatment is triggered. Treatment practice then impacts the level of resistance in the surrounding sea lice population. Our simulation finds that treatment drives selection for resistance and coordination between managers is key. We also find that position in the hydrologically-derived network of farms can impact individual farm infestation levels and the topology of this network can impact overall infestation and resistance. We show how coordination and triggering of treatment alongside varying hydrological topology of farm connections affects the evolution of lice resistance, and thus optimise salmon quality within socio-economic and environmental constraints. Network topology drives infestation levels in cages, treatments, and hence treatment-driven resistance. Thus farmer behaviour may be highly dependent on hydrologically position and local level of infestation.

Patterns of asexual reproduction of the soybean aphid, Aphis glycines (Matsumura), with and without the secondary symbionts Wolbachia and Arsenophonus, on susceptible and resistant soybean genotypes

Plant breeding is used to develop crops with host resistance to aphids, however, virulent biotypes often develop that overcome host resistance genes. We tested whether the symbionts, Arsenophonus (A) and Wolbachia (W), affect virulence and fecundity in soybean aphid biotypes Bt1 and Bt3 cultured on whole plants and detached leaves of three resistant, Rag1, Rag2 and Rag1 + 2, and one susceptible, W82, soybean genotypes. Whole plants and individual aphid experiments of A. glycines with and without Arsenophonus and Wolbachia did not show differences in overall fecundity. Differences were observed in peak fecundity, first day of deposition, and day of maximum nymph deposition of individual aphids on detached leaves. Bt3 had higher fecundity than Bt1 on detached leaves of all plant genotypes regardless of bacterial profile. Symbionts did not affect peak fecundity of Bt1 but increased it in Bt3 (A+W+) and all Bt3 strains began to deposit nymphs earlier than the Bt1 (A+W-). Arsenophonus in Bt1 delayed the first day of nymph deposition in comparison to aposymbiotic Bt1 except when reared on Rag1 + 2. For the Bt1 and Bt3 strains, symbionts did not result in a significant difference in the day they deposited the maximum number of nymphs nor was there a difference in survival or variability in number of nymphs deposited. Variability of number of aphids deposited was higher in aphids feeding on resistant plant genotypes. The impact of Arsenophonus on soybean aphid patterns of fecundity was dependent on the aphid biotype and plant genotype. Wolbachia alone had no detectable impact but may have contributed to the increased fecundity of Bt3 (A+W+). An individual based model, using data from the detached leaves experiment and with intraspecific competition removed, found patterns similar to those observed in the greenhouse and growth chamber experiments including a significant interaction between soybean genotype and aphid strain. Combining individual data with the individual based model of population growth isolated the impact of fecundity and host resistance from intraspecific competition and host health. Changes to patterns of fecundity, influenced by the composition and concentration of symbionts, may contribute to competitive interactions among aphid genotypes and influence selection on virulent aphid populations.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Parasite management in aquaculture exerts selection on salmon louse behaviour

The evolution of pest resistance to management strategies is a major challenge for farmed systems. Mitigating the effects of pest adaptation requires identifying the selective pressures imposed by these strategies. In Atlantic salmon (Salmo salar) aquaculture, barriers are used to prevent salmon louse (Lepeophtheirus salmonis) larvae (copepodids) from entering salmon cages. These barriers are effective against shallow-swimming copepodids, but those swimming deeper can pass underneath and infest salmon. Laboratory experiments suggest that depth regulation in copepodids is a variable behavioural trait with a genetic basis. We used biological-hydrodynamic dispersal models to assess how this trait variation alters the dispersion of lice through the ocean environment and into farms. The dispersal of copepodids with 3 behavioural phenotypes (deep, mean or shallow) was modelled over winter-spring and spring-summer periods in a Norwegian fjord system with intensive aquaculture. The infestation pressure of each phenotype on barrier cages was estimated from their modelled depth distributions: copepodids deeper than 10 m were predicted to successfully pass underneath barriers. The deep phenotype was the most abundant below 10 m and reached infestation pressures 3 times higher than that of the mean phenotype. In contrast, the shallow phenotype infestation pressure reached less than half that of the mean phenotype. These differences in relative fitness indicate that barriers can impose strong directional selection on the swimming behaviour of copepodids. The strength of this selection varied seasonally and geographically, with selection for the deep phenotype stronger in winter-spring and at coastal locations than in spring-summer and within fjords. These findings can be applied across farms to slow louse adaptation, by limiting barriers during situations of strong selection, although this must be balanced against trade-offs to short-term efficacy. More broadly, our study highlights new ways in which dispersal models can address evolutionary questions crucial for sustainable parasite management in aquaculture.

The potential for cleaner fish-driven evolution in the salmon louse Lepeophtheirus salmonis: Genetic or environmental control of pigmentation?

The parasitic salmon louse represents one of the biggest challenges to environmentally sustainable salmonid aquaculture across the globe. This species also displays a high evolutionary potential, as demonstrated by its rapid development of resistance to delousing chemicals. In response, farms now use a range of non-chemical delousing methods, including cleaner fish that eat lice from salmon. Anecdotal reports suggest that in regions where cleaner fish are extensively used on farms, lice have begun to appear less pigmented and therefore putatively less visible to cleaner fish. However, it remains an open question whether these observations reflect a plastic (environmental) or adaptive (genetic) response. To investigate this, we developed a pigment scoring system and conducted complimentary experiments which collectively demonstrate that, a) louse pigmentation is strongly influenced by environmental conditions, most likely light, and b) the presence of modest but significant differences in pigmentation between two strains of lice reared under identical conditions. Based on these data, we conclude that pigmentation in the salmon louse is strongly influenced by environmental conditions, yet there are also indications of underlying genetic control. Therefore, lice could display both plastic and adaptive responses to extensive cleaner fish usage where visual appearance is likely to influence survival of lice.

Migratory hosts can maintain the high-dose/refuge effect in a structured host-parasite system: The case of sea lice and salmon

Migration can reduce parasite burdens in migratory hosts, but it connects populations and can drive disease dynamics in domestic species. Farmed salmon are infested by sea louse parasites, often carried by migratory wild salmonids, resulting in a costly problem for industry and risk to wild populations when farms amplify louse numbers. Chemical treatment can control lice, but resistance has evolved in many salmon-farming regions. Resistance has, however, been slow to evolve in the north-east Pacific Ocean, where large wild-salmon populations harbour large sea louse populations. Using a mathematical model of host-macroparasite dynamics, we explored the roles of domestic, wild oceanic and connective migratory host populations in maintaining treatment susceptibility in associated sea lice. Our results show that a large wild salmon population, unexposed to direct infestation by lice from farms; high levels of on-farm treatment; and a healthy migratory host population are all critical to slowing or stopping the evolution of treatment resistance. Our results reproduce the "high-dose/refuge effect," from the agricultural literature, with the added requirement of a migratory host population to maintain treatment susceptibility. This work highlights the role that migratory hosts may play in shared wildlife/livestock disease, where evolution can occur in ecological time.

Experimental transmission of Caligus rogercresseyi between two different fish species

Caligus rogercresseyi is the dominant sea louse parasite affecting the salmon and trout industry in southern Chile. This parasite has a wide range of native and endemic fish hosts. The Patagonian blenny Eleginops maclovinus, which is parasitized mostly by the caligid species Lepeophtheirus spp. and C. rogercresseyi, is presumably responsible for the transmission of C. rogercresseyi to salmonids. The aim of this study was to characterize the transmission of parasites between different fish species and parasite cohort development under laboratory conditions. Parasite abundances and intensities were quantified. Transmission of parasites from Patagonian blenny to Atlantic salmon Salmo salar was lower (~9%, mainly corresponding to C. rogercresseyi) than from salmon to Patagonian blenny (14.7-26.9%, where only C. rogercresseyi were observed). This suggests that the transmission of C. rogercresseyi from salmon individuals is higher than the transmission from a native fish. Parasite cohorts developed successfully on both fish species, but apparently under different developmental rates. Water temperature, oxygen, and juvenile abundances were the variables that better explained cohort development success and variation in C. rogercresseyi adult abundances over time.

Sea lice exposure to non-lethal levels of emamectin benzoate after treatments: a potential risk factor for drug resistance

The avermectin derivative emamectin benzoate (EMB) has been widely used by salmon industries around the world to control sea lice infestations. Resistance to this anti-parasitic drug is also commonly reported in these industries. The objective of this study was to quantify the number of sea lice potentially exposed to sub-lethal concentrations of EMB while fish clear the drug after treatments. We assessed juvenile sea lice abundance after 38 EMB treatments on six Atlantic salmon farms, in a small archipelago in British Colombia, Canada, between 2007 and 2018. We fitted a standard EMB pharmacokinetic curve to determine the time when fish treated with this product would have EMB tissue concentrations below the recommended target therapeutic level. During the study, we estimated that for each sea lice treatment there was, on average, an abundance of 0.12 juvenile sea lice per fish during the time period when the concentrations of EMB would have been lower than 60ppb, the recommended therapeutic treatment level for sea lice. The findings from this study on metaphylactic anti-parasitic treatments identify a potential driver for drug resistance in sea lice that should be further explored.

Susceptibility of Caligus rogercresseyi collected from the native fish species Eleginops maclovinus (Cuvier) to antiparasitics applied by immersion

A major challenge for Chilean salmon farming is infestation by the ectoparasite Caligus rogercresseyi. In addition, there is evidence that a loss of chemotherapeutic treatment efficacy against important fish pathogens is occurring in salmon farming, including antiparasitic efficacy. Currently, there are known techniques that allow the determination of the susceptibility profile of parasites to antiparasitic treatment. However, there is scarce information about both threshold values and categorization of antiparasitic susceptibility for C. rogercresseyi. Bioassay technique allowed the determination of both mean values and the natural variation of EC50%, which were contrasted with available susceptibility thresholds. Results allowed to determine that parasites from the native fish host, Eleginops maclovinus, are susceptible to azamethiphos, deltamethrin and cypermethrin treatments, showing a high susceptibility profile to antiparasitics.

Modelling sea lice control by lumpfish on Atlantic salmon farms: interactions with mate limitation, temperature and treatment rules

Atlantic salmon farming is one of the largest aquaculture sectors in the world. A major impact on farm economics, fish welfare and, potentially, nearby wild salmonid populations, is the sea louse ectoparasite Lepeophtheirus salmonis. Sea louse infestations are most often controlled through application of chemicals, but in most farming regions, sea lice have evolved resistance to the small set of available chemicals. Therefore, alternative treatment methodologies are becoming more widely used. One increasingly common alternative treatment involves the co-culture of farmed salmon with cleaner fish, which prey on sea lice. However, despite their wide use, little is understood about the situations in which cleaner fish are most effective. For example, previous work suggests that a low parasite density results in sea lice finding it difficult to acquire mates, reducing fecundity and population growth. Other work suggests that environmental conditions such as temperature and external sea louse pressure have substantial impact on this mate limitation threshold and may even remove the effect entirely. We used an Agent-Based Model (ABM) to simulate cleaner fish on a salmon farm to explore interactions between sea louse mating behaviour, cleaner fish feeding rate, temperature and external sea louse pressure. We found that sea louse mating has a substantial effect on sea louse infestations under a variety of environmental conditions. Our results suggest that cleaner fish can control sea louse infestations most effectively by maintaining the population below critical density thresholds.

The Use of Kernel Density Estimation With a Bio-Physical Model Provides a Method to Quantify Connectivity Among Salmon Farms: Spatial Planning and Management With Epidemiological Relevance

Connectivity in an aquatic setting is determined by a combination of hydrodynamic circulation and the biology of the organisms driving linkages. These complex processes can be simulated in coupled biological-physical models. The physical model refers to an underlying circulation model defined by spatially-explicit nodes, often incorporating a particle-tracking model. The particles can then be given biological parameters or behaviors (such as maturity and/or survivability rates, diel vertical migrations, avoidance, or seeking behaviors). The output of the bio-physical models can then be used to quantify connectivity among the nodes emitting and/or receiving the particles. Here we propose a method that makes use of kernel density estimation (KDE) on the output of a particle-tracking model, to quantify the infection or infestation pressure (IP) that each node causes on the surrounding area. Because IP is the product of both exposure time and the concentration of infectious agent particles, using KDE (which also combine elements of time and space), more accurately captures IP. This method is especially useful for those interested in infectious agent networks, a situation where IP is a superior measure of connectivity than the probability of particles from each node reaching other nodes. Here we illustrate the method by modeling the connectivity of salmon farms via sea lice larvae in the Broughton Archipelago, British Columbia, Canada. Analysis revealed evidence of two sub-networks of farms connected via a single farm, and evidence that the highest IP from a given emitting farm was often tens of kilometers or more away from that farm. We also classified farms as net emitters, receivers, or balanced, based on their structural role within the network. By better understanding how these salmon farms are connected to each other via their sea lice larvae, we can effectively focus management efforts to minimize the spread of sea lice between farms, advise on future site locations and coordinated treatment efforts, and minimize any impact of farms on juvenile wild salmon. The method has wide applicability for any system where capturing infectious agent networks can provide useful guidance for management or preventative planning decisions.

Statistical modelling of sea lice count data from salmon farms in the Faroe Islands

Fiskaaling regularly counts the number of sea lice in the attached development stages (chalimus, mobiles and adult) for the salmon farms in the Faroe Islands. A statistical model of the data is developed. In the model, the sea-lice infection is represented by the chalimus (or mobile) lice developing into adult lice and is used to simulate past and current levels of adult lice-including treatments-as well as to predict the adult sea lice level 1-2 months into the future. Time series of the chalimus and adult lice show cross-correlations that shift in time and grow in size with temperature. This implies in situ the temperature-dependent development times of about 56 down to 42 days and the inverted development times (growth rates) of 0.018 up to 0.024 lice/day at 8-10°C. The temperature dependence DT=α1T+α2α3=17,840T+7.439-2.128is approximated byD1T=105.2-6.578T≈49 days at the mean temperature 8.5°C-similar to DchaT=100.6-6.507T≈45 days from EWOS data. The observed development times at four sites for a year (2010-11) were 49, 50, 51 and 52 days, respectively. Finally, we estimate the sea lice production from fish farms to discuss approaches to control the sea lice epidemics-preferably by natural means. This study is useful for understanding sea lice levels and treatments, and for in situ analysis of the sea-lice development times and growth rates.

The Phe362Tyr mutation conveying resistance to organophosphates occurs in high frequencies in salmon lice collected from wild salmon and trout

The parasitic salmon louse, and its resistance to chemical delousing agents, represents one of the largest challenges to the salmon aquaculture industry. We genotyped lice sampled from wild salmon and sea trout throughout Norway with the recently identified Phe362Tyr mutation that conveys resistance to organophosphates. These results were compared to data from lice sampled on farmed salmon in the same regions. The resistant (R) allele was observed in salmon lice from wild salmon and sea trout throughout Norway, although its frequency was highest in farming-intense regions. In most regions, the frequency of the R allele was higher in lice collected from wild sea trout than wild Atlantic salmon, and in all regions, the frequency of the R allele was similar in lice collected from wild sea trout and farmed Atlantic salmon. The R allele is only selected for in fish-farms where organophosphates are used for delousing. Therefore, our results suggest extensive exchange of lice between farmed and wild hosts, and indicate that in farming-dense regions in Norway, aquaculture represents a major driver of salmon louse population structure. Finally, these data suggest that the wild hosts within the regions studied will not delay the spread of resistance when organophosphates are used.

Managing aquatic parasites for reduced drug resistance: lessons from the land

Atlantic salmon farming is one of the largest aquaculture industries in the world. A major problem in salmon farms is the sea louse ectoparasite Lepeophtheirus salmonis, which can cause stress, secondary infection and sometimes mortality in the salmon host. Sea lice have substantial impacts on farm economics and potentially nearby wild salmonid populations. The most common method of controlling sea louse infestations is application of chemicals. However, most farming regions worldwide have observed resistance to the small set of treatment chemicals that are available. Despite this, there has been little investigation of treatment strategies for managing resistance in aquaculture. In this article, we compare four archetypical treatment strategies inspired by agriculture, where the topic has a rich history of study, and add a fifth strategy common in aquaculture. We use an agent-based model (ABM) to simulate these strategies and their varying applications of chemicals over time and space. We analyse the ABM output to compare how the strategies perform in controlling louse abundance, number of treatments required and levels of resistance in the sea louse population. Our results indicated that among the approaches considered applying chemicals in combination was the most effective over the long term.

A selection study on a laboratory-designed population of salmon lic (Lepeophtheirus salmonis) using organophosphate and pyrethroid pesticides

Resistance towards antiparasitic agents in the salmon louse (Lepeophtheirus salmonis) is a widespread problem along the Norwegian coast, reducing treatments efficacies and slowing down the envisioned expansion of Norwegian salmon production. The present study was conducted in order to assess the efficacies of two of the most widely used anti-parasitic substances–azamethiphos and deltamethrin–as well as assessing the benefit of having a resistant genotype compared to being fully sensitive when exposed to one of these substances. Atlantic salmon were exposed to a mix of salmon lice copepodids from a fully sensitive, a double resistant and a multi-resistant strain. Once the lice reached pre-adult stages, one group was exposed to 100 μg/L azamethiphos for 60 minutes, the other to 2 μg/L deltamethrin for 30 minutes, and the last was kept in a seawater control. Detached lice were collected at a series of time points following exposure, and all lice (immobilized and surviving) were analysed for both pyrethroid (sensitive “S” and resistant “R”) and azamethiphos (fully sensitive “SS”, heterozygous resistant “RS” and fully resistant “RR”) resistance markers. We found that the efficacies of deltamethrin on parasites with genotype S and R were 70.3 and 13.2%, respectively. The overall efficacy of the deltamethrin treatment was 32.3%. The efficacies of azamethiphos on parasites with genotype SS, RS and RR were 100, 80 and 19.1%, respectively. The overall efficacy of the azamethiphos treatment was 80.4%. Survival analyses revealed that the median survival time in deltamethrin-sensitive and–resistant parasites were 16.8 and >172 hours, respectively. The differences were even more pronounced in the azamethiphos-treated group, where SS, RS and RR parasites survived for 0.26, 6.6 and >172 hours, respectively. The substantial differences in survival between sensitive and resistant lice following treatment demonstrate the ability of medicinal treatments to drive genetic selection towards a much more resistant salmon lice population within a very short time span if there is no influx of sensitive genotypes.

Evosystem Services: Rapid Evolution and the Provision of Ecosystem Services

Evolution is recognized as the source of all organisms, and hence many ecosystem services. However, the role that contemporary evolution might play in maintaining and enhancing specific ecosystem services has largely been overlooked. Recent advances at the interface of ecology and evolution have demonstrated how contemporary evolution can shape ecological communities and ecosystem functions. We propose a definition and quantitative criteria to study how rapid evolution affects ecosystem services (here termed contemporary evosystem services) and present plausible scenarios where such services might exist. We advocate for the direct measurement of contemporary evosystem services to improve understanding of how changing environments will alter resource availability and human well-being, and highlight the potential utility of managing rapid evolution for future ecosystem services.

Lessons from sea louse and salmon epidemiology

Effective disease management can benefit from mathematical models that identify drivers of epidemiological change and guide decision-making. This is well illustrated in the host-parasite system of sea lice and salmon, which has been modelled extensively due to the economic costs associated with sea louse infections on salmon farms and the conservation concerns associated with sea louse infections on wild salmon. Consequently, a rich modelling literature devoted to sea louse and salmon epidemiology has been developed. We provide a synthesis of the mathematical and statistical models that have been used to study the epidemiology of sea lice and salmon. These studies span both conceptual and tactical models to quantify the effects of infections on host populations and communities, describe and predict patterns of transmission and dispersal, and guide evidence-based management of wild and farmed salmon. As aquaculture production continues to increase, advances made in modelling sea louse and salmon epidemiology should inform the sustainable management of marine resources.

Surveillance of the Sensitivity towards Antiparasitic Bath-Treatments in the Salmon Louse (Lepeophtheirus salmonis)

The evolution of drug resistant parasitic sea lice is of major concern to the salmon farming industry worldwide and challenges sustainable growth of this enterprise. To assess current status and development of L. salmonis sensitivity towards different pesticides used for parasite control in Norwegian salmon farming, a national surveillance programme was implemented in 2013. The programme aims to summarize data on the use of different pesticides applied to control L. salmonis and to test L. salmonis sensitivity to different pesticides in farms along the Norwegian coast. Here we analyse two years of test-data from biological assays designed to detect sensitivity-levels towards the pesticides azamethiphos and deltamethrin, both among the most common pesticides used in bath-treatments of farmed salmon in Norway in later years. The focus of the analysis is on how different variables predict the binomial outcome of the bioassay tests, being whether L. salmonis are immobilized/die or survive pesticide exposure. We found that local kernel densities of bath treatments, along with a spatial geographic index of test-farm locations, were significant predictors of the binomial outcome of the tests. Furthermore, the probability of L. salmonis being immobilized/dead after test-exposure was reduced by odds-ratios of 0.60 (95% CI: 0.42–0.86) for 2014 compared to 2013 and 0.39 (95% CI: 0.36–0.42) for low concentration compared to high concentration exposure. There were also significant but more marginal effects of parasite gender and developmental stage, and a relatively large random effect of test-farm. We conclude that the present data support an association between local intensities of bath treatments along the coast and the outcome of bioassay tests where salmon lice are exposed to azamethiphos or deltamethrin. Furthermore, there is a predictable structure of L. salmonis phenotypes along the coast in the data, characterized by high susceptibility to pesticides in the far north and far south, but low susceptibility in mid Norway. The study emphasizes the need to address local susceptibility to pesticides and the need for restrictive use of pesticides to preserve treatment efficacy.