BioRisk | A | 69-1 78 (2022) A Cras gr ei Ea rcs doi: 10.3897/biorisk.| 7.77384 RESEARCH ARTICLE & B lO R IS k https://biorisk.pensoft.net State of the antioxidant defense system in wedge clams from Bulgarian Black Sea as a measure of resistance to environmental impacts Almira Georgieva', Elina Tsvetanova', Nesho Chipev', Albena Alexandrova' | Laboratory of Free Radical Processes, Institute of Neurobiology, Bulgarian Academy of Sciences, Sofia, Bulgaria Corresponding author: Albena Alexandrova (aalexandrova@abv.bg) Academic editor:Roumiana Metcheva | Received 31 October 2021 | Accepted 29 November 2021 | Published 21 April 2022 Citation: Georgieva A, Tsvetanova E, Chipev N, Alexandrova A (2022) State of the antioxidant defense system in wedge clams from Bulgarian Black Sea as a measure of resistance to environmental impacts. In: Chankova S, Peneva V, Metcheva R, Beltcheva M, Vassilev K, Radeva G, Danova K (Eds) Current trends of ecology. BioRisk 17: 169-178. https://doi.org/10.3897/biorisk. 17.77384 Abstract Pollution and climate change can induce oxidative stress (OS) in aquatic organisms. Reduced activity or incoordination between antioxidant enzymes in marine bivalves may cause cellular impairment with ef- fects on higher levels of ecological organization. The present study aims to assess the condition factor and the activity of antioxidant enzymes catalase (CAT), superoxide dismutase (SOD), glutathione peroxidase (GPx), glutathione-S-transferase (GST), and glutathione concentration (GSH) in soft tissues of Donax trunculus (Linnaeus, 1758) as indicators of the state of marine ecosystems along the Bulgarian Black Sea coast. The wedge clams were sampled manually from different localities in June and in September. The activity of antioxidant enzymes and GSH concentrations varied both seasonally and among localities. Higher values of GSH, SOD and GPx were registered in wedge clams collected in autumn compared to those collected in summer. In wedge clams higher activity of the major phase II detoxification enzyme GST was observed in summer at localities with intensive tourism, suggesting an activation of the cell detoxification processes, presumably in response to increased environmental pressure. In conclusion, the observed presence of elevated enzyme activities suggested activation of the antioxidant protection system of the wedge clams in response to environmental pressure, indicating their ability to cope with induced OS and adapt to local conditions, and thus maintain ecosystem health. Keywords Antioxidant enzymes, Bulgarian Black Sea, condition factor, Donax trunculus, glutathione Copyright Almira Georgieva et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 170 Almira Georgieva et al. / BioRisk 17: 169-178 (2022) Introduction Pollution of aquatic ecosystems with heavy metals (Bat et al. 2018), polycyclic aro- matic hydrocarbons (PAHs) (Georgieva et al. 2016; Honda and Suzuki 2020), mi- croplastics (Andrady 2011), and pesticides (Mojiri et al. 2020) is becoming a serious environmental problem. Pollutants enter the marine environment by industrial and municipal wastewater, agricultural activities, tourist flow and atmospheric inputs (Maranho et al. 2015). As a result, marine biotas are subjected to an increasing and complex mixture of chemicals. Marine organisms respond to the impact of these pollutants by activating different cellular pathways and defense systems aiming to overcome the negative consequences, which may affect higher levels of organization. Thus, stress responses are considered an ecological driving force and trigger of evolu- tion (Steinberg 2012). Oxidative stress (OS) is the result of the disturbance of cellular pro/antioxidant balance and is a general reaction of aerobic organisms to endogenous and exogenous stimuli, including environmental contamination. OS biomarkers are considered as ear- ly-warning indicators of pollution, proposed to be used in biomonitoring programs. The activation of the antioxidant defense system is a factor determining the survival of the organism in conditions of increased environmental pressure. Lower antioxidant en- zymes activities or lower concentrations of non-enzymatic antioxidants suggest higher susceptibility to stressors, and higher possibility of tissue damages in marine organisms (Steinberg 2012). Comparative studies showed that bivalves are characterized by a highly active antioxidant enzyme complex and a much greater pool of non-enzymatic antioxidants than vertebrates (Soldatov et al. 2007). Wedge clams Donax trunculus (Linnaeus, 1758) are the most common species inhabiting the shallow sandy sea bottom of the Black Sea, where most of the pollutants are deposited. Being filter-feeders they are capable of bioaccumulating and concentrat- ing pollutants in their tissues. Furthermore, as a basic component of marine ecosys- tems, they are essential for maintaining ecosystem health. This study aimed to assess the state and changes of the antioxidant defense system of wedge clams D. trunculus from representative localities along the Bulgarian Black Sea coastal area as a measure of their resistance and adaptation to variable environ- mental pressures. Materials and methods Sampling Adult clams D. trunculus (length 23-35 mm) were sampled from their natural sandy habitats along the Bulgarian Black Sea coast (Fig. 1) in two seasons: summer (June) and autumn (September) of 2020. The samples were gathered manually or were obtained from commercial providers. The clams were stored at -80 °C until biochemical analyses. Antioxidant defense system in wedge clams from Bulgarian Black Sea 171 Condition Index (Cl) calculation Before the biochemical analyses, the weight and length of random clam individuals from the same sites were measured with an analytical scale (accuracy of 0.01 g) and a caliper. The Condition index (CI) was computed according to Riascos et al. (2012) as Cl=DW/ SL*100, where DW is the dry whole soft tissue weight (g) and SL is shell length (mm). Tissue preparation On the day of the analyses, the clams were thawed, their soft tissues extracted and homogenized with 0.1 M K-PO F buffer, pH 7.4. The resulting homogenates were centrifuged, first at 3000 rpm for 10 min. to obtain a post-nuclear fraction to measure glutathione concentration and then at 12 000 rpm for 20 min. to obtain a post-mi- tochondrial supernatant where the activities of the antioxidant enzymes were assayed. Biochemical analysis The protein concentrations of the tissue post-nuclear and post mitochondrial super- natant were measured according to Lowry et al. (1951). The concentrations of total glutathione (GSH) were assayed using the method described by Rahman et al. (2006) and the results were expressed as ng GSH/mg protein. The antioxidant enzyme ac- tivities were determined as follow: catalase (CAT) activity by Aebi (1984); superoxide dismutase (SOD) activity by Peskin and Winterbourn (2017); glutathione peroxidase (GRx) and glutathione-S-transferase (GST) activities by using kits, Glutathione Per- oxidase Cellular Activity Assay Kit (Cat. No CGP1) and GST Assay Kit (Cat. No Be ae Code Locality N E Sl Shkorpilovtsi 42.9564 28.0345 S2 Slanchev Bryag 42.6902 27.7171 «Kabakum Yarang _*Sy. Konstantin & Elena S3 Dyuni 42.3728 27.7163 “oct S4 Primorsko 42.2773 27.7321 S5 Arkutino 42.3445 27.7322 “sl Bryag Burgas _.§ S6 Ahtopol 42.1033 27.9259 gine 87 Kabakum 43.2656 28.0345 - ePrimorsko ee Ste 7) : Dabivpol ws Yai 4 $8 Sv. Konstantin and Elena 43.2332 28.0111 SS wk “ uv ba $9 Tsarevo 42.1611 27.8656 Figure |. Sampling localities along the Bulgarian Black Sea coast with geographical coordinates. 12 Almira Georgieva et al. / BioRisk 17: 169-178 (2022) CS0410), respectively, both purchased from Sigma-Aldrich Co. LLC (USA). The ob- tained values of all measured enzymes were expressed in U/mg protein. Statistical analysis Significance of differences of means between groups was determined using Student’s t-statistic. Factorial ANOVA was applied to analyze meaningful factors of variable de- pendences. In order to identify meaningful predictors of dependent variables multiple regression analysis was carried out. Results In the samples gathered in summer, no statistically significant differences were observed in the GSH levels and SOD and CAT activities of the clams by localities (Table 1). However, significant differences between the clams were observed in their GPx and GST activities. The lowest GPx activity was measured in samples from Skorpilovtsi (S1). The highest GST activity was measured in samples from Slanchev Bryag (S2) and the lowest — from Ahtopol (S6) (Table 1), which resulted in marked statistical differences between them. There were no significant differences in the measured OS biomarkers among local- ities in the autumn samples. ‘The same pattern was also present in the summer samples. The only exceptions were the GPx and GST activities (Table 1). In autumn, statisti- cally significant differences of GPx activities were observed between wedge clams from Arkutino (S5), Ahtopol (S6) and Kabakum (S7). There were also significant differ- ences between clams from Kabakum (S7), where the lowest GPx was measured, and Table |. Assessment of OS biomarkers (mean+sd) in clams from selected sites along the Bulgarian Black Sea coast in summer and autumn with significance of the differences of indicators at the same site between seasons (*p<0.05; **p<0=001; ***p<0.0001) and between localities (‘p<0.05; “p<0.001; “*p<0.0001 as the S(n) indicates the particular locality). Code Indicators GSH SOD CAT GPx GST (U/mg protein) Locality (ng/mg protein) (U/mg protein) (U/mg protein) (U/mg protein) Summer S1 Shkorpilovtsi 369.33 +80 3,24 41.3 2.56 +0.51 0.99726 +0.4 110.85? +32.5 $2 Slanchev Bryag 427.33 £44.8 4.56 +1.1 1.95 £0.40 1.89") 40.5 164.20'%° +28.0 S3 Dyuni 331.00 14.8 3.84 +1.0 2.48 0.634 1.8875! 41.1 112.03 +38.5 S4 Primorsko 396.55 +46.9 4.67*** +0.8 2.03 +0.27 Z.22HTTSI +0.8 149.94 437.2 S5 Arkutino 332.33 +20.2 4.51** +0.6 2.39 +0.33 2.028715! 40.7 151.68 +37.0 S6 Ahtopol 279.16** +51.1 5.52 +2.0 2.04 £0.36 1.89775! +0.5 85.83'? 445.4 Autumn S4 Primorsko 434.17 130.3 12w.06*** +3.7 2.47 +0.47 5.7207 Rots] 38 133.97 165758 35.9 $5 Arkutino 357.98 £124.5 10.30*** 44.3 2.74 £0.46 TASES? +16 100.2575 24.0 S6 Ahtopol 396.08** +49.8 9.58 +3.6 2.01 £0.35 3.117 41.6 60.027S455.99 415.1 S7 Kabakum 425.73 +136.8 10.52 +1.6 3.07 £0.17 1.327 +0.5 43.787154 +14 $8 Sv. Konstantin & Elena 378.80 +139.2 7.41 +1.0 2.95 £1.07 4.23 $1.23 61.244 +12.3 S9 Tsarevo 322.63 134.74 11.37 +2.0 2.32 £0.66 6.22757 +1.3 112.1775 +23.1 Antioxidant defense system in wedge clams from Bulgarian Black Sea N73 those from Tsarevo where the activity of GPx was significantly higher (Table 1). GST activities measured in samples from Primorsko (S4) differed significantly from those in the samples from the localities Ahtopol (S6), Kabakum (S7) and Sv. Konstantin and Elena (S8). Statistically significant differences were also present between the clams from Ahtopol (S6) and Tsarevo (S9) (Table 1). The measured antioxidant markers of clams were compared between the two sea- sons studied for the locations Primorsko (S4), Arkutino (S5), and Ahtopol (S6) (Table 1). No seasonal differences were present in the activities of CAT and GST. The level of GSH was raised in clams sampled in autumn from Ahtopol (S6). Similarly, the activi- ties of SOD and GPx were higher in the autumn samples from Primorsko (S4) and Arkutino (S5), compared to those in summer. In the present study, we calculated the CI of clams from the samples gathered in the two seasons from selected representative localities (Fig. 2). The CI of wedge clams gathered in summer were much more similar among the studied localities, than in the wedge clams collected in autumn. In June, significantly lower CI values were found in samples from Slanchev Bryag (0.65+0.04) and Duni (0.75+0.06). In wedge clams from the autumn samples CI values showed higher dis- similarity among localities, than in the summer (Fig. 2). Significantly lower CI values were calculated for wedge clams in the autumn samples from Ahtopol (0.43+£0.06) re — $1, 51, Dyuni 83 - Sl. Bryag $2 Arkutino S45 Arkutino S45 Kabakum 87 Shkorpilovtsi $1 Condition index = = = La | = Gi > Hh & * Primorsko $4 = Ahtopol $6 tk Primorsko S4 i Ahtopol S6 Hi BE: ee v. Konstantin & Elena $8 ri Tsarevo 89 I Summer Autumn Figure 2. Calculated condition indexes (meantSE) of wedge clams, gathered from localities along the Bulgarian Black Sea coast in summer and autumn with significance of the differences *p<0.05; ***<0.0001; S(n) indicates the particular locality from which the difference was significant). 174 Almira Georgieva et al. / BioRisk 17: 169-178 (2022) Table 2. Factorial ANOVA on dependence of wedge clams’ Condition Index (CI) on the factors Season and Locality. Sum Squares Degr. of Freedom Mean SS F P Intercept 0.52175 1 0.52171 1144.71 0.000000 Locality 0.01684 2 0.00842 18.47 0.000001 Season 0.00101 1 0.00101 2.22 0.140660 Locality*Season 0.02995 2 0.01497 32.86 0.000000 Error 0.02734 60 0.00045 Table 3. Results of regression analysis of the Condition Index (CI) of wedge clams (dependent variable) on the five studied antioxidants as predictor variables. b Std.Err. b t (86) p-value Intercept 0.04764 0.01831 2.60166 0.01092 GSH -0.00002 0.00002 -0.84815 0.39870 SOD 0.00048 0.00103 0.46778 0.64112 CAT 0.00977 0.00404 2.41578 0.01781 GPx 0.00017 0.00189 0.09154 0.92727 GST 0.00015 0.00007 2.01366 0.04717 and Kabakum (0.48+0.03). The CI values of the autumn samples from the remain- ing localities were significantly higher, even compared to the samples from all summer localities. Factorial ANOVA was applied to analyse the overall dependence of CI (i.e. health status) of the wedge clams on the two factors Season and Locality (Table 2). As can be seen from Table 2 the individual effect of Locality on CI was highly sig- nificant. The individual effect of Season alone was not significant. However, the joined effect of the two factors was highly significant, thus indicating that the state of the envi- ronment at a locality has a leading role in the condition and health status of the wedge clams, the seasonal effects being subordinate and specific for the particular locality. In order to identify meaningful predictors of the CI of the wedge clams among the studied antioxidant biomarkers, multiple regression analysis was carried out (Table 3) where CI was set as the dependent variable and the GSH concentration and the activ- ity of the antioxidant enzymes as predictor variables. The obtained results (Table 3) indicated the presence of significant correlative rela- tions between CI and the activity of CAT and GST. Thus, the two enzymes seem to be good predictor variables of the state and health of the wedge clams under changing environmental pressures. Discussion The results obtained in this study demonstrated the presence of significant variability in the pro/antioxidant processes in wedge clams sampled from representative loca- tions of the Bulgarian Black Sea coast. The variations were most probably the result of Antioxidant defense system in wedge clams from Bulgarian Black Sea 175 changes in the marine environmental conditions at the locations such as the geologi- cal characteristics of the area, proximity to large settlements and industrial sites, river inflows, etc. (Jenderedjian et al. 2007). In addition, seasonal variability was also pre- sent, which can depend on a large number of interrelated factors such as temperature, oxygen saturation, local hydrological cycle (which regulates productivity in an area), food availability, the metabolic status of the clams themselves, gonadal ripening and spawning (Soldatov et al. 2007; Nogueira et al. 2017), as well as the concentrations of metals (Reis et al. 2017) and PAHs (Koudryashova et al. 2019) in seawater. Data re- ported here indicated that glutathione-related enzymes GPx and GST were the mainly affected biomarkers. GSH itself, which is a co-substrate of both enzymes, did not show significant changes among localities. Although GSH plays a central role in antioxidant systems, this marker alone is not sufficient to indicate the OS rate due to the complex regulation of its homeostasis (Manduzio et al. 2005). Hence, changes in GSH level may be transient, and GSH should be measured in parallel with GPx and glutathione reductase, which are involved in its turnover, so that environmental stressors can be more precisely detected (Kulinskii and Kolesnichenko 2009). GPx, which is catalyzing the reduction of both hydrogen peroxide and organic hydroperoxides, provides an ef- fective protection against oxidative damage and free radicals (Ighodaro and Akinloye 2018). Increased activity of GPx has been found in bivalves in response to marine pol- lution with heavy metals (Mansour et al. 2020) and PAHs (Moreira et al. 2020). In this study, increased GPx activities were found in clams only from Primorsko and Arkutino (both seasons), which was probably due to local environmental characteristics. This assumption was also supported by the finding that in the wedge clams from these loca- tions an increased GST activity was measured, compared to the other locations. GST plays a principal role in the metabolism of xenobiotics by catalyzing their conjugation with GSH, making them more water-soluble and easily eliminated from the body. GST also plays a role in the protection against OS having peroxidase activity (Hayes et al. 2005). Furthermore, GST can regenerate S-thiolated proteins, which are the consequence of OS development in cells. Studies have reported increased activity and expression of various isoforms of GST in bivalves caused by xenobiotics such as PAHs, polychlorinated biphenyls, DDT, dioxins, metals and microplastics (Vidal-Lifan et al. 2010, 2014; O'Donovan et al. 2018). In our study, the highest GST activity among all samples was found in the soft tissue of clams collected from Slanchev Bryag in summer. This was most probably due to intensive touristic activities and pollution at this local- ity, which led to activation of phase II metabolism, and antioxidant defense system to overcome pollution. Furthermore, the wedge clams from this locality had low CI, which suggested decline in their physiological condition and proved the presence of significant biorisk for the marine biota. The activity of GST is known as a biomarker of exposure to pollutants since the enzyme has exhibited highly significant positive correlations with toxicants’ amount and the use of GST as indicator is an advantage because of the regularity of its amendment and low seasonal variability (Vidal-Lifan et al. 2010). Indeed, in the present research we did not find seasonal changes in the activity of GST in the clams from one and the same locality. In our study, statistically 176 Almira Georgieva et al. / BioRisk 17: 169-178 (2022) significant differences in summer and autumn were found for GSH, SOD and GPx for a number of localities (Table 1). In general, the activity of the measured enzymes as well as the concentration of GSH were higher in autumn. This was most probably due to the activation of the antioxidant defense system of wedge clams after the summer touristic flow in locations known for their high touristic activity. In the present study, CI of clams was calculated in samples from selected locali- ties in the two seasons, as an indicator of the physiological condition and health of individual wedge clams. ‘The effects of the Locality and Season, as general factors, were assessed using Factorial ANOVA. ‘The result indicated highly significant individual ef- fect of Locality on CI, while the Season had not such effect. However, the joint effect of Locality and Season was highly significant. This result can be interpreted as a proof of the leading role of the state of the environment in a particular locality on the health of wedge clams, while seasonal effects seemed to be subordinate and dependent on the presence of different environmental factors and changes of a seasonal character at the particular locality. It was also proved that activity of the antioxidant enzymes CAT and GST was correlated with CI and they seemed to be good predictor variables of the presence of risks for the state and health of the wedge clams under variations of local environmental pressures. Conclusion The established variability in the antioxidant enzyme complex of D. trunculus clearly reflected the state of the particular environmental conditions in the studied locations of the Bulgarian Black Sea coast. The results indicated that the wedge clams demon- strated good and effective antioxidant defense, which was activated by increased envi- ronmental pressures at different locations, indicating their ability to cope with induced oxidative stress and adapt to local conditions, thus contributing to ecosystem health maintenance. Further work is obviously needed to prove the capacity of antioxidant enzymes as reliable bioindicators of ecosystem changes. Acknowledgements This work was supported by grant KI1-06-H31/6 of National Science Fund, Bulgaria. We thank Clam Shells Ltd. for assistance in providing material from wedge clams. References Aebi H (1984) Catalase in vitro. Methods in Enzymology 105: 121-126. https://doi. org/10.1016/S0076-6879(84)05016-3 Andrady AL (2011) Microplastics in the marine environment. Marine Pollution Bulletin 62(8): 1596-1605. https://doi.org/10.1016/j.marpolbul.2011.05.030 Antioxidant defense system in wedge clams from Bulgarian Black Sea LA Bat L, Arici E, Oztekin A (2018) Human health risk assessment of heavy metals in the Black Sea: Evaluating mussels. Current World Environment 13(1): e17. https://doi.org/10.12944/ CWE.13.1.03 Georgieva S, Stancheva M, Makedonski L (2016) Investigation about the presence of organo- chlorine pollutants in mussels from the Black Sea Bulgaria. Analele Universitatii Ovidius Constanta. Seria Chimie 27(1): 8-12. https://doi.org/10.1515/auoc-2016-0006 Hayes JD, Flanagan JU, Jowsey IR (2005) Glutathione transferases. The Annual Review of Pharmacology and Toxicology 45: 51-88. https://doi.org/10.1146/annurev.pharm- tox.45.120403.095857 Honda M, Suzuki N (2020) Toxicities of Polycyclic Aromatic Hydrocarbons for Aquatic Ani- mals. International Journal of Environmental Research and Public Health 17(4): e1363. https://doi.org/10.3390/ijerph17041363 Ighodaro OM, Akinloye OA (2018) First line defence antioxidants-superoxide dismutase (SOD), catalase (CAT) and glutathione peroxidase (GPX): Their fundamental role in the entire antioxidant defence grid. Alexandria Journal of Medicine 54(4): 287-293. https:// doi.org/10.1016/j.ajme.2017.09.001 Jenderedjian K, Hakobyan S, Jenderedjian A (2007) Use of benthic invertebrates as indica- tors of pollution origin in agricultural and urban areas. In: Ebel A, Davitashvili T (Eds) Air, water and soil quality modeling for risk and impact assessment, 21-220. https://doi. org/10.1007/978-1-4020-5877-6_20 Koudryashova Y, Chizhova T, Tishchenko P, Hayakawa K (2019) Seasonal Variability of Polycy- clic Aromatic Hydrocarbons (PAHs) in a Coastal Marine Area in the Northwestern Region of the Sea of Japan/East Sea (Possiet Bay). Ocean Science Journal 54(4): 635-655. https:// doi.org/10.1007/s12601-019-0031-9 Kulinskii VI, Kolesnichenko LS (2009) Glutathione system. I. Synthesis, transport, glutathione transferase, glutathione peroxidase. Biomeditsinskaia Khimiia 55(3): 255-277. https:// doi.org/10.1134/S1990750809020036 Lowry O, Rosenbrough N, Farr A, Randal R (1951) Protein Measurement with the FolinPhenol Reagent. The Journal of Biological Chemistry 193(1): 265-275. https://doi.org/10.1016/ S$0021-9258(19)52451-6 Manduzio H, Rocher B, Durand F, Galap C, Leboulenge F (2005) The point about oxidative stress in mollusks. Invertebrate Survival Journal 91(2): 91-104. https://hal-normandie- univ.archives-ouvertes.fr/hal-01738557 Mansour C, Guibbolini M, Rouane Hacene O, Saidane Mosbahi D, Risso-de Faverney C (2020) Oxidative stress and damage biomarkers in clam Ruditapes decussatus exposed to a polluted site: The reliable biomonitoring tools in hot and cold seasons. Archives of Environmental Con- tamination and Toxicology 78(3): 478-494. https://doi.org/10.1007/s00244-020-007 13-2 Maranho LA, DelValls TA, Martin-Diaz ML (2015) Assessing potential risks of wastewater discharges to benthic biota: An integrated approach to biomarker responses in clams (Ru- ditapes philippinarum) exposed under controlled conditions. Marine Pollution Bulletin 92(1—-2): 11-24. https://doi.org/10.1016/j.marpolbul.2015.01.009 Mojiri A, Zhou JL, Robinson B, Ohashi A, Ozaki N, Kindaichi T, Farraji H, Vakili M (2020) Pesticides in aquatic environments and their removal by adsorption methods. Chemos- phere 253: e126646. https://doi.org/10.1016/j.chemosphere.2020.126646 178 Almira Georgieva et al. / BioRisk 17: 169-178 (2022) Moreira LB, Sasaki ST, Taniguchi S, Peres TE, Figueira RCL, Bicego MC, Marins RV, Costa- Lotufo LV, Abessa DMS (2020) Biomarkers responses of the clam Anomalocardia flexuosa in sediment toxicity bioassays using dredged materials from a semi-arid coastal system. Heliyon 6(5): e04030. https://doi.org/10.1016/j.heliyon.2020.e04030 Nogueira L, Mello DF, Trevisan R, Garcia D, da Silva Acosta D, Dafre AL, de Almeida EA (2017) Hypoxia effects on oxidative stress and immunocompetence biomarkers in the mussel Pernaperna (Mytilidae, Bivalvia). Marine Environmental Research 126: 109-115. https://doi.org/10.1016/j.marenvres.2017.02.009 O’Donovan S, Mestre NC, Abel S, Fonseca TG, Carteny CC, Cormier B, Bebianno MJ (2018) Ecotoxicological effects of chemical contaminants adsorbed to microplastics in the clam Scrobicularia plana. Frontiers in Marine Science 5: e143. https://doi.org/10.3389/ fmars.2018.00143 Peskin AV, Winterbourn CC (2017) Assay of superoxide dismutase activity in a plate assay us- ing WST-1. Free Radical Biology & Medicine 103: 188-191. https://doi.org/10.1016/j. freeradbiomed.2016.12.033 Rahman I, Kode A, Biswas SK (2006) Assay for quantitative determination of glutathione and glutathione disulfide levels using enzymatic recycling method. Nature Protocols 1(6): 3159-3165. https://doi.org/10.1038/nprot.2006.378 Reis PA, Salgado MA, Vasconcelos V (2017) The spatial and seasonal variation of trace metals in coastal seawater and soft tissue of Chthamalus montagui around the northwest coast of Portugal. Ocean Science Journal 52(2): 207-219. https://doi.org/10.1007/s12601-017- 0013-8 Riascos JM, Avalos CM, Pacheco AS, Hrilmayer O (2012) Testing stress response of the bivalve Ptotothaca thaca to El Nifo- La Nina thermal conditions. Marine Biology Research 8(7): 654-666. https://doi.org/10.1080/17451000.201 1.653367 Soldatov AA, Gostyukhina OL, Golovina IV (2007) Antioxidant enzyme complex of tissues of the bivalve Mytilus galloprovincialis lam. under normal and oxidative-stress conditions: A review. Applied Biochemistry and Microbiology 43(5): 556-562. https://doi.org/10.1134/ S$0003683807050092 Steinberg CEW (2012) Environmental Stresses: Ecological Driving Force and Key Player in Evolution. Stress Ecology: Environmental stress as ecological driving force and key player in evolution. Springer, Dordrecht, 480 pp. https://doi.org/10.1007/978-94-007-2072-5 Vidal-Lifan L, Bellas J, Campillo JA, Beiras R (2010) Integrated use of antioxidant enzymes in mussels, Mytilus galloprovincialis, for monitoring pollution in highly productive coastal areas of Galicia (NW Spain). Chemosphere 78(3): 265-272. https://doi.org/10.1016/j. chemosphere.2009.10.060 Vidal-Lifan L, Bellas J, Etxebarria N, Nieto O, Beiras R (2014) Glutathione-S-transferase, glu- tathione peroxidase and acetylcholinesterase activities in mussels transplanted to harbour areas. The Science of the Total Environment 470-47 1: 107-116. https://doi.org/10.1016/j. scitotenv.2013.09.073