Next Article in Journal
Using Single-Species and Whole Community Stream Mesocosm Exposures for Identifying Major Ion Effects in Doses Mimicking Resource Extraction Wastewaters
Previous Article in Journal
Water Quality Modelling for Nitrate Nitrogen Control Using HEC-RAS: Case Study of Nakdong River in South Korea
Previous Article in Special Issue
Active Biomonitoring of Heavy Metal Concentrations in Aquatic Environment Using Mosses and Algae
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Environmental Variables on Biochemical Biomarkers in the Amphipod Monoporeia affinis from the Gulf of Riga (Baltic Sea)

1
Latvian Institute of Aquatic Ecology, Agency of Daugavpils University, Voleru Str. 4, LV-1007 Riga, Latvia
2
Department of Marine Systems, Tallinn University of Technology, Akadeemia Tee 15a, 12618 Tallinn, Estonia
3
Marine Research Centre, Finnish Environment Institute SYKE, Latokartanonkaari 11, FI-00790 Helsinki, Finland
*
Author to whom correspondence should be addressed.
Water 2023, 15(2), 248; https://doi.org/10.3390/w15020248
Submission received: 14 December 2022 / Revised: 29 December 2022 / Accepted: 4 January 2023 / Published: 6 January 2023
(This article belongs to the Special Issue Biological Monitoring of Pollution in Water Environment)

Abstract

:
The complexity of the marine environment and the increasing anthropogenic pressure create a necessity to expand existing monitoring approaches. The main goal of this study was to depict the effects of selected, seasonally varying environmental factors on a battery of biomarkers in the benthic amphipod Monoporeia affinis from the Gulf of Riga (GoR). Seasonal variability in acetylcholinesterase (AChE), catalase (CAT), glutathione reductase (GR), and glutathione S-transferase (GST) activities was investigated at six coastal stations (20–30 m) in August and November in 2020 and 2021. In addition, the biomarkers were measured at seven deep-water stations (>30 m) in November 2021. In general, the results indicated no significant influence of the measured environmental variables on the biomarker activities, except for deep-water stations, where chlorophyll a significantly affected enzymatic activity. The current study indicated that M. affinis has a higher GST, CAT and GR activity in summer compared to autumn in coastal stations, showing seasonal variability of these biomarkers. However, summarizing the biomarker levels recorded at each station and season, the integrated biomarker response (IBR) index showed the most stressed health status of the M. affinis populations in the deep-water stations 135 and 107 and coastal regions in the north-eastern part of the GoR (station 160B). This suggests that the impact on enzymatic responses of benthic organisms could be due to port activities leading to the accumulation of pollutants in muddy sediments regionally. Moreover, for the monitoring of biological effects of contaminants there is a need to establish the background levels of biomarkers, i.e., responses to the different natural environmental factors in the GoR region.

1. Introduction

Environmental quality control of aquatic ecosystems receives more and more attention due to increasing economic development and rising pollution. Anthropogenic activities are the main factors in increasing the level of contaminants in aquatic environments [1]. As a result, long-term exposure of aquatic organisms to contaminants leads to changes in physiological responses, accumulation of pollutants in tissues, changes in population level, and, finally, changes in species diversity [2]. Applying molecular, biochemical and physiological biomarkers in environmental monitoring programs helps to detect and identify the sub-lethal early-warning effects of pollution on organisms concerning their health, fitness, growth, and reproductive capacity related to ecosystem health in general [3].
Biochemical biomarkers offer information regarding the potential impact of toxic pollutants on the health of organisms and they represent different types of biological responses to different stressors [4]. A large number of chemical contaminants undergo oxidative reactions in cells, leading to the excess production of reactive oxygen species (ROS), which may cause a common phenomenon called oxidative stress [5]. Antioxidant enzyme activity is an important part of the antioxidant defense system (ADS) targeted to prevent oxidative damage to cellular components. In ecotoxicology, the most commonly used enzymes to detect oxidative stress include catalase (CAT), glutathione reductase (GR) as well as the phase II biotransformation enzyme glutathione S-transferase (GST) [6,7]. The inhibition of acetylcholinesterase activity (AChE), a key enzyme in the cholinergic nervous system, is a classic biomarker of exposure to neurotoxic compounds, particularly to organophosphates and carbamate pesticides [8], but has been recorded to occur also in connection to exposure to other pollutants such as trace metals, detergents and cyanobacterial toxins [9,10]. These common biomarkers are often used to detect exposure to pollution in aquatic organisms in laboratory and field studies [11].
The selection of suitable organisms for a biomonitoring program is an essential step. Amphipods are considered to be excellent bioindicator organisms as they are widespread over large salinity and habitat ranges [12] and respond to various types of environmental contaminants [6]. They are often an important food source for many fish and invertebrate species, and are thus considered highly relevant organisms for ecotoxicological studies [13]. The amphipod Monoporeia affinis (Lindström 1855) is an ecological keystone species of the soft-bottom macrozoobenthic communities in the Baltic Sea and thus a relevant bioindicator organism to be used to monitor and assess anthropogenic impacts [14,15,16,17]. Malformed embryos in amphipods are currently recommended as a supplementary indicator of contamination effects on aquatic organisms in the Baltic Sea and are regularly monitored by some countries [13,18,19].
For the correct interpretation of biological responses to chemical contamination, it is necessary to recognize the variability in the applied parameters under different natural environmental conditions such as temperature, dissolved oxygen level, salinity, photoperiod and food availability [20,21,22]. Variability in some responses is also a natural feature of the annual physiological cycle of the species and can be caused by intrinsic confounding factors such as reproductive status [21,23,24]. Regarding M. affinis, the combined effects of pollution in the sediment and oxygen deficiency on ADS responses have been recorded using biochemical biomarkers and reproductive disorders both in laboratory experiments [14] and in field studies [18]. In the Gulf of Riga (GoR), the Baltic Sea subregion of the current study, previous studies regarding the effects of pollution on the growth, reproduction and survival of amphipods have been carried out under standardized laboratory conditions [25,26,27,28]. However, biochemical biomarker studies from field-collected samples in the area have been carried out only on the soft-bottom clam Macoma balthica (Linnaeus, 1758) and the mussel Mytilus spp. [10,20,22,26]. Information on the biological effects of contaminants on amphipods from field studies in the Baltic Sea is relatively scarce and, partly therefore, the application of biological effects methods in monitoring in the region has been underdeveloped until very recently. The research on biochemical biomarkers in M. affinis would be the first information in the GoR region and could therefore be used as the basis of regular biomonitoring of contamination effects in the area for future.
The aim of this study was to assess seasonal and spatial variability in selected biochemical biomarkers in M. affinis in the GoR as a possible starting point for regular monitoring activities in the area. The above-mentioned ADS response biomarkers CAT, GR and GST were selected to detect possible contaminant-induced oxidative stress while AChE inhibition represents exposure to directly neurotoxic substances or indirect inhibitory effects caused by other types of compounds. The samplings were carried out in 2020 and 2021 in summer and autumn in both years, using a network of sampling sites consisting of shallow coastal and deep-water areas of the GoR. A common set of environmental variables were recorded simultaneously to examine the relationships of these factors with the variability observed in the biomarker levels in the study area.

