5urn:lsid:arphahub.com:pub:8D1BC1DD-8175-5933-B147-C839B202D5BAurn:lsid:zoobank.org:pub:C132FA7C-A22D-4474-AAAF-B7A6889D0FF9NeoBiotaNB1619-00331314-2488Pensoft Publishers10.3897/neobiota.45.3311633116Research ArticleBivalviaBiological InvasionsHabitats Ecosystems & Natural SpacesMarine & Freshwater ecologyEuropeRelative impacts of the invasive Pacific oyster, Crassostreagigas, over the native blue mussel, Mytilusedulis, are mediated by flow velocity and food concentrationJoycePatrick W.S.patrick.joyce.254@gmail.comhttps://orcid.org/0000-0003-1058-79011KregtingLouise1DickJaimie T.A.1School of Biological Sciences, Queen’s University Marine Laboratory, 12–13 The Strand, Portaferry, BT22 1PF, UKQueen’s UniversityBelfastUnited KingdomInstitute for Global Food Security, School of Biological Sciences, Queen’s University Belfast, MBC, 97 Lisburn Road, Belfast, BT9 7BL, UKQueen’s UniversityBelfastUnited KingdomSchool of Natural and Built Environment, Queen’s University Marine Laboratory, 12–13 The Strand, Portaferry, BT22 1PF, UKQueen’s UniversityBelfastUnited KingdomSchool of Natural and Built Environment, Queens University Belfast, Belfast, BT9 5AG, UKQueen’s UniversityBelfastUnited Kingdom
Corresponding author: Patrick W. S. Joyce (pjoyce04@qub.ac.uk)
Academic editor: A. Ricciardi
20191504201945193705881962-0984-55B2-BC88-52B181F9EA3026488721501201927032019Patrick W.S. Joyce, Louise Kregting, Jaimie T.A. DickThis is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
The ecological impacts of invasive species can be severe, but are generally viewed as highly unpredictable. Recent methods combining per capita feeding rates, population abundances and environmental contexts have shown great utility in predicting invader impacts. Here, clearance rates of the invasive Pacific oyster, Crassostreagigas, and native mussel, Mytilusedulis, were investigated in a laboratory experiment where oscillatory water flow and algal food concentrations were manipulated. Invasive oysters had lower clearance rates than native mussels in all experimental groups and did not differ among flow velocities or food concentrations. Native mussel clearance rates were higher at 5 cm s-1 compared to 0 and 15 cm s-1 flow velocities and increased with increasing food concentration. The Relative Impact Potential (RIP) metric was used to assess (i) the influence of flow velocity and food concentration on potential impacts of C.gigas on plankton resources and, (ii) the impacts of coexisting reefs, containing both species, on resources compared to monospecific native mussel beds. Greatest Relative Impact Potential of invasive oysters was seen at the lowest flow velocity, but became reduced with increasing flow velocity and food concentration. Relative Impact Potentials of coexisting reefs were generally greater than monospecific native mussel beds, with greatest impacts predicted at lowest flow velocity. We suggest that the greatest ecological impacts and competition potential of C.gigas will occur in areas with low flow velocity, but that increased flow will mediate co-existence between the two species.
Joyce PWS, Kregting LT, Dick JTA (2019) Relative impacts of the invasive Pacific oyster, Crassostrea gigas, over the native blue mussel, Mytilus edulis, are mediated by flow velocity and food concentration. NeoBiota 45: 19–37. https://doi.org/10.3897/neobiota.45.33116
Introduction
The ecological impacts of invasive species can be severe (Simberloff et al. 2013, Dick et al. 2017b) and the rate of invasive species introductions shows no sign of decline (Seebens et al. 2017, 2018). It is therefore necessary to improve our abilities to predict the ecological impacts of current, emerging and future invasive species (Dick et al. 2013, 2017b, Alexander et al. 2014). Recently, the quantification and comparison of per capita effects of invasive species, through for example comparative functional responses (inter- and intraspecific comparisons of consumption in relation to resource densities), have been revealed as a strong predictor of invasive species impacts (Dick et al. 2014, 2017b). Invasive species are often associated with higher consumption rates than comparative native species, with these higher per capita metrics predicting ecological impact (Dick et al. 2013). Further, the recent incorporation of contexts such as temperature (South and Dick 2017; South et al. 2017) or habitat complexity (Alexander et al. 2012; Wasserman et al. 2016) allow such experiments to increase our understanding and prediction of impacts under various environmental contexts (Dick et al. 2014, Paterson et al. 2015).
Recently, the Relative Impact Potential (RIP), a metric proposed by Dick et al. (2017b), has combined per capita resource use with population abundances to better predict the ecological impacts of invasive species. Although Dick et al. (2017b) primarily use functional responses combined with abundance data to produce RIP scores, they suggest that suitable, relevant proxies for such measures can be used in their place. For example, in the case of filter feeders, algal uptake or clearance rates as a measure of per capita resource use are more common than functional responses (e.g. Alexander et al. 2015). Similarly, biomass may be an equally relevant metric as a proxy for population abundance depending upon the species in question (Dick et al. 2017b). The RIP combines facets of the ‘Total Response’ and the ‘Parker-Lonsdale’ equations (see Dick et al. 2017b and Parker et al. 1999, respectively) into one metric that produces absolute values of species impact on a resource. These absolute values can then be used to compare impact, for example, the baseline impacts of native consumers in relation to invasive species. This approach was highly successful in identifying high impact invaders and indeed RIP scores are correlated tightly with independent measures of the degree of ecological impact of such invaders (Dick et al. 2017b). Here, “impact” is defined as a documented effect on a native population, whereas invasion “success” is best defined as the rate of establishment or spread of a species. The semantic distinction between the two terms should be clarified as no link between the two has been found (Ricciardi and Cohen 2007).
Per capita resource use has traditionally been used in animal ecology to investigate impacts on resources (Holling 1959, 1966, Dick et al. 2017b), whereas plant ecologists use the same method to explicitly study interspecific resource competition (Tilman 1977, Dick et al. 2017a). In ecosystems containing sessile animals that cannot move and search for different resources, such comparative resource use may reveal patterns of interspecific competition as well as impacts on resources. Investigations of interactions between sessile organisms often consider space as the only limiting resource worth studying (Connell 1961). However, local seston depletion can occur above bivalves beds (Wildish and Kristmanson 1984, Dolmer 2000a, b) leading to resource limitation (e.g. Vismann et al. 2016), and hence potential competition between filter feeding species due to their limited ability to actively search for new resources.
While comparative per capita resource use has been applied to a range of taxonomic groups, the method has only recently been applied to filter feeders (Alexander et al. 2015, Kemp and Aldridge 2018). However, the incorporation of environmental parameters such as water motion, a fundamental process for filter feeders, have not been included. Sessile suspension feeders rely greatly upon water motion as it is the bulk water column flow that supplies them with fresh food (Genin et al. 1986). Although water motion is necessary for replenishing plankton resources, hydrodynamic forces exerted by the movement of water can also exert destructive forces on organisms, therefore trade-offs between food provision and dislodgement pressures occur (Denny 2006).
The clear majority of work and our understanding of bivalve feeding with regards to water motion has been conducted in uni-directional currents characteristic of estuaries, inland bays and harbours. It is unknown if these studies provide a reasonable basis for the prediction of responses of bivalves to oscillatory water motion characteristic of wind and swell-driven open coasts (Denny and Gaylord 2002). Continuous reversals in flow direction increase turbulence within the water column (Denny et al. 1998), thus it is less likely that seston depletion above bivalve beds would occur. Several studies into the effects of flow velocity on clearance rates of bivalves have been conducted using mussels (Mytilus spp.) but with conflicting results. Some studies have found clearance rates to be unaffected by increasing flow velocity while others show significant reductions in clearance rates with increases in flow (Denis 1999; Newell et al. 2001; Widdows et al. 2002; Ackerman and Nishizaki 2004; Nielsen and Vismann 2014). To date, there have been no studies on the influence of water motion on clearance rates of the invasive Pacific oyster, Crassostreagigas.
The Pacific oyster, C.gigas, is one of the most ‘globalised’ marine invertebrates, dominating shellfish production in many regions (Ruesink 2007, Herbert et al. 2016), and is considered invasive in several countries. For example, the Wadden Sea has seen C.gigas settle onto beds of the native blue mussel, Mytilusedulis, on such a scale that there has been a shift in dominance from native mussels to non-native oysters (Kochmann et al. 2008), suggesting that C.gigas can compete with native M.edulis for resources and potentially impact those resources to the detriment of the wider community.
The present study thus examined the Relative Impact Potentials of the invasive Pacific oyster, C.gigas, and the native blue mussel, M.edulis, in relation to effects of oscillatory flow velocity and algal food concentration on their clearance rates. The experimental treatments simulated environmental conditions experienced on inshore coasts. The main objectives were to: (i) assess the influence of oscillatory flow velocity and food concentration on the clearance rates of the two species; (ii) combine per capita resource use with field biomass, using the RIP metric to identify conditions that may lead to impacts on plankton resources; and (iii) use the RIP metric to compare the impacts on plankton resources of coexisting bivalve beds with those of monospecific native mussel beds on plankton resources.
MethodsBivalve collection
In August 2016, adult Pacific oysters, Crassostreagigas, with a shell length 65–105 mm, were obtained from a local commercial oyster farm, Killough Oysters Ltd. Adult native mussels, Mytilusedulis, with a shell length of 45–50 mm, were collected from an intertidal rocky shore in Strangford Lough, County Down, Northern Ireland (54°28'11.2"N, 5°32'25.4"W). Animals of these sizes were used as they are representative of adult organisms, thus results from the experiments would provide data for mature populations. Animals were housed at Queen’s University Marine Laboratory, Portaferry in large holding tanks (~500L) with through-flowing, sand filtered seawater pumped directly from the adjacent Strangford Lough. Prior to experimental testing, shells were cleaned of any mud and epibionts and returned to the holding tanks for at least 48 hours prior to testing.
Experimental tank system
Clearance rates of the bivalves were determined in an aerated experimental tank system designed to simulate oscillatory water motion (full details of the design can be found in Kregting et al. 2015). The tank system consisted of four tanks where the bivalves were moved back and forth through a stationary body of aerated water to simulate water motion representative of the horizontal oscillatory water motion benthic animals experience at the seabed on shallow inshore coasts. The horizontal oscillatory water motion was simulated in the four laboratory experimental tanks by two horizontal rods mounted above the tanks on a steel frame allowing free oscillatory movement of the rods. The rods were attached to a rotating arm driven by a 12 V car windscreen wiper motor. Two detachable, vertical polypropylene arms with perpendicular base plates were fixed to each rod (arms = 4). Each arm was suspended over a 65 L polypropylene container (60 × 40 × 32 cm). The driving motor was powered using a regulated power supply (Skytronic 0–30 V) which could be altered to control the horizontal velocity of the arms. Three flow scenarios were selected; static (0 cm s-1) and two which oscillated over a distance of 21 cm with amplitudes of 5 and 15 cm s-1.
Either 4 oysters or 10 mussels were attached to the experimental base plates. Different numbers of each species were used in the experiment to keep the area covered by the animals the same, with clearance rates then corrected by biomass (see below). Oysters were attached to the baseplates using cyanoacrylate glue. Mussels were placed onto baseplates and covered with plastic mesh netting to hold them in place allowing natural byssus attachment. Plates with animals attached were placed into 1 µm filtered, UV sterilised seawater for 22 hours to standardise starvation. After the starvation period which allowed sufficient byssus attachment from mussels, the mesh netting was removed from the mussels prior to testing.
Microalgal culture
The microalga Tetraselmissuecica was chosen for the experiment, as plankton of this size (~6–10 µm) (Chrétiennot-Dinet et al. 1986, Hansen et al. 1996) are retained with high efficiency by both species (Bougrier et al. 1997; Ward and Shumway 2004). Algae were cultured in 1 µm filtered, UV sterilised seawater using f/2 media and were on-grown until sufficient stock could be maintained in a 150L bag culture.
Clearance rate experiment
Experimental tanks were filled with 30 L of 1 µm filtered, UV sterilised seawater and aerated at one end to ensure the water was well mixed, but not interfering with the oscillatory movement, allowing use of the clearance equation (see below). The selected animals were subjected to a randomly selected flow velocity for 30 minutes before the addition and mixing of a randomly selected, pre-defined volume of algal monoculture (Table 1). To measure algal depletion in the tanks, two 3 ml water samples were taken immediately after the algal monoculture was sufficiently mixed within the experimental tanks and again after 1 hour. Cell concentration of the water samples were analysed using an electronic particle counter (Coulter Z1). Experimental tanks were emptied, cleaned with freshwater, and rinsed with 1 µm filtered, UV sterilised seawater after each trial. This process was carried out for both species at the three flow velocities (0, 5, 15 cm s-1) and five algal culture volumes (4, 8, 16, 32, 64 ml; corresponding cell concentrations in Table 1) with four replicates per experimental group. Due to the experimental setup, accurate measurement of animal valve gape was not achievable however, all animals were visually inspected during feeding trials. After experimentation, the soft tissue of each animal was removed from the shell and dried at 70 °C for 24 hours to determine the shell-free dry weight (SFDW) of each replicate. Control trials without animals in the experimental tanks (n = 2) were conducted to identify any natural reductions in algal concentration over the feeding period due to sinking.
Volumes of Tetraselmissuecica added to experimental tanks with corresponding initial cell concentrations within experimental tanks for clearance trials (mean ± S.E.).
ml of T.suecica
Cell concentration (cells ml-1) ± S.E.
4
5954 ± 188
8
8198 ± 265
16
13567 ± 342
32
22221 ± 381
64
42003 ± 664
Clearance rate calculations
Due to adequate water mixing within experimental tanks, the ‘clearance method’ (Riisgård et al. 2013) was used to measure the rates of algal consumption of the bivalves. Clearance rates (CR), measured as the volume of water cleared of particles per hour (h) per gram of SFDW (g), were calculated as:
where V is the volume of water in the experimental tank, C0 and Ct are algal concentrations at time 0 and time t, SFDW is the shell-free dry weight of animal flesh in each replicate. SFDW was used to standardise clearance rates between species as, although the area occupied by both species was kept constant, differences in biomass occurred between the two species.
Data analyses
All analyses were performed in R 3.3.1 (R Core Team 2012). One replicate from two separate experimental groups were removed from the analysis due to mussel detachment during the feeding period. A three-factor analysis of variance (ANOVA) compared clearance rates between species (2 levels; C.gigas and M.edulis), among flow velocities (3 levels; 0, 5, 15 cm s-1), and among food concentrations (5 levels; 4, 8, 16, 32, 64 ml of Tetraselmissuecica). Levene’s test for homogeneity of variance (F29,88 = 0.9, p > 0.05) and Shapiro-Wilk’s test for normality (p > 0.05) ensured ANOVA assumptions were met. Significant differences between treatments were compared with Tukey’s honest significant difference post hoc test.
Species biomass and the Relative Impact Potential (RIP) metric
A systematic search of the on-line scientific databases Scopus, Web of Science and Google Scholar was used to collect field biomass data for both Crassostreagigas and Mytilusedulis. All searches were performed in October 2017 using the search terms (Crassostreagigas OR Magallanagigas OR Mytilusedulis) AND (biomass OR abundance OR density) AND (invasive OR non-native OR native). References from retrieved articles were screened for other relevant publications. Literature was selected (Table 2) if biomass estimates were given as total wet weight (WW), shell-free dry weight (SFDW) or ash-free dry weight (AFDW). Data given as WW or AFDW were converted to SFDW using published weight conversion factors for bivalves (Ricciardi and Bourget 1998), as SFDW was used in the clearance rate calculations. Clearance rates averaged across food concentrations, as well as those at the lowest and highest food concentrations for each species from this study, were combined with biomasses for each species to create RIP biplots (Laverty et al. 2017). Biplots represent the Relative Impact Potential of C.gigas compared to M.edulis to under the contexts of ‘flow velocity’ and ‘food concentration’. Biomass data from the Wadden Sea were also available for reefs where the two species are coexisting, dubbed ‘oyssel reefs’ (Reise et al. 2017). In these cases, clearance rates averaged across food concentrations, as well as at the lowest and highest food concentrations of each species, were multiplied by their proportional contribution to the overall reef biomass. The proportionally adjusted clearance rates were then combined to give an overall clearance rate for the mixed species reef. Clearance rates combined with biomass data were used to create RIP biplots to represent Relative Impact Potentials of coexisting ‘oyssel reefs’ compared to monospecific M.edulis beds. RIP biplots combine biomass and clearance rate data to give a visualisation of ecological impact with greater impacts being shifted towards the top and right of the plot, and lesser impacts being shifted towards the bottom and left of the plot (Laverty et al. 2017).
Biomass data for Crassostreagigas, Mytilusedulis, and coexisting ‘oyssel’ reefs from the Wadden Sea with corresponding references. Symbols denote separate species biomass contributions to coexisting ‘oyssel’ reefs.
Species
Biomass (g SFDW m-2)
Reference
C.gigas
508†
(Markert et al. 2010)
348‡
(Markert et al. 2013)
201*
(Markert et al. 2013)
118
(Fey et al. 2010)
M.edulis
328
(Markert et al. 2010)
247†
(Markert et al. 2010)
85‡
(Markert et al. 2013)
71*
(Markert et al. 2013)
166
(Munch-Petersen and Kristensen 2001)
Coexisting reef
755†
(Markert et al. 2010)
433‡
(Markert et al. 2013)
273*
(Markert et al. 2013)
ResultsClearance rate experiment
Visual inspection found that all animals were open and appeared to be feeding during experimental trials. Control groups saw changes in algal concentrations < 2% of the changes that occurred in treatments with animals, thus any changes in algal concentration over the feeding period with animals present were attributed to intake by the animals and not sinking.
Overall, clearance rates of Crassostreagigas were significantly lower than those of Mytilusedulis (Table 3; Figure 1). There was a significant main effect of flow velocity on clearance rate (Table 3), with clearance rates at 5 cm s-1 significantly higher than both other velocities tested (Tukey’s HSD, p < 0.05). However, the significant ‘species’ × ‘flow’ interaction (Table 3) reflects the lack of change in the C.gigas clearance rate but increase in M.edulis clearance rate at 5 cm s-1.
Overall, clearance rate increased with food concentration (Table 3). The significant ‘species’ × ‘food concentration’ interaction (Table 3) reflects that increasing clearance rates over increasing algal food concentrations occurred only for the native M.edulis (Tukey’s HSD, p < 0.05; Fig. 1).
Clearance rates of the native mussel, Mytilusedulis (blue circles), and invasive Pacific oyster, Crassostreagigas (red triangles), as a function of algal food concentration at flow velocities of 0, 5 and 15 cm s-1.
https://binary.pensoft.net/fig/293892
Three-way ANOVA of the effects of species (2 levels; Crassostreagigas and Mytilusedulis), flow velocity (3 levels; 0, 5, 15 cm s-1), and food concentration (5 levels; 4, 8, 16, 32, 64 ml of algal monoculture) on clearance rates.
Df
Mean Sq
F value
Pr(>F)
Species
1
151.24
275.725
< 0.001
Flow
2
3.28
5.985
<0.01
Food
4
5.83
10.636
< 0.001
Species × flow
2
4.3
7.836
< 0.001
Species × food
4
3.83
6.989
< 0.001
Flow × food
8
1.09
1.985
0.0575
Species × flow × food
8
0.71
1.303
0.2523
Residuals
88
0.55
Relative impact potentials
Biomass data from the Wadden Sea show that C.gigas generally has a higher biomass than M.edulis. Combined with average clearance rates from this study, C.gigas is shown to have similar Relative Impact Potential to M.edulis at 0 and 15 cm s-1 flow velocities (Fig. 2). The increased average clearance rate of M.edulis at 5 cm s-1 reveals that impacts of M.edulis at this flow velocity are higher than C.gigas (i.e. shifted further to the right, Fig. 2). At low food concentrations, the RIP of C.gigas is higher than that of M.edulis at both 0 and 5 cm s-1 due to reduced M.edulis clearance rates (i.e. M.edulis toward bottom and shifted left, Fig. 2). At 15 cm s-1, impacts of the two species under low food conditions are similar. High food concentrations indicate greater impacts of M.edulis under all flow scenarios due to increased clearance rates.
RIP biplots of the native mussel, Mytilusedulis (blue), and invasive Pacific oyster, Crassostreagigas (red), using biomass data from the Wadden Sea (mean ± S.E.). Squares indicate clearance rate (CR; L h-1 g-1) at minimum food level, circles indicate average CR over all food levels, triangles indicate CR at maximum food level (mean ± S.E.). Impact increases from bottom left to top right.
https://binary.pensoft.net/fig/293893
The total biomass of coexisting reefs was higher than monospecific M.edulis beds (Table 2). The elevated average clearance rate of M.edulis at 5 cm s-1 leads to the impacts of coexisting reefs and monospecific M.edulis beds to be similar. At 0 and 15 cm s-1, however, the impacts of coexisting reefs are higher than those of monospecific M.edulis beds (Fig. 3). At low food concentrations, impacts of coexisting reefs are shown to be higher than monospecific M.edulis beds at all flow velocities (Fig. 3). High food concentrations lead to similar impacts of coexisting reefs and monospecific M.edulis beds at 0 and 15 cm s-1 but greater impacts of monospecific M.edulis beds at 5 cm s-1 (Fig. 3).
RIP biplots of the native mussel, Mytilusedulis (blue), and coexisting ‘oyssel’ reefs (orange) using biomass data from the Wadden Sea (mean ± S.E.). Squares indicate clearance rate (CR; L h-1 g-1) at minimum food level, circles indicate average CR over all food levels, triangles indicate CR at maximum food level (mean ± S.E.). Impact increases from bottom left to top right.
https://binary.pensoft.net/fig/293894Discussion
Comparative resource use and Relative Impact Potential studies involving native and invasive bivalves to investigate species interactions have not been explored in depth (Alexander et al. 2015, Kemp and Aldridge 2018). This study has thus coupled a comparative resource use concept, with population abundance and the environmental context of oscillatory flow, to examine whether invasion impact on resources and competitive effects by the invasive Pacific oyster, Crassostreagigas, could be predicted based on algal uptake in the different hydrodynamic conditions tested. We found that per capita resource use of the invasive Pacific oyster, C.gigas was lower than that of native Mytilusedulis. However, when accounting for field densities, sites that have seen large invasions of C.gigas may experience ecological impacts on resource communities, especially in areas with little water motion.
The flow velocities chosen in this study are within the range that mussels and oysters are likely to experience regularly in open coastal areas, for example, the Wadden Sea (Janssen-Stelder 2000). Changes in flow velocity and food concentration had no effect on invasive C.gigas but significantly altered native M.edulis clearance rates. Flow velocity and food concentration mediated the Relative Impact Potential of the invasive C.gigas over the native M.edulis. Low flow velocities and food concentrations led to the RIP of C.gigas being higher than that of M.edulis, suggesting a greater impact on resources (i.e. plankton) by C.gigas under such conditions. Increases in flow and food however, increased the RIP of M.edulis due to the increases in clearance rate, thus suggesting a lower comparative impact of C.gigas when flow velocity and food concentrations increase. The RIP biplots also show that in reefs where coexistence between the two species occurs, impacts on plankton resources are likely to be greater than monospecific M.edulis beds in the majority of flow velocity and food concentration contexts tested.
This is the first study investigating the effects of oscillatory water flow on bivalve clearance rates thus we cannot compare the results found to other studies. Previous investigation into M.edulis clearance rates in uni-directional currents have provided mixed results (Ackerman 1999, Denis 1999, Widdows et al. 2002, Ackerman and Nishizaki 2004). Here, we found a uni-modal response of M.edulis clearance rates with increasing water velocity. Although unclear, it may be predicted that clearance rates may decrease at lower oscillatory flow velocities compared to uni-directional current velocities due to increased turbulence created by oscillating motion. Such turbulence may inhibit feeding as well as the fact that in oscillatory flows, inhalant siphons would face into the flow 50% of the time, which has been linked with decreased clearance rates (Newell et al. 2001). The lack of influence of flow velocity on the clearance rate of C.gigas differed from responses shown by M.edulis. No previous studies have investigated the influence of flow velocity on C.gigas clearance rates although this species can be found in environments with a wide range of hydrodynamic conditions from high energy to sheltered environments (Wrange et al. 2010, Strand et al. 2012, Dolmer et al. 2014).
Food concentration only significantly increased M.edulis clearance rate. This is consistent with patterns observed whereby, at a lower threshold, bivalves can cease filtering (Denis 1999, Riisgård et al. 2013, Sarnelle et al. 2015). This is not shown by C.gigas which, although it has a lower clearance rate, appears to maximise its feeding capability even at low food levels.
Here, our measured clearance rates for C.gigas were <1 L h-1 g-1, which is lower than other studies ranging from 2–11.8 L h-1 g-1 (Walne et al. 1972, Gerdes et al. 1983, Bougrier et al. 1995, Dupuy et al. 2000). These studies, however, used a range of oyster sizes generally smaller than those used here, which may lead to higher body weight specific clearance rates. Previous studies also measured clearance rates in static systems which provide unrealistic, idealised conditions for filtration. Although our measured clearance rates are lower than other laboratory studies, they are similar to field observations of < 1 L h-1 g-1 (Wheat and Ruesink 2013). Our measured clearance rates for M.edulis are comparable to those of other laboratory studies (examples in Troost et al. 2010).
Although C.gigas is a successful invader, it produces varied ecological impacts, both positive and negative depending on context (Padilla 2010, Herbert et al. 2012, 2016). Invasion “success” and “impact” should be distinguished as the two are not necessarily correlated (Ricciardi and Cohen 2007). Common misinterpretation leads to the incorrect use of the terms whereby success should be defined as the rate of establishment and spread whereas impact is a documented effect on native populations. Here, we show that the relatively low per capita clearance rates of C.gigas found are in line with theory that high relative per capita rates are associated with high impact, with the corollary being that low impact should be associated with low per capita rates (Dick et al. 2013, 2014, 2017b). By combing per capita clearance rates with field biomass data into the RIP metric, we show again that, even with the higher biomass of C.gigas, because this is mitigated by lower per capita feeding rates, the invader is predicted to have relatively low impact on native resources.
Further, although the RIPs would usually be used to assess or predict species impacts on a resource, we contend that it may also be useful in understanding interspecific competition (Dick et al. 2017a), due to the inability of bivalves to move and search for new food resources. Although plankton resources are not thought of to be limiting, it has been shown that seston depletion can occur above bivalve beds (Wildish and Kristmanson 1984, Dolmer 2000a), and that bivalve beds can become food limited (Vismann et al. 2016), which may result in interspecific competition for limited resources. The RIPs shown here suggest that only under low flow and low food conditions may the invader, C.gigas, exhibit ecological impacts over, and compete effectively with the native, M.edulis. Such areas are therefore most likely to see resource based ecological impacts from C.gigas. However, due to daily fluctuations in wind/storm driven changes in flow velocity as well as seasonal variations in plankton abundance, M.edulis will not always be outcompeted for resources which is reflected by the regular coexistence found between the two species (Holm et al. 2016; Reise et al. 2017). Although these species may compete for other resources (i.e. space), our data, using measured clearance rates, are fully in line with field patterns of invasive impacts and coexistence, and indeed, the per capita, RIP and context-dependency approach could be used more predictively for emerging and potential invasive species impacts (Dick et al. 2017b).
The RIP metric also revealed that the impacts of mixed species ‘oyssel reefs’ (Reise et al. 2017) on resources are predicted to be greater than baseline impacts of monospecific M.edulis beds. The magnitude of these impacts is shown to be greatest in areas with low flow velocities and food concentrations. A lack of investigation into differences in plankton resources over native monospecific M.edulis beds compared to those that have been invaded by C.gigas cannot allow statements of whether such field impacts have occurred. However, with the RIP clearly a successful predictive tool for invasive species impacts on resources, studies into plankton abundances above native and non-native bivalve beds clearly need further investigation.
Conclusion
Due to the differential effects of flow velocity on C.gigas and M.edulis, additional investigations into the effects of flow type (i.e. currents vs. waves) may be required to further understand differences in feeding and growth in situ, especially as these species are sensitive to interspecific competition. As growth is directly related to feeding, investigation into growth rates in different hydrodynamic conditions, which vary both naturally and due to anthropogenic influences, should be conducted to validate the results of this study. Based on our results, we suggest that areas with little water motion and those where local food limitation may occur are likely to be most at risk of impacts from C.gigas. Areas with increased water motion are unlikely to be resource limited due to increased food replenishment however, destructive forces exerted by the water itself may affect species growth and success. The field patterns of low invasion impact and coexistence with the native analogue, M.edulis, further highlight the excellent explanatory and predictive power of coupling per capita resource use with field abundances for invasion ecology, however, investigation into growth rates under these environmental contexts would achieve an even better understanding of competition between the two species.
Acknowledgements
All authors conceived the work and designed the experiment. PWSJ carried out data collection and analysis. All authors contributed to writing the final manuscript and gave approval for publication. The authors thank Aidan Flaherty for assistance in adapting the experimental setup, Emma Gorman for providing the algal culture and Daniel Barrios-O’Neill for advice regarding statistical analyses. JTAD also thanks the Natural Environment Research Council (NERC). This work was funded by the Department for the Economy (DfE) N. Ireland and part of SuperGen Marine Energy Research Consortium II, which was funded by the UK Engineering and Physical Science Research Council. The authors extend their gratitude to Andrea Anton Gamazo and an anonymous reviewer for their helpful comments which have improved the manuscript.
ReferencesAckermanJD (1999) Effect of velocity on the filter feeding of dreissenid mussels (Dreissenapolymorpha and Dreissenabugensis): implications for trophic dynamics.56: 1551–1561. https://doi.org/10.1139/f99-079AckermanJDNishizakiMT (2004) The effect of velocity on the suspension feeding and growth of the marine mussels Mytilustrossulus and M.californianus: implications for niche separation.49: 195–207. https://doi.org/10.1016/j.jmarsys.2003.06.004AlexanderMEAdamsRDickJTARobinsonTB (2015) Forecasting Invasions: resource use by mussels informs invasion patterns along the South African coast.162: 2493–2500. https://doi.org/10.1007/s00227-015-2742-5AlexanderMMEDickJO’ConnorNHaddawayNFarnsworthK (2012) Functional responses of the intertidal amphipod Echinogammarusmarinus: effects of prey supply, model selection and habitat complexity.468: 191–202. https://doi.org/10.3354/meps09978AlexanderMEDickJTAWeylOLFRobinsonTBRichardsonDM (2014) Existing and emerging high impact invasive species are characterized by higher functional responses than natives.10(2): 20130946. https://doi.org/10.1098/rsbl.2013.0946BougrierSGeaironPDeslous-PaoliJMBacherCJonquièresG (1995) Allometric relationships and effects of temperature on clearance and oxygen consumption rates of Crassostreagigas (Thunberg).134: 143–154. https://doi.org/10.1016/0044-8486(95)00036-2BougrierSHawkinsAJSHéralM (1997) Preingestive selection of different microalgal mixtures in Crassostreagigas and Mytilusedulis, analysed by flow cytometry.150: 123–134. https://doi.org/10.1016/S0044-8486(96)01457-3ConnellJH (1961) The Influence of interspecific competition and other factors on the distribution of the barnacle Chthamalus Stellatus.42: 710–723. https://doi.org/10.2307/1933500DenisL (1999) Clearance rate responses of Mediterranean mussels, Mytilusgalloprovincialis, to variations in the flow, water temperature, food quality and quantity.12: 279–288. https://doi.org/10.1016/S0990-7440(00)86639-5DennyMWGaylordB (2002) The mechanics of wave-swept algae.205: 1355–1362DennyMWGaylordBHelmuthBDanielT (1998) The menace of momentum: Dynamic forces on flexible organisms.43: 955–968. https://doi.org/10.4319/lo.1998.43.5.0955DickJTAAlexanderMEJeschkeJMRicciardiAMacIsaacHJRobinsonTBKumschickSWeylOLFDunnAMHatcherMJPatersonRAFarnsworthKDRichardsonDM (2014) Advancing impact prediction and hypothesis testing in invasion ecology using a comparative functional response approach.16: 735–753. https://doi.org/10.1007/s10530-013-0550-8DickJTAAlexanderMERicciardiALavertyCDowneyPOXuMJeschkeJMSaulW-CHillMPWassermanRBarrios-O’NeillDWeylOLFShawRH (2017a) Functional responses can unify invasion ecology.19: 1667–1672. https://doi.org/10.1007/s10530-016-1355-3DickJTAGallagherKAvlijasSClarkeHCLewisSELeungSMinchinDCaffreyJAlexanderMEMaguireCHarrodCReidNHaddawayNRFarnsworthKDPenkMRicciardiA (2013) Ecological impacts of an invasive predator explained and predicted by comparative functional responses.15: 837–846. https://doi.org/10.1007/s10530-012-0332-8DickJTALavertyCLennonJJBarrios-O’NeillDMensinkPJRobert BrittonJMédocVBoetsPAlexanderMETaylorNGDunnAMHatcherMJRosewarnePJCrookesSMacIsaacHJXuMRicciardiAWassermanRJEllenderBRWeylOLFLucyFEBanksPBDoddJAMacNeilCPenkMRAldridgeDCCaffreyJM (2017b) Invader Relative Impact Potential: a new metric to understand and predict the ecological impacts of existing, emerging and future invasive alien species.54: 1259–1267. https://doi.org/10.1111/1365-2664.12849DolmerP (2000a) Feeding activity of mussels Mytilusedulis related to near-bed currents and phytoplankton biomass.44: 221–231. https://doi.org/10.1016/S1385-1101(00)00052-6DolmerP (2000b) Algal concentration profiles above mussel beds.43: 113–119. https://doi.org/10.1016/S1385-1101(00)00005-8DolmerPHolmMWStrandÅLindegarthSBodvinTNorlingPMortensenS (2014) The invasive Pacific oyster, Crassostreagigas, in Scandinavia coastal waters: A risk assessment on the impact in different habitats and climate conditions. Institute of Marine Research (Fisken og Havet, Vol 2).DupuyCVaquerALamHoai TRougierCMazouniNLautierJCollosYGallS Le (2000) Feeding rate of the oyster Crassostreagigas in a natural planktonic community of the Mediterranean Thau Lagoon.205: 171–184. https://doi.org/10.3354/meps205171FeyFDankersNSteenbergenJGoudswaardK (2010) Development and distribution of the non-indigenous Pacific oyster (Crassostreagigas) in the Dutch Wadden Sea.18: 45–59. https://doi.org/10.1007/s10499-009-9268-0GerdesD (1983) The Pacific oyster Crassostreagigas: Part I. Feeding behaviour of larvae and adults.31: 195–219. https://doi.org/10.1016/0044-8486(83)90313-7HerbertRJHHumphreysJDaviesCJRobertsCFletcherSCroweTP (2016) Ecological impacts of non-native Pacific oysters (Crassostreagigas) and management measures for protected areas in Europe.25: 2835–2865. https://doi.org/10.1007/s10531-016-1209-4HerbertRJHRobertsCHumphreysJFletcherS (2012) The Pacific oyster (Crassostreagigas) in the UK: Economic, legal and environmental issues associated with its cultivation, wild establishment and exploitation. Report for the Shellfish Association of Great Britain.HollingCS (1959) The components of predation as revealed by a study of small-mammal predation of the European pine sawfly.91: 293–320. https://doi.org/10.4039/Ent91293-5HollingCS (1966) The functional response of invertebrate predators to prey density.98: 1–86. https://doi.org/10.4039/entm9848fvHolmMWDavidsJKDolmerPHolmesETheis NielsenTVismannBWinding HansenB (2016) Coexistence of Pacific oyster, Crassostreagigas (Thunberg, 1793), and blue mussels, Mytilusedulis (Linnaeus, 1758), on a sheltered intertidal bivalve bed? Aquatic Invasions 11: 155–165. https://doi.org/10.3391/ai.2016.11.2.05Janssen-StelderB (2000) The effect of different hydrodynamic conditions on the morphodynamics of a tidal mudflat in the Dutch Wadden Sea.20: 1461–1478. https://doi.org/10.1016/S0278-4343(00)00032-7KempJSAldridgeDC (2018) Comparative functional responses to explain the impact of sympatric invasive bivalves (Dreissena spp.) under different thermal regimes.84: 175–181. https://doi.org/10.1093/mollus/eyy006KochmannJBuschbaumCVolkenbornNReiseK (2008) Shift from native mussels to alien oysters: Differential effects of ecosystem engineers.364: 1–10. https://doi.org/10.1016/j.jembe.2008.05.015KregtingLTHepburnCDSavidgeG (2015) Seasonal differences in the effects of oscillatory and uni-directional flow on the growth and nitrate-uptake rates of juvenile Laminariadigitata (Phaeophyceae).51: 1116–1126. https://doi.org/10.1111/jpy.12348LavertyCGreenKDDickJTABarrios-O’NeillDMensinkPJMédocVSpataroTCaffreyJMLucyFEBoetsPBrittonJRPeggJGallagherC (2017) Assessing the ecological impacts of invasive species based on their functional responses and abundances.19: 1653–1665. https://doi.org/10.1007/s10530-017-1378-4MarkertAEsserWFrankDWehrmannAExoKM (2013) Habitat change by the formation of alien Crassostrea-reefs in the Wadden Sea and its role as feeding sites for waterbirds.131: 41–51. https://doi.org/10.1016/j.ecss.2013.08.003MarkertAWehrmannAKrönckeI (2010) Recently established Crassostrea-reefs versus native Mytilus-beds: differences in ecosystem engineering affects the macrofaunal communities (Wadden Sea of Lower Saxony, southern German Bight).12: 15–32. https://doi.org/10.1007/s10530-009-9425-4Munch-PetersenSKristensenS (2001) On the dynamics of the stocks of blue mussels (Mytilusedulis L.) in the Danish Wadden Sea.465: 31–43. https://doi.org/10.1023/A:1014539414345NewellCRWildishDMacDonaldB (2001) The effects of velocity and seston concentration on the exhalant siphon area, valve gape and filtration rate of the mussel Mytilusedulis.262: 91–111. https://doi.org/10.1016/S0022-0981(01)00285-4NielsenPVismannB (2014) Clearance rate of Mytilusedulis (L.) as a function of current velocity and mussel aggregation.33: 457–463. https://doi.org/10.2983/035.033.0214PadillaDK (2010) Context-dependent impacts of a non-native ecosystem engineer, the Pacific oyster Crassostreagigas.50: 213–225. https://doi.org/10.1093/icb/icq080ParkerIMSimberloffDLonsdaleWMGoodellKWonhamMKareivaPMWilliamsonMHHolleB VonMoylePBByersJEGoldwasserL (1999) Impact: Toward a framework for understanding the ecological effects of invaders.1: 3–19. https://doi.org/10.1023/A:1010034312781PatersonRADickJTAPritchardDWEnnisMHatcherMJDunnAM (2015) Predicting invasive species impacts: a community module functional response approach reveals context dependencies.84: 453–463. https://doi.org/10.1111/1365-2656.12292R Core Team (2012) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna.ReiseKBuschbaumCBüttgerHWegnerKM (2017) Invading oysters and native mussels: from hostile takeover to compatible bedfellows. Ecosphere 8: e01949. https://doi.org/10.1002/ecs2.1949RicciardiABourgetE (1998) Weight-to-weight conversion factors for marine benthic macroinvertebrates.163: 245–251. https://doi.org/10.3354/meps163245RicciardiACohenJ (2007) The invasiveness of an introduced species does not predict its impact.9: 309–315. https://doi.org/10.1007/s10530-006-9034-4RiisgårdHUPleissnerDLundgreenKLarsenPS (2013) Growth of mussels Mytilusedulis at algal (Rhodomonassalina) concentrations below and above saturation levels for reduced filtration rate.9: 1005–1017. https://doi.org/10.1080/17451000.2012.742549RuesinkJ (2007) Biotic resistance and facilitation of a non-native oyster on rocky shores.331: 1–9. https://doi.org/10.3354/meps331001SarnelleOWhiteJDGeelhoedTEKozelCL (2015) Type III functional response in the zebra mussel, Dreissenapolymorpha.72: 1202–1207. https://doi.org/10.1139/cjfas-2015-0076SeebensHBlackburnTMDyerEEGenovesiPHulmePEJeschkeJMPagadSPyšekPKleunenM vanWinterMAnsongMArianoutsouMBacherSBlasiusBBrockerhoffEGBrunduGCapinhaCCaustonCECelesti-GrapowLDawsonWDullingerSEconomoEPFuentesNGuénardBJägerHKarteszJKenisMKühnILenznerBLiebholdAMMosenaAMoserDNentwigWNishinoMPearmanDPerglJRabitschWRojas-SandovalJRoquesARorkeSRossinelliSRoyHEScaleraRSchindlerSŠtajerováKTokarska-GuzikBWalkerKWardDFYamanakaTEsslF (2018) Global rise in emerging alien species results from increased accessibility of new source pools. Proceedings of the National Academy of Sciences of the USA. https://doi.org/10.1073/pnas.1719429115SeebensHBlackburnTMDyerEEGenovesiPHulmePEJeschkeJMPagadSPyšekPWinterMArianoutsouMBacherSBlasiusBBrunduGCapinhaCCelesti-GrapowLDawsonWDullingerSFuentesNJägerHKarteszJKenisMKreftHKühnILenznerBLiebholdAMosenaAMoserDNishinoMPearmanDPerglJRabitschWRojas-SandovalJRoquesARorkeSRossinelliSRoyHEScaleraRSchindlerSŠtajerováKTokarska-GuzikBKleunenM vanWalkerKWeigeltPYamanakaTEsslF (2017) No saturation in the accumulation of alien species worldwide. Nature Communications 8: 14435. https://doi.org/10.1038/ncomms14435SimberloffDMartinJ-LGenovesiPMarisVWardleDAAronsonJCourchampFGalilBGarcía-BerthouEPascalMPyšekPSousaRTabacchiEVilàM (2013) Impacts of biological Invasions: what’s what and the way forward.28: 58–66. https://doi.org/10.1016/j.tree.2012.07.013SouthJDickJTA (2017) Effects of acute and chronic temperature changes on the functional responses of the dogfish Scyliorhinuscanicula (Linnaeus, 1758) towards amphipod prey Echinogammarusmarinus (Leach, 1815). Environmental Biology of Fishes: 1–13. https://doi.org/10.1007/s10641-017-0633-ySouthJDickJTAMcCardMBarrios-O’NeillDAntonA (2017) Predicting predatory impact of juvenile invasive lionfish (Pteroisvolitans) on a crustacean prey using functional response analysis: effects of temperature, habitat complexity and light regimes. Environmental Biology of Fishes: 1–11. https://doi.org/10.1007/s10641-017-0633-yStrandÅBlandaEBodvinTDavidsJKFastJensen LHolm-HansenTHJelmertALindegarthSMortensenSMoyFENielsenPNorlingPNybergCChristensenHTVismannBHolmMWWinding HansenBDolmerP (2012) Impact of an icy winter on the Pacific oyster (Crassostreagigas Thunberg, 1793) populations in Scandinavia.7: 433–440. https://doi.org/10.3391/ai.2012.7.3.014TilmanD (1977) Resource competition between plankton algae: An experimental and theoretical approach.58: 338–348. https://doi.org/10.2307/1935608TroostK (2010) Causes and effects of a highly successful marine invasion: Case-study of the introduced Pacific oyster Crassostreagigas in continental NW European estuaries.64: 145–165. https://doi.org/10.1016/j.seares.2010.02.004VismannBHolmMWDavidsJKDolmerPPedersenMFBlandaEChristensenHTNielsenPHansenBW (2016) Field clearance of an intertidal bivalve bed: relative significance of the co-occurring blue mussel Mytilusedulis and Pacific oyster Crassostreagigas.25: 107–119. https://doi.org/10.3354/ab00661WalnePR (1972) The influence of current speed, body size and water temperature on the filtration rate of five species of bivalves.52: 345–374. https://doi.org/10.1017/S0025315400018737WardJEShumwaySE (2004) Separating the grain from the chaff: particle selection in suspension- and deposit-feeding bivalves.300: 83–130. https://doi.org/10.1016/j.jembe.2004.03.002WassermanRJAlexanderMEWeylOLFBarrios-O’NeillDFronemanPWDaluT (2016) Emergent effects of structural complexity and temperature on predator-prey interactions. Ecosphere 7: e01239. https://doi.org/10.1002/ecs2.1239WheatERuesinkJL (2013) Commercially-cultured oysters (Crassostreagigas) exert top-down control on intertidal pelagic resources in Willapa Bay, Washington, USA.81: 33–39. https://doi.org/10.1016/j.seares.2013.04.006WiddowsJLucasJSBrinsleyMDSalkeldPNStaffFJ (2002) Investigation of the effects of current velocity on mussel feeding and mussel bed stability using an annular flume.56: 3–12. https://doi.org/10.1007/s10152-001-0100-0WildishDJKristmansonDD (1984) Importance to mussels of the benthic boundary layer.41: 1618–1625. https://doi.org/10.1139/f84-200WrangeA-LValeroJHarkestadLSStrandØLindegarthSChristensenHTDolmerPKristensenPSMortensenS (2010) Massive settlements of the Pacific oyster, Crassostreagigas, in Scandinavia.12: 1145–1152. https://doi.org/10.1007/s10530-009-9535-z