Skip to main content

Survival, healing, and swim performance of juvenile migratory sea lamprey (Petromyzon marinus) implanted with a new acoustic microtransmitter designed for small eel-like fishes

Abstract

Background

Little is known about the transformer stage of the parasitic lampreys, a brief but critical period that encompasses juvenile out-migration from rivers to lakes or oceans to begin parasitic feeding. Information about this life stage could have significant conservation implications for both imperiled and invasive lampreys. We investigated tag retention, survival, wound healing, and swim performance of newly transformed sea lamprey (Petromyzon marinus) implanted with a new micro-acoustic transmitter, the eel–lamprey acoustic transmitter (ELAT), in a controlled laboratory environment.

Results

The 61-day survival of our tagged subjects was 71%, within the range reported in similar studies of juvenile lampreys. However, survival was significantly lower in the tagged animals (vs control), with no effect statistically attributable to measures of animal length, mass, condition, or population of origin (Great Lakes vs. Atlantic drainage). Mortality in tagged fish was concentrated in the first four days post-surgery, suggesting injury from the surgical process. An unusually long recovery time from anesthesia may have contributed to the increased mortality. In a simple burst swim assay, tagged animals swam significantly slower (− 22.5%) than untagged animals, but were not significantly different in endurance swim tests. A composite wound healing score at day four was a significant predictor of maximum burst swim speed at day 20, and wound condition was related to animal mass, but not length, at the time of tagging.

Conclusions

Impairments to survival and swim performance of juvenile sea lamprey implanted with the ELAT transmitter were within currently reported ranges for telemetry studies with small, difficult to observe fishes. Our results could be improved with more refined anesthesia and surgical techniques. The ability to track migratory movements of imperiled and pest populations of parasitic lampreys will improve our ability to estimate vital rates that underlie recruitment to the adult population (growth, survival) and to investigate the environmental factors that regulate the timing and rates of movement, in wild populations.

Background

Migration is a time of heightened threat and uncertainty regarding the risk of predation and the localization of critical resources [60]. This is particularly so when migration occurs early in life, when variation in daily growth and survival may be substantial and influence recruitment to the adult population [24, 37]. A number of economically and ecologically valued fishes exhibit juvenile feeding migrations from rivers into lakes or oceans, including salmon (Oncorhynchus spp.) and lake sturgeon (Acipenser fulvescens). Threats to the viability of these populations, and to control invasive species with similar life histories, have prompted considerable interest in ascertaining the behavior, timing, and survival of juvenile fishes as they move downstream [35, 56].

Several lamprey species are among the most poorly understood fishes that exhibit juvenile migration, including species that range from imperiled and ecologically or commercially valued (e.g., Pacific lamprey Entosphenus tridentatus, European river lamprey, Lampetra fluviatilis, sea lamprey Petromyzon marinus in select European regions) to pestilential and targeted for reduction (e.g., invasive sea lamprey in the Laurentian Great Lakes). These lampreys spawn in rivers, after which newly hatched larvae (ammocoetes) bury in stream sediments for typically four to seven years feeding on organic detritus and algae [14]. At the end of the larval stage, they transform and migrate downstream into estuaries or large lakes to commence parasitic feeding on fishes [1, 14, 21]. This period of physiological and geographical transition (aka the transformer stage) [9] is of short duration, yet is considered critical to the development of effective conservation and management practices [21]. For example, out-migrating lamprey within their native ranges typically must transit through a gauntlet of dams and their affiliated water intake structures, resulting in physical injury, increased predation, and direct or delayed mortality [27, 48]. Providing for safe passage through these structures is necessary for the protection of the species [26].

In the Laurentian Great Lakes, where a single sea lamprey may consume upwards of 21 kg of fish biomass, control is achieved through the application of selective pesticides (lampricides) to kill larvae in streams prior to transformation [5, 43]. Rivers are selected for lampricide treatment based on estimates of larval abundance, informed by population models and expert opinion [34]. The ability to estimate the fate of out-migrating sea lampreys would allow treatment decisions to be made based on probable parasite production vs. larval abundance, incorporating system-specific differences in rates of survival through migration and first feeding. For example, in Lake Michigan, the lower reaches of rivers systems often contain drowned river mouth lake–wetland complexes that are rich in piscivorous predators [38], and piscivores are important predators of out-migrating juvenile lampreys [48, 58, 62]. Life stage-specific spatial population models for sea lamprey currently assume all out-migrating sea lamprey have an equal likelihood of survival until arrival in the lakes, regardless of differences in migratory distance, river size, habitat types, or predator populations, each of which contribute to differences in mortality and growth among watersheds [31, 59].

Survival estimation during migration requires tracking an animal’s status over considerable distance and/or time. The most frequently used tool to monitor the movements of fishes and other aquatic organisms over large distances is telemetry, involving the implantation or attachment of a transmitter that may be detected when the fish moves into the range of a receiver [32]. Until recently, small (≤ 12 mm) passive integrated transponder (PIT) tags were the only transmitters available for the study of out-migrating lamprey. PIT tags lack an internal power source, requiring tagged fishes to move within 1–2 m of the antenna to be detected [3, 53]. The proportion of PIT-tagged larval and out-migrating lamprey detected by PIT antennas is relatively low, with reported detection rates ranging from 5 to 14% [13, 47]. Internally powered telemetry tags (acoustic and radio transmitters) offer an alternative that substantially increases the detection range (tens to hundreds of meters) and detection probability (as high as 80–100% on a well-designed array covering the entirety of the stream channel) at a single receiver [45]. Acoustic transmitters have several desirable properties for use in small fishes, in particular, they do not require an external trailing antenna (as radio telemetry typically does), that may result in negative impacts to swim performance, lower likelihood of predator attack, and higher survival in tagged individuals (reviewed by Crossin et al. [11] and references therein).

A primary assumption of telemetry analysis is the fish’s movement behavior (e.g., swim speed, timing) is not substantially altered by bearing the tag, or because of the implantation process. Historically, acoustic transmitters have proven too large to implant into small fishes with narrow body cavities. Transmitters with sufficient battery life to support demographic studies of fishes require large batteries that would impede movement and survival in small fishes (Liedtke 2019). This situation is compounded in anguilliform swimming fishes, where tags may physically impede propagation of the propulsive wave or compress organs along the length of the very narrow body. The fabrication of a novel microbattery and transmitter suitable for use in juvenile sea lamprey—the Eel–Lamprey Acoustic Tag or ELAT—offers a potential solution to this problem [15, 39]. This tag is compatible with the JSATS (Juvenile Salmon Acoustic Telemetry System) receiver, with a detection range of 80–140 m, a distance 100 × greater than the approximate range of a 12 mm PIT tag [17]. The transmitter has a source level of a 147 dB and is in a cylindrical encasing with dimensions 12 mm × 2 mm, 0.08 g dry-mass and transmits at a frequency of 416.7 kHz (± 0.5%). Currently, the standard PRI (Pulse Rate Interval) of five seconds allots approximately 30 days of use [16]. Pilot field studies using the ELAT tag have demonstrated the high detection rates (> 95%) in juveniles of both Pacific lamprey and American eel (Anguilla rostrata) [17, 39]. At approximately the same size as the PIT tags used in previous studies of juvenile lampreys (12 mm x 2 mm, 0.08 g), this is the first acoustic telemetry tag that does not persistently violate the two percent of body mass standard for use without impediment of movement in juvenile anguilliform fishes due to its cylindrical shape [46, 53, 63]. Here, we report an assessment of tag retention, survival, wound healing, and impacts to swim performance of the ELAT transmitter surgically implanted into newly transformed sea lampreys. Our goal was to perform a comprehensive evaluation of the impacts of ELAT transmitter implantation to support the establishment of criteria for use in field studies designed to enumerate critical demographic parameters (e.g., stage-specific mortality rates) and the movement ecology of the poorly understood transformer stage.

Methods:

Study design

