Effects of hypoxia on survival, behavior, metabolism and cellular damage of Manila clam (Ruditapes philippinarum)

Qiao Li; Song Sun; Fang Zhang; Minxiao Wang; Mengna Li

doi:10.1371/journal.pone.0215158

Abstract

The Manila clam Ruditapes philippinarum has become a common and dominant macrobenthic species in coastal areas of the northwestern Pacific and temperate waters of Europe; it is also a major cultured shellfish, with annual worldwide production exceeding 3.3 million tonnes. This species faces greater risk of exposure to hypoxia as eutrophication worsens throughout its coastal habitats; however, its tolerance to hypoxia remains unclear, and the toxicological indicators including LC₅₀ and LT₅₀ have not yet been assessed. Previous studies on the effects of hypoxia on marine benthos have focused largely on functional responses, such as metabolism and gene expression, leaving potential structural damage to the mitochondria or the cells unknown. In this study we assessed the effects of hypoxia on Manila clam in terms of survival, behavior, metabolism and cellular damage, using a newly designed automated hypoxia simulation device that features exceptional accuracy and good stability. The clams exhibited strong tolerance to hypoxia as the 20-day LC₅₀ for dissolved oxygen (DO) was estimated to be 0.57 mg L^-1, and the LT₅₀ at 0.5 mg L^-1 DO was 422 hours. Adaptations included fewer buried clams and a depressed metabolism, while the unexpected rise in the activities of key enzymes involved in glycolysis may indicate a diverse strategy of shellfish under hypoxia. Cellular damage was observed as collapse of the mitochondrial cristae and both cellular and mitochondrial vacuolization. This multi-level study complements and updates our knowledge of the effects of hypoxia on marine benthos, by improving our understanding of the potential for marine ecological transformation under hypoxic conditions and providing useful information for Manila clam farming.

Citation: Li Q, Sun S, Zhang F, Wang M, Li M (2019) Effects of hypoxia on survival, behavior, metabolism and cellular damage of Manila clam (Ruditapes philippinarum). PLoS ONE 14(4): e0215158. https://doi.org/10.1371/journal.pone.0215158

Editor: Ismael Aaron Kimirei, Tanzania Fisheries Research Institute, UNITED REPUBLIC OF TANZANIA

Received: October 23, 2018; Accepted: March 27, 2019; Published: April 18, 2019

Copyright: © 2019 Li et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the manuscript and its Supporting Information files.

Funding: This study was supported by the NSFC—Shandong Joint Fund for Marine Ecology and Environmental Sciences [grant number U1606404]; the Strategic Priority Research Program of the Chinese Academy of Sciences [grant number XDA19060201, XDA23050501]; National Key R&D Program of China [grant number 2017YFC1404405]; and Taishan Scholars Project to Song Sun. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

As the eutrophication of coastal waters intensifies, both the extent and strength of oxygen depletion (hypoxia) increase [1,2]. Though certain mobile species can flee hypoxic zones, less-tolerant organisms with poor mobility often suffer mass mortalities under severe hypoxia [2–4], often leading to further harmful ecological changes, such as jellyfish blooms [5,6] and other forms of ecosystem degradation [7,8].

The commercially important Manila clam Ruditapes philippinarum has strong adaptability to a wide range of salinities and water temperatures; the species has a widespread natural distribution in coastal areas of the northwestern Pacific and has become established on the Atlantic coast of Europe. As a common and dominant species in many coastal regions, it functions as a key part of food webs and the biogeochemical cycle in these marine and brackish ecosystems [9–11]. Owing to its rapid growth and ample productivity, with more than 3.3 million tons produced globally per year, Manila clam supports aquaculture industries in China, Korea and Japan, and is among the top-25 aquatic species introduced in European countries like Spain, France, Italy and Ireland [12,13]. However, worsening hypoxia in shallow coastal waters and estuaries increasingly challenges the survival of Manila clam and mariculture industry development. In China, the severity of summer hypoxia has been reportedly rising in coastal areas of Shandong Peninsula, an important natural habitat and mariculture region for Manila clam. Dissolved oxygen (DO) in Laizhou Bay and Rushan Bay dropped below 2.0 mg L^-1, the situation was worse in the Xiaoqing River estuary, with DO concentrations of less than 0.5 mg L^-1 [14–16]. In Japan, mass mortalities of Manila clam caused by hypoxia have occurred many times during summer, bringing tremendous economic loss to the fishery [17,18].

