AUTHOR’S PAGE PROOFS: NOT FOR CIRCULATION
CSIRO PUBLISHING Marine and Freshwater Research, 2018, 69, 1–7 Short Communication https://doi.org/10.1071/MF17092
Inter-population variability in growth and reproduction of invasive bleak Alburnus alburnus (Linnaeus, 1758) across the Iberian Peninsula
D. LatorreA, G. Maso´ A,J, A. HinckleyB, D. Verdiell-CubedoA, A. S. TarkanC,K, A. Vila-GispertA, G. H. CoppD, J. CucheroussetE, E. da SilvaF, C. Ferna´ndez-DelgadoG, E. Garcı´a-BerthouA, R. MirandaH, F. J. Oliva-PaternaI, A. Ruiz-NavarroI, J. M. SerranoB and D. AlmeidaA,L,M Y
AGRECO, Institute of Aquatic Ecology, University of Girona, C/ M. Aure`lia Capmany 69, E-17003 Girona, Spain. BDepartment of Ecology, Complutense University of Madrid, C/ Jose´ Antonio Novais 12, E-28040 Madrid, Spain. CFaculty of Fisheries, Mug˘la Sıtkı Koc¸man University, TR-48000 Mentese, Mug˘la, Turkey. DSalmon & Freshwater Team, Cefas, Pakefield Road, Lowestoft, Suffolk, NR33 OHT, UK. ECNRS, Universite´ Toulouse III Paul Sabatier, ENFA, UMR 5174 EDB (Laboratoire E´volution & Diversite´ Biologique), 118 Rue Narbonne, F-31062 Toulouse, France. FDepartment of Anatomy, Cell Biology & Zoology, University of Extremadura, Avenida De Elvas, s/n, E-06006 Badajoz, Spain. GDepartment of Zoology, University of Co´rdoba, Campus de Rabanales, Edificio Charles Darwin, E-14014 Co´rdoba, Spain. HDepartment of Environmental Biology, School of Sciences, University of Navarra, C/ Irunlarrea 1, E-31080 Pamplona, Spain. IDepartment of Zoology and Physical Anthropology, University of Murcia, Avenida Teniente Flomesta 5, E-30100 Murcia, Spain. AQ1 JDepartment of Biodiversity and Ecological Restoration, Instituto Pirenaico de Ecologı´a (IPE-CSIC), Avenida Nuestra Sen˜ora de la Victoria 16, E-22700 Jaca, Spain. KDepartment of Ecology and Vertebrate Zoology, Faculty of Biology and Environmental AQ2 Protection, 12/16 Banacha Str., PL-90-237, University of Ło´dz´, Poland. LDepartamento de Ciencias Me´dicas Ba´sicas, Facultad de Medicina, Universidad San Pablo-CEU, CEU Universities, Campus de Monteprı´ncipe, E-28925 Alcorco´n, Madrid, Spain. MCorresponding author. Email: [email protected]
Abstract. The native European freshwater cyprinid fish, common bleak Alburnus alburnus (Linnaeus, 1758), is an invasive non-native fish in the Iberian Peninsula, where it threatens the valuable endemic fish fauna. Despite the bleak’s invasive status, there is a dearth of studies on the biological traits of non-native populations in Europe’s Mediterranean region, and the presentROOF study aimed to compare bleak growth and reproductive ONL traits across the principal rivers of Iberia 5 with those from a native ‘reference’ bleak population in France. Non-native bleak from the River Tagus had the highest back-calculated total lengths (TLs), growth rate and body condition, but the lowest reproductive investment and smallest egg size of all studied populations. Whereas, these latter two traits were the highest in the River Ebro, where fecundity was the lowest for all Iberian rivers. The youngest age and the smallest TL at maturity were observed in the River Ebro. The River SeguraP had the lowest back-calculated TLs and growth rate, the highest fecundity and proportion of females, and the 10 oldest age at maturity. Population traits of bleak in the River Saoˆne (the native ‘reference’ population) were generally different from those in Iberian rivers, except for the River Ebro. This wide inter-population variability contributes to the species’ successful establishment in novel habitats.
Additional keywords: back-calculation, body condition, freshwater fish, non-native species, sex ratio.
Received 6 April 2017, accepted 5 December 2017, published online xx xxxxx xxxx
Journal compilation Ó CSIRO 2018 www.publish.csiro.au/journals/mfr 2 Marine and Freshwater Research D. Latorre et al.
Introduction as well as with regard to a native bleak population of the River Invasive fishes can be major drivers of biodiversity loss in Saoˆne (eastern France), which was used as the ‘reference’ freshwater ecosystems at the global scale (Mooney and Hobbs population because France has historically been the ‘donor 2000). The potentially disruptive effects of non-native fishes are region’ for non-native fish introductions to Iberia along the 5 of particular conservation concern in the Iberian Peninsula, so-called ‘Perpignan–Barcelona corridor’ (see Clavero and 5 where endemism in freshwater fishes is high (.50% of native Garcı´a-Berthou 2006). species) and many of the endemic species are seriously threat- ened (Reyjol et al. 2007). Indeed, introduced fishes currently represent .30% of the total fish species and this continues to Materials and methods 10 increase, which renders the Iberian Peninsula a hotspot for fish Each Iberian sampling site was surveyed for two consecutive invasions (Leunda 2010; Almeida et al. 2013). days by different methods (see details below). In particular, The common bleak Alburnus alburnus (Linnaeus, 1758) is sampling dates were on 25–28 May for Ebro (two sites), 28–31 10 native to most of Europe, from the Pyrenees Mountains to the May for Tagus (two sites), 29 May–1 June for GuadianaY (two Ural Mountains, and including the south-eastern part of Great sites), 22–25 May for Segura (two sites) and 7–10 May for 15 Britain, where this cyprinid species inhabits lakes and medium- Guadalquivir (two sites) and 5–6 June for Saoˆne (one site). The to-large rivers (Keith et al. 2011). In the Iberian Peninsula, following five main Iberian rivers, which are all of typical reservoirs are important locations for the introduction of non- Mediterranean climate, were sampled for bleak along a latitu- 15 native fishes, such as the following trophy species: northern pike dinal gradient (from north to south): Ebro (4184703000N, Esox lucius, largemouth bass Micropterus salmoides, Eurasian 180502400W and 4182603900N, 082802600W), Tagus (3984900700N, 20 perch Perca fluviatilis, European catfish (also known as sheat- 482002500W and 3985801200N, 484201700W), Guadiana fish) Silurus glanis and pikeperch Sander lucioperca. Common (3885905800N, 585105200W and 3885001800N, 681305300W), Segura bleak was first introduced to Spanish reservoirs from southern (3880600500N, 181705000W and 3880405400N, 085303700W) and 20 France during the 1990s as a ‘forage fish’ for the above- Guadalquivir (3783700100N, 583504500Wand3783004600N, mentioned piscivorous fishes (Vinyoles et al. 2007), including 585603800W). These geographic coordinates (given in upstream to 25 the rivers Ebro (in 1992), Guadiana (in 1999), Segura (in 2004) downstream order) correspond to two sampling sites per river. (Andreu-Soler et al. 2004), Tagus (in 2005) and Guadalquivir (in The River Saoˆne reference site (4780201300N, 580605200E), where 2006). Since its introduction, the bleak has displayed a strong the climate is temperate oceanic, was then surveyed in early June. 25 invasive character throughout Iberian fresh waters by natural To provide an appropriate assessment of reproductive traits, and human-assisted (mainly stocking) dispersal, becoming fish were collected from May to June 2012, just before the 30 abundant within Iberian fish assemblages (e.g. Maceda-Veiga spawning period of bleak for each study area (authors, pers. AQ3 et al. 2010). Nevertheless, Iberian fishes are poorly adapted to obs.). Sampling sites (n ¼ 11, two sites per Iberian river plus one sustain the many introductions of non-native fishes, which on the River Saoˆne) were selected to encompass similar envi- 30 usually results in native species being displaced from available ronmental conditions. In particular, land use was mainly for resources. In particular, bleak threatens serval Iberian fishes, agricultural exploitation, sites were located in well-regulated 35 mainly through trophic and interference (i.e. aggression) com- middle reaches of the main channel, and far from the influences petition (Leunda 2010). As an example, direct observations of of main tributaries and towns. Also, sites were separated bleak by snorkelling have clearly shown that bleak can disturb .50 km within each Iberian river. A variety of consistent 35 the foraging behaviour of endemic Ebro nase Parachondrostoma sampling protocols was followed by wading and from boats, miegii (Almeida and Grossman 2012). according to the European legislation (CEN/ISO Standards, EC 40 From a conservation perspective, there is a clear need for Directive 2014/101/UE, European Commission 2014). This information on bleak population traits across the Iberian Penin- allowed obtaining a representative sample of bleak across the sula with which to assess the risks posed by this species to native broadest possible TL-size range from each river. Catch methods 40 species and ecosystems (Almeida et al. 2013). In particular, consisted of following a zig-zagging and upstream direction in growth and reproduction are key aspects for assessing the both banks at each site (100-m river length) by electrofishing 45 phenotypic variability of non-native fishes, which can facilitate (2000 W-pulse DC generator at 200–250 V, 2–3 A, 30 min per the invasion process (e.g.ROOF Copp and Fox 2007). However, bank), dip nets ONL (1.5-m-long pole, 30-cm-diameter net, 10-mm variations in bleak population traits have received relatively mesh size), seine nets (20 2 m, 10-mm mesh size) and gill-nets 45 little study across Europe in general and in the Mediterranean (20 1 m, 10-mm mesh size, 50% hanging ratio, 1.5 m deep). region in particular. The available literature on bleak autoecology Electrofishing was performed the first sampling day and nets 50 are simply relatedP to species’ spatial distribution (Vinyoles et al. were used the second day. All Iberian surveys followed the same 2007), habitat preferences (Mun˜oz-Mas et al. 2016) and a few sampling protocols (e.g. proportional effort in terms of people and population traits restricted to local scales (Almeida et al. 2014; time) to ensure comparability among the study rivers, with 50 Maso´ et al. 2016). methods being agreed and applied by each research team per The aim of the present study was to assess the inter-population river). A professional fisherman sampled for bleak in France. To 55 variability of traits invasive bleak across the Iberian Peninsula, encompass the existing environmental variability, fish were including back-calculated total lengths (TLs) at age, growth rate, collected from all meso-habitats present in the study rivers (e.g. body condition, reproductive investment, fecundity, egg size, age runs, pools, shallows). Fish were stored on ice during transport to 55 and TL at maturity and sex ratio. For this purpose, bleak the laboratory, before being frozen ( 208C). All field procedures populations were compared across the principal Iberian rivers complied with animal use and care regulations of Europe. Invasive bleak across the Iberian Peninsula Marine and Freshwater Research 3
Bleak were defrosted and measured for TL ( 1 mm), scale were expressed as a percentage of each eTLA; and (4) the mean samples taken from the dorsal area on both flanks, and then percentage from different ages was used to calculate a growth dissected and weighed for eviscerated mass (eW, 0.01 g). Sex index (GI, %) for each bleak individual. ANOVAs were used to was determined, as well as the maturity status in females, which test for significant differences of back-calculated TLs and GI 5 were classified as sexually mature if their ovaries contained yolked between populations. 5 oocytes (e.g. Tarkan et al. 2016). Reproductive traits shared by To provide an integrated quantification of body condition both sexes (i.e. gonad investment, age and TL at maturity) were (BC), eW was used to avoid bias from gonad mass and gut analysed in mature females only (Copp and Fox 2007; Maso´ et al. content. To assess reproductive investment (RI) in mature 2016). Therefore, mature females were also weighed for gonad females, gW was used (Maso´ et al. 2016). Female fecundity 10 mass (gW, 0.1 mg). Sub-samples (0.01 g) were taken from (FEC) was estimated according to the formula 10 anterior, middle and posterior portions of each mature ovarian lobe. Yolked oocytes were counted in every sub-sample and FEC ¼ gW D measured for egg diameter (ED, 0.01 mm) by using a PC-based Y image analysis (Aphelion, ADCIS, Saint-Contest, France). where FEC is the number of yolked oocytes in the entire ovary, 15 Age was determined from scales using a micro-projector gW is the gonad (ovary) mass and D is the density of yolked (magnification 48 ) by counting true annuli from acetate oocytes (number of oocytes g 1 of ovarian tissue). ANCOVAs impressions (Tarkan et al. 2016). The total scale radius and that were used to test for significant differences in BC (covariate: 15 of annual increments were measured from the focus to the TL), RI, FEC and ED (covariate: eW). posterior edge along the anterior–posterior axis. Age determina- Mean age-at-maturity (AaM) of each population was calcu- 20 tions were completed independently by two readers and when lated from the percentage of mature females in each age class, the interpretations were different, an additional reading was using the following formula (adapted by Fox and Crivelli 2001): taken. Age determination was validated by opercula examina- tion from sub-samples (10% of individuals per river) under a XAmax ¼ binocular microscope (40 ). AaM AMA MðA 1Þ 25 Non-linear and linear equations were fitted to determine which 20
model best described the relationship between TL and scale where A is the age in years, MA is the proportion (from 0 to 1) of radius. As the resulting models were linear equations for the six mature fish at Age A, and Amax is the maximum age in the fish bleak samples, back-calculation of TLs at age were estimated by sample. A modified version of this formula (10-mm-TL inter- the Fraser–Lee equation (adapted from Francis 1990), as follows: vals in place of age classes) was used to calculate mean TL at maturity (LaM), as per Fox and Crivelli (2001). 25 ¼ þð Þ ð : Þ TLA TL0 TLC TL0 RA RC General linear models (GLMs) were used to test for differ- 30 ences between males and females on the examined growth traits where TLA is TL when Growth mark (i.e. annulus) A was formed, (i.e. back-calculated TLs, GI and BC), as per Almeida et al. TLC is TL at the time of capture, RA is the distance (i.e. scale (2014). Because no difference was found, this categorical factor radius) from centre to the Annulus A, RC is the total scale radius at (i.e. sex) was not included in subsequent models, so as to 30 the time of capture and TL0 is the intercept on the length axis from increase the statistical power of the remaining sources of 35 linear regression between all TLC and RC values. Thus, the overall variation (see the same procedure in Alcaraz and Garcı´a- intercept TL0 acts as a ‘weighting factor’ to reduce bias resulting Berthou 2007). Data were pooled per Iberian river because the from differences in the size distribution of the examined popula- effect of ‘sampling site’ within each river in previous general- tions (see Maso´ et al. 2016 and Tarkan et al. 2016 for the same ised, linear, mixed models (GLMMs) was not significant (using 35 procedure). The particular TL0 value was 14 mm, which was used ‘site’ as the random factor) – a comprehensive review of this 40 as a fixed TL–scale intercept. This is in accordance with Maso´ statistical technique is given in Johnson et al. (2015). ANOVAs et al. (2016), where an intercept of 12 mm was found for four and ANCOVAs were followed by post hoc Tukey–Kramer Iberian bleak populations in north-eastern Spain. honestly significant difference (HSD) tests. Sex ratio (SR, For comparisons of growth rates, the Hickley and Dexter male : female) was compared among rivers using pairwise 40 (1979) procedure was followedROOF by using back-calculated data, Chi-Square (x2 ONL) tests (with Yates’ correction). Data were 45 as follows: (1) mean TLs from all populations were calculated at transformed into log (x þ 1) except percentage data, which þ 10 every age A, which were plotted against each mean TL at age A 1 were logit-transformed (Warton and Hui 2011). Assumptions of to obtain a regression line for the Walford (1946) method (used normality of distributions and homogeneity of variances were below); (2) expected TLs at every age (eTLA) were obtained from verified through Shapiro–Wilks and Levene’s tests respectively. 45 the formulae (adaptedP from the von Bertalanffy growth equation) Statistical analyses were performed with SPSS v19 (SYSTAT A Software Inc., Chicago, IL, USA). The significance level was set eTL ¼ TL1 ð1 S Þ A at P ¼ 0.05. Sequential Bonferroni corrections were performed 50 and for every set of multiple tests.
TL1 ¼ TL Cð1 SÞ Y Results and discussion 50
where TLY and S are the Y-intercept and the slope of the Walford In total, 1200 bleak (n ¼ 100 individuals 2 sites 5 Iberian plot respectively; (3) all back-calculated TLs at every age A rivers þ 200 individuals 1 site from the River Saoˆne) were 4 Marine and Freshwater Research D. Latorre et al.
a ab b b a b collected and processed, the minimum and maximum TLs being 44 and 199 mm. Significant differences were found between the six study populations (Table 1) for back-calculated TL at Age 1 (F5,1194 ¼ 32.89, P , 0.001), Age 2 (F5,896 ¼ 41.41, P , 0.001) and Age 3 (F5,517 ¼ 27.66, P , 0.001). Only the Tagus popu- 5 lation differed significantly in TL at Age 4 (F5,189 ¼ 12.47, P , 0.001). Significant differences were found among popu- 2.291.44 95.82.18 80.42.14 1.53 85.02.87 2.08 90.82.19 3.17 92.0 3.35 83.3 1.41 3.00 AQ4 lations for GI (F5,1194 ¼ 20.41, P , 0.001) and BC/eW (F5,1193 ¼ 7.12, P , 0.001). Body mass ranged from 0.65 to c c a ab bc ab 53.74 g. Among reproductive traits, ovary mass ranged from 10 0.02 0.02 0.04 0.02 0.03 0.03 0.22 to 6.69 g. Significant differences were found among
), growth index (GI, %), body condition ¼ , i populations for RI (F5,328 21.51, P 0.001). Absolute FEC A 1.14 1.18 0.95 1.05 1.11 1.02 (
i ranged from 732 to 12 380 eggs, with significant differencesY at maturity (AaM, years), length at maturity found among populations (F5,328 ¼ 16.89, P , 0.001). Oocyte a a b ab d c size ranged from 0.73 to 1.41 mm ED, again with significant 15 382 413 484 496 364 527 differences found among populations (F5,328 ¼ 14.34, P , 0.001) -value from Bonferroni correction) 1829 1864 3162 2800 8069 5085 Bleak from the River Tagus had the fastest growth rate, P highest body condition value, and the lowest gonad mass and c c a b b b egg size values. Whereas Ebro bleak has the highest values for 20 critical
0.09 0.10 0.15 0.12 0.09 0.13 mean gonad mass and mean egg diameter, but the lowest ,
P fecundity of all Iberian populations studied. Segura bleak had
2.48 2.54 1.02 1.44 1.68 1.69 the lowest growth rate and the highest fecundity. This variety of population strategies is probably displayed to compensate the
b ab c ab b a particular environmental pressures of novel habitats and facili- 25 0.70 0.70 0.77 0.68 0.67 0.69 tate bleak establishment (e.g. see a comprehensive review in Sakai et al. 2001). (only for SR) tests ( 8.31 2 The youngest and oldest mean AaMs were found in the rivers 12.59 11.19 17.58 10.38 12.80 x Ebro and Segura respectively, whereas the shortest and longest mean LaMs were observed in the river Ebro and Saoˆne 30 ab a c b a ab respectively (Table 1). All of the bleak populations were 0.8 0.5 0.9 0.5 0.7 0.6
GI BC RI FECdominated by ED males, AaM with LaM SR SR being lower in the rivers Segura and Saoˆne, whereas the rivers Guadalquivir, Tagus and Guadiana 93.3 92.4 100.9 120.1 106.3 100.0 has the second highest SR values, and the River Ebro an intermediate value (Table 1). 35 a a b a a a Bleak from the rivers Saoˆne and the closest Iberian river, the, 4.3 3.1 3.6 3.2 3.4 4.9 4 (LaM, mm TL) and sex ratio (SR, male : female) had similar population responses (i.e. growth and reproductive s.e., after ANOVAs or ANCOVAs (see details in Materials and methods section). Significant differences among river populations A