ABSTRACT

Objective

The use of captive breeding programs for the conservation of freshwater fishes, be it for research or reintroduction, is becoming more common. Interspecific differences in rearing ecology make applying common rearing techniques difficult for threatened species. This becomes increasingly difficult for species in remote areas, as accessing fish for hatchery propagation can take longer. Bull Trout Salvelinus confluentus has recently been listed as threatened in Alberta and is garnering greater consideration for artificial propagation. Currently, there are no active hatchery programs for Bull Trout, and historically, very few attempts have occurred with limited success. The goal of this research was to establish a laboratory population of Bull Trout and compare the effectiveness of delayed and same day fertilization when transporting gametes or embryos over large distances.

Methods

Here, eggs were either fertilized in the field prior to transport (field) or 1 d after transport (green). Both treatments were collected from Smith-Dorrien Creek, Alberta, and transported to Winnipeg, Manitoba, by airplane (approximately 1,300 km). Differences between treatments were tested using a binomial generalized linear model.

Results

Fertilization rates were similar and high (>99%) for both treatments, and our model did not detect a significant treatment effect. Although a significant treatment effect was not observed, the green treatment had 10% lower survival over the duration of the experiment. The greatest mortality rates were observed at the swim-up stage in both treatments.

Conclusions

Both methods had similar fertilization rates and high survival to the fry stage (green = 36% and field = 46%). The results suggest that both methods are viable options when transporting Bull Trout gametes or embryos long distances for species recovery and research purposes.

Lay Summary

Use of wild source fish for restoration and research can require transport of fish or gametes long distances. Assessment of methods for rearing and transport indicated that survival was similar and high for both fertilized and unfertilized eggs transported ~1,300 km by plane.

INTRODUCTION

There is a long history of fish propagation, initially for consumptive purposes and recreation and more recently for the conservation of fish species and experimentation (Greenberg, 2010). Research studying the rearing of at-risk fish species for conservation purposes has become more common, particularly to recover and reestablish native populations (i.e., fish that were not introduced by humans or allowed to invade due to anthropogenic activities) that have declined due to various anthropogenic effects on their habitats (Halverson, 2008; Jackson et al., 2004). The creation and use of laboratory populations of fish that are propagated from native sources rather than domestic hatchery strains gives researchers the opportunity to better understand the ecology and physiology of the species for conservation purposes, as differences are frequently observed between domestic and wild populations (Scott et al., 2015). When selecting source populations for establishing native broodstocks in hatcheries for reintroduction, it is important to select sources that are similar (genetics and life history) to past populations (Dunham et al., 2011). In some circumstances, source populations for experimentation or restocking may be located in areas that are logistically challenging to access and/or have highly variable spawning times that do not allow gametes to be crossed among populations on the same date.

Salmonidae is a family of fishes that has a long history of propagation for food, recreation, and conservation purposes (Greenberg, 2010; Halverson, 2010). Much of this work has been focused on a few commercially and recreationally important species: most notably Atlantic Salmon Salmo salar and Rainbow Trout Oncorhynchus mykiss, resulting in their establishment in various coldwater habitats throughout the world (Greenberg, 2010; Halverson, 2010). An assessment of stocking rates suggested that Rainbow Trout represent the most dominant fish species among those that have been stocked by biomass in the United States between 1931 and 2004 (Halverson, 2008). Unlike Rainbow Trout, Bull Trout Salvelinus confluentus, which is a culturally and ecologically important salmonid, has received limited attention for artificial rearing (not present in a review of U.S. stocking from 1931 to 2004; Halverson, 2008). To date, few documented attempts have occurred to rear Bull Trout in captivity (Hayes & Banish, 2017; Kissinger et al., 2024) and even fewer studies have been published on the practices of egg-take and early rearing of Bull Trout (Fredenberg et al., 1995). Bull Trout are not typically considered for captive breeding for recreational or commercial purposes (Hayes & Banish, 2017; Kissinger et al., 2024), probably because they spawn in the fall in cold headwater streams that can be difficult to access and require a longer incubation period than spring spawners (Nelson & Partz, 1992). Bull Trout exhibit slow growth because they occupy cold low-productivity habitats (Selong et al., 2001), they are more vulnerable to stress from high temperatures (Selong et al., 2001) and lower oxygen than other salmonids, and larger juveniles and adults are territorial and piscivorous (Nelson & Partz, 1992), resulting in potential harassment and cannibalism when they are held in captivity. Additionally, within parts of their natural distribution, Bull Trout have historically been considered a less-desirable fish and in some cases have suffered intentional eradication measures (Alberta Sustainable Resource Development, 2012), which may further explain their absence as a commonly propagated fish species.

Although Bull Trout are not an ideal candidate for propagation in commercial or recreational hatchery facilities, their listing as a species at risk (Canada Gazette, 2019; U.S. Fish & Wildlife Service, 1999) in parts of their southern distribution in the USA and Canada has heightened the need for laboratory research and the potential for captive breeding (Alberta Sustainable Resource Development, 2012; Dunham et al., 2011). Obtaining gametes from Bull Trout often requires transport over large distances and long periods. Typical fish husbandry practices encourage the fertilization of eggs immediately or soon after removal, but in cases where access to hatchery facilities is limited, fertilized eggs or gametes can be transported to rearing facilities (Bobe & Labbé, 2009). Although the success of fertilization declines the longer gametes are held, no change in fertilization success was observed between 0 and 7 d after the removal of eggs from Rainbow Trout in a study by Carpentier and Billard (1978). Successful fertilization has been documented in Rainbow Trout up to 18 d following the removal of eggs, but fertilization success was very low (near 0, Babiak & Dabrowski, 2003). Although success was observed following long periods of holding at controlled temperatures for Rainbow Trout, differences in fertilization success have been observed among species (Carpentier & Billard, 1978; Jeuthe et al., 2016) and the gametes in most studies that have assessed long-term holding were kept in controlled environments that had minimal exposure to a range of possible transport-related effects (e.g., jostling, low dissolved oxygen, pressure changes, and temperature change; Bobe & Labbé, 2009).

