Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Intraspecific variation in oxidative stress tolerance in a model cnidarian: Differences in peroxide sensitivity between and within populations of Nematostella vectensis

  • Lauren E. Friedman,

    Roles Formal analysis, Investigation, Methodology

    Current address: Harvard T. H. Chan School of Public Health, Boston, MA

    Affiliation Department of Biology, Boston University, Boston, MA, United States of America

  • Thomas D. Gilmore,

    Roles Funding acquisition, Supervision, Writing – review & editing

    Affiliation Department of Biology, Boston University, Boston, MA, United States of America

  • John R. Finnerty

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Resources, Supervision, Writing – original draft, Writing – review & editing

    jrf3@bu.edu

    Affiliation Department of Biology, Boston University, Boston, MA, United States of America

Abstract

Nematostella vectensis is a member of the phylum Cnidaria, a lineage that includes anemones, corals, hydras, and jellyfishes. This estuarine anemone is an excellent model system for investigating the evolution of stress tolerance because it is easy to collect in its natural habitat and to culture in the laboratory, and it has a sequenced genome. Additionally, there is evidence of local adaptation to environmental stress in different N. vectensis populations, and abundant protein-coding polymorphisms have been identified, including polymorphisms in proteins that are implicated in stress responses. N. vectensis can tolerate a wide range of environmental parameters, and has recently been shown to have substantial intraspecific variation in temperature preference. We investigated whether different clonal lines of anemones also exhibit differential tolerance to oxidative stress. N. vectensis populations are continually exposed to reactive oxygen species (ROS) generated during cellular metabolism and by other environmental factors. Fifteen clonal lines of N. vectensis collected from four different estuaries were exposed to hydrogen peroxide. Pronounced differences in survival and regeneration were apparent between clonal lines collected from Meadowlands, NJ, Baruch, SC, and Kingsport, NS, as well as among 12 clonal lines collected from a single Cape Cod marsh. To our knowledge, this is the first example of intraspecific variability in oxidative stress resistance in cnidarians or in any marine animal. As oxidative stress often accompanies heat stress in marine organisms, resistance to oxidative stress could strongly influence survival in warming oceans. For example, while elevated temperatures trigger bleaching in corals, oxidative stress is thought to be the proximal trigger of bleaching at the cellular level.

Introduction

The starlet sea anemone, Nematostella vectensis, is a leading laboratory model organism in the phylum Cnidaria due to its ease of culture, amenability to a wide range of experimental protocols, and the availability of abundant genomic resources [119]. In its natural habitat—estuaries along the Atlantic coast of North America—this animal experiences dramatic variation in key environmental variables over short temporal and spatial scales [20]. Furthermore, there is extensive genetic and phenotypic variation between populations of N. vectensis [2127]. Atlantic coast populations have a high degree of genetic differences and regional population genetic structure that is not correlated to relative geographic distances [21, 25]. These differences within local populations are hypothesized to be the result of limited gene flow within and between N. vectensis populations, adaptation to local environmental conditions, and episodic anthropogenic dispersal over large distances. The estuarine habitats inhabited by N. vectensis are often physically isolated, and these habitats can differ dramatically with respect to key environmental parameters, such as temperature, pH, and salinity [3, 20]. Moreover, the early stages of the N. vectensis lifecycle exhibit limited dispersal ability, due in part to the egg mass being negatively buoyant [3]. The combination of diverse environments and limited dispersal will tend to promote adaptive genetic differentiation between populations. For this reason, N. vectensis has been advocated as a useful model for studying adaptive microevolution in a coastal invertebrate beset by rapid environmental change [6, 24].

One environmental stressor that varies dramatically in salt marshes and other estuarine habitats is hydrogen peroxide (H2O2) concentration. For example, in an intertidal salt-flat off the coast of Germany, H2O2 concentrations in the water were found to vary almost 100-fold, from 0.05 μM to approximately 4.5 μM, with concentrations at the sediment-water boundary layer being 2–4 times higher than concentrations in the adjacent bottom water [28]. Peroxide levels can accumulate in shallow coastal habitats due to natural inputs such as photochemical production and atmospheric wet deposition as well as anthropogenic inputs such as runoff from industrial processes, e.g., mining and industrial manufacturing. Salt marshes are particularly susceptible to peroxide accumulation because large amounts of organic matter can fuel peroxide production through photo-degradation [29], and isolated pools may seldom be flushed.

N. vectensis is a benthic salt-marsh organism with relatively low mobility, so this species is readily exposed to compounds that accumulate in sediments of estuaries. In its natural habitat, N. vectensis can be exposed to a wide range of exogenously produced reactive oxygen species (ROS) including H2O2. This anemone must, therefore, be under strong selection to counteract oxidative stress. N. vectensis is known to tolerate a wide range of other highly variable environmental parameters, e.g., salinities from ~0–50 ppt and temperatures from 2–28.5°C [3, 20]. The ability to tolerate rapid and dramatic variations in these environmental conditions suggests that individual N. vectensis have wide environmental tolerances. At the same time, the pronounced environmental differences between and within estuaries inhabited by N. vectensis suggest the potential for local adaptation to differing selection pressures. The potential for local adaptation is increased by the fact that N. vectensis’ natural dispersal ability appears limited, as evidenced by significant genetic structure between and within estuaries [21, 25]. Specifically with reference to peroxide tolerance, one might expect that all N. vectensis would exhibit a wide tolerance to peroxide, with no observable differences between or within populations. However, given that a recent study identified substantial differences in temperature-specific growth and regeneration rates in N. vectensis [24], we sought to determine whether different clonally derived lines of anemones would also exhibit differential peroxide tolerance. Indeed, as elevated temperatures promote oxidative stress (e.g., [30]), we hypothesized that the same clonal lines that are best adapted to high temperatures would also be most resistant to oxidative stress.

In this study, we exposed fifteen clonal lines of N. vectensis collected from four geographically isolated estuaries to different levels of H2O2. Our study included three clonal lines (one each from Baruch, South Carolina, Meadowlands, New Jersey, and Kingsport, Nova Scotia) known to vary in their growth and/or regeneration rates at 21 and 29°C [24], as well as 12 previously uncharacterized clonal lines established from founder individuals collected in Great Sippewissett Marsh, Falmouth, MA. We monitored the effects of H2O2 on adult survival and on the ability of animals to regenerate a new “head” following bisection through the body column. We identified substantial differences in the effects of peroxide on survival and regeneration between the SC, NJ, and NS clonal lines and among clonal lines originating from the single estuary in Sippewissett.

