Hydrobiologia (2018) 819:259–273 https://doi.org/10.1007/s10750-018-3644-6 PRIMARY RESEARCH PAPER Trophic interactions in a lowland river ﬁsh community invaded by European barbel Barbus barbus (Actinopterygii, Cyprinidae) Catherine Gutmann Roberts J. Robert Britton Received: 6 March 2018 / Revised: 2 May 2018 / Accepted: 8 May 2018 / Published online: 22 May 2018 The Author(s) 2018 Abstract Following their invasion, non-indigenous driven by inter-speciﬁc differences in functional ﬁsh can potentially exclude native ﬁshes from their morphology and habitat use. original niches through competition, or can partition in their resource use with native species to facilitate co- Keywords Biological invasion Isotopic niche existence. Here, using stable isotope analysis, the Non-indigenous Stable isotope analysis Trophic trophic interactions of invasive European barbel niche Barbus barbus and other ﬁshes were tested in an invaded river of relatively low ﬁsh species diversity and where no other Barbus species were present. Testing was over three distinct life stages: age Introduction 0 ? (\ 38 mm), juveniles (86–231 mm) and adults ([ 386 mm). There were strong patterns of isotopic Invasions of non-indigenous ﬁshes can result in niche partitioning between the juvenile ﬁshes, with adverse impacts in native ﬁsh communities, including some inter-speciﬁc niche differences also apparent in competitive displacement and exclusion (Gozlan 0 ? ﬁshes. For adult B. barbus and chub Squalius et al., 2010). Understanding how an invasive ﬁsh can cephalus, however, niche convergence was evident. impact native species requires knowledge on their Within the B. barbus population, the niches of the trophic interactions, such as whether they share prey adult ﬁsh differed signiﬁcantly from the 0? and resources, resulting in niche convergence, or exploit juvenile ﬁsh, indicating considerable dietary changes different resources, resulting in niche partitioning with development. These results suggested that niche (Cucherousset et al., 2012; Tran et al., 2015; Copp partitioning at the most abundant life stages were et al., 2017). Quantifying the feeding relationships of facilitating the co-existence of invasive B. barbus with introduced and native ﬁshes is thus important for other ﬁshes in the community, with this most likely understanding their ecological risks to the native communities (Cucherousset & Olden, 2011) and facilitates assessment of the ecological impacts that might develop (Gozlan et al., 2010; Tran et al., 2015; Handling editor: Michael Power Copp et al., 2017). C. G. Roberts J. R. Britton (&) European barbel Barbus barbus (Linneaus, 1758) Department of Life and Environmental Sciences, Faculty of the Cyprinidae family is now invasive in many of Science and Technology, Bournemouth University, European rivers outside of its native range, especially Poole BH12 5BB, Dorset, UK e-mail: firstname.lastname@example.org 123 260 Hydrobiologia (2018) 819:259–273 rivers in Italy and Western Britain (Britton & Pegg, indigenous populations, with populations being rela- 2011). Attaining lengths to approximately 800 mm tively abundant and widespread through basins such as and weights in excess of 8 kg (Amat-Trigo et al., the River Severn (Amat-Trigo et al., 2017). Although 2017), they are generally valued for sport angling, knowledge of invasive B. barbus trophic interactions with this the primary driver of introductions (Britton & with indigenous ﬁshes is limited in these rivers, both Pegg, 2011). A highly vagile species, they can disperse Basˇic´ et al. (2015) and Gutmann Roberts et al. (2017) relatively quickly through river systems (Hunt & revealed that in rivers in both their native and invasive Jones, 1974), often leading to rapid colonisation ranges, angling baits based on marine-derived nutri- (Carosi et al., 2017). In invaded rivers where they ents can provide a strong trophic subsidy. This results are sympatric with endemic Barbus ﬁshes, such as in in some individual B. barbus and S. cephalus (gener- Italian rivers (e.g. the Tiber basin), long-term data ally [ 400 mm) specialising on this allochthonous suggest populations of endemic Barbus tyberinus resource. However, in rivers where angling is less Bonaparte, 1839 are being displaced by B. barbus, intense and so where this subsidy is lower, and in body with the mechanism suggested to involve asymmetric sizes that rarely consume these baits (\ 400 mm), competition between the ﬁshes (Carosi et al., 2017). there remains a distinct knowledge gap on the trophic This displacement is in addition to genetic impacts relationships of invasive B. barbus with other species. caused by introgression that results in a loss of genetic Moreover, there is also minimal knowledge on how integrity in the endemic Barbus ﬁshes (Meraner et al., their diet and trophic niche sizes change with increas- 2013; Zaccara et al., 2014). Elsewhere in Europe, ing body size, and in relation to these changes in the invasive B. barbus populations are often present in indigenous ﬁshes. This is despite the diet of ﬁsh communities where other Barbus ﬁshes are absent. usually being gape limited, where gape size is a Thus, whilst they are sympatric with indigenous function of body length (e.g. Persson et al., 1996), cyprinid ﬁshes such as chub Squalius cephalus suggesting considerable dietary shifts will occur with (Linnaeus, 1758), they have lower functionally sim- increasing body length. Data on the inter- and intra- ilarity with these ﬁshes than with congeners. Conse- speciﬁc trophic relationships of B. barbus are also quently, the strength of their interactions might be less missing in their invasive range more generally, where intense and their invasion might be less likely to incur competitive interactions between invasive and ende- negative ecological impacts. mic Barbus ﬁshes have, to date, been inferred from Examples of systems invaded by B. barbus and relative body condition data (e.g. Carosi et al., 2017). where endemic Barbus ﬁshes are absent are rivers in The aim of this study was to quantify the trophic Western England. In Britain, B. barbus is only interactions of a population of invasive B. barbus with indigenous to eastern ﬂowing rivers in England due other ﬁshes in a river where no other Barbus ﬁshes to their previous connections with mainland Europe at were present. The focus was on determining the extent the end of the last glacial period (Wheeler & Jordan, of trophic niche sharing within and between species, 1990). Research on these indigenous B. barbus and how this altered across a range of life stages (as suggests many of these populations are imperilled inferred from body sizes). The River Teme, western due to losses of habitat and river connectivity (Basˇic´ England, was the study river, where non-indigenous B. et al., 2017). Consequently, enhancement stocking barbus have been present since the 1970s (Antognazza often supports these populations, with hatchery-reared et al., 2016). The objective was to determine the individuals released at lengths between 120 and trophic niche sizes and overlaps between invasive B. 250 mm and age 1? and 2? years (Britton et al., barbus and native ﬁshes at three different life stages: 2004; Antognazza et al., 2016). Studies on the trophic age 0? (young-of-the-year), juveniles and adults. It interactions of these stocked ﬁsh suggest substantial was hypothesised that (1) due to the consistent patterns partitioning in their trophic niches with S. cephalus, of inter-speciﬁc trophic partitioning between B. bar- the species that has the most similar functional traits bus and native cyprinid ﬁshes in their indigenous and body sizes as B. barbus in these rivers (Basˇic & range (Basˇic & Britton, 2016), these patterns of inter- Britton, 2016). speciﬁc partitioning are present in their non-indige- In their invasive range in Western England, pop- nous, invasive range; and (2) within the ﬁshes, there ulations tend to be more successful than many were signiﬁcant shifts in the position of the trophic 123 Hydrobiologia (2018) 819:259–273 261 niches across the three life stages, with populations occasionally present in samples but not included in having a relatively large niche comprising smaller analyses were bullhead Cottus gobio Linnaeus, 1758 sub-sets. and stone loach Barbatula barbatula (Linnaeus, As the B. barbus population of the River Teme (and 1758). Brown trout Salmo trutta Linnaeus, 1758 and the River Severn basin generally) is an important juvenile Atlantic salmon Salmo salar Linnaeus, 1758 angling resource (Amat-Trigo et al., 2017), the use of are also present in the river but are more prevalent stomach contents analysis via destructive sampling of upstream of the town of Ludlow, outside of the study the juvenile and adult ﬁsh was not possible. Conse- reach (Fig. 1). Compared with the area of river located quently, trophic analyses were based on stable isotope close to the conﬂuence with the River Severn and that analysis (SIA), where the ecological application of was used by Gutmann Roberts et al. (2017), angling 13 15 carbon (as d C) and nitrogen (as d N) stable isotopes pressure was relatively light in the study reach, and is based on the predictable relationship between the thus inputs of angling baits containing high propor- isotope composition of a consumer and its prey. It thus tions of pelletized ﬁshmeal were considered as rela- provides a temporally integrative and powerful tool to tively low. analyse trophic interactions between native and non- The 0? ﬁsh were sampled from a single area of native ﬁshes (Cucherousset et al., 2012). For compar- nursery habitat located close to Bransford (Fig. 1). isons of SI data within and between the ﬁshes, two They were sampled using a micromesh seine net metrics were used: the signiﬁcance of differences in (25 9 2 m) on 12 September 2016. The ﬁsh were 15 13 d N-d C centroids, and core isotopic niche sizes (as euthanised via anaesthetic overdose (MS-222) and a proxy of the trophic niche) and overlaps calculated transported back to the laboratory on ice. In the using standard ellipse areas (Jackson et al., laboratory, they were identiﬁed to species, measured 2011, 2012). (standard length, nearest mm) and a sample of dorsal muscle tissue removed and dried to constant weight at 50C. The juvenile and adult ﬁsh samples were Methods collected using angling and electric ﬁshing between July and September 2015 and 2016, with SIA based on Sampling details and stable isotope analysis scales (Busst & Britton, 2016, 2017). Correspond- ingly, for each captured ﬁsh, identiﬁcation was to The study was conducted on the middle reaches of the species level, followed by measuring (fork length, River Teme, from the town of Tenbury Wells nearest mm) and the collection of between three and 0 0 (5219 N, - 224 W) to the village of Bransford ﬁve scales from the area between the base of the dorsal 0 0 (5210 N, - 216 W) (Fig. 1). Across the Teme ﬁn and above the lateral line. As scales grow as ﬁsh catchment, altitude varies between 24.3 and 544.5 length increases, only the outer portion of scales mAOD, and land-use is primarily grassland (59%), reﬂects their most recent growth (Hutchinson & with some horticulture (24%) (CEH, 2018). In the Trueman, 2006; Basˇic´ et al., 2015). Consequently, study reach, the river generally comprised sequences only the very outer portion of the sampled scales was of pool and rifﬂes, where maximum depths rarely used in SIA. One scale was prepared per ﬁsh, with this exceeded 2 m and widths rarely exceeded 15 m. A involving their thorough washing with distilled water, ﬂow gauging station towards the downstream end of removal of the scale outer edge using dissection scissors and then drying to constant weight as per the the reach near Bransford had a long-term Q95 of 3 -1 3 -1 3 2.0 m s , Q50 of 10.2 m s and Q10 of 42.4 m 0? ﬁsh samples. The other scale samples were used to -1 s (CEH, 2018). In the study reach, the cyprinid ﬁsh age the ﬁsh (Amat-Trigo et al., 2017). Scale decalci- community had relatively limited diversity, with only ﬁcation was not performed prior to SIA, since the invasive B. barbus, and S. cephalus, dace Leuciscus removal of inorganic carbonates has no signiﬁcant 13 15 leuciscus (Linnaeus, 1758) and minnow Phoxinus effect on scale d C and d N values (Sinnatamby phoxinus (Linnaeus, 1758) present. Grayling Thymal- et al., 2007, 2010; Woodcock & Walther, 2014). lus thymallus (Linnaeus, 1758) were also present in Concomitantly, qualitative samples for SIA of samples at the upper end of the reach and so were also macroinvertebrates were collected from two areas of 0 0 included in analyses. Other species that were the river, ‘Area 1 and ‘Area 2 . Samples in Area 1 123 262 Hydrobiologia (2018) 819:259–273 Fig. 1 Inset: location of the River Teme in Great Britain. Main study area was the stretch of the river between the two dashed map: The River Teme catchment showing its conﬂuence with lines. Macroinvertebrates were collected at locations marked the River Severn. Arrows mark the direction of river ﬂow. The with asterisks were collected from Tenbury Wells (5219 N, spp. larvae (n = 3 per Area) were also taken in 2016. 0 0 0 - 224 W) and Lindridge (5232 N, - 251 W), and All samples were taken back to the laboratory where 0 0 from Area 2 at Bransford (5210 N, - 216 W) in they were washed in distilled water and dried to June and September 2015 and 2016. Samples were constant weight as per the ﬁsh samples; note that in collected using kick sampling. Macro-invertebrate each case, one sample comprised between three and samples collected in 2015 contained very high six individuals. proportions of the amphipod Gammarus pulex (Lin- The dried muscle, scale and invertebrate samples naeus, 1758) in both sampling areas ([ 50%). Gam- were then submitted to the Cornell Isotope Laboratory marus spp. are common prey items for riverine ﬁshes in New York, USA, for SIA. This involved the generally (e.g. MacNeil et al., 1999), as well as the samples being ground to powder, weighed in tin ﬁshes analysed here more speciﬁcally (e.g. Mann, capsules (nearest 1,000 lg) and analysed on a Thermo 1974; Basˇic´ et al., 2015). Thus, samples were taken to Delta V isotope ratio mass spectrometer (Thermo describe the stable isotope data of ﬁsh putative prey in Scientiﬁc, USA) interfaced to a NC2500 elemental 2015. This sampling was repeated in 2016, with G. analyser (CE Elantach Inc. USA). Standards were pulex samples taken for SIA to enable consistent veriﬁed against international reference materials and temporal and spatial testing of differences in ﬁsh calibrated against the primary reference scales for 13 15 putative prey resources. However, to increase the d C and d N. The accuracy and precision were diversity of these baseline samples, samples of checked every ten samples using a standard animal Chironomid larvae (n = 6 per Area) and Trichoptera sample (mink). The outputs were values of d C and 123 Hydrobiologia (2018) 819:259–273 263 d N(%) for each sample. As C:N ratios were below correction on the stable isotope data, differences were 3.5, indicating low lipid content, there was no need for tested in the temporal data in the uncorrected and the d C to be lipid corrected (Post et al., 2007; Skinner corrected data of the juvenile ﬁshes using ANOVA. et al., 2016). The juvenile ﬁshes were used in preference to the adult ﬁshes for this, as the diet of the latter was also likely to Data analysis have had some inﬂuence from angling baits containing ˇ ´ marine-derived ﬁshmeal (Basic et al., 2015; Gutmann The 0? ﬁsh utilised in the analysis were all between Roberts et al., 2017). This testing was not completed 17 and 38 mm and in their ﬁrst year of life. The for the 0? ﬁshes, as their samples were taken from a juvenile ﬁsh were between 86 and 231 mm and single site in 2016. However, their data were also between ages 1? and 4? years; note that in this corrected to enable their results to be compared with length range, some L. leuciscus would have been the juvenile and adult ﬁshes. The stable isotope sexually mature, but with B. barbus and S. cephalus correction equations were being immature. The adult ﬁsh, comprising only B. 15 15 d N d N i base barbus and S. cephalus, were all C 386 mm TP ¼ þ 2 ð1Þ 3:4 (Table 2). The ﬁsh ages were derived by scale ageing using a projecting microscope and accounting annual 13 13 d C d C i meaninv d C ¼ ð2Þ marks as per Amat-Trigo et al. (2017). For inter- corr CR inv speciﬁc data analyses, these length classes were considered separately. This was because the habitat where TP is the trophic position of the ﬁsh, d N is i i the isotopic ratio of the ﬁsh, d N is the isotopic use of these species tended to be quite different, with base the 0? ﬁshes all sampled from marginal areas of the ratio of primary consumers, 3.4 is the fractionation between trophic levels and 2 is the trophic position of river where ﬂows were minimal, the juvenile ﬁshes were generally captured from relatively shallow and the baseline organism (Post, 2002); and d C is the corr corrected carbon isotope ratio of the ﬁsh, d C is the fast-ﬂowing rifﬂe habitats, and the adult ﬁshes have uncorrected isotope ratio of the ﬁsh, d relatively large home ranges in the basin, often C is the meaninv mean invertebrate isotope ratio and CR is the exceeding 5 km (Hunt & Jones, 1974). By only inv 13 13 invertebrate carbon range (d C - d C ) (Ols- completing inter-speciﬁc analyses within these groups max min of lengths, then the data were being tested between son et al., 2009). The initial analyses using the corrected SI data tested ﬁshes of relatively similar body sizes. This meant that these analyses would be more ecologically relevant for differences in TP and C between species and Corr between different life stages of the same species using testing the hypothesis than comparing data between species of very different length ranges (Basˇic & ANOVA or Welch’s test, with the latter used where the data were normally distributed but violated the assump- Britton, 2015). Prior to analysing the stable isotope analysis of the tion of homogeneity of variance. For each life stage, the corrected SI data were then used to test the signiﬁcance ﬁshes, the stable isotope data of the macro-inverte- 15 13 of differences in their d N-d C centroids, and differ- brate samples were tested for spatial (Area 1 versus ences in the positions and overlaps of their core trophic Area 2) and temporal (2015 vs. 2016) differences. 15 13 niches. For testing differences in the d N-d C Testing used generalized linear models (GLM) due to centroids per life stage and species, the SIA data were the relatively low sample sizes that were not normally normalised by square root transformation and a resem- distributed. The GLM revealed some signiﬁcant differences (cf. Results). Thus, to enable the stable iso- blance matrix computed using Euclidean distances (Dethier et al., 2013). A PERMANOVA model was tope data of the juvenile and adult ﬁsh to be combined for use across the entire study reach, their isotopic data then ﬁtted to this distance matrix using the adonis function in the vegan package in R. This calculated the required ‘correction’ (Jackson & Britton, 2014). 15 13 Correspondingly, the d N data were converted to signiﬁcance of the differences in d N-d C centroids per group (Oksanen et al., 2007; R Core Team, 2017). trophic position (TP; Eq. 1) and the d C data were As the adonis function is similar to traditional corrected to C (Eq. 2) (Olsson et al., 2009; Jackson Corr ANOVA, it provided a pseudo F-statistic and P value & Britton, 2014). To identify the effect of this 123 264 Hydrobiologia (2018) 819:259–273 based on 999 permutations of the data (Dixon, 2003). revealed minimal differences in mean values, with Using the same method, it was then determined whether overlaps in their 95% conﬁdence limits (Chironomid: different life stages within B. barbus and S. cephalus d C: Area 1: - 31.36 ± 0.39, Area 2: 15 13 15 had signiﬁcant differences in theird N-d C centroids. 31.76 ± 0.42%; d N: Area 1: 9.88 ± 0.35, Area 2: With more than two life stages of ﬁsh being used per 9.74 ± 0.12%; Trichoptera: d C: Area 1: test, pairwise comparisons tested the signiﬁcance of 32.36 ± 0.45, Area 2: 32.25 ± 0.38%; d N: Area differences between the groups, with Bonferroni 1: 9.20 ± 0.31, Area 2: 8.86 ± 0.42%). In G. pulex, adjustment for multiple comparisons. however, some spatial and temporal differences in To compare ‘core’ trophic niche size and overlaps their SI data were apparent that, when tested in GLMs, within and between species, the isotopic niche was revealed signiﬁcant differences (d C: Wald 2 15 2 used, where the isotopic niche is an approximation of v = 12.05, P \ 0.01; d N: Wald v = 23.5, the trophic niche. It is acknowledged that the isotopic P \ 0.01; Table 1A). Pairwise comparisons of the niche varies slightly from the trophic niche due to it mean SI values from both models revealed these being inﬂuenced by factors other than diet (Jackson signiﬁcant differences were both spatial and temporal et al., 2011), such as growth and metabolic rate of for both stable isotopes (Table 1B, C). Consequently, individuals (Busst & Britton, 2017). It was calculated the use of Eqs. 1 and 2 to correct the ﬁsh SI data used using the metric ‘standard ellipse area’ (SEA), a the G. pulex SI data only (Table 1A). Prior to data bivariate measure of the distribution of individuals in correction, there were signiﬁcant differences in the isotopic space (Jackson et al., 2012). To examine the juvenile ﬁsh stable isotope data between years 13 15 size and overlap of the ‘core’ isotopic niches of each (ANOVA: d C F = 5.05, P = 0.03; d N 1,45 size group by species, ellipses were plotted that F = 11.56, P \ 0.01; Fig. 2A). However, these 1,45 enclosed 40% of the predicted data and thus the typical signiﬁcant differences were no longer apparent fol- resource use of that life stage of ﬁsh. The ellipses were lowing isotopic correction (ANOVA: d C calculated within the R package SIBER v2.1.3 (Jack- F = 0.10, P = 0.75; d NF = 1.11, P = 0.30; 1,45 1,45 son et al., 2011, 2012) and, due to some relatively small Fig. 2B). Note that in the tests, SI data from L. sample sizes, a corrected Bayesian estimate of Stan- leuciscus were not included as they were only present dard Ellipse Area (SEA ) was calculated. This was in samples in 2016. followed by a calculation utilising a Markov chain Monte Carlo simulation with 10 iterations for each Intra- and inter-speciﬁc stable isotope analysed group that provided 95% conﬁdence limits relationships (SEA ) of the isotopic niche size (Jackson et al., 2011; R Core Team, 2017). Using SEA , the extent of niche The lengths of each species were very similar across overlap (%) between species and life stages was then the 0? ﬁshes. There was greater natural variation also estimated. This was determined using the maxi- between the lengths of ﬁshes as juveniles (T. thymallus mum likelihood ﬁtted standard ellipses, with the extent were smaller than other ﬁshes) and adults (S. cephalus of the overlap between two groups thus represented by were generally smaller than B. barbus) (Table 2). The the overlap of their core niches. This was calculated only species with all life stages represented in analyses using Bayesian modelling in the SIBER package, with were B. barbus and S. cephalus. For B. barbus,C corr the denominator being the sum of non-overlapping was signiﬁcantly higher in adults than the 0? ﬁsh and area of the two ellipses (Jackson et al., 2011). juveniles (P \ 0.01), whilst TP was signiﬁcantly lower for adults versus the 0? ﬁsh (P \ 0.01) (Table 3). For S. cephalus, the 0? ﬁsh had signiﬁ- Results cantly lower C than juveniles and adults (P \ 0.01) corr and signiﬁcantly higher TP (P \ 0.01) (Table 3). Stable isotope correction for macro-invertebrate Between the species, differences in C between corr temporal and spatial differences 0? B. barbus and 0? S. cephalus were not signiﬁ- cant, but was between both of these 0? ﬁshes and Comparison of spatial differences in SI data of the 0? P. phoxinus (P \ 0.01; Table 3). The TP of 0? B. Chironomid larvae and Trichoptera spp. in 2016 barbus was signiﬁcantly higher than both S. cephalus 123 Hydrobiologia (2018) 819:259–273 265 Table 1 (A) Mean stable isotope data of Gammarus pulex 14 (A) from the two sampled areas of the River Teme in 2015 and 2016, where values are estimates (± 95% CL) from the gen- eralized linear model (GLM), and their signiﬁcance of differ- ences (as P values) according to pairwise comparisons with 13 15 Bonferroni adjustment in the GLM for (B) d C and (C) d N 12 (A) 13 15 Area Year n d C(%) d N(%) 1 2015 6 - 30.68 ± 0.58 10.78 ± 0.57 2016 4 - 29.44 ± 0.82 8.73 ± 0.71 2 2015 3 - 29.10 ± 0.82 10.22 ± 0.82 2016 6 - 29.86 ± 0.82 9.16 ± 0.82 (B) -31.0 -29.0 -27.0 -25.0 δ C(‰) d C Area 1, Area 1, Area 2, Area 2, 3.5 2015 2016 2015 2016 (B) Area 1, – 2015 3.0 Area 1, 0.05 – 2.5 Area 2, 0.01 1.00 – Area 2, 0.67 1.00 1.00 – 2.0 (C) 1.5 d N Area 1, Area 1, Area 2, Area 2, 2015 2016 2015 2016 1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 Area 1, – Ccorr Area 1, < 0.01 – Fig. 2 Uncorrected (A) and corrected (B) stable isotope data of ‘juvenile’ ﬁsh (except Leuciscus leuciscus) sampled in 2015 Area 2, 1.00 0.04 – (clear circle) and 2016 (ﬁlled circle) Area 2, < 0.01 1.00 0.43 – differences in C and TP between adult B. barbus corr and S. cephalus (P [ 0.05; Table 3). Sample number (n) represents the number of samples analysed, 15 13 where one sample = 3–6 individual gammarids Differences in d N-d C centroids Values in bold are signiﬁcant at P B 0.05 The overall test for differences in the positions of the 15 13 d N-d C centroids of the 0? ﬁshes was signiﬁcant and P. phoxinus (P \ 0.01), although it was not (F = 15.72, P \ 0.01); pairwise comparisons indi- 2,58 signiﬁcantly different between S. cephalus and P. cated the signiﬁcant differences were between P. phoxinus (P [ 0.05; Table 3). In juveniles, C of B. corr phoxinus and the other ﬁshes (P \ 0.01 in both cases), barbus was signiﬁcantly lower than all other ﬁshes but with differences between B. barbus and S. (P \ 0.01; Table 3). The TP of juvenile B. barbus was cephalus being not signiﬁcant (P [ 0.05; Table 4). signiﬁcantly higher than T. thymallus, signiﬁcantly For the juvenile ﬁshes, the overall model was signif- lower than L. leuciscus, but not signiﬁcantly different icant (F = 18.41, P \ 0.01), with signiﬁcant dif- 3,76 to S. cephalus (Table 3). There were no signiﬁcant ferences between all species (P \ 0.01; Table 4). This TP δ N (‰) 266 Hydrobiologia (2018) 819:259–273 Table 2 The number, ﬁsh length ranges, mean lengths (95% CL) and measurement type (SL standard length, FL fork length) of each life stage of ﬁsh analysed for their stable isotopes across the two sampling areas Species n Length range (mm) Mean length (mm) (± 95% CL) Length measurement 0? B. barbus 30 18–34 25.2 ± 1.8 SL 0? S. cephalus 15 17–36 27.3 ± 2.4 SL 0? P. phoxinus 16 17–38 27.3 ± 2.8 SL Juvenile B. barbus 16 105–231 158 ± 15 FL Juvenile S. cephalus 16 112–207 153 ± 11 FL Juvenile L. leuciscus 30 102–214 167 ± 11 FL Juvenile T. thymallus 15 86–205 122 ± 16 FL Adult B. barbus 21 540–690 584 ± 17 FL Adult S. cephalus 21 386–570 466 ± 22 FL Table 3 Outputs of Species Test Testing df FP ANOVA/Welch’s test of corrected carbon (C ) and corr C corr trophic position (TP) for Barbus barbus ANOVA Life stage 2,64 32.76 \0.01 comparisons between life Squalius cephalus Welch’s Life stage 2,31 17.66 \0.01 stages and species for Barbus barbus and Squalius 0? ANOVA Species 3,60 10.01 \0.01 cephalus Juvenile Welch’s Species 3,36 16.61 \0.01 Adult ANOVA Species 1,40 2.09 0.16 Species/ life stage Test Testing df FP TP Barbus barbus ANOVA Life stage 2,64 47.17 \0.01 Squalius cephalus ANOVA Life stage 2,49 6.52 \0.01 0? ANOVA Species 3,60 12.36 \0.01 Juvenile Welch’s Species 3,30 60.53 \0.01 Note data for all sites and Adult ANOVA Species 1,40 0.02 0.90 years are combined was in contrast to the adult B. barbus and S. cephalus, Core isotopic niches (standard ellipse areas) which were not signiﬁcantly different (F = 1.77, 1,41 The 95% conﬁdence intervals of the core isotopic P = 0.18). 15 13 The model testing differences in d N-d C cen- niches (as standard ellipse areas) at each life stage suggested they were similar in size between the troids between the different life stages of B. barbus was signiﬁcant (F = 28.89, P \ 0.01), with pair- species (Table 5). In general, the core isotopic niches 2,64 of the 0? ﬁshes had low overlap (maximum 7% wise comparisons indicating the signiﬁcant differ- ences were between adults and the other life stages between B. barbus and S. cephalus), the juvenile ﬁshes (P \ 0.01; Table 4). Whilst the overall model was had no niche overlap, but in adult B. barbus and S. also signiﬁcant in S. cephalus (F = 17.31, cephalus, their niches overlapped by 55% (Fig. 3). 2,49 P \ 0.01), pairwise comparisons indicated the signif- Within B. barbus, there was no overlap in their core icant differences were only between the 0? ﬁshes and niches between the 0?, juvenile and adult ﬁsh (Fig. 3). the other life stages (P \ 0.01; Table 4). In S. cephalus, there was no niche overlap between 0? and juveniles, with this increasing to 2% between 123 Hydrobiologia (2018) 819:259–273 267 15 13 Table 4 The signiﬁcance of differences in d N-d C centroids between the different life stages of the ﬁshes, as represented by P values (with Bonferroni adjustment for multiple comparisons) derived in PERMANOVA 0? B. 0? S. 0? P. J B. J S. J L. J T. A B. A S. barbus cephalus phoxinus barbus cephalus leuciscus thymallus barbus cephalus 0? B. barbus - 0? S. cephalus 0.74 – 0? P. phoxinus \ 0.01* \ 0.01* – J B. barbus 0.22 – – – J S. cephalus – \ 0.01* – 0.01* – J L. leuciscus – – – 0.02* 0.01* – J T. thymallus – – – 0.01* 0.01* 0.02* – A B. barbus \ 0.01* – – \ 0.01* – – – – A S. cephalus – \ 0.01* – – 0.62 – – 0.18 – Table 5 Mean Lifestage and species C TP SEAc (± 95% CL) corr stable isotope data (± 95% CL) and standard ellipse 0? B. barbus 0.06 ± 0.38 3.30 ± 0.08 0.77 ± 0.28 areas (as SEAc; ± 95% CI 0? S. cephalus 0.49 ± 0.65 2.92 ± 0.13 0.96 ± 0.51 SEAb) for the sampled 0? P. phoxinus - 1.49 ± 0.49 3.02 ± 0.11 0.73 ± 0.36 ﬁshes in the study river and across the three life stages Juvenile B. barbus 0.39 ± 0.58 2.70 ± 0.13 0.54 ± 0.28 (0?, juvenile and adult) Juvenile S. cephalus 2.57 ± 0.75 2.63 ± 0.14 0.59 ± 0.30 Juvenile L. leuciscus 1.42 ± 0.39 3.03 ± 0.05 0.44 ± 0.16 Juvenile T. thymallus 1.57 ± 0.20 2.13 ± 0.14 0.28 ± 0.14 Adult B. barbus 2.52 ± 0.40 2.62 ± 0.15 1.34 ± 0.58 Adult S. cephalus 3.22 ± 0.45 2.61 ± 0.15 1.89 ± 0.83 0 ? and adults, and then the juvenile niche sitting and core isotopic niches. The centroids were calcu- entirely within the adult niche (Fig. 4). lated using all SI data per life stage and species, whereas cores niches are based on a predicted 40% of the SI data to indicate typical resource use (Jackson Discussion et al., 2011, 2012). There were some consistent results from these analyses that aligned with Hypothesis 1, Hypothesis 1 tested whether there were consistent especially in the juvenile ﬁshes where there were 15 13 inter-speciﬁc patterns of trophic partitioning between signiﬁcant differences in d N-d C centroids between B. barbus and the other ﬁshes. It was formulated due to all species and no overlaps in their core niches. In the these patterns of niche partitioning being evident 0? ﬁshes, there was less consistency in the results of between the ﬁshes in the B. barbus indigenous range both analyses, with isotopic niches showing low inter- 15 13 (Basic & Britton, 2016). An alternative to this speciﬁc overlap, but with d N-d C centroids show- hypothesis would be the ﬁshes having high niche ing signiﬁcant differences only between P. phoxinus overlap, as has been suggested between invasive and and the other ﬁshes. Whilst there was poor alignment endemic Barbus ﬁshes in Italian rivers, where it of the results in the adult ﬁshes with Hypothesis 1, both 15 13 appears to have resulted in the competitive displace- analyses provided consistent results; the d N-d C ment of the endemics (Carosi et al., 2017). Hypothesis centroids of the adult B. barbus and S. cephalus were 15 13 1 was tested using two analyses, d N-d C centroids not signiﬁcantly different and their core niches had 123 268 Hydrobiologia (2018) 819:259–273 This pattern was evident in rivers that had been (a) 4.0 stocked with hatchery-reared B. barbus at sizes below 3.5 ˇ ´ 250 mm and remained evident in adult ﬁshes (Basic 3.0 and Britton, 2016). In contrast, Gutmann Roberts et al. 2.5 (2017) revealed that in the lower reaches of the study 2.0 river, there was high overlap in the core isotopic 1.5 niches of adult B. barbus and S. cephalus, primarily 1.0 the result of individual ﬁshes specialising in the -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 consumption of pelletized marine ﬁshmeal utilised by (b) 4.0 anglers. This niche overlap was also evident in the 3.5 adult ﬁshes here, where inter-speciﬁc niche differ- 3.0 ences were only signiﬁcant in the 0? and juvenile ﬁshes. In B. barbus, there were strong and signiﬁcant 2.5 patterns in niche partitioning between their different 2.0 life stages, suggesting considerable ontogenetic shifts 1.5 in their diet that resulted in their population having a 1.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 relatively large core isotopic niche that was composed of at least three distinct sub-sets. This was consistent (c) 4.0 15 13 with Hypothesis 2, although the d N-d C centroids 3.5 suggested differences were not signiﬁcant between the 3.0 0? and juvenile ﬁsh. In contrast, the isotopic niches of 2.5 S. cephalus were more similar over their three studied 2.0 life stages, with only the niche of the 0? ﬁsh being 1.5 distinct from the other life stages, with the juvenile and 1.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 adult niches overlapping completely, contrary to Ccorr Hypothesis 2. Fig. 3 Corrected Carbon (C ) versus trophic position (TP) Stable isotope data of 0? ﬁshes can be confounded corr for Barbus barbus (ﬁlled triangle), Squalius cephalus (ﬁlled by issues of their data still showing a strong parental square), Leuciscus leuciscus (delta) and Thymallus thymallus signal. For example, in anadromous brown trout Salmo (open square), and the positions of their core isotopic niches (as trutta, newly emerged fry retained a strong parental, SEA ), where solid line: B. barbus, small dashed line: S. marine-based isotopic signal that enabled their differ- cephalus, long dashed line: L. leuciscus, and dash/dot line: T. thymallus.(a)0? ﬁshes; (b) juvenile ﬁshes; and (c) adult ﬁshes entiation from fry produced from non-anadromous parents, but this difference was much reduced after relatively high overlap (55%). Across all life stages four months of feeding in freshwater (Briers et al., and analyses, there was no evidence to support the 2013). In 0? smallmouth bass Micropterus dolomieu, alternative hypothesis that invasion by B. barbus had post-hatch embryos had elevated d N values that resulted in competitive displacement of native ﬁshes, were associated with their parental origin, but these as suggested by Carosi et al. (2017) for endemic values subsequently decreased rapidly due to their Barbus. It is, however, acknowledged that this was not exogenous feeding during their metamorphosis from tested implicitly here, given the absence of data from larvae into juveniles (Vander Zanden et al., 1998). the pre-invaded period or from sites with B. barbus Here, the 0? ﬁshes utilised were all of lengths above absent. 17 mm, were all fully formed juveniles rather than The pattern of isotopic niche partitioning between larvae, and were up to 10 weeks old. Their stable iso- B. barbus and other ﬁshes was thus consistent with a tope data were very distinct from those of the adult number of isotopic studies completed on populations ﬁshes; in terms of uncorrected data, the 0? ﬁshes were in their indigenous range (Basˇic´ et al., 2015;Basˇic´ and depleted in d Cbyupto8% compared to adult Britton, 2015, 2016). These studies all suggested that conspeciﬁcs. Consequently, the strong patterns of core B. barbus and S. cephalus have distinct core isotopic isotopic niche partitioning detected in these 0? ﬁshes niches, with minimal inter-speciﬁc resource sharing. were interpreted as resulting from their dietary TP Hydrobiologia (2018) 819:259–273 269 Fig. 4 Intra-speciﬁc (a) 4.0 comparisons of Corrected Carbon (C ) versus corr trophic position (TP) and 3.6 positions of core isotopic niches (as SEA ) for (a) Barbus barbus and (b) Squalius cephalus, and 3.2 where ﬁlled diamond, solid line: 0? ﬁsh; open circles, dashed line: juvenile ﬁsh; 2.8 ﬁlled triangle, dotted line: adult ﬁsh 2.4 2.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 4.0 (b) 3.6 3.2 2.8 2.4 2.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 Ccorr differences formed by their exogenous feeding within all samples from all areas sampled (in contrast to other the river, rather than being a legacy of their parental macro-invertebrate taxa), there were signiﬁcant dif- isotopes. ferences in their isotopic data between years and sites In this study, the ﬁsh SI data were ‘corrected’ to (whereas data from other macro-invertebrate taxa enable data to be combined across two sampling years revealed similar values, at least spatially), and liter- and relatively long stretch of river. The use of standard ature suggests Gammarus spp. are an important prey equations to ‘correct’ SI data in this manner is well item for many ﬁshes (e.g. MacNeil et al., 1999). established (e.g. Olsson et al., 2009; Jackson & Correction of the juvenile ﬁsh SI data using the G. Britton, 2014). However, this correction relies upon pulex data was also demonstrated as removing signif- an adequate description of the SI data of the ﬁsh prey icant temporal differences. To enable the SI data of the resources. Here, correction focused on use of G. pulex 0? ﬁshes to be comparable to the other ﬁshes, they SI data. The rationale for this was they were present in were also corrected using the G. pulex data, although TP 270 Hydrobiologia (2018) 819:259–273 stomach contents analyses (SCA) for the 0? ﬁshes barbus presence, it is acknowledged that this is had suggested that these amphipods were a minor prey speculative given that isotopic niche sizes of the item (Gutmann Roberts & Britton, 2018). However, native ﬁshes were unable to be measured in B. barbus individuals of [ 20 mm standard length in all the absence. Notwithstanding, the inter-speciﬁc isotopic 0? ﬁsh species analysed had infections of the intesti- niche partitioning evident in the study suggests that nal parasite Pomphorhynchus tereticollis (C. Gutmann despite their similar ecological guilds and sharing Roberts, unpublished data). This parasite has gam- similar habitats (especially the 0? ﬁshes), there were marids as its intermediate host (Kaldonski et al., sufﬁcient differences between the ﬁshes in their 2008), suggesting G. pulex might have been consumed functional traits and/or habitat utilisation to enable in greater proportions in the 0? ﬁshes than suggested substantial differentiation in their core isotopic niches by SCA. Moreover, in the River Severn basin, to occur (Robinson et al., 1993; Borcherding et al., Pomphorhynchus spp. has been reported as prevalent 2013; Negus & Hoffman, 2013). in all the ﬁshes studied here (Brown 1984), suggesting On one hand, the results here could suggest that the gammarids are a common and important prey item of ecological impacts of invasive B. barbus are relatively ﬁsh in the river. Therefore, although it is acknowl- minor in the river, as there was little evidence to edged that the macro-invertebrate baseline SI data suggest there was high diet similarity in the ﬁshes at used to correct the data here could have utilised a their most abundant life stages (0? and juveniles). wider range of taxa, especially in 2015, it is strongly This inference is supported by other recent studies on argued that the use of G. pulex SI data to correct the native B. barbus that have revealed strong patterns of ﬁsh SI data was justiﬁed and appropriate. inter-speciﬁc core isotopic niche partitioning (e.g. Following introductions of non-native ﬁshes, Basic & Britton, 2015, 2016). However, these studies adverse ecological impacts often develop through were all limited to assessing trophic interactions via increased inter-speciﬁc competition for food resources stable isotope analysis, with studies suggesting that between invasive and sympatric native ﬁshes (Gozlan when compared with other dietary analysis methods, et al., 2010; Cucherousset et al., 2012). Given the such as stomach contents analyses, different results relatively similar size ranges of the invasive B. barbus can occur, resulting from differences between items with other cyprinid ﬁshes across the different life ingested (stomach contents) and assimilated (SIA) stages (albeit with some inter-speciﬁc length differ- (e.g. Locke et al., 2013). Consequently, some caution ences within life stages), this suggests there was is necessary if isotopic niche overlaps are to be used to considerable potential for inter-speciﬁc competitive infer the strength of competitive interactions. More- interactions, especially given the ﬁshes were from over, invasive B. barbus can potentially result in other relatively similar functional guilds (Basic & Britton, ecological concerns, such as causing habitat alter- 2016). However, the lack of overlap in the isotopic ations, given recent work has demonstrated that in niches of the 0? and juvenile ﬁshes—the life stages their native range, B. barbus act as ‘zoogeomorphic when their abundances tend to be highest—suggested agents’ (Pledger et al., 2014, 2016). This is where their low dietary overlap, with only the isotopic niches of benthic foraging activities can reduce bed material the adult ﬁshes indicating some dietary overlap. stability, increase bedload transport, and impact Schulze et al. (2012) suggested that species within microtopographic roughness and sediment structure the same ecological guild can only coexist when they (Pledger et al., 2014, 2016). This benthic foraging respond differently to limited resource availability could then also impact upon aspects of the macro- with, for example, specialised species only persisting invertebrate communities, such as decreased abun- if their competitors are generalists. Evidence in dance via predation or reduced species richness via literature supports this, with reduced trophic niche disturbance. However, these aspects were beyond the sizes in many co-existing ﬁshes when compared to scope of this study and so further research is required allopatry (Bolnick et al., 2010; Tran et al., 2015). In to provide increased understandings of how B. barbus the study river, however, even where the isotopic invasions affect macro-invertebrate communities and niches of the ﬁshes were partitioned, the niches were sediment structure. Finally, it should also be noted that similarly sized. Although this suggests there had not the study river was low in ﬁsh species richness, with a been any niche constriction in the native ﬁshes in B. cyprinid ﬁsh community comprising only four species 123 Hydrobiologia (2018) 819:259–273 271 assessments of trophic relationships in riverine ﬁsh com- (including B. barbus), with other ﬁsh taxa being very munities. Journal of Applied Ichthyology 31: 296–300. limited in diversity and abundance in the study reach, Basˇic´, T. & J. R. Britton, 2016. Characterizing the trophic niches and with no other Barbus species present. Conse- of stocked and resident cyprinid ﬁshes: consistency in quently, if B. barbus are introduced into a river with partitioning over time, space and body sizes. Ecology and Evolution 6: 5093–5104. considerably higher native ﬁsh species richness, Basˇic´, T., J. R. Britton, M. C. Jackson, P. Reading & J. Grey, irrespective of the presence of any other Barbus 2015. Angling baits and invasive crayﬁsh as important species, there is the possibility that there will be a trophic subsidies for a large cyprinid ﬁsh. Aquatic Sciences greater probability of higher niche overlaps within 77: 153–160. Basˇic´, T., J. R. Britton, S. P. Rice & A. G. Pledger, 2017. species in the ﬁsh community and thus higher potential Impacts of gravel jetting on the composition of ﬁsh for ecological impacts to result. It is therefore spawning substrates: Implications for river restoration and recommended that such introductions proceed only ﬁsheries management. Ecological Engineering 107: 71–81. with caution and full risk assessment (Roy et al., Bolnick, D. I., T. Ingram, W. E. Stutz, L. K. Snowberg, O. L. Lau & J. S. Paull, 2010. Ecological release from inter- 2018). speciﬁc competition leads to decoupled changes in popu- In summary, across three life stages of invasive B. lation and individual niche width. Proceedings of the Royal barbus, there were some strong patterns of isotopic Society—Biological Sciences 277: 1789–1797. niche partitioning with native ﬁshes, with this parti- Borcherding, J., M. Dolina, L. Heermann, P. Knutzen, S. Kru- ger, S. Matern, R. van Treeck & S. Gertsen, 2013. Feeding tioning initially evident between some ﬁshes during and niche differentiation in three invasive gobies in the their ﬁrst year of life that became strongly apparent at Lower Rhine, Germany. Limnologica 43: 49–58. juvenile life stages. These invasive B. barbus thus Briers, R. A., J. O. Waterman, K. Galt & R. N. Campbell, 2013. integrated into this riverine food web via exploiting Population differentiation and temporal changes of car- otenoid pigments and stable isotope ratios in the offspring different food resources to the native ﬁshes that of anadromous and non-anadromous trout Salmo trutta. facilitated their co-existence. Ecology of Freshwater Fish 22: 137–144. Britton, J. R. & J. Pegg, 2011. Ecology of European barbel Acknowledgements We thank the Environment Agency for Barbus barbus: implications for river, ﬁshery, and con- their assistance with collection of scales from the juvenile ﬁshes. servation management. Reviews in Fisheries Science 19: CGR was supported by a studentship funded by the Severn 321–330. Rivers Trust and Bournemouth University, with additional Britton, J. R., I. G. Cowx & G. Peirson, 2004. Sources of error in support from the Barbel Society. the ageing of stocked cyprinids. Fisheries Management and Ecology 11(6): 415–417. Open Access This article is distributed under the terms of the Brown, A.F., 1984. The ecology of Pomphorhynchus laevis Creative Commons Attribution 4.0 International License (http:// (Muller, 1776) (Acanthocephala) in the River Severn, creativecommons.org/licenses/by/4.0/), which permits unre- England. Ph.D. thesis, University of Liverpool. stricted use, distribution, and reproduction in any medium, Busst, G. M. & J. R. Britton, 2016. High variability in provided you give appropriate credit to the original stable isotope diet–tissue discrimination factors of two author(s) and the source, provide a link to the Creative Com- omnivorous freshwater ﬁshes in controlled ex situ condi- mons license, and indicate if changes were made. tions. Journal of Experimental Biology 219: 1060–1068. Busst, G. M. & J. R. Britton, 2017. Tissue-speciﬁc turnover rates of the nitrogen stable isotope as functions of time and growth in a cyprinid ﬁsh. Hydrobiologia. https://doi.org/ 10.1007/s10750-017-3276-2. References Carosi, A., L. Ghetti, G. La Porta & M. Lorenzoni, 2017. Eco- logical effects of the European barbel Barbus barbus (L., Amat-Trigo, F., C. Gutmann Roberts & J. R. Britton, 2017. 1758) (Cyprinidae) invasion on native barbel populations Spatial variability in the growth of invasive European in the Tiber River basin (Italy). The European Zoological barbel Barbus barbus in the River Severn basin, revealed Journal 84: 420–435. using anglers as citizen scientists. Knowledge and Man- CEH (2018) 54029—Teme at Knghtsford Bridge. http://nrfa. agement of Aquatic Ecosystems 418: 17–23. ceh.ac.uk/data/station/meanﬂow/54029. Accessed on 15 Antognazza, C. M., D. Andreou, S. Zaccara & J. R. Britton, Feb, 18. 2016. Loss of genetic integrity and biological invasions Copp, G. H., J. R. Britton, Z. Guo, V. R. Edmonds-Brown, J. result from stocking and introductions of Barbus barbus: Pegg, L. Vilizzi & P. I. Davison, 2017. Trophic conse- insights from rivers in England. Ecology and Evolution 6: quences of non-native pumpkinseed Lepomis gibbosus for 1280–1292. native pond ﬁshes. Biological Invasions 19: 25–41. Basˇic´, T. & J. R. Britton, 2015. Utility of ﬁsh scales from stock Cucherousset, J. & J. D. Olden, 2011. Ecological impacts of assessment surveys in stable isotope analysis for initial nonnative freshwater ﬁshes. Fisheries 36: 215–230. 123 272 Hydrobiologia (2018) 819:259–273 Cucherousset, J., S. Bouletreau, A. Martino, J. M. Roussel & F. Negus, M. T. & J. C. Hoffman, 2013. Habitat and diet differ- Santoul, 2012. Using stable isotope analyses to determine entiation by two strains of rainbow trout in Lake Superior the ecological effects of non-native ﬁshes. Fisheries based on archival tags, stable isotopes, and bioenergetics. Management and Ecology 19: 111–119. Journal of Great Lakes Research 39: 578–590. Dethier, M. N., E. Sosik, A. W. Galloway, D. O. Duggins & C. Oksanen, J., R. Kindt, P. Legendre, B. O’Hara, M. H. H. Ste- A. Simenstad, 2013. Addressing assumptions: variation in vens, M. J. Oksanen & M. Suggests, 2007. The vegan stable isotopes and fatty acids of marine macrophytes can package. Community ecology package 10: 631–637. confound conclusions of food web studies. Marine Ecology Olsson, K., P. Stenroth, P. E. R. Nystro¨m & W. Grane´li, 2009. Progress Series 478: 1–14. Invasions and niche width: does niche width of an intro- Dixon, P., 2003. VEGAN, a package of R functions for com- duced crayﬁsh differ from a native crayﬁsh? Freshwater munity ecology. Journal of Vegetation Science 14(6): Biology 54: 1731–1740. 927–930. Persson, L., J. Andersson, E. Wahlstrom & P. Eklov, 1996. Size- Gozlan, R. E., J. R. Britton, I. Cowx & G. H. Copp, 2010. speciﬁc interactions in lake systems: predator gape limi- Current knowledge on non-native freshwater ﬁsh intro- tation and prey growth rate and mortality. Ecology 77: ductions. Journal of Fish Biology 76: 751–786. 900–911. Gutmann Roberts, C., T. Basˇic´, F. Amat Trigo & J. R. Britton, Pledger, A. G., S. P. Rice & J. Millett, 2014. Reduced bed 2017. Trophic consequences for riverine cyprinid ﬁshes of material stability and increased bedload transport caused angler subsidies based on marine-derived nutrients. by foraging ﬁsh: a ﬂume study with juvenile barbel (Barbus Freshwater Biology 62: 894–905. barbus). Earth Surface Processes and Landforms 39: Gutmann Roberts, C. & Britton, J.R. 2018. Quantifying trophic 1500–1513. interactions and niche sizes of juvenile ﬁshes in an invaded Pledger, A. G., S. P. Rice & J. Millett, 2016. Bed disturbance via riverine cyprinid ﬁsh community. Ecology of Freshwater foraging ﬁsh increases bedload transport during subsequent Fish (in press). high ﬂows and is controlled by ﬁsh size and species. Hunt, P. C. & J. W. Jones, 1974. A population study of Barbus Geomorphology 253: 83–93. barbus (L.) in the River Severn. England. Journal of Fish Post, D. M., 2002. Using stable isotopes to estimate trophic Biology 6: 255–267. position: models, methods, and assumptions. Ecology 83: Hutchinson, J. J. & C. N. Trueman, 2006. Stable isotope anal- 703–718. yses of collagen in ﬁsh scales: Limitations set by scale Post, D. M., C. A. Layman, D. A. Arrington, G. Takimoto, J. architecture. Journal of Fish Biology 69: 1874–1880. Quattrochi & C. G. Montana, 2007. Getting to the fat of the Jackson, M. C. & J. R. Britton, 2014. Divergence in the trophic matter: models, methods and assumptions for dealing with niche of sympatric freshwater invaders. Biological Inva- lipids in stable isotope analyses. Oecologia 152: 179–189. sions 16: 1095–1103. R Core Team, 2017. R: A language and environment for statis- Jackson, A. L., R. Inger, A. C. Parnell & S. Bearhop, 2011. tical computing. R Foundation for Statistical Computing, Comparing isotopic niche widths among and within com- Vienna. http://www.R-project.org/. munities: SIBER–Stable Isotope Bayesian Ellipses in R. Robinson, B. W., D. S. Wilson, A. S. Margosian & P. T. Lotito, Journal of Animal Ecology 80: 595–602. 1993. Ecological and morphological differentiation of Jackson, M. C., I. Donohue, A. L. Jackson, J. R. Britton, D. pumpkinseed sunﬁsh in lakes without bluegill sunﬁsh. M. Harper & J. Grey, 2012. Population-level metrics of Evolutionary Ecology 7: 451–464. trophic structure based on stable isotopes and their appli- Roy, H. E., W. Rabitsch, R. Scalera, A. Stewart, B. Gallardo, P. cation to invasion ecology. PLoS ONE 7: e31757. Genovesi, F. Essl, T. Adriaens, S. Bacher, O. Booy & E. Kaldonski, N., M. J. Perrot-Minnot, S. Motreuil & F. Ce´zilly, Branquart, 2018. Developing a framework of minimum 2008. Infection with acanthocephalans increases the vul- standards for the risk assessment of alien species. Journal nerability of Gammarus pulex (Crustacea, Amphipoda) to of Applied Ecology 55: 526–538. non-host invertebrate predators. Parasitology 135: Schulze, T., H. Do¨rner, U. Baade & F. Ho¨lker, 2012. Dietary 627–632. niche partitioning in a piscivorous ﬁsh guild in response to Locke, S. A., G. Bulte, M. R. Forbes & D. J. Marcogliese, 2013. stocking of an additional competitor–The role of diet spe- Estimating diet in individual pumpkinseed sunﬁsh Lepomis cialisation. Limnologica-Ecology and Management of gibbosus using stomach contents, stable isotopes and par- Inland Waters 42(1): 56–64. asites. Journal of Fish Biology 82: 522–537. Sinnatamby, R. N., J. E. Bowman, J. B. Dempson & M. Power, Mann, R. H. K., 1974. Observations on the age, growth, repro- 2007. An assessment of de-calciﬁcation procedures for 13 15 duction and food of the dace, Leuciscus leuciscus (L.), in d C and d N analysis of yellow perch, walleye and two rivers in southern England. Journal of Fish Biology 6: Atlantic salmon scales. Journal of Fish Biology 70: 237–253. 1630–1635. MacNeil, C., J. T. Dick & R. W. Elwood, 1999. The dynamics of Skinner, M. M., A. A. Martin & B. C. Moore, 2016. Is lipid predation on Gammarus spp. (Crustacea: Amphipoda). correction necessary in the stable isotope analysis of ﬁsh Biological Reviews 74: 375–395. tissues? Rapid Communications in Mass Spectrometry 30: Meraner, A., A. Venturi, G. F. Ficetola, S. Rossi, A. Candiotto & 881–889. A. Gandolﬁ, 2013. Massive invasion of exotic Barbus Tran, T. N. Q., M. C. Jackson, D. Sheath, H. Verreycken & J. barbus and introgressive hybridization with endemic R. Britton, 2015. Patterns of trophic niche divergence Barbus plebejus in Northern Italy: where, how and why? between invasive and native ﬁshes in wild communities are Molecular Ecology 22: 5295–5312. 123 Hydrobiologia (2018) 819:259–273 273 predictable from mesocosm studies. Journal of Animal productivity and CO concentration. Hydrobiologia 644: Ecology 84: 1071–1080. 245–259. Trigo, F. A., C. Gutmann Roberts & R. J. Britton, 2017. Spatial Wheeler, A. & D. R. Jordan, 1990. The status of the barbel, variability in the growth of invasive European barbel Barbus barbus (L.) (Teleostie, Cyprinidae) in the United Barbus barbus in the River Severn basin, revealed using Kingdom. Jouranl of Fish Biology 37: 393–399. anglers as citizen scientists. Knowledge & Management of Woodcock, S. H. & B. D. Walther, 2014. Trace elements and Aquatic Ecosystems 418(17): 1–6. stable isotopes in Atlantic tarpon scales reveal movements Vander Zanden, M. J., M. Hulshof, M. S. Ridgway & J. across estuarine gradients. Fisheries Research 153: 9–17. B. Rasmussen, 1998. Application of stable isotope tech- Zaccara, S., C. M. Antognazza, L. Buonerba, J. R. Britton & G. niques to trophic studies of age-0 smallmouth bass. Crosa, 2014. Human-mediated contact zone between Transactions of the American Fisheries Society 127: endemic and invasive Barbus species (Osteichthyes: 729–739. Cyprinidae) in a regulated lowland river: genetic infer- Ventura, M. & E. Jeppesen, 2010. Evaluating the need for acid ences and conservation implications. Italian Journal of 13 15 treatment prior to d C and d N analysis of freshwater ﬁsh Zoology 81: 571–583. scales: effects of varying scale mineral content, lake
Hydrobiologia – Springer Journals
Published: May 22, 2018
It’s your single place to instantly
discover and read the research
that matters to you.
Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.
All for just $49/month
Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly
Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.
Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.
Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.
All the latest content is available, no embargo periods.
“Hi guys, I cannot tell you how much I love this resource. Incredible. I really believe you've hit the nail on the head with this site in regards to solving the research-purchase issue.”Daniel C.
“Whoa! It’s like Spotify but for academic articles.”@Phil_Robichaud
“I must say, @deepdyve is a fabulous solution to the independent researcher's problem of #access to #information.”@deepthiw
“My last article couldn't be possible without the platform @deepdyve that makes journal papers cheaper.”@JoseServera