To better understand the physiology of threatened Bull Trout, we set out to establish a laboratory population at the Department of Fisheries and Oceans in Winnipeg, Manitoba, by retrieving gametes from Smith-Dorrien Creek, Alberta (approximately 1,300-km straight line and approximately 20-h transport time). We also sought to determine whether there is a significant difference in survival rate between two propagation methods: (1) transporting field-fertilized hardened eggs (field) and (2) transporting unfertilized eggs and milt separately (green) over large periods and distance followed by fertilization in the laboratory. We assessed survival rate to multiple critical developmental stages, including fertilization (1 d postfertilization [dpf]), 50% eyed, 50% hatch, 50% swim-up, and fry at 342 dpf. We also compared accumulated thermal units (ATU) among treatments to examine potential differences in developmental rates at these developmental stages. Finally, we documented the husbandry practices and timing of the critical developmental stages for both treatments because this information is limited in the literature for Bull Trout. Understanding the process and potential differences in the effectiveness of these two propagation methods will provide staff with a greater suite of options when collecting gametes should additional time be needed for transport or fish from one location become ripe prior to another, where crosses would be beneficial (i.e., broodstock development).

METHODS

Experimental design

All Bull Trout were captured on September 15, 2021, and eggs and milt were harvested from six females and six males from Smith-Dorrien Creek. The gametes were fertilized via two methods (field and green, Figure 1).

Graphical representation of the experimental design. Bull Trout eggs were collected on two separate days to be transported to a lab and fertilized on the same day, resulting in two treatments (green treatment: delayed fertilization of 4,048 eggs postcollection; field treatment: same-day fertilization of 6,111 eggs). Fertilized eggs from each treatment group were divided among incubation trays (three trays for green; two trays for field) At 180 d postfertilization, the alevins were moved to larger tanks to continue growing (alevins from the green group were divided among two tanks and alevins from the field group were divided among three tanks). Created with BioReader.com.
Figure 1.

Graphical representation of the experimental design. Bull Trout eggs were collected on two separate days to be transported to a lab and fertilized on the same day, resulting in two treatments (green treatment: delayed fertilization of 4,048 eggs postcollection; field treatment: same-day fertilization of 6,111 eggs). Fertilized eggs from each treatment group were divided among incubation trays (three trays for green; two trays for field) At 180 d postfertilization, the alevins were moved to larger tanks to continue growing (alevins from the green group were divided among two tanks and alevins from the field group were divided among three tanks). Created with BioReader.com.

For the green treatment, gametes were taken from three females (4,048 eggs) and three males on September 15. Unfertilized eggs were placed in 3.78-L Ziploc bags and approximately 3–5 mL of milt was collected from each male and shipped in 15-mL falcon tubes. The gamete containers were filled with pure O2 and placed in a cooler with ice packs. The gametes were then shipped by airplane in a pressurized compartment (flight time approximately 2 h) and received at Fisheries and Oceans Canada Freshwater Institute (Winnipeg, Manitoba) on September 16 and fertilized that day, resulting in an 18–22-h delay from egg-take to fertilization. Temperature was monitored inside the cooler using hobo temperature loggers (Onset HOBO TidbiT v2 Water Temperature Data Logger, 124 UTBI-001) and gradually increased during shipment from 1.9°C to 6.0°C. Due to concerns with milt from one of the males (i.e., water in the vial), milt from only two males was used to fertilize the eggs from the three females. For both treatments, milt was combined prior to fertilization to increase the chances of fertilization and to randomize the within-treatment influence of parental effects. The eggs and milt were measured and fertilized using a ratio of 200 mL of eggs to 100 µL of milt cocktail.

For the field treatment, gametes were taken from three females (6,111 eggs) and three males on September 16. By taking and fertilizing the eggs on September 16, both treatments (field and green) were fertilized on the same day. The eggs were then water-hardened in water from Smith-Dorrien Creek that had been boiled and cooled to eliminate the chance of disease transfer through contaminated water (Wagner, 2002). The water was saturated to 100% dissolved oxygen before the eggs were water hardened by bubbling them in water from a cylinder of compressed oxygen. The percentage of saturation was monitored using a handheld multiparameter water quality instrument (YSI ProDigital ODO/CT, SKU 627150-1). The fertilized eggs were water-hardened for 2 h prior to movement and were placed in a 22.7-L collapsible container and then in a cooler with ice packs. The cooler was then shipped to Fisheries and Oceans Canada Freshwater Institute by air and received September 16, a duration of 10–11 h after the gametes were collected and fertilized in the field. Temperature was monitored inside the cooler using hobo temperature loggers (Onset HOBO TidbiT v2 Water Temperature Data Logger, UTBI-001) and gradually increased during shipment from 2.9°C to 3.5°C upon arrival at the Freshwater Institute.