The tolerance to hypoxia by marine benthos varies among taxa. Bivalves, such as Manila clam, are believed to be more tolerant than fishes and crustaceans. The mean LC₅₀ for several common bivalves was estimated to be 1.42 ± 0.14 mg L^-1, which is much lower than that estimated for some crustaceans (2.45 ± 0.14 mg L^-1) [4]. In addition, different species of bivalves have different levels of tolerance; the LC₅₀ for Atlantic surf clam Spisula solidissima was measured as 0.5 mg L^-1, while for Baltic clam Macoma balthica it was 1.7 mg L^-1 [19,20]. Previous studies have indicated that Manila clams are able to survive more than a week in severe oxygen-deficient conditions [17,18]. However, most experiments with this species have been conducted at a single extremely low DO level, often < 0.5 mg L^-1, and there are no quantitative results for LC₅₀ or LT₅₀.

Animals have adopted a variety of strategies to survive hypoxia, from behavioral to metabolic. Typically, a buried bivalve will first ensure its oxygen supply by emerging from the sediment, to then enhance its respiratory rate and the oxygen-carrying capability of the blood. But if the hypoxia lasts for a relatively long time, a metabolic response may be crucial [21,22]. Farrer’s scallop Chlamys farreri and blue mussel Mytilus edulis, whose distribution overlaps with that of Manila clam, were found to have depressed their oxygen consumption rate (OCR), from 0.003 μg g^-1 h^-1 and 0.14 mg g^-1 h^-1 to 0.0005 μg g^-1 h^-1 and 0.05 mg g^-1 h^-1, respectively, when DO levels dropped from 6.0 mg L^-1 to less than 2.0 mg L^-1 [23]. The activities of enzymes, the driving force of metabolism, also react to hypoxic stress. In fishes, crustaceans and mollusks, antioxidant enzymes like superoxide dismutase and catalase are widely reported as being affected by hypoxia. As a natural response to protect an organism from physiological damage, these enzymes are usually immediately activated when DO level starts to drop; yet under long-term hypoxia their activities have a greater chance of becoming depressed, exposing the benthos to other potential damage via environmental stress [24–26]. However, respiratory enzyme activities, which are closely related to energy supply, have received little attention, especially in bivalves. Moreover, recent studies of Manila clam as well as other marine benthos have focused on functional responses, while structural responses at the cellular level are still unclear.

In the present research, long-term laboratory experiments were carried out to determine the effects of hypoxia on the survival, behavior, metabolism and cellular damage of R. philippinarum. A novel hypoxia simulation device, designed following the principle of hysteresis and with a few other innovations, was adopted for the first time, enabling us to precisely reduce the DO increment in the experimental design to as little as 0.2 mg L^-1.

Materials and methods

Experimental animals and sediments