The objective of this study was to examine the effects of the ELAT tag on sea lamprey survival, condition, and swim performance. To achieve this, we surgically implanted ELAT tags into 59 newly transformed juvenile parasitic sea lamprey and examined tag retention, survival, wound healing, and swim performance during a 61-day study period, corresponding to anticipated tag life at a PRI of 10 s. Wound healing assessment involved two scoring metrics, wound closure and wound inflammation, and a composite score that combined the two measures into an overall measure of wound condition. Swim performance was analyzed by comparing two performance metrics, maximum burst swim speed (20 days post-surgery) and time-to-exhaustion swimming (29–31 days post-surgery), to untagged control fish. We focused on the animal’s burst swim ability, as it is a facultative behavior suggestive of ability to escape a perceived threat, and swim-to-exhaustion, as it is characteristic of long stretches of active swimming during a migration [19]. The swim measures were separated by ten days to ensure the stress of the burst swim test did not influence the results of the time-to-exhaustion test. All measures were statistically compared to a group of 54 control animals that did not experience the anesthesia, surgery, or tag implantation. Due to the difficulty in acquiring this life history stage in sea lamprey, subjects were collected from both Great Lakes and Atlantic basin drainages. Where appropriate, source was included as a covariate in the statistical models.

Collection and housing of subjects

Sea lamprey used in this study were collected from tributaries in the Great Lakes Basin (N = 64), and a hydropower diversion canal on the Connecticut River that flows into the Atlantic Ocean (N = 49).

We selected 113 sea lamprey ranging 140 mm–160 mm total length (TL) based on recommendations from previous studies using similarly sized PIT tags in juvenile sea lamprey [63] and Pacific lamprey [46]. We chose to include animals across a range of sizes near the median of the size distributions in the Great Lakes and the Atlantic basin (Fig. 1) to examine the relationship of body size to wound healing and swim performance.

Fig. 1
figure 1

Density plots of total length (mm) for historic collections of out-migrating juvenile sea lamprey for animals captured from Great Lakes streams (Gray, N = 267) and Atlantic drainages in the northern United States (red, N = 653) are overlaid with a stacked histogram representing the two harvest locales and illustrating the experimental subjects total length (Great Lakes: N = 59, Anadromous: N = 54). Data for historical collections were compiled by J. Hume, Michigan State University

Subjects were housed at the U.S. Geological Survey Hammond Bay Biological Station from February 11, 2020 to April 13, 2020. Lampreys were separated by source and kept in eight separate 23L 40 cm x 22 cm x 26 cm tanks that were supplied with constant flow-through water from Lake Huron and continuous aeration. As our findings are intended to inform use of the ELAT tag in natural rivers, the holding tanks were constructed to include certain natural features the animal would experience during out-migration. First, previous tagging effects studies observed fungal growth and infection on experimental subjects held in the laboratory [7, 49, 50]. To ensure the microbiome the sea lamprey were exposed to was akin to the natural flora, we lined the tanks with a 4 cm layer of mixed substrate (cobble, sand, fine-grained sediment) collected from nearby Schmidt’s Creek. In addition, unfiltered inflowing water from Lake Huron was heated to a temperature typical of central Lake Michigan tributaries during the fall migration (~ 4.5° C).

Surgical implantation of the transmitter:

ELAT transmitters were fabricated at the Pacific Northwest National Laboratory by Daniel Deng. This experiment used a modified ELAT where the microbattery was replaced with a PIT transmitter to allow for individual identification of the subjects via a hand scanner. The measured dimensions of the modified tags used in this study were as follows (mean ± 1 SD): weight in air, 0.0809 ± 0.003 gm; length, 12.07 ± 0.10 mm; and width, 0.20 ± 0.003 mm. On February 12, 2020, 59 fully transformed sea lamprey (per [76]) were surgically implanted with modified ELAT transmitters (30 Great Lakes and 29 Atlantic drainage), and 54 entered the experiment as control animals (34 Great Lakes and 20 Atlantic drainage). Animals in the tagged group measured (mean ± 1 SD) 150.16 ± 5.69 (total length, mm) and weighed 4.37 ± 0.64 g. Animals in the control group measured (mean ± 1 SD) 148.30 ± 6.08 (total length, mm) and weighed 4.18 ± 0.50 g.

The implantation surgery followed protocols established by Moser et al. [49] and Christiansen [7]. An anesthetic bath of AQUI-S-20E (10% eugenol) and lake water was prepared to sedate the animals undergoing the tagging procedure. Anesthesia concentrations from prior studies (0.02–0.06 ml l−1, USFWS 2013, Simard 2017) proved ineffective at inducing stage-IV anesthesia in pre-surgery screening [66]. Similar occurrences have been noted for other anesthetic agents applied to juvenile lamprey when using doses recommended for teleost fishes. Stage-IV anesthesia was induced in (mean ± 1 SD) 17.0 ± 4.5 min via immersion in a solution of 0.7 ml l−1 AQUI-S-20E (eugenol concentration = 70 ppm). When each animal reached stage-IV anesthesia, it was placed on a damp surgery board and the surgeon (T.F. Haas) made a lateral 3 mm incision into the body cavity adjacent to the 13th myomere using a Premier Edge Restricted Depth 3 mm microscalpel (OASIS® Medical). A disinfected (immersion in 70% ethanol for ≥ 2 min) then rinsed transmitter was inserted posteriorly into the opening until it was completely enveloped in the animal’s body cavity (mean ± 1 SD = 1.4 ± 0.93 min, range = 0–4 min); no suturing or other method of artificially closing was used [13, 14, 49]. After the transmitter was fully inserted, each animal was placed into an oxygenated recovery tank until active swimming or suction attachment to the side of the tank was observed, indicating initial recovery (mean ± 1 SD time to recovery = 102.08 ± 67.13 min, range = 18–414 min). After recovery, the subjects were returned to their holding tanks for the experimental phase.

Mortality, wound healing, and condition:

Mortalities for both tagged and control treatment groups were enumerated daily. Wounds were visually evaluated by a single individual (TFH) at days 4, 20, 29, 44, and 61 for all animals alive and retaining a transmitter at the time of observation. For evaluation, animals were removed from their tank and individually placed in a small, clear, tank where a photo was taken when the animal ceased movement. No anesthesia was used. Control animals were not examined assessment criteria were modified from that of Wagner et al. [72] and Moser et al. [49], where wound closure and wound inflammation were scored separately on a scale of 1–4 (Table 1). An additional metric, a composite wound score, was produced from the wound closure and wound inflammation scores. This was formulated by plotting the wound closure and wound inflammation scores on an X–Y plane (axes range from 0 to 4) and measuring the Euclidean distance from the origin to the observed scores (Eq. 1). To ease interpretation, that distance was then rescaled to a range of 0–10, based on the minimum (1.41, i.e., wound closure and inflammation scores of 1) and maximum (5.657, i.e., wound closure and inflammation scores of 4) possible distances from the origin:

$$\left( {\frac{Euclidean \,Distance - \sqrt 2 }{{\left( {\surd 32 - \surd 2} \right)/10}}} \right).$$
(1)
Table 1 Criteria for assessing wound closure and wound inflammation. Modified from Moser et al. [49] and Wagner [72]

Swim performance

Twenty days after implantation and after the second wound healing assessment, we measured burst swim velocity in all implanted and control animals remaining in the study. Animals were placed in a 155 × 13x10cm plastic trough filled with 8 cm of water at ± 1 °C of the holding tank temperature. A plastic mesh grid marked in 5 cm intervals was positioned in the bottom of the trough to measure distance moved, and all interior surfaces were covered with plastic mesh to prevent animals from attaching to the trough surface. Each animal was placed in the lower end of the trough and allowed a five-minute acclimation period prior to testing. After the acclimation period, maximum speed was measured by inducing a fast-start (“startle”) response. This was done by squirting approximately 3 ml of water from a 5 ml pipette at the water surface above the animal’s head [12, 50]. Each animal underwent three burst trials with a 3-min recovery period in between trials. A Go-Pro camera (60 frames/sec.) was mounted above the trough to record trials, and video footage was processed and analyzed using Kinovea® motion analysis software. Burst speed was measured as the distance the animal traveled in the initial 30 frames (0.5 s) immediately following the first frame showing ripples produced by the pipette ejection. The maximum observed speed (cm s-1) of the three trials was used in analyses below, hereafter referred to as an individual’s burst speed.

On days 30 through 32 post-surgery, all implanted animals remaining in the study (n = 30) and 32 control animals chosen haphazardly from approximate same source proportions were subjected to an endurance swimming test (time-to-exhaustion) in a swim tunnel. Animals were individually placed in a 48.2 L plastic mesh-lined chamber of a 121 L Blazka-type swim respirator, with a Leeson Washguard Adjustable Speed AC Motor and Controller. The impeller was powered to 4.5 Hz. This speed was chosen as a result of preliminary trials where lamprey outside of the experiment were placed in the swim chamber and subjected to various velocities, with the optimal velocity being chosen as swift enough to (A) induce swimming and (B) prompt the cessation of swimming within several minutes. The inner-mesh lining of the swim tunnel was necessary to prevent subjects from attaching to the side of the tube, but rendered the manufacturer’s regression equation linking power units of the swim-tunnel impeller to water velocity imprecise. COVID-19 restrictions prevented empirical water velocity calibration; however, passive particulates were moving through the tunnel at ~ 15 cm s−1 (approximately 1 body lengths/s for our subjects). Consequently, this assay was used to compare differences between groups, not to precisely estimate time-to-exhaustion in tagged and untagged sea lamprey.