Once in the lab, the eggs from both treatments were placed in a Marisource Heath Tray flow-through stacked incubation system with five trays dedicated to eggs (green: trays 1, 2, and 3; field: trays 4 and 5). For the green treatment, a single tray was used for each female to keep the lineages separate, and for the field, two trays were used to hold the embryos from the fertilized treatment, as they were mixed during shipment. Each tray was 82 × 60× 64 cm, and the eggs were placed in the trays between 1 and 2 eggs deep on a trout screen (mesh size 1.4 × 1.8 mm). Fertilized eggs were approximately 5.5 mm in diameter after water hardening. The system was flow through (freshwater inflow rate 5 L/min), but some water (flow rate 3 L/min) was recirculated to increase flow and act as a redundancy should power be lost to the main water source. Initial fertilization was not enumerated until 24 h postfertilization, where staff considered successful fertilization as eggs with the presence of embryo development. Mortalities from each tray were enumerated daily and were identified as opaque eggs. To assess mortalities once hatched, alevins were gently probed to assess their status if they were not actively moving. The timing of different developmental stages (e.g., fertilization, 50% eyed, 50% hatched, and 50% swim-up) were documented. Alevins from each treatment were kept separate, but trays within treatments were pooled together and transferred to larger flow-through tanks (volume 188 L/tank, green: tanks A and B; field: tanks C, D, and E; Figure 1) on March 15, 2022. An extra tank was used for the field treatment to accommodate the higher number of fish, and each tank was held at approximately 1,900 fish/tank. Mortalities from each treatment were monitored until August 23, 2022, when the tanks were combined and treatments were no longer discernible.

Accumulated thermal units

The temperatures for all treatments were increased to mimic natural changes that are observed in the wild. In general, the eggs were fertilized around 4°C and held at this temperature until 121 dpf, when both treatments were 50% hatched, at which time temperature was increased gradually to approximately 10°C over the next 70 d (Table 1; Supplement 1 [see online Supplementary Material]).

Table 1.

Accumulated thermal units (ATUs) calculated for each developmental life stage milestone for both treatment groups. dpf = days postfertilization.

Life stageTreatmentDatedpfATUs (˚C)Tank
50% eyedGreen2021-11-2875276.08Tank A—tray 1
2021-11-2875277.623Tank B—tray 2
2021-11-2774279.548Tank C—tray 3
Field2021-11-2774286.018Tank D—tray 4
2021-11-2673289.28Tank E—tray 5
50% hatchGreen2022-01-13121446.98Tank A—tray 1
2022-01-12120444.294Tank B—tray 2
2022-01-12120453.210Tank C—tray 3
Field2022-01-10118451.159Tank D—tray 4
2022-01-09117458.474Tank E—tray 5
50% emergedGreen2022-04-152131,081.920Tank A—tray 1
2022-04-152131,078.664Tank B—tray 2
2022-04-152131,089.186Tank A—tray 3
Field2022-04-142121,084.698Tank C—tray 4
2022-04-142121,091.722Tank D—tray 4
2022-04-142121,099.766Tank D—tray 5
2022-04-132111,101.043Tank E—tray 5
Life stageTreatmentDatedpfATUs (˚C)Tank
50% eyedGreen2021-11-2875276.08Tank A—tray 1
2021-11-2875277.623Tank B—tray 2
2021-11-2774279.548Tank C—tray 3
Field2021-11-2774286.018Tank D—tray 4
2021-11-2673289.28Tank E—tray 5
50% hatchGreen2022-01-13121446.98Tank A—tray 1
2022-01-12120444.294Tank B—tray 2
2022-01-12120453.210Tank C—tray 3
Field2022-01-10118451.159Tank D—tray 4
2022-01-09117458.474Tank E—tray 5
50% emergedGreen2022-04-152131,081.920Tank A—tray 1
2022-04-152131,078.664Tank B—tray 2
2022-04-152131,089.186Tank A—tray 3
Field2022-04-142121,084.698Tank C—tray 4
2022-04-142121,091.722Tank D—tray 4
2022-04-142121,099.766Tank D—tray 5
2022-04-132111,101.043Tank E—tray 5
Table 1.

Accumulated thermal units (ATUs) calculated for each developmental life stage milestone for both treatment groups. dpf = days postfertilization.

Life stageTreatmentDatedpfATUs (˚C)Tank
50% eyedGreen2021-11-2875276.08Tank A—tray 1
2021-11-2875277.623Tank B—tray 2
2021-11-2774279.548Tank C—tray 3
Field2021-11-2774286.018Tank D—tray 4
2021-11-2673289.28Tank E—tray 5
50% hatchGreen2022-01-13121446.98Tank A—tray 1
2022-01-12120444.294Tank B—tray 2
2022-01-12120453.210Tank C—tray 3
Field2022-01-10118451.159Tank D—tray 4
2022-01-09117458.474Tank E—tray 5
50% emergedGreen2022-04-152131,081.920Tank A—tray 1
2022-04-152131,078.664Tank B—tray 2
2022-04-152131,089.186Tank A—tray 3
Field2022-04-142121,084.698Tank C—tray 4
2022-04-142121,091.722Tank D—tray 4
2022-04-142121,099.766Tank D—tray 5
2022-04-132111,101.043Tank E—tray 5
Life stageTreatmentDatedpfATUs (˚C)Tank
50% eyedGreen2021-11-2875276.08Tank A—tray 1
2021-11-2875277.623Tank B—tray 2
2021-11-2774279.548Tank C—tray 3
Field2021-11-2774286.018Tank D—tray 4
2021-11-2673289.28Tank E—tray 5
50% hatchGreen2022-01-13121446.98Tank A—tray 1
2022-01-12120444.294Tank B—tray 2
2022-01-12120453.210Tank C—tray 3
Field2022-01-10118451.159Tank D—tray 4
2022-01-09117458.474Tank E—tray 5
50% emergedGreen2022-04-152131,081.920Tank A—tray 1
2022-04-152131,078.664Tank B—tray 2
2022-04-152131,089.186Tank A—tray 3
Field2022-04-142121,084.698Tank C—tray 4
2022-04-142121,091.722Tank D—tray 4
2022-04-142121,099.766Tank D—tray 5
2022-04-132111,101.043Tank E—tray 5

