Abstract

Marine organisms are exposed to a pH decrease and to alteration of carbonate chemistry due to ocean acidification (OA) that can represent a source of oxidative stress which can significantly affect their antioxidant defence systems efficiency. The polychaetes Platynereis dumerilii and P. massiliensis (Nereididae) are key species of the benthic community to investigate the effect of OA due to their physiological and ecological characteristics that enable them to persist even in naturally acidified CO2 vent systems. Previous studies have documented the ability of these species to adapt to OA after short- and long-term translocation experiments, but no one has ever evaluated the basal antioxidant system efficiency comparing populations permanently living in habitat characterized by different pH conditions (acidified vs. control). Here, individuals of both Platynereis species, sampled from a natural CO2 vent system and from a nonventing “control” site in three different periods (April 2016, October 2016, and February 2017), were compared highlighting signals which suggested the ability of both species to acclimatize to high pCO2–low pH with slight seasonal variations of their antioxidant efficiency and the absence of disturbances of the oxidative status of Platynereis spp. tissues.

1. Introduction

World climatic alterations are mainly driven by atmospheric CO2 partial pressure (pCO2) increase, as a consequence of the anthropogenic activity, which is predicted to reach 800 ppm by the end of the current century [15]. This change at the atmospheric level also affects the ocean surface through the phenomenon of ocean acidification (OA): an increase of the dissolved CO2, alteration of seawater carbonate chemistry, and the consequent reduction of the ocean pH. Ocean surface pH has fallen by about 0.1 units since the beginning of the industrial era and is expected to further decrease by 0.3-0.5 units by the end of the current century [1, 4, 5]. Marine organisms are continuously exposed to a range of environmental parameters, such as pH, salinity, and temperature, varying over temporal and spatial scales, which may represent a source of oxidative stress that entails Reactive Oxygen Species (ROS) additional production. ROS, endogenous and highly reactive oxygen-bearing molecules, are commonly produced at low concentrations during several natural cellular pathways of aerobic metabolism, and under basal and stable conditions their adverse effects are prevented by antioxidant defence systems (i.e., low molecular weight scavengers and enzymes). In stressful conditions, this balance may be altered leading to uncontrolled ROS formation that translates into cellular oxidative damage against biological macromolecules including lipids, proteins, and DNA, impairing normal cellular functions. This unbalance in favour of oxidants is termed “oxidative stress”. In the last years, scientific literature has provided evidence that global climate change, especially OA, affects antioxidant systems efficiency of several marine organisms [620]. The biological effects of low pH–high pCO2 have been investigated through not only laboratory/mesocosm experiments, but also studying natural volcanic CO2 vent systems which occur in different parts of the world (e.g., [2128]). Such “natural laboratories” have so far provided environmentally realistic overviews of the conditions with which marine organisms will interface according to the near-future OA predictions [5, 29] and represent an important tool to detect information about putatively tolerant species to OA and their capability to modulate the antioxidant system in response to pH variations. Despite the growing interest of the scientific community in this research field, there is still little knowledge about how the acidification processes affect antioxidant defence systems of the benthic biota.