2. Materials and Methods

2.1. Study Area

The Gulf of Riga is located in the northeastern part of the Baltic Sea and is a semi-enclosed basin with a relatively low salinity of 0.5–7.7 ppt due to its isolation from the Baltic Proper and significant freshwater discharges [29,30]. The GoR surface area is 16.330 km2 (3.9% of the Baltic Sea area), water volume 424 km3 (2.1% of the Baltic Sea volume), average depth 26 m, and maximum depth 62 m [31]. Near-bottom water temperature remains low all year round (4–6 °C) [32,33] and dissolved oxygen concentrations range from 2 to 6 mg/L [31,34]. At the depths of 20–30 m, soft bottom substrates or sandy sediments and boulders with well-developed benthic communities prevail with richness up to 13 species, while silty and muddy sediments mostly dominate in the deeper parts (>30 m) of the gulf consisting of only a few amphipod and polychaete species [27,32,35,36].

2.2. Sampling

The M. affinis sampling stations were chosen based on the species abundance data available from the Latvian national marine monitoring. The amphipods were collected aboard the R/V “Salme” in August and November in 2020 and 2021 with a Van Veen grab (0.1 m2) at the depths of 20–43 m (Figure 1) and sieved through a 0.5 mm mesh. The amphipods were immediately frozen in liquid nitrogen, and later placed in a −80 °C freezer for the later biomarker analyses. Salinity, temperature, chlorophyll a, and dissolved oxygen concentration, in the near-bottom layer were measured at each site with a multiparameter water quality probe CTD profiler SBE 19 plus SeaCAT (Bellevue, WA, USA).

2.3. Sample Preparation and Biomarker Analyses

Whole bodies of five amphipods were pooled for each of the 12 replicate biomarker samples per station. Six replicates were homogenized for 2 × 45 s (Retsch MM400 homogenizer, Haan, Germany) in cold (4 °C) 0.1 M phosphate buffer with 0.15 M KCl (pH 7.4) (1:4 w/v ratio) to measure the ADS enzyme activities. The remaining six replicates were used to determine AChE with homogenization using cold 0.02 M phosphate buffer with 0.1% Triton X-100 (pH 7.0). During homogenization, the vials were kept on ice. The homogenate was centrifuged at 12,000× g at 4 °C for 20 min and the supernatant was used for the assays. Four measurements were carried out per replicate. The enzyme activities were measured with a microplate reader (Spark® Multimode Microplate Reader, TECAN, Grödig, Austria) and analyzed using Magellan software (software version 2.2). The reaction rate was evaluated according to the best linear range of the obtained curve. The specific activities of all the enzymatic biomarkers were calculated against the total protein content of the sample.
The CAT activity was determined according to Claiborne [37]. The supernatant was diluted with the phosphate buffer in a 1:10 (v/v) ratio. The CAT activity was measured recording the decrease of 30 mM H2O2 at 240 nm. The measurement of GST activity is based on the conjugation of reduced glutathione (GSH) to 1-chloro-2,4-dinitrobenzene (CDNB) using a modification of the method based on Habig et al. [38]. The supernatant in a reaction mix (Dulbecco’s buffer PBS, 0.1 M GSH and 0.1 M CDNB) were used to measure GST activity at 340 nm. The GR activity level was evaluated as the oxidation rate of NADPH. The supernatant in a reaction mix (0.1 M phosphate buffer + 2 mM EDTA [pH 7.5], 2 mM GSSG, 2 mM NADPH, 3 mM DTNB) was analyzed for GR activity at 412 nm [39]. For the determination of AChE activity according to Ellman et al. [40], the supernatant in 0.02 M phosphate buffer (pH 7.0), 0.1 M acetylcholine iodide (ACTC) and 0.01 M 5,5′-dithiobis (2-nitrobenzoic acid) (DTNB) was measured at a 412 nm.
The total amount of protein in the homogenate of each sample was determined with the Bradford assay using bovine serum albumin as the standard [41].

2.4. Statistical Analysis