After a 3-min acclimation period with no water flow, animals were induced by water flow to freely swim against a current. As the test progressed, animals would become impinged on the mesh barrier at the downstream end of the tunnel. While the animal maintained position in the flow, actively swimming against the current, elapsed time was recorded. If the animal became impinged, but maintained active swimming motions (tail-beating) resisting the current, elapsed time continued to be recorded. In the first instance an animal became impinged with a continued absence (10 s) of anguilliform movement resisting the water current, a Smith-Root backpack electro-shocking unit sent a brief 12 V, 1 Hz, and 5% duty cycle electrical current into the rear of the swim tunnel for up to ten seconds. This mild voltage was intended to irritate rather than stun the animal so swimming could recommence with full musculature control. Trials continued for animals induced to resume swimming within ten seconds of continuous electrical current. We deemed exhaustion occurred (trial completed) when either of two events occurred: (1) an absence of anguilliform swimming for ten seconds after the second impingement, or (2) no resumption of active swimming within ten seconds of electrical current after the initial impingement. If the animal continued to swim for 60 min without exhaustion, the trial was ended.

Data analysis (tag retention and survival)

Tanks were monitored daily for shedding of transmitters and mortalities by visually searching the substrate for transmitters and immobile animals. No burrowing activity was observed during this experiment; therefore, a visual inspection of the surface ensured that any mortalities or shed tags were seen and that the animals were minimally disturbed. Immobile animals were tapped with a net to determine if the subject was non-responsive. Dead animals and shed tags were removed from tanks and recorded daily until the end of this study.

Daily mortality data recorded from a 61-day holding period were used to generate Kaplan–Meier survival curves for control and tagged animals; a log-rank test (Mantel-Cox) was used to assess the difference in survival among the two groups. The p-value resulting from a Log-Rank test determines the level of significant difference in survival between tagged and control groups for the entirety of the experiment. This potentially creates “blind spots” in the analysis, which would prevent the identification of critical periods when large changes in mortality occurred between groups (e.g., immediate post-surgical or post-swim testing mortality). These periods have high informative value for future studies, especially for field studies seeking to use this technology (e.g., setting the post-surgical holding period if surgery related mortality is delayed by a few days). We further explored the data to determine if and when any blind spots occurred. To do this, we repeatedly subset the survival data, creating a new dataset for each day (dayi) including data from dayi until the end of the experiment. This nullified deaths prior to dayi, thereby establishing dayi as day 0 in the analysis for each dataset (i.e., the first subset equals the whole dataset, from day 1 onward, the second contains only day 2 onward, etc.). A Log-Rank test of each dataset produced p-values for the period following dayi that were evaluated by generating Kaplan–Meier survival curves. This progressive p-value provided insight by illustrating the pattern of the significance level in the difference between the tagged and control groups survival through time. These p-values were used for qualitative visualization only, with statistical significance hinging on the full dataset log-rank test.

Upon reviewing the survival curves of the tagged and control groups, it was apparent there was a time-dependent component in the survival curves. Mortalities appeared to occur in punctuated events within a group (i.e., were not homogeneous throughout the 61-day holding period). Kaplan–Meier-based survival analyses restricted to only one covariate are unable to capture this effect. Therefore, additional analyses were performed using Cox-Proportional Hazard (Cox model) regression (R package “survival” v.2.44–1.1, [69, 70]) with a time-dependent covariate (TDC) to determine if tag implantation was associated with an increased risk of mortality, and to what extent differences changed over time. The time-dependent covariate was a binomial categorization that was set to one for tagged and zero for control animals surviving to Day 32, which gave the model an explicit value to measure the significance and magnitude of a tag effect from days 32 to 61. This model was compared using Akaike Information Criterion (AIC) and heteroscedasticity of the Schoenfeld residuals to a null model of one covariate (Group) to assess the role of the TDC in the model’s fit. We used a counting process form data frame to build the TDC into the dataset, constructed using methodology similar to that of Zhang et al. [78] and Therneau [70]. An individual animal is represented by one row if mortality was experienced prior to Day 32 and two rows if mortality was experienced after Day 32 or the animal survived the experiment in entirety (covariate set to 1 if tagged, 0 if control). Day 32 was chosen as the break point as it coincided with the first mortality event recorded in the control group. Schoenfeld residuals were examined to determine if the assumptions of a proportional hazard model were met. In both the Kaplan–Meier and Cox-Proportional Hazard analyses, tagged animals that shed their tag were not right-censored at the time of the shed event to preserve statistical power through sample size. Finally, animals experiencing early mortalities (< 5 days) were compared to animals surviving > 5 days through separate logistic regression with the three size metrics as predictors. Source (Great Lakes vs. Anadromous) was the additional covariate.

Data analysis (wound healing and condition)

The two wound healing scores, and the composite score, were assessed through time via separate non-parametric Friedman tests, with the assessment score as the dependent variable and numeric assessment (1–5) as the predictor. Post hoc Wilcoxon pairwise rank sum tests using Bonferroni correction explored the significance of relationships between wound scores and numeric assessment [55]. Relationships between each of three size measurements (Total Length in mm, mass in g, and Fulton’s Condition Factor; Eq. 2; Ricker [57]) and each wound score were explored through simple linear regression [55]:

$$\frac{mass(g)}{{Total \,Length(mm)}^{3}* {10}^{6}}.$$
(2)

Data analysis (swim performance)

A two-way unbalanced ANOVA [22] was performed to determine the relationship between maximum observed burst speed and treatment group (Control vs. Tagged) with source (Great Lakes vs. Anadromous) as a covariate. Separate simple linear regression models determined the relationship of size and composite wound score to maximum burst speed for the tagged animals. A time-to-event analysis was used to compare the tagged and control group’s time-to-exhaustion probabilities by evaluating their respective Kaplan–Meier survival curves via log-rank test (R package “survival” v.2.44–1.1, [70], survival = continuing to swim). Plotted through time (log-seconds) on the x-axis, survival probabilities dropped on the y-axis as members of the respective groups reached their time of exhaustion. The log-rank test assessed the statistical difference between the cumulative exhaustion probabilities of each group. To assess significance, the log-rank test statistic, where the expected value is the product of the risk of event (number of exhaustion events/number unexhausted) and the number of unexhausted animals at the given time, is compared to the critical value (chi-square) for one degree of freedom [25]. Additionally, an unbalanced two-way ANOVA tested the null hypothesis that there was no significant difference between swim-to-exhaustion time and covariates group and source.

Results

Four days into the experiment one holding tank containing 14 tagged and 1 control animal from the anadromous tank experienced 24 h without inflowing water or aeration. These animals were used in analyses prior to that point and censored from any analysis after that date [41].

Tag retention and survival

Eleven of the 59 tagged animals (19%) shed their tags throughout the experiment. Two shedding events occurred when handling animals for endurance swim performance tests. The majority (5/9) of the remaining shedding events occurred within the first four days post-implantation, and all animals that shed their tags had wounds that were completely open with no apposition in the initial wound assessment on Day 4 (wound closure score = 4). The four additional shedding events were recorded on days 19, 20, and 29 (two individuals).

Twenty-five sea lamprey died during the experiment (8 in the control group, 17 implanted with ELAT transmitters). The tagged group experienced increased mortality at the commencement of the experiment, with five animals dying within four days post-surgery (Days 1, 1, 2, 4, 4). Tag burden (tag weight/body weight, %) averaged 1.87 ± 0.036 across the entire tagged group, with mortalities distributed across the range (Fig. 2), but higher in animals with a tag burden above 2% (6/15, 40% mortality as a group) vs. those at or below 2% (11/44, 25% mortality as a group). However, no size metric (TL, Mass, Condition Factor) nor source was a significant predictor of early mortality (logistic regression, all p-values > 0.1). No mortalities were observed in the control group until Day 31, with a total 8 mortalities accumulating gradually over the following 25 days.