Within benthic community, polychaetes represent a key group in marine habitat and are often used as bioindicators in monitoring programs for their high sensitivity to metal exposure [3033] and anthropogenic pressure [34, 35]. Due to their physiological and ecological characteristics, polychaetes are among the most abundant invertebrates under low pH conditions, such as along the naturally acidified CO2 vent system of the Castello at the Ischia Island [22, 3639], and various studies investigated the role of antioxidant systems in response to low pH–high pCO2 [12, 1618, 24, 40, 41]. The polychaete Platynereis dumerilii (Audouin and Milne-Edwards, 1834) (Nereididae) represents a key species: it showed high tolerance to environmental stress [42, 43], including low pH [38], and for this reason it was employed as model organism in two recent in situ transplant experiments, carried out to evaluate the effects of OA on the oxidative sensitivity of populations living inside and outside the vent area of Ischia [12, 24]. The short-term translocation (only 5 days) displayed a true local adaptation to low pH conditions, with an exclusive genotype apparently restricted to the acidified areas of Castello Aragonese (Ischia) characterized by higher metabolic rate, measured as oxygen consumption [24]. The genotype of the acidified areas was later identified as the only known sibling species of Platynereis dumerilii, P. massiliensis (Moquin-Tandon, 1869): morphologically indistinguishable species in its adult (nonreproductive) stage, but characterized by a completely different reproductive biology [4446]. A recently published study that combined together genetic and reproductive biology analyses revealed that both Platynereis species actually represented two different complexes of siblings [47]. Based on preliminary genetic results, the Castello vent site of Ischia appeared dominated by the brooding P. massiliensis sibling, while the control site by the broadcasting species P. dumerilii [24, 46, 47]. In the long-term translocation (30 days) the antioxidant sensitivity of different polychaete species (Platynereis dumerilii, Polyophthalmus pictus, and Syllis prolifera) and their antioxidant capacity to counteract oxyradicals formation in control and low pH conditions was evaluated, highlighting as the population of Platynereis originating from the vent showed higher constitutive antioxidant efficiency [12], which may allow them to cope with short-term and chronic exposure to higher oxidative pressure without further enhancement of antioxidant defences [12]. From these results, a hypothesis of long-term adaptation of the Platynereis vent-inhabiting population emerged, suggesting the need of this species for greater antioxidant protection in conditions of chronic oxidative exposure to low pH–high pCO2 [12]. Recent laboratory experiments on both species also highlighted a differentiation on the expression of some target genes involved in the oxidative metabolism as a result of the exposure to different pH conditions [48, 49]. Specimens of Platynereis dumerilii from a control site near the vents, Sant’Anna, showed significant lower levels of NADH dehydrogenase mRNA expression compared to P. cfr massiliensis from Castello acidified sites confirming, as already stated by Calosi et al. [24], that living under acidified conditions entails a higher energetic consumption and metabolic rate [48]. In line with Calosi et al. [24], translocation experiment did not show significant effect for this gene expression in Platynereis dumerilii but, in contrast, a significant downregulation of P. cfr massiliensis NADH dehydrogenase from low pH to control conditions was observed indicating a reduction in the oxidative metabolism of this species [48]. Differently to what was previously asserted by Lucey et al. [46], recent phylogenetic analysis carried out in the frame of a PhD thesis on Platynereis spp. [50], samples collected in the south-acidified areas of Castello Aragonese and the control zone of Sant’Anna rocks highlighted a less evident spatial segregation of the two Platynereis species between the two sites. For this reason, it was not possible to consider the Castello vent area as an exclusive domain zone for P. massiliensis, as well as the control area of Sant’Anna for P. dumerilii.

Considering the preliminary results on antioxidant efficiency of putative P. dumerilii and the presence of both species in control and acidified sites, the aim of this study was to provide new insights about the basal levels of the antioxidant system in Platynereis spp., comparing populations living at different pH conditions (acidified vs. control). This comparison between populations was carried out in three different periods: April 2016, October 2016, and February 2017, in relation to different temperature conditions. The oxidative effects of different pH levels were evaluated thorough analysis of single antioxidant activities, such as catalase (CAT), glutathione S-transferases (GST), glutathione reductase (GR), and Se-dependent and Se-independent glutathione peroxidases (GPx), that can be very sensitive in revealing a prooxidant condition [51] and thorough their integration with total oxyradical scavenging capacity (TOSC), which quantify the capability of Platynereis spp. to neutralize different forms of oxyradicals including peroxyl radicals (ROO), hydroxyl radicals (HO), and peroxynitrite (HOONO) [52].

2. Materials and Methods

2.1. Study Areas, Sample Collection, and Processing

The study was conducted at the Ischia island (Gulf of Naples, Italy), a volcanic island well known for the presence of numerous submarine CO2 vent systems [29, 53, 54], including the area of Castello Aragonese on the north-eastern side (40° 43.84 N, 13° 57.08 E) as the first vent system studied in the world [29]. Gas bubbles, composed by 90-95% CO2, 3-6% N2, 0.6-0.8% O2, 0.2-0.8% CH4, with no sulphur, are released at ambient seawater temperature at about 1.4 x 106 l d−1 [21]. The salinity of the water (38) and total alkalinity (2.5 mequiv Kg−1) are relatively uniform [21, 38]. Gas emissions, which occur between 0.5 and 3.0 m depth, create a gradient of pH on both the north and south sides of the Castello islet. According to Kroeker et al. [22], three different pH zones can be identified along a rocky reef of approximately 150 m in length on each side of the islet: a control area with normal pH conditions and no venting activity (N1 and S1), an intermediate area with moderate vent activity and low pH conditions (N2 and S2), and a high venting activity area characterized by extreme low pH conditions (N3 and S3) (Figure 1). The acidified S2–S3 stations where Platynereis collection was performed showed a wide range of pH variability as reported by Ricevuto et al. [38, Supplement material]. The control site, called Sant’Anna rocks, characterized by a very stable pH value, which ranged around a mean of 8.01, is located within the Cartaromana Bay, approximately 600 m from the south side of Castello Aragonese [12] (Figure 1).

