Royal Society Publishing

Transmission dynamics of parasitic sea lice from farm to wild salmon

Martin Krkošek, Mark A Lewis, John P Volpe

Abstract

Marine salmon farming has been correlated with parasitic sea lice infestations and concurrent declines of wild salmonids. Here, we report a quantitative analysis of how a single salmon farm altered the natural transmission dynamics of sea lice to juvenile Pacific salmon. We studied infections of sea lice (Lepeophtheirus salmonis and Caligus clemensi ) on juvenile pink salmon (Oncorhynchus gorbuscha) and chum salmon (Oncorhynchus keta) as they passed an isolated salmon farm during their seaward migration down two long and narrow corridors. Our calculations suggest the infection pressure imposed by the farm was four orders of magnitude greater than ambient levels, resulting in a maximum infection pressure near the farm that was 73 times greater than ambient levels and exceeded ambient levels for 30 km along the two wild salmon migration corridors. The farm-produced cohort of lice parasitizing the wild juvenile hosts reached reproductive maturity and produced a second generation of lice that re-infected the juvenile salmon. This raises the infection pressure from the farm by an additional order of magnitude, with a composite infection pressure that exceeds ambient levels for 75 km of the two migration routes. Amplified sea lice infestations due to salmon farms are a potential limiting factor to wild salmonid conservation.

Keywords:

1. Introduction

Diseases can threaten wildlife when reservoir host populations are created, from which diseases can spill over into wildlife populations (McCallum & Dobson 1995; Daszak et al. 2000). Marine salmon farms located along wild salmon migration routes are spatially concentrated host populations that may serve as reservoirs and perturb the dynamics of any sympatric natural salmonid host–parasite system. Industrial-scale salmon farming has grown 55-fold over the last two decades (Porter 2003), and has spread globally, throughout and beyond the native ranges of wild salmonids. Sea lice (Lepeophtheirus salmonis Krøyer and Caligus spp.) are directly transmitted marine ectoparasites of salmonids: their life cycles have obligate parasitic and free-swimming stages, but no obligate intermediate host. Lice frequently infest farm salmon, and many studies have linked planktonic lice and lice parasitizing wild salmonids with the presence of farms (Tully & Whelan 1993; Costelloe et al. 1996, 1998a,b; Todd et al. 1997; Mackenzie et al. 1998; Tully et al. 1999; Bjørn et al. 2001; Bjørn & Finstad 2002; Marshall 2003; Morton & Williams 2004; Morton et al. 2004; McKibben & Hay 2004; Penston et al. 2004; Carr & Whoriskey 2004). Sympatric wild salmonid populations may then be affected: farms have been implicated in the infestation and collapse of pink salmon (Oncorhynchus gorbuscha) cohorts in Pacific Canada (PFRCC 2002; Morton & Williams 2004; Morton et al. 2004), wild sea trout (Salmo trutta) and Atlantic salmon (Salmo salar) populations in Europe (McVicar 1997, 2004). Assessing the role of salmon farms in these declines first requires a comprehensive understanding of how farms alter the natural transmission dynamics in this system, particularly involving the vertical transmission of lice from adult salmonid cohorts to susceptible juveniles. Here, we use mechanistic spatial models coupled with extensive field data to analyze how an isolated salmon farm perturbed the natural transmission dynamics of sea lice to juvenile pink and chum (Oncorhynchus keta) salmon along two narrow and restricted migration corridors in British Columbia, Canada.

Two native sea louse species coexist on salmonids in Pacific waters off North America: L. salmonis and C. clemensi (Parker & Margolis 1964). Both species have planktonic larval stages, and juvenile and adult parasitic stages. Planktonic nauplii hatch from gravid parasitic females and develop into infective copepodids. After settling on a host fish, copepodids develop through distinct chalimus and motile pre-adult and adult stages (Kabata 1972; Johnson & Albright 1991a,b). Attached stages feed on the mucus, scales and blood of the host fish leading to osmotic stress and emaciation of sufficiently infected hosts (Pike & Wadsworth 2000). The ecologies of these species differ in that L. salmonis are salmonid specialists, whereas C. clemensi are generalists occurring on members of several piscine families (Parker & Margolis 1964).