The biomarker levels were expressed as the mean and standard error (mean ± SE) recorded for amphipods collected from each station by each month, year, and the coastal/deep-water grouping. To ensure homogeneity of data while assessing the seasonal and interannual variability based on the grouped data of the biomarker activity of all coastal stations together, only stations 101A, 167B, VAD2, 163B, 162B were considered in the calculations. The environmental variables recorded were expressed as the mean and standard deviation (mean ± SD). The normality of the biomarker data distribution was checked via the Shapiro–Wilk test. As most of the data were not normally distributed, non-parametric tests were applied. Due to failed data normality, the non-parametric Kruskal–Wallis one-way analysis of variance in combination with the Wilcoxon rank sum test was performed to investigate seasonal and annual variability in medians of biomarker activities across the entire GoR and at the individual sampling stations. Correlations between biomarker levels and physicochemical variables were examined by means of the Spearman’s rank correlation and visualized as a correlation matrix using the “corrplot” package of the R software. Data exploration, artworks, and statistical analyses were performed using the R software for Windows, release 4.0.3.
All the measured biomarkers were combined into one general “stress index” known as the Integrated Biomarker Response (IBR) index [42]. The procedure used for the IBR calculation was based on the original paper by Beliaeff and Burgeot [42], modified by Broeg and Lehtonen [43]. The AChE, GST, CAT, and GR data from the two seasons (August and November) in the year 2020 were used for the IBR calculation for the coastal stations while the November data from the year 2021 were used to compare IBR at the coastal and deep-water stations.

3. Results

3.1. Environmental Factors

The measured physicochemical parameters (Table 1) varied moderately between the two seasons (August and November) and the station depth group. The data presented in Table 1 shows that both compared years, 2020 and 2021, differed significantly in terms of water temperature: in August the median in the group of coastal water sites (20–30 m depth) was respectively 8.1 °C and 3.3 °C, while in November it was 10.8 °C and 8.4 °C, respectively. At the coastal stations, higher average values in temperature, chlorophyll a and oxygen level were observed in November (9.6 ± 0.4 °C, 2.6 ± 0.2 mg/m3, and 10.4 ± 0.2 mg/L, respectively) compared to those measured in August (7.5 ± 1.4 °C, 1.3 ± 0.3 mg/m3 and 5.5 ± 0.5 mg/L, respectively). In addition, in 2021, higher concentrations of dissolved oxygen and chlorophyll a in water were recorded in all coastal stations compared to the year 2020. A similar trend was recorded between the depth groups in November, with higher values in temperature, chlorophyll a and oxygen level being detected at the coastal stations compared to the deep-water stations (8.1 ± 0.5 °C, 2.2 ± 0.2 mg/m3, and 9.1 ± 1.2 mg/L, respectively). The average water temperature was slightly lower (by 0.4 °C) in deep-water sites in comparison to coastal sites. At the same time, the median in both groups in (November 2021) of stations (coastal and deep-water) was identical (8.4 °C). However, in November 2021 the higher concentrations of chlorophyll a and oxygen content were recorded at coastal stations.

3.2. Seasonal Variability in Biomarkers at the Coastal Stations

In M. affinis, the median AChE activity of the pooled data of coastal stations showed no significant variability between August and November (p > 0.05) in both year (Figure 2). Although, looking at each station separately, the AChE activity level at station 160B was lower in August than in November in 2020 (p = 0.020) and at station 163B in 2021 (p = 0.01). Interannual variability could be seen both in August (227.5 ± 7.6 and 191.9 ± 8.3 nmol/min/mg/protein, in 2020 and 2021, respectively, p = 0.005) and November (242.1 ± 11.2 and 179.3 ± 7.4 nmol/min/mg/protein, p < 0.001).
Grouped coastal stations CAT activity did not show any significant variability between August and November in 2020, however it was significant (p = 0.017) in 2021 (Figure 2). In 2020, a higher activity was recorded in August at some individual stations, namely 101A (p = 0.009), 160B (p = 0.043), and 167B (p = 0.008). No significant differences between the years could be found in the August samples but in November 2020 the activity was significantly lower (p < 0.001) compared to November 2021 (126.9 ± 3.9 and 154.4 ± 3.6 µmol/min/mg/protein, respectively).
In GST, significant interannual (2020 versus 2021) variability has been observed for both August (p = 0.019) and November (p = 0.003) (Figure 2). Seasonal variability was significant (p < 0.001) both for 2020 (133.9 ± 5.1 and 105.1 ± 3.2 nmol/min/mg/protein in August and November, respectively) and 2021 (148.8 ± 3.9 and 123.8 ± 2.4 nmol/min/mg/protein). Seasonal differences in GST activity were significant at most of the coastal stations including 160B (p = 0.021), 162B (p = 0.014), and 163B (p = 0.027) in 2020, and 101A (p = 0.004), 162B (p = 0.010), 167B (p = 0.004) in 2021.
Monoporeia affinis at the individual stations 101A and 167B in 2020 (p = 0.028 and 0.037, respectively) and 101A and 163B in 2021 (p = 0.024 and 0.016) showed higher GR activity in August (Figure 2). Interannual variability was significant in the November samples where the GR activity across the stations was higher in 2021 (3.9 ± 0.2 nmol/min/mg/protein) compared to 2020 (3.3 ± 0.1 nmol/min/mg/protein) (p < 0.001).

3.3. Biochemical Biomarkers at the Deep-Water Stations

Biomarker levels in M. affinis at the deep-water stations were measured only in November 2021 (Figure 3). The highest values of AChE activity were observed at station 135 (262.2 ± 24.6 nmol/min/mg/protein), being significantly higher (p < 0.05) than at stations 114A, 137A and 111 (115.0 ± 7.8, 140.3 ± 7.1 and 163.9 ± 12.4 nmol/min/mg/protein, respectively).
For CAT, the activity levels were higher at stations 107, 135 and 142 (from 201.3 ± 7.1 to 214.7 ± 17.2 µmol/min/mg/protein) compared to the other four stations (p < 0.01).
The GST activity was significantly higher (p < 0.05) in samples collected from stations 107 and 135 (142.7 ± 6.9 and 152.1 ± 13.7 nmol/min/mg/protein, respectively), compared to stations 102A and 111 (105.7 ± 4.30 and 105.2 ± 5.9 nmol/min/mg/protein) (Figure 3).
The GR activity was significantly lower (p = 0.011 to 0.031) in samples collected from station 111 (3.0 ± 0.2 nmol/min/mg/protein) compared to the other stations (5.2 ± 0.4, 5.7 ± 0.4, 4.6 ± 0.4, 6.7 ± 0.6, 5.1 ± 0.3 nmol/min/mg/protein for stations 102A, 107, 114A, 135 and 142, respectively), except 137A (4.1 ± 0.5 nmol/min/mg/protein).
Moreover, the activity of CAT and GR among the deep-water stations was significantly higher (p < 0.001) than it was among the coastal stations in November 2021.

