Evaluating dispersal potential of an invasive fish by the use of aerobic scope and osmoregulation capacity

Non-indigenous species (NIS) can impact marine biodiversity and ecosystem structure and function. Once introduced into a new region, secondary dispersal is limited by the physiology of the organism in relation to the ambient environment and by complex interactions between a suite of ecological factors such as presence of predators, competitors, and parasites. Early prediction of dispersal potential and future ‘area of impact’ is challenging, but also a great asset in taking appropriate management actions. Aerobic scope (AS) in fish has been linked to various fitness-related parameters, and may be valuable in determining dispersal potential of aquatic invasive species in novel environments. Round goby, Neogobius melanostomus, one of the most wide-ranging invasive fish species in Europe and North America, currently thrives in brackish and fresh water, but its ability to survive in high salinity waters is unknown to date. We show that AS in round goby is reduced by 30% and blood plasma osmolality increased (indicating reduced capacity for osmoregulation) at salinities approaching oceanic conditions, following slow ramping (5 PSU per week) and subsequent long-term acclimation to salinities ranging between 0 and 30 PSU (8 days at final treatment salinities before blood plasma osmolality measurements, 12–20 additional days before respirometry). Survival was also reduced at the highest salinities yet a significant proportion (61%) of the fish survived at 30 PSU. Reduced physiological performance at the highest salinities may affect growth and competitive ability under oceanic conditions, but to what extent reduced AS and osmoregulatory capacity will slow the current 30 km year-1 rate of advance of the species through the steep salinity gradient from the brackish Baltic Sea and into the oceanic North Sea remains speculative. An unintended natural experiment is in progress to test whether the rate of advance slows down. At the current rate of advance the species will reach the oceanic North Sea by 2018/2019, therefore time for taking preventative action is short.


Introduction
factor underpinning the success of this particular NIS. Indeed, round goby is currently expanding into areas of still higher salinities; from the most western part of the Baltic Sea secondary dispersal is taking place northward along the coastline at a rate of approximate 30 km year -1 [48]. This area comprises a transition zone between the brackish Baltic Sea and the full oceanic North Sea, where at present the range of the round goby covers salinities from approximately 10 to 17 PSU (Fig 1). An unintended natural experiment is thus in progress to test whether the rate of advance slows down.
By investigating the effect of salinity (0-30 PSU) on two physiological traits, namely AS and osmoregulation, the present study provides a mechanistic assessment of the secondary dispersal potential of the species within the transition zone between the brackish Baltic Sea and the oceanic North Sea, which is characterised by a steep salinity gradient. We expected that the highest AS for round goby would occur at near isoosmotic salinities (10PSU, the salinity at the location of origin of the experimental fish), with decreasing AS at salinities deviating from these values. We also expected that round goby would have limited ability to osmoregulate, that is, to maintain stable blood plasma osmolality, at higher salinities. ] were on average 0.05 and 0.1 mg L -1 , respectively, from at least two week before initiation of experiments and throughout the experimental period, and pH had stabilized at 7.9.
The equations were as follows: There was a significant effect of salinity on blood plasma osmolality (Kruskal-Wallis test, p = 0.001). The blood plasma osmolality was not statistically different at salinities between 0 and 25 PSU, and the average within these treatments ranged from 332 and 352 mOsm kg -1 . However at 30 PSU, the blood plasma osmolality had increased to 372 mOsm kg -1 , significantly higher than in the 0-20 PSU treatments (Mann-Whitney U tests, p < 9.743 Ã 10 −4 ) (Fig 3).
There was a significantly negative relationship between osmotic potential (blood plasma osmolality minus ambient osmolality) and AS (Spearman Rank Correlation, r s = -1, n = 6, p < 0.01).