Fig. 2
figure 2

Frequency histogram of tag burden in a group of 56 sea lamprey fitted with an ELAT transmitter. Animals that died (n = 17) are indicated in orange. The upper box plot reports the median (vertical line), 25/75th percentiles (box), and 10/90th percentiles (whiskers) with outliers (filled circles)

The cumulative survival probabilities (mean ± 95% CI) were 0.849 ± 0.092 and 0.647 ± 0.124 for the control and tagged groups, respectively. The observed difference in survival was significant (Mantel-Cox log-rank test, p = 0.011). Examination of the Kaplan–Meier survival curves and the progressive p-value (Fig. 3) suggests distinct periods of elevated mortality in both groups. The principal differences in survival between tagged and control groups appear in the four-day period immediately post-surgery (whereafter, the progressive p-value loses statistical significance), and Days 18–30, roughly corresponding to the period following the burst swim tests (Day 20, Kaplan–Meier curves, Fig. 3). Prior to date of the swim-to-exhaustion test (Day 30), no mortalities were recorded in the control group. Following this test, both groups exhibited similar reductions in survival. The Cox-Proportional hazard model suggested tagged animals were 7.4 times more likely to experience mortality in the first 32 days of the experiment vs. control animals (coefficient = 2.00, Hazard ratio = 7.4, p = 0.008). The time-dependent covariate coefficient was − 1.80 (p = 0.065, 95% CI -3.7—0.11, se = 0.97) which implies that after 32 days, the hazard ratio in the tagged group was reduced to of 1.22 \(\left({e}^{2.00-1.80}\right)\). The ΔAIC (1.86) relative to the null model supports the assertion that the difference of associated risk through time between the two groups is not homogeneous. Moreover, the global Schoenfeld residuals of the time-dependent covariate model (p = 0.63) present less heteroscedasticity than those of the null model (p = 0.10).

Fig. 3
figure 3

Kaplan–Meier survival curves with 95% confidence intervals for control (red) and tagged (blue) groups with the overall significance (Mantel-Cox log-rank test). A (x) symbol denotes a right-censorship event in the respective group. The progressive P-value (black line, right Y-axis) portrays the significance level associated with Log-Rank tests throughout each day of the experimental period. Each point on the line corresponds to the p-value for a Log-Rank test, adjusting the experimental period to begin on that day. The purple dashed line represents a significance level of 0.05

Wound healing and condition

No animals received a composite score of 0 indicating a completely healed wound by the end of the study period (Fig. 4). One animal had a completely closed wound (wound closure = 1) on the final wound assessment, but mild inflammation was present (wound inflammation = 2). The initial composite wound score (Day 4) was unrelated to total length (separate univariate linear regressions: p = 0.07, Adjusted R2 = 0.04, SE = 0.06) or condition factor (p = 0.27, Adjusted R2 = 0.005, SE = 2.5) but negatively associated with mass (p = 0.02, Adjusted R2 = 0.08, SE = 0.5). Each additional gram of mass was associated with an improvement of 1.2 composite score units (i.e., a 12% improvement in wound condition). Wound closure scores improved throughout the experimental period, with the greatest improvement occurring between the first and second observation dates (Days 4 and 20, post-surgery, mean change = 0.47 score units ± 0.24). Animals surviving the full experiment showed a slight decline on average in wound inflammation from Days 4 to 61 (0.17 score units ± 0.15), but inflammation scores improved between Days 4 and 20, with the same mean as wound healing. A Friedman rank sum test revealed time was a significant predictor for each of the wound scoring metrics (Table 2). Post hoc pairwise Wilcoxon Rank Sum Tests revealed a significant improvement in wound scores performed after Day 20 (Table 2).

Fig. 4
figure 4

Wound assessment scores. A Distribution of wound closure scores for animals remaining in the study at the time of assessment. B Distribution of wound inflammation scores for animals remaining in the study at the time of assessment. A & B 1 = Violet, 2 = Yellow, 3 = Blue, 4 = Green C Box plots of composite wound scores for animals remaining in the study at the time of assessment, showing median (black line), and maximum and minimum values (lines)

Table 2 Cumulative wound scores (closure, inflammation, composite) Friedman Rank Sum comparison through time, with total p-values and those of pairwise Wilcoxon Rank Sum Tests between successive assessments

Swim performance

Burst speed differed between control and tagged groups (ANOVA: p = 0.003, F = 9.51), but was not related to origin of the animals (Great Lakes vs Atlantic drainages; p = 0.33, F = 0.9). Mean maximum burst speed in cm/s (± 1 SE) was 26.22 ± 1.81 for the tagged group and 33.83 ± 1.42 for the untagged group (Fig. 5). The greatest burst speed recorded in each group was 57.02 cm/s for tagged animals and 56.88 cm/s for control animals.

Fig. 5
figure 5

A Violin plots of recorded maximum burst speeds for control and tagged sea lamprey, overlaid with box plots exhibiting median (red, black bar), 25% and 75% quartiles (boxes) and non-outlying maximum and minimum speeds (lines, within 1.5 × interquartile range). The raw data are jittered on the X-axis to ease interpretation. B Scatterplot and regression of the maximum burst swim speed vs. Day 4 composite wound score. For reference, the mean (± 95% CI) maximum burst speed for the control group is plotted at a value of 0 on the x-axis (black), but was not included in the regression analysis

In the tagged group, condition and body size were not related to maximum burst swim speed (linear regression: condition factor, p = 0.36; total length, p = 0.36; mass, p = 0.97). There was a significant negative correlation between observed maximum burst speed and each of the wound scores recorded on Day 4 post-surgery (linear regression: wound closure, p = 0.031, adjusted R2 = 0.12; wound inflammation, p = 0.022, adjusted R2 = 0.12; composite wound score, p = 0.007, adjusted R2 = 0.18) (Fig. 5B). However, the composite wound score recorded on the day of the burst swim test (Day 20) was not a significant predictor of maximum burst swim speed (linear regression: p = 0.22, adjusted R2 = 0.02).

Three animals in each of the tagged and control groups could not be induced to swim in the swim tunnel and were not evaluated for swim-to-exhaustion. Two additional animals in each group exhibited lethargy upon retrieval from their holding tank and were censored from the analysis. Sample sizes used in analysis were n = 25 tagged and n = 27 control. One animal from the control group completed the full trial without exhausting, whereas all other animals in both groups exhausted prior to 60 min. Time-to-exhaustion did not significantly differ between control and tagged groups (ANOVA: p = 0.32, F = 1.02) and was not related to origin of the animals (Great Lakes vs Atlantic drainages; p = 0.17, F = 1.96). Median time-to-exhaustion in min. for tagged and control groups were 1.45 (Range = 0.73 – 26.28) and 2.35 (Range = 0.53–60), respectively (Fig. 6A). Similarly, no significant difference in the Kaplan–Meier survival curves was detected (log-rank test, p = 0.18; Fig. 6B).

Fig. 6
figure 6

A Violin plots of recorded time-to-exhaustion for control and tagged sea lamprey, overlaid with box plots exhibiting median (red, black bar), 25% and 25% quartiles (boxes) and maximum and minimum values for non-outlying results (lines). Data are jittered on the X-axis to ease interpretation. B Survival (exhaustion) curves for control (red) and tagged (blue) animals in Swim-to-Exhaustion tests (p = 0.18)

Discussion

In telemetry studies of juvenile or small fishes, negative impacts of the transmitter on behavior, swim performance, and physiology may result in erroneous conclusions regarding movement trends, behavioral tendencies, and survival rates. Accurate estimation of these components for out-migrating juveniles is critical to understanding the efficacy of conservation measures, such as the installation of fish passage devices and the restoration of spawning and rearing habitats. The results of this study suggest implantation of an ELAT transmitter into newly transformed sea lamprey results in survival and tag retention akin to that reported for similar sized PIT tags. The 61-day survival of our tagged subjects was 71%, well within the range reported in previous studies with larval and juvenile lampreys 25–100%; [14, 29, 49,50,51, 63]. However, survival was significantly lower in the tagged animals (vs control), with no effect statistically attributable to measures of animal length, mass, condition, or population of origin (Great Lakes vs. Atlantic drainage). There was some evidence that individuals with higher tag burdens were more likely to experience mortality, but given the small sample sizes, that finding should be interpreted with caution.

