Oxygen Depletion Affects Kinematics and Shoaling Cohesion of Cyprinid Fish

Abstract

Numerous anthropogenic stressors impact rivers worldwide. Hypoxia, resulting from organic waste releases and eutrophication, occurs very commonly in Mediterranean rivers. Nonetheless, little is known about the effects of deoxygenation on the behavior of Mediterranean freshwater fish. To fill this knowledge gap, we assessed the impact of three different dissolved oxygen levels (normoxia, 48.4%, 16.5% saturation) on kinematics indicators (swimming velocity, acceleration, distance traveled) and shoaling cohesion of adult Iberian barbel, Luciobarbus bocagei, a widespread cyprinid species inhabiting a broad range of lotic and lentic habitats. We conducted flume experiments and video-tracked individual swimming movements of shoals of five fish. Our results reveal significant differences between the treatments regarding kinematics. Swimming velocity, acceleration, and total distance traveled decreased stepwise from the control to each of the two oxygen depletion treatments, whereby the difference between the control and both depletion levels was significant, respectively, but not between the depletion levels themselves. Shoaling cohesion showed dissimilarities between the treatments regarding the maximum distance between fish, as the high depletion treatment differed from each of the other two, indicating that under severe oxygen depletion some individuals move away from the shoal. Overall, our results show how oxygen depletion changes fish behavior, which may entail ecological responses, highlighting the need to maintain an unfragmented river network to ensure movement dispersal among habitats, thus providing conditions for species escapement from hypoxia.

1. Introduction

Mediterranean rivers are characterized by frequent hydrological disturbances, including seasonally predictable drought events, exposing biota to harsh environmental conditions [1]. An extensive period of low summer flows naturally leads to an increase in temperature and primary production, and subsequently to a reduction of oxygen levels [2]. Mediterranean fish populations are adapted to naturally occurring seasonal oxygen fluctuations [1]. However, human activities, such as water abstraction, the discharge of organic wastes, and the diffusion of agricultural fertilizers into aquatic systems, are common problems in Mediterranean rivers, causing further diminishment of oxygen levels of already stressed freshwater systems and possibly leading to hypoxic conditions [3,4] that affect fish and subsequently the entire aquatic community [3,4,5,6].

Hypoxic conditions can cause behavioral and physiological stress in fish [7], which can entail many ecological consequences, including impaired reproduction, feeding, and predator avoidance [8,9]. Regarding swimming activity, fish experience a trade-off between activity increase as they try to find more oxygenated water and activity decrease due to lower oxygen availability [10]. When exposed to hypoxia, fish groups may reduce their school density and size as they tend to increase their horizontal spacing to one another in a tentative manner to counteract the limiting effect of oxygen depletion on schools [5]. However, this may cause disadvantages, such as lower coordination or higher energetic costs of locomotion [11]. Both activity changes, individual fish movements and shoaling cohesion, are, however, species-dependent and are influenced by the degree of hypoxia [5,9]. Too low oxygen saturation will at some point be lethal for fish, whereby the lethal concentration also differs between species [12].

The Mediterranean region has been identified as one of the most prominent hot spots of climate change [13]. According to future predictions of climate change for Southern Europe, precipitation will decrease by up to 20% [14], which will also entail significant changes in seasonal runoff asymmetry [15]. Hence, exacerbated and extended summer low flow conditions, in combination with existing human pressures, will likely increase the frequency and duration of hypoxic events in rivers [1,15,16]. To protect and manage rivers, it is vital to understand the impact of deoxygenation on individual fish movements, schooling behavior, and the implications thereof for fish populations, which are a pivotal component of aquatic ecosystems.

To our knowledge, no study has assessed the effects of oxygen depletion on the movements of Mediterranean freshwater fish. Therefore, this work analyzed the impact of three different dissolved oxygen levels on the behavior of a common fish species of the Iberian Peninsula, the Iberian barbel, Luciobarbus bocagei. We hypothesized that increasing levels of oxygen depletion would (1) reduce fish movement and (2) decrease the cohesion of fish shoals.

2. Materials and Methods

We conducted flume experiments in a section (2.0 × 1.0 m) of an artificial channel, set at a slope of 3% and a discharge of 28 L∙s⁻¹, resulting in a water depth of 0.7 m. We caught 75 wild adult Iberian barbel, Luciobarbus bocagei (150–250 mm total length), via electrofishing in the Lizandro River, Portugal (33 m a.s.l.), and brought them to the lab facilities (92 m a.s.l.), where they were held in three 800 L tanks (25 fish per tank) in a controlled environment for 48 h before experimentation. After the trials, fish were returned alive to the river. Further details on fish sampling, holding, and the testing facility are given in Branco et al. [17,18] and Santos et al. [19].