Materials and methods

Animal collection and husbandry

The clone lines used in this study were derived from individual field-caught anemones collected prior to 2008 for a population genetic study from estuaries in four locations: Kingsport, Nova Scotia; Great Sippewissett Marsh, Falmouth, MA; Meadowlands, NJ; and Baruch, SC [25]. Within Sippewissett Marsh, founder individuals were collected from three distinct pools designated Pool A, B and C. Shortly after capture, individual anemones were isolated to separate culture containers and then used to establish clonal lines by repeatedly bisecting the animals and allowing them to regenerate. All clonal lines were housed in the laboratory using standard culture conditions [1] for at least five years before being used in this study. The clone lines used in this study were previously tested for their temperature-specific growth and regeneration rates [24].

Preliminary peroxide exposure tests

The levels of peroxide used in our experiments were determined empirically. In a pilot experiment conducted by J. C. Sullivan in the Finnerty lab [31], multiple strains of N. vectensis, including the Nova Scotia, New Jersey and South Carolina clone lines tested in this study, were challenged to regenerate under a range of H2O2 concentrations from 0.0001–0.01% (32.6–3260 μM). At concentrations exceeding 0.0005% (163 μM) H2O2, anemone mortality reached 100%. At 0.0005% H2O2, some anemones survived the exposure, albeit while exhibiting a pronounced delay in the regeneration rate relative to control conditions. Thus, for the current study, we characterized adult survivorship at 0.0005% H2O2, while we evaluated regeneration rate at half that concentration (0.00025% or 82 μM).

Peroxide survival experiments

Several individual anemones from each clone line (n = 5–10) were maintained in 0.0005% (163 μM) H2O2/11 ppt artificial seawater (ASW; Instant Ocean™) or under control conditions (11 ppt artificial seawater) for a period of two weeks. The ASW/H2O2 was made fresh daily, and both ASW/H2O2 and ASW media were changed daily. Over the course of these experiments, each individual anemone was housed in a single well of a 12-well tissue-culture plate (Thermo Scientific Nunc) and maintained at 20°C on a 12 h-light/12 h-dark diel cycle. Once per day, individual anemones were scored for tentacle number and for the ability to contract the tentacles and body column in response to mechanical stimulation. Single factor ANOVAs were performed to identify significant differences in average tentacle number between peroxide-exposed and control individuals on each day of the study (p<0.05). Animals that had lost the ability to respond to mechanical stimulation were scored as deceased. All surviving anemones that had been subjected to peroxide exposure were retired from future stress-tolerance studies due to concerns over potential lingering epigenetic effects of exposure.

Peroxide regeneration experiments

Individual anemones from each clonal line (n = 10–20) were bisected along the transverse axis at the midpoint of the body column, approximately equidistant from the tip of the foot and the opening of the mouth [1]. The aboral halves of 50% of the bisected anemones (n = 5–10) from each clonal line were then transferred to 0.00025% (82 μM) H2O2/11ppt ASW (i.e., half of the H2O2 concentration used in the survival experiments). These animals were then allowed to regenerate for 14 days, while the other 50% were allowed to regenerate for the same period under control conditions (11 ppt ASW). The H2O2/ASW or ASW media were changed daily. The anemones were sorted between H2O2 and control treatments so that mean body length was approximately equal between treatments. Individuals used in these experiments measured 0.5–3.0 cm in length from the tip of the foot (or physa) to the oral opening. Regenerating anemones were scored daily for the presence and number of tentacles, and after 14 days, they were scored for their ability to feed. Previous studies have shown that the regeneration process is slowed when the bisection is not clean or the mesenteries are exposed [32]; therefore, if the bisection site of the regenerating anemones was not a clean cut, then the individual was not included in the study. ANOVAs were performed to determine significant differences in average tentacle number during regeneration between H2O2-exposed and control individuals. As with the H2O2 survival experiments, all surviving anemones that successfully regenerated when exposed to H2O2 were retired from any future stress-tolerance studies due to concerns over potential lingering epigenetic effects of exposure.

Results

Differences between genets in adult survivorship upon exposure to H2O2

To determine whether anemones varied in their tolerance to H2O2, over a period of 14 days, we compared the number of tentacles in adult polyps incubated in various concentrations of H2O2 to those in anemones in control conditions (Fig 1). We used this metric as an indicator of the anemones’ condition because tentacle number is easy to score, and the tentacles are particularly sensitive to tissue damage from elevated H2O2. Furthermore, it is necessary for the tentacles to be extended during feeding, and tentacles are typically extended in unstressed animals. In all clone lines tested, there were no changes in mean tentacle numbers in control conditions. However, relative to controls, there was a significant decrease in tentacle number in several clonal lines under conditions of H2O2 exposure (Figs 2 and 3). In the clonal line from Nova Scotia (NS), mean tentacle number remained generally steady over the first ten days of peroxide exposure, but declined on days 11–12, and remained lower than in control anemones, though the difference was not statistically significant. In the clonal line from New Jersey (NJ), H2O2-exposed and control animals began to diverge in tentacle number at day 12, and the difference between treatments became statistically significant by day 14. In the South Carolina clonal line, neither a decline in tentacle number nor any visible tissue disfigurement was observed over the course of the treatment period.

thumbnail
Fig 1. Representative adult tolerance to H2O2 exposure.

Representative individuals from Sippewissett MA (top row: clone line MA B1; bottom row: clone line MA A1) during exposure to 0.0005% H2O2 (A-D) or control conditions (E-F) over 7 days. After a week, the majority of the clonal lines showed severe disfigurement and loss of tentacles (as in B), but some clonal lines appeared relatively resistant to H2O2 exposure (as in D).

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

thumbnail
Fig 2. Tentacle number during 14-day exposure to H2O2 or control conditions.