Pink and chum salmon are unique among salmonids in their precocious entry into marine waters: juveniles emerge from gravel and immediately leave their natal streams at approximately 28–35 and 30–40 mm fork length, respectively (Groot & Margolis 1991; personal observation). This makes them the smallest salmonids to contend with marine parasites. In this early marine phase, these species form mixed schools and occupy the same near-shore (often intertidal) habitats during their seaward migration (Groot & Margolis 1991). Marine entry begins in March (Groot & Margolis 1991) with juveniles initially free of lice. Direct vertical transmission of sea lice from returning conspecifics does not begin until late July, but transmission can occur indirectly immediately upon marine entry of juveniles via alternate wild hosts and farm salmon. The magnitudes of these transmission routes are unknown, but if the primary natural route is direct then fish farms may provide a significant transmission pathway for lice that operates much earlier in the salmon life cycle then would otherwise naturally occur.

2. Field methods

We sampled juvenile pink and chum salmon at 1–4 km intervals for 40 and 60 km along two narrow and restricted migration routes relative to an isolated salmon farm, and quantified the abundance of copepodid, chalimus and motile stages. The migration routes are labelled route I and II, respectively, and were located in marine fjordic habitats in British Columbia, Canada. An isolated salmon farm (farm A) was situated midway along both migration routes. A second salmon farm (farm B) was situated such that migratory salmon indirectly passed within 7 km of it toward the end of route I. We did not sample for approximately 20–60 km of the migration routes between the landward end of the study area and the various natal streams of the studied populations. Two replicate sets of samples (79–237 juvenile pink and chum salmon per sample) were obtained from each site in the spring of 2003 (17–27 April and 9–23 May). Four datasets result: two replicates of the two migration routes, each representing a spatially structured snapshot of louse population structure. Datasets are labelled: I-April, I-May, II-April, II-May. Route I and route II datasets share the same data landward of farm A. Further details of the field site are withheld to maintain industry anonymity.

At each site, juvenile pink and chum salmon (measuring 2.8–10 cm fork length) were captured by beach seine (30 m long, 4 mm mesh size). The beach seine was drawn in to approximately 1 m2×30 cm and a live subset was retained in 30 l buckets of seawater using a 15×15 cm2 dipnet. Care was taken to maintain randomness by varying the location, depth and speed of the subsampling procedure. Individual fish were removed from a bucket with a 15×15 cm2 dipnet and placed in a 15×27 cm2 clear plastic envelope for analysis. Fins and lateral, dorsal and ventral surfaces were scanned with a 10× magnification hand lens and the number of copepodid, chalimus and motile lice was recorded. Ovigerous females were classified as either L. salmonis or C. clemensi. Fish were then returned to buckets where more than 99% recovered and were subsequently released at the location of capture. This nonlethal assay is known to produce slight underestimates of copepodid and chalimus abundances and is only feasible for fish less than 10 cm fork length (Krkošek et al. in press). Temperature readings were taken at most sites and salinity readings were taken at a subset of sample sites around farm A with a Hydrolab Quanta electronic water quality meter.

3. Model

The large-scale movement of louse larvae in long and narrow fjordic habitats is limited to movements up and down the habitat length, and we model it with a one-dimensional domain. Juvenile salmon migrate down this domain, initially free of lice, and first encounter infective copepodids that originate from two primary host populations: farm salmon and sympatric wild hosts. As the infection progresses, juvenile salmon become a secondary source of louse larvae themselves. A distinction between the primary sources is their spatial distributions: a farm is a point source of lice whereas wild hosts are a distributed source. Each source corresponds to distinct spatial profiles in the dynamics of free-swimming and parasitic stages, which form the basis for distinguishing between farm-origin and natural-origin lice in spatially structured data.

Three models were developed (see Appendix A): Ψ0 includes only natural primary sources of lice; Ψ1 includes only farm sources of lice; and Ψ2 includes both farm and natural primary sources of lice. Larvae from wild hosts are described by a uniform spatial distribution; a diversity of free-swimming hosts combined with the planktonic stages of lice suggests these larvae are well mixed. Advection–diffusion submodels describe the movements of nauplii and copepodids from the farm and the juvenile salmon. Infection dynamics on migratory juvenile salmon are linked to these larval distributions by a spatially dependent Poisson process (Papoulis 1963) that assigns a probability to each datum. In contrast to the common aggregated distributions of macroparasites on their host populations (Shaw & Dobson 1995; Shaw et al. 1998), the sea lice data conformed well to the characteristics of a Poisson process (see figure 2). The mean abundances of copepodid, chalimus and motile stages areEmbedded Imagewhere v is the mean seaward migration velocity of salmon, β is the transmission coefficient and sc and sh are the proportions of lice surviving copepodid and chalimus stages. The λs are the distances salmon travel in the cumulative mean durations of copepodid (c), chalimus (h) and motile (m) stages. The spatial distribution of infective larvae is L(x): the number of planktonic copepodids within the unit of volume defined by the detection radius of lice centred at x. Only relative values become important, so the detection radius is not explicitly required. Maximum likelihood was used to fit models, likelihood ratios (Hilborn & Mangel 1997) to test if farms do not infect wild juvenile salmon and if secondary infection from the juvenile salmon is significant, and Akaike information criteria (AIC; Burnham & Anderson 2002) to select the best model among Ψ0, Ψ1 and Ψ2. Both species of lice and juvenile salmon were included in the same analysis due to similarities in host behaviour and parasite life cycles.