The significant time-dependent covariate in the Cox model and visualization of the progressive p-values suggest the statistical difference between the two groups was heavily weighted by mortalities occurring in the first four days of the experiment. Three lines of evidence suggest the immediate post-surgical mortalities and tag losses were likely caused by imperfect surgical technique, the acute stress of the surgery, or both. First, mortalities and tag losses were distributed across a range of body size, a pattern consistent with a surgeon/surgery effect rather than the consequences of tag internalization [4, 10]. Prior work noted nicking the gastrointestinal tract or other internal organs, and/or protrusion of the intestine through the incision, was associated with rapid mortality (< ~ 7 days) in juvenile lampreys [13, 49, 63] as observed here. Although we did not measure it explicitly, we observed visual evidence of discoloration around the wound site that is consistent with descriptions of internal hemorrhaging provided by Dawson et al. [13, 14]. Two of 5 animals experiencing early mortality in this study exhibited protruding intestines and visual evidence of internal hemorrhaging near the wound site. Second, all tagged animals that survived the four days post-surgery also survived the following two weeks leading up to the first swim test and exhibited clear improvement in wound condition. Third, the time to stage-IV anesthesia was long relative to other anesthetics used with lampreys, and longer than recommendations of all anesthetics (Summerfelt & Smith 1990). Data reported by Christiansen et al. [7] regarding Pacific lamprey transformers suggest a 70 mg/L dose of AQUI-S should result in full sedation in six to nine minutes, we observed a mean time to Stage-IV sedation of 17 min. Strikingly, the time to recovery reported by Christiansen et al. (3–8 min) was considerably faster than we observed (a mean of 102 min). Because our goal was to investigate the use of the ELAT in field studies, our holding tanks and anesthetic bath were maintained at the typical water temperature experienced by an out-migrating transformer in Great Lakes streams (4.5 °C), a temperature substantially colder than reported by Christiansen et al. (12° C). In fishes, colder water temperature can substantially increase the time to uptake and recovery from various anesthetic agents, leading to increases in the rate of post-surgical complications and death [64, 77]. Our results suggest that use of alternative anesthetics (e.g., MS-222) should be considered when performing surgeries in cold water on juvenile sea lamprey. In addition, our study was required to follow Federal Investigational New Animal Drug (INAD) guidelines for the use of AQUI-S (i.e., dosage and exposure time) contained in INAD #11-741 [71]. At a minimum, anesthesia management procedures for cold-water use of AQUI-S need to be better developed for use in juvenile sea lamprey. Future studies with greater availability of subjects should also include anesthesia-only and anesthesia+surgery sham controls to better disentangle the relative effects of each factor in survival and performance tests.

It is noteworthy that wound healing was slower in our subjects when compared to the two available studies of PIT tagging transforming or transformed juvenile lampreys despite a similar incision size (3 mm). In the present study, one surviving tagged fish exhibited a fully closed wound at the conclusion of the 61-day study period (2.2%), with the majority of improvement in the wounds occurring between Days 4 and 20. Simard et al. [63], studying metamorphosing anadromous (Fort River, MA) and land-locked (Lake Champlain) sea lamprey, report post-surgical skin closure rates of approximately 15–20% after 60 days of healing, though rates varied across populations (12 mm PIT tag surgically inserted, ~ 5 °C). Mesa et al. [46] report only 2.2% of transformed Pacific lamprey held at 9–18 °C exhibited signs of poorly healed or infected incisions 40 days following implantation of a 9 mm PIT tag. At Day 20, our results were akin to those reported by Moser et al. [49] for Pacific lamprey ammocoetes tagged with an 8.4 mm PIT tag and held in the field for 15 days at 6– 11 °C (i.e., no healed wounds). Despite frequent reports of fungal infections in other laboratory studies using juvenile lampreys [46, 50, 61, 63], infections did not arise in our subjects. The water temperature used here was predominantly lower than that of the studies mentioned, in most cases by > 4 °C (~ 5–18 °C); thus, the use of metamorphosed animals held at low water temperature may have both slowed the rate of wound healing and inhibited fungal infection. However, it is likely that use of this transmitter in out-migrating lamprey within Great Lakes watersheds will involve animals experiencing similar water temperatures. Slow healing may be a significant concern in field studies, where more vigorous swimming and burial behaviors may further slow wound closure and may also increase the likelihood of tag loss (vs. that observed in tank studies). Moser et al. [49] report slower healing for juvenile Pacific lamprey held in a natural stream vs. those held in laboratory tanks. Options for defending against this outcome in field studies include the addition of a suture to hasten wound closure and better retain the transmitter, performing surgeries under warmer water in the laboratory (with appropriate acclimation periods), or simply targeting studies to warmer times of year.

Implantation with an ELAT transmitter negatively affected swim performance. Maximum burst swim speed in tagged animals was 7.6 cm/s (22%) less than untagged animals, regardless of source population, suggesting a substantial reduction in the ability to perform high-speed acceleration. Interestingly, the composite wound score taken four days post-surgery was a significant predictor of burst swim speed measured more than two weeks later, whereas the score on the day of the test was not. Further, the predicted burst swim speed for a wound score of zero (completely healed) was 36.9 cm/s, similar to the mean burst speed observed in untagged animals (33.8 ± 1.4 cm/s). Mueller et al. [50] report a 6 cm/s reduction in the burst swim speed of actively out-migrating Pacific lamprey, although the relative reduction vs. untagged animals (7%) was substantially lower than in our study, despite using similarly sized subjects. The difference is largely attributable to that the recorded burst swimming speeds of Pacific lamprey are much higher which appears consistent across studies (means of 51–76 cm/s), e.g., [12, 67]. It is important to note that these studies were also performed at higher temperatures (12 and 20 °C, respectively) and that swim performance may be strongly affected by temperature in fishes [6, 74, 75]. Conversely, and despite a 38% decrease in median value, we failed to detect a significant negative effect of the transmitter in the time-to-exhaustion swim test. Why the effect of the tag manifested more strongly during burst vs. steady swimming is of interest, as burst or “escape” swimming is more associated with successful passage through hydraulic challenges and the avoidance of predators during out-migration [21]. When a juvenile sea lamprey accelerates from rest, it generates a high amplitude body bend (double that of steady swimming) that propagates over 75% of the body length [20]. It is plausible there is an amplitude threshold where the transmitter becomes impinged on the viscera, restricting the ability to generate thrust, and aggravating any internal injury associated with the tag implantation surgery. That threshold is more likely to be breached during high acceleration swimming, and conceivably may affect the estimation of natural mortality rates in field studies with small, narrow-bodied anguilliform fishes.

A core assumption in telemetry studies is that tagged individuals are characteristic of the natural population in terms of size, sex, movement characteristics, and fate [30]. Our experimental animals were somewhat smaller than the historic mean size observed in the Great Lakes, and did not encompass the full size range of available fish. Although we did not observe significant statistical effects of animal size on the measured outcomes, the surgeon reported difficulty inserting the transmitter into subjects with smaller body cavities, and there was some evidence that individuals with tag burdens in excess of two percent of body mass suffered greater mortality. Several lamprey also exhibited light discoloration posterior to the implanted tag soon after surgery, a symptom of blood impediment through the dorsal aorta caused by pressure from the tag [13] that appeared prior to recovery from anesthesia. However, it was not a consistent indicator of impending mortality. Two of 5 animals displaying this symptom no longer exhibited discoloration by the first wound assessment on Day 4 and survived the length of the study, perhaps due to a posterior shift in tag position as posited by Moser et al. [49]. It was also notable that in both swim tests there was substantial variation in performance among individuals that was unrelated to measures of body size or population of origin. A number of studies demonstrate both burst and steady swimming speeds, and the underlying aerobic scope for activity, are repeatable traits within individual fish during swim performance tests, e.g., [8, 23, 36, 52]. Understanding consistent individual differences in tag effects on swim performance is likely an important metric to understanding whether and how the effects of the transmitter on swimming will manifest in studies designed to measure movement rates and survival in the wild [73]. Further, out-migrating juvenile lampreys exhibit strong nocturnal patterning in downstream movements, burrowing into sediment or seeking shelter during daylight hours [12, 47, 54]. In a study with smaller (means of < 120 mm) larval sea lamprey implanted with 8 or 9 mm PIT tags, the tag significantly increased both the time (2.3 × longer) and effort (defined as number of “stops” while burrowing, 1 additional stop on average) necessary for the animal to successfully burrow (vs. untagged fish; [13]).

