Open Access
Issue
Knowl. Manag. Aquat. Ecosyst.
Number 420, 2019
Article Number 50
Number of page(s) 12
DOI https://doi.org/10.1051/kmae/2019040
Published online 04 December 2019

© M.L. Božanić et al., Published by EDP Sciences 2019

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License CC-BY-ND (https://creativecommons.org/licenses/by-nd/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. If you remix, transform, or build upon the material, you may not distribute the modified material.

1 Introduction

Strong anthropogenic influence, discharge of industrial and communal wastewaters and runoff of artificial fertilizers from agricultural land are factors that cause accumulation of various substances, including metals, in rivers (Schueler, 2000; Çevik et al., 2009; Triquet-Amiard et al., 2013), which leads to deterioration of the quality of aquatic environments and threatens the survival of aquatic organisms. Together with other animal organisms, insects are constantly exposed to action of the indicated anthropogenic stressors, which cause daily and seasonal fluctuations of water temperature, changes of dissolved oxygen concentrations and food availability and the presence of pollutants (De Almeida et al., 2007; Despotović, 2013). Metals are among the most significant pollutants, in view of their stability in the ecosystem and transport through the food chain (Shahbaz et al., 2013) and considering the fact that they enter aquatic ecosystems via runoff from land subjected to drainage by man or as a consequence of other forms of human activity (Sakan et al., 2011; Waykar and Deshmukh, 2012). Owing to their toxicity and increased concentration in aquatic ecosystems, metals represent a growing problem on the global level, since once they enter an ecosystem they constitute a potential danger to living organisms for a number of years (Geffard et al., 2003; Fernandes et al., 2007; Snodgrass et al., 2008; Besser et al., 2009; Abdel-Baki et al., 2011), in as much as they are deposited both in the sediment (Fernandes et al., 2007; Abdel-Baki et al., 2011) and in tissues of aquatic organisms, where in elevated concentrations they have various toxic effects (Shahbaz et al., 2013: Ezejiofor et al., 2013). It is known that river sediments can receive high concentrations of heavy metals, which remain as difficultly degradable matter in the rivers (Fernandes et al., 2007; Abdel-Baki et al., 2011), and if the sediment contains a large amount of organic material it can bind metals in even greater measure (Zhang et al., 2014; Hsu et al., 2016; Strom et al., 2011). Increase in the amount of organic matter in a riveŕs sediment can occur under influence of the outlet water from a trout farm, which introduces remains of uneaten food and trout faeces into the recipient stream (Liao, 1970; Weston et al., 1996), with the result that metals accumulate in the sediment directly below the trout farm, enabling aquatic organisms to accumulate them in some measure (Burridge et al., 2010).

The list of significant polluters of aquatic environments includes Zn, Cu, Pb, Cd, Hg, Ni and Cr. Some of them like Cu, Zn and Cr are essential elements because they are co-factors of enzymes or a constituent of the respiratory pigment (haemocyanin). However, certain metals such as As, Hg and Cd in higher concentrations have a toxic effect. Their ions, even in low concentrations, can have serious effects on the health of living beings (Despotović, 2013). The degree of accumulation of metals in the body of organisms is affected both by physiological processes in the organisms themselves and by changes in physical and chemical parameters of the environment (Vicente-Martorell et al., 2009; Shulkin et al., 2003). Feeding also exerts considerable influence on metal accumulation in the body of vertebrates and represents one of the main sources of metals for many aquatic invertebrates (Wang, 2002). Thus, differences in the feeding of organisms of the benthofauna determine differences in the bioaccumulation of metals in the body of aquatic invertebrates (Clements, 1991; Pourang, 1996; Corbi et al., 2010).

Because of characteristics of their life cycle [confinement to the bottom, long duration of the life cycle, limited mobility (Hellawell, 1986; Bonada et al., 2006; De Pauw et al., 2006)], aquatic invertebrates are suitable bioindicators for monitoring the degree of accumulation of heavy metals both from the water and from the sediment (Fialkowski et al., 2003). Moreover, certain responses of aquatic organisms can indicate a danger of contamination by heavy metals in the earliest phases of pollution (Gremyatchikh et al., 2009), and monitoring changes at the individual level is more informative than monitoring changes at the community level for purposes of biomonitoring of aquatic ecosystems (Al-Shami et al., 2010). In view of what has been said and considering the fact that larvae of Ephemeroptera accumulate certain heavy metals in higher concentrations than some other taxa (Burrows and Whitton, 1983; Fialkowski et al., 2003) thanks to their specific morphological and functional characteristics (Cain et al., 1992), we chose the species Ephemera danica as a test organism. It is a good candidate for a model organism owing to its burrowing way of life, mode of feeding (an active gatherer of detritus filterers) and relatively large body dimensions (Clements and Kiffney, 1994), as well as because it can be collected in numbers sufficient for analysis (Nummelin et al., 2007). Also, its larval phase lasts one or two years, which makes it possible to monitor accumulation of heavy metals over a certain period of time (Svensson, 1977; Winkelmann and Koop, 2007; Bennett, 2007).

The present work is one piece of a comprehensive investigation of the trout farm effects on the various aspects of E. danica larvae biology. The first part of this investigation, dealing with the effects on the larvae antioxidative defense, has already being published (Božanić et al., 2018). In as much as river sediment burdened with organic pollution binds greater amounts of heavy metals than an unburdened recipient stream (Hsu et al., 2016; Burridge et al., 2010), the purpose of the present work was to establish whether an increase in the concentration of certain heavy metals occurs in water and sediment of the Skrapež River, a stream under influence of the outlet water from a trout farm, and ascertain the extent of metal bioaccumulation in the body of Ephemera danica larvae so as to determine their potential usefulness as a bioindicator of the pollution of river systems by heavy metals.

2 Material and methods

2.1 Description of investigated localities

Detailed description of the study area (Fig. 1 in Božanić et al., 2018), the trout farm and investigated localities is given in our previous publication (Božanić et al., 2018) therefore, only the basic data are given here. Located in western Serbia, the Skrapež River is 47.7 km long and has a drainage area of 647.65 km2 (Gavrilović and Dukić, 2002). The Kraj Vodenice (=Millside) Trout Farm is located along its course. Trout are reared in four independent basins on this fish farm, where they are given extruded feed (Skreting Optiline he 3p gal) with 0.9% phosphate content, 40% protein content and 28% fat content. Sampling was conducted at a control locality about 180 m upstream from the fish farm (locality SK1) and at three downstream localities (SK2, 30 m downstream from the fish farm; SK3, 300 m beyond locality SK2; and SK4, 300 m away from locality SK3).

2.2 Sampling of larvae and determination of heavy metal content

As already described in our previous publication (Božanić et al., 2018) specimens of E. danica larvae were collected seasonally (in April, July and October of 2015 and in January of 2016) at four localities on the Skrapež River by using a benthos sieve and tweezers. At each locality, 20 specimens of larvae were collected and fixed in the field with 96% alcohol in glass bottles. The samples were transported to a laboratory of the Environmental Protection Agency, where they were analyzed for the presence of heavy metals according to U.S. EPA Method 200.3. This method is applicable for analysis of the following elements: arsenic (As), cadmium (Cd), chromium (Cr), copper (Cu), iron (Fe), lead (Pb), nickel (Ni) and manganese (Mn).

Samples of E. danica larvae were first homogenized and a weighed amount (5 g) of a frozen tissue sample was then transferred to a 125-ml flask. A measured volume (10 ml) of concentrated nitric acid was added and the mixture warmed with gentle swirling of the samples on a hot plate until the tissue was solubilized (U.S. EPA Method 200.3). The temperature was increased to near boiling until the solution began to turn brown. The mixture was cooled, another batch (5 ml) of concentrated nitric acid was added to it and the sample was returned to the hot plate until the solution once again began to turn brown. It was then cooled and yet another batch (2 ml) of concentrated nitric acid was added, after which it was returned to the hot plate until its volume was reduced to 5–10 ml. After this, the mixture was again cooled, a measured volume (2 ml) of 30% hydrogen peroxide was added to it and the sample was returned to the hot plate until its volume was reduced to 5–10 ml. It was then cooled once more, a measured volume (2 ml) of concentrated hydrochloric acid was added to it and the sample was returned to the hot plate until its volume was reduced to 5 ml. The sample was then allowed to cool, after which it was quantitatively transferred to a 100-ml volumetric flask. Finally, the sample was diluted to volume with ASTM type I water and mixed so as to allow any insoluble material to separate (U.S. EPA Method 200.3). The concentration of metals in samples prepared in this way was determined by atomic absorption spectrophotmetry (AAS) on a Perkin Elmer Analyst 600 instrument.

2.3 Analysis of heavy metals in the water and sediment

Water samples for analysis of heavy metal content were taken 30 cm below the water surface using plastic bottles with a volume of 500 ml. Detailed procedure for sediment collection is described in Vranković et al. (2018). Analysis of heavy metals content in water and sediment was conducted at the Institute of Chemistry, Technology and Metallurgy of Belgrade University's Chemistry Centre. Samples of water and sediment were analysed for the presence of the following elements: As, Cd, Cr, Cu, Fe, Pb, Ni and Mn.

Water samples were stabilized by the addition of 0.5 mL HNO3 (70 wt.%, ACS reagent, Sigma Aldrich) and stored at 4 °C. Sediments were first homogenized, and the samples were then dried at 105 °C overnight. After this, sediment samples were powdered and passed through a sieve (<500 µm). Batches of samples weighing 0.5 g (accuracy ± 0.1 mg) were measured out into Teflon™ cuvettes and 15 ml of HNO3 (70 wt.%, ACS reagent, Sigma Aldrich) and HCl (37 wt.%, ACS reagent, Sigma Aldrich) (1:3, v/v) were added. Microwave digestion was done using the ETHOS 1 Advanced Microwave Digestion System (Milestone, Italy) for a period of 20 min at 200 °C. After digestion, the contents of the cuvette were diluted to 25 ml with rinsing of the precipitate (Božanić et al., 2019). Element concentrations were determined by the analytical technique of inductively coupled plasma optical emission spectrometry (ICP-OES) using a Thermo Scientific iCAP 6500 Duo ICP instrument (Thermo Fisher Scientific, Cambridge, UK). For calibration of the instrument, the following plasma standard solutions were used: SS-Low Level Elements ICV Stock (10 mg/L) and ILM 05.2 ICS Stock 1 (500 mg/L) (VHG Labs, Inc., LGC Standards, Manchester, NH 03103 USA); and Multi-Element Plasma Standard Solution 4, Specpure®, 1000 µg/ml (Alfa Aesar GmbH & Co., Germany). Quality control of the analytical process, performed using U.S. EPA Method 200.7 and LPC Solution Certified Reference Material (ULTRA Scientific, USA), indicated that the resulting concentrations were within 98–103%., The measurements were performed in triplicate (n = 3) for each digested sample.

2.4 Biota-sediment accumulation factor

Biota-sediment accumulation factor (BSAF) was used to quantify the capacity of E. danica larvae to accumulate heavy metals from the sediment (Mackay and Fraser, 2000). BSAF was calculated as the ratio of concentration of a metal in E. danica larvae (M E) and its concentration in the sediment (M s): BSAF = M E/M s.

2.5 Statistical analyses

Data were presented as mean ± standard error. Values of metal concentrations in the sediment, their concentrations in larvae of E. danica and BSAF values at the investigated localities were compared statistically using one-way ANOVA followed by the post hoc Fisher least significant difference (LSD) test (P < 0.05). The t-test was used to compare metal concentrations in the sediment and E. danica at the same locality. The existence of correlation between pairs of variables was tested with Pearson product moment correlation (P < 0.05). Statistical tests were performed with the aid of the Sigma Plot 11 software (Systat Software Inc., San Jose, CA, USA).

Co-inertia analysis (CIA) was used to describe the influence of metal concentrations in the sediment on their concentrations in E. danica larvae and appropriate BSAF values (Dolédec and Chessel, 1994). CIA was computed using the ADE-4 software (Thioulouse et al., 1997).

3 Results

In the present study, we investigated larvae of Ephemera danica for their ability to biologically accumulate in their bodies selected metals from the water and sediment in their characteristic aquatic environment, which is conspicuous for high concentrations of chromium and nickel in the sediment. In addition to this, we examined the influence of wastewater from a trout farm on concentrations of the tested heavy metals, both in water and sediment of the recipient stream and in larvae of E. danica dwelling in it.

Co-inertia analysis (CIA) was used to investigate the relationship between concentrations of heavy metals in the sediment on the one hand and their concentrations and capacity for accumulation in larvae of E. danica (measured by the biota sediment accumulation factor or BSAF) on the other, as well the influence of wastewater from the trout farm on these processes (Fig. 1). This analysis indicated the existence of a statistically highly significant (P < 0.001) co-structure between the PCA correlation matrix of heavy metal concentrations in the sediment and PCA correlation matrix of heavy metal concentrations in E. danica on the one hand, and BSAF on the other. In order to graphically represent the observed co-structure (Fig. 1), two factorial axes were retained, with the F1 axis explaining most of the co-structure (85.7%) and the F2 axis only small part of it (7.7%). Great strength of the co-structure was indicated by high values of the coefficients of correlation along the F1 axis (0.83) and F2 axis (0.87), as well as by short length of the arrows representing position of the localities.

thumbnail Fig. 1

Triplot presentation of the co-inertia analysis (CIA) of eight metal concentrations in the sediment and E. danica and their biota sediment accumulation factor performed on 16 samples [four sites, SK1-SK4, sampled in spring (sp), summer (su), autumn (a), and winter (w)]. The ordination diagram of eight normalized metal concentrations in the sediment (Ass, Cds, Crs,Cus, Fes, Mns, Nis, and Pbs) in the CIA is represented by grey arrows projected onto the F1 × F21 factorial map. Positions of metal concentrations in E. danica (AsE, CdE, Crs,CuE, FeE, MnE, NiE, and PbE) and their BSAF (AsB, CdB, CrB,CuB, FeB, MnB, NiB, and PbB) are also plotted onto the CIA F1 × F21 plane. Standardized co-inertia scores of sediment and E. danica data for each sample (black arrows) are projected onto the F1 × F22 factorial map. Origin of the arrow indicate the position of the sample as ordinated by metal concentrations in the sediment, while the arrowhead indicate its position as ordinated by the values of metal concentrations in E. danica and their BSAF.

3.1 Metal concentrations in the sediment

Analysis of metal concentrations in the sediment indicates that the vectors representing the concentrations of iron (Fes), copper (Cus), chromium (Crs), nickel (Nis), lead (Pbs) and cadmium (Cds) are aligned along the F1 axis and are oriented toward its negative end. Situated at the negative end of the F1 axis, the positions of all seasons of locality SK2 (origins of black arrows) which lies directly below the inflow of wastewater from the trout farm stand out clearly from positions of the other localities, which are placed either near the middle of the F1 axis (all seasons of the SK1 control locality and winter and summer of locality SK3) or at its positive end (all seasons of locality SK4 and autumn and spring of locality SK3). Along the F1 axis, which explains by far the greatest part of variation of the data, seasonal samples for each of the localities (with the exception of SK3) are very close together, which indicates that seasonal variations in concentrations of the aforementioned metals in the sediment are much less pronounced than longitudinal variations. This means that a clear increase in concentrations of the given metals occurs at the locality directly below the trout farm (SK2), the indicated increase being more significant than eventual natural longitudinal changes, but it is strictly localized and complete recovery occurs already at the next locality (SK3). For this reason, the F1 axis can be considered the axis of the trout farm's effect. One-way ANOVA indicated that Fes, Cus, Crs, Nis, Pbs and Cds increase statistically significantly at SK2 in relation to SK1 and then decline statistically significantly at SK3 and SK4, except in the case of the Pbs concentration, where a statistically significant decrease occurs just at S4 (Fig. 2B–E, G, H). Also, the concentrations of these metals in the sediment of SK4 are without exception lower than at S1, although that decrease is statistically significant only in the cases of Cus and Crs (Fig. 2B, D). There are no statistically significant differences between localities in the cases of the concentrations of manganese and arsenic in the sediment (Mns and Ass, respectively), and the value of Ass is even somewhat lower at SK2 than at the other localities. This is a consequence of seasonal variations greater than longitudinal ones, as is also indicated by CIA, where the vectors of these two metals are linked with the F2 axis, which represents the axis of seasonal variations.

thumbnail Fig. 2

Changes in average metal concentrations in the sediment (black bars) and E. danica (grey bars) and their BSAF (dark-grey bars) for arsenic (A), chromium (B), cadmium (C), copper (D), iron (E), manganese (F), nickel (G) and lead (H) at four investigated localities (SK1, SK2, SK3, and SK4) along the Skrapež River. Samples at different localities for each parameter that are significantly different are marked with different letters (a, b, c). A statistically significant difference between metal concentrations in the sediment and E. danica at each locality is marked with an asterisk (*).

3.2 Metal concentrations in larvae of E. danica

Analysis of the co-inertia of metal concentrations in larvae of E. danica indicates that the concentrations of chromium (CrE), copper (CuE), iron (FeE) and lead (PbE) just like the concentrations of these metals in the sediment are linked with the negative end of the F1 axis and are positively correlated with concentrations in the sediment (Fig. 1). As indicated by ANOVA, these metals are characterized by a statistically greater concentration at SK2 in relation to the other localities. Also, the concentrations of all four metals in E. danica show a statistically significant positive correlation with the concentration of the corresponding metal in the sediment: for CrE, R = 0.80 and P > 0.001 (Fig. 3A); for CuE, R = 0.72 and P = 0.002, (Fig. 3C); for FeE, R = 0.634 and P = 0.008; and for PbE, R = 0.57 and P = 0.021.

In co-inertia analysis, the concentration of cadmium in E. danica (CdE) is aligned counter to the vector of Cds, which indicates that their values change in opposite directions. Thus, ANOVA showed that CdE at SK2 is statistically greater than at SK1 and SK3 (Fig. 2C), while Pearson's test indicated that Cds and CdE are negatively correlated (R = –0.54, P = 0.030). On the CIA graph, the positions of concentrations of arsenic, nickel and manganese in E. danica (AsE, NiE and MnE, respectively) are situated virtually at a right angle to the vectors of concentrations of the given metals in the sediment (Fig. 1), and there is no statistically significant correlation between them. Only in the case of NiE does a statistically significant increase occur at SK2 in relation to SK1, while there are no significant differences between localities in the cases of AsE and MnE. If we compare the concentrations of a given metal in E. danica and in the sediment, the t-test indicates that a statistically significant difference exists in the cases of all metals except Pb and at each of the investigated localities, whereas in the case of Pb none exists at any of the localities (Fig. 2).

3.3 Biota sediment accumulation factor − BSAF

The values of BSAF indicate the direction of these differences. They are positive in the cases of arsenic (AsB), cadmium (CdB) and copper (CuB), which indicates that E. danica accumulates these metals in relation to the sediment (Fig. 2A, C, D). Accumulation is most pronounced in the case of arsenic, where AsB = 27 ± 7, whereas somewhat lower (but still significant) values were obtained for cadmium (CdB = 9.3 ± 1.5) and copper (CuB = 6.5 ± 0.9). In the case of lead, PbB is close to one (Fig. 2H), which means that its concentrations are the same in the sediment and in E. danica. As for the other metals, the value of this index is less than one (Fig. 2B, E, F, G), which means that they do not accumulate in larvae of E. danica. This is expressed especially strongly in the case of Ni, where the concentrations in E. danica are 60 times lower than in the sediment (NiB = 0.016 ± 0.001).

In analysis of the relationship between the concentrations of metals in the sediment and their BSAF, CIA indicates that CuB and FeB are linked with the negative end of the F1 axis (Fig. 1) and that they are characterized by a positive correlation with concentrations in the sediment (R = 0.53 and P = 0.036 for CuB, Fig. 3C; R = 0.51 and P = 0.044 for FeB). Moreover, ANOVA showed that the values of CuB and FeB attain a statistically significant maximum at SK2 (Fig. 2D, E). As for AsB, CrB and CdB, they are aligned counter to the vectors of their concentrations in the sediment and are characterized by a negative correlation (R = –0.80 and P < 0.001 for AsB; R = –0.77 and P < 0.001 for CrB, Fig. 3A; and R = –0.88 and P < 0.001 for Cdb, Fig. 3B).

3.4 Metal concentrations in the water

Co-inertia analysis (CIA) was used to analyze concentrations of the tested metals in the water and estimate their influence on the concentration and capacity for accumulation in E. danica, although this time it was shown that no significant (P = 0.7) co-structure exists between the PCA correlation matrix of heavy metal concentrations in the water and the PCA correlation matrix of heavy metal concentrations in E. danica and BSAF, so the analysis could not be completed (Tab. 1). This indicates that there is no significant influence of the concentration of heavy metals in the water and E. danica, a conclusion additionally confirmed by the absence of any statistically significant correlation between the concentration of metals in the water and E. danica [except in the case of iron (R = 0.52 and P = 0.037)], as well as by the absence of statistically significant variation of any of the metals at the investigated localities (Tab. 1).

Table 1

Average water concentrations of heavy metals (in mg/kg) at the localities of the Skrapež river.

4 Discussion

Determination of the total concentration of metals in the sediment and water does not give a precise estimate of their influence on the environment, and it is therefore necessary to monitor how much a certain element accumulates in organisms, i.e., how great is its bioavailabilty. Absorption of metals by aquatic invertebrates leads to accumulation of the metals in their bodies, which through the food chain further affects all consumers (including man) with potential long-term implications for human health (Notten et al., 2005; Fernandes et al., 2007; Corbi et al., 2008; Agah et al., 2009) and causes pollution of aquatic ecosystems (Malik et al., 2010). In order to understand the influence of a trout farm on accumulation of heavy metals in aquatic invertebrates, we investigated larvae of E. danica as a model organism.

4.1 Analysis of heavy metal concentrations in the sediment and in the body of E. danica larvae

Very high concentrations of Cr and Ni are a basic characteristic of the sediment at the investigated localities on the Skrapež River. According to published reference values for quality of the sediment in freshwater ecosystems (MacDonald et al., 2000), the lowest concentration of Ni is 7.4 times greater and that of Cr 2.45 times greater than PEC (the probable effect concentration, 48.6 and 111 mg/kg, respectively). A similar conclusion can be drawn if we compare the concentrations of these two metals with the maximum permissible concentrations (MPC) (Official Gazette No. 50/2012), which are significantly lower than the concentrations of these elements in the Skrapež sediment (the MPC for nickel is 44 mg/kg, while that for chromium is 240 mg/kg). Chromium is rarely found in natural waters (Mertz, 1986), but it can enter them and be deposited as a pollutant in the sediment as a result of various human activities. It is sometimes included as a micronutrient in fertilizers used in growing different cereals. According to the results obtained by Santos et al. (2009), soils are polluted with different concentrations of chromium depending on agricultural activity. Because the valley of the Skrapež is under agricultural crops, that is probably one of the reasons for the higher concentration of chromium in the riveŕs sediment. Chromium concentrations are higher at locality SK2 than at the other localities. There is more organic matter at that locality than at the other localities due to effluents from the trout farm, and the content of organic matter can have decisive influence on binding of chromium in the substrate, as was demonstrated by Santos et al. (2009) in their research. With respect to the concentration of nickel in the sediment, it is most likely of geochemical origin: the elevated concentration of this element in the sediment can be attributed to the fact that geological composition of the Skrapež basin includes magmatic rocks, metamorphic rocks and schists in which nickel is concentrated (Cannon et al., 1978; Kovačević-Majkić, 2009). Moreover, there are veins of nickel on Mt. Povlen, on which the Skrapež River arises (Kovačević-Majkić, 2009), and this is an additional reason for its high concentration. As in the case of Cr, an increase in the concentration of Ni occurs at the locality directly below the trout farm. This is probably a consequence of its elevated concentration in the outlet water due to the use of chemical substances in agricultural processes, which causes residues of the metal to be found in the raw materials used in production of trout feed.

For the other investigated metals in the Skrapež River, even the highest concentrations registered in the sediment are below their TEC (threshold effect concentration, MacDonald et al., 2000) values. That is especially clearly expressed in the cases of As and Cd, whose contamination factors are the lowest and whose highest registered concentrations are nine and 10 times lower, respectively, than the TEC (9.79 and 0.99 mg/kg respectively). Moreover, according to No. 50/2012 of the Official Gazette, the concentrations of all metals with the exceptions of Ni and Cr are far below the maximum permissible concentrations in surface waters.

In regard to the trout farḿs effect on heavy metal concentrations in the sediment, CIA (Fig. 1) and ANOVA (Fig. 2) clearly showed that the concentrations of Fes, Cus, Crs, Nis, Pbs and Cds at locality SK2 just below the trout farm attain a distinct and statistically significant maximum with very little seasonal variation, which indicates that the trout farm is the direct cause of this effect. This can be explained by the fact that release of effluents, which contain remains of trout feed and faeces (Kronvang et al., 1993; Iversen, 1995), from the trout farm into the recipient stream causes accumulation of the organic matter in the sediment immediately below the trout farm thus increasing its metal binding capacity (Zhang et al., 2014; Hsu et al., 2016; Burridge et al., 2010). Increase in the concentration of biogenic metals like Cu and Fe below the trout farm can also be attributed to their presence in trout feed (which is rich in trace elements of the iron and copper type), whose uneaten remains enter the recipient stream with the outlet water and accumulate in the sediment directly below the fish farm (Miller et al., 1993; Kamunde et al., 2002). Also, pollution of the environment with copper can arise through the use of copper-based pesticides, and copper salts in this way enter the water and then accumulate in the sediment (deOliviera-Filho et al., 2004; Flemming and Trevors, 1989).

As for concentrations of the tested metals in the body of E. danica larvae, co-inertia analysis and ANOVA showed that the concentrations of chromium, copper, iron, nickel and lead in the body of E. danica larvae (CrE, CuE, FeE, NiE, PbE) attain a statistically significant maximum at locality SK2. Except in the case of NiE, they exhibit a statistically significant positive correlation with concentrations of the same metals in the sediment (Fig. 3A and 3C). In contrast to them, CdE attains a minimum at SK2 and is characterized by a negative correlation with Cds (Fig. 3B). Previously published data on the relationship of heavy metal concentrations in the sediment and aquatic invertebrates are fairly contradictory. Thus, Harrahy and Clements (1997) noticed a significant increase of concentrations in the body of Chironomus tentans with increase of concentrations in the sediment for all metals (Cd, Zn and Pb) except copper. In keeping with this, in rivers receiving organic pollution it was established that the concentration of copper in Epobdella octoculata and Asellus aquaticus declines with an increase in its concentration in the sediment, while that of lead increases (Eyres and Pugh-Thomas, 1978). Meanwhile, Goodyear and McNeill (1999) assert that all invertebrates have a body concentration of copper in positive correlation with concentrations in the sediment. Hare et al. (1991) found Cu and Pb concentrations in the body of Hexagenia to be in proportion with those in the sediment. Corbi et al. (2010) state that iron is accumulated in insects in correlation with its concentrations in the sediment. Jop (1991) monitored the concentrations of cadmium in E. danica, E. vulgata, Leptophlebia vespertina and Baetis vernus and contrary to our results concluded that the concentration of cadmium in them is proportional to concentrations in the sediment. Contrary to this, Corbi et al. (2010) report that Cd was not registered in insect larvae even when present in the sediment.

thumbnail Fig. 3

Effects of metal concentration in the sediment on metal concentration in E. danica and its BSAF. (A) Effect of changes in Cr concentration in the sediment (Crs) during the investigation period on its concentration in E. danica (CrE, black circles) and BSAF (CrB, open triangles) at the investigated localities on the Skrapež River. The black solid line represents the linear fit of the data in the case of CrE (R = 0.80, P < 0.001) and the grey one in the case of CrB (R = –0.77, P < 0.001). (B) Effect of changes in Cd concentration in the sediment (Cds) during the investigation period on its concentration in E. danica (CdE, black circles) and BSAF (CdB, open triangles) at the investigated localities on the Skrapež River. The black solid line represents the linear fit of the data in the case of CdE (R = –0.54, P = 0.030) and the grey one in the case of CdB (R = –0.88, P < 0.001). (C) Effect of changes in Cu concentration in sediment (Cus) during the investigation period on its concentration in E. danica (CuE, black circles) and BSAF (CuB, open triangles) at the investigated localities on the Skrapež River. The black solid line represents the linear fit of the data in the case of CuE (R = 0.72, P = 0.002) and the grey one in the case of CdB (R = 0.53, P = 0.036).

4.2 Biota sediment accumulation factor − BSAF

Every aquatic invertebrate assimilates metals from its surroundings, such assimilation depending upon the species of invertebrate and the relative biological capabilities of the metal (Rainbow and Wang, 2001). Factors that influence the fate of a metal are physiology of the invertebrate, whether the metal is used for an essential physiological function, whether it is excreted or stored in the body, and whether it gains access to a “wrong” biomolecule and thereby has a toxic effect (Rainbow, 2002). Our research indicates that metals characterized by a high concentration in the sediment (Ni, Cr, Fe and Mn) have a BASF value of less than one, that is their bioaccumulation is absent in E. danica, whereas elements characterized by a low concentration in the sediment have a BASF value greater than one, which is to say that their bioaccumulation occurs in E. danica. Meanwhile, if the relationship of concentrations of these elements in the sediment and their BSAF values is examined within the framework of these two groups, it is possible to discern further distinctions, namely that Fe and Cu are characterized by a positive correlation of these two values, whereas Cr, As and Cd are characterized by a negative one. Moreover, CuB and FeB attain a statistically significant maximum at SK2 (Fig. 2D, E). Only Pb, whose concentration in the sediment is also relatively low, is characterized by a BSAF value of about one. In the cases of Ni and Cr, which are toxic for living beings (Poulton et al., 1989; Eisler, 1998; Kotaś and Stasicka, 2000) and which are present in the sediment in exceptionally high concentrations, such a result could be expected and probably represents a defense mechanism of E. danica against assimilation of these metals in high concentrations lethal for it. This is confirmed by the aforementioned negative correlation between CrB and Crs (Fig. 3A). Furthermore, in comparison with the other metals, nickel and chromium are less mobile, which means that they are less available to organisms of the bottom fauna and have smaller chances of entering the food chain (Leonard and Wood, 2013). Our results are in agreement with the findings of Shulkin et al. (2003), who reported an increase in concentration of all metals except nickel in the body of the shellfish Crenomitllis grayanus. Also, De Jonge et al. (2014) noted that burdening of the body in Leuctra sp., Rhithrogena sp., Simuliidae and Perlodidae is not linked with the level of total Ni. Watras et al. (1985) measured the concentration of nickel in Daphnia magna and noticed that ingested nickel either is not incorporated or else is quickly eliminated. Muyssen et al. (2004) concluded that the BSAF value for Ni is almost always less than one, i.e., the concentrations of nickel in benthic organisms are lower than concentrations in the sediment. In the case of chromium, Hamidian et al. (2016) noticed that slight bioaccumulation of this metal occurs in chironomids in the Shoor River (Iran), but its concentrations in the sediment were 5–10 times lower than in the Skrapež River, which in view of the negative correlation between CrB and Crs observed in our study can explain such a result in the Shoor River. All of this supports our conclusion that aquatic invertebrates probably develop defense mechanisms against this toxic metal, and that the greater its content in the sediment, the less it is assimilated by benthic organisms.

The positive correlation between accumulation of copper and iron in the sediment and in the body of E. danica can be attributed to the fact that these metals are essential and necessary for metabolic activity of the organism (Cohen et al., 2001). As has already been said, Goodyear and McNeill (1999) in their study also found a positive correlation between the concentration of copper in the investigated invertebrates and its concentration in the sediment. Moreover, they concluded that copper concentrations in the organisms increase in proportion with increase of copper concentration in the sediment, without signs of biological limitation. These results were corroborated by Gundacker (2000), who compared the concentration of metals in the sediment and in the bodies of shellfish and concluded that bioaccumulation of copper occurred in them, in contrast to lead. It can be asserted that the situation is similar with iron, since Corbi et al. (2010) found the concentrations of iron in the sediment to be very high and in positive correlation with its accumulation in insects, especially in species of the genus Chironomus.

As for bioaccumulation of Cd and As, a negative correlation is evident between the levels of these elements in the sediment and in the body of E. danica. To be specific, although the level of bioaccumulation of Cd and especially As is very high, its value in both cases decreases with increase in concentrations of the elements in the sediment. In view of the fact that we are dealing here with exceptionally toxic substances, this negative correlation probably represents a defense mechanism of E. danica against their excessive accumulation. In the case of the level of Cd bioaccumulation, our results corroborate those obtained in the investigations Girgin et al. (2010), where it is stated that representatives of the order Ephemeroptera accumulate greater concentrations of cadmium than do members of other orders of aquatic invertebrates. In addition, Gundacker (2000) reported a significant level of Cd bioaccumulation in shellfish. On the other hand, research on bioaccumulation of heavy metals in chironomids in the Shoor (Iran) indicated that it is lacking in the cases of As and Cd (Hamidian et al., 2016). In the case of Cd, that can be explained by the fact that its concentration in the sediment of the Shoor River was below the limit of detection, while the concentration of As in the sediment was more than 10 greater than in our investigation, which in view of the observed negative correlation between AsB and Ass can explain the absence of bioaccumulation in the Shoor River.

4.3 Analysis of concentrations of heavy metals in the water in relation to concentrations in E. danica

All of the analyses performed in the present study unambiguously indicate that concentrations of the investigated metals in the water have no significant influence on their concentrations in larvae of E. danica. The absence of influence could be a consequence of the way of life of these organisms. As larvae they burrow into the substrate making tunnels in the process, and feed by filtering organic detritus (Clements and Kiffney, 1994). For this reason, there exists a positive correlation between concentrations in the sediment and in the body of the larvae, but not between the latter and metal concentrations in the water. Similarly, Bervoets et al. (1997) indicated that the concentrations of lead and cadmium in bodies of the organisms they investigated are in correlation with total lead concentration in the sediment, but not with their concentration in poor water. In contrast to our results, in the study of Brown (1977) who investigated a river polluted by waste materials from an abandoned mine a significant positive correlation was recorded between the concentration of copper in the body of Trichoptera larvae and in the water. Also, a correlation was found to exist between the concentrations of one or more of the studied metals in the bodies of Ecdyonurus venosus, Brachyptera risi, Leuctra spp., Perla bipunctata, Rhyacophila dorsalis and Dicranota sp. and their concentrations in the water (Burrows and Whitton, 1983).

4.4 Effects of metals on the activity of antioxidative enzymes in E. danica larvae

In parallel with this investigation we used the same system (Skrapež River and The Kraj Vodenice Trout Farm) to monitor the effects of changes in water chemistry on the activity of enzymes superoxide dismutase (SOD) and glutathione peroxidase (GPx) and concentration of total glutathione (GSH) in E. danica larvae (Božanić et al., 2018). SOD, GPx and GSH are important constitutens of the anioxidant defence system (Valko et al., 2016). It is well documented that one of the mechanisms of metal-induced toxicity is the generation of reactive oxygen species (ROS) (Valko et al., 2016). For example, metals, including Fe, Cu, Cr, and V undergo redox cycling, while As, Cd, Hg, Ni, and Pb, deplete glutathione and protein-bound sulfhydryl groups, resulting in the production of ROS and altering the activity of antioxidative enzymes responsible for ROS degradation (Stohs and Bagchi, 1995; Valko et al., 2005; Gomes et al., 2013).

Comparison of metal concentrations in the sediment with the activities of SOD and GPx and GSH concentration in E. danica larvae showed no statisticaly significant correlation (data not shown) which is in agreement with results obtained on Gammarus dulensis in the Crnica River (Vranković et al., 2018). However, GPx activity in E. danica larvae was stimulated by the increase in NiE (R = 0.746, P < 0.001) and PbE (R = 0.617, P = 0.011) and the decrease in CdE (R = –0.702, P = 0.002), indicating that apart from the concentrations of NO3 , NH4 + and value of DO% in the water of the Skrapež River, NiE PbE and CdE could olso influence GPx activity. As expected, metals accumulated in the larval body exerted greater influence on the antioxidant defence system compared to those in sediment.

5 Conclusions

This is a rear study investigating the influence of the trout farm effluents on the accumulation of heavy metals in the sediment of the recipient stream, and macrozoobenthos organisms. We clearly showed that trout farm discharge induced localized and statistically significant increase of concentrations of Fe, Cu, Cr, Ni, Pb and Cd in sediment and Ephemera danica larvae, but not in water indicating that sediment is the main source of heavy metals for this organism. In addition, evaluation of E. danica as a potential bioindicator of heavy metal pollution showed that only metals with low values of contamination factor (As, Cd and Cu) accumulated in the larvae. Moreover, values of Biota sediment accumulation factor for toxic metals (As, Cd and Cr) showed negative correlation with their concentrations in sediment indicating existence of defense mechanisms in E. danica against assimilation of these metals in high concentrations. In conclusion E. danica larvae proved to be good bioindicators of the pollution of river systems by As, Cd and Cu.

Acknowledgements

The study was supported by projects of the Ministry of Education, Science and Technological Development of the Republic of Serbia Project No. TR31075 and OI172030.

References

  • Abdel-Baki AS, Dkhil MA, Al-Quraishy S. 2011. Bioaccumulation of some heavy metals in tilapia fish relevant to their concentration in water and sediment of Wadi Hanifah, Saudi Arabia. Afr J Biotechnol 10: 2541–2547. [Google Scholar]
  • Agah H, Leermakers M, Elskens M, Fatemi SMR, Baeyens W. 2009. Accumulation of trace metals in the muscles and liver tissues of five fish species from the Persian Gulf. Environ Monit Assess 157: 499–514. [CrossRef] [PubMed] [Google Scholar]
  • Al-Shami S, Rawi C, Nor S, Ahmad A, Ali A. 2010. Morphological Deformities in Chironomus spp. (Diptera:Chironomidae) Larvae as a Tool for Impact Assessment of Anthropogenic and Environmental Stresses on Three Rivers in the Juru River System, Penang, Malaysia. Environ Entomol 39: 210–222. [CrossRef] [PubMed] [Google Scholar]
  • Bennett C. 2007. A seven year of the life cycle of the mayfly Ephemera danica . Freshwater Forum 27: 3–14. [Google Scholar]
  • Bervoets L, Blust R, De Wit M, Verheyen R. 1997. Relationships between river sediments characteristics and trace metal concentrations in tubificid worms and chironomid larvae. Environ Pollut 95: 345–356. [Google Scholar]
  • Besser J, Brumbaugh W, Allert A, Poulton B, Schmitt C, Ingersoll CG. 2009. Ecological impacts of lead mining on Ozark streams: toxicity of sediment and pore water. Ecotoxicol Environ Saf 72: 516–526. [CrossRef] [PubMed] [Google Scholar]
  • Bonada N, Prat N, Resh HV, Statzner B. 2006. Developments in aquatic insect biomonitoring: a comparative analysis of recent approaches. Annu Rev Entomol 51: 495–523. [CrossRef] [PubMed] [Google Scholar]
  • Božanić M, Todorović D, Živić M, Perić-Mataruga V, Marković Z, Živić I. 2018. Influence of a trout farm on antioxidant defense in larvae of Ephemera danica (Insecta: Ephemeroptera). Knowl Manag Aquat Ecol 47: 419–447. [Google Scholar]
  • Božanić M, Marković Z, Živić M, et al. 2019. Mouthpart deformities of Chironomus plumosus larvae caused by increased concentrations of copper in sediment from carp fish pond. Turk J Fish Aquat Sci 19: 251–259. [Google Scholar]
  • Brown BE. 1977. Effects of mine drainage on the River Hayle, Cornwall, factors affecting concentrations of copper, zinc and iron in water, sediments and dominant invertebrate fauna. Hydrobiologia 52: 221–233. [Google Scholar]
  • Burridge L, Weis JS, Cabello F, Pizarro J, Bostick K. 2010. Chemical use in salmon aquaculture: a review of current practices and possible environmental effects. Aquaculture 306: 7–23. [Google Scholar]
  • Burrows G, Whitton BA. 1983. Heavy metals in water, sediments and invertebrates from a metal-contaminated river free of organic pollution. Hydrobiologia 106: 263–273. [Google Scholar]
  • Cain DJ, Luoma SN, Carter JL, Fend SV. 1992. Aquatic insects as bioindicators of trace element contamination in cobble-bottom rivers and streams. Can J Fish Aquat Sci 49: 2141–2154. [Google Scholar]
  • Cannon HL, Connally GG, Epstein JB, Parker JG, Thorton I, Wixson G. 1978. Rocks: geological sources of most trace elements. In: Report to the workshop at southscas plantation Captiva Island, FL, US. Geochemistry and the Environment 3: 17–31. [Google Scholar]
  • Çevik F, Göksu M, Derici O, Findik O. 2009. An assessment of metal pollution in surface sediments of Seyhan dam by using enrichment factor, geoaccumulation index and statistical analyses. Environ Monit Assess 152: 309–317. [Google Scholar]
  • Clements W, Kiffney P. 1994. Integrated laboratory and field approach for assessing impacts of heavy metals at the Arkansas river, Colorado. Environ Toxicol Chem 13: 397–404. [Google Scholar]
  • Clements WH. 1991. Community responses of stream organisms to heavy metals: a review of descriptive and experimental approaches. In: Newman MC, McIntosh AW, ed. Ecotoxicology of metals: current con-cepts and applications, Lewis Publishers, Chelsea, 363–391. [Google Scholar]
  • Cohen T, Hee S, Ambrose R. 2001. Trace metals in fish and invertebrates of three California Coastal Wetlands. Mar Pollut Bull 42: 232–242. [Google Scholar]
  • Corbi JJ, Froehlich CG, Trivinho-Strixino S, Dos Santos A. 2010. Bioaccumulation of metals in aquatic insects of streams located in areas with sugar cane cultivation. Químca Nova 33:644–648. [CrossRef] [Google Scholar]
  • Corbi JJ, Trivinho-Strixino S, Dos Santos A. 2008. Environmental evaluation of metals in sediments and dragonflies due to sugar cane cultivation in Neotropical streams. Water Air Soil Pollut 195: 325–333. [Google Scholar]
  • De Almeida EA, Bainy ACD, De Melo Loureiro AP, et al. 2007. Oxidative stress in Perna perna and other bivalves as indicators of environmental stress in the Brazilian marine environment: antioxidants, lipid peroxidation and DNA damage. Comp Biochem Physiol Part A 146: 588–600. [CrossRef] [Google Scholar]
  • De Jonge M, Lofts S, Bervoets L, Blust R. 2014. Relating metal exposure and chemical speciation to trace metal accumulation in aquatic insects under natural field conditions. Sci Total Environ 496: 11–21. [CrossRef] [PubMed] [Google Scholar]
  • De Pauw N, Gabriëls W, Goethals P. 2006. River monitoring and assessment methods based on macroinvertebrates. In Ziglio G, Siligardi M, Flaim G, eds. Biological monitoring of rivers: applications and perspectives. Chichester, UK: John Wiley & Sons, 113–134. [Google Scholar]
  • Despotović S. 2013. Parameters of the antioxidant defence system and heavy metal concentrations in the visceral mass of selected snail and mussel species from the Danube, Tisa and Velika Morava rivers. Doctoral dissertation, University of Belgrade, 137 p. [Google Scholar]
  • Dolédec S, Chessel D. 1994. Co-inertia analysis: an alternative method for studying species-environment relationships. Freshw Biol 31: 277–294. [Google Scholar]
  • Eisler R. 1998. Nickel Hazards to Fish, Wildlife, and Invertebrates: A Synoptic Review. U.S. Patuxent Wildlife Research Center, Washington DC, 76 p. [Google Scholar]
  • Eyres JP, Pugh-Thomas M. 1978. Heavy metal pollution of the River Irwell (Lancashire, UK) demonstrated by analysis of substrate material and macroinvertebrate tissue. Environ Pollut 16: 129–136. [Google Scholar]
  • Ezejiofor TIN, Ezejiofor AN, Udebuani AC, et al. 2013. Environmental metals pollutants load of a densely populated and heavily industrialized commercial city of Aba, Nigeria. J Toxicol Environ Health Sci 5: 1–11. [CrossRef] [Google Scholar]
  • Fernandes C, Fontainhas-Fernandes A, Peixoto F, Salgado MA. 2007. Bioaccumulation of heavy metals in Liza saliens from the Esomriz-Paramos coastal lagoon, Portugal. Ecotoxicol Environ Saf 66: 426–431. [CrossRef] [PubMed] [Google Scholar]
  • Fialkowski W, Klonowska-Olejnika M, Smith BD, Rainbow PS. 2003. Mayfly larvae (Baetis rhodani and B. vernus) as biomonitors of trace metal pollution in streams of a catchment draining a zinc and lead mining area of Upper Silesia, Poland. Environ Pollut 121: 253–267. [Google Scholar]
  • Flemming CA, Trevors JT. 1989. Copper toxicity and chemistry in the environment: a review. Water Air Soil Pollut 44: 143–158. [Google Scholar]
  • Gavrilović Lj, Dukić D. 2002. River of Serbia. Institute for textbook publishing and teaching aids, Belgrade, 218 p. (in Serbian). [Google Scholar]
  • Geffard O, Geffard A, His E, Budzinski H. 2003. Assesment of the bioavaliability and toxicity of sediment associated polycyclic aromatic hydrocarbons and heavy metals applied to Crassostrea gigas embryos and larvae. Mar Pollut Bull 46: 481–490. [Google Scholar]
  • Girgin S, Kazancı N, Dügel M. 2010. Relationship between aquatic insects and heavy metals in an urbanstream using multivariate techniques. Int J Environ Sci Technol 7: 653–664. [CrossRef] [Google Scholar]
  • Gomes T, Gonzalez-Rey M, Rodríguez-Romero A, et al. 2013. Biomarkers in Nereis diversicolor (Polychaeta: Nereididae) as management tools for environmental assessment on the southwest Iberian coast. Sci Mar 77: 69–78. [Google Scholar]
  • Goodyear KL, McNeill S. 1999. Bioaccumulation of heavy metals by aquatic macro-invertebrates of different feeding guilds: a review. Sci Total Environ 229: 1–19. [Google Scholar]
  • Gremyatchikh V, Tomilina II, Grebenyuk LP. 2009. The effect of mercury chloride on morphofunctional parameters in Chironomus riparius Meigen (Diptera, Chironomidae) larvae. Inland Water Biol 1: 89–95. [CrossRef] [Google Scholar]
  • Gundacker C. 2000. Comparison of heavy metal bioaccumulation in freshwater molluscs of urban river habitats in Vienna. Environ Pollut 110: 61–71. [Google Scholar]
  • Hamidian AH, Zareh M, Poorbagher H, Vaziri L, Ashrafi S. 2016. Heavy metal bioaccumulation in sediment, common reed, algae, and blood worm from the Shoor river, Iran. Toxicol Ind Health 32: 398–409. [CrossRef] [PubMed] [Google Scholar]
  • Hare L, Tessier A, Campbell PGC. 1991. Trace element distribution in aquatic insects: variations among genera, elements, and lakes. Can J Fish Aquat Sci 48: 1481–1491. [Google Scholar]
  • Harrahy EA, Clements WH. 1997. Toxicity and bioaccumulation of a mixture of heavy metals in Chironomus tentans (Dipter: Chironomidae) in synthetic sediment. Environ Toxicol Chem 16: 317–327. [Google Scholar]
  • Hellawell JM. 1986. Biological Indicators of Freshwater Pollution and Environmental Management. London & New York: Elsevier Applied Science Publishers, 518 p. [Google Scholar]
  • Hsu L-C, Huang C-Y, Chuang Y-H, et al. 2016. Accumulation of heavy metals and trace elements in fluvial sediments received effluents from traditional and semiconductor industries. Sci Rep 6: 34250. [CrossRef] [PubMed] [Google Scholar]
  • Iversen TM. 1995. Fish farming in Denmark: Environmental impact of regulative legislation. Water Sci Technol 31: 73–84. [Google Scholar]
  • Jop KM. 1991. Concentration of metals in various larval stages of four Ephemeroptera species. Bull Environ Contam Toxicol 46: 901–905. [CrossRef] [PubMed] [Google Scholar]
  • Kamunde C, Grosell M, Higgs D, Wood C. 2002. Copper metabolism in actively growing rainbow trout (Onchorhynchus mykiss): interactions between dietary and waterborne copper uptake. J Exp Biol 205: 279–290. [PubMed] [Google Scholar]
  • Kotaś J, Stasicka Z. 2000. Chromium occurence in the environment and methods of it speciation. Environ Pollut 107: 263–283. [Google Scholar]
  • Kovačević-Majkić J. 2009. Hydrogeographic study of the Skrapež River. Geographical Institute “Jovan Cvijić” SANU, Serbian Geographical Society, Belgrade, Special Issues 74: 133 p. (in Serbian) [Google Scholar]
  • Kronvang B, Ertebjerg G, Grant R, Kristensen P, Hovmand M, Kirkegard J. 1993. Nationwide monitoring of nutrients and their ecological effects: state of the Danish aquatic environmental. Ambio 22: 176–187. [Google Scholar]
  • Leonard EM, Wood CM. 2013. Acute toxicity, critical body residues, Michaelis–Menten analysis of bioaccumulation, and ionoregulatory disturbance in response to waterborne nickel in four invertebrates: Chironomus riparius, Lymnaea stagnalis, Lumbriculus variegatus and Daphnia pulex . Comp Biochem Physiol Part C 158: 10–21. [Google Scholar]
  • Liao PB. 1970. Pollution potential of salmonid fish hatcheries. Water Sewage Works 117: 291–297. [Google Scholar]
  • MacDonald DD, Ingersoll CG, Berger TA. 2000. Development and Evaluation of Consensus-Based Sediment Quality Guidelines for Freshwater Ecosystems. Arch Environ Contam Toxicol 39: 20–31. [CrossRef] [PubMed] [Google Scholar]
  • Mackay D, Fraser A. 2000. Bioaccumulation of persistent organic chemicals: mechanisms and models. Environ Pollut 110: 375–391. [Google Scholar]
  • Malik N, Biswas AK, Qureshi TA, Borana K, Virha R. 2010. Bioaccumulation of heavy metals in fish tissues of a freshwater lake of Bhopal. Environ Monit Assess 160: 267–276. [CrossRef] [PubMed] [Google Scholar]
  • Mertz W. 1986. Trace Elements in Human and Animal Nutrition. London: Academic Press, 560 p. [Google Scholar]
  • Miller P, Lanno R, McMaster M, Dixon D. 1993. Relative contributions of dietary and waterborne copper to tissue copper burdens and waterborne copper tolerance in rainbow trout (Onchorhynchus mykiss). Can J Fish Aquat Sci 50: 1683–1689. [Google Scholar]
  • Muyssen BTA, Brix KV, DeForest DK, Janssen CR. 2004. Nickel essentiality and homeostasis in aquatic organisms. Environ Rev 12: 113–131. [Google Scholar]
  • Notten MJM, Oosthoek AJP, Rozema J, Aerts R. 2005. Heavy metal concentrations in a soil-plant-snail food chain along a terrestrial soil pollution gradient. Environ Pollut 138: 178–190. [Google Scholar]
  • Nummelin M, Lodenius M, Tulisalo E, Hirvonen H, Alanko T. 2007. Predatory insects as bioindicators of heavy metal pollution. Environ Pollut 145: 339–347. [Google Scholar]
  • Official Gazette of the Republic of Serbia 50/2012 (2012) Regulation on limit values of pollutants in surface and groundwaters and sediments and deadlines for achieving them. Accessed 18 May 2012. (in Serbian) [Google Scholar]
  • Poulton BC, Beitinger TL, Stewart KW. 1989. The Effect of Hexavalent Chromium on the Critical Thermal Maximum and Body Burden of Clopper Clio (Plecoptera: Perlodidae). Arch Environ Contam Toxicol 18: 594–600. [Google Scholar]
  • Pourang N. 1996. Heavy metal concentrations in surficial sediments and benthic macroinvertebrates from Anzali wetland, Iran. Hydrobiologia 331: 53–61. [Google Scholar]
  • Rainbow PS, Wang WX. 2001. Comparative assimilation of Cr, Cr, Se, and Zn by the barnacle Elminius modestus from phytoplankton and zooplankton diets. Mar Ecol Prog Ser 218: 239–248. [Google Scholar]
  • Rainbow PS. 2002. Trace metal concentrations in aquatic invertebrates: why and so what? Environ Pollut 120: 497–507. [Google Scholar]
  • Sakan S, Đorđević D, Dević G, Relić D, Anđelković I, Đuričić J. 2011. A study of trace element contamination in river sediments in Serbia using microwave-assisted aqua regia digestion and multivariate statistical analysis. Microchem J 99: 492–502. [Google Scholar]
  • Santos A, Oliveira LC, Botero WG, et al. 2009. Distribuição e biodisponibilidade de crômio em solos contaminados por resíduos de couro. Quimica Nova 7: 1693–1697. [Google Scholar]
  • Schueler T. 2000. Cars are leading source of metal loads in California. The Practice of Watershed Protection, Center for Watershed Protection, Ellicott City, MD. pp. 44. [Google Scholar]
  • Shahbaz M, Hashmi MZ, Malik RN, Yasmin A. 2013. Relationship between heavy metals concentrations in egret species, their environment and food chain differences from two Headworks of Pakistan. Chemosphere 93: 274–282. [CrossRef] [PubMed] [Google Scholar]
  • Shulkin VM, Presleyb BJ, Kavunc VI. 2003. Metal concentrations in mussel Crenomytilus grayanus and oyster Crassostrea gigas in relation to contamination of ambient sediments. Environ Int 29: 493–502. [CrossRef] [PubMed] [Google Scholar]
  • Snodgrass J, Casey R, Joseph D, Simon J. 2008. Microcosm investigationsof storm water pond sediment toxicity to embryonic and larval amphibians: variation in sensitivity among species. Environ Pollut 154: 291–297. [Google Scholar]
  • Stohs S, Bagchi D. 1995. Oxidative mechanisms in the toxicity of metal ions. Free Radic Biol Med 18: 321–336. [CrossRef] [PubMed] [Google Scholar]
  • Strom D, Simpson SL, Batley GE, Jolley DF. 2011. The influence of sediment particle size and organic carbon on toxicity of copper to benthic invertebrates in oxic/suboxic surface sediments. Environ Toxicol Chem 30: 1599–1610. [CrossRef] [PubMed] [Google Scholar]
  • Svensson B. 1977. Life cycle, energy fluctuations and sexual differentiation in Ephemera danica (Ephemeroptera), a stream-living mayfly. Oikos 29: 78–86. [Google Scholar]
  • Thioulouse J, Chessel D, Dolédec S, Olivier JM. 1997. ADE-4: a multivariate analysis and graphical display software. Stat Comput 7: 75–83. [Google Scholar]
  • Triquet-Amiard C, Amiard J-C, Rainbow PS. 2013. Ecological Biomarkers: Indicators of Ecotoxicological Effects. CRC Press, Taylor & Francis Group. 450 p. [Google Scholar]
  • USEPA. 1994. Method 200.7. Determination of Metals and Trace Elements in Water and Wastes by Inductively Coupled Plasma-Atomic Emission Spectrometry. Washington D.C. [Google Scholar]
  • USEPA. 1996. Method 200.3. Sample Preparation Procedure for Spectrochemical Determination of Total Recoverable Elements in Biological Tissues. Washington D.C. [Google Scholar]
  • Valko M, Jomova K, Rhodes CJ, Kuča K, Musílek K. 2016. Redox- and non-redox-metal-induced formation of free radicals and their role in human disease. Arch Toxicol 90: 1–37. [CrossRef] [PubMed] [Google Scholar]
  • Vicente-Martorell JJ, Galindo-Riaño MD, García-Vargas M, Granado-Castro MD. 2009. Bioavailability of heavy metals monitoring water, sediments and fish species from a polluted estuary. J Hazard Mater 162: 823–836. [Google Scholar]
  • Vranković J, Živić M, Radojević A, et al. 2018. Evaluation of oxidative stress biomarkers in the freshwater gammarid Gammarus dulensis exposed to trout farm outputs. Ecotox Environ Safe 163: 84–95. [CrossRef] [PubMed] [Google Scholar]
  • Wang WX. 2002. Interactions of trace metals and different marine food chains. Mar Ecol Prog Ser 243: 295–309. [Google Scholar]
  • Watras CJ, MacFarlane J, Morel FMM. 1985. Nickel Accumulation by Scenedesmus and Daphnia: Food-Chain Transport and Geochemical Implications. Can J Fish Aquat Sci 42: 724–730. [Google Scholar]
  • Waykar B, Deshmukh G. 2012. Evaluation of bivalves as bioindicators of metal pollution in freshwater. Bull Environ Contam Toxicol 88: 48–53. [CrossRef] [PubMed] [Google Scholar]
  • Weston DP, Phillips MJ, Kelly LA. 1996. Environmental impacts of salmonid culture. Dev Aquacult Fish Sci 29: 919–967. [Google Scholar]
  • Winkelmann C, Koop JHE. 2007. The management of metabolic energy storage during the life cycle of mayflies: a comparative field investigation of the collector-gatherer Ephemera danica and the scraper Rhithrogena semicolorata . J Comp Physiol B 177: 119–128. [CrossRef] [PubMed] [Google Scholar]
  • Zhang C, Yu Z, Zeng G, et al. 2014. Effects of sediment geochemical properties on heavy metal bioavailability. Environ Internat 73: 270–281. [CrossRef] [Google Scholar]

Cite this article as: Božanić ML, Dojčinović BP, Živić MŽ, Marković ZZ, Manojlović DD, Živić IM. 2019. Bioaccumulation of heavy metals in Ephemera danica larvae under influence of a trout farm outlet waters. Knowl. Manag. Aquat. Ecosyst., 420, 50.

All Tables

Table 1

Average water concentrations of heavy metals (in mg/kg) at the localities of the Skrapež river.

All Figures

thumbnail Fig. 1

Triplot presentation of the co-inertia analysis (CIA) of eight metal concentrations in the sediment and E. danica and their biota sediment accumulation factor performed on 16 samples [four sites, SK1-SK4, sampled in spring (sp), summer (su), autumn (a), and winter (w)]. The ordination diagram of eight normalized metal concentrations in the sediment (Ass, Cds, Crs,Cus, Fes, Mns, Nis, and Pbs) in the CIA is represented by grey arrows projected onto the F1 × F21 factorial map. Positions of metal concentrations in E. danica (AsE, CdE, Crs,CuE, FeE, MnE, NiE, and PbE) and their BSAF (AsB, CdB, CrB,CuB, FeB, MnB, NiB, and PbB) are also plotted onto the CIA F1 × F21 plane. Standardized co-inertia scores of sediment and E. danica data for each sample (black arrows) are projected onto the F1 × F22 factorial map. Origin of the arrow indicate the position of the sample as ordinated by metal concentrations in the sediment, while the arrowhead indicate its position as ordinated by the values of metal concentrations in E. danica and their BSAF.

In the text
thumbnail Fig. 2

Changes in average metal concentrations in the sediment (black bars) and E. danica (grey bars) and their BSAF (dark-grey bars) for arsenic (A), chromium (B), cadmium (C), copper (D), iron (E), manganese (F), nickel (G) and lead (H) at four investigated localities (SK1, SK2, SK3, and SK4) along the Skrapež River. Samples at different localities for each parameter that are significantly different are marked with different letters (a, b, c). A statistically significant difference between metal concentrations in the sediment and E. danica at each locality is marked with an asterisk (*).

In the text
thumbnail Fig. 3

Effects of metal concentration in the sediment on metal concentration in E. danica and its BSAF. (A) Effect of changes in Cr concentration in the sediment (Crs) during the investigation period on its concentration in E. danica (CrE, black circles) and BSAF (CrB, open triangles) at the investigated localities on the Skrapež River. The black solid line represents the linear fit of the data in the case of CrE (R = 0.80, P < 0.001) and the grey one in the case of CrB (R = –0.77, P < 0.001). (B) Effect of changes in Cd concentration in the sediment (Cds) during the investigation period on its concentration in E. danica (CdE, black circles) and BSAF (CdB, open triangles) at the investigated localities on the Skrapež River. The black solid line represents the linear fit of the data in the case of CdE (R = –0.54, P = 0.030) and the grey one in the case of CdB (R = –0.88, P < 0.001). (C) Effect of changes in Cu concentration in sediment (Cus) during the investigation period on its concentration in E. danica (CuE, black circles) and BSAF (CuB, open triangles) at the investigated localities on the Skrapež River. The black solid line represents the linear fit of the data in the case of CuE (R = 0.72, P = 0.002) and the grey one in the case of CdB (R = 0.53, P = 0.036).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.