In each experiment, we monitored the movement of a shoal of five fish [20] under three levels of dissolved oxygen (DO) conditions: (1) a control level (normoxic condition with DO saturations >80% [21]), (2) a mild depletion level of 48.4% DO saturation (4.4 ± 0.2 mg∙L⁻¹), and (3) a high depletion level of 16.5% DO saturation (1.4 ± 0.1 mg∙L⁻¹) (

Table 1. Physiochemical and biological description of the experiments (average values ± SD), and number of analyzed videos and frames per treatment.

	Variable	Treatments
		Control	Mild Depletion	High Depletion
Water quality	O2 (%)	(Normoxia)	48.4 ± 1.9	16.5 ± 1.2
	O2 (mg∙L−1)	(Normoxia)	4.4 ± 0.2	1.4 ± 0.1
	Temperature (°C)	19.4 ± 0.4	19.6 ± 0.3	22.7 ± 0.1
	Conductivity (μS∙cm−1)	152 ± 4.5	492 ± 16.4	738 ± 33.4
	pH	8.1 ± 0.1	8.2 ± 0.1	8.4 ± 0.1
Fish	Total length (mm)	202 ± 23.9	213 ± 29.2	191 ± 33.0
	Body mass (g)	81 ± 27.5	94 ± 37.5	70 ± 37.7
Video	Analyzed videos	5	5	4 1
	Analyzed frames	30,025	30,025	24,020 1

¹ Due to a technical video failure, one replicate of the high depletion treatment level had to be excluded from further analysis.

). These levels allowed for testing along a decreasing gradient, whereby the latter represents a level of extreme oxygen depletion that may induce stress or even death in fish [6,22]. We reduced oxygen values by adding sodium sulfite (Na₂SO₃). This compound is a recognized oxygen scavenger [23] that has been used to create oxygen-deficient conditions in fish research [18,24], as well as in aquatic research more generally [25,26]. Each experiment lasted for eight minutes (six min acclimation [27] plus two min video-tracking) and was replicated five times, whereby a new subset of fish was used for each trial, giving a total of 75 tested fish.

The experiments were conducted in strict accordance with ethical provisions on the welfare of experimental animals enforced by the European Union and were coordinated by J.M. Santos, who holds FELASA level C certification (www.felasa.eu) to direct animal experiments. Fish sampling permits were obtained from the Institute for Nature Conservation and Forests (ICNF). Fish experiments and maintenance in the laboratory and experimental facility were authorized by the Department for Health and Animal Protection (Direcção de Serviços de Saúde e Protecção Animal; permit issued on 6 October 2011) in accordance with the recommendations of the “protection of animal use for experimental and scientific work”.