Temperature was monitored in the heath trays and tanks using hobo temperature loggers that were set for 10-min measurements throughout the duration of the experiment. The temperature measures were averaged to obtain the daily temperature of each tray. The sum of average daily temperatures from each tray and tank was used to calculate ATU for survival to different developmental stages for each treatment. The number of surviving individuals out of the initial number of inseminated oocytes and ATU were calculated to 342 dpf and compared among treatments.

Statistical analysis

The proportion of individuals surviving at the end of each developmental stage allowed us to explicitly model survival within each treatment group. Due to the small sample size, we were unable to employ more complex modeling approaches. We assumed independence of observations, although we recognized the presence of unaccounted batch effects between developmental stages.

To model the probability, πi, of the proportion of surviving individuals at the end of each discrete developmental stage as a function of the covariates, we applied a binomial generalized linear model with a logit link function:

where “survived” was the number of inseminated oocytes that survived at the end of a developmental stage out of the initial number, ni, at the beginning of the stage. The covariates “treatment” (categorical with two levels: field and green) and developmental “stage” (categorical with five levels [1 dpf, 50% eyed, 50% hatched, 50% swim-up, and fry]) were included as fixed effects. Interactions between these factors were included to obtain treatment- and stage-specific survival estimates. The 95% point-wise confidence intervals were estimated on the linear predictor scale and transformed back to the response scale using the inverse of the link function. The post hoc analysis included pairwise comparisons between treatment groups for each stage using the “emmeans” R package (Lenth, 2024). All statistics were completed in R Statistics (R Core Team, 2023).

RESULTS

Development

The fertilized eggs from the field treatment reached each respective developmental stage slightly earlier (1–3 d) than the eggs from the green treatment (Table 1). Eyed development was first observed at 60 dpf, and less than 3 weeks later, all eggs were eyed. Hatching was first observed in the field treatment 106 dpf, and by 120 dpf, both treatments had achieved 50% hatching success (Table 1). Hatching was completed by 124 dpf in both treatments, and 50% swim-up was achieved by 212 dpf and completed after 234 dpf (Table 1).

Accumulated thermal units

Temperature ranged from 3.5°C to 13.2°C, with an average of 7.6°C throughout the duration of the experiment (Supplement 1). At all developmental stages, the field treatment had higher average ATU (Table 1). However, we observed no statistical differences between the field treatment and green treatment for ATU with 50% eyed (p = 0.098), 50% hatched (p = 0.427), and 50% swim-up (p = 0.321). At the time of 50% swim-up, values for ATU were slightly greater in the field treatment (1,082 ATU green and 1,093 ATU field).

Survival

Both treatments successfully reared fish to the fry stage (Figure 2), but linear trends in mortality significantly differed between each developmental stage, suggesting a cumulative 10% difference among treatments by the end of this experiment (Table 2; Figure 2B). Despite a 10% difference at the end of the experiment, pairwise comparisons of fertilization rates were similar and high in both treatments (>99%; Table 3) and no significant linear treatment effect was observed (p = 0.074; Table 2), suggesting that the transportation and timing of fertilization method did not significantly influence initial fertilization and overall survival (Figure 2). The greatest decrease in survival was observed at the 50% swim-up stage, and the pairwise comparisons indicate that survival of the green treatment was significantly lower at all developmental stages except fertilization and 50% hatched (Figure 2A; Table 3).

The figure displays (A) pairwise comparisons of predicted survival probability at each developmental stage and corresponding 95% point-wise confidence intervals and (B) the proportion of fertilized oocytes surviving to the end of each developmental stage in each treatment group (field: red; green: black; *p < 0.05; **p < 0.01; ***p < 0.001).
Figure 2.

The figure displays (A) pairwise comparisons of predicted survival probability at each developmental stage and corresponding 95% point-wise confidence intervals and (B) the proportion of fertilized oocytes surviving to the end of each developmental stage in each treatment group (field: red; green: black; *p < 0.05; **p < 0.01; ***p < 0.001).

Table 2.

Model coefficient estimates and associated confidence intervals and p-values for the proportion of surviving fertilized oocytes. Bold indicates significance p < 0.05.

PredictorsEstimatesCIp
(Intercept)5.425.06 to 5.82<0.001
Stage [50% eyed]−1.22−1.67 to −0.81<0.001
Stage [50% hatch]−2.44−2.86 to −2.07<0.001
Stage [50% swim-up]−4.87−5.28 to −4.51<0.001
Stage [fry]−3.96−4.37 to −3.59<0.001
Treatment [Green]0.69−0.03 to 1.500.074
Stage [50% eyed] × Treatment [green]−1.86−2.71 to −1.10<0.001
Stage [50% hatch] × Treatment [green]−0.54−1.37 to 0.210.176
Stage [50% swim-up] × Treatment [green]−1.04−1.85 to 0.310.008
Stage [fry] × Treatment [green]−1.06−1.88 to 0.330.007
PredictorsEstimatesCIp
(Intercept)5.425.06 to 5.82<0.001
Stage [50% eyed]−1.22−1.67 to −0.81<0.001
Stage [50% hatch]−2.44−2.86 to −2.07<0.001
Stage [50% swim-up]−4.87−5.28 to −4.51<0.001
Stage [fry]−3.96−4.37 to −3.59<0.001
Treatment [Green]0.69−0.03 to 1.500.074
Stage [50% eyed] × Treatment [green]−1.86−2.71 to −1.10<0.001
Stage [50% hatch] × Treatment [green]−0.54−1.37 to 0.210.176
Stage [50% swim-up] × Treatment [green]−1.04−1.85 to 0.310.008
Stage [fry] × Treatment [green]−1.06−1.88 to 0.330.007
Table 2.

