ABSTRACT
Animals that inhabit environments that fluctuate in oxygen must not only contend with disruptions to aerobic metabolism, but also the potential effects of reactive oxygen species (ROS). The goal of this study was to compare aspects of ROS metabolism in response to O2 variability (6 h hypoxia or hyperoxia, with subsequent normoxic recovery) in two species of intertidal sculpin fishes (Cottidae, Actinopterygii) that can experience O2 fluctuations in their natural environment and differ in whole-animal hypoxia tolerance. To assess ROS metabolism, we measured the ratio of glutathione to glutathione disulfide as an indicator of tissue redox environment, MitoP/MitoB ratio to assess in vivo mitochondrial ROS generation, thiobarbituric acid reactive substances (TBARS) for lipid peroxidation, and total oxidative scavenging capacity (TOSC) in the liver, brain and gill. In the brain, the more hypoxia-tolerant Oligocottus maculosus showed large increases in TBARS levels following hypoxia and hyperoxia exposure that were generally not associated with large changes in mitochondrial H2O2. In contrast, the less-tolerant Scorpaenichthys marmoratus showed no significant changes in TBARS or mitochondrial H2O2 in the brain. More moderate increases were observed in the liver and gill of O. maculosus exposed to hypoxia and hyperoxia with normoxic recovery, whereas S. marmoratus had a greater response to O2 variability in these tissues compared with the brain. Our results show a species- and tissue-specific relationship between hypoxia tolerance and ROS metabolism.
INTRODUCTION
The majority of reactive oxygen species (ROS) are generated at sites along the mitochondrial electron transport system (ETS) and the rate of ROS generation is dependent upon the interactive effects of mitochondrial activity, proton motive force, ETS reduction/oxidation (redox) state and O2 availability (Quinlan et al., 2013). Mitochondria continuously produce ROS at low rates for cell signaling (D'Autréaux and Toledano, 2007), but ROS generation can increase when ETS flux is low, redox environment in mitochondria is highly reduced and proton motive force is high (Brand, 2000; Niknahad et al., 1995). ROS generation also increases when O2 availability is high, e.g. during normoxic recovery after hypoxia or during hyperoxia, when the redox environment in mitochondria is relatively oxidized (Aon et al., 2010; Yusa et al., 1987). High rates of mitochondrial ROS production can overwhelm total cellular scavenging capacities, resulting in oxidative damage to DNA, proteins and lipids (Cadet, 2003; Gutteridge, 1995; Reznick and Packer, 1994). In most vertebrates, particularly mammals, high rates of ROS production and the associated tissue damage typically occurs during periods of O2 fluctuation such as those associated with ischemia-reperfusion injury (Chouchani et al., 2014).
Animals inhabiting environments that naturally fluctuate in O2, particularly those that experience hypoxia, often display adaptations that enhance whole-animal O2 uptake and use under O2-limiting conditions, to sustain aerobic metabolism (Scott et al., 2009; Storz et al., 2009). In contrast, much less is known about whether animals that inhabit these environments have also evolved strategies to minimize ROS production or mitigate the damaging effects of ROS. Several studies support the idea that hypoxia tolerance is associated with lowered mitochondrial ROS generation (shown in elasmobranchs, mollusks and hypoxia-acclimated killifish Fundulus heterolitus; Du et al., 2016; Hickey et al., 2012; Ivanina and Sokolova, 2016), whereas there does not appear to be a consistent link between hypoxia tolerance and antioxidant defenses in fish (Leveelahti et al., 2014). However, a previous in vitro study of brain mitochondria from multiple species of intertidal sculpin (Cottidae, Actinopterygii) showed that mitochondria from more hypoxia-tolerant species emitted more H2O2 per O2 molecule consumed under oligomycin-induced state IV conditions compared with a less hypoxia-tolerant species, and H2O2 emitted per mg mitochondrial protein measured did not differ (Lau and Richards, 2019). Together, these in vitro studies suggest that at the level of the mitochondria, the relationship between hypoxia tolerance and ROS metabolism is nuanced, possibly due to variation in the role of ROS in cell signaling and the role of mitochondria as regulators of ROS metabolism in vivo (Munro and Treberg, 2017). Thus, the goal of this study was to determine if in vivo ROS metabolism differs between two species of sculpin, Oligocottus maculosus Girard 1856 and Scorpaenichthys marmoratus (Ayres 1854), which vary in hypoxia tolerance and experience different O2 regimes in their natural environment. O. maculosus resides predominately in rock pools in the marine intertidal zone, where O2 can fluctuate daily from near anoxia to hyperoxia (Richards, 2011), whereas S. marmoratus primarily inhabits the subtidal environment, which is near normoxia and far more O2 stable. Previous work has demonstrated that O. maculosus is far more hypoxia tolerant than S. marmoratus and showed they differ in brain mitochondrial function, but not tissue mitochondrial content (Lau and Richards, 2018; Lau et al., 2017; Mandic et al., 2009a).
In order to manipulate and gain insight into how ROS metabolism differs between hypoxia-tolerant and -intolerant fish, each species of sculpins was exposed to hypoxia (3.5 kPa O2) or hyperoxia (64 kPa O2) for 6 h, with or without 1 h recovery in normoxic water. These exposure times and PO2 levels mimic those observed in the marine intertidal zone (Richards, 2011). Following O2 manipulations, in vivo ROS metabolism was assessed in liver, gill and brain. Tissue redox status was assessed as the relative ratio of glutathione (GSH) and glutathione disulfide (GSSH), which is one of the dominant redox couples within the cell (Rahman et al., 2006). To estimate in vivo mitochondrial H2O2 production, we used the ratiometric mitochondria-targeted mass spectrometry probe MitoB, which is phenylboronic acid conjugated to a triphenylphosphonium ion (TPP+; Cochemé et al., 2011; Logan et al., 2014). MitoB can be injected into an animal where it accumulates in the mitochondrial matrix and reacts with H2O2 to form a stable phenol product, MitoP, which can then be assessed in tissue samples via mass spectrometry. To gain insight into the effects of our O2 manipulations on oxidative damage, we assessed lipid peroxidation by measuring thiobarbuturic acid reactive substance (TBARS) levels (Gutteridge, 1995). We also assessed tissue total oxidative scavenging capacity (TOSC). We hypothesized that the more hypoxia-tolerant sculpin O. maculosus would be more resistant to changes in environmental O2 and show less disturbance of in vivo ROS metabolism compared with the less hypoxia-tolerant sculpin S. marmoratus.
MATERIALS AND METHODS
Chemicals
MitoB compounds used for this study were synthesized by Prof. Richard Hartley (University of Glasgow). The details of the synthesis of MitoB and associated compounds are described in Cochemé et al. (2011). Chemicals for biochemical analysis were purchased from Sigma-Aldrich unless otherwise specified.
Animals
Adult O. maculosus (3.95±0.17 g; mean±s.e.m.) and S. marmoratus (21.71±1.64 g) were collected near Bamfield Marine Sciences Centre (British Columbia, Canada) at Ross Islets (48°52.4′N, 125°9.7′W) and Wizard's Rock (48°51.5′N, 125°9.4′W) using either handheld nets or pole seines during the lowest point in the tidal cycle. Animals were transported to the University of British Columbia (UBC) and housed in a recirculating aquaculture system (RAS) with artificial seawater at 12°C and maintained on a diet of shrimp, Atlantic krill and bloodworms for at least 3 weeks before experimentation. All experimental procedures were reviewed and approved by the UBC Animal Care Committee under animal use protocol number A13-0309.
Validation of MitoB use in sculpins
Individual O. maculosus were retrieved from stock RAS tanks and injected intraperitoneally with 30 nmol MitoB (50 µl injection volume in phosphate-buffered saline containing a small amount of food coloring to facilitate visibility) using a BD Ultra-Fine™ syringe (6 mm needle). Following injection, fish were transferred to separate 1.1 liter plastic mesh baskets (with three fish each) and held in a wet table with recirculating seawater at 12°C. At 0.5, 1, 4, 8, 24 and 72 h post-injection, three animals were removed from the mesh baskets, anesthetized with 0.5 g l−1 benzocaine, and gills, brain, muscle and liver were dissected, frozen in liquid nitrogen, and stored at −80°C until processing (see below).
To further validate whether MitoB behaves similarly in marine teleosts as in murine models (Cochemé et al., 2011), we used an ion-selective electrode for TPP+ to monitor the uptake of MitoB into isolated liver mitochondria (see Lau et al., 2017 for isolation protocol) from another species of marine sculpin, Artedius lateralis (Girard 1854) (Fig. S1). After MitoB titrations, the addition of NADH-generating substrates for complex I (5 mmol l−1 pyruvate, 10 mmol l−1 malate and 10 mmol l−1 glutamate) and 10 mmol l−1 succinate for complex II substrate to stimulate state II respiration rate led to a decrease in MitoB detected in the extra-mitochondrial surroundings, indicating an influx of MitoB into the matrix (Fig. S1). The addition of 0.5 µmol l−1 of the uncoupler carbonyl cyanide 4-(trifluoromethoxy)phenylhydrazone (FCCP) caused an increase in levels of MitoB detected in the extra-mitochondrial surroundings, indicating an efflux of MitoB from the matrix.
Experimental protocol and sampling
Both species of sculpin were exposed to one of five treatments after which tissue sampling occurred: (1) 6 h normoxia, (2) 6 h hypoxia, (3) 6 h hypoxia followed by 1 h normoxic recovery, (4) 6 h hyperoxia and (5) 6 h hyperoxia followed by 1 h normoxic recovery. For the hypoxia treatment, both species were exposed to 3.5 kPa, which is at the critical O2 tension for O2 consumption rate (Pcrit) for more hypoxia-tolerant O. maculosus and is equivalent to exposure to 65% of Pcrit in the less hypoxia-tolerant S. marmoratus (Mandic et al., 2009b). Thus, the hypoxia exposure was more severe relative to Pcrit for S. marmoratus than for O. maculosus. For hyperoxia treatment, both species were exposed to 64 kPa, which is a PO2 that is frequently experienced in the higher intertidal and less in the subtidal region (Richards, 2011).
The treatments were set up in two covered aquaria, one for normoxic control animals and the other for O2 manipulations (hypoxia and hyperoxia exposures were carried out on two separate days), which were placed in a wet table with recirculating water chilled to 12°C. Small circulating water pumps inside each aquarium ensured adequate mixing throughout the experiment. To initiate the experiment, individuals of both species were taken from their stock tank and injected intraperitoneally with MitoB (approximately 30 nmol per O. maculosus and 60 nmol per S. marmoratus in 50 µl phosphate-buffered saline, yielding ∼7.6 and 2.8 nmol g−1 in the two species, respectively which is within the range used in a previous study on brown trout; Salin et al., 2017) and placed into 1.1 liter plastic mesh baskets weighted with gravel (4–6 O. maculosus and 2–3 S. marmoratus in each basket) in the treatment aquaria. Although there were slight differences in the amount of MitoB injected in the two species, these differences had no impact on the study outcomes because both MitoB and MitoP were well above the detection limit of the mass spectrometer and the ratiometric nature of MitoP/MitoB analysis is insensitive to minor differences in the amount of mitoB injected. One hour after MitoB injection, fish for the first normoxic control time point were sampled as described below. The PO2 of the seawater was then adjusted to the desired level with either nitrogen gas (for hypoxia), 100% O2 (for hyperoxia), or aerated with compressed air (normoxia) and maintained at these levels for 6 h. Oxygen levels were monitored with a hand-held O2 probe (Oakton Instruments; for hypoxia) or a FOXY fluorescent O2 probe (Ocean Optics; for hyperoxia) and O2 level never deviated from the set point by more than 1.5%. After 6 h of hypoxia/hyperoxia treatment, a total of 12 individuals of O. maculosus and 8 S. marmoratus were sampled as described below. Following sampling at 6 h, both the hypoxic and hyperoxic aquaria were quickly returned to normoxia (within 30 min) by aeration and a total of 12 individuals of O. maculosus and 8 S. marmoratus were sampled at 1 h recovery once normoxia was established. Paired normoxic controls of 12 individuals of O. maculosus and 8 S. marmoratus were also sampled at 6 h exposure and 1 h recovery. For tissue sampling, the mesh basket containing fish was gently moved into a submerged 3.7 l container and both the fish and container were removed from the treatment aquarium and a high dose of anesthetic (0.5 g l−1 benzocaine) was introduced to the container. This procedure ensured that the animals were never exposed to air and that water PO2 was maintained at the desired level during anesthesia. We elected to use a high dose of benzocaine to facilitate rapid anesthesia (within 30 s). No signs of distress were noted during fish transfer or anesthesia. Once the fish were unresponsive to touch, they were removed from the anesthetic bath and brain, gill and liver were dissected, frozen in liquid N2, and stored at −80°C until analyses of MitoB, MitoP, TBARS, GSH/GSSG and TOSC. Owing to the small amount of tissue available, we did not perform all analysis on all tissues sampled from individuals, but instead prioritized analysis of MitoP/MitoB, TBARS, TOSC and performed analysis of GSH/GSSG only on available samples. As a result, some of the n values for GSH/GSSG are low, but still sufficient for basic statistical analysis.
Extraction and purification of MitoP/MitoB
Approximately 30 mg of frozen tissue was homogenized in 210 µl ice-cold 60% acetonitrile/0.1% formic acid spiked with internal standards (10 µmol l−1 d15-MitoB/5 µmol l−1 d15-MitoP) in a bullet bead homogenizer (Next Advance, Troy, NY, USA). Homogenates were centrifuged at 16,000 g for 10 min at 4°C, after which the supernatant was transferred to a new 1.5 ml microcentrifuge tube. The tissue was then resuspended in 200 µl 60% acetonitrile/0.1% formic acid, homogenized again and centrifuged at 16,000 g for 10 min at 4°C. The resulting supernatant was combined with the previous supernatant, vortexed for 10 s, and incubated at 4°C for 30 min. Samples were then centrifuged at 16,000 g for 10 min, and the supernatant was filtered with the Millipore centrifugal filter plate (0.45 µm hydrophilic, low protein binding Durapore membrane), centrifuged at 3000 g for 10 min. The filtrate was collected and the samples were dried in a vacuum speed centrifuge (CentriVap Concentrator; Labconco, Kansas City, MO, USA). The dried sample was then resuspended in 250 µl 20% acetonitrile/0.1% formic acid, vortexed for 5 min and centrifuged at 16,000 g for 10 min. 200 µl of the suspension was used for LC-MS/MS analysis. Standard curves for MitoB (0 to 1000 pmol) and MitoP (0 to 1000 pmol) were prepared with salmon tissue. The standards were prepared for analysis with the same procedure as described above (standard curves in Fig. S2).
LC-MS/MS analysis for MitoP/MitoB
The LC-MS/MS analysis was carried out according to the general protocol described elsewhere (Cochemé et al., 2012) with a few modifications. For LC-MS/MS analyses, the mass spectrometer was connected in series to an I-class Aquity LC system (Waters). Samples were stored in an autosampler at 4°C and 2 µl of the sample was injected into a 15 µl flow-through needle and RP-UPLC at 40°C using an Acquity UPLC® BEH C18 column (1×50 mm, 1.7 µm; Waters) with a Waters UPLC filter (0.2 µm). MS buffers A (95% water/5% acetonitrile/0.1% formic acid) and B (90% acetonitrile/10% water/0.1% formic acid) were infused at 200 µl min−1 using the following gradients: 0–0.3 min, 5% B; 0.3–3.0 min, 5–100% B; 3–4 min, 100% B; 4.0–4.1 min, 100–5% B; 4.1–4.6 min, 5% B. Eluant was diverted to waste at 0–1 min and 4–4.6 min. The compounds were detected in multiple reactions monitoring in positive ion mode. With this method selected, precursor ions are fragmented to product ions, which are used for the precursor ion identification. For quantification, the following transitions were used: MitoB, 397>183; d15–MitoB, 412>191; MitoP, 369>183; d15–MitoP, 384>191. The peak area of MitoB, MitoP and internal standard of samples and standard curves were quantified using MassLynx 4.1 software. The amount of MitoB and MitoP in the samples was determined using the standard curves. The deuterated internal standards were used to normalize samples to account for variation in sample volume.
Tissue glutathione redox potential
Glutathione (GSH) and glutathione disulfide (GSSG) were assayed with the enzymatic recycling method described in Rahman et al. (2006). Briefly, approximately 25 mg of each tissue was homogenized in 180 µl buffer with 0.1 mol l−1 KH2PO4 and 5 mmol l−1 EDTA at pH 7.5 (KPE buffer) with 0.1% Triton X-100 and 0.6% sulfosalicylic acid. The homogenized sample was then centrifuged at 3000 g for 10 min at 4°C. The supernatant was then divided for assessment of total GSH and GSSG. To determine GSSG, 50 µl of the samples or GSSG standards were incubated for 1 h at room temperature with 1 µl vinylpyridine (1:10 v/v in KPE buffer) to derivatize GSH. After 1 h, 3 µl triethanolamine (1:6 in KPE buffer) was added to the samples/standards and they were incubated for 10 min at room temperature, which was followed by the addition of 3 µl of 1 mol l−1 HCl to neutralize the sample. The GSSG and total GSH samples were then assayed using the same protocol in which glutathione reductase (GR) converts GSSG into GSH. For the assay, the buffer was prepared with equal volumes of 1.7 mmol l−1 5,5′-dithio-bis(2-nitrobenzoic acid) (DTNB) and glutathione reductase (3.33 U ml−1 KPE), of which 125 µl was added to 20 µl of standard/sample. 60 µl of 0.8 mmol l−1 β-NADPH was added to start the reaction and the rate of TNB formation was monitored at 412 nm for 5 min. GSH was calculated as the difference between total glutathione and GSSG. At a pH of 7.5 and standard/assay temperature of 25°C, the standard reduction potential of GSH/GSSG is −270 mV {=−240 mV+[(7.5−7.0)×−59] mV; Schafer and Buettner, 2001}. In order to calculate the tissue GSH reduction potential, we used E=−270 mV–30×log([GSH]2/[GSSG]), where [GSH] and [GSSG] are the molar concentrations of GSH and GSSG.
Thiobarbuturic acid reactive substances (TBARS)
Tissue TBARS was assessed using a commercially available kit (TBARS Parameter™ kit; R&D Systems). Approximately 30 mg of each tissue was homogenized in a buffer consisting of 100 mmol l−1 Tris-HCl, 2 mmol l−1 EDTA, and 5 mmol l−1 MgCl2·6H2O at pH 7.75 followed by centrifugation at 1600 g for 10 min at 4°C. A part of the supernatant was reserved and stored at −80°C for the analysis of total scavenging oxidative capacity (TOSC; described below) and protein concentration (Bradford, 1976), while the remaining supernatant was used immediately for the determination of TBARS. For the analysis of TBARS, equal volumes of sample and acid reagent were combined and incubated for 15 min at room temperature. The sample was then centrifuged twice at 12,000 g for 4 min at 4°C and the supernatant was used for the assay. TBARS standards were prepared according to the manufacturer's instructions. The samples and standards were then incubated with TBA reagent (2:1 v/v) for 3 h at 50°C, and absorbance was measured at 532 nm. TBARS levels were normalized to tissue protein concentration.
Total oxidative scavenging capacity (TOSC)
TOSC was determined as the rate of H2O2 removal, where H2O2 in the sample was monitored with Amplex Ultrared (Invitrogen). 100 µl of assay buffer (0.1 mmol l−1 Amplex Ultrared, 1 U ml−1 horseradish peroxidase in 50 mmol l−1 sodium citrate at pH 6.0) was added to 50 µl catalase standard (standards between 0 to 12 U of catalase activity were prepared) or sample (same homogenate prepared for the TBARS assay) in a spectrophotometer plate and pre-read at 565 nm. 50 µl of 160 µmol l−1 H2O2 was then added to each well and the endpoint absorbance was measured after 5 min. The final TOSC was expressed as µmole H2O2 min−1 (equivalent to 1 unit of catalase activity) and was normalized to tissue protein concentration.
Calculations and statistical analyses
Tissue GSH redox potential, MitoP/MitoB, TBARS and TOSC were measured in the four separate normoxic control samples taken during the exposures for both species (1 h after MitoB injection, paired with the 6 h of hypoxia, 6 h of hyperoxia, and after an additional 1 h of recovery). As the normoxic values were not significantly different between these samples within a species (one-way ANOVA; data not shown), we grouped them into a single normoxia sample to compare with our O2 treatment and recovery samples. For all measurements, hypoxic and hyperoxic groups had separate normoxic controls except for the S. marmoratus GSH redox potential, where the same normoxic control was used for hypoxic and hyperoxic groups because of the small sample size.
Two-way ANOVAs were performed to compare the effect of species (O. maculosus and S. marmoratus) and treatment (normoxia, 6 h hypoxia/hyperoxia treatment and recovery from hypoxia/hyperoxia) on tissue GSH redox potential, MitoP/MitoB, TBARS and TOSC. TBARS and TOSC were normalized to normoxic control values in order to investigate species-specific responses to hypoxia and hyperoxia exposure. Significant species by treatment interactions were followed with post hoc Holm–Sidak multiple comparisons to investigate within-species treatment effects (see Tables S1 and S2).
RESULTS
Validation of use of MitoB in marine sculpins
Within 30 min of injection into resting normoxic O. maculosus, MitoB was detected in brain, gill and liver, after which there was an exponential loss of MitoB over 72 h after injection (Fig. 1A). MitoP, which is formed when MitoB is oxidized by mitochondrial H2O2, appeared quickly in the brain, gill and liver, but at much lower concentrations than MitoB (Fig. 1B). MitoP also showed an exponential loss over the 72 h after injection (Fig. 1B). In white muscle, MitoB and MitoP levels peaked at 1 h after injection, but tissue concentrations were much lower compared with the other tissues sampled (highest levels at 1 h with 0.167 pmol MitoB and 0.06 pmol MitoP mg−1 tissue), so white muscle was not analysed for this study (see Fig. S3). Owing to similar patterns of MitoP and MitoB loss from tissues, MitoP/MitoB was generally stable over 72 h of consistent normoxic treatment (Fig. 1C), but the higher levels of MitoP and Mito B detected over the first 8 h after injection resulted in a high signal to noise ratio; thus, all experiments were completed within 8 h of MitoB injection.
Brain
Hypoxia and hypoxia+recovery
In both O. maculosus and S. marmoratus there were no significant effects of hypoxia or hypoxia followed by 1 h of normoxic recovery (hereafter ‘hypoxia-recovery’) on GSH redox potential or MitoP/MitoB ratio (Fig. 2A,B, see Table 1 for ANOVA results). There was a significant effect of treatment and species on brain TBARS as well as a significant treatment by species interaction (Fig. 2C; Table 1). Post hoc analysis revealed that O. maculosus had significantly higher TBARS during hypoxia exposure than in normoxia (P=0.0012) or when exposed to hypoxia+recovery (P=0.017), whereas in S. marmoratus there were no significant effects of hypoxia or hypoxia+recovery on brain TBARS. There were no significant effects of treatment or species on brain TOSC (Fig. 2D; Table 1).
Hyperoxia and hyperoxia+recovery
There was no effect of species, or hyperoxia and hyperoxia followed by 1 h of normoxic recovery (hereafter ‘hyperoxia+recovery’) on GSH redox potential (Fig. 2E; Table 2 for ANOVA results). There was a significant species effect on MitoP/MitoB, where O. maculosus had a higher MitoP/MitoB compared with S. marmoratus across treatments (Fig. 2F; Table 2). There was a significant effect of species and a significant treatment by species interaction on TBARS (Fig. 2G; Table 2). Post hoc analysis revealed that O. maculosus had higher TBARS after hyperoxia exposure than in normoxia (P=0.0013), whereas there were no changes in TBARS in the brain of S. marmoratus in any treatment. There was a significant species effect on TOSC levels: O. maculosus had higher TOSC levels than S. marmoratus (Fig. 2H). There was no significant treatment effect on brain TOSC.
Liver
Hypoxia and hypoxia+recovery
There was a significant effect of species, but no effect of treatment on liver GSH redox potential (Fig. 3A). The liver from O. maculosus had a more oxidized GSH redox potential than S. marmoratus (Fig. 3A). There was also a significant species effect on MitoP/MitoB, where O. maculosus had higher MitoP/MitoB than S. marmoratus (Fig. 3B). There was no effect of treatment on MitoP/MitoB. Neither species showed changes in liver TBARS and TOSC levels (Fig. 3C,D; Table 1).
Hyperoxia and hyperoxia+recovery
There were no effects of species or hyperoxia treatment on liver GSH redox potential (Fig. 3E). There were significant species and treatment effects on MitoP/MitoB: O. maculosus had higher MitoP/MitoB than S. marmoratus, which increased during hyperoxia in both species (Fig. 3F). There were no effects of species or treatment on liver TBARS (Fig. 3G). There was a significant species effect and a treatment by species interaction on TOSC levels. Post hoc analysis revealed that S. marmoratus had significantly higher TOSC levels in the hyperoxia treatment compared with normoxia (P=0.0089) and O. maculosus showed no changes with treatment (Fig. 3H).
Gill
Hypoxia and hypoxia+recovery
In the gills, there were no species or treatment effects on GSH redox potential (Fig. 4A). There was a significant effect of treatment on MitoP/MitoB, where MitoP/MitoB was lower in the hypoxia and hypoxia+recovery treatments compared with the normoxic controls for both species (Fig. 4B). There was a significant species effect on TBARS, where S. marmoratus generally showed higher TBARS than O. maculosus, especially in the hypoxia and hypoxia+recovery treatments (Fig. 4C). There were no changes in TOSC levels with treatment in either species (Fig. 4D).
Hyperoxia and hyperoxia+recovery
There was no species or treatment effect on gill GSH redox potential (Fig. 4E; Table 2). There was a significant treatment effect on MitoP/MitoB, but no species effect. MitoP/MitoB progressively increased in both species following hyperoxia and hyperoxia+recovery (Fig. 4F). There were no changes to TBARS and TOSC levels in either species exposed to hyperoxia and hyperoxia+recovery treatments (Fig. 4G,H).
DISCUSSION
The goal of the present study was to determine whether there were differences in in vivo ROS metabolism during exposure to hypoxia, hyperoxia and subsequent normoxic recovery between two species of sculpins that differ in hypoxia tolerance and inhabit environments that differ in the frequency and severity of O2 fluctuations. We hypothesized that the more hypoxia-tolerant species would be able to minimize the impacts of environmental O2 variability on ROS metabolism such that oxidative stress would be lower than in the less hypoxia-tolerant species. Our analysis demonstrated that the more hypoxia-tolerant O. maculosus exhibited more oxidative damage in the brain following exposure to both hypoxia and hyperoxia (which recovered in normoxia) compared with the less tolerant S. marmoratus. In comparison, the responses in the liver and gills were relatively subtle compared with those in the brain of O. maculosus, and S. marmoratus showed more changes in response to O2 fluctuations in these tissues. Overall, these data do not support our original hypothesis, which is likely because of the complex tissue- and species-specific responses in ROS metabolism during environmental O2 fluctuations. The relationship between ROS metabolism and hypoxia tolerance may also be further complicated by the role that ROS can play in cell signaling.
Our validation study demonstrated that MitoB uptake into isolated mitochondria from O. maculosus was similar to that observed in Drosophila (Fig. S1; Cochemé et al., 2011). When MitoB was injected into O. maculosus, MitoB and MitoP were detected in brain, gill and liver within 30 min of injections. The timeline of MitoB loss in marine sculpins, however, differed from that of the freshwater brown trout (Salin et al., 2015), possibly because of differences in osmoregulatory strategies. In brown trout, MitoB was still at high levels in liver at 72 h post-injection (94% of concentration at first measured time point of 3 h; Salin et al., 2017), whereas the levels in the liver of O. maculosus were much lower at 72 h (0.5% of concentration at 0.5 h; Fig. 1). Despite the decreases in MitoB in tissues over time, the levels of MitoB present in the tissues of sculpins allowed us to reliably detect MitoB and MitoP for the required 8 h duration of our study. There were also differences in tissue MitoP/MitoB between sculpins and brown trout, with brown trout white muscle showing a much higher ratio than in liver (0.27 compared with 0.073, respectively), whereas in the white muscle of sculpins, there was a low uptake of MitoB, which rendered the calculation of MitoP/MitoB unreliable (Fig. S3; Salin et al., 2017). Although we cannot explain the species-specific differences in white muscle MitoB uptake, the low uptake of MitoB in sculpin white muscle is consistent with its low mitochondrial content. It is also interesting to note that MitoB uptake was observed in the sculpin brain, which was not observed in mammals, likely due to the presence of a tight blood–brain barrier. Teleost fish are believed to also have a tight endothelial-based blood–brain barrier analogous to that of other vertebrates, but the uptake of MitoB in sculpin brain suggests that there may be functional differences in the blood–brain barrier between teleosts and mammals that allows MitoB uptake (Kniesel and Wolburg, 2000; Wolburg et al., 1983). These differences in MitoB uptake emphasize the importance of validating its use in different organisms and tissues, but overall, the results of our validation study suggest that MitoB can be used to reliably detect in vivo H2O2 accumulation in brain, liver and gill of sculpins over an 8 h exposure window.
Net in vivo ROS accumulation is a product of ROS generation, which is affected by cellular redox environment, and ROS scavenging, which is also affected by the redox environment and the antioxidants present (Aon et al., 2010). Tissue-specific differences in ROS-induced oxidative damage can be ascribed to variation in ROS generation capacity and susceptibility to oxidation (Murphy et al., 2011). In order to draw meaningful comparisons of ROS metabolism across tissues and between species, it is important to consider the drivers of ROS generation and the subsequent downstream effects of ROS accumulation, in addition to measuring the amount of ROS generated. Therefore, in this comparative study, we chose to quantify representative measures of redox environment, mitochondrial ROS, ROS damage and scavenging capacity. Although the indices we chose are not all-encompassing (e.g. there are numerous cellular redox couples and we only measured one), they provide insight into the important components of ROS metabolism across tissues and reveal interspecific differences in whole-animal responses to O2 fluctuations.
In both species of sculpin and in all tissues examined, there were no O2 treatment-associated changes in redox environment (as assessed by GSH redox potential), and no large changes in mitochondrial H2O2 in brain and liver. Our previous work in isolated mitochondria from these species demonstrated an association between experimentally manipulated GSH redox potential and H2O2 production, whereby additions of extramitochondrial GSH (to induce reductive stress) resulted in increased in H2O2 production (Lau and Richards, 2019). We also demonstrated that mitochondria from these two species possessed similar capacities to buffer extramitochondrial reductive stress. In combination, these studies suggest that environmentally relevant changes in O2 do not induce large changes in tissue GSH redox state in either species, which could explain why there are no large and consistent changes in H2O2 in vivo. It is important to point out, however, that our sample size for the GSH redox analysis is low as a result of tissue sample limitations, so the precise relationship between in vivo GSH redox potential and H2O2 production should be viewed with caution. Furthermore, mitochondrial H2O2 emission did not differ between species except under oligomycin-induced state 4 respiration state (Lau and Richards, 2019), which represents a very specific set of conditions (i.e. where ATP synthase is completely inhibited) that do not occur in vivo. Thus, our observations at the whole-animal level appear to be consistent with our previous in vitro study.
Despite the modest changes in mitochondrial H2O2 detected in all tissues using MitoB, we observed signs of oxidative stress in some tissues. For example, the large increases in TBARS in the brain of O. maculosus following hypoxia and hyperoxia exposure were not accompanied by significant changes in MitoP/MitoB during hypoxia exposure and recovery (Fig. 2B), and the changes of MitoP/MitoB during hyperoxia and recovery were subtle and only differentiated by species (Fig. 2F). There are several possible explanations for this discrepancy. First, the oxidative damage observed in the brain of O. maculosus could have been caused by ROS released by ETS complex sites oriented toward the cytoplasmic side (i.e. complex III and glycerol-3-phosphate dehydrogenase; Quinlan et al., 2013) and not detected by MitoB, which is localized to the mitochondrial matrix. Second, ROS other than H2O2 (e.g. superoxide and hydroxyl radicals which are normally quickly scavenged and converted to less reactive H2O2) and reactive species other than ROS (e.g. reactive sulfur and nitrogen species; DeLeon et al., 2016) may be causing oxidative damage in response to O2 stress. Third, non-mitochondrial sites of ROS production, e.g. microsomes and peroxisomes, may be responsible for the observed oxidative damage (Brown and Borutaite, 2012; Chance et al., 1979). Despite these caveats that come with using MitoB to detect mitochondrial H2O2, the lack of a direct connection between ROS generation and presumed downstream effects of ROS is not unique to this study and has been observed in other organisms. For example, in molluscs exposed to anoxia+recovery and hypoxia+recovery, changes in aconitase activity (the inhibition of which is directly related to the presence of superoxide) did not mirror changes in total antioxidant capacity nor measures of oxidative damage (Ivanina and Sokolova, 2016). These results illustrate the importance of taking a multifaceted approach to understanding ROS metabolism since measuring only one aspect, e.g. oxidative damage, may potentially be misleading.
While we predicted that the more O2-sensitive brain would show less oxidative damage compared with liver and gills, especially in the hypoxia-tolerant O. maculosus, in fact, we observed the opposite, where the brain of O. maculosus accumulated high levels of TBARS during hypoxia and hyperoxia exposure. In contrast, the same O2 stressors did not result in any measurable TBARS in the brain of S. marmoratus. Tissue specificity in oxidative damage has been observed in other studies (Johannsson et al., 2018; Lushchak et al., 2005; Lushchak and Bagnyukova, 2007), with no consistent pattern of tissues being similarly impacted with changes to environmental O2. In addition, there were no consistent patterns among the types of oxidative damage detected, e.g. protein carbonyls, lipid peroxides and TBARS, within a tissue during the hypoxia+recovery exposures in these studies. Species and tissue differences in the extent and type of oxidative damage observed could be related to inherent variation in their susceptibility to oxidative damage, e.g. due to differences in membrane lipid composition (Tribble et al., 1992) or level of misfolded proteins (Dukan et al., 2000). Our assessment of oxidative damage was limited to TBARS and it is possible that other markers of lipid peroxidation (e.g. 4-hydroxynonenal and lipid peroxides) or markers of other types of oxidative damage (e.g. protein carbonyls, mtDNA damage) would offer additional insights into the downstream effects of ROS. Differences in antioxidant mechanisms could also contribute to explaining species and tissue variation in the levels of oxidative damage measured (Leveelahti et al., 2014; Welker et al., 2012). In the present study, only the liver of S. marmoratus experienced changes in TOSC in response to hyperoxia. The lack of change in TOSC in the other tissues in both species could indicate that the existing levels of antioxidants were sufficient to minimize H2O2 accumulation and oxidative damage, and/or that other antioxidants mechanisms targeting non-H2O2 reactive species were involved, since our assessment of TOSC only captures H2O2 metabolism.
Interpreting the species-specific effects of changes in environmental O2 on ROS metabolism is further complicated by the fact that ROS and the products of oxidative damage can play important roles in cell signaling (Ayala et al., 2014; Sies, 2017; Wong et al., 2008). For example, malondialdehyde (the lipid peroxidation product quantified by the TBARS assay) and other lipid peroxidation adducts can act as secondary messengers of free radicals that subsequently affect protein function and structure, after which they can be enzymatically metabolized (Ayala et al., 2014; Zarkovic et al., 2013). As a result, the increases in TBARS observed in the brain of O. maculosus during hypoxia and hyperoxia exposure, could be part of a signaling process to coordinate physiological responses to changes in environmental O2 levels. In organisms that inhabit the environmentally dynamic higher intertidal, such as O. maculosus, these putative signaling mechanisms may be important and possibly confer a survival advantage. Consistent with the idea that TBARS may serve as a signaling molecule in the brain of O. maculosus is the fact that TBARS levels quickly recovered to control values after 1 h of normoxic recovery from hypoxia, possibly because this tissue possesses a higher capacity to repair lipid peroxidation. Also, O. maculosus displayed this transient response despite being exposed to a relatively less severe level of hypoxia compared with S. marmoratus. Further studies are necessary to directly test whether lipid peroxidation products and other products of oxidative damage play an important role in cell signaling during exposure to environmental O2 variability.
One of the main challenges associated with interpreting tissue responses to whole-animal hypoxia or hyperoxia exposure across species is taking into account their species-specific physiological modifications to the O2 transport cascade. For example, when species are exposed to the same level of environmental hypoxia and hyperoxia, physiological adjustments to the O2 transport cascade may not yield similar tissue PO2. Indeed, it must be noted here that our choice of hypoxic PO2 used herein represents a more severe hypoxic stress for S. marmoratus than O. maculosus, but it is interesting that the observed responses are generally greater in O. maculosus than S. marmoratus, perhaps due to a hypoxia-induced inhibition of metabolism in the latter species. Attempts to scale hypoxia exposures to species-specific indices of hypoxia tolerance (e.g. Pcrit) do not appear to fully normalize tissue O2 exposures or physiological responses (Mandic et al., 2014; Speers-Roesch et al., 2012). The impact of these species-specific physiological adjustments on tissue PO2 should, however, be minimized at the respiratory surface that interacts directly with the environment. Indeed, in both species examined here, gill mitochondrial H2O2 production during hypoxia and hyperoxia exposures were similar, suggesting similar levels of O2 stress. During hypoxia exposure, the decrease in ROS production observed in both species is likely due to the lower availability of O2 substrate for ROS production (Fig. 4B); however, recovery from hypoxic stress is generally thought to lead to an increase of H2O2 generation. This was not observed in the sculpin gills and in fact both species decreased mitochondrial H2O2 in response to hypoxia+recovery (Fig. 4B). The lower H2O2 after hypoxia+recovery could be due to an active suppression of mitochondrial H2O2 generation in both species exposed to a similar hypoxic PO2 at the gills. In S. marmoratus, however, this decrease in mitochondrial H2O2 appears to be concomitant with higher TBARS levels (although there was only a significant species effect, with treatment effect approaching significance, P=0.071; Table 1) during both hypoxia and hypoxia+recovery compared with O. maculosus, suggesting that S. marmoratus still experienced oxidative damage. Hyperoxia and hyperoxia+recovery yielded the expected increases in mitochondrial H2O2 in both species, as O2 is abundant and it further increases in normoxic recovery, suggesting that ROS production can be sensitive to O2 fluctuations even at supraphysiological O2 tensions (Fig. 4F).
Conclusions
Overall, we observed that the more hypoxia-tolerant sculpin O. maculosus differed in tissue responses to O2 variability when compared with the less hypoxia-tolerant S. marmoratus. We observed in vivo responses in the brain that suggests that hypoxia tolerance may not necessarily be related to an overall decrease in ROS generation to minimize oxidative stress, but these responses were also dependent on the tissue and the PO2 in question. While lowering ROS emission and minimizing oxidative damage has been suggested as a putatively adaptive response to O2 limitation in other hypoxia-tolerant animals [assessed in gastrocnemius muscle in high-altitude deer mice (Peromyscus maniculatus) and liver in killifish; Du et al., 2016; Mahalingam et al., 2017], our in vivo and in vitro observations (Lau and Richards, 2019) of the hypoxia-tolerant O. maculosus suggests that minimizing ROS and oxidative damage may not be a universal adaptive strategy to inhabiting O2 variable environments. The species and tissue specificity observed in this study and others is likely due to the multi-dimensional nature of environmental O2 variation that animals experience in diverse habitats (Mandic and Regan, 2018), and potentially complicated by the fact that ROS can also act as a signaling molecule (Dickinson and Chang, 2011).
Acknowledgements
We would like to thank Joshua Emerman, Monica Goldade, Dr Tammy Rodela, Derek Somo and Andrew Thompson for their help with fish collections, and to Dr Rodela and D. Somo for helping with tissue sampling.
Footnotes
Author contributions
Conceptualization: G.Y.L., J.G.R.; Methodology: G.Y.L., S.A.; Validation: G.Y.L., S.A.; Formal analysis: G.Y.L., S.A.; Investigation: G.Y.L.; Resources: M.P.M., J.G.R.; Writing - original draft: G.Y.L., J.G.R.; Writing - review & editing: G.Y.L., S.A., M.P.M., J.G.R.; Visualization: G.Y.L.; Supervision: M.P.M., J.G.R.; Funding acquisition: M.P.M., J.G.R.
Funding
This work was funded by a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery grant to J.G.R., and by a Wellcome Trust Investigator Award (110159/Z/15/Z) to M.P.M. G.Y.L. was supported by a NSERC Canada Graduate Scholarship, Michael Smith Foreign Supplement, and a travel grant from The Company of Biologists.
References
Competing interests
The authors declare no competing or financial interests.