Platynereis spp. samples were collected in the south-acidified sites of Castello (named as S3 and S2 in previous papers e.g., [36], or as low pH and extreme low pH in [22]) and in the normal/control pH area of Sant’Anna rocks. Samplings were carried out in four different periods: April 2016, October 2016, and February 2017 (mean monthly seawater temperatures 16.5°C, 21.9°C, and 14.6°C, respectively). Worms were collected in each sampling site and period by sampling local macroalgae belonging to the species Halopteris scoparia, Jania rubens, Dictyota spp., and Cladophora spp. where these polychaetes live associated. Macroalgal thalli were collected in cotton fabric bags by snorkelling and SCUBA diving at 0.5–2 m depth. After collection, samples were transported to the Villa Dohrn-Benthic Ecology Center (approx. 4 km from the Castello area) inside cool boxes within one hour. Once in the lab, algal thalli were sorted and the Platynereis spp. species was identified thanks to the typical sinuous swimming movement, immediately transferred into separated 1.5 ml microcentrifuge tubes (pooling approximately 5-10 individuals per Eppendorf, according to the body mass of the collected samples), frozen, and temporarily stored at -80°C until the transport to the Laboratory of Ecotoxicology and Environmental Chemistry of Ancona (Italy) for the antioxidant analyses.

2.2. Analyses of Antioxidants and Total Oxyradical Scavenging Capacity

For the analysis of the antioxidant enzyme activities, pools were homogenized (1:10 w:v) in 100 mM of potassium phosphate buffer (pH 7.5) containing NaCl (1.5%), 0.1 mg mL−1 phenylmethylsulphonyl fluoride (PMSF), 0.1 mg mL−1 bacitracin and 0.008 TIU mL−1 aprotinin as protease inhibitors. After centrifuging at 100,000 x g for 70 min at 4°C, supernatants were collected and used for the subsequent analyses. Enzymatic activity measurements including catalase (CAT), glutathione S-transferases (GST), glutathione reductase (GR), and glutathione peroxidases (GPx) were carried out using a Varian (Model Cary 3) spectrophotometer at the constant temperature of 18°C according to Bocchetti et al. [30]. CAT activity was determined by the decrease in absorbance at 240 nm (ε = 0.04 mM−1 cm−1) due to H2O2 consumption (12 mM) in 100 mM K-phosphate buffer (pH 7.0). GST activity was quantified at 340 nm using 1-chloro-2,4 dinitrobenzene (CDNB) as substrate (ε = 9.6 mM−1 cm−1). The assay conditions were 100 mM potassium phosphate buffer (pH 6.5), 1.5 mM CDNB, and 1.5 mM GSH. GR activity, also known as glutathione-disulphide reductase (GSR), was measured spectrophotometrically at 340 nm following the oxidation of NADPH during the reduction of GSSG (extinction coefficient, ε = 6.22 mM−1 cm−1). The assay was carried out in 100 mM potassium phosphate buffer (pH 7.0), 1 mM GSSG and 60 μM NADPH. The activity of Se-dependent and Se-independent GPx forms was determined in two enzymatic assays in which GSSG is converted to the reduced form GSH. The consumption of NADPH was quantified as decrease of absorbance at 340 nm (ε = 6.22 mM−1 cm−1) in 100 mM K-phosphate buffer pH 7.5, 1 mM EDTA, 1 mM dithiothreitol (DTT), 2 mM GSH, 1-unit glutathione reductase, 0.24 mM NADPH, and 0.8 mM cumene hydroperoxide as substrate.

For the total oxyradical scavenging capacity (TOSC) analysis, polychaetes were homogenized following the same protocol reported above, with 0.5 μg mL−1 pepstatin as additional protease inhibitor and without phenylmethylsulphonyl fluoride (PMSF). This assay was based on the capability of cellular antioxidants to reduce the oxidation of α-keto-γ-methiolbutyric acid (KMBA), and the consequent formation of ethylene gas, in presence of artificially generated oxyradicals. The ethylene formation was monitored at 12 min time intervals by gas-chromatographic analyses and the TOSC values were calculated from the equation: TOSC = 100 – (/ x 100), where and were the integrated areas calculated under the kinetic curve produced during the reaction course for sample (SA) and control (CA) reactions, respectively [52].

In order to obtain the specific antioxidant activity and TOSC values, data were normalized with the relative protein concentration according to Lowry method [55] by using Bovine Serum Albumin (BSA) as standard.