Model coefficient estimates and associated confidence intervals and p-values for the proportion of surviving fertilized oocytes. Bold indicates significance p < 0.05.

PredictorsEstimatesCIp
(Intercept)5.425.06 to 5.82<0.001
Stage [50% eyed]−1.22−1.67 to −0.81<0.001
Stage [50% hatch]−2.44−2.86 to −2.07<0.001
Stage [50% swim-up]−4.87−5.28 to −4.51<0.001
Stage [fry]−3.96−4.37 to −3.59<0.001
Treatment [Green]0.69−0.03 to 1.500.074
Stage [50% eyed] × Treatment [green]−1.86−2.71 to −1.10<0.001
Stage [50% hatch] × Treatment [green]−0.54−1.37 to 0.210.176
Stage [50% swim-up] × Treatment [green]−1.04−1.85 to 0.310.008
Stage [fry] × Treatment [green]−1.06−1.88 to 0.330.007
PredictorsEstimatesCIp
(Intercept)5.425.06 to 5.82<0.001
Stage [50% eyed]−1.22−1.67 to −0.81<0.001
Stage [50% hatch]−2.44−2.86 to −2.07<0.001
Stage [50% swim-up]−4.87−5.28 to −4.51<0.001
Stage [fry]−3.96−4.37 to −3.59<0.001
Treatment [Green]0.69−0.03 to 1.500.074
Stage [50% eyed] × Treatment [green]−1.86−2.71 to −1.10<0.001
Stage [50% hatch] × Treatment [green]−0.54−1.37 to 0.210.176
Stage [50% swim-up] × Treatment [green]−1.04−1.85 to 0.310.008
Stage [fry] × Treatment [green]−1.06−1.88 to 0.330.007
Table 3.

Pairwise comparisons of fertilized oocyte survival for two treatment groups across key developmental stages. dpf = days postfertilization; Inf = infinite. Bold indicates significance p < 0.05.

ContrastStageEstimateSEdfZp
Green/field1 dpf1.9920.768Inf1.7870.074
Green/field50% eyed0.3090.04Inf−9.053<0.001
Green/field50% hatch1.1620.117Inf1.4950.135
Green/field50% swim-up0.7070.03Inf−8.049<0.001
Green/fieldfry0.6920.047Inf−5.425<0.001
ContrastStageEstimateSEdfZp
Green/field1 dpf1.9920.768Inf1.7870.074
Green/field50% eyed0.3090.04Inf−9.053<0.001
Green/field50% hatch1.1620.117Inf1.4950.135
Green/field50% swim-up0.7070.03Inf−8.049<0.001
Green/fieldfry0.6920.047Inf−5.425<0.001
Table 3.

Pairwise comparisons of fertilized oocyte survival for two treatment groups across key developmental stages. dpf = days postfertilization; Inf = infinite. Bold indicates significance p < 0.05.

ContrastStageEstimateSEdfZp
Green/field1 dpf1.9920.768Inf1.7870.074
Green/field50% eyed0.3090.04Inf−9.053<0.001
Green/field50% hatch1.1620.117Inf1.4950.135
Green/field50% swim-up0.7070.03Inf−8.049<0.001
Green/fieldfry0.6920.047Inf−5.425<0.001
ContrastStageEstimateSEdfZp
Green/field1 dpf1.9920.768Inf1.7870.074
Green/field50% eyed0.3090.04Inf−9.053<0.001
Green/field50% hatch1.1620.117Inf1.4950.135
Green/field50% swim-up0.7070.03Inf−8.049<0.001
Green/fieldfry0.6920.047Inf−5.425<0.001

DISCUSSION

As far as we know, this is the first study to show that Bull Trout gametes and embryos can be successfully raised to hatch following long-distance transportation by air and up to 22 h of holding. The data suggest that similar amounts of ATU are required for each developmental stage among treatments but also show that mortality was higher in the green treatment with 10% less survival by the end of the experiment than in the field treatment. Although higher mortality was observed, both treatments successfully reared large numbers of Bull Trout to the fry stage and our model did not detect a significant treatment effect. The high and similar fertilization rates between treatments suggest that both techniques are equally effective.

Development