3.4. Correlation Analysis

Only some of the physicochemical variables showed strong or moderate correlations within a season at the coastal stations (Figure 4). In general, more correlations were found in November compared to August. Near-bottom oxygen concentrations showed a negative correlation with AChE activity of the amphipods both in August and November while being positively correlated with CAT and GR activity in November. Moreover, chlorophyll a concentration exhibited a strong negative correlation with AChE and was positive with CAT in November, while no significant relationships between the variables could be observed in August. Regarding temperature, only AChE activity in November showed a significant positive correlation. In the deeper areas of the gulf in November, chlorophyll a concentration showed a significant positive correlation with AChE, CAT, and GST (Figure 4).

3.5. Integrated Biomarker Response

The IBR calculations showed variability in the integrated response at the different stations, seasons, and depth ranges (Figure 5). Across the coastal stations, higher index values were detected in August at stations 160B and 167B (Figure 5a). In November, the most impacted amphipod population was observed at the same station 160B while the least impacted was 162B in both seasons. In the comparison of coastal and deep-water station data in November, higher IBR index values were commonly observed at the deep-water stations, especially at 135, 107 and 142, with the notable exception of station 111 (Figure 5b).

4. Discussion

This study focused on the spatial and temporal differences in selected biomarker responses in a key species in the Gulf of Riga study area and the associations of the responses with some environmental factors. The results obtained provide the first information on biochemical biomarkers in M. affinis in this sea area and can therefore be used as the basis of regular biomonitoring of contamination effects in the area.
Seasonal variability has been recognized as an important factor influencing the baseline levels of biomarkers and also the responsiveness of organisms to pollution stress [21,44,45], e.g., in previous studies on M. balthica from the GoR [25]. In the present study, seasonal variability in GST and CAT was detected in the amphipods at all the study stations. Of the effects of specific environmental variables measured, the impact of chlorophyll a on AChE, CAT, and GST was observed in the samples collected in November in the deeper parts of the GoR. At the same time, a significant negative correlation between chlorophyll a in the bottom layer and AChE activity of the amphipods was detected at the coastal stations in November. Normally, a high chlorophyll a level in water would indicate good feeding opportunities for filter and deposit feeders and result in an improved physiological condition; however, if associated with toxic phytoplankton blooms such as those regularly occurring in the Baltic Sea in late summer it can cause the inhibition of AChE due to residual cyanobacterial toxins [9,10,14]. Among other abiotic factors, Löf et al. [46] found that both salinity and temperature significantly modified some biomarker responses to contaminants in M. affinis. Negative effects of increased temperature during contaminant exposure on M. affinis, both on a biochemical and organism level, were recorded also in laboratory conditions [47]. However, in our study, temperature positively and oxygen negatively correlated with AChE activity in November at the coastal stations. A rise in temperature can increase the oxygen consumption of the amphipods, thus increasing energy cost and inhibiting quality of reproductivity [48].
Apart from abiotic variables, biotic factors such as reproductive status, body size and food availability can affect enzymatic responses in amphipods [14,49,50,51]. Therefore, the physiological status of the organisms should be taken into account in the interpretation of biomarker responses in field studies [52]. Monoporeia affinis has a long reproductive cycle; oogenesis starts in late summer and mating takes place in November [53], and the physiological processes during the breeding period lead to the mobilization of energy stores within the organism, potentially increasing their sensitivity to environmental stressors [52]. Significant seasonal differences were also detected in the activity of enzymatic biomarkers between summer and autumn in the clam M. balthica from the GoR [20] and the Gulf of Finland [22]. The current study indicated that M. affinis has a higher GST, CAT and GR activity in summer compared to autumn. For the amphipod Hyalella kaingang, Braghirolli et al. [21] reported a significantly higher CAT activity in summer compared to autumn; however, opposite results were obtained regarding GST activity. Conclusively, the observations above emphasize the importance of being aware of possible different seasonal baselines in biomarkers that can depend on various abiotic and biotic factors [54].
Various biological processes provide information on the types of stressors affecting benthic soft-bottom communities and the reproduction of amphipods possibly related to their observed population decreases in the Baltic Sea [17,55]. Parallel studies in the GoR (unpublished data) indicate that the reproductive success of M. affinis is moderately associated with the measured biomarker levels under natural conditions and this information can be utilized when designing early warning monitoring schemes using these parameters. Although the observations on the variability in biomarkers presented in this study are not temporally broad enough to draw conclusions on their long-term dynamics, the significant differences observed in all four biomarker responses between 2020 and 2021 highlight the necessity for longer-term monitoring to reveal the drivers of the observed population declines.
The central part of the GoR (i.e., the deep-water stations) is characterized by a low-diversity benthic community, muddy sediments and relatively high trace metal concentrations [27], where pollutants are relocated and accumulated by currents. The integrated stress response (IBR) in amphipods, assessed by combining the information obtained from all four biomarkers selected for this study, shows a peak in autumn and is more pronounced in the deep-water areas. The elevated biomarker IBR responses were generally observed in M. affinis collected from the central region of the GoR (deep-water stations 135, 142) and the north-eastern part of the GoR (stations 160B and 107) as well as at stations situated in the river estuaries. In many cases the biomarker levels were higher at the deep-water stations compared to the coastal area. Elevated trace metal concentrations, particularly of Cd, have been recorded in sediments at the deep-water areas as well as at one sublittoral station (160B) [27]. Obviously, there are also many other contaminants present in the GoR marine environment. Using trace metal levels in sediments as a rough proxy of anthropogenic chemical pollution and their apparent association with the biomarker responses, including the IBR, gives an indication of the usefulness of the biomarker approach when assessing the biological effects of environmental contamination in the area. However, as shown in this study, natural variability has also to be carefully considered when interpreting the results.