4. Results

A total of 41 samples were collected across the four datasets yielding a total of 5514 juvenile salmon that were sampled for sea lice infections. From these fish, we counted a total of 552 copepodids, 2078 chalimi and 1015 motiles. Of these motiles, there were 12 ovigerous L. salmonis females and 53 ovigerous C. clemensi females. The distributions of lice on juvenile salmon ranged from slightly underdispersed to slightly overdispersed, and all samples closely followed the 1 : 1 variance to mean line, which accounted for 96% of the variability in those data (figure 1). Sea temperatures were 8–10 °C and salinity readings were in the range 27.6–30.3‰.

Figure 1

Log variance versus log mean for all 41 samples of copepodids (black circles) chalimi (grey squares) and motiles (clear triangles). There are 79–237 fish per sample. The solid line is the variance=mean line which accounts for 96% of the variation. The dashed line is the best-fit linear model, which has slope 1.08 and accounts for 97% of the variation. Compare with fig. 5 in Shaw & Dobson (1995).

The analysis revealed that lice from farm salmon infected wild juvenile salmon (likelihood ratio test, R>242.4, d.f.=8; p<0.001 for all datasets) and that juvenile salmon were a secondary source of infection (likelihood ratio test, R>19.8, d.f.=1, p<0.001 for all datasets). Model selection statistics (AIC) indicated Ψ2 was the superior model for all datasets; the minimum ΔAIC between all models and Ψ2 was 22.4 and the probability Ψ2 was the best model approached unity in each dataset (minimum Akaike weight =0.9998). Statistical details are provided in the Electronic Appendix. These results strongly suggest that both wild and farm primary hosts were important sources of sea lice infecting juvenile salmon and in addition that juvenile salmon themselves were an important secondary source of sea lice.

The sea lice data were spatially structured: most lice were observed on juvenile salmon after they passed farm A (figure 2). Juvenile salmon carried low burdens of lice prior to their encounter with farm A. Near farm A, a large abundance of parasitic copepodids was observed followed by subsequent peaks in copepodids, chalimi and motiles further down the migration routes. The fit of model Ψ2 agrees well with these data (figure 2) and explains these patterns. Prior to passing farm A, louse abundances were at natural ambient levels determined by a balance between immigration and emigration/death rates through each stage. Near farm A, a large cohort of lice colonized the juvenile salmon. These lice developed through subsequent chalimus and motile stages as their hosts migrated, producing the spatially displaced peaks in chalimi and motiles. Larvae were subsequently produced from the motile population on the juvenile hosts and produced the secondary infection waves of copepodids and chalimi apparent in the data.

Figure 2

Abundance of parasitic louse stages on juvenile pink and chum salmon at points along their migration routes relative to a salmon farm located at x=0 (farm A). Salmon migrate in the rightward (seaward) direction. Columns correspond to datasets and rows correspond to louse stages. Error bars are bootstrapped 95% confidence intervals. Solid lines are the maximum likelihood best fits of model Ψ2 to each dataset. Sample sizes are in the range 79–237 fish per sample.

Due to the mechanistic structure of the model, it was possible to analyse the transmission dynamics across each component in this system using parameter estimates from the maximum likelihood fits of model Ψ2 to each dataset (figure 3). The analysis reveals that larvae originating from the farm and from the farm-produced cohort of lice on the juvenile salmon were responsible for the majority of the infection dynamics observed in the data. These differences can be quantified (table 1). Assuming the farm is 0.2 km in length, then the production of infective copepodids by the farm was on average 3.14×104 times greater than natural production in this spatial interval. This corresponds to an infection pressure near the farm that was on average 73 times greater than ambient levels and exceeded ambient levels for an average of 30 km. Inclusion of the dynamics of the farm origin cohort of lice on the juvenile salmon suggests that the production of larvae due to the farm was 2.08×105 times greater than ambient levels per unit space, with a composite spatial profile of infection pressure that exceeded ambient levels for 75 km (figure 3).

Figure 3