To date, little is known about the development of Bull Trout in hatchery systems, as there are only a few examples of hatchery and laboratory rearing (Gould, 1987; Hayes & Banish, 2017; Kissinger et al., 2024). The extant information about Bull Trout development suggests that average ATU in our study, for both treatments from fertilization to 50% emergence, was approximately 300 ATU higher than has been described previously (Bowerman et al., 2014; Gould, 1987). Variability in ATU within a species can be large (e.g., 650–1,440 ATU, for Chum Salmon Oncorhynchus keta; Merritt & Raymond, 1983), indicating that some level of plasticity or local adaptation is present within a species. The finding that Bull Trout from Smith-Dorrien Creek require higher ATU than populations that are found at more southern latitudes may be driven by regional differences in climate, where winters are longer and colder relative to southern range extents (Bowerman et al., 2014; Gould, 1987). Spawning and swim-up phenology can vary among populations (Austin et al., 2019; Bowerman et al., 2014), which is thought to be governed by differences in water temperature, and subsequently, ATU, where fish occupying colder streams typically spawn earlier (Austin et al., 2019). When four genetically differentiated strains of Brook Trout were brought into a laboratory and reared at the same temperatures, significant differences in degree-days to hatch were observed among strains, which was believed to be linked to local adaptation (Baird et al., 2002). Based on the geographic differences observed between the Bull Trout population in Smith-Dorrien Creek and other Bull Trout populations referenced prior, it is likely that there are genetic differences, as population differentiation within Bull Trout has been observed at much smaller scales (Warnock et al., 2010). Although they were not tested in this study, these genetic differences may be linked to observed differences in ATU among populations and warrant further investigation.

The presence of high mortality rates in both treatments during swim-up was not surprising because transitioning from endogenous to exogenous nutrient sources can be challenging for fish, especially when they are artificially reared (Jeuthe et al., 2015). The survival to swim-up that we observed in our experiment was lower than survival rates in hatchery systems (50–80%, Essington et al., 2000; Hilborn, 1992). However, this is not surprising, as Arctic Char Salvelinus alpinus, which are close relatives to Bull Trout, have similarly low survival in hatchery systems (Jeuthe et al., 2016), which could be linked to a preference for development at colder temperatures or difficulties transitioning to exogenous food sources (Jeuthe et al., 2015). Although survival to the fry stage in this Bull Trout population was lower than that of many hatchery-strain salmonids, an average survival of 41% for these two methods to 342 dpf suggests that both methods that were employed here represent a viable option for rearing Bull Trout in captivity. Should this become a common approach for species recovery, improvements to the swim-up and fry stage would have the greatest benefits to overall survival.

Managing at-risk species

Demonstrating that Bull Trout populations can be reared in captivity following long-distance transport opens the door for not only laboratory experiments, like those described here, but also the possibility of creating broodstocks for restoration purposes. Transport and holding of gametes for up to 22 h with nearly 100% fertilization success for the green treatment is highly encouraging because it allows users to keep milt and eggs separate if in situ fertilization is not possible. The ability to keep gametes separate until arrival in a laboratory facility and hold gametes for extended periods offers users the flexibility to conduct various parental crosses, even in situations where local populations spawn at different times (e.g., Austin et al., 2019). This can be especially important when attempting to incorporate as much genetic diversity as possible into a recovery composite broodstock from various sources (Whiteley et al., 2015). In addition, as many Bull Trout populations receive higher levels of disruption in areas that are near population centers, many of the high-density stocks that would be suitable donors are found in remote locations that are challenging to access (requiring long hikes and/or helicopter access).

In Alberta, there are four provincial hatcheries and three are within the historic distribution of Bull Trout. All three hatcheries within the historic distribution of Bull Trout are in the southwestern corner of the province and would be approximately 550 linear km from the most northern extent of Bull Trout within the province of Alberta, creating challenges for potential collection and transport of gametes or embryos. Alternatively, Alberta Provincial Fisheries biologists are testing different methods for rearing native salmonids at sites (i.e., Westslope Cutthroat Trout Oncorhynchus lewisi), including remote streamside incubation chambers and portable incubation trailers (B. Meagher, Alberta Environment and Protected Areas, personal communication). However, with fall-spawning species like Bull Trout, winter freezing presents challenges for remote streamside incubation chambers (as are used for Westslope Cutthroat Trout) and portable incubation trailers. Therefore, Alberta Environment and Protected Area staff are currently supporting a research initiative at the University of Calgary to assess the feasibility of egg incubation chambers for Bull Trout that are buried in the sediment (methods similar to those of Baxter & McPhail, 1999). The results from this research indicate high levels of survival to hatch (71%) and that this method could be a viable alternative to hatchery propagation for fall-spawning species in cold climates (Lepine, 2024).

Another risk when transporting fish is disease transfer through transporting live fish and/or water. This is of particular concern with the recent documentation of whirling disease in Alberta (Canadian Food Inspection Agency, 2019). Consequently, the water for the field treatment from this study had to be decontaminated prior to use and transport (Wagner, 2002). Transporting gametes without the need for water decreases the probability of spreading whirling disease because gametes are not susceptible to infection (Markiw, 1991). Showing that both methods had similar fertilization rates suggests that the green method can be used, which reduces the risk of transferring aquatic diseases to fish and between locations. Although risks of disease transfer are reduced using the green method, other factors must be considered prior to moving fish or gametes and practitioners should consider existing guidance when undertaking relocations (e.g., Dunham et al., 2011).

Conclusion

We show that the green and field methods for transporting Bull Trout embryos and gametes over long distances were equally effective at initial fertilization, no overall treatment effect was observed, and both approaches successfully reared large numbers of fish to the fry stage. Improvements to rearing techniques after the swim-up stage would increase overall survival rates, and we suggest that further research and considerations focus on these stages of development. We show that both approaches provide viable options for research and species recovery actions for at-risk Bull Trout populations.

SUPPLEMENTARY MATERIAL

Supplementary material is available at North American Journal of Aquaculture online.

DATA AVAILABILITY

The data that support the findings of this study are available from the corresponding author upon reasonable request.

ETHICS STATEMENT

