Research Article |
Corresponding author: Emma M. DeRoy ( emma.m.deroy@gmail.com ) Academic editor: Jonathan Jeschke
© 2022 Emma M. DeRoy, Steven Crookes, Kyle Matheson, Ryan Scott, Cynthia H. McKenzie, Mhairi E. Alexander, Jaimie T. A. Dick, Hugh J. MacIsaac.
This is an open access article distributed under the terms of the CC0 Public Domain Dedication.
Citation:
DeRoy EM, Crookes S, Matheson K, Scott R, McKenzie CH, Alexander ME, Dick JTA, MacIsaac HJ (2022) Predatory ability and abundance forecast the ecological impacts of two aquatic invasive species. NeoBiota 71: 91-112. https://doi.org/10.3897/neobiota.71.75711
|
Characterising interspecific interaction strengths, combined with population abundances of prey and their novel predators, is critical to develop predictive invasion ecology. This is especially true of aquatic invasive species, which can pose a significant threat to the structure and stability of the ecosystems to which they are introduced. Here, we investigated consumer-resource dynamics of two globally-established aquatic invasive species, European green crab (Carcinus maenas) and brown trout (Salmo trutta). We explored the mediating effect of prey density on predatory impact in these invaders relative to functionally analogous native rock crab (Cancer irroratus) and Atlantic salmon (Salmo salar), respectively, feeding on shared prey (Mytilus sp. and Tenebrio molitor, respectively). We subsequently combined feeding rates with each predator’s regional abundance to forecast relative ecological impacts. All predators demonstrated potentially destabilising Type II functional responses towards prey, with native rock crab and invasive brown trout exhibiting greater per capita impacts relative to their trophic analogues. Functional Response Ratios (attack rates divided by handling times) were higher for both invasive species, reflecting greater overall per capita effects compared to natives. Impact projections that incorporated predator abundances with per capita effects predicted severe impacts by European green crabs. However, brown trout, despite possessing higher per capita effects than Atlantic salmon, are projected to have low impact owing to currently low abundances in the sampled watershed. Should brown trout density increase sixfold, we predict it would exert higher impact than Atlantic salmon. Such impact-forecasting metrics and methods are thus vital tools to assist in the determination of current and future adverse impacts associated with aquatic invasive species.
Aquatic invasive species, consumption rate, feeding, freshwater, functional response, Functional Response Ratio, impact, invasion, marine, predation, Relative Impact Potential
Invasive species exert measurable and often catastrophic changes in recipient communities (
Analysis of a predator’s density-dependent consumption rates [i.e. its functional response (FR)] can provide insights into its per capita effect (
While species’ resource consumption can provide insights into their projected ecological impacts (
The objective of our study was to discern whether per capita and overall impacts differed between aquatic invasive species and respective native analogues, using two globally-established invasive species. We utilised the aforementioned trio of metrics (i.e. FR, FRR and RIP) to quantify the predatory impacts of two aquatic invasive species – the marine European green crab (Carcinus maenas) (hereafter, green crab) and the freshwater brown trout (Salmo trutta) – each of which are established in Canada and other regions globally. Both are listed amongst the 100 of the worst invasive species (
To accurately direct management efforts of invasive species, researchers must understand their projected effect across and within regions to which the species has spread, relative to native analogues. Such predictions provide essential information to possible management interventions of invasive species.
Invasive green crab (Carcinus maenas) and native rock crab (Cancer irroratus) (N = 30 each) were collected during the summer of 2015 using Fukui traps (baited with herring) from the upper subtidal zone at North Harbour within Placentia Bay, Newfoundland (NL). Green crab was first detected in this region in 2007 (
Only male crabs with all appendages intact were selected to avoid potential variation in foraging that could result from morphological or behavioural differences between the sexes (
Mytilus sp. mussel prey (25 ± 3 mm) – on which both crab species are known to feed (
Crabs and mussels were transported in containers with seawater to the Northwest Atlantic Fisheries Centre in St. John’s, NL. Species were held separately in holding tanks (275 l) equipped with a flow-through seawater system (11.8 ± 1.5 °C) and fed ad libitum mussels and scallops. The photoperiod (13 h light:11 h dark) was kept constant throughout the experiment. Crabs and mussels were allowed to acclimatise to the system and monitored at least one week prior to and post use in FR trials.
Rock crabs were significantly larger [carapace width (notch to notch) ± SE: rock crab: 104.3 ± 1.57 mm; green crab: 64.2 ± 0.62 mm] and heavier (mass: rock crab: 201.4 ± 7.55 g; green crab: 82.8 ± 2.45 g) than green crab (Wilcoxon rank sum: W = 0, P < 0.0001). Cheliped size, which can be a proxy for crushing strength, for the rock crab was also larger (22.9 ± 0.36 mm) than the green crab crusher cheliped (19.1 ± 0.37 mm) (Wilcoxon rank sum: W = 81, P < 0.0001). This difference resulted from the intentional selection of typical full-sized adult rock and green crabs found in the same habitats, which further allowed comparisons with other studies that used the same approach. Use of both invasive and native adult crabs, therefore, permitted us to discern maximum potential impact of these species.
Experimental trials with invasive brown trout (Salmo trutta) (N = 31: mean ± SE wet weight: 49.4 ± 2.1 g) and native Atlantic salmon (Salmo salar) (N = 18: 91.7 ± 5.4 g) were conducted at the University of Windsor’s Freshwater Restoration Ecology Centre (FREC, LaSalle, ON Canada). Brown trout were purchased from Kolapore Springs Fish Hatchery (Thornbury, ON, Canada) in the summer of 2015 and transported to FREC in insulated tanks with continuously aerated water. Atlantic salmon were reared at FREC.
We selected brown trout as our focal invader given its cosmopolitan distribution and long invasion history (
All fish were acclimatised for one week during which time they were fed mealworms (Tenebrio molitor) ad libitum. Animals were housed in climate-controlled facilities prior to and during experiments (15–17 °C air temperature; 10 h light:14 h dark regime). Fish from individual species were held communally in recirculating housing tanks (800 l; 5% turnover per day), in accordance with University of Windsor’s Animal Care guidelines.
FR trials were run across six circular opaque fibre-glass tanks (275 l; 100 cm diameter and ~ 50 cm water depth) configured in rows of two. Each tank was set up with its own individual light source and inflow to standardise environmental conditions (10.25 °C ± 0.04; ~ 5–10 l/minute flow rate). All tanks were covered with mesh (1.3 cm opening) to prevent potential escape. A random number generator allotted predators and prey density treatments to individual trial tanks.
Trials were conducted between 7am and 3pm. Individual crabs were selected at random and held in experimental tanks supplied with flow-through seawater 48 hours prior to experimental trial to acclimatise and standardise hunger. To initiate a trial, mussels (free of epibionts) were presented haphazardly throughout the tank at six densities (2, 4, 8, 16, 32 and 64 mussels per tank). Each feeding trial lasted five hours, after which we examined prey capture, defined as any crab-mussel interaction that resulted in the crushing or opening the shells of a mussel. We conducted five replicates at each prey density and one control trial for each prey density in the absence of a predator to quantify background mortality rates. Each crab was only used once. We excluded any trial in which the foraging crab moulted in the week following the experiment to further minimise potential variation in crab behaviour during the feeding trial.
Fisheries and Oceans Canada provided regional abundance estimates (CPUE ± SE) for both green and rock crab in North Harbour, Placentia Bay, NL. An average multi-year estimate (2015–2019) was used to account for spatiotemporal variability in population densities. Each yearly estimate was based on 12 traps (four lines of three traps set perpendicular to the shore in the shallow subtidal) set during each of five monthly surveys (June through to October). Trapping estimates recorded the number of crabs obtained per trap per day. The soak time during each deployment was approximately 24 hours with traps set at low tide.
Fish were starved for 24 hours to standardise hunger levels and acclimatised to experimental tanks prior to trial onset. Fish were randomly selected and assigned to one of two flow-through 50 l trial tanks (Mean ± SE: 10.24 °C ± 0.16, flow rate: 1 l/ minute) containing aquarium water. Species were alternated between trials. Tanks were wrapped in black plastic to mitigate observer influences.
To initiate the start of a trial, mealworms (1 cm, cut using a razor) were introduced to the water surface at one of six prey densities (8, 16, 32, 64, 128 and 175 prey per tank). Due to limited stock of Atlantic salmon, three repetitions were conducted per prey density with no re-use. Five replicates were performed per density for brown trout, with the exception of the prey density of 175, for which six replicates were conducted. Mealworms were launched via a weigh boat from the same point across trials. In this regard, they mimicked drifting invertebrates on which salmonids commonly feed (
To compute an estimate of the relative impact for our study species, we procured abundance estimates for Atlantic salmon and brown trout within the Credit River watershed (2015–2019) from Credit Valley Conservation Authority. The Credit River watershed is an important system for juvenile salmonids, including Atlantic salmon, of which both naturally and hatchery-reared individuals are present. Atlantic salmon is native to this region and is currently the subject of restoration efforts (
Abundance data were obtained using single pass electrofishing [see Credit Valley Conservation Authority (2019) for a detailed overview of their methodology]. Abundance estimates were procured in the summer, several months after stocking. Estimates were calculated as the number of individual fish divided by stream area (m2).
Data analyses were performed in R, version 4.0.2 (
We tested for effects of species (factor, two levels), prey density (factor, six levels) and their interaction on consumption rate (continuous) using a GLM (glmmTMB,
We assessed differences in per capita feeding rates via FR curves. We fitted both Type II and III FR models to consumption rate data, using maximum likelihood estimation (bbmle, Bolker and R Development Core Team 2016) and compared fit via AIC. To account for prey depletion over the duration of the experiment, we modelled the resultant Type II FRs using Rogers’ random predator equation (
Ne = N0 (1–exp (a (Neh–T))) (Eqn. 1)
where Ne is the number of prey eaten, N0 is the initial density of prey, a is attack constant, h is handling time and T is the total experimental period (5 hours).
Predator consumption rates – as well as consumer-resource interaction variables, such as search rate, detection distance and handling time on which such rates depend – often vary with individual mass (
a = a0m0.65 (Eqn. 2)
and handling time scaled negatively with predator body mass:
h = h0m-0.65 (Eqn. 3)
In both Eqns. 2 and 3, a0 and h0 are constants and m is predator mass (g).
We fitted the allometrically-scaled FR models using all data for a given species to obtain initial parameter estimates for bootstrapping. We then bootstrapped (N = 100) the data to construct 95% confidence intervals around the fitted curves and extract median values for model parameters. Convergence in FR confidence intervals indicated a lack of significant difference between species’ consumption rates.
We computed FRRs (attack rate a divided by handling time h, i.e. a/h) for each species using median attack rate and handling time parameters. The FRR is a novel metric that has successfully differentiated ecologically-damaging invasive species (
Finally, we determined the maximum feeding rate of each predator (1/h) and combined these values with field abundance estimates to derive a Relative Impact Potential (RIP) estimate according to
This allowed us to discern the relative impact of introduced green crab to native rock crab.
Using the same methodology as described above, differences in overall prey consumption amongst species (factor, two levels), prey density (factor, six levels) and wet weight (continuous) were assessed using a GLM with negative binomial error distribution (glmmTMB,
FR type was confirmed following the protocol outlined above. We subsequently incorporated allometric functions in FR models to account for size discrepancies between salmonids (
Allometrically-scaled FR models were fitted following the aforementioned methodology to obtain median estimates of attack rate and handling time. We then computed FRRs (a/h) for each species as well as corresponding maximum feeding rates (1/h). We subsequently used both species’ maximum feeding rates and abundances to compute the RIP estimate. Stocking effort and abundance were both greater for the native Atlantic salmon. Given field abundance disparities between our focal salmonid predators, we projected impact potential of brown trout in increments of 0.01 ind/m2 to determine the point at which RIP would exceed a value of 1. That is, we determined when brown trout’s projected ecological impact may exceed that of native Atlantic salmon.
In control trials, we experienced no prey mortality and thus ascribed all prey death to predation, which was also directly observed. Predator consumption rates were best described by Type II FRs (Fig.
Functional responses of invasive green crab and native rock crab towards mussel prey. Lines represent initial functional response fits from the random predator equation; shaded areas are 95% confidence intervals (n = 100 non-parametric bootstraps).
On average, rock crabs consumed more mussels than green crabs, though the difference was not significant (Wilcoxon: W = 415, P = 0.61). Analysis of species' per capita effects revealed more nuanced differences in consumptive impact. Rock crab consumed more mussels than green crab, both with and without correcting for size differences between crab species (Table
Relative Impact Potential (RIP) and Functional Response Ratio (FRR) scores, as well as mean ± standard error (SE) estimates of maximum feeding rate, recorded for both invasive-native species pairs. RIP > 1 are predicted to be high impact invaders, those < 1 are low impact relative to native predators.
System | Species | Maximum feeding rate (1/h) (± SE) | RIP | FRR (a/h) (± SE) |
Marine | Green crab | 0.82 (0.01) | 71 | 0.17 (0.01) |
Rock crab | 2.02 (0.19) | 0.12 (0.01) | ||
Freshwater | Brown trout | 0.48 (0.01) | 0.20 | 0.04 (0.002) |
Atlantic salmon | 0.41 (0.02) | 0.001 (<0.001) |
Green crab had a higher FRR, reflecting a steeper FR curve at low prey densities (i.e. larger attack rate, a) that compensated for a higher handling time (h; and, hence, lower maximum feeding rate), indicating European green crab will potentially impact prey populations more than rock crab. Further, green crab abundance (mean ± SE: 29.44 ± 6.91) was orders of magnitude greater than that of rock crab (0.17 ± 0.12), driving a large RIP value (Table
RIP biplot comparing invasive green crab and rock crab feeding upon native mussel prey. Biplots generated using mean ± standard error (SE) estimates for FRs (allometrically-scaled maximum feeding rate, prey/5 hour) and field abundances (CPUE). Ecological impact increases from bottom left to top right.
These results corroborate two independent Ecological Impact Scores used by
Both salmonids exhibited Type II FRs (Fig.
Functional responses of invasive brown trout and native Atlantic salmon towards dried mealworm prey. Lines represent initial functional response fits from the random predator equation; shaded areas are 95% confidence intervals (n = 100 non-parametric bootstraps).
Consumption rates increased significantly with increasing prey density (χ2 = 32.40, df = 5, P < 0.0001) and by predator mass (χ2 = 16.60, df = 1, P < 0.0001). Brown trout was more voracious than Atlantic salmon across all levels of prey availability (χ2 = 46.17, df = 1, P < 0.0001) (Fig.
Brown trout exhibited a higher maximum feeding rate and FRR relative to Atlantic salmon (Table
RIP biplot comparing invasive brown trout and Atlantic salmon feeding upon mealworm prey. Biplots generated using mean ± standard error (SE) estimates for FRs (allometrically-scaled maximum feeding rate, prey/hour) and field abundances (ind/m2). Ecological impact increases from bottom left to top right.
Understanding differences in resource consumption by invasive and native species can provide meaningful insights into potential impacts of invaders in colonised ranges (
Our study highlights strong density-dependence of both per capita and total estimated population effects. All species demonstrated inverse density-dependent prey mortality and potentially destabilising Type II FRs for prey populations. While invasive species often exhibit higher FR curves relative to functionally analogous native species (
Invasive species’ impacts are often context-dependent, in part mediated by abundance (but also per capita differences; see
The high attack rate of green crab at low prey densities can potentially drive mutual prey species to become increasingly rare or even extinct (
Analysis of freshwater salmonids revealed greater levels of consumption by brown trout across all levels of prey availability, despite their smaller size. These findings are consistent with FRs of other high impact invasive species (
Field abundance provides an estimate or proxy of predator numerical response (
Abundance discrepancies may dampen the potential for interspecific competition and produce limited, but strong interactions between the two species, as evidenced by their overlapping isotopic niches (
The current risk of brown trout appears low, based on analyses of population-level impact. However, numerical estimates of abundance, as reported herein, may not represent abundances of the focal species in other systems. It is possible that brown trout impact could be far more substantial in areas where the numerical difference between these species is lower (and in systems without external manipulation), as corroborated by our model (Fig.
Our findings have important implications. If left unchecked, they suggest that burgeoning brown trout populations are likely to produce significant ecological impacts, potentially to the detriment of both native Atlantic salmon and prey populations. Despite the invader’s high per capita effect, management interventions can suppress its potential population-level impact on recipient systems by keeping relative densities low. Sustaining native species’ populations while ensuring productive fisheries – like that of brown trout – therefore depends on balanced management (
Future research should consider the ecological impact of our focal species across a wider variety of prey types in field and laboratory settings. A growing body of literature reinforces the resource-specific nature of invasive species’ per capita effects (e.g.
Additionally, investigation into non-consumptive effects of both green crab and brown trout are needed to ascertain implications for native predators and ultimately consequences to prey. For example, habitat use and depth distribution overlap in shallow waters in areas where our focal crab species co-occur (
Understanding the synergistic influence of co-occurring stressors on invasive species’ impacts is a priority area for invasion science (
Functional and numerical response methodology provides meaningful insights into assessing invader impact and has become especially robust when used in conjunction with the FRR and RIP metrics. These results imply that, if the per capita impacts and relative abundance of non-indigenous species are well-known, its potential relative impact can be predicted and appropriate management actions devised, if needed. Our findings further underscore the importance of population suppression to effectively manage invasive species and promote co-existence with native analogues and prey populations. While our results provide novel insights into the implications of our focal predators, further work is required that incorporates environmental change scenarios.
This work was supported by a Natural Sciences and Engineering Research Council of Canada Discovery Grant and a Canada Research Chair in Aquatic Invasive Species to H.J.M and the Fisheries and Oceans Canada Aquatic Invasive Species Science Program (NL Region) to C.H.M. Experiments were conducted and animals were handled in accordance with the CCAC and the University of Windsor’s Animal Care Committee. The authors thank Dr. Trevor Pitcher for use of FREC facilities. Rock and green crab were procured under a DFO scientific collection licence. The authors thank the following individuals for their assistance in crab collections and laboratory experiments: Terri Wells, Vanessa Reid and Haley Lambert. The authors thank Drs. James Dickey and Ross Cuthbert for their modelling assistance. Abundance data for Lake Ontario were used with permission of Credit Valley Conservation Authority (2020).
Table S1
Data type: docx file
Explanation note: Unscaled mean ± standard error (SE) estimates of maximum feeding rate and Functional Response Ratio (FRR) scores, recorded for both invasive-native species pairs.