We recorded fish movements with a GoPro Hero 3 Black edition camera mounted above the channel. The camera was set to record at 1280 × 720 pixels. We processed the videos with GoPro Studio to exclude the fisheye effect and then analyzed them at 10 frames per second, resulting in 1200 data points per single fish (Table 1) with the software Tracker (http://physlets.org/tracker/), whereby we calibrated the measurement settings of each video with the size of the channel’s bottom tiles (20 × 20 cm).

We calculated kinematics indicators, i.e., mean swimming velocity (cm∙s⁻¹), acceleration (cm∙s⁻²), and distance traveled (cm) for every single fish to investigate spontaneous locomotion activity changes between DO treatments [27]. To test for differences in shoaling cohesion between treatments, we measured the horizontal (XY) inter-individual distance between all fish within a shoal, based on center of mass [27,28], whereby we used the average of the minimum, mean, and maximum values as a proxy for the shoal cohesion. In all analyses, we first tested for overall statistical significance between the treatment levels. If this was detected, it was followed by pair-wise post-hoc tests. For most of the tests, data did not fulfill all assumptions required for parametric tests. Therefore, omnibus tests were conducted with the non-parametric Kruskal-Wallis H test, followed by Dunn’s pairwise tests (adjusted using the Bonferroni correction) (α = 0.05) using SPSS 24. The dataset to this article can be found online in the Supplementary Material Section (Datasheet S1).

3. Results

3.1. Kinematics

We detected significantly different swimming velocities between the treatments (χ² = 21.37, p < 0.001, Figure 1a). Pairwise comparison found that fish in the control treatment (T0) swam faster than those in the mild depletion treatment (T1) (χ² = 18.92, p < 0.01), as well as those in the high depletion treatment (T2) (χ² = 26.96, p < 0.001), whereas the mild and high depletion treatment (T1 and T2) were not distinct from each other. The average swimming velocity was 10.8 (±0.3 Standard Error (SE)), 8.6 (±0.5 SE), and 7.9 (± 0.4 SE) cm∙s⁻¹ for T0, T1, and T2, respectively.

For acceleration, we found a similar pattern (χ² = 34.26, p < 0.001, Figure 1b), whereby the pairwise comparison was highly significant for both T0 and T1 (χ² = 27.40, p < 0.001), as well as T0 and T2 (χ² = 32.04, p < 0.001). Average fish acceleration was 27.1 (± 0.6 SE), 20.7 (± 1.0 SE), and 20.3 (± 0.7 SE) cm∙s⁻² for T0, T1, and T2, respectively.

We found significant differences in the total distance traveled by each fish among the treatments (χ² = 13.10, p < 0.001, Figure 1c), though pairwise comparison revealed that fish only traveled longer distances in T0 (1247 cm ± 44.0 cm SE) when compared to T2 (987 cm ± 45.6 cm SE) (χ² = 21.51, p < 0.001). T0 and T1 (1055 cm ± 57.2 cm SE) were marginally significant (χ² = 13.76, p = 0.050).

3.2. Shoaling Cohesion

Average minimum distances between fish ranged from 8.5–14.4 cm, mean distances from 46.6–67.8 cm, and maximum distances from 82.2–124.5 cm (

Table 2. Average minimum, mean, and maximum distances (cm ± SE) between fish within a shoal during the three treatments.

	Distance (cm ± SE)
Treatments	Minimum	Mean	Maximum
T0—Control (normoxia)	14.4 ± 2.7 a	48.3 ± 8.5 a	82.2 ± 16.0 a
T1—Mild depletion (48.4% DO)	8.5 ± 0.9 a	46.6 ± 6.7 a	88.4 ± 13.0 a
T2—High depletion (16.5% DO)	10.5 ± 0.4 a	67.8 ± 15.4 a	124.5 ± 26.0 b
Kruskal-Wallis test	χ2 = 1.16, p = 0.56	χ2 = 4.77, p = 0.09	χ2 = 8.89, p < 0.01

“a” and “b” indicate a significant difference between two treatments, whereby the same letter indicates no statistically significant difference by Kruskal-Wallis test.

). Regarding the minimum and mean distance between the fish, we did not detect differences between the three treatments. However, for the maximum distance between fish, we found significant dissimilarities between the treatments (Table 2), as the distance between fish was higher in T2 when compared to T0 (χ² = −15.02, p < 0.05) as well as T1 (χ² = −16.90, p < 0.05).

4. Discussion

Hypoxic events affect Mediterranean rivers worldwide [2,29], and fish, in particular, are highly sensitive to reduced oxygen concentrations [9,30]. The ecological impact of such events ranges from beneficial to mortality [31], but common responses include behavioral alterations [5,9,18].

4.1. Kinematics

Our results showed that swimming velocity, acceleration, and distance traveled significantly differed between the treatment levels. Pair-wise tests revealed that fish in the control treatment exhibited higher spontaneous activity than those in the mild depletion as well as high depletion treatment, which implies that barbel reduce their activity if DO values are lowered. A similar behavioral response is documented for several other species [9,32]. Such lowered activity can be interpreted as an adaptive energy-saving response [33]. Carassius carassius, e.g., reduces its locomotor activity to 50% of that during normoxia when exposed to prolonged periods of anoxia, thereby saving 35–40% of total energy consumption [34]. Swimming speed of Ammodytes tobianus decreased by 95% after 40 min at the lowest stepwise reduced critical partial pressure of oxygen (3.1 kPa) [32]. Even though the acclimation phase and analysis time in our study were short, fish responded quickly to decreased DO levels by significantly reducing their movement activity.

4.2. Shoaling Cohesion

Under low DO (20%) conditions, school measurements of Clupea harengus significantly increased [35]. We found a similar pattern regarding the maximum distance between individuals at the high depletion treatment in comparison to both the control and the mild one, whereas this was not the case for the minimum and average distance. Nevertheless, this might indicate that at approximately 15% DO, some barbel start to depart from the shoal, while others still seem to prefer the close vicinity of the shoal. In such high oxygen depleted situations, individual movement appears to override shoal movements as individuals risk leaving shoal protection to find more oxygenated areas further away, thus reducing shoal integrity by splitting the shoal into two or more groups [11].

4.3. Ecological Implications

In natural situations, fish might try to escape unfavorable conditions. In a hypoxic environment, Oncorhynchus mykiss selected lower temperature ranges, enhancing their chance of survival [36]. In the confined flume section of the present study, barbel seemed to reduce their swimming activity, possibly adopting a “sit-and-wait” strategy, which has also been described for other species [5]. In connected pool habitat experiments, however, barbel continued to move despite reduced oxygen levels, possibly searching for higher-oxygenated areas. Nevertheless, the number of movements significantly decreased under low (15%) DO concentrations [18]. Under such oxygen-depleted conditions, isolated habitats may, therefore, cause severe threats to stream fish, highlighting the need to maintain an unfragmented river network to ensure movement dispersal among habitats, and thus providing conditions for species escapement from hypoxia.

In our experiments, the addition of sodium sulfite not only reduced DO levels but also led to an increase in water temperature, conductivity, and pH, making it difficult to completely disentangle the effects of lowered DO from the increase of the other parameters. Nevertheless, similar changes in these variables can often occur with the release of organic stressors, such as sewage, into rivers [18,25]. To increase understanding of fish resistance to oxygen-depleted conditions and its implications for aquatic ecosystems, future studies should investigate the effects of long-term DO depletion on riverine fish species.