All fish were treated in accordance with the guidelines of the following permits as required by the Government of Canada to transport and conduct research on aquatic species at risk. Alberta Environment and Parks (Research License #21-0217RL; Research and Collection Permit #21-384), DFO Freshwater Institute Animal Care Committee (Animal Use Protocol #FWI-ACC-2021-66), Canadian Food Inspection Agency (Domestic Movement Permit #D-2021-119874), Government of Manitoba (Live Fish Handling Permit #34097479 LFH-ITP 25-2021), and Species at Risk program (Species at Risk Act Section 73 Permit #21-PCAA-00048).

FUNDING

Funding was provided by Fisheries and Oceans Canada (DFO) through the Program and the Forestry Resource Improvement Association of Alberta (FOOMOD-01-030).

ACKNOWLEDGMENTS

The authors thank T. Rudolfsen, M. Teillet (DFO Edmonton), G. Eisler, A. Robinson, and H. Keyes (AJM Environmental Inc., Calgary) for assistance with planning logistics, fieldwork, and fish capture. John Enns (Calgary) provided valuable knowledge and experience collecting and shipping gametes/fertilized eggs. A. Spiteri, C. Scotland, and J. Symes (Alberta Environment and Parks, Kananaskis District) assisted with field navigation, planning logistics, and safety. L. Peterson and multiple volunteers from Trout Unlimited Canada (Calgary) provided field assistance with redd surveys. K. Wautier and T. Durhack (DFO Winnipeg) received and fertilized gametes in Winnipeg and provided support with the design and preparation of incubation units. M. Panciera and B. Romaniuk (DFO Winnipeg) helped care and raise fish after hatching. Thanks to T. Winter (Alberta Environment and Protected Areas) and R. Bennett (DFO Winnipeg) for help with review.

REFERENCES

Alberta Sustainable Resource Development
. (
2012
).
Bull Trout conservation management plan 2012–17
(Species at Risk Conservation Management Plan No. 8).
Alberta Sustainable Resource Development
.

Austin
,
C. S.
,
Essington
,
T. E.
, &
Quinn
,
T. P.
(
2019
).
Spawning and emergence phenology of Bull Trout Salvelinus confluentus under differing thermal regimes
.
Journal of Fish Biology
,
94
,
191
195
.

Babiak
,
I.
, &
Dabrowski
,
K.
(
2003
).
Refrigeration of Rainbow Trout gametes and embryos
.
Journal of Experimental Zoology Part A: Comparative Experimental Biology
,
300A
,
140
151
.

Baird
,
H. B.
,
Krueger
,
C. C.
, &
Josephson
,
D. C.
(
2002
).
Differences in incubation period and survival of embryos among Brook Trout strains
.
North American Journal of Aquaculture
,
64
,
233
241
.

Baxter
,
J. S.
, &
McPhail
,
J. D.
(
1999
).
The influence of redd site selection, groundwater upwelling, and over-winter incubation temperature on survival of Bull Trout (Salvelinus confluentus) from egg to alevin
.
Canadian Journal of Zoology
,
77
,
1233
1239
.

Bobe
,
J.
, &
Labbé
,
C.
(
2009
). Chilled storage of sperm and eggs. In
E.
 
Cabrita
,
V.
 
Robles
, &
P.
 
Herráez
(Eds.),
Methods in reproductive aquaculture: Marine and freshwater species
(pp.
219
236
).
CRC Press
.

Bowerman
,
T.
,
Neilson
,
B. T.
, &
Budy
,
P.
(
2014
).
Effects of fine sediment, hyporheic flow, and spawning site characteristics on survival and development of Bull Trout embryos
.
Canadian Journal of Fisheries and Aquatic Sciences
,
71
,
1059
1071
.

Canada Gazette
. (
2019
). Exemption order for certain licenses, authorizations and documents (Bull Trout [Saskatchewan–Nelson rivers populations]) (SOR/2019-288, pp. 5824). Canada Gazette, Part II, 153(17). Retrieved July 2024, from https://www.gazette.gc.ca/rp-pr/p2/2019/2019-08-21/pdf/g2-15317.pdf

Canadian Food Inspection Agency
. (
2019
). Confirmed detections of whirling disease–Alberta. Canadian Food Inspection Agency. https://inspection.canada.ca/en/animal-health/aquatic-­animals/diseases/reportable-diseases/whirling-disease/alberta

Carpentier
,
P.
, &
Billard
,
R.
(
1978
).
Conservation à court terme des gamètes de salmonidés à des températures voisines de 0°C
.
Annales de Biologie Animale Biochimie Biophysique
,
18
,
1083
1088
.

Dunham
,
J. B.
,
Gallo
,
K.
,
Shively
,
D.
,
Allen
,
C.
, &
Goehring
,
B.
(
2011
).
Assessing the feasibility of native fish reintroductions: A framework applied to threatened Bull Trout
.
North American Journal of Fisheries Management
,
31
,
106
115
.

Essington
,
T. E.
,
Quinn
,
T. P.
, &
Ewert
,
V. E.
(
2000
).
Intra- and interspecific competition and the reproductive success of sympatric Pacific salmon
.
Canadian Journal of Fisheries and Aquatic Sciences
,
57
,
205
213
.

Fredenberg
,
W.
,
Dwyer
,
P.
, &
Barrows
,
R.
(
1995
).
Experimental Bull Trout hatchery progress report 1993–1994
.
Montana Department of Fish, Wildlife and Parks
.

Gould
,
W. R.
(
1987
).
Features in the early development of Bull Trout (Salvelinus confiuentus)
.
Northwest Science
,
61
,
264
268
.

Greenberg
,
P.
(
2010
).
Four fish: The future of the last wild food
.
Penguin Press
.

Halverson
,
A.
(
2008
).
Stocking trends: A quantitative review of governmental fish stocking in the United States, 1931 to 2004
.
Fisheries
,
33
,
69
75
.

Halverson
,
A.
(
2010
).
An entirely synthetic fish: How Rainbow Trout beguiled America and overran the world
.
Yale University Press
.

Hayes
,
M. F.
, &
Banish
,
N. P.
(
2017
).
Translocation and reintroduction of native fishes: A review of Bull Trout Salvelinus confluentus with applications for future reintroductions
.
Endangered Species Research
,
34
,
191
209
.

Hilborn
,
R.
(
1992
).
Institutional learning and spawning channels for Sockeye Salmon (Oncorhynchus nerka)
.
Canadian Journal of Fisheries and Aquatic Sciences
,
49
,
1126
1136
.

Jackson
,
J. R.
,
Boxrucker
,
J. C.
, &
Willis
,
D. W.
(
2004
). Trends in agency use of propagated fishes as a management tool in inland fisheries. In
M. J.
 
Nickum
,
P. M.
 
Mazik
,
J. G.
 
Nickum
, &
D. D.
 
MacKinlay
(Eds.),
Propagated fish in resource management
(Symposium 44, pp.
121
138
).
American Fisheries Society
.

Jeuthe
,
H.
,
Brännäs
,
E.
, &
Nilsson
,
J.
(
2015
).
Thermal stress in Arctic Charr Salvelinus alpinus broodstock: A 28 year case study
.
Journal of Fish Biology
,
86
,
1139
1152
.

Jeuthe
,
H.
,
Brännäs
,
E.
, &
Nilsson
,
J.
(
2016
).
Effects of variable egg incubation temperatures on the embryonic development in Arctic Charr Salvelinus alpinus
.
Aquaculture Research
,
47
,
3753
3764
.

Kissinger
,
B. C.
,
Sullivan
,
M. G.
,
Paul
,
A. J.
,
Meinke
,
A.
, &
Post
,
J.
(
2024
).
Establishment of Bull Trout in a previously fishless subalpine lake by translocation
.
North American Journal of Fisheries Management
,
44
,
520
531
.

Lenth
,
R
. (
2024
).
emmeans: Estimated marginal means, aka least-squares means
. R package version 1.10.6-090003. https://rvlenth.github.io/emmeans/

Lepine
,
T. M.
(
2024
).
An experimental test of Bull Trout egg to alevin survival in the wild, with implications for conservation translocations in Alberta
[Master's thesis, University of Calgary].
PRISM
. https://hdl.handle.net/1880/119598

Markiw
,
M. E.
(
1991
).
Whirling disease: Earliest susceptible age of Rainbow Trout to the triactinomyxid of Myxobulus cerebralis
.
Aquaculture
,
92
,
1
6
.

Merritt
,
M. F.
, &
Raymond
,
J. A.
(
1983
).
Early life-history of Chum Salmon in the Noatak River and Kotzebue Sound
 
(Report 1)
.
Alaska Department of Fish and Game
.

Nelson
,
J.
, &
Partz
,
M.
(
1992
).
The fishes of Alberta
.
The University of Alberta Press
.

R Core Team
. (
2023
).
R: A language and environment for statistical computing
.
R Foundation for Statistical Computing
.

Scott
,
M. A.
,
Dhillon
,
R. S.
,
Schulte
,
P. M.
, &
Richards
,
J. G.
(
2015
).
Physiology and performance of wild and domestic strains of diploid and triploid Rainbow Trout (Oncorhynchus mykiss) in response to environmental challenges
.
Canadian Journal of Fisheries and Aquatic Sciences
,
72
,
125
134
.

Selong
,
J. H.
,
McMahon
,
T. E.
,
Zale
,
A. V.
, &
Barrows
,
F. T.
(
2001
).
Effect of temperature on growth and survival of Bull Trout, with application of an improved method for determining thermal tolerance in fishes
.
Transactions of the American Fisheries Society
,
130
,
1026
1037
.

U.S. Fish and Wildlife Service
. (
1999
).
Endangered and threatened wildlife and plants; determination of threatened status for Bull Trout in the coterminous United States
.
Federal Register
,
64
,
58910
58933
. https://www.govinfo.gov/content/pkg/FR-1999-11-01/pdf/99-28295.pdf

Wagner
,
E. J.
(
2002
). Whirling disease prevention, control, and management: A review. In
Jerri L.
 
Bartholomew
&
J.
 
Christopher Wilson
(Eds.),
Whirling disease: Reviews and current topics
(Symposium 29, pp.
217
225
).
American Fisheries Society
.

Warnock
,
W. G.
,
Rasmussen
,
J. B.
, &
Taylor
,
E. B.
(
2010
).
Genetic clustering methods reveal Bull Trout (Salvelinus confluentus) fine-scale population structure as a spatially nested hierarchy
.
Conservation Genetics
,
11
,
1421
1433
.

Whiteley
,
A. R.
,
Fitzpatrick
,
S. W.
,
Funk
,
W. C.
, &
Tallmon
,
D. A.
(
2015
).
Genetic rescue to the rescue
.
Trends in Ecology and Evolution
,
30
,
42
49
.

Author notes

CONFLICTS OF INTEREST: The authors declare that they have no conflict of interest.

This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial License (https://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact [email protected] for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact [email protected].

Supplementary data