Discussion
Maintaining physiological performance in novel environments is the most fundamental constraint potentially limiting dispersal of invasive species. The results support our prior expectation and show that the AS of round goby was significantly reduced at the highest salinities (25 and 30 PSU). The same pattern was seen for the osmoregulatory capacity, where round goby were able to maintain stable blood plasma osmolality up to 25 PSU, but had increased blood plasma osmolality at 30 PSU, indicating that 25 PSU is the species' limit for unperturbed osmoregulation. Nevertheless, there was fairly high survival (i.e. 61%) of round goby at 30 PSU (Fig 4), and some individuals were able to maintain low blood plasma osmolality at the high salinities.
The present reduced AS is the outcome of the combined effects of increasing SMR at salinities departing from 10 PSU, and MMR following the inverse pattern, being reduced at 25 and 30 PSU. The augmented SMR likely reflects the additional costs of maintaining internal homeostasis, governing not only ion transport per se, but all factors related to this, including hormone level change and control and also potential underlying stress responses at the more extreme salinities [31, [49][50][51]. The gills are the principal organs of osmoregulation, but because of their involvement in respiration, they are at the same time the main site of water and ion leakage [30]. The reduced MMR observed at the two highest salinities is thus likely caused by reduced diffusion potential of the gill epithelium, initiated to reduce the osmotic water loss, which then simultaneously reduces the ability to take up oxygen (the osmo-respiratory compromise) [29,30]. In combination, this agrees well with the present significant negative correlation between AS and osmotic potential, suggesting the latter to be one (of likely several) determining factor on AS.
Using AS to predict future species distribution during climate warming has gained much attention due to the oxygen and capacity limiting thermal tolerance (OCLTT) hypothesis [25,52], which is an extension of the Fry Paradigm [20,27]. The OCLTT hypothesis claims that AS is highest at an individual's thermal optimum and drops as individuals approach their thermal limit, the latter caused by a failure of the fishes capacity to supply oxygen to tissues at temperature extremes [25,52]. In spite of recent questioning of the general applicability of the OCLTT hypothesis [53][54][55], the results of the present study adhere to the Fry Paradigm and demonstrate impairment of oxygen transport systems when salinities exceed those which are optimal for AS. Thus, AS might be useful when evaluating distribution patterns of invasive fish species though salinity gradients. With the present substantial survival (61%-and the survivors obviously tolerating an increased plasma osmolality) and a mere 30% reduction of AS at the highest salinity (30 PSU), the question is whether the repressed physiological performance will hinder successful northward invasion into the full oceanic North Sea. Although growth potential and ability to compete with native species for resources could be affected in this high-saline environment [50,56], it may be that 70% of the round goby's current AS is sufficient to allow a population to establish there.
Previous studies on adult round goby with fresh water origin (Laurentian Great Lakes) have shown poor or no survival at salinities of 20 PSU or higher. In the study by Karsiotis et al. [57] all individuals died within 1 month at salinities of 25 and 30 PSU, and only 5% survived at 20 PSU, while Ellis and MacIsaac [58] showed a mortality of 100% after 48h at 30 PSU (only salinity tested). This is in stark contrast to the present study showing extensive survival of fish after three months at 30 PSU. The round goby in our study originated from Baltic brackish water, and thus the origin of the fish (fresh versus brackish waters) seemingly affects the species' ability to cope with salinity stress. Furthermore, Ellis and MacIsaac [58] and Karsiotis et al. [57] used acute exposure or very rapid salinity acclimation (up to 10 PSU h -1 ) in contrast to our study, where the salinity increase was 5 PSU week -1 . In Adriatic sturgeon acclimation to brackish water prior to an increasing salinity challenge test enabled the fish to maintain stable blood plasma osmolality [51], apparently by stimulating a beneficial increase in the activity of ion-translocating enzymes in the gills [59].
Gradual activation of various physiological functions (cardiac, metabolic and locomotor) in extreme environments has been shown in both aquatic and terrestrial animals [14,[60][61][62][63][64][65], and, together with the present results, point to the need to take acclimation into account when trying to understand and predict NIS performance in novel environments.
In addition to acclimation, several mechanisms exist by which physiological constraints may be overcome, hence enabling dispersal of individuals into and through unfavorable environments. It is now widely recognized that epigenetically induced changes in gene expression (i.e. changes that do not involve modification of DNA sequences) can evolve from a number of environmental stressors which then again determine how susceptible or resilient the individual is to specific stressors later in life ( [66] and references herein). Likewise, maternal effects may act as adaptive mechanisms by buffering effects of a changing climate and thereby helping to alleviate some of the related fitness consequences [67,68]. Furthermore, variations in the early rearing environment of genetically similar individuals can result in distinctively different physiological and behavioural responses when subject to stressors, thus enabling some individuals to cope better than others [69,70]. Notably, in our study, some individuals were able to maintain unperturbed blood plasma osmolality levels at 25 and 30 PSU, resulting in high betweenindividual variation within these groups. Although our study was not designed to test for consistency in individual responses, these between-individuals differences may reflect variation in plasticity of the salinity tolerance, suggesting the existence of individuals with superior performance in high saline environments. Recent studies on fish have demonstrated how interindividual differences in AS are associated with fitness-related components like dominance, position of individuals within schools, feeding capacity and growth [71][72][73][74]. Future studies targeting potential differences between individuals in salinity tolerance are needed to shed light on why some round goby appear to be more robust to environmental changes as these may be successful in relation to invasion of the North Sea.

Mechanistic Species Distribution Models (SDMs)
The potential value of statistical species distribution models (SDMs) to describe species distributions has been demonstrated (e.g. [75]), but the predictive capability of these models has also be disputed, as they merely describe patterns in past time [76]. Mechanistic SDMs, parametised with physiological traits, are more robust when extrapolating distributions to future environmental conditions [76][77][78][79]. The critical issue is, however, which traits to incorporate into the SDM. Behavioral traits related to habitat preference [75], aggressiveness [80], or generalist vs. specialist strategies [81] may be relevant. However, the importance of behavioral traits may vary even on a short time scale (i.e. year-to-year), if the magnitude of competition or predation pressure fluctuates. Suboptimal habitats may not prevent colonization. For example, the round goby prefers relatively shallow protected areas of rocky substrate, but has been observed at water depth >20 m (own observations) and also appears to be thriving on mud bottom with only few rocks [82] and in tributaries [75]. The scope and applicability of the OCLTT hypothesis is under debate, yet many studies suggest that AS provides a more robust predictor of the ultimate geographical distribution [25,83], as it links environmental properties directly to the physiological performance of the individual, including growth [73,74,84] and swimming performance [71,85,86] and indirectly via biotic interactions [87]. For example, Marras et al. [87] demonstrated that the relationship between temperature and AS provides a trait-based tool to evaluate the potential invasiveness of a NIS across temperature gradients.
Physiological constraints can only be overcome through evolutionary adaptation, epigenetics, maternal effects and/or effects of the early rearing environment (see above), and this could potentially be incorporated into mechanistic SDMs as researchers are beginning to develop an understanding of the processes underpinning local adaptations in non-static, non-equilibrium scenarios such as invasions across environmental gradients [79,88,89].

Concluding remarks
In summary, our results suggest that physiological traits should be considered when evaluating the potential for dispersal of invasive species in novel environments. We found that salinity above 25 PSU causes a decrement in AS as well as a concomitant perturbation of plasma osmolality. This reduced physiological performance at the highest salinities may affect growth and competitive ability under oceanic conditions [50,56], but to what extent reduced AS and osmoregulatory capacity will slow the current 30 km year -1 rate of advance of the species [48] through the steep salinity gradient from the brackish Baltic Sea and into the oceanic North Sea remains speculative. An unintended natural experiment is in progress to test whether the rate of advance slows down. At the current rate of advance the species will reach the oceanic North Sea by 2018/2019, therefore time for taking preventative action is short.

Materials and methods
Ethical permit 2015-15-0201-00546 from the Danish Animal Ethics Committee covered all experiments reported here. Fish were observed twice daily by animal caretakers, the authors or the veterinarian at the aquarium. Humane endpoints were used, i.e. fish that showed signs of suffering or distress (ceasing of feeding, loss of equilibrium, increased ventilation frequency, and lack of normal movement and social interactions) were euthanized by a sharp blow to the head where after the head was cut off. Three weeks into the acclimatization protocol (i.e. when one group of fish was at 0 PSU, one group at 10 PSU and the remaining four groups at 15 PSU) necropsy of two euthanized fish revealed heavy infection (skin and gills) with Trichodina sp. Individual fish in all tanks were subsequently given a 20min freshwater bath [90,91] after which they were returned to their acclimation tank. Regular samplings of skin scrapes throughout the remaining experimental period did not reveal any further Trichodina sp. infections, and no signs of gill lesions/gill abnormalities were seen. From this, and also considering the short period of infection, it was assumed that Trichodina sp. had not affected the capacity for either oxygen uptake or ion exchange at the time when data collection commenced. For blood sampling for measurements of blood plasma osmolality, the fish were anaesthetized (see below for details).

Fish and experimental setup
Round gobies (approximately 200 fish) were caught with fyke nets in the brackish water estuary Guldborgsund (salinity 10 PSU [92]) in the western Baltic Sea (54˚42'N, 11˚51'E). They were transported to Den Blå Planet National Aquarium Denmark, Kastrup, Denmark, and held in freshwater (0 PSU) for 3 months until the experimental acclimation began. Fish were fed to satiation three times a week with 1.5-4 mm commercial high-nutrition fish feed pellets on which they fed well after 1-2 weeks of customization. Fresh-instead of brackish water was chosen due to practical reasons, as the species thrives well in fresh water (own previous experience; [93]). Altogether 100 adult fish of comparable size were allocated into five separate acclimation tanks containing 450 L 10 PSU water, with 20 fish in each tank. After one week at 10 PSU, the salinity was increased gradually by 5 PSU per week to obtain treatment salinities of 10, 15, 20, 25 and 30. Thus, the 10 PSU group were already at their final treatment salinity, the 15 PSU had reached their final treatment salinity after one week, the 20 PSU group after two weeks, the 25 PSU group after three weeks, while fish in the 30 PSU group had reached their final salinity after four weeks. A remaining batch of 20 fish, of comparable size to the abovementioned 100 fish, was kept in the fresh water, i.e. constituting the 0 PSU group. Once each group of fish reached their experimental salinity, they were kept at this salinity for 8 days and then blood was sampled for determination of blood plasma osmolality. Subsequently, the fish were left for a between 12 and 20 days before respirometry (for further details of experiments see below). Survival was determined by summing the number of dead fish in each treatment group over the entire period, i.e. from the start, when fish were separated into the tanks and acclimated to the different salinities until the end of all experiments. Plastic pipes and artificial seaweed were used as shelter structures for the fish in all tanks. The feeding regime continued as described above, and there were no signs of feeding ceasing at any point of time in any group of fish. Brackish water was obtained by mixing filtered oceanic sea water (30 PSU) and non-chlorinated tap water (0 PSU) and salinity measured with a digital refractometer (Pocket refractometer PAL-06S from Atago). Water temperature was 15-17˚C, and the daylight period held at 9 h. To maintain adequate water quality, the water in each tank was recirculated through 20 L canister filters with a capacity of 3500 L h -1 and 50% exchanged twice a week. Water quality was measured once a week, ammonium (NH 4 + ) and nitrite (NO 2 -) levels with a spectrophotometer (7500 Photometer, Palintest instruments Ltd, Gateshead, UK) and pH with a hand held pH meter (HQ30d Portable pH Meter, Hach, Loveland, US).

Respirometry
Forty-eight fish (n = 8 in each group) were used for respirometry (Table 1). They were kept for a minimum of 20 days at their final treatment salinity before respirometry experiments were initiated, a time period considered sufficient for their blood plasma osmolality level to be restored after the salinity change [49,94,95]. Individuals were fasted for two days before being placed in the respirometer. Four respirometers with volumes between 1126 and 1157 mL were held in an 80 L aerated, temperature-regulated water bath (18˚C, the range throughout the study was± 0.5˚C, and within one trial ± 0.25˚C), placed at an undisturbed site. The water was exchanged once a week. To minimize visual disturbance, the respirometers were shielded off from each other and from above with non-transparent polyethylene plates. Oxygen consumption rate (Ṁ O 2 ) was measured with an intermittent-flow respirometry set-up [96,97]. The partial pressure of oxygen (pO 2 ) was measured every second with fiber-optic sensors, connected to a Witrox 4 oxymeter (Loligo Systems, Tjele, Denmark). Using the software AutoResp™ (Loligo Systems, Tjele, Denmark) dpO2 dt (kPaO 2 h -1 ) was determined by linear regression and Ṁ O 2 (mg O 2 kg -1 h -1 ) calculated according to the equation: where α is the dpO2 dt , V resp is the total volume of the respirometer minus the volume of the fish (L), β is oxygen solubility at the given temperature and salinity and BM is the body mass of the fish (kg). The animal volume was assumed equal to its mass. Table 1. Size (mean ± SD) of round goby used in respirometry and for blood plasma osmolality experiments. Length is total length, and Fulton's condition factor (K) was calculated as K = 100×W×L -3 (where W is in g and L is total length in cm). Whole animal oxygen consumption rate (Ṁ O 2,whole animal ; mg O 2 h -1 ) scales allometrically with BM over very large BM ranges [98], yet with narrow BM ranges one can assume isometric scaling. To investigate whether BM in the present study was taken appropriately into account in the calculation of Ṁ O 2 (Eq 1), simple linear regressions analyses (Ṁ O 2,whole animal = a Ã BM + b) were conducted, one for MMR yielding an R 2 of 0.567 (F(1,46) = 62.554, p = 4.028 Ã 10 −10 ) and one for SMR yielding and R 2 of 0.698 (F(1,46) = 106.212, p = 1.547 Ã 10 −13 ). This supported that Ṁ O 2,whole animal here could be assumed to scale isometrically with BM, justifying using Eq 4.
Before and after each trial, background Ṁ O 2 was measured in the empty respirometers using a wait period of 400 s and a measurement period of 3600 s. To yield maximum metabolic rate (MMR), individual fish were manually exercised to exhaustion (between 3-5 min) in a circular tank (40 cm diameter, water depth 10 cm) [99], after which the fish was placed in the respirometer. The time until the experiment began was approximately 60 s, and the first Ṁ O 2 measurement taken with wait period of 60 s and a measuring period of 150 s. The background Ṁ O 2 from before the trial was subtracted from this initial Ṁ O 2 value, which was then considered MMR [99,100].
After measurement of MMR, fish remained in the respirometers overnight to allow for Ṁ O 2 measurement of the individual at a resting, non-digesting state (standard metabolic rate; SMR [81]) (flush period of 240 s, wait period of 60 s, measuring period of 300 s). A single experimental trial lasted approximately 22 h. The applied measurement periods yielded almost exclusively regression with R 2 > 0.99, and the oxygen content never fell below 70%. Ṁ O 2 values with R 2 < 0.95 were removed before data processing, and a double Gaussian distribution was fitted to the frequency distribution of the Ṁ O 2 values from the individual fish (Fig 5).
This algorithm separates the Ṁ O 2 associated with spontaneous activity and handling stress from those of the resting animal, and SMR was determined as the lowest mean value of the two Gaussian distributions [101].
Background Ṁ O 2 from before the trial was subtracted from MMR measurements, whereas the value after the trial was subtracted from the SMR measurements. Aerobic scope (AS) was expressed as MMR minus SMR. The background Ṁ O 2 was 1% of MMR measurements and 6% of SMR measurements (calculated on global averages).

Blood plasma osmolality
The fish (n = 9-10 in each group, Table 1) were kept for eight days at their final treatment salinity before blood sampling, a time period considered sufficient for their blood plasma osmolality level to be restored after the salinity change [49,94,95]. Following blood sampling fish were returned to their tank, and a randomised sample of these fish was used a minimum of 12 days later for respirometry. The 12 day period was deemed long enough to restore blood hematocrit, hemoglobin and osmolality [102,103]. The fish were anaesthetized in a 0.133 g L -1 3-aminobenzoicacid ethyl ester (MS-222) solution buffered with sodium bicarbonate (0.30 g l -1 ) [104] and blood subsequently sampled by caudal puncture with lithium-heparinized 27 gauge hypodermic needle [105]. The samples (approximately 200 μL) were immediately spun at 6000 rpm for 5 min, and plasma isolated and stored in a -18˚C fridge. Osmolality was subsequently measured in a vapour pressure point osmometer (Vapor Pressure Osmometer 5600, Wescor, South Logan, USA).

Statistical analyses
Statistics were performed in the software SPSS statistics 23 (IBM, Armonk, USA) and R [106] using a significance level (α) of 0.05. There was variance homogeneity for MMR, SMR, and AS (Levene's test, all p > 0.147), and normality for the data in all treatments (Shapiro-Wilk's tests, p > 0.102) except for SMR at a salinity of 25 (Shapiro-Wilk's test, p = 0.048). This single occurrence of non-normality was ignored since ANOVAs are relatively robust towards non-normality [107]. A general linear model was used to conduct one-way ANOVAs to test the effect of salinity on SMR, MMR and AS. The models also tested if the 0, 15, 20, 25, and 30 PSU treatments produced SMR, MMR, and AS significantly different from 10 PSU (the reference level), which is the condition nearest isoosmotic levels, as well as the environmental salinity in Guldborgsund where the round goby were caught. It is well known that some individuals are extremely difficult to exercise to exhaustion (i.e. to achieve true measurements of MMR) using standard tests [21]. Hence, to avoid confounding of the analysis by abnormal behaviour, we excluded data if the MMR value was found to be outside 2 standard deviations of the mean calculated for the remaining 7 data points in the treatment group (i.e. outside the 95% confidence interval). This resulted in the removal of two outliers (one from 10 PSU and one from 15 PSU). Since MMR constitutes a part of AS, the two individuals identified as outliers on MMR, were also disqualified for AS, but their SMR were still included in the analysis. The model also parameterized the response in relation to a predefined set of parameters, given as: where a,b and c are constants to be estimated and SAL is salinity. The second order term was included, as local maxima and minima were expected. The osmolality data was log transformed to obtain variance homogeneity (Levene's test, p = 0.114), but the transformation did not yield normality at a salinity of 25 (Shapiro-Wilk's tests, p = 0.025) and 30 (Shapiro-Wilk's tests, p = 0.002). Therefore, a Kruskal-Wallis' test was applied, followed by pairwise Mann-Whitney U tests with a Bonferroni correction for multiple comparisons (Bonferroni's α = 0.003).
A Spearman Rank Correlation was carried out between the average osmotic potential for each treatment ([blood plasma osmolality minus ambient osmolality]) and average for each treatment AS.
Supporting information S1