It is routine for tag effect studies to recommend minimum size limits for use of animals in telemetry work, based on the desire to satisfy the criteria of minimal tag effect on the animal’s performance. However, this practice conflicts with the need to use representative animals, as demographic measures, such as mortality, are highly size dependent in juvenile fishes [2, 65]. In other words, preferential use of larger subjects that are less likely to suffer tag effects (i.e., non-probability sampling) may bias estimates of vital rates made using those animals in field studies. Conceivably, this may be compensated in the modeling of mortality rates from telemetry studies by using size-dependent estimates of tag-induced mortality. Similar studies to ours in larval and juvenile Pacific and sea lamprey offer minimum size recommendations of 120–150 mm total length, based on the relationship between mortality and length observed in the lab [46, 50, 61, 63]. While there was no significant correlation of any size measurement and mortality in our tagged animals, larger mass, but not total length, was associated with an improved composite wound score, suggesting lesser injury during the surgery. Post-metamorphosis, the intestines of a sea lamprey range laterally throughout the body of the animal, lessoning the available space between the viscera and the body wall [42]. We suggest that body girth at the tag location, or a similar metric (e.g., ratio of body depth to length), or a measure of condition factor, may prove more suitable metrics to set the minimum threshold for implantation. These metrics require a larger dataset, perhaps integrated across previous and future studies to solidify a concrete threshold. Our findings support a ratio of 3.33:1 (length (mm) to mass (mg)), above a minimum length and mass of 150 mm and 45 mg, should ensure high survival. We further recommend a brief holding period (2–4 days) after recovery from anesthesia and prior to release to assess wound condition and potential early tag loss. If the wound does not appear to have inflammation over a moderate level and/or a complete lack of apposed tissue around the incision, it is likely that the animal will experience limited physical incapacitation compared to an untagged conspecific. Additionally, although not experienced in our study, prior work with juvenile lampreys suggest fungal infection is a risk whenever the fish are held in tanks, especially at warmer temperatures. A prophylactic anti-fungal regime may be needed in those conditions.

Conclusions

The ELAT transmitter has now proven feasible in laboratory and field settings across a range of anguilliform species with relevant behavioral unknowns, such as the regionally endangered American brook lamprey (Lethenteron appendix) of the United States Eastern Seaboard and the critically endangered European eel (Anguilla Anguilla) (Department, N.H. Fish and Game [18], The IUCN Red List of Threatened Species [68]) [16, 17, 39]. Our findings demonstrate its utility for juvenile sea lamprey. The ability to gain insight into the fate of migrants as a function of the features of the rivers that regulate mortality (e.g., predator density, habitat type, etc.) is likely to lead to the development of new control methods to target sea lamprey in their last phase before the initiation of parasitism [33]. In addition, it allows for the investigation of the effects of habitat degradation and barriers to downstream migration (dams and weirs, water intake structures), and overharvesting; the principal threats to several migratory lampreys in their native range [28, 35, 40, 44, 48]. Ultimately, the ELAT transmitter has the potential to open the door to a greater understanding of the out-migration parasitic lampreys, a critical and poorly understood stanza in the lives of these consequential and enigmatic creatures.

Availability of data and materials

All datasets used in this publication are available upon request to the corresponding author.

Abbreviations

ELAT:

Eel–lamprey acoustic transmitter

PIT:

Passive Integrated Transponder

JSATS:

Juvenile Salmon Acoustic Telemetry System

PRI:

Pulse Rate Interval

USFWS:

United States Fish & Wildlife Service

INAD:

Investigational New Animal Drug

TDC:

Time-Dependent Covariate

AIC:

Akaike Information Criterion