2.3. Statistical Analysis of Data

Permutational multivariate analysis of variance (PERMANOVA) and pairwise nonparametric tests with square root transformation, Euclidean distance, and 9,999 number of permutations were conducted to test the differences (α = 0.05). “Population origin” and “sampling period” were considered as fixed factors with two and three levels, respectively (sites: Castello and Sant’Anna; periods: April 2016, October 2016, February 2017), to test the response temporal variation between population. Multivariate principal component analysis (PCA) was applied to visualize the relationships among the different populations/sampling periods and all statistical analyses were performed using PRIMER/PERMANOVA v 6 [56].

3. Results

Results are showed in both Table 1 and Figure 2; the activity of GR significantly differed between the two Platynereis spp. populations in April and February (pairwise comparison p < 0.05, see asterisks Table 1). Sant’Anna population showed significant differences between April–October and April–February (pairwise comparison p < 0.05, see letters Table 1). The highest enzyme activity was recorded in April (GR = 54.53 nmol/min/mg prt) and the lowest one in February (GR = 6.98 nmol/min/mg prt). On the contrary, GR activity was higher in October in organisms sampled from Castello with values comparable to those observed in Sant’Anna population (SA = 30.43 nmol/min/mg prt; CA = 32.11 nmol/min/mg prt), while constant values were observed for this population in the other months. Concerning the capability to neutralize HOONO, the Platynereis spp. populations showed significant differences in October and February between the investigated sites (Figure 2). While the total oxyradical scavenging capacity toward HOONO in organisms from Castello was constant over different sampling periods, conversely, more marked variation was observed for Platynereis specimens from Sant’Anna site, with the highest TOSC efficiency measured in October (HOONO = 786.5 UTosc/mg prt) and the lowest one in February (336.42 UTosc/mg prt). Concerning other antioxidant parameters, in worms sampled from Sant’Anna, CAT exhibited the highest activity in April, followed by a significant decrease over the examined periods; on the contrary, Castello specimens showed constant values with some differences (p < 0.05) between specimens in October and February (Figure 2). Similar trend, with the highest enzymatic activity in April, was also observed for the GST in both populations. On the contrary, the activity of Se-dependent and Se-independent GPx was lower in April and increased in October and February in organisms sampled from both sites. TOSC values toward ROO did not show significant differences, while against HO limited variations were observed in organisms sampled from Castello with a significant decrease in February.

PCA analysis provided a two-dimensional pattern explaining 82.5% of the total variance (56.2% and 26.3% in the first and second axes, respectively) (Figure 3). However, despite some differences occurred mainly for single antioxidants activity between populations and periods, no clear groups’ separation occurred (Figure 3).

4. Discussion

This study represents the first attempt of a background analysis of the antioxidant parameters of Platynereis spp. (putative different sibling species) living in different pH conditions (normal and acidified conditions). Specimens morphologically identified as Platynereis dumerilii were already employed as model organisms for some transplant experiments to investigate the prooxidant effect of ocean acidification [12, 24, 48, 49], but the basal level efficiency of the antioxidant defence systems in terms of species sensitivities and seasonality of populations submitted to natural acidified and normal pH conditions was never evaluated.

Our overall results indicated that the two populations showed few differences, probably based on seasonal-related conditions in their habitat, in particular for glutathione reductase activity and the total oxyradical capability to counteract HOONO radical. GR catalyses the reduction of glutathione disulphide (GSSG) to the active form of glutathione (GSH) which is a fundamental molecule for preventing oxidative stress and maintaining the reduced environment of cell, using NADPH as cofactor. The GR activity of specimens sampled from control site (Sant’Anna population) showed the highest activity in April and a rapid decrease until February. Conversely, the Castello vent population showed similar trend and values during the whole year, suggesting more stable environmental conditions during the examined periods. Similar considerations were also supposed for TOSC toward HOONO, which showed for specimens collected in Sant’Anna a maximum capability to counteract radical species in October and a minimum one in February, while the vent population displayed constant activity during the whole sampling period.