The spatial distributions of planktonic copepodids inferred by model Ψ2 on a relative scale. Juvenile salmon migrate in the rightward (seaward) direction. The thick grey line is the total abundance of copepodids produced by all sources. The horizontal lines near zero are the ambient infection pressures. The thin dark curves oriented about x=0 are the distributions of copepodids produced directly by farm A and the second curves to the right are the distributions of copepodids produced by the farm-origin cohort of lice on the juvenile salmon. The latter distribution was found by solving the model with κ=0 to eliminate any contribution of lice from natural sources. Corresponding datasets are I-April (a), I-May (b), II-April (c), and II-May (d).

View this table:
Table 1

Summary of the spatial infection pressures caused by the salmon farm relative to ambient levels. (Larval production per unit space was calculated based on a farm length of 0.2 km. Secondary larval production from the farm-origin cohort of lice on the juvenile salmon was calculated by solving the model with κ=0 to eliminate any contribution from natural origin lice. Infection pressure is the local abundance of planktonic copepoidids.)

5. Discussion

Our analysis identifies a signal amidst noise in empirical data. The signal is the ensemble of spatial distributions of lice predicted by the models. Any source of error that confounds the model predictions would detract from the results, while any unaccounted sources of variation that add noise to the data reduce the statistical power of the analysis. There were many unaccounted sources of noise in the sea lice data: variation in temperature and salinity; inter- and intra-specific variations in lice life-history parameters; inter- and intra-specific variations in juvenile salmon behaviour and host–parasite interactions; mixing of juvenile salmon with different immunological histories; deviations from a mean seaward migration velocity; deviations from a uniform juvenile salmon spatial distribution; temporal variations in the infestation levels on the farm; infection originating from farm B; density-dependent effects on louse survivorship and/or host mortality, potential patchiness in planktonic louse distributions; and generally, the 1D mathematical representation of a dynamic 3D biological system. The only confounding source of error would occur if a population of natural hosts was aggregated around the farm, producing the spatial distributions we attributed to the salmon farm. Were this to occur, our calculations indicate such an aggregated wild population would have to be either four orders of magnitude more dense than anywhere else in the study area or four orders of magnitude more infested than other wild hosts. Furthermore, the spatial distributions in the data require a point source that is stationary for at least two louse life cycles. For L. salmonis, this is at least 100 days at 10 °C (Johnson & Albright 1991a,b). Therefore, not only would such a population be unrealistically dense or infested, it would also be unrealistically stationary.

Given the paucity of confounding factors and abundant sources of noise, the strength of our results suggests the unaccounted sources of variation must be trivial relative to the effect of the salmon farm. Indeed, our calculations suggest the farm raised infection levels by four orders of magnitude, and it is unlikely that other sources of error could vary by such a magnitude. However, a degree of noise is reflected in the variability of some parameter estimates (see Electronic Appendix) and in the deviation of some data from the best-fit models in figure 3. The statistical results are strong because those deviations are small relative to the overall dynamic pattern observed in the data and that pattern could only occur were a salmon farm the primary driver of sea lice dynamics on wild juvenile salmon. This is in agreement with Morton et al. (2004), who found virtually no lice on juvenile pink and chum salmon in several regions of British Columbia without salmon farms. However, a natural infestation of juvenile pink salmon with C. clemensi has been reported (Parker & Margolis 1964). Other European studies have found external oceanographic co-variates more important in determining louse dynamics (e.g. Marshall 2003), but such studies taken over seasonal time-scales import large temporal variations in temperature and salinity that affect the dynamics of all lice, but tell little about the interactions across wild and farm host populations. We have avoided such confounds by focusing on smaller temporal scales and explicitly examining the interactions among all host populations.

The models predict a close proximity of louse larvae to their source (figure 4a) and this agrees well with the abundance of copepodids around the salmon farm and around the peaks in the motile stages in the data. These results also agree with other investigations of spatial distributions of louse larvae that found consistent high levels near farms (Costelloe et al. 1996, 1998a,b). In contrast, one may expect most larvae to be transported long distances from their source owing to the one to two week lifespan of planktonic stages. This may be correct for a temporal pulse of larvae (figure 4b), but a farm is a continuous source, and it is the continuous production of larvae that produces the sustained maximum density of larvae near a farm. However, while a farm is a continuous source, population dynamics of lice on farms tend to oscillate due to population growth between intermittent treatment events (Revie et al. 2002) and this may affect the distribution of larvae around the source. Using biologically realistic parameter values, we found that the spatial profile of larvae produced by an oscillating source with amplitude half the mean tends to the predicted distribution of the constant mean source (figure 4c and d).