Data for 15 N. vectensis clone lines are shown. In all graphs, days of exposure are given along the X-axis, and tentacle numbers are shown on the Y-axis. Data from the control condition are shown in black, while data from H2O2 exposure are shown in colored lines. The number of individuals tested under control conditions and H2O2 exposure is provided for each clonal line (n = #/#). Statistically significant differences between H2O2 and control conditions on a given day are indicated using asterisks (* = p<0.05; ** = p<0.005). Error bars depict the standard deviations for each mean.

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

thumbnail
Fig 3. Mean tentacle counts during H2O2 exposure.

For each of 15 N. vectensis clonal lines, mean tentacle number for 4–11 individuals is plotted against days of exposure to 0.0005% H2O2. Thin vertical lines depict standard errors.

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

Differences between the NJ, NS, and SC clonal lines in their response to elevated H2O2 were not surprising given previously demonstrated differences in their response to elevated temperatures [24]. However, we observed more substantial differences in H2O2 tolerance among 12 N. vectensis clonal lines generated from individuals collected from three separate pools in Sippewissett Marsh, MA in 2008 (Figs 2 and 3). In particular, three Sippewissett clonal lines appeared minimally affected by H2O2 exposure (MA-A1, MA-C2, and MA-C4), while the other nine clonal lines displayed varying levels of disfigurement and loss of tentacles over time (Figs 2 and 3). The three clonal lines minimally affected by H2O2 exposure had been collected from two different pools in Sippewissett Marsh. Two of the resistant clonal lines came from Pool C, which also contained three clonal lines that were sensitive to H2O2 exposure. The other resistant clonal line came from Pool A (MA-A1), and there were four sensitive clonal lines also from that pool (MA-A2, A3, A4 and A5). Unfortunately, we are unable to directly compare the performance of MA-A1 under H2O2-exposed versus control conditions because this clonal line perished, apparently due to disease, before the control could be performed.

Regeneration is differentially affected by H2O2 in different genets

To determine whether H2O2 had an effect on regeneration, we bisected animals, isolated the aboral halves, and then counted the number of tentacles in these regenerating animals once per day over a 14-day period under control conditions or in the presence of H2O2 (Fig 4). In some clonal lines, the rate of regeneration was significantly affected by H2O2 exposure (Fig 5). The NS clonal line showed a significant effect of H2O2 exposure on average number of tentacles present on days 5–12 (p<0.05), with highly significant effects on days 5, 11, and 12 (p<0.005). The H2O2-exposed NJ individuals showed significantly lower average tentacle numbers on day 7 as well as days 11–14 of regeneration (p<0.05), with highly significant differences apparent by day 14 (p<0.005). The SC clone line showed no statistically significant effects of H2O2 exposure on tentacle number throughout the 14-day regeneration period.

thumbnail
Fig 4. Representative regenerating N. vectensis during H2O2 exposure.

An adult N. vectensis from New Jersey, approximately 1 cm in length (A), was bisected through the body column (at the site of the arrow) using a scalpel. The cut divided the anemone into (B) an oral fragment bearing the mouth and tentacles (tn) and an aboral fragment bearing the foot, or physa (phy). The scalpel scored the bottom of the plastic petri dish in which the bisection was performed (arrow). Following bisection, the regenerating aboral fragment was exposed to 0.00025% H2O2 over a period of 14 days. Seven days after bisection (C), the aboral fragment had developed 10 tentacles. The tentacle tips are numbered in the photo. By day 14 (D), 12 tentacles had developed. In panel D, the pharynx (pha) and several of the eight complete mesenteries (mes) are visible.

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

thumbnail
Fig 5. Tentacle counts in regenerating N. vectensis.

Regeneration ability as measured by mean tentacle number over time (n = 5–11) during exposure to 0.00025% H2O2 or control conditions. In all graphs, the number of days since the anemones were bisected is indicated along the X-axis, and tentacle number is shown on the Y-axis. Data from the control condition are shown in black, while data from H2O2 exposure are shown in colored lines. The number of individuals tested under control conditions and H2O2 exposure is provided for each clonal line (n = #/#). Statistically significant differences between H2O2 and control conditions on a given day are indicated using asterisks (* = p<0.05; ** = p<0.005).

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

Continual H2O2 exposure had varying effects on regeneration in 11 clonal lines from Sippewissett Marsh (Fig 5). In general, the ability to regenerate in 0.00025% H2O2 was correlated with the ability to survive in 0.0005% H2O2. From MA pool A, the four clonal lines that were sensitive to H2O2 in the adult survival assay (MA-A2, A3, A4, A5) also showed delayed regeneration in the presence of H2O2 at some time points. In clonal lines MA-A4 and A5, a significant detrimental effect of H2O2 was observed until the end of the study, becoming increasingly severe over the 14-day period. In MA-A2 and A3, the tentacle counts of individual anemones undergoing H2O2 treatment were similar to the mean tentacle numbers of the control specimens by the end of the study. Unfortunately, the clonal line from pool A that exhibited the greatest tolerance to H2O2 in the adult survival assay (MA-A1) was lost due to an apparent disease before it could be subjected to the H2O2 /regeneration assay. While neither of the two clonal lines from pool B was able to tolerate 0.0005% H2O2 as an adult, both were able to regenerate without any detectable effect in 0.00025% H2O2. From pool C, in the two clonal lines that demonstrated the greatest H2O2 tolerance as adults (MA-C4 and C2), there were also negligible or minor effects of H2O2 on regeneration. In MA-C2, regeneration was not affected by H2O2, whereas in MA-C4, regeneration was affected by H2O2 on days 3 and 4, but H2O2-treated individuals were indistinguishable from control anemones on days 5–14. Some of the less resistant clones (MA-C1 and C3) showed delayed rates of tentacle development early in regeneration (e.g., day 4), but were similar to control animals by the end of the study. By contrast, clonal line MA C5 was the most sensitive to H2O2, as it showed the most persistent and severe detrimental effects of H2O2 on regeneration, with a highly significant (p<0.005) decrease in average tentacle number versus control animals from days 8–14.

In a comparison of the performance of all clonal lines in the presence of H2O2 (Fig 6), the greatest variation in tentacle number was observed on days 4–6. For example, on day 4, the clonal line with the lowest mean tentacle count (MA-C1) averaged fewer than two tentacles per individual anemone, whereas the clonal line with the highest mean tentacle count (MA-C4) averaged approximately ten tentacles per individual. By day 14, many clonal lines appeared to achieve a steady-state number of tentacles ranging between ~9 (MA-B1) and ~14 (NJ).

thumbnail
Fig 6. Mean tentacle counts during regeneration under H2O2 exposure.

For each of 14 N. vectensis clonal lines, mean tentacle number for 5–11 individuals (see Fig 6) is plotted against days of exposure to 0.00025% H2O2. Thin vertical lines depict standard deviations.

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

Discussion

In this study, we applied exogenous H2O2 to assess the oxidative stress response of the estuarine sea anemone Nematostella vectensis. This is a straightforward and biologically relevant way to challenge the oxidative stress response of a coastal marine invertebrate because organisms in shallow coastal waters regularly encounter elevated H2O2 levels that can trigger oxidative stress (e.g., [28]). Furthermore, peroxide is the most abundant and long-lived ROS in marine environments, and unlike other ROS, H2O2 is uncharged and can diffuse freely across membranes (reviewed in [33]).

The peroxide levels used in our experiments (82 and 163 μM) are substantially higher than those recorded in marine and freshwater environments (Table 1). Across numerous ocean basins, the peroxide concentration of sea surface waters have been found to range from 10–420 nM [3436], while peroxide concentrations in the surface waters of lakes have been found to exceed 1 μM (e.g., [37]). At a site more comparable to the habitat of N. vectensis, peroxide levels in a shallow tidal pool located on a sand flat, ~500 m off the coast of Germany were found to vary between 0.5 and 4.5 μM [28]. However, in contrast to the high marsh habitat in which the anemones used in the current study were collected, this sand flat is inundated by the tides for ~19 h each day [28]. Of note, the levels of H2O2 in rainwater can be substantially higher than the surface waters of oceans or lakes, with continental rainwater ranging from 0.1–247 μM (reviewed in [38, 39]) and oceanic rainwater ranging from 3.5 to 82 μM [40, 41]. As a result, precipitation events can substantially increase H2O2 levels in surface waters [4042].

thumbnail
Table 1. Hydrogen peroxide levels measured in aquatic environments.

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

While H2O2 levels have not yet been reported for any saltmarsh, we suspect that the pools inhabited by N. vectensis regularly experience higher concentrations than those typically found in oceans or lakes. The saltmarsh pools where the anemones in this study were collected are small, shallow bodies of water. As such, they receive proportionately large inputs of H2O2 through precipitation. Snow and rain can both have much higher peroxide concentrations than seawater (Table 1), and in shallow, stratified layers of marine surface waters, rain events have been shown to increase the H2O2 concentration by more than 10-fold [40]. Furthermore, saltmarsh pools receive abundant direct sunlight and are rich in organic matter. As a result, the production of H2O2 via the photochemical decay of organic matter, which is thought to be the principal source of H2O2 in surface waters [29, 43], should be particularly high. Furthermore, saltmarsh pools, particularly those at higher elevation, may only rarely be inundated by the tides, so their H2O2 content would not often be diluted by an influx of seawater. Finally, in the upper layer of the sediment, where the anemones burrow, photosynthetic algae and bacteria can generate locally high levels of oxygen [44] and concomitantly high levels of ROS.

The effects of H2O2 exposure have been studied in a handful of marine invertebrates, including polychaete worms, molluscs, crustaceans, and corals, but the concentrations of peroxide that were used and the effects that were monitored have varied substantially among these studies (Table 2). In whole-organism studies, micromolar concentrations of H2O2 (0.5–20 μM) were shown to decrease oxygen consumption and aerobic metabolism in three intertidal species, a burrowing polychaete [45], a clam [46], and a limpet [47]. Five micromolar H2O2 also caused a substantial drop in water filtration rate in an intertidal clam [46]. A higher level of H2O2 (50 μM) was used in experiments on the estuarine polychaete, Laeonereis acuta. In this infaunal worm, peroxide exposure caused lipid peroxidation and DNA damage [48]; there was also a biphasic effect on aerobic metabolism, with O2 consumption falling after 4 d exposure and rising again after 7 d [49]. In the tropical reef coral, Galaxea fascicularis, exposure to 0.3 or 3 μM H2O2 was accompanied by an increase in catalase activity, but no increase in superoxide dismutase activity [50], and this pattern was observed at both 27°C and 31°C [51]. In comparison to our study, all of these organismal investigations utilized lower levels of H2O2, and they did not examine tissue integrity, mortality, or a developmental process (such as regeneration). Significantly higher levels of peroxide have been applied to an isolated body-wall preparation obtained from an intertidal polychaete (42–664 μM; [52]) and to a mitochondrial suspension from a sipunculid (90.9 μM; [53]). Finally, millimolar levels of H2O2 have been shown to induce synchronous spawning in molluscs in a study that did not investigate subsequent impacts of H2O2 exposure on survivorship or physiological function [54]. In studies on cultured cells, H2O2 levels comparable to or greater than those employed here have been used to induce “mild” oxidative stress; for example, Huang et al. exposed HeLa cells to 1000 μM H2O2 for 1 h or 100 μM H2O2 for 6 hours to investigate its effects on protein localization [55]. Of note, our study demonstrates the highest level of survivable peroxide exposure that has been found to date for any marine or estuarine invertebrate.

thumbnail
Table 2. Effects of hydrogen peroxide exposure on marine invertebrates.

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

While previous studies on oxidative stress tolerance in marine animals have generally sought to characterize a species’ average or typical response to a given challenge, our study is the first to examine intraspecific variation in the oxidative stress response of a marine animal. Indeed, given their differing thermal preferences, we suspected that these clonal lines might differ in peroxide tolerance, because thermal stress is known to induce oxidative stress and because ROS such as H2O2 can work synergistically with acute temperature changes to cause oxidative stress. For example, in the Antarctic limpet Nacella concina, H2O2 exposure and elevated temperatures contribute additively to oxidative stress [47]. Similarly, in corals thermal stress has been shown to result in oxidative DNA damage (e.g., [56]) and to induce antioxidant enzyme activity [51]. Although our study did not directly determine whether there is a correlation between sensitivity to thermal stress and to oxidative stress, the SC clonal line, which was found to be most tolerant of high temperature (29°C; [24]), was among the least affected by H2O2, and the NS clonal line, which was most sensitive to high temperature, was more affected by H2O2. However, the Sippewisset clone lines varied extensively in H2O2 tolerance, and the previous temperature study did not sample enough individuals from this single estuary to determine whether thermal tolerance varies correspondingly within this population.

The most surprising result of our study is the significant differences that we identified in survival and regeneration during H2O2 exposure among clonal lines taken from a single marsh in Sippewissett, MA. The Sippewissett clonal lines used in our experiments were derived from wild-caught individuals that lived within a hundred meters of each other, some collected from the same small pool where they would have presumably encountered similar environmental stressors; i.e., each of the genets collected from the same pool would have been subjected to the same daily and seasonal variation in temperatures. Given these similar temperature exposures, one would not expect that animals from a single pool would differ widely in their temperature preferences. Nevertheless, the differences among clonal lines from this single pool in their ability to tolerate H2O2 exposure actually exceeded the variation exhibited between the SC and NS clonal lines. The clonal lines from Sippewisset were collected at the same time and maintained in the laboratory for an equivalent length of time prior to the conduct of this study. Therefore, the differences we observed cannot be explained by differing lengths of time under laboratory culture.

There are multiple possible explanations for the variation we observed within and between pools at Sippewissett. For example, there may be adaptation or long-term acclimation to different microhabitats. Spatially- and temporally-variable factors that can affect peroxide levels include sediment composition and shifting sediment distribution, the amount of organic matter, water depth and tidal action, the amount of pollution and heavy metal contamination, primary production, flooding events (and the associated interconnectedness of tidal pools), aerial exposure, periods of anoxia, solar radiation, and other microbial plant and algal species present [33]. Comparable fine-scale variation in other environmental variables has been shown to affect organismal physiology. For example, substantial environmental micro-variation in temperature has been shown to exist in tide pool habitats [57], and mussels originating only 24 m apart were found to differ substantially in their response to thermal variation [58]. Similarly, changes in vertical gradients of oxygen, H2S, temperature, and pH have been shown to effect changes in molecular response to ROS in the burrowing marine polychaete Heteromastus filiformis [59]. Gradients of pH have also been shown to affect both oxygen consumption and H2O2 levels. Previous exposure to sublethal stressors has also been shown to increase stress resistance in a number of model systems an in response to a number of stressors (reviewed in [60]).

Relatively few studies have investigated intraspecific variation in the response to H2O2, and the bulk of these studies have been performed on plants and algae. For example, intertidal individuals of the green macroalga Ulva lactuca exhibited a greater ability to scavenge exogenous H2O2 than subtidal individuals of the same species [61]. To our knowledge, the only experimental study documenting intraspecific variation in H2O2 tolerance in an animal model was conducted on the soil nematode Caenorhabditis remanei [62]. In this study, multiple strains of C. remanei were exposed to 3.5 mM H2O2, an oxidative challenge that kills all individuals within 4-8h of exposure. The offspring of different sires were found to differ in their ability to withstand this level of H2O2 exposure, and the heritability of this trait was found to be 2x greater than that for heat stress resistance and 11x greater than that for lifespan [62]. This study used a qualitatively different challenge than used in the current study, where the anemones were tested at a level of H2O2 exposure they could survive for two weeks.

Intraspecific variability in oxidative stress response has been studied in a number of marine invertebrates, including corals and sea anemones. ROS-levels and antioxidant enzyme activities exhibit seasonal variation in the coral Pocillopora capitata [63]. In a population of the sea anemone Anemonia viridis spanning a natural pCO2 gradient, levels of the algal secondary metabolite dimethylsulfoniopropionate (DMSP) as well as the activity of the enzyme superoxide dismutase (SOD) were found to be lower where CO2 concentrations were higher [64]. The authors suggest that the symbiotic algae inhabiting the sea anemone’s tissues experience greater photosynthetic efficiency at higher CO2 levels, and they therefore generate fewer ROS. In both of these naturalistic studies, the role of ROS could not be isolated from other correlated environmental factors. Nor could these studies discriminate environmental versus genetic causes of the observed phenotypic changes, because it was not known whether the individuals experiencing different environments were part of the same genet or different genets.

Moving forward, an important goal will be to understand the genetic mechanisms (or stably inherited epigenetic mechanisms) that explain the differences in peroxide tolerance observed here. Importantly, according to population genetic studies using AFLPs, the Sippewisset population is genetically distinct from any of the other twenty-three populations that have been studied, including the Baruch, Kingsport, and New Jersey populations [21, 24]. All 36 Sippewisset individuals characterized in these studies appear far more similar to each other than to individuals from any other population [25]. Thus, the most highly resistant clone line from Sippewisset, MA-C2, is much more similar genetically to the least resistant clones line from Sippewissett, MA-A2, than it is to the highly resistant SC clone line. This suggests that differences in oxidative stress tolerance could have evolved independently, and perhaps with a different genetic mechanism, in different populations.

The current study illustrates the importance of incorporating organismal assays into studies of oxidative stress. That an organism is experiencing oxidative stress is often inferred from elevated levels of antioxidant enzymes or non-enzymatic antioxidant defense mechanisms. However, resistant individuals may exhibit greater stress tolerance precisely because they express stress-response genes at higher levels, either constitutively or conditionally (i.e., in response to the relevant stressor) [65]. What ultimately determines the fate of an organism under stress is its ability to maintain physiological function in the face of stressful environmental conditions.

If the degree of fine-scale variability in oxidative stress tolerance observed here in N. vectensis is widespread among natural populations, it could have profound importance for the resilience of species and ecological communities in this era of rapidly increasing global temperatures because thermal stress is so often coupled with oxidative stress. This is particularly true for reef-building corals. The loss of photosynthetic endosymbiotic zooxanthellae (“bleaching”) is a leading cause of coral mortality. Bleaching is most often associated with elevated temperatures, but the proximate mechanism underlying bleaching is thought to be oxidative stress [66, 67]. As a result, those corals that prove most resistant to thermal bleaching may be those with the most robust anti-oxidant defenses. We suggest that future studies should assess the natural variability in oxidative stress resistance that exists among corals (and other relatively heat-sensitive taxa), and determine if the molecular mechanisms responsible for oxidative stress resistance also confer greater resistance to thermal stress. Such knowledge could prove valuable in developing predictive models for the fate of natural populations in the face of climate change.

Acknowledgments

The research performed here was based upon preliminary experiments conducted by J. C. Sullivan.

References

  1. 1. Stefanik DJ, Friedman LE, Finnerty JR. Collecting, rearing, spawning and inducing regeneration of the starlet sea anemone, Nematostella vectensis. Nat Protoc. 2013;8(5):916–23. pmid:23579780.
  2. 2. Stefanik DJ, Wolenski FS, Friedman LE, Gilmore TD, Finnerty JR. Isolation of DNA, RNA and protein from the starlet sea anemone Nematostella vectensis. Nat Protoc. 2013;8(5):892–9. pmid:23579778.
  3. 3. Hand C, Uhlinger KR. The culture, sexual and asexual reproduction, and growth of the sea anemone Nematostella vectensis. Biol Bull. 1992;182(2):169–76. PubMed PMID: WOS:A1992KC97400005. pmid:29303672
  4. 4. Ambrosone A, Marchesano V, Mazzarella V, Tortiglione C. Nanotoxicology using the sea anemone Nematostella vectensis: from developmental toxicity to genotoxicology. Nanotoxicology. 2014;8(5):508–20. pmid:23641943.
  5. 5. Botman D, Kaandorp JA. Spatial gene expression quantification: a tool for analysis of in situ hybridizations in sea anemone Nematostella vectensis. BMC Res Notes. 2012;5:555. pmid:23039089; PubMed Central PMCID: PMC3532226.
  6. 6. Darling JA, Reitzel AR, Burton PM, Mazza ME, Ryan JF, Sullivan JC, et al. Rising starlet: the starlet sea anemone, Nematostella vectensis. Bioessays. 2005;27(2):211–21. pmid:15666346.
  7. 7. DuBuc TQ, Dattoli AA, Babonis LS, Salinas-Saavedra M, Rottinger E, Martindale MQ, et al. In vivo imaging of Nematostella vectensis embryogenesis and late development using fluorescent probes. BMC Cell Biol. 2014;15:44. pmid:25433655; PubMed Central PMCID: PMC4264334.
  8. 8. Fischer AH, Tulin S, Fredman D, Smith J. Employing BAC-reporter constructs in the sea anemone Nematostella vectensis. Integr Comp Biol. 2013;53(5):832–46. pmid:23956207.
  9. 9. Fritzenwanker JH, Technau U. Induction of gametogenesis in the basal cnidarian Nematostella vectensis (Anthozoa). Dev Genes Evol. 2002;212(2):99–103. pmid:11914942.
  10. 10. Genikhovich G, Technau U. Bromodeoxyuridine labeling of S-phase nuclei in the starlet sea anemone Nematostella vectensis. Cold Spring Harb Protoc. 2009;2009(9):pdb prot5284. pmid:20147269.
  11. 11. Genikhovich G, Technau U. Anti-acetylated tubulin antibody staining and phalloidin staining in the starlet sea anemone Nematostella vectensis. Cold Spring Harb Protoc. 2009;2009(9):pdb prot5283. pmid:20147268.
  12. 12. Genikhovich G, Technau U. In situ hybridization of starlet sea anemone (Nematostella vectensis) embryos, larvae, and polyps. Cold Spring Harb Protoc. 2009;2009(9):pdb prot5282. pmid:20147267.
  13. 13. Genikhovich G, Technau U. Induction of spawning in the starlet sea anemone Nematostella vectensis, in vitro fertilization of gametes, and dejellying of zygotes. Cold Spring Harb Protoc. 2009;2009(9):pdb prot5281. pmid:20147266.
  14. 14. Genikhovich G, Technau U. The starlet sea anemone Nematostella vectensis: an anthozoan model organism for studies in comparative genomics and functional evolutionary developmental biology. Cold Spring Harb Protoc. 2009;2009(9):pdb emo129. pmid:20147257.
  15. 15. Layden MJ, Rottinger E, Wolenski FS, Gilmore TD, Martindale MQ. Microinjection of mRNA or morpholinos for reverse genetic analysis in the starlet sea anemone, Nematostella vectensis. Nat Protoc. 2013;8(5):924–34. pmid:23579781.
  16. 16. Putnam NH, Srivastava M, Hellsten U, Dirks B, Chapman J, Salamov A, et al. Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization. Science. 2007;317(5834):86–94. pmid:17615350.
  17. 17. Sullivan JC, Ryan JF, Watson JA, Webb J, Mullikin JC, Rokhsar D, et al. StellaBase: the Nematostella vectensis Genomics Database. Nucleic Acids Res. 2006;34(Database issue):D495–9. pmid:16381919; PubMed Central PMCID: PMC1347383.
  18. 18. Tulin S, Aguiar D, Istrail S, Smith J. A quantitative reference transcriptome for Nematostella vectensis early embryonic development: a pipeline for de novo assembly in emerging model systems. Evodevo. 2013;4:16. pmid:23731568; PubMed Central PMCID: PMC3748831.
  19. 19. Wolenski FS, Layden MJ, Martindale MQ, Gilmore TD, Finnerty JR. Characterizing the spatiotemporal expression of RNAs and proteins in the starlet sea anemone, Nematostella vectensis. Nat Protoc. 2013;8(5):900–15. pmid:23579779.
  20. 20. Hand C, Uhlinger KR. The unique, widely distributed, estuarine sea anemone, Nematostella vectensis Stephenson—a Review, new facts, and questions. Estuaries. 1994;17(2):501–8. PubMed PMID: WOS:A1994PK14200015.
  21. 21. Darling JA, Reitzel AM, Finnerty JR. Regional population structure of a widely introduced estuarine invertebrate: Nematostella vectensis Stephenson in New England. Mol Ecol. 2004;13(10):2969–81. pmid:15367113.
  22. 22. Sullivan JC, Wolenski FS, Reitzel AM, French CE, Traylor-Knowles N, Gilmore TD, et al. Two alleles of NF-κB in the sea anemone Nematostella vectensis are widely dispersed in nature and encode proteins with distinct activities. PLoS One. 2009;4(10):e7311. pmid:19806194; PubMed Central PMCID: PMC2751831.
  23. 23. Darling JA, Reitzel AM, Finnerty JR. Characterization of microsatellite loci in the widely introduced estuarine anemone Nematostella vectensis. Mol Ecol Notes. 2006;6(3):803–5. PubMed PMID: WOS:000239778600067.
  24. 24. Reitzel AM, Chu T, Edquist S, Genovese C, Church C, Tarrant AM, et al. Physiological and developmental responses to temperature by the sea anemone Nematostella vectensis. Mar Ecol Prog Ser. 2013;484:115–30. PubMed PMID: WOS:000320321300008.
  25. 25. Reitzel AM, Darling JA, Sullivan JC, Finnerty JR. Global population genetic structure of the starlet anemone Nematostella vectensis: multiple introductions and implications for conservation policy. Biol Invasions. 2008;10(8):1197–213. PubMed PMID: WOS:000260306700002.
  26. 26. Reitzel AM, Herrera S, Layden MJ, Martindale MQ, Shank TM. Going where traditional markers have not gone before: utility of and promise for RAD sequencing in marine invertebrate phylogeography and population genomics. Mol Ecol. 2013;22(11):2953–70. pmid:23473066; PubMed Central PMCID: PMC3669247.
  27. 27. Reitzel AM, Sullivan JC, Finnerty JR. Discovering SNPs in protein coding regions with StellaSNP: Illustrating the characterization and geographic distribution of polymorphisms in the estuarine anemone Nematostella vectensis. Estuar Coast. 2010;33(4):930–43. PubMed PMID: WOS:000278468300011.
  28. 28. Abele-Oeschger D, Tug H, Rottgers R. Dynamics of UV-driven hydrogen peroxide formation on an intertidal sandflat. Limnol Oceanogr. 1997;42(6):1406–15. PubMed PMID: WOS:000071910200010.
  29. 29. Cooper WJ, Zika RG. Photochemical formation of hydrogen peroxide in surface and ground waters exposed to sunlight. Science. 1983;220(4598):711–2. pmid:17813875.
  30. 30. Madeira D, Narciso L, Cabral HN, Vinagre C, Diniz MS. Influence of temperature in thermal and oxidative stress responses in estuarine fish. Comp Biochem Physiol A. 2013;166(2):237–43. pmid:23774589.
  31. 31. Sullivan JC. Bioinformatic and molecular analsyses of the Nematostella genome reveal a mixture of ancient conservation and rapid recent evolution. Boston, MA: Boston University; 2008.
  32. 32. Reitzel AM, Burton PM, Krone C, Finnerty JR. Comparison of developmental trajectories in the starlet sea anemone Nematostella vectensis: embryogenesis, regeneration, and two forms of asexual fission. Invertebr Biol. 2007;126(2):99–112. PubMed PMID: WOS:000248912400001.
  33. 33. Lesser MP. Oxidative stress in marine environments: biochemistry and physiological ecology. Annu Rev Physiol. 2006;68:253–78. pmid:16460273.
  34. 34. Yuan JC, Shiller AM. Distribution of hydrogen peroxide in the northwest Pacific Ocean. Geochem Geophy Geosy. 2005;6. doi: Artn Q09m02 PubMed PMID: WOS:000232109700001.
  35. 35. Yuan JC, Shiller AM. The distribution of hydrogen peroxide in the Southern and central Atlantic Ocean. Deep-Sea Res Pt Ii. 2001;48(13):2947–70. PubMed PMID: WOS:000169483800011.
  36. 36. Moore CA, Farmer CT, Zika RG. Influence of the Orinoco River on hydrogen peroxide distribution and production in the eastern Caribbean. J Geophys Res-Oceans. 1993;98(C2):2289–98. PubMed PMID: WOS:A1993KM44300007.
  37. 37. Hakkinen PJ, Anesio AM, Graneli W. Hydrogen peroxide distribution, production, and decay in boreal lakes. Can J Fish Aquat Sci. 2004;61(8):1520–7. PubMed PMID: WOS:000224976800017.
  38. 38. Sakugawa H, Kaplan IR, Tsai WT, Cohen Y. Atmospheric hydrogen peroxide. Environ Sci Technol. 1990;24(10):1452–62. PubMed PMID: WOS:A1990EA64600005.
  39. 39. Vione D, Maurino V, Minero C, Pelizzetti E. The atmospheric chemistry of hydrogen peroxide: A review. Ann Chim-Rome. 2003;93(4):477–88. PubMed PMID: WOS:000183659900020.
  40. 40. Cooper WJ, Saltzman ES, Zika RG. The contribution of rainwater to variability in surface hydrogen peroxide. J Geophys Res-Oceans. 1987;92(C3):2970–80. PubMed PMID: WOS:A1987G535100014.
  41. 41. Yuan JC, Shiller AM. The variation of hydrogen peroxide in rainwater over the South and Central Atlantic Ocean. Atmos Environ. 2000;34(23):3973–80. PubMed PMID: WOS:000088749000012.
  42. 42. Miller WL, Kester DR. Peroxide variations in the Sargasso Sea. Mar Chem. 1994;48(1):17–29. PubMed PMID: WOS:A1994QB63400002.
  43. 43. Zika RG, Moffett JW, Petasne RG, Cooper WJ, Saltzman ES. Spatial and temporal variations of hydrogen peroxide in Gulf of Mexico waters. Geochim Cosmochim Acta. 1985;49:1173–84.
  44. 44. De Wit R, Jonkers HM, Van Den Ende FP, Van Gemerden H. In situ fluctuations of oxygen and sulfide in marine microbial sediment ecosystems. Neth J Sea Res. 1989;23(3):271–81.
  45. 45. Abele-Oeschger D, Oeschger R, Theede H. Biochemical adaptations of Nereis diversicolor (Polychaeta) to temporarily increased hydrogen peroxide levels in intertidal sandflats. Mar Ecol Prog Ser. 1994;106(1–2):101–10. PubMed PMID: WOS:A1994NB76600009.
  46. 46. Abele-Oeschger D, Buchner T. Effect of environmental hydrogen peroxide accumulation on filtration rates of the intertidal bivalve clam Cerastoderma edule Verh Dtsch Zool Ges. 1995;88:93.
  47. 47. Abele D, Burlando B, Viarengo A, Portner HO. Exposure to elevated temperatures and hydrogen peroxide elicits oxidative stress and antioxidant response in the Antarctic intertidal limpet Nacella concinna. Comp Biochem Phys B. 1998;120(2):425–35. PubMed PMID: WOS:000075742600022.
  48. 48. da Rosa CE, Bianchini A, Monserrat JM. Antioxidant responses of Laeonereis acuta (Polychaeta) after exposure to hydrogen peroxide. Braz J Med Biol Res. 2008;41(2):117–21. PubMed PMID: WOS:000253337700007. pmid:18297192
  49. 49. da Rosa CE, Iurman MG, Abreu PC, Geracitano LA, Monserrat JM. Antioxidant mechanisms of the Nereidid Laeonereis acuta (Anelida: Polychaeta) to cope with environmental hydrogen peroxide. Physiol Biochem Zool. 2005;78(4):641–9. pmid:15957118.
  50. 50. Higuchi T, Fujimura H, Arakaki T, Oomori T. Activities of antioxidant enzymes (SOD and CAT) in the coral Galaxea fascicularis against increased hydrogen peroxide concentrations in seawater. Proc 11th Intl Coral Reef Symp. 2009:926–30.
  51. 51. Higuchi T, Suzuki Y, Fujimura H. Multiple effects of hydrogen peroxide and temperature on antioxidants and bleaching. Proc 12th Intl Coral Reef Symp. 2012:4C_3.
  52. 52. Storch D, Abele D, Portner HO. The effect of hydrogen peroxide on isolated body wall of the lugworm Arenicola marina (L.) at different extracellular pH levels. Comp Biochem Physiol C. 2001;128(3):391–9. pmid:11255112.
  53. 53. Buchner T, Abele D, Portner HO. Oxyconformity in the intertidal worm Sipunculus nudus: the mitochondrial background and energetic consequences. Comp Biochem Physiol B. 2001;129(1):109–20. pmid:11337254.
  54. 54. Morse DE, Duncan H, Hooker N, Morse A. Hydrogen peroxide induces spawning in mollusks, with activation of prostaglandin endoperoxide synthetase. Science. 1977;196(4287):298–300. pmid:403609.
  55. 55. Huang C, Han Y, Wang Y, Sun X, Yan S, Yeh ET, et al. SENP3 is responsible for HIF-1 transactivation under mild oxidative stress via p300 de-SUMOylation. EMBO J. 2009;28(18):2748–62. pmid:19680224; PubMed Central PMCID: PMCPMC2750016.
  56. 56. Nesa B, Hidaka M, editors. Thermal stress increases oxidative DNA damage in coral cell aggregates. 11th Int Coral Reef Symp; 2008; Ft. Lauderdale, Florida.
  57. 57. Dowd WW, King FA, Denny MW. Thermal variation, thermal extremes and the physiological performance of individuals. J Exp Biol. 2015;218(Pt 12):1956–67. pmid:26085672.
  58. 58. Jimenez AG, Jayawardene S, Alves S, Dallmer J, Dowd WW. Micro-scale environmental variation amplifies physiological variation among individual mussels. Proc Biol Sci. 2015;282(1820):20152273. pmid:26645201; PubMed Central PMCID: PMCPMC4685784.
  59. 59. Abele D, Grosspietsch H, Portner HO. Temporal fluctuations and spatial gradients of environmental P-O2, temperature, H2O2 and H2S in its intertidal habitat trigger enzymatic antioxidant protection in the capitellid worm Heteromastus filiformis. Mar Ecol Prog Ser. 1998;163:179–91. PubMed PMID: WOS:000073032800017.
  60. 60. Freire CA, Welker AF, Storey JM, Storey KB, Hermes-Lima M. Oxidative stress in estuarine and intertidal environments (temperate and tropical). In: Abele D, Vázquez-Medina JP, Zenteno-Savin T, editors. Oxidative Stress in Aquatic Ecosystems: Blackwell Publishing Ltd.; 2011. p. 41–57.
  61. 61. Ross C, Van Alstyne KL. Intraspecific variation in stress-induced hydrogen peroxide scavenging by the ulvoid macroalga Ulva lactuca. J Phycol. 2007;43(3):466–74. PubMed PMID: WOS:000246686100006.
  62. 62. Reynolds RM, Phillips PC. Natural variation for lifespan and stress response in the nematode Caenorhabditis remanei. PLoS One. 2013;8(4):e58212. pmid:23658604; PubMed Central PMCID: PMC3637273.
  63. 63. Linan-Cabello MA, Flores-Ramirez LA, Zenteno-Savin T, Olguin-Monroy NO, Sosa-Avalos R, Patino-Barragan M, et al. Seasonal changes of antioxidant and oxidative parameters in the coral Pocillopora capitata on the Pacific coast of Mexico. Mar Ecol-Evol Persp. 2010;31(3):407–17. PubMed PMID: WOS:000280671300002.
  64. 64. Borell EM, Steinke M, Horwitz R, Fine M. Increasing pCO(2) correlates with low concentrations of intracellular dimethylsulfoniopropionate in the sea anemone Anemonia viridis. Ecol Evol. 2014;4(4):441–9. PubMed PMID: WOS:000331354500011. pmid:24634728
  65. 65. Koenigstein S, Pohlmann K, Held C, Abele D. Ecological comparison of cellular stress responses among populations—normalizing RT-qPCR values to investigate differential environmental adaptations. Bmc Ecol. 2013;13. doi: Unsp 21 PubMed PMID: WOS:000319456000001.
  66. 66. Downs CA, Fauth JE, Halas JC, Dustan P, Bemiss J, Woodley CM. Oxidative stress and seasonal coral bleaching. Free Radic Biol Med. 2002;33(4):533–43. pmid:12160935.
  67. 67. Lesser MP. Oxidative stress causes coral bleaching during exposure to elevated temperatures. Coral Reefs. 1997;16(3):187–92. PubMed PMID: WOS:A1997XN81300008.
  68. 68. Abele D, Ferreyra GA, Schloss I. H2O2 accumulation from photochemical production and atmospheric wet deposition in Antarctic coastal and off-shore waters of Potter Cove, King George Island, South Shetland Islands. Antarct Sci. 1999;11(2):131–9. PubMed PMID: WOS:000084414100002.
  69. 69. Cooper WJ, Lean DRS. Hydrogen peroxide concentration in a northern lake—Photochemical formation and diel variability. Environ Sci Technol. 1989;23(11):1425–8. PubMed PMID: WOS:A1989AX23600021.
  70. 70. Abele-Oeschger D, Sartoris FJ, Portner HO. Hydrogen peroxide causes a decrease in aerobic metabolic rate and in intracellular pH in the shrimp Crangon crangon. Comp Biochem Phys C. 1997;117(2):123–9. PubMed PMID: WOS:A1997XC74300002.