Based on limited variations, our findings did not show a clear correlation between pH conditions and alteration of the oxidative status in Platynereis spp. populations, which was instead already observed in previous ocean acidification experiments with several invertebrate species [6, 8, 1012, 1416, 41, 48, 49, 57]. An in situ transplant experiment of 30 days into naturally acidified conditions carried out with the fan worm Sabella spallanzanii (Sabellidae) highlighted a significant decrease of enzymatic activities of CAT and GPx and the impairment of the overall capability to neutralize hydroxyl radicals (HO) [41]. The effect of pH decrease was also investigated in the polychaete Diopatra neapolitana; after 28 days of exposure to low pH levels higher enzymatic activity (CAT, SOD, GSTs) and oxidative alterations were recorded [16]. The results of TOSC assay performed after an in situ reciprocal transplant experiment in the Castello vents on specimens morphologically identified as Platynereis dumerilii displayed insights of long-term adaptation of the vent population to greater prooxidant challenge [12]. Conversely to our findings, Ricevuto et al. [12] observed that Platynereis vent specimens showed more elevated basal antioxidant efficiency toward ROO and HOONO when compared to the population collected in control pH conditions, suggesting the need of a greater antioxidant protection in conditions of chronic oxidative exposure [12]. A different response to pH conditions was also highlighted by Wäge et al. [48] in whose study P. cfr massiliensis from acidified sites showed a marked upregulation of gene expression involved in the energy metabolism compared to P. dumerilii from control site. This inconsistency with our results might be due to a mixing of the two species (or complexes) in the studied areas, Castello Aragonese and Sant’Anna rocks, only 600 meters away from each other. The two species would therefore be separated exclusively by a chemical barrier (pH) for which nonsignificant effects were detected at least in P. dumerilii [48].

The PCA analysis suggested that the antioxidant response of Sant’Anna population was more differentiated when compared with specimens collected in the acidified areas of Castello, without a clear trend in the different periods (Figure 3). Both statistical analyses supported that the two populations seemed to have different temporal-related trends of the antioxidant defence systems, even if the lowest mean values of antioxidant capacity were mainly recorded in the period characterized by the lowest water temperature (February). This phenomenon was most evident in the Sant’Anna population, which showed a marked decrease of antioxidant activities of CAT, GST, and GR and a lower scavenging capacity toward ROO and HOONO, in February. The antioxidant defence systems’ efficiency of marine organisms can be influenced by several environmental factors, such as annual fluctuations in solar irradiance, changes in water temperature, fluctuating oxygen concentration, and exposure to chemical pollutants [5863]. The lowest temperature of the winter season could entail lower prooxidant pressure and, as a consequence, the need of a decreased antioxidant efficiency, as reported for the European eel, Anguilla anguilla, [64], in order to counteract the increase of environmental ROS formation during the summer. High temperature increased ROS production and the consequent enhancement of antioxidant enzymes’ activity also in the mussel species Mytilus coruscus [57].

In this study, Platynereis spp. population from control pH conditions displayed a higher temporal variability, showing the need to modulate the redox response to keep the oxidative stress level of the tissues under control during different periods of the year. Conversely, the overall ability of the vent population to maintain stable levels of antioxidant defences, regardless the period of the year and seasonal-related trends, suggested that natural enhancement of environmental prooxidant conditions was balanced with slight changes of individual antioxidants. The different native pH conditions, over three sampling periods, were not translated into significant differences between populations in the other antioxidant biomarkers analysed, and the similar antioxidant responses highlighted between the two populations confirmed the high tolerance of these species (or complex of sibling species, [47]). The long-term exposure to moderately elevated pCO2 conditions, such as those expected in global climate change scenarios at the end of this century [1], could minimally affect the cellular redox status, as already observed in two marine bivalve species, Crassostrea virginica and Mercenaria mercenaria [9].

5. Conclusions

In conclusion, this study provided the first baseline for a direct characterization of antioxidant responses in Platynereis spp. from naturally acidified and control pH conditions. The different pCO2 and pH levels of the studied habitats did not seem to strongly affect Platynereis spp. defence systems’ efficiency and no one of the two populations stood out for a stronger or lower antioxidant capacity. Platynereis spp. specimens appeared able to acclimatize to low pH conditions with slight seasonal variations of the antioxidant defence systems, with enzymatic activity and TOSC kept constant throughout the year. The inconsistency with previous studies suggested the need to further investigate the seasonal variation of the basal antioxidant systems efficiency of genetically characterized Platynereis specimens collected in vent and nonvent sites.

Data Availability

The individual raw data used to support the findings of this study are available from the corresponding author upon request.

Ethical Approval

All applicable international, national, and/or institutional guidelines for the care and use of animals were followed.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

We thank Ilaria Lanzoni for the assistance with the analyses at the Laboratory of Ecotoxicology and Environmental Chemistry of Polytechnic University of Marche, Ancona. This work was part of the PhD thesis of Giulia Valvassori, supported by the Stazione Zoologica Anton Dohrn (Napoli, Italy) fellowship.