5. Conclusions

The current study represents the first comprehensive report on a battery of biochemical biomarkers measured in M. affinis collected from the GoR at different seasons and years. The results showed that biological effects methods can be used in environmental quality assessment in this Baltic Sea area. Higher response levels were recorded in populations inhabiting the central (deeper) and northeastern regions of the GoR. Furthermore, fewer differences in biomarker levels were observed in amphipods collected in August compared to November, indicating apparent seasonal variability. The observations highlight the importance of understanding the physiology of a species in its natural habitat, including the periods of greater vulnerability to environmental stressors due to biotic and abiotic factors. Using a specifically designed battery of biomarkers together with measuring a variety of contaminants and key physicochemical variables in monitoring improves our understanding of biochemical responses to environmental stress. This research shows that the use of amphipods for integrated chemical-biological monitoring as well as for the evaluation of sediment quality of the GoR is possible. Finally, the results contribute to the knowledge on the effects of environmental factors on biomarker responses in amphipods and to clarify their dependency on the season that will lead to the development of their natural baseline/threshold values needed in biomonitoring.

Author Contributions

E.S.: conceptualization, data curation, resources, project administration, funding acquisition, formal analysis, validation, visualization, writing—original draft, writing—review and editing; I.B.: conceptualization, resources, validation, visualization, writing—original draft, writing—review and editing; N.S.: data curation, software, visualization, writing—original draft, writing—review and editing. N.K.: validation, writing—review and editing; R.T.: methodology, writing—review and editing; K.K.L.: investigation, methodology, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Regional Development Fund, 1.1.1.2/16/l/001 Post-doctoral project No. 1.1.1.2/VIAA/3/19/465.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.

Acknowledgments