Manila clams used for the experiments were collected from Xidayang Cun Farm in Jiaozhou Bay (36°11'15"N, 120°16' 34"E) and the sampling was carried out by fishing boats using sucking pumps. The sampling was approved by the Xidayang Cun Council. All clams were transported to the laboratory within three hours. During transportation, the clams were held in storage boxes with seawater and with compressed air continuously pumped in. Before the experiments, all clams were brushed and kept in aquariums without sediment. Two-thirds of the water (salinity 30.0 ± 1.0) was replaced every two days, and Chlorella vulgaris (9.13 ± 2.48 mg L^-1) was provided as food once a day. The water temperature was maintained at 22°C using heaters. Compressed air was pumped in for 30 minutes every 8 hours, and the DO concentration was closely monitored to ensure that it did not drop below 6.0 mg L^-1. All aquariums were covered with 2-mm-thick black canvas to avoid light stimulus.

Natural sediments were collected from the intertidal zone at the same site where the clams were collected. During sampling, several surveying rods were inserted into the sediment to specify the depth so that only the top 12 ± 3 cm of sediment inhabited by Manila clams was collected. After transportation to the laboratory, the collected sediment was sieved through 2-mm mesh to remove other benthos and then placed in the experimental tanks to achieve a 12-cm sediment layer. The water above the sediment was 45 cm deep and there was no headspace.

Before being placed in the experimental tanks, the clams were carefully examined for their health and measured for weight and length. First, as a sign of health, the shells needed to be strictly intact and in a half-open state when there was no stimulus. Next, when the siphons were touched with a glass rod two times, the clam had to withdraw its siphons quickly. Among the clams satisfying these requirements, the ones with a similar weight and length were selected for the experiments (Table 1).

Download:

Table 1. Mean (± SD) shell length and wet weight of the Ruditapes philippinarum selected for the two hypoxia experiments.

https://doi.org/10.1371/journal.pone.0215158.t001

A novel hypoxia simulation device

To ensure more accurate and stable DO concentrations and to eliminate the disturbance introduced by bubbling in nitrogen or changing the tank water, and thereby to maintain good water quality during the long-term experiments, we developed an automated hypoxia simulation device which also featured a temperature-control function. The device consisted of four identical environmental simulation systems, a water-exchange system, and a central control system (Fig 1). For each environmental simulation system, a supply tank served as a chamber to allow adjustment of the DO concentration and water temperature; this was placed below four replicate experimental tanks. This design prevented the animals from being shocked by the environmental oscillations. The supply tank was equipped with both a chiller and a heater for temperature regulation, while DO adjustment was achieved by bubbling in nitrogen or compressed air under the control of solenoid valves. To measure DO and water temperature, sensors were installed inside the supply tank. Of these, the DO sensors were based on the principle of luminescence quenching, which offers a variety of advantages over traditional electro-chemical probes (e.g. no need for stirring; less potential for fouling; long-term stability without calibration; and faster response time) [27]. Once the system was functioning, seawater was pumped from the supply tank to the experimental tanks and then back again. Between them was a purifier to filter out the metabolic wastes produced by the clams. The environmental simulation system contained 1,160 L of seawater, with 800 L in the supply tank and 90 L in each experimental tank. The huge water capacity together with the recirculating arrangement ensured a fairly stable environment.

Download:

Fig 1. Design of the automated hypoxia simulation device.

(A) The sketch showed the main structure of the device. The structure labeled with ‘a’ was the environmental simulation system; the structure labeled with ‘b’ was the water-exchange system; the structure labeled with ‘c’ was the central control system. (B) A logic diagram indicated how the device was running, including the water flow and gas flow.

https://doi.org/10.1371/journal.pone.0215158.g001

The water-exchange system featured a 900-L reserve tank which had a similar structure as the supply tank. When the system was activated, the DO concentration and temperature of the seawater in the reserve tank could be adjusted to the same levels as in the target supply tank. After that, a certain volume of seawater would be discharged from the target supply tank, and then an equal amount of seawater would be supplemented under the control of the water solenoid valves.

The central control system comprised a DO meter, a relay and a computer, which automatically actuated the compressed air or nitrogen solenoid valves for adjusting the DO level and also controlled the chiller or heater.

Since all sensors require time to output readings, the DO data received was often several minutes old. However, this design would lead to an excessive input of nitrogen or oxygen once the desired DO concentration had been reached. Worse, it would also lead to an excessive input of the gas since it took time for oxygen to spread evenly in the water. To effectively solve this problem, we employed the principle of hysteresis in designing the central control system. The hysteresis principle involves a control for switching logic; this took the time-delay described above into consideration, based on the relationship between the current DO reading, the actual DO value, and the DO spreading time. Thus, the actuation and closing of the solenoid valves could be performed before the DO reading exceeded the control range. In addition, the water exchange could also proceed automatically by the actuation of a corresponding water solenoid valve, after instructions were relayed.

The water supply and reserve tanks were made of polypropylene, while the experimental tanks were made of acrylic. A DO meter (OX11875) and DO sensors (OX11250) were bought from Loligo Systems (Viborg, Denmark). The pumps (MP-30R) were produced by Wenzhou Xinxishan Co., Ltd (Wenzhou, China). The solenoid valves (Q22MD-L8) were produced by Dalian Fangda Pneumatic Component Co., Ltd (Dalian, China). The water chillers, heaters (HSLS1000R Series, 660 kW), and custom-made purifiers were provided by Dalian Huixin Titanium Equipment Development Co., Ltd (Dalian, China).

Experiment 1: Responses under moderate and severe hypoxia

For the first 20-day hypoxia tolerance experiment, four DO concentrations were set: 0.5, 1.0 and 2.0 mg L^-1 for the treatment groups, and 6.0 mg L^-1 for the control group, corresponding to the four environmental simulation systems of the device (salinity 30.0 ± 1.0; water temperature 22 ± 0.5°C; flow rate 250 L h^-1). Each clam selected from the holding aquarium was marked with an individual number on the shell before being placed into an experimental tank. Every 96 clams were randomly assigned to each of the four groups, with 24 clams assigned to each of four replicate tanks. Three tanks (tanks A, B and C) were intended for observations of survival and behavior, as well as the analyses of metabolism and cellular damage performed at the end of the experiment; one tank (tank D) was intended for metabolism research on the 10^th day of the experiment. The experimental tanks were covered with 2-mm-thick black canvas to block out sunlight. All clams could bury themselves and adapt to the environment for 96 hours. During the first 72 hours, we replaced any clams that died or did not bury, while during the last 24 hours we only removed the dead or unburied ones without making replacements. To check if a clam is dead, we touched the siphons or foot three times with a glass rod. If the clam showed no response, it was judged to be dead. All handling was completed at least 8 hours before the experiment commenced. The animals were not fed during the acclimation period and throughout the whole experiment.

The DO concentration was maintained at 6.0 mg L^-1 during the 96-hour adaptation period and thereafter reduced to the target (treatment) values within 48 hours. During the experiment, mortalities and the number of clams emerging from the sediment in tanks A, B and C were monitored every 8 hours, while DO was monitored daily using the Winkler method, and pH was measured at the beginning and end of the experiment. In addition, to maintain quality of the seawater, one-third was replaced every 3 days.

On the 10^th day of the experiment, 8 clams were chosen randomly from tank D for measurements of OCR and ammonia-N excretion rate (AER) as well as for determination of the O:N ratio, and another 8 clams were chosen for enzyme activities analysis.

At the end of the experiment, 19 clams were chosen randomly from tanks A, B and C. Of these, 8 clams were used for the measurement of OCR and AER as well as determination of the O:N ratio; 8 clams were used for enzyme activity analysis; and 3 clams were prepared for study under a transmission electron microscope (TEM). To ensure randomness of sampling and mitigate potential bias, a random-number software was used, and the clams were removed according to the number on their shell. To collect clams still buried in the sediment and to avoid excessive disturbance to them, a low-power pump was used to suck the sediment away, layer by layer. After the buried clams were exposed, they were carefully removed for the given analysis.

To measure the OCR, AER and O:N ratio, each chosen clam was immediately moved from the experimental tank to a 1-L sealed beaker of water. The beaker was warmed by a 22°C water bath and covered with 2-mm-thick black canvas. Seawater stored in the supply tank was used as the source water for the corresponding group of clams, except that the clams exposed to 0.5 mg L^-1 DO were analyzed using the source water from the 1.0 mg L^-1 DO supply tank, out of concern that a starting DO concentration of 0.5 mg L^-1 might be too low to control for error in the OCR measurement. The source water was sterilized by ultraviolet lamp for 24 hours before use. The concentrations of DO and ammonia-N were measured at the beginning and 1 hour later. The calculation formulas were as follows: where DO₀ and N₀, respectively, refer to the concentrations of DO and ammonia-N at the beginning of the experiment, and DO_t and N_t are the concentrations at the end; W is the dry weight of the clam, and T is the duration of the measurement; the O:N ratio reflects the composition of metabolic fuel [28]. An O:N ratio with a value of 7 indicates that energy will be provided by protein; the higher the value, the more the energy supply is derived from fats and carbohydrates instead of protein [29].

For measures of enzyme activity, the activities of lactate dehydrogenase (LDH), phosphofructokinase (PFK), and pyruvate kinase (PK) in adductor muscle were determined spectrophotometrically by monitoring the absorbance change at 340 nm caused by the oxidation of NADH. The adductor muscles were cut into roughly 3 × 3 × 3 mm sections and immediately stored in liquid nitrogen. Commercial assay kits (the total protein assay kit A045-4, lactate dehydrogenase assay kit A020-2, phosphofructokinase test kit A129 and pyruvate kinase assay kit A076-1; Nanjing Jiancheng Bioengineering Institute, Nanjing, China) were used according to the manufacturer’s instructions. The experimenter responsible for the measurements was blinded to the treatment origin of the clams as the centrifuge tubes were renumbered after sampling by another experimenter.

For observations with the TEM (JEM-1200EX; JEOL, Tokyo, Japan), foot muscles were sectioned into roughly 2 × 2 × 3 mm pieces and immediately fixed with 1 mL of 2.5% glutaraldehyde fixative before being stored in a refrigerator at 4°C. All samples were renumbered and sent to the Medical College of Qingdao University (Qingdao, China) for cell and organelle structure analysis within 2 hours; the experimenters were blinded from any information about the experiment design as well as the clam treatment groups.

Before Experiment 1, we had carried out another two pre-experiments on survival and behavioral response with the same experimental design as Experiment 1. One of the pre-experiments was terminated on the 15^th day for accidental electricity cut off.

Experiment 2: Survival rates within the key DO concentration range

To further explore the survival rate against DO concentration, especially the key DO concentration range that corresponded to a drastic change in the survival rate, another set of 20-day tolerance experiments was carried out. The DO concentrations were set based on the survival results of Experiment 1, and 6.0 mg L^-1 was adopted for all control groups. Each experiment involved 96 randomly chosen clams in each group, with 24 clams in each of four replicate tanks (Table 2). Other environmental conditions were kept the same as in Experiment 1 (salinity 30.0 ± 1.0; water temperature 22 ± 0.5°C; flow rate 250 L h^-1). All tanks were covered with 2-mm-thick black canvas. Before the start of the experiments, all clams were allowed to adapt for 96 hours. As before, we replaced any dead or unburied clams found during the first 72 hours, but only removed these without planting new ones during the following 24 hours. The animals were not fed during the adaption periods and the whole experiments. Mortalities were checked daily and DO was monitored by the Winkler method. One-third of the seawater was replaced every three days, as in Experiment 1.

Download:

Table 2. Dissolved oxygen (DO) concentrations used in Experiment 2.

https://doi.org/10.1371/journal.pone.0215158.t002

Statistical analyses

OCR, AER, O: N and all measures of enzyme activity were tested for statistical significance. To ensure preciseness of the conclusions, all data were first tested for normality and homogeneity of variances. If the results passed all two tests, analysis of variance (ANOVA) would be adopted, and a Tukey test was used for multiple comparisons. If the results failed to pass normality tests or homogeneity of variances tests, the Kruskal–Wallis H-test was adopted, and multiple comparisons were conducted with the Mann–Whitney U-test. The numbers of surviving and buried clams were analyzed using the latter method, as nonparametric statistical methods are more reliable (although less sensitive) than ANOVA when values for replicates are not large enough (Supporting information). No difference being introduced by hypoxia was given as the null hypothesis and a p-value of < 0.05 defined statistical significance of the results.

To better control the potential deviation caused by other factors, we made a comparison across Pre-experiment 1, Pre-experiment 2 and Experiment 1. The calibrated values (CV) of survival rates and burial rates for each group in each experiment was defined as follows: where n₁ and n₂, respectively, refer to the initial number of surviving clams in one specific group and the number of surviving clams on the 15^th day in the same group; N₁ and N₂, respectively, refer to the initial number of surviving clams in the corresponding control group and the number of surviving clams on the 15^th day of the same control group; n₃ and N₃, respectively, refer to the number of buried clams in the specific group and the corresponding control group mentioned above. The CV were then compared across different DO concentrations, using Kruskal–Wallis H-test.

The 20-day LC₅₀ for 0.5 mg L^-1 DO and the LT₅₀ were determined using probit analysis [30]. All statistical analyses were carried out with SPSS 16.0 (IBM, Armonk, New York State, United States of America).

Results

Survival and behavioral response

The experimental DO concentrations were accurately and stably controlled by the automated hypoxia simulation device during both experiments, which ensured reliability of the results (Fig 2). The measures of DO concentration levels in all the experiments were attached to the Supporting information.

Download:

Fig 2. Dissolved oxygen (DO) treatments used in Experiments 1 and 2; the control group in each experiment was maintained at 6.0 mg L^-1 DO.

(A) Experiment 1 consisted of three treatment groups, set at 0.5, 1.0 and 2.0 mg L^-1 DO, and a control group. (B) Experiment 2 comprised four independent experiments with different sets of treatments: Experiment 2.1 –treatments 0.40, 0.60 and 0.80 mg L^-1 DO, and a control group; Experiment 2.2 –treatments 0.45, 0.65 and 0.85 mg L^-1 DO, and a control group; Experiment 2.3 –treatments 0.50, 0.70 and 0.90 mg L^-1, and a control group; Experiment 2.4 –treatments 0.55, 0.75 and 0.95 mg L^-1, and a control group.

https://doi.org/10.1371/journal.pone.0215158.g002

Hypoxia poses a great threat to the survival of Manila clam. However, mortalities only occurred under the low DO concentrations of 0.5 and 1.0 mg L^-1, while a moderately hypoxic environment of 2.0 mg L^-1 did not cause any clams to die (Fig 3). During the first 7 days of exposure, at all DO concentrations, the clams displayed strong tolerance to hypoxia as no mortalities were recorded. During the next 3 to 5 days, a few of the clams exposed to 0.5 mg L^-1 DO began to die at a slow rate, followed by deaths among those exposed to 1.0 mg L^-1 DO. The survival rate decreased sharply among clams exposed to the most severe hypoxia (0.5 mg L^-1 DO), resulting in much lower survival among these than the clams exposed to 1.0 mg L^-1 DO (p = 0.05). The LT₅₀ at 0.5 mg L^-1 DO was calculated as 422 hours, and the 20-day LC₅₀ was estimated to be 0.57 mg L^-1 DO when the experiments using more DO concentrations were conducted as part of Experiment 2. When the DO concentration was decreased from 1.0 to 0.4 mg L^-1, survival dropped from nearly 90% to approximately 20%; survival began to drop more suddenly around 0.8 mg L^-1 DO and dropped most dramatically at about 0.6 mg L^-1 DO (Fig 4).

Download:

Fig 3. Changes in the number of surviving and buried Ruditapes philippinarum in Experiment 1.

(A) Treatment group exposed to 0.5 mg L^-1 DO. The dashed lines denote the estimated critical-time thresholds dividing the curve into different phases; the first dashed line marks the beginning of the mortalities, and the second cuts the rest of the curve into two parts that could each be linearly fitted with best fitness (the largest sum of R values in Excel). (B) Treatment group exposed to 1.0 mg L^-1 DO. (C) Treatment group exposed to 2.0 mg L^-1 DO. (D) Control group maintained at 6.0 mg L^-1 DO. N = 72 per treatment.

https://doi.org/10.1371/journal.pone.0215158.g003

Download:

Fig 4. Survival rates of the Ruditapes philippinarum in Experiment 2.

Each circle represents one replicate in one specific experiment. The 20-day LC₅₀ for DO was estimated by probit analysis, based on the combined results of Experiments 2.1–2.4. N = 96 per treatment.

https://doi.org/10.1371/journal.pone.0215158.g004

As shown in Table 3, the variances between the results of Experiment 1 and the two Pre-experiments were small after calibration by control group, and the effects of hypoxia on survival and behavior were still significant, suggesting that there was little chance that Experiment 1 was affected by a potential factor.

Download:

Table 3. Calibrated values (CV) for survival rates and burial rates of the Ruditapes philippinarum (* indicates p < 0.05).

https://doi.org/10.1371/journal.pone.0215158.t003

As a major behavioral response towards hypoxia, the Manila clams emerged from the sediment in their tanks to obtain more oxygen; at the lower DO concentrations, the less the clams buried (Fig 3). However, the clams that were successively exposed to all DO concentrations seemed to reach a relatively stable condition around the 9^th day, even after dramatic changes in burial rates. All clams exposed to 0.5 mg L^-1 DO emerged at the surface of the sediment, which suggests this level of hypoxia is an extremely challenging environment for the species.

Metabolic response

OCR, AER, and the O:N ratio were significantly influenced by hypoxia. However, these values indicated a response only to certain DO concentrations, especially the lower concentrations (Fig 5). By the end of the experiment, no difference in OCR was detected between the clams exposed to 1.0, 2.0 and 6.0 mg L^-1 DO (p = 0.878), whereas the clams exposed to 0.5 mg L^-1 DO exhibited distinct respiratory depression (p = 0.022). For AER, a significant suppression effect under hypoxia could be observed when the DO concentration dropped to 1.0 mg L^-1 and below. The clams exposed to 2.0 and 6.0 mg L^-1 DO did not differ significantly in AER (p = 0.85). The O:N ratio rose sharply when the DO dropped from 2.0 to 1.0 mg L^-1 and below. There was no significant difference between the results for clams removed on the 10^th day and at the end of the experiment, whether for values of OCR, AER or O:N (p > 0.05).

Download:

Fig 5. Metabolic response to hypoxia by Ruditapes philippinarum on the 10^th day and at the end of the experiment, under different DO concentrations.

(A) Oxygen consumption rate on the tenth day; (B) ammonia-N excretion rate on the tenth day; (C) O:N ratio on the tenth day. (D) oxygen consumption rate at the end of experiment; (E) ammonia-N excretion rate at the end of experiment; (F) O:N ratio at the end of experiment. Error bars indicate mean ± SE (** indicates p < 0.01, * indicates p < 0.05).

https://doi.org/10.1371/journal.pone.0215158.g005

The effects of hypoxia on LDH, PFK, and PK activities in adductor muscle were significant. Compared with the response mode of PK and PFK, whose activities peaked under the most severe hypoxia (0.5 mg L^-1 DO); LDH reached its highest activity at the DO concentration of 2 mg L^-1 (Fig 6). The depression of LDH activity when the DO dropped from 2.0 to 1.0 mg L^-1 may suggest diversified pathways of pyruvate. The results on the 10^th day did not differ from those at the end of the experiment (p > 0.05), except that LDH activity did not peak at the DO concentration of 2.0 mg L^-1.

Download:

Fig 6. Biochemical responses to hypoxia measured in adductor muscle of Ruditapes philippinarum on the 10^th day and at the end of the experiment, under different DO concentrations.

(A) Lactate dehydrogenase activity on the tenth day; (B) phosphofructokinase activity on the tenth day; (C) pyruvate kinase activity on the tenth day. (D) lactate dehydrogenase activity at the end of experiment; (E) phosphofructokinase activity at the end of experiment; (F) pyruvate kinase activity at the end of experiment. Error bars indicate mean ± SE (** indicates p < 0.01, * indicates p < 0.05).

https://doi.org/10.1371/journal.pone.0215158.g006

Cellular damage

Under hypoxic stress, damage was observed at the cellular level in the Manila clams (Fig 7). Mitochondria, the main site of aerobic respiration for cells, were damaged to various degrees, with collapsed cristae, shriveled membranes and induced cell inclusion. Though the cell remained intact with no membrane rupture, vacuolization was apparent from the lighter staining of chromatin. Other organelles seemed to present a normal appearance except that the myofilaments tended to dissolve.

Download:

Fig 7. Transmission electron microscope images of the foot muscle of Ruditapes philippinarum.

(A), (B), (C) and (D) represent the DO groups of 6.0, 2.0, 1.0 and 0.5 mg L^-1, respectively. White arrows point to normal mitochondria while black arrows point to mitochondria with collapsed cristae (judged by sparse cristae) and vacuolization (judged by shallower staining). Cell vacuolization was apparent according to the shallower cell-staining. Some other cellular structures for point of reference are listed as follows: ‘a’ refer to glycogen granules; ‘b’ refer to myofilaments; ‘c’ refer to cell membrane.

https://doi.org/10.1371/journal.pone.0215158.g007