References

  1. Baer J, Hartmann F, Brinker A. Abiotic triggers for sea and river lamprey spawning migration and juvenile outmigration in the River Rhine Germany. Ecol Freshwater Fish. 2018;27(4):988–98.

    Article  Google Scholar 

  2. Bailey KM, Houde ED. Predation on Eggs and Larvae of Marine Fishes and the Recruitment Problem. In: Blaxter JHS, Southward AJ, editors. Advances in Marine Biology. Cambridge: Academic Press; 1989. p. 1–83.

    Google Scholar 

  3. Baker CF, Reeve K, Baars D, Jellyman D, Franklin P. Efficacy of 12-mm half-duplex passive integrated transponder tags in monitoring fish movements through stationary antenna systems. North Am J Fish Manag. 2017;37(6):1289–98.

    Article  Google Scholar 

  4. Brown RS, Oldenburg EW, Seaburg AG, Cook KV, Skalski JR, Eppard MB, et al. Survival of seaward-migrating PIT and acoustic-tagged juvenile Chinook salmon in the snake and Columbia rivers: an evaluation of length-specific tagging effects. Anim Biotelem. 2013;1(1):1–13.

    Article  Google Scholar 

  5. Burkett DP, Barber JM, Steeves TB, Siefkes MJ. Sea lamprey control 2020–2040: charting a course through dynamic waters. J Great Lakes Res. 2021;2021(47):S809–14.

    Article  Google Scholar 

  6. Cano-Barbacil C, Radinger J, Argudo M, Rubio-Gracia F, Vila-Gispert A, García-Berthou E. Key factors explaining critical swimming speed in freshwater fish: a review and statistical analysis for Iberian species. Sci Rep. 2020;10(1):18947.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Christiansen HE, Gee LP, Mesa MG. Anesthesia of juvenile pacific lampreys with MS-222, BENZOAK, AQUI-S 20E, and Aquacalm. North Am J Fish Manag. 2013;33(2):269–76.

    Article  Google Scholar 

  8. Claireaux G, Handelsman C, Standen E, Nelson JA. Thermal and temporal stability of swimming performance in the European Sea Bass. Physiol Biochem Zoology. 2007;80(2):186–96.

    Article  Google Scholar 

  9. Clemens B. A call for standard terminology for lamprey life stages. Fisheries. 2019;44(5):243–5.

    Article  Google Scholar 

  10. Cooke SJ, Graeb BDS, Suski CD, Ostrand KG. Effects of suture material on incision healing, growth and survival of juvenile largemouth bass implanted with miniature radio transmitters: case study of a novice and experienced fish surgeon. J Fish Biol. 2003;62(6):1366–80.

    Article  Google Scholar 

  11. Crossin GT, Heupel MR, Holbrook CM, Hussey NE, Lowerre-Barbieri SK, Nguyen VM, et al. Acoustic telemetry and fisheries management. Ecol Appl. 2017;27(4):1031–49.

    Article  PubMed  Google Scholar 

  12. Dauble DD, Moursund RA, Bleich MD. Swimming behaviour of juvenile pacific lamprey Lampetra tridentata. Environ Biol Fish. 2006;75(2):167–71.

    Article  Google Scholar 

  13. Dawson HA, Potts DD, Maguffee AC, O’Connor LM. Feasibility of passive integrated transponder technology to study in situ movements of Larval Sea lamprey. J Fish Wildlife Manag. 2015;6(1):71–82.

    Article  Google Scholar 

  14. Dawson HA, Quintella BR, Almeida PR, Treble AJ, Jolley JC. The ecology of larval and metamorphosing lampreys. In: Docker MF, editor. Lampreys: biology, conservation and control: volume 1. Springer: Dordrecht; 2015. p. 75–137.

    Chapter  Google Scholar 

  15. Deng ZD, Carlson TJ, Li H, Xiao J, Myjak MJ, Lu J, et al. An injectable acoustic transmitter for juvenile salmon. Sci Rep. 2015;5(1):8111.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Deng ZD, Li H, Lu J, Xiao J, Myjak MJ, Martinez JJ, et al. An acoustic micro-transmitter enabling tracking of sensitive aquatic species in riverine and estuarine environments. Cell Rep Phys Sci. 2021;2(5):100411.

    Article  CAS  Google Scholar 

  17. Deng, Z.D., Titzler, P.S., Mueller, R.P., Lu, J., Martinez, J.J., Fu, T., Li, H., Cable, J.S., Deters, K.A., Hubbard, J.M., and Colotelo, A.H. Pilot field trial of the juvenile lamprey/eel tag and RME plan to guide future juvenile pacific lamprey acoustic telemetry studies. 2018

  18. Department, N.H Fish and Game. Endangered and Threatened Wildlife of NH | Nongame | New Hampshire Fish and Game Department. https://www.wildlife.state.nh.us/nongame/endangered-list.html. Accessed 17 Mar 2021.

  19. Domenici P, Kapoor BG. Fish Locomotion: An Eco-ethological Perspective. Science Publishers; 2010.

  20. Du Clos KT, Dabiri JO, Costello JH, Colin SP, Morgan JR, Fogerson SM, et al. Thrust generation during steady swimming and acceleration from rest in anguilliform swimmers. J Exp Biol. 2019;22:jeb212464.

    Article  Google Scholar 

  21. Evans TM, Wagner CM, Miehls SM, Johnson NS, Haas TF, Dunlop E, et al. Before the first meal: the elusive pre-feeding juvenile stage of the sea lamprey. J Great Lakes Res. 2021;47:580–9.

    Article  Google Scholar 

  22. Fox, J, Sanford Weisberg. An {R} Companion to Applied Regression, Third Edition. Thousand Oaks CA: 2019; Sage. https://socialsciences.mcmaster.ca/jfox/Books/Companion/

  23. Fuiman LA, Cowan JH Jr. Behavior and recruitment success in fish larvae: repeatability and covariation of survival skills. Ecology. 2003;84(1):53–67.

    Article  Google Scholar 

  24. Furey NB, Armstrong JB, Beauchamp DA, Hinch SG. Migratory coupling between predators and prey. Nat Ecol Evol. 2018;2(12):1846–53.

    Article  PubMed  Google Scholar 

  25. Goel MK, Khanna P, Kishore J. Understanding survival analysis: Kaplan-Meier estimate. Int J Ayurveda Res. 2010;1(4):274–8.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Goodman DH, Reid SB, Som NA, Poytress WR. The punctuated seaward migration of Pacific lamprey (Entosphenus tridentatus): environmental cues and implications for streamflow management. Can J Fish Aquat Sci. 2015;72(12):1817–28.

    Article  Google Scholar 

  27. Griffioen AB, van Keeken OA, Hamer AL, Winter HV. Passage efficiency and behaviour of sea lampreys (Petromyzon marinus, Linnaeus 1758) at a large marine–freshwater barrier. River Res Appl. 2022;38(5):906–16.

    Article  Google Scholar 

  28. Hansen MJ, Madenjian CP, Slade JW, Steeves TB, Almeida PR, Quintella BR. Population ecology of the sea lamprey (Petromyzon marinus) as an invasive species in the Laurentian Great Lakes and an imperiled species in Europe. Rev Fish Biol Fisheries. 2016;26(3):509–35.

    Article  Google Scholar 

  29. Hanson KC, Barron JM. Evaluation of the effects of marking pacific lamprey ammocoetes with visual implant elastomer, coded wire tags, and passive integrated transponders. Trans Am Fish Soc. 2017;146(4):626–33.

    Article  Google Scholar 

  30. Hondorp DW, Holbrook C, Krueger CC. Effects of acoustic tag implantation on lake sturgeon Acipenser fulvescens: lack of evidence for changes in behavior. Anim Biotelemetry. 2015;3(1):1–3.

    Article  Google Scholar 

  31. Howe EA, Marsden JE, Donovan TM, Lamberson RH. A life cycle approach to modeling sea lamprey population dynamics in the Lake Champlain basin to evaluate alternative control strategies. J Great Lakes Res. 2012;38:101–14.

    Article  Google Scholar 

  32. Hussey NE, Kessel ST, Aarestrup K, Cooke SJ, Cowley PD, Fisk AT, et al. Aquatic animal telemetry: a panoramic window into the underwater world. Science. 2015;348(6240):1255642.

    Article  PubMed  Google Scholar 

  33. Johnson NS, Miehls S. Guiding out-migrating Juvenile Sea Lamprey (petromyzon Marinus) with Pulsed Direct Current. River Res Appl. 2014;30(9):1146–56.

    Article  Google Scholar 

  34. Jones LM, Irwin BJ, Hansen GJA, Dawson HA, Treble AJ, Liu W, et al. An Operating Model for the Integrated Pest Management of Great Lakes Sea Lampreys. Open Fish Sci J. 2009;2:59–73.

  35. Jones ML. Toward improved assessment of sea lamprey population dynamics in support of cost-effective sea lamprey management. J Great Lakes Res. 2007;3:35–47.

    Article  Google Scholar 

  36. Kolok AS, Farrell AP. Individual variation in the swimming performance and cardiac performance of Northern Squawfish Ptychocheilus oregonensis. Physiol Zool. 1994;67(3):706–22.

    Article  Google Scholar 

  37. Krueger DM, Rutherford ES, Mason DM. Modeling the influence of parr predation by walleyes and brown trout on the long-term population dynamics of chinook salmon in Lake Michigan: a stage matrix approach. Trans Am Fish Soc. 2013;142(4):1101–13.

    Article  Google Scholar 

  38. Larson JH, Trebitz AS, Steinman AD, Wiley MJ, Mazur MC, Pebbles V, et al. Great Lakes Rivermouth ecosystems: scientific synthesis and management implications. J Great Lakes Res. 2013;39(3):513–24.

    Article  Google Scholar 

  39. Liedtke TL, Lampman RT, Deng ZD, Beals TE, Porter MS, Hansen AC, et al. Movements of juvenile Pacific Lamprey (Entosphenus tridentatus) in the Yakima and Columbia Rivers, Washington, 2018—A pilot study using acoustic telemetry. Reston, VA: U.S. Geological Survey; 2019; p 40. (Open-File Report; vols. 2019–1058). Report No.: 2019–1058.

  40. Maitland PS. Ecology of the River, Brook and Sea Lamprey. Conserving Natura 2000 Rivers Ecology Series No. 5. Peterborough: English Nature; 2003.

  41. Mantel N. Evaluation of survival data and two new rank order statistics arising in its consideration. Cancer Chemother Rep. 1966;50(3):163–70.

    CAS  PubMed  Google Scholar 

  42. Manzon RG, Youson JH, Holmes JA. Lamprey metamorphosis. In: Lampreys: biology, conservation and control. 139–214 M.F. Docker. Fish and Fisheries Series. Springer; 2015.

  43. Marsden JE, Siefkes MJ. Control of Invasive Sea Lamprey in the Great Lakes, Lake Champlain, and Finger Lakes of New York. In: Docker MF, editor. Lampreys: biology, conservation and control: volume 2. Dordrecht: Springer; 2019. p. 411–79.

    Chapter  Google Scholar 

  44. Mateus C, Rodríguez-Muñoz R, Quintella B, Alves M, Almeida P. Lampreys of the Iberian Peninsula: distribution, population status and conservation. Endanger Species Res. 2012;16:183–98.

    Article  Google Scholar 

  45. McMichael GA, Eppard MB, Carlson TJ, Carter JA, Ebberts BD, Brown RS, et al. The Juvenile Salmon Acoustic telemetry system: a new tool. Fisheries. 2010;35(1):9–22.

    Article  Google Scholar 

  46. Mesa MG, Copeland ES, Christiansen HE, Gregg JL, Roon SR, Hershberger PK. Survival and growth of juvenile Pacific lampreys tagged with passive integrated transponders (PIT) in freshwater and seawater. Trans Am Fish Soc. 2012;141(5):12601268.

    Article  Google Scholar 

  47. Miehls SM, Holbrook CM, Marsden JE. Diel activity of newly metamorphosed juvenile sea lamprey (Petromyzon marinus). PLoS ONE. 2019;14(2):e0211687.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Moser ML, Jackson AD, Lucas MC, Mueller RP. Behavior and potential threats to survival of migrating lamprey ammocoetes and macrophthalmia. Rev Fish Biol Fishe. 2015;25(1):103–16.

    Article  Google Scholar 

  49. Moser ML, Jackson AD, Mueller RP, Maine AN, Davisson M. Effects of passive integrated transponder (PIT) implantation on Pacific lamprey ammocoetes. Anim Biotelemetry. 2017;5(1):1

    Article  Google Scholar 

  50. Mueller RP, Moursund RA, Bleich MD. Tagging Juvenile Pacific lamprey with passive integrated transponders: methodology, short-term mortality, and influence on swimming performance. North Am J Fish Manag. 2006;26(2):361–6.

    Article  Google Scholar 

  51. Mueller R, Liss S, Deng ZD. Implantation of a New Micro Acoustic Tag in Juvenile Pacific Lamprey and American Eel. JoVE. 2019;(145):e59274.

    Google Scholar 

  52. Nelson JA, Gotwalt PS, Reidy SP, Webber DM. Beyond Ucrit: matching swimming performance tests to the physiological ecology of the animal, including a new fish ‘drag strip.’ Comp Biochem Physiol Part A: Mol Integr Physiol. 2002;133(2):289–302.

    Article  CAS  Google Scholar 

  53. Nightingale J, Stebbing P, Taylor N, McCabe G, Jones G. The long-term effects and detection ranges of passive integrated transponders in white-clawed crayfish Austropotamobius pallipes. Knowl Manag Aquat Ecosyst. 2018;419:20.

    Article  Google Scholar 

  54. Potter IC, Huggins RJ. Observations on the morphology, behaviour and salinity tolerance of downstream migrating River lampreys (Lampetra fluviatitis). J Zool. 1973;169(3):365–79.

    Article  Google Scholar 

  55. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, 2021.

  56. Rechisky EL, Welch DW, Porter AD, Jacobs-Scott MC, Winchell PM. Influence of multiple dam passage on survival of juvenile Chinook salmon in the Columbia River estuary and coastal ocean. PNAS. 2013;110(17):6883–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Ricker WE. Computation and interpretation of biological statistics of fish populations. Bulletin of the Fisheries Research Board of Canada, Bulletin 191, Ottawa, 1975.

  58. Riley WD, Ibbotson AT, Beaumont WRC, Pawson MG, Cook AC, Davison PI. Predation of the juvenile stages of diadromous fish by sea bass (Dicentrarchus labrax) in the tidal reaches of an English chalk stream. Aquat Conserv Mar Freshwat Ecosyst. 2011;21(3):307–12.

    Article  Google Scholar 

  59. Robinson JM, Wilberg MJ, Adams JV, Jones ML. A spatial age-structured model for describing sea lamprey (Petromyzon marinus) population dynamics. Can J Fish Aquat Sci. 2013;70(12):1709.

    Article  Google Scholar 

  60. Sabal MC, Boyce MS, Charpentier CL, Furey NB, Luhring TM, Martin HW, et al. Predation landscapes influence migratory prey ecology and evolution. Trends Ecol Evol. 2021;36(8):737–49.

    Article  PubMed  Google Scholar 

  61. Schreck CB, Fitzpatrick MS, Lerner DT. Determination of passage of juvenile lamprey: Development of tagging protocol. Annual Report, U.S. Army Corps of Engineers BPS-P-00-15b, 1999.

  62. Schultz LD, Heck MP, Kowalski BM, Eagles-Smith CA, Coates K, Dunham JB. Bioenergetics Models to Estimate Numbers of Larval Lampreys Consumed by Smallmouth Bass in Elk Creek, Oregon. North Am J Fish Manag. 2017;37(4):714–23.

    Article  Google Scholar 

  63. Simard LG, Sotola VA, Marsden JE, Miehls S. Assessment of PIT tag retention and post-tagging survival in metamorphosing juvenile sea lamprey. Anim Biotelemetry. 2017;5(1):18.

    Article  Google Scholar 

  64. Sneddon LU. Clinical anesthesia and analgesia in fish. J Exot Pet Med. 2012;21(1):32–43.

    Article  Google Scholar 

  65. Sogard SM. Size-selective mortality in the juvenile stage of teleost fishes: a review. Bull Mar Sci. 1997;60:1129–57.

    Google Scholar 

  66. Summerfelt RC, Smith LS. Anaesthesia, surgery and related techniques. In: Schreck CB, Moyle PB, editors. Methods for fish biology. Bethesda MD: American Fisheries Society; 1990. p. 213–272.

  67. Sutphin ZA, Hueth CD. Swimming performance of Larval Pacific Lamprey (Lampetra tridentata). Northwest Sci. 2010;84(1):196–200.

    Article  Google Scholar 

  68. The IUCN Red List of Threatened Species [Internet]. IUCN Red List of Threatened Species. https://www.iucnredlist.org/en. Accessed Mar 17 2021.

  69. Therneau T, Crowson C. Using time dependent covariates and time dependent coefficients in the cox model. Survival Package. 2013.

  70. Therneau T. A Package for Survival Analysis in R. R package version 2.44–11. 2020.

  71. Johnson B. Study Protocol For A Compassionate Aquaculture Investigational New Animal Drug (INAD) Exemption For AQUI-S®20E (eugenol) (INAD #11-741). 2013. https://www.fws.gov/media/study-protocol-compassionate-aquaculture-investigational-new-animal-drug-inad-exemption-aqui. Accessed 9 July 2022

  72. Wagner GN, Stevens ED, Byrne P. Effects of suture type and patterns on surgical wound healing in rainbow trout. Trans Am Fish Soc. 2000;129(5):1196–205.

    Article  Google Scholar 

  73. Ward TD, Algera DA, Gallagher AJ, Hawkins E, Horodysky A, Jørgensen C, Killen SS, McKenzie DJ, Metcalfe JD, Peck MA, Vu M, Cooke SJ. Understanding the individual to implement the ecosystem approach to fisheries management. Conserv Physiol. 2016;4(1):cow005. https://doi.org/10.1093/conphys/cow005

    Article  PubMed  PubMed Central  Google Scholar 

  74. Wardle CS. Effects of temperature on the maximum swimming speed of fishes. In: Ali MA, editor. Environmental physiology of fishes. Boston: Springer; 1980. p. 519–31.

    Chapter  Google Scholar 

  75. Winger PD, He P, Walsh SJ. Factors affecting the swimming endurance and catchability of Atlantic cod (Gadus morhua). Can J Fish Aquat Sci. 2000. https://doi.org/10.1139/f00-049.

    Article  Google Scholar 

  76. Youson JH, Potter IC. A description of the stages in the metamorphosis of the anadromous sea lamprey, Petromyzon marinus L. Can J Zool. 1979;57(9):1808–17.

    Article  Google Scholar 

  77. Zahl IH, Kiessling A, Samuelsen OB, Hansen MK. Anaesthesia of Atlantic cod (Gadus morhua) — Effect of pre-anaesthetic sedation, and importance of body weight, temperature and stress. Aquaculture. 2009;295(1):52–9.

    Article  CAS  Google Scholar 

  78. Zhang Z, Reinikainen J, Adeleke KA, Pieterse ME, Groothuis-Oudshoorn CGM. Time-varying covariates and coefficients in Cox regression models. Ann Transl Med. 2018;6(7):11–11.

    Article  Google Scholar 

Download references

Acknowledgements

The authors would like to acknowledge several people and entities that were key to the completion of this project. Those include Hammond Bay Biological Station and it’s staff who implemented COVID-19 protocols that allowed us to complete this experiment in a safe manor at the onset of the pandemic of 2020, the Great Lakes Fisheries Commission for funding, staff at Hammond Bay Biological station for assistance during burst swim tests and surgeries, Mary Moser for review and advice on swim tests, and the reviewers at Animal Biotelemetry. Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

Funding

This work was funded by the Great Lakes Fisheries Commission.

Author information

Authors and Affiliations

Authors

Contributions

TFH and CMW performed literature review of previous experiments similar in nature and determined appropriate methodology to meet the desired specific outputs. TFH performed all experiments and was the primary writer of the publication. CMW was the primary reviewer and secondary contributor of the publication. TCS was chiefly responsible for the survival analysis and assisted in the interpretation of endurance swimming results. SMM was responsible for aiding and setup of the holding design and guidance of swim tunnel use. ZDD fabricated the tags used in the experiment and TB was responsible for day-to-day mortality and tag shed checks and the latter wound assessments. All authors contributed to the review of this manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Taylor F. Haas.

Ethics declarations

Ethics approval and consent to participate

This study was performed under Institutional Animal Care and Use Committee (IACUC) ID: PROTO201900108, under the New Protocol Review and these regulations: USDA regulations (9 CFR Parts 1, 2, & 3), the Public Health Service Policy on Humane Care and Use of Laboratory Animals (PHS Policy), the Guide for the Care and Use of Laboratory Animals, 8th Edition (the Guide), and the Guide for the Care and Use of Agricultural Animals in Research and Teaching, 3rd Edition (Ag Guide).

Consent for publication

All the authors have consented to the publication of this manuscript as written, in the scientific journal Animal Biotelemetry. John Hume has expressed written consent for the use of the length frequency data of historical catch in Fig. 1.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Haas, T.F., Castro-Santos, T., Miehls, S.M. et al. Survival, healing, and swim performance of juvenile migratory sea lamprey (Petromyzon marinus) implanted with a new acoustic microtransmitter designed for small eel-like fishes. Anim Biotelemetry 11, 9 (2023). https://doi.org/10.1186/s40317-023-00318-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40317-023-00318-1

Keywords