Our special thanks to the Latvian Institute of Aquatic Ecology laboratory assistants and researchers who participated in the sediment and amphipods sampling, as well as colleagues who made physicochemical analyses.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bashir, I.; Lone, F.A.; Bhat, R.A.; Mir, S.A.; Dar, Z.A.; Dar, S.A. Concerns and Threats of Contamination on Aquatic Ecosystems. In Bioremediation and Biotechnology; Springer: Cham, Switzerland, 2020; pp. 1–26. [Google Scholar] [CrossRef] [Green Version]
  2. Hook, S.E.; Gallagher, E.P.; Batley, G.E. The Role of Biomarkers in the Assessment of Aquatic Ecosystem Health. Integr. Environ. Assess. Manag. 2014, 10, 327–341. [Google Scholar] [CrossRef] [Green Version]
  3. Smit, M.G.D.; Bechmann, R.K.; Hendriks, A.J.; Skadsheim, A.; Larsen, B.K.; Baussant, T.; Bamber, S.; Sanni, S. Relating biomarkers to whole-organism effects using species sensitivity distributions: A pilot study for marine species exposed to oil. Environ. Toxicol. Chem. 2009, 28, 1104–1109. [Google Scholar] [CrossRef]
  4. Lionetto, M.G.; Caricato, R.; Giordano, M.E. Pollution Biomarkers in Environmental and Human Biomonitoring. Open Biomark. J. 2019, 9, 1–9. [Google Scholar] [CrossRef] [Green Version]
  5. Stohs, S.J.; Bagchi, D. Oxidative mechanisms in the toxicity of metal ions. Free Radic. Biol. Med. 1995, 18, 321–336. [Google Scholar] [CrossRef] [Green Version]
  6. Turja, R.; Sanni, S.; Stankevičiūtė, M.; Andreikėnaitė Butrimavičienė, L.; Devier, M.-H.; Budzinski, H.; Lehtonen, K. Biomarker responses and accumulation of polycyclic aromatic hydrocarbons in Mytilus trossulus and Gammarus oceanicus during exposure to crude oil. Environ. Sci. Pollut. Res. 2020, 27, 15498–15514. [Google Scholar] [CrossRef] [Green Version]
  7. Vranković, J.; Živić, M.; Radojević, A.; Perić-Mataruga, V.; Todorović, D.; Marković, Z.; Živić, I. Evaluation of oxidative stress biomarkers in the freshwater gammarid Gammarus dulensis exposed to trout farm outputs. Ecotoxicol. Environ. Saf. 2018, 163, 84–95. [Google Scholar] [CrossRef]
  8. Umar, A.M.; Aisami, A. Acetylcholinesterase Enzyme (AChE) as a Biosensor and Biomarker for Pesticides: A Mini Review. Bull. Environ. Sci. Sustain. Manag. 2020, 4, 7–12. [Google Scholar] [CrossRef]
  9. Kankaanpää, H.; Leiniö, S.; Olin, M.; Sjövall, O.; Meriluoto, J.; Lehtonen, K.K. Accumulation and depuration of cyanobacterial toxin nodularin and biomarker responses in the mussel Mytilus edulis. Chemosphere 2007, 68, 1210–1217. [Google Scholar] [CrossRef]
  10. Lehtonen, K.K.; Kankaanpää, H.; Leiniö, S.; Sipiä, V.O.; Pflugmacher, S.; Sandberg-Kilpi, E. Accumulation of nodularin-like compounds from the cyanobacterium Nodularia spumigena and changes in acetylcholinesterase activity in the clam Macoma balthica during short-term laboratory exposure. Aquat. Toxicol. 2003, 64, 461–476. [Google Scholar] [CrossRef]
  11. Lionetto, M.G.; Caricato, R.; Giordano, M.E. Pollution Biomarkers in the Framework of Marine Biodiversity Conservation: State of Art and Perspectives. Water 2021, 13, 1847. [Google Scholar] [CrossRef]
  12. Whiteley, N.M.; Rastrick, S.P.S.; Lunt, D.H.; Rock, J. Latitudinal variations in the physiology of marine gammarid amphipods. J. Exp. Mar. Biol. Ecol. 2011, 400, 70–77. [Google Scholar] [CrossRef]
  13. Podlesińska, W.; Dąbrowska, H. Amphipods in estuarine and marine quality assessment—A review. Oceanologia 2019, 61, 179–196. [Google Scholar] [CrossRef]
  14. Gorokhova, E.; Löf, M.; Reutgard, M.; Lindström, M.; Sundelin, B. Exposure to contaminants exacerbates oxidative stress in amphipod Monoporeia affinis subjected to fluctuating hypoxia. Aquat. Toxicol. 2013, 127, 46–53. [Google Scholar] [CrossRef]
  15. Gorokhova, E.; Löf, M.; Halldórsson, H.P.; Tjärnlund, U.; Lindström, M.; Elfwing, T.; Sundelin, B. Single and combined effects of hypoxia and contaminated sediments on the amphipod Monoporeia affinis in laboratory toxicity bioassays based on multiple biomarkers. Aquat. Toxicol. 2010, 99, 263–274. [Google Scholar] [CrossRef]
  16. Lehtonen, K.K. Seasonal variations in the physiological condition of the benthic amphipods Monoporeia affinis and Pontoporeia femorata in the Gulf of Riga (Baltic Sea). Aquat. Ecol. 2004, 38, 441–456. [Google Scholar] [CrossRef]
  17. Sundelin, B.; Eriksson Wiklund, A.-K. Malformations in embryos of the deposit-feeding amphipod Monoporeia affinis in the Baltic Sea. Mar. Ecol. Prog. Ser. 1998, 171, 165–180. [Google Scholar] [CrossRef]
  18. Löf, M.; Sundelin, B.; Bandh, C.; Gorokhova, E. Embryo aberrations in the amphipod Monoporeia affinis as indicators of toxic pollutants in sediments: A field evaluation. Ecol. Indic. 2016, 60, 18–30. [Google Scholar] [CrossRef]
  19. HELCOM. Reproductive Disorders: Malformed Embryos of Amphipods. HELCOM Supplementaryindicator Report. 2018, pp. 1–23, ISSN 2343-2543. Available online: https://helcom.fi/wp-content/uploads/2019/08/Reproductive-disorders-malformed-embryos-of-amphipods-HELCOM-supplementary-indicator-2018.pdf (accessed on 3 January 2023). [CrossRef]
  20. Barda, I.; Purina, I.; Rimsa, E.; Balode, M. Seasonal dynamics of biomarkers in infaunal clam Macoma balthica from the Gulf of Riga (Baltic Sea). J. Mar. Syst. 2014, 129, 150–156. [Google Scholar] [CrossRef]
  21. Braghirolli, F.M.; Oliveira, M.R.; Oliveira, G.T. Seasonal variability of metabolic markers and oxidative balance in freshwater amphipod Hyalella kaingang (Crustacea, Amphipoda). Ecotoxicol. Environ. Saf. 2016, 130, 177–184. [Google Scholar] [CrossRef]
  22. Leiniö, S.; Lehtonen, K.K. Seasonal variability in biomarkers in the bivalves Mytilus edulis and Macoma balthica from the northern Baltic Sea. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 2005, 140, 408–421. [Google Scholar] [CrossRef]
  23. Benito, D.; Ahvo, A.; Nuutinen, J.; Bilbao, D.; Saenz, J.; Etxebarria, N.; Lekube, X.; Izagirre, U.; Lehtonen, K.K.; Marigómez, I.; et al. Influence of season-depending ecological variables on biomarker baseline levels in mussels (Mytilus trossulus) from two Baltic Sea subregions. Sci. Total Environ. 2019, 689, 1087–1103. [Google Scholar] [CrossRef]
  24. Sheehan, D.; Power, A. Effects of seasonality on xenobiotic and antioxidant defence mechanisms of bivalve molluscs. Comp. Biochem. Physiol. Part C Pharmacol. Toxicol. Endocrinol. 1999, 123, 193–199. [Google Scholar] [CrossRef]
  25. Berezina, N.A.; Strode, E.; Lehtonen, K.K.; Balode, M.; Golubkov, S.M. Sediment quality assessment using Gmelinoides fasciatus and Monoporeia affinis (Amphipoda, Gammaridea) in the northeastern Baltic Sea. Crustaceana 2013, 86, 780–801. [Google Scholar] [CrossRef]
  26. Putna, I.; Strode, E.; Barda, I.; Purina, I.; Rimša, E.; Jansons, M.; Balode, M.; Strake, S. Sediment quality of the ecoregion Engure, Gulf of Riga, assessed by using ecotoxicity tests and biomarker responses. Proc. Latv. Acad. Sciences. Sect. B. Nat. Exact Appl. Sci. 2014, 68, 20–30. [Google Scholar] [CrossRef] [Green Version]
  27. Strode, E.; Jansons, M.; Purina, I.; Balode, M.; Berezina, N.A. Sediment quality assessment using survival and embryo malformation tests in amphipod crustaceans: The Gulf of Riga, Baltic Sea AS case study. J. Mar. Syst. 2017, 172, 93–103. [Google Scholar] [CrossRef]
  28. Strode, E.; Balode, M. Toxico-resistance of Baltic amphipod species to heavy metals. Crustaceana 2013, 86, 1007–1024. [Google Scholar] [CrossRef]
  29. Ojaveer, E. Large-scale processes in the ecosystem of the Gulf of Riga. In Ecosystem of the Gulf of Riga between 1920 and 1990; Ojaveer, E., Ed.; Estonian Academy Publishers: Tallinn, Estonia, 1995; pp. 268–277. [Google Scholar]
  30. Skudra, M.; Lips, U. Characteristics and inter-annual changes in temperature, salinity and density distribution in the Gulf of Riga. Oceanologia 2017, 59, 37–48. [Google Scholar] [CrossRef] [Green Version]
  31. Berzinsh, V. Dynamics of hydrological parameters of the Gulf of Riga. In Ecosystem of the Gulf of Riga between 1920 and 1990; Ojaveer, E., Ed.; Estonian Academy: Tallinn, Estonia, 1995; pp. 8–32. [Google Scholar]
  32. Raudsepp, U. Interannual and Seasonal Temperature and Salinity Variations in the Gulf of Riga and Corresponding Saline Water Inflow From the Baltic Proper. Hydrol. Res. 2001, 32, 135–160. [Google Scholar] [CrossRef]
  33. Yurkovskis, A. Long-term land-based and internal forcing of the nutrient state of the Gulf of Riga (Baltic Sea). J. Mar. Syst. 2004, 50, 181–197. [Google Scholar] [CrossRef]
  34. Eglīte, E.; Lavrinovičs, A.; Müller-Karulis, B.; Aigars, J.; Poikāne, R. Nutrient turnover at the hypoxic boundary: Flux measurements and model representation for the bottom water environment of the Gulf of Riga, Baltic Sea. Oceanologia 2014, 56, 711–735. [Google Scholar] [CrossRef] [Green Version]
  35. Kotta, J.; Lauringson, V.; Martin, G.; Simm, M.; Kotta, I.; Herkül, K.; Ojaveer, H. Gulf of Riga and Pärnu Bay. In Ecology of Baltic Coastal Waters; Schiewer, U., Ed.; Ecological Studies; Springer: Berlin/Heidelberg, Germany, 2008; pp. 217–243. ISBN 978-3-540-73523-6. [Google Scholar]
  36. Viška, M.; Soomere, T. Simulated and observed reversals of wave-driven alongshore sediment transport at the eastern Baltic Sea coast. Baltica 2013, 26, 145–156. [Google Scholar] [CrossRef] [Green Version]
  37. Claiborne, A. Catalase Activity. In CRC Handbook of Methods for Oxygen Radical Research; CRC Press: Boca Raton, FL, USA, 1985; ISBN 978-1-351-07292-2. [Google Scholar]
  38. Habig, W.H.; Pabst, M.J.; Jakoby, W.B. Glutathione S-transferases. The first enzymatic step in mercapturic acid formation. J. Biol. Chem. 1974, 249, 7130–7139. [Google Scholar] [CrossRef]
  39. Smith, I.K.; Vierheller, T.L.; Thorne, C.A. Assay of glutathione reductase in crude tissue homogenates using 5,5’-dithiobis(2-nitrobenzoic acid). Anal. Biochem. 1988, 175, 408–413. [Google Scholar] [CrossRef]
  40. Ellman, G.L.; Courtney, K.D.; Andres, V.; Featherstone, R.M. A new and rapid colorimetric determination of acetylcholinesterase activity. Biochem. Pharmacol. 1961, 7, 88–95. [Google Scholar] [CrossRef]
  41. Bradford, M.M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 1976, 72, 248–254. [Google Scholar] [CrossRef] [PubMed]
  42. Beliaeff, B.; Burgeot, T. Integrated biomarker response: A useful tool for ecological risk assessment. Environ. Toxicol. Chem. 2002, 21, 1316–1322. [Google Scholar] [CrossRef] [PubMed]
  43. Broeg, K.; Lehtonen, K.K. Indices for the assessment of environmental pollution of the Baltic Sea coasts: Integrated assessment of a multi-biomarker approach. Mar. Pollut. Bull. 2006, 53, 508–522. [Google Scholar] [CrossRef] [Green Version]
  44. Jemec, A.; Drobne, D.; Tisler, T.; Sepcic, K. Biochemical biomarkers in environmental studies-lessons learnt from enzymes catalase, glutathione S-transferase and cholinesterase in two crustacean species. Environ. Sci. Pollut. Res. 2010, 17, 571–581. [Google Scholar] [CrossRef]
  45. Verlecar, X.N.; Jena, K.B.; Chainy, G.B.N. Seasonal variation of oxidative biomarkers in gills and digestive gland of green-lipped mussel Perna viridis from Arabian Sea. Estuar. Coast. Shelf Sci. 2008, 76, 745–752. [Google Scholar] [CrossRef]
  46. Löf, M.; Sundelin, B.; Liewenborg, B.; Bandh, C.; Broeg, K.; Schatz, S.; Gorokhova, E. Biomarker-enhanced assessment of reproductive disorders in Monoporeia affinis exposed to contaminated sediment in the Baltic Sea. Ecol. Indic. 2016, 63, 187–195. [Google Scholar] [CrossRef]
  47. Jacobson, T.; Prevodnik, A.; Sundelin, B. Combined effects of temperature and a pesticide on the Baltic amphipod Monoporeia affinis. Aquat. Biol. 2008, 1, 269–276. [Google Scholar] [CrossRef] [Green Version]
  48. Eriksson Wiklund, A.K.; Sundelin, B. Impaired reproduction in the amphipods Monoporeia affinis and Pontoporeia femorata as a result of moderate hypoxia and increased temperature. Mar. Ecol. Prog. Ser. 2001, 222, 131–141. [Google Scholar] [CrossRef] [Green Version]
  49. Chainy, G.B.N.; Paital, B.; Dandapat, J. An Overview of Seasonal Changes in Oxidative Stress and Antioxidant Defence Parameters in Some Invertebrate and Vertebrate Species. Scientifica 2016, 2016, 6126570. [Google Scholar] [CrossRef] [Green Version]
  50. Glippa, O.; Engström-Öst, J.; Kanerva, M.; Rein, A.; Vuori, K. Oxidative stress and antioxidant defense responses in Acartia copepods in relation to environmental factors. PLoS ONE 2018, 13, e0195981. [Google Scholar] [CrossRef] [Green Version]
  51. Taddei, A.; Räsänen, K.; Burdon, F.J. Size-dependent sensitivity of stream amphipods indicates population-level responses to chemical pollution. Freshw. Biol. 2021, 66, 765–784. [Google Scholar] [CrossRef]
  52. Sroda, S.; Cossu-Leguille, C. Seasonal variability of antioxidant biomarkers and energy reserves in the freshwater gammarid Gammarus roeseli. Chemosphere 2011, 83, 538–544. [Google Scholar] [CrossRef]
  53. Wiklund, A.-K.E.; Sundelin, B.; Rosa, R. Population decline of amphipod Monoporeia affinis in Northern Europe: Consequence of food shortage and competition? J. Exp. Mar. Biol. Ecol. 2008, 367, 81–90. [Google Scholar] [CrossRef]
  54. Luna Acosta, A.; Bustamante, P.; Godefroy, J.; Fruitier-Arnaudin, I.; Thomas-Guyon, H. Seasonal variation of pollution biomarkers to assess the impact on health status of juvenile Pacific oysters Crassostrea gigas exposed in situ. Environ. Sci. Pollut. Res. 2010, 17, 999–1008. [Google Scholar] [CrossRef]
  55. Rousi, H.; Laine, A.O.; Peltonen, H.; Kangas, P.; Andersin, A.-B.; Rissanen, J.; Sandberg-Kilpi, E.; Bonsdorff, E. Long-term changes in coastal zoobenthos in the northern Baltic Sea: The role of abiotic environmental factors. ICES J. Mar. Sci. 2013, 70, 440–451. [Google Scholar] [CrossRef]