Figure 4

The spatial distributions of sea lice larvae around a point source. (a) The expected differences in the spatial distributions between nauplii (dashed line) and copepodids (solid line). (b) Differences in the spatial distributions of nauplii produced by a pulse release (dashed lines at successive time intervals) and by a continuous constant source (solid line) under strong advection. The dynamics of an oscillating source of amplitude half the mean (c). The frequency of oscillations is 14.8 per year, calculated from Revie et al. (2002). Three solutions for the spatial distributions of nauplii about the oscillating source are plotted in (d) for the time points indicated by the vertical dashed lines in (c). The expected distribution for the constant mean source is shown by the filled circles.

While no general conclusions can be made on the transmission dynamics of lice from farm to wild salmon based on this study alone (we only considered a single isolated salmon farm), the summation of these results with those of the current literature (listed in §1) strongly suggest these processes pervade most systems where wild and farm salmonids occur in sympatry. The data and analysis presented here extend our current understanding of lice interactions between wild and farm salmon by providing quantitative estimates of transmission rates across host populations while also considering subsequent transmission dynamics of lice within the wild population; initial transmission from farm to wild juvenile salmon may be minor compared with the subsequent population dynamics of the farm-produced lineage of lice. The spatial profiles of copepodid abundances inferred by the models are strikingly consistent across all four datasets, lending support to this result.

Louse transmission from farm to wild salmon occurs during the most vulnerable salmon life‐history phase (because of their small size and recent marine transition). For pink and chum salmon, this is a period prior to the return of adult conspecifics when near-shore native adult salmon abundances are lowest within annual variations (Groot & Margolis 1991). Presumably then, the natural rates of transmission to these juvenile fish are also at their lowest during this time. This agrees with the low rates of natural transmission shown here and by Morton et al. (2004) and suggests the primary natural route of vertical transmission to juveniles is direct: from returning conspecifics months later in the salmon life cycle after significant body growth in the juveniles has occurred. The effect of a single salmon farm was to raise infection pressure by four to five orders of magnitude over natural rates, exceeding ambient levels for 75 km along two wild salmon migration routes. In some areas, juvenile salmon must migrate past several salmon farms, compounding these effects. This may have severe consequences, as these species may not have evolved mechanisms to contend with such high infestation pressures so early in their life‐history. In fact, the near absence of natural vertical transmission of lice during this time would facilitate the evolution of the precocious marine phase evidenced by these species. Currently, the lethal infection levels for juvenile pink and chum salmon are unknown, but given their small size on marine entry (0.2–0.5 g; Groot & Margolis 1991, personal observation) and the known lethal limits for other species (1 motile louse per 0.75–1.6 g body weight for European salmonids; Grimnes & Jakobsen 1996; Bjørn & Finstad 1997), any amplification of louse transmission during this life-history phase may threaten the viability of the affected populations. The European experience is that louse management on farms has not resulted in positive effects on depressed wild salmonid populations (McVicar 2004). The rate of industry expansion in British Columbia is much faster than evolutionary time-scales, and therefore, an adaptive response by native salmonids is unlikely. Declines of wild salmonid populations correlated with sea lice infestations and salmon farms in Canada (PFRCC 2002; Morton & Williams 2004; Morton et al. 2004) and Europe (McVicar 1997, 2004) demand a cautionary approach to continued industrial expansion.

Acknowledgments

We are indebted to Alexandra Morton and Eric Nelson for their generous hospitality and logistical support of the fieldwork. We thank Marjorie Wonham for comments on the manuscript and Chris Jerde for discussions on the analysis. Special thanks to Chris Williamson, Brendan Connors, Michelle Nelson and M.W. for their hard work in the field. This work was supported by The David Suzuki Foundation; Raincoast Conservation Society; Raincoast Research Society; an NSERC Discovery Grant to J.P.V.; NSERC and a Canada Research Chair to M.A.L.; and a Walter H. Johns Graduate Fellowship and NSERC Industrial Postgraduate Scholarship to M.K.

Footnotes

  • Present address: School of Environmental Studies, University of Victoria, Victoria, BC, Canada V8W 2Y2.

  • The supplementary Electronic Appendix is available at http://dx.doi.org/10.1098/rspb.2004.3027 or via http://www.journals.royalsoc.ac.uk.

  • As this paper exceeds the maximum length normally permitted, the authors have agreed to contribute to production costs.

    • Received September 19, 2004.
    • Accepted November 24, 2004.

References

View Abstract