Figure 1. Map of the study area in the Gulf of Riga with markings of coastal (20–30 m) and deep-water (>30 m) stations.
Figure 1. Map of the study area in the Gulf of Riga with markings of coastal (20–30 m) and deep-water (>30 m) stations.
Water 15 00248 g001
Figure 2. Activity of AChE, CAT, GST, and GR in M. affinis at the coastal stations (20–30 m) in the GoR in August and November 2020 and 2021. The red asterisk indicates statistically significant (p < 0.05) seasonal differences at the respective stations. Black dots denote outliers calculated according to the Interquartile range (IQR) criterion.
Figure 2. Activity of AChE, CAT, GST, and GR in M. affinis at the coastal stations (20–30 m) in the GoR in August and November 2020 and 2021. The red asterisk indicates statistically significant (p < 0.05) seasonal differences at the respective stations. Black dots denote outliers calculated according to the Interquartile range (IQR) criterion.
Water 15 00248 g002
Figure 3. Activity (mean ± SE) of AChE, CAT, GST, and GR in M. affinis at the deep-water stations in November 2021.
Figure 3. Activity (mean ± SE) of AChE, CAT, GST, and GR in M. affinis at the deep-water stations in November 2021.
Water 15 00248 g003
Figure 4. Spearman’s rank correlation coefficients between biomarker levels and physicochemical variables were calculated for six coastal stations in August and November, and for seven deep-water stations in November. Correlation coefficients (Spearman’s rho) with a p < 0.05 are shown and marked based on their sign (+ or −) and strength (continuous color scale).
Figure 4. Spearman’s rank correlation coefficients between biomarker levels and physicochemical variables were calculated for six coastal stations in August and November, and for seven deep-water stations in November. Correlation coefficients (Spearman’s rho) with a p < 0.05 are shown and marked based on their sign (+ or −) and strength (continuous color scale).
Water 15 00248 g004
Figure 5. Integrated biomarker response index (IBR) in the amphipod M. affinis collected from (a) coastal, year 2020 and (b) coastal and deep-water stations (marked *), year 2021 from the GoR.
Figure 5. Integrated biomarker response index (IBR) in the amphipod M. affinis collected from (a) coastal, year 2020 and (b) coastal and deep-water stations (marked *), year 2021 from the GoR.
Water 15 00248 g005
Table 1. Physicochemical parameters measured at the different study stations in the near-bottom layer in August and November 2020 and 2021.
Table 1. Physicochemical parameters measured at the different study stations in the near-bottom layer in August and November 2020 and 2021.
YearMonthStationDepth [m]Temp. [°C]Absolute Salinity [ppt]Chlorophyll a [mg/m3]Oxygen [mg/L]
2020August101A227.25.80.94.0
2021223.36.11.85.9
2020November2210.96.01.69.5
2021238.46.03.711.1
2020August167B216.45.80.84.7
2021213.26.21.66.4
2020November2110.85.92.09.7
2021228.46.03.011.1
2020August163B2212.15.90.96.4
20212218.45.92.26.9
2020November2210.85.92.19.7
2021218.55.93.011.1
2020August162B258.96.00.84.0
2021253.96.11.56.8
2020November2510.75.91.99.6
2021248.35.94.011.3
2020AugustVAD2265.55.90.94.6
2021263.06.21.58.3
2020November2610.35.92.39.8
2021257.76.13.011.1
2020August160B2210.15.90.92.9
2020November2210.45.72.49.9
2021November227.95.93.011.5
2021November107328.46.02.711.2
2021November111388.86.52.010.7
2021November114A338.66.81.710.9
2021November142428.26.32.310.7
2021November135458.26.12.65.3
2021November102A425.36.21.64.0
2021November137A429.06.02.210.7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Strode, E.; Barda, I.; Suhareva, N.; Kolesova, N.; Turja, R.; Lehtonen, K.K. Influence of Environmental Variables on Biochemical Biomarkers in the Amphipod Monoporeia affinis from the Gulf of Riga (Baltic Sea). Water 2023, 15, 248. https://doi.org/10.3390/w15020248

AMA Style

Strode E, Barda I, Suhareva N, Kolesova N, Turja R, Lehtonen KK. Influence of Environmental Variables on Biochemical Biomarkers in the Amphipod Monoporeia affinis from the Gulf of Riga (Baltic Sea). Water. 2023; 15(2):248. https://doi.org/10.3390/w15020248

Chicago/Turabian Style

Strode, Evita, Ieva Barda, Natalija Suhareva, Natalja Kolesova, Raisa Turja, and Kari K. Lehtonen. 2023. "Influence of Environmental Variables on Biochemical Biomarkers in the Amphipod Monoporeia affinis from the Gulf of Riga (Baltic Sea)" Water 15, no. 2: 248. https://doi.org/10.3390/w15020248

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop