The physiological mechanisms that limit thermal tolerance are broadly relevant to comparative biology and global change. Species differences in macromolecular stability play important roles in evolved patterns of heat tolerance, but other mechanisms such as oxidative stress have also been hypothesized to contribute. For example, mussels in the genus Mytilus exhibit evolved physiological differences at several levels of organization that have been linked with interspecific differences in whole-organism heat tolerance. Both omics and behavioral studies suggested that variation in resistance to oxidative stress plays a role in these differences. Functional data are needed to test this hypothesis. Here, we compared three Mytilus congeners to examine whether susceptibility to oxidative stress contributes to acute heat tolerance. We assayed the activity of two antioxidant enzymes (catalase, superoxide dismutase), as well as levels of oxidative damage to lipids, DNA and individual proteins (using gel-based proteomics methods). In addition, we assessed these oxidative stress responses after repeated episodes of heat stress experienced in air or while immersed in seawater, given that survival and competitive outcomes between Mytilus congeners differ in these two contexts. The results are generally inconsistent with patterns that would be expected if oxidative stress contributes to thermal sensitivity. Rather, the more heat-tolerant congeners suffer comparable or even elevated levels of oxidative damage. As predicted, different treatment contexts led to distinct changes in proteome-wide abundance patterns and, to a lesser extent, protein carbonylation profiles. Overall, the results question the relevance of oxidative damage as a mediator of heat tolerance in this genus.

Elucidating the mechanisms that establish thermal optima and limits can shed light on the factors governing organisms' distribution patterns and their sensitivities to global change (Dong et al., 2022; Pörtner, 2002; Somero et al., 2017; Williams et al., 2016). Adaptive evolved differences in thermal sensitivities of the major classes of macromolecules – proteins, nucleic acids and membranes – play a central role in thermal adaptation (Dong and Somero, 2009; Somero et al., 2017), but it remains less clear to what extent other mechanisms might also contribute to adaptive differences among species. Oxidative stress figures prominently among the other mechanistic factors that have been hypothesized to contribute to species differences in thermal sensitivity (Abele et al., 1998; Feidantsis et al., 2020; Pörtner, 2010; Tomanek, 2011), although the available evidence does not uniformly support a direct link between oxidative stress and thermal tolerance (e.g. Mueller et al., 2014). Oxidative stress results from the imbalance between the production of highly reactive metabolic byproducts (e.g. reactive oxygen species, ROS) and scavenging of these harmful byproducts by antioxidant systems before they inflict macromolecular damage to nucleic acids, lipids or proteins (Halliwell and Gutteridge, 2007). If left unrepaired, oxidative damage can greatly impair cellular function. It is often observed that ROS production and/or oxidative damage markers increase at elevated body temperatures (Abele et al., 1998, 2001; Bagnyukova et al., 2007; Gleason et al., 2017; Heise et al., 2006, 2003; Jimenez et al., 2016; Wang et al., 2018). Because future environmental patterns are expected to include more unpredictable and more extreme temperatures (e.g. Harris et al., 2018; Stillman, 2019), it is imperative to understand more fully the mechanistic underpinnings of thermal tolerance, including the role of oxidative stress.

The rocky intertidal environment has proven to be an excellent system for studying adaptive differences in environmental stress tolerance. The likelihood that oxidative stress also may exert a strong impact on intertidal zone animals stems from patterns of environmental variation experienced by these species. Dictated by the rhythms of the tides and their interactions with solar and atmospheric factors (Denny and Harley, 2006; Helmuth, 1998), this habitat exposes organisms to episodic extremes of temperature, as well as wide variation in oxygen availability and food supply. Variations in these key environmental factors are found across latitude and at much smaller scales, such as across vertical gradients. Not surprisingly, even closely related inhabitants of this dynamic habitat display sometimes striking differences in their susceptibility to abiotic stressors including temperature (e.g. Dong et al., 2022, 2008; Dowd and Somero, 2013; Stillman and Somero, 2000; Tomanek and Helmuth, 2002; Tomanek and Somero, 2000).

The three congeners of Mytilus mussels found on the Pacific coast of North America vary dramatically in several aspects of thermal physiology, from the molecular to the organismal level (reviewed in Lockwood and Somero, 2011), which makes them a particularly useful system for studying the mechanistic basis of inter-specific differences in thermal tolerance. Mytilus trossulus Gould is less heat tolerant than Mytilus galloprovincialis Lamarck or Mytilus californianus Conrad, as illustrated by differences in survival after acute, high-temperature episodes (Dowd and Somero, 2013) or following chronic exposure to moderate temperatures (Schneider, 2008). The heat-tolerant M. galloprovincialis has a lower resting heart rate and a higher critical temperature for cessation of cardiac activity than M. trossulus (Braby and Somero, 2006b); the critical temperature for cessation of cardiac activity is even higher for M. californianus (Logan et al., 2012). At the level of single macromolecules, a handful of metabolic enzymes of M. galloprovincialis, as well as cytosolic malate dehydrogenase for M. californianus, have been shown to maintain the appropriate configuration for substrate binding affinity at higher temperatures, unlike the same enzymes from M. trossulus (Fields et al., 2006; Lockwood and Somero, 2012). The species also vary in their molecular responses to heat stress, such as the threshold for induction of synthesis of heat-shock proteins and patterns of protein ubiquitination (Buckley et al., 2001; Hofmann and Somero, 1996). Comparative transcriptomic and proteomic analyses of M. trossulus and M. galloprovincialis have further characterized the temperature-dependent expression profiles that might underpin their differences in heat tolerance. Specifically, these two species differ in their relative induction of small heat shock proteins that stabilize the cytoskeleton and, most relevant to this study, in their induction of constituents of the antioxidant system (Fields et al., 2012; Lockwood et al., 2010; Tomanek and Zuzow, 2010). Genomic sequences for chaperones and antioxidant system proteins also appear to be under strong selection in the more heat-tolerant M. galloprovincialis relative to M. trossulus (Popovic and Riginos, 2020). These specific examples of differences in expression patterns and biochemistry are reflected in a higher frequency of stabilizing amino acid substitutions across the proteome of M. galloprovincialis than in M. trossulus (Popovic and Riginos, 2020).

The evolutionary and biogeographic history of the genus Mytilus has also been reconstructed (Koehn, 1991; Seed, 1992), providing important context for the interpretation of inter-specific physiological differences (Lockwood and Somero, 2011). The California mussel, M. californianus, arose in the relatively cool waters of the eastern Pacific. In this region, M. californianus experiences episodic exposure to high body temperatures only at low tide (Denny et al., 2011). The blue mussel M. trossulus also originated in the eastern Pacific. Finally, the blue mussel M. galloprovincialis diverged more recently from the Mytilus edulis lineage in the more saline and warmer Mediterranean Sea (Seed, 1992). These evolutionary events have been followed by a period of human-mediated spread, through both intentional (e.g. introductions for aquaculture) and unintentional means (e.g. transport in ballast water; Carlton and Geller, 1993; Wilkins et al., 2008). Mytilus galloprovincialis has been a particularly successful invader around the globe (Branch and Steffani, 2004; Daguin and Borsa, 1999), and in the eastern Pacific it has displaced native M. trossulus throughout the southern, warmer portion of its former range (McDonald and Koehn, 1988). Enhanced heat tolerance of M. galloprovincialis is a key factor in its invasive success (Lockwood and Somero, 2011).

Whereas differences in thermal responses at the level of cardiac function, gene expression and protein thermal stability have been shown to correlate with whole-organism heat tolerance and with global (re)distribution patterns of Mytilus congeners, the putative contribution of oxidative stress resistance to these inter-specific patterns has not been functionally evaluated. The factors influencing oxidative stress in many intertidal animals are complex, however, because episodes of elevated body temperature tend to coincide with times when oxygen availability is low. For example, in intertidal mussels, high body temperatures occur when the tide is low, air temperature and/or solar irradiance is high, and wave action is minimal. To counter the detrimental consequences of desiccation, many intertidal animals close the shell valves (Miller and Dowd, 2017) or otherwise seal up during aerial episodes of heat stress. This behavioral response limits the availability of oxygen to the tissues and forces reliance on less efficient anaerobic pathways in the cytosol (Bayne et al., 1976; Shick et al., 1988; Widdows and Shick, 1985). Hypoxia is known to coincide with ROS production in a variety of aquatic organisms (Welker et al., 2013), and the combination of hypoxia and acute temperature increases might be expected to inflate levels of oxidative stress (Castro et al., 2020; Sun et al., 2020). For example, M. galloprovincialis gill cells produced more ROS upon exposure to 30°C in air than to 15°C in air (Wang et al., 2018). Other toxic byproducts of glycolysis such as methylglyoxal (derived from glucose) might also be produced. In addition, mussels and other intertidal organisms quickly reopen the shell valves upon reimmersion (Dowd and Somero, 2013; Miller and Dowd, 2017). This behavior coincides with a drop in body temperature and leads to a rapid influx of oxygen that is analogous to reperfusion after ischemia. This rapid reperfusion has been identified as a key contributor to oxidative macromolecular damage in stroke and disease models (Li and Jackson, 2002), as well as in comparative physiology models (Bickler and Buck, 2007; Lushchak et al., 2005), and it could represent a key period of oxidative stress in mussels (Rivera-Ingraham et al., 2013). For example, Gracey et al. (2008) identified a temporal pattern of gene expression in M. californianus that may be linked to cycles of oxidative stress. Furthermore, in behavioral studies, we have observed that all three species studied here close up during acute heat stress, but the more heat-sensitive congener M. trossulus appears to avoid further cycles of hypoxia–reoxygenation by remaining open during recovery periods between stressful temperature episodes. This behavioral evidence, combined with the omics results described above, led us to hypothesize that M. trossulus suffers more oxidative damage than its more heat-tolerant congeners during episodes of heat stress (Dowd and Somero, 2013).

Building on this substantial background knowledge, here we addressed two outstanding questions regarding inter-specific differences in acute heat tolerance among Mytilus congeners. First, what role do oxidative stress and oxidative damage to macromolecules play in these observed differences? To more directly address the hypothesis that Mytilus congeners vary in resistance to oxidative damage during/after exposure to episodes of heat stress, we compared several functional measures across the three Mytilus species: activities of two antioxidant defense enzymes (catalase and superoxide dismutase); and the amount of lipid oxidation and DNA oxidation, and two-dimensional gel-based proteomics measures of oxidative protein damage (amino acid side-chain carbonylation). We predicted that more heat-tolerant congeners would have higher antioxidant capacities and/or accumulate less oxidative damage from acute heat stress episodes. We also expected that the heat-sensitive M. trossulus might accumulate oxidative damage to key proteins that would otherwise play protective roles during/after heat stress; carbonylation is known to impede protein function (e.g. Grimsrud et al., 2007; Pierce et al., 2008). Second, how do metrics of oxidative stress in these mussels change as a result of repeated exposure to thermal stress in air versus seawater? Earlier Mytilus omics studies examined responses to single episodes of heat stress in seawater only (Lockwood et al., 2010; Tomanek and Zuzow, 2010); we extended our analyses to include one and three consecutive daily episodes of thermal stress in seawater or in air. The experimental design employed here better reflects conditions as they are experienced by mussels in the field (Denny et al., 2011; Miller and Dowd, 2017), and it attempts to capture species differences that might not appear after a single episode or when mussels are immersed in seawater throughout the exposure. Recent monitoring and modeling studies show that thermal stress events often cluster into multi-day windows for intertidal species; these relatively rare ‘heat waves’ coincide with midday low tides when air temperature is high and wave action is low (Harley, 2008; Mislan and Wethey, 2015; Seuront et al., 2019). Moreover, both behavior and survival of these Mytilus congeners exhibit clear differences between thermal stress in air and thermal stress in seawater (Dowd and Somero, 2013), and ecological and genetic studies suggest that the outcomes of competition among Mytilus congeners can vary between these environmental contexts (McDonald et al., 1991; Schneider and Helmuth, 2007; Wonham, 2004).

The results reported here highlight numerous differences in the physiological consequences of heat stress in air versus heat stress in seawater; these differences may help to explain idiosyncrasies of Mytilus distribution and invasion patterns. However, the results largely contradict predictions of the hypothesis that oxidative stress mediates acute heat tolerance, calling into question its relevance in this study system.

Animal collection and maintenance

We studied three Mytilus congeners from the Pacific coast of North America that vary in their heat tolerance. Mytilus trossulus (heat sensitive) and M. galloprovincialis (heat tolerant) have been the subject of extensive comparative studies outlined above; we included heat-tolerant M. californianus to allow stronger inferences and to avoid spurious patterns in a two-species comparison (Garland and Adolph, 1994). Adult M. californianus were collected from a single location in the middle rocky intertidal zone at Hopkins Marine Station (HMS) in Pacific Grove, CA, USA (36.6216°N, 121.9042°W). In order to avoid hybrid blue mussels – HMS sits within the hybrid zone for the native and invasive congeners (Braby and Somero, 2006a; Hilbish et al., 2010) – M. trossulus and M. galloprovincialis were collected from docks outside the hybrid zone in Newport Harbor, OR, USA (44.6226°N, 124.0520°W), and Santa Barbara Harbor, CA, USA (34.4089°N, 119.6893°W), respectively, and transported to HMS over ice and in air-filled coolers. Based on the collection site features, the blue mussels had been acclimatized in the wild to constant immersion, whereas M. californianus had been acclimatized to rhythmic immersion and emersion according to the lunar cycle at HMS. All three species were subsequently held under common conditions of constant immersion in flow-through seawater tanks at HMS (salinity 34±1 ppt; temperature 12–15°C) for at least 4 weeks – and no longer than 12 weeks – prior to the start of any experiment. This allowed sufficient time for acclimation to common-garden conditions (Buckley et al., 2001). Although field studies have documented differences in mean size among these species where they co-occur (Braby and Somero, 2006a; Dutton and Hofmann, 2008), experiments herein were carried out on individuals of comparable size (M. californianus 61–89 mm; M. galloprovincialis 55–78 mm; M. trossulus 50–78 mm). Sex of the animals was not determined. Mussels were fed 3 times a week with a mixture of coccolithophorids (Isochrysis spp., Pavlova spp.), green algae (Tetraselmis spp.) and diatoms (Thalassiosira weissflogii) (Shellfish Diet 1800©, Reed Mariculture). Mussels were not fed during the 3 days of the episodic heat stress exposures.

Because the two blue mussels cannot be distinguished morphologically, we verified the putative species identifications for those species (M. trossulus from Newport, OR, USA; M. galloprovincialis from Santa Barbara, CA, USA) using the polymorphic molecular markers Glu-5′ and ITS-1. These marker regions were amplified with established primers (Heath et al., 1995; Rawson et al., 1996). Only one putative M. trossulus was mis-identified; its molecular markers were consistent with M. galloprovincialis. This individual was excluded from analyses.

Heat ramp experiments

Heat ramp experiments were conducted in the same manner for each species in both seawater and air (as in Dowd and Somero, 2013). Seawater heat ramps were conducted in an insulated cooler. Seawater temperature was controlled with a programmable laboratory water bath temperature controller (Cole-Parmer Polystat®). The outflow freshwater circulation from the water bath was passed through a countercurrent heat exchanger, while a pump in the cooler passed seawater in the opposite direction through the heat exchanger and back to the cooler. Air heat ramps were carried out in insulated polystyrene boxes in a humidified, walk-in chamber using established methods (Logan et al., 2012). Air temperature was regulated with a programmable benchtop temperature controller (Omega CSC32) connected to a mat heating element (Omega Model SRMU020230, 120 V, 1.25 A, 5.08×76.2 cm). A fan inside the box ensured adequate air mixing, and ice packs in the boxes were used to dampen oscillations in the air temperature otherwise caused by cycling of the heating element. In all cases, the animals were held at a control temperature (13–14°C) within the acclimation range until just before the start of the heat ramp. Temperatures were increased and decreased gradually from 13–14°C at 16°C h−1 for all species; heating rates as high as 20.2°C h−1 have been recorded for intertidal M. californianus at HMS on the warmest days, but they are more typically less than 8°C h−1 (Denny et al., 2011; Miller and Dowd, 2017). In each experiment, the animals were held at the peak temperature (seawater 33±1°C; air 32±1°C) for 45 min before cooling began. The peak temperature of 32–33°C approximates the highest temperature measured at a low-shore site at HMS (Miller and Dowd, 2017). Mussels were returned to control temperature seawater immediately at the conclusion of the heat ramp. In each experiment, individuals were sampled after exposure to one or three heat ramps. Heat waves of 3 days are typical of the shoreline at HMS (Denny et al., 2011). Successive ramps started at the same time (between 12:00 h and 13:00 h) on subsequent days. For simplicity, we refer to treatments as SWX and AirX, where X is the number of heat stress episodes (1 or 3).

Individuals in each treatment group were sampled immediately prior to the beginning of the next day's heat ramp. Thus, they were given ∼20 h of recovery in control conditions. There were two reasons for this design. First, if oxidative damage does contribute to species differences in survival, which play out over days to weeks (Dowd and Somero, 2013), we would expect metrics of damage to accumulate and persist in the more heat-sensitive congener. Second, the design captures mussels' physiological states at the time when the next stressful period of heat/immersion would begin.

Tissue sampling

We focused our analyses on gill tissue. Gill has high ciliary activity for suspension feeding (Jørgensen, 1974) and is the primary site of oxygen uptake, so we expected that consequences of changes in oxygen availability might be most pronounced in this tissue. Furthermore, gill function is essential for nitrogenous waste excretion (Thomsen et al., 2016). Thus, damage to gill tissue could have pronounced effects on the whole organism. Gill has also been the focus of earlier Mytilus transcriptomics and proteomics studies (Gracey et al., 2008; Lockwood et al., 2010; Tomanek and Zuzow, 2010), providing important background for comparison with the present results.

At the appropriate time in each experiment, gill tissue of each individual was dissected out and flash-frozen in liquid nitrogen, after which tissues were stored at −80°C until analysis. Aside from the DNA oxidation assay, all subsequent assays maintained the individual mussel as the unit of analysis (i.e. samples were not pooled).

Antioxidant enzyme activity assays

All enzymatic assays were conducted at a standard cuvette temperature of 25±0.2°C on a UV-1601 spectrophotometer (Shimadzu, Kyoto, Japan), equipped with a temperature-controlled cell that was regulated by a Lauda RM6 recirculating water bath (Brinkmann, Westbury, NY, USA). Tissues were homogenized in a bead homogenizer at a constant 1:4 ratio of tissue to potassium phosphate buffer (50 mmol l−1, pH 6.8 at 4°C), then centrifuged at 18,000 g and 4°C for 30 min. The supernatant was removed, centrifuged a second time for 5 min, and then separated into aliquots for the different enzyme activity assays.

Total superoxide dismutase (SOD; EC 1.15.1.1) activity (including both the cytosolic Cu-Zn and mitochondrial Mn-dependent forms) and catalase (CAT; EC 1.11.1.6) activity were assayed on the same gill homogenates. SOD was assayed in triplicate for each sample using the xanthine oxidase/cytochrome c method (modified from Dowd et al., 2013; McCord and Fridovich, 1969). The reaction mixtures contained 10 µmol l−1 cytochrome c and 50 µmol l−1 xanthine in 50 mmol l−1 potassium phosphate buffer containing 0.1 mmol l−1 EDTA, at pH 7.6. We included 16 U CAT enzyme in each reaction to prevent reoxidation of cytochrome c by the H2O2 produced by the SOD reaction (Crapo et al., 1978). Baseline reactions were started by adding xanthine oxidase to generate sufficient superoxide ion (the substrate of SOD) to cause an increase in A550 of ∼0.020 min−1 due to the reduction of cytochrome c. One unit of SOD activity was defined as the quantity needed to achieve 50% competitive inhibition of this baseline rate, and activity in mussel samples was standardized per gram fresh mass. CAT activity was measured in triplicate for each sample by following the rate of disappearance of hydrogen peroxide (decrease in absorbance at 240 nm) in the same phosphate buffer solution (Beers and Sizer, 1952; Dowd et al., 2013). One unit of CAT activity was defined as the quantity decomposing 1 µmol of H2O2 (molecular extinction coefficient 0.0394 l mmol−1 cm−1) in 1 min at pH 7.6 (note that this is not the standard pH of 7.0 in most of the literature). CAT specific activity was standardized in units of per gram fresh mass.

Lipid oxidative damage (TBARS) assay

Lipid peroxidation represents one of the most common forms of macromolecular damage under exposure to oxidative stress. The thiobarbituric acid reactive substances (TBARS) assay was run following established methods (Hermes-Lima and Storey, 1995). Lipids were extracted in 1.1% phosphoric acid and then incubated with 1% thiobarbituric acid (TBA) in 50 mmol l−1 NaOH. Blanks were treated with 3 mmol l−1 hydrochloric acid instead of TBA. Butylated hydroxytoluene (BHT) was included in the assay buffer at 10 µmol l−1 because of its chain-breaking capacity (Pikul et al., 1983); this modification reduces artefactual increases in peroxides induced by heating samples during the preparation. After heating for 15 min at 100°C, lipids were extracted in isobutanol. We measured blank (replacing TBA with 3 mmol l−1 HCl) and TBA-treated absorbance at 532 nm and 600 nm on a microplate reader in duplicate for each sample. Sample equivalents of malondialdehyde (MDA) concentrations were determined from a standard curve of 1,1,3,3-tetramethoxypropane treated in the same fashion; this chemical is converted to the detectable standard MDA via acid hydrolysis during the reaction.

DNA oxidative damage assay

Oxidative damage to DNA was quantified using the New 8-OHdG Check kit (Cosmobio USA), which measures the amount of 8-hydroxy-2′-deoxyguanosine (8-OHdG) in tissues based on competition with a known standard. 8-OHdG is a commonly measured and relatively stable product of DNA oxidation (Baek et al., 2000; Halliwell and Gutteridge, 2007; Liu et al., 2004). Briefly, DNA was extracted from tissues, pooled and quantified into equal total amounts (33 µg) from 5 individuals in each treatment, digested to single nucleotides with P1 nuclease, eluted in buffer and then assayed in duplicate on a microplate with bound anti-8-OHdG antibody, all as per the kit manufacturer's instructions. Concentrations were estimated from a logit standard curve using kit standards. Because of limited nucleic acid material, there was only one replicate per treatment/species combination.

Statistical analyses for enzymatic activity and oxidative damage

Data for the preceding assays were analyzed in R (http://www.R-project.org/) using linear models that included species, treatment and a species-by-treatment interaction term. For the enzymes, we initially included shell length as a covariate, but likelihood ratio tests showed that this term did not significantly improve the model (SOD P=0.599; CAT P=0.389; TBARS P=0.319). Therefore, shell length was removed from final analyses. SOD activity data were log10 transformed to meet assumptions.

Protein abundance and protein oxidative damage (gel-based proteomics)

We first attempted to quantify total oxidative protein damage using immunochemical dot-blots for protein carbonylation, a stable marker of protein oxidation (Stadtman, 1993), with an antibody reactive with dinitrophenylhydrazine (DNPH)-derivatized amino acid side-chains (Dalle-Donne et al., 2003). However, we encountered substantial background signal. The DNPH derivatization reaction targets all carbonyl groups regardless of their molecular source (DNA, protein, etc.). We could not abolish this background by treatment of the samples with streptomycin sulfate to precipitate DNA (Luo and Wehr, 2009) or by protein precipitation in 20% trichloroacetic acid (Colombo et al., 2016). Furthermore, these artefacts remained when we excluded the detergents Triton X-100 and Nonidet P-40, which can contain DNPH-reactive impurities and promote the release of contaminating DNA from the nucleus (Luo and Wehr, 2009), from the homogenization buffer. We also removed dithiothreitol (DTT), which via the Fenton reaction with iron can cause artefactual increases in carbonyl measurements (Luo and Wehr, 2009; Stadtman, 1993), with similar lack of effect. We suspect this persistent background results from the accumulation of carbonyl-containing pigments such as fucoxanthins from the algal food. Numerous pigments, many of which contain carbonyl groups, have been identified in Mytilus mussels (Campbell, 1970; Hertzberg et al., 1988), and gill is a heavily pigmented tissue.

We were interested as much in variation in the specific targets of oxidative protein damage among Mytilus congeners as in variation in the total amount of oxidative damage. Therefore, we used a two-dimensional gel proteomics approach with subsequent Western blotting to detect specific carbonylated proteins (Chaudhuri et al., 2006; McDonagh et al., 2005). This technique allowed us to separate and identify individual proteins with stronger (or weaker) carbonyl signal following thermal stress events relative to the control conditions, while excluding background signal from pigments or other impurities. For these analyses, we focused on a comparison of control animals with those subjected to three episodes of heat stress in air or in water (n=4 mussels per treatment per species for proteomics).

2D gel electrophoresis

The 2D gel electrophoresis protocol followed established methods (Dowd et al., 2010a,b). Briefly, mussel gill proteins were extracted by homogenization in RIPA buffer. Gill protein extracts from each individual were processed in duplicate for all subsequent steps for M. galloprovincialis and M. trossulus; one replicate was used for quantifying protein abundance while the other was used for western blotting for protein carbonylation. For M. californianus, we only processed gels for identifying protein carbonylation. Salts were removed via precipitation in 20% trichloroacetic acid (w/v) in acetone; we processed aliquots of homogenate carrying 750 µg of protein each. Proteins were then resuspended in urea/thiourea rehydration buffer and loaded onto 17 cm immobilized pH gradient (IPG) strips (pH 3–10 non-linear) by overnight passive rehydration at 4°C, prior to isoelectric focusing (IEF) as previously described (Dowd et al., 2010a,b). The two IPG strips for each mussel were run on the same day. In total, two 12-lane IEF runs were needed per species, with treatments evenly distributed across days. IPG strips were frozen at −80°C until further analysis. IPG strip equilibration and reduction (using DTT and IAA), second dimension SDS-PAGE gel runs, and colloidal Coomassie Blue G250 staining were all carried out as previously described for one of the two IPG strips for each individual (Dowd et al., 2010a,b). Prior to equilibration and reduction, the second IPG strip per mussel was washed in 10 mmol l−1 DNPH in 2 mol l−1 HCl for 15 min (McDonagh and Sheehan, 2006); this reaction derivatized protein carbonyl groups to a stable dinitrophenylhydrazone for later detection. The DNPH derivatization reaction was halted by addition of 2 mol l−1 Tris buffer twice for 7 min each. Derivatization at this stage of the process prevents differential IEF of proteins that might be caused by the DNP adduct (Reinheckel et al., 2000), facilitating downstream image processing and comparative analysis of total protein gels with anti-DNPH blots. Proteins were separated in the second dimension on 16 cm (12% acrylamide) SDS-PAGE gels in Bio-Rad Protean II cells at 200 mA constant current for ∼6.5 h. Protein abundance gels were submerged in a colloidal Coomassie Blue G-250 staining solution and gently rotated in the dark overnight, followed by several rounds of destaining in ultrapure water. Gels were scanned on an Epson V750 Professional scanner using the film/transparency adapter at 600 dpi resolution and with automatic gain control turned off.

2D western blotting for protein carbonylation

For each individual mussel, the DNPH-derivatized protein gel was immediately used for large-format western blotting. Proteins were transferred to PVDF membrane overnight at 15 V and 4°C. Membranes were subsequently blocked with 4% blotting grade milk in phosphate-buffered saline containing 0.05% Tween-20 (PBS-T), incubated with primary antibody at 1:10,000 dilution in the blocking solution (rabbit anti-DNP keyhole limpet hemocyanin; Molecular Probes catalog number A-6430), washed 4 times in PBS-T, incubated with secondary antibody at 1:5000 dilution in the blocking solution (goat anti-rabbit HRP-conjugated; Santa Cruz Biotech catalog number sc-2004), and washed again 4 times in PBS-T. Blots were visualized by addition of ECL Reagent (Amersham Pharmacia) for 1 min, followed by exposure to blue light autoradiography film for 10 min (M. galloprovincialis and M. trossulus) or 2 min (M. californianus). Films were immediately scanned for densitometry and image analysis using the same methods as for the gels.

Image analysis and statistics

Protein abundance gel and carbonylation blot image analyses were carried out using Delta 2D software v4.0 (Decodon, Germany) as previously described (Dowd et al., 2010a,b). Images for the protein abundance gels and anti-DNP protein carbonylation blots were analyzed separately, because of considerable differences in the number of protein spots visible using the two methods. Images were overlaid manually using the Delta 2D ‘exact’ warping strategy, with extensive manual confirmation of alignment (typically 200–300 manual warping assignments per gel pair). In each type of analysis, spot volumes (intensity×area) were normalized by the total volume of all spots visible in the union image. We replaced the few instances of missing data with a small value of 0.001 to avoid statistical issues introduced by missing values in the data matrix. The normalized data were then arcsine-square root transformed to better meet the assumptions of downstream statistical tests (Dowd, 2012). To enable spot-picking for proteins of interest in the carbonylation analysis from abundance gels, we manually aligned each species' union carbonylation blot image to the union protein abundance image using Delta2D.

The proteomics data were analyzed spot-wise in the Matlab R2017b bioinformatics package using a one-way ANOVA across the three treatments, with a strict alpha value of 0.01 as the threshold statistical significance to avoid false positives. Each species was analyzed separately for each type of data (protein abundance and carbonylation signal). Post hoc comparisons were computed with a least significant difference (LSD) test. As a complementary approach to account for the false discovery rate when conducting numerous statistical tests, we applied the q-value method (Storey and Tibshirani, 2003) with a threshold q<0.1 within each analysis. The carbonylation blots had fewer protein spots than the abundance gels, and no protein carbonylation spots met the q-value threshold for any of the three species. Lastly, to visualize proteome-wide treatment differences in abundance and carbonylation, we performed species-specific, unsupervised principal components analyses (PCA) and plotted individual mussels' PCA scores in 3D PC space.

Protein identification via mass spectrometry

We excised protein spots of interest using 1 mm biopsy punches, focusing on clearly identifiable spots that were most differentially expressed or contributed most to treatment differences in carbonylation. One set of protein spots were analyzed on an Applied Biosystems 4700 Proteomics Analyzer MALDI-TOF-TOF mass spectrometer at the University of California, Davis, USA, using previously described methods (Dowd et al., 2010a,b). Bioinformatics searches used both MS and MS/MS spectra to search against the NCBI and SwissProt Mollusca entries in Mascot (Perkins et al., 1999), with mass error tolerances of 50 ppm for MS fragments and 0.5 Da for MS/MS ions, and a maximum of 1 missed cleavage.

There were some examples of constitutive or inducible differences among the three mussel species in their antioxidant capacities and/or levels of oxidative damage, but these were not systematically different between M. trossulus and its more heat-tolerant congeners. Furthermore, repeated episodes of heat stress in seawater or in air had only muted effects on the antioxidant and oxidative damage profiles of the three congeners. However, there were several instances in which responses to heat stress in air differed dramatically from responses to heat stress in seawater.

CAT activity

We predicted a lower or less inducible activity of the antioxidant enzyme CAT in the heat-sensitive M. trossulus relative to either M. californianus or M. galloprovincialis. The activity of catalase was significantly impacted by species (P=0.002; Table S1), but not entirely in the predicted manner. As predicted, CAT activity of M. galloprovincialis exceeded that of M. trossulus in the control temperature seawater, SW1, Air1 and Air3 treatments. However, CAT activity of M. galloprovincialis also exceeded that of M. californianus in the control, SW1 and Air3 treatments (Fig. 1). CAT activity did not differ between heat-tolerant M. californianus and heat-sensitive M. trossulus. There was no overall treatment effect on CAT activity (P=0.933), but the interaction term was only marginally non-significant (P=0.064). For M. galloprovincialis, CAT activity after three episodes of heat stress in seawater was significantly reduced relative to that in all but the Air1 treatment (Fig. 1B; Table S1).

Fig. 1.

Catalase (CAT) antioxidant enzyme activity in three Mytilus mussel congeners sampled in control conditions and after one or three episodic exposures to heat stress in air or in seawater. (A) Mytilus californianus, (B) Mytilus galloprovincialis and (C) Mytilus trossulus. SWC, control temperature seawater; SWX and AirX indicate the number (X=1 or 3) of heat stress episodes in seawater or in air (33±1°C seawater; 32±1°C air). Full details of treatments are given in Materials and Methods. The number of points (corresponding to individual mussels) equals the sample size. In this and all subsequent boxplots, the central line indicates the median, the lower and upper bounds of the box indicate the 25th and 75th percentiles, respectively, and the whiskers indicate 1.5× the interquartile range.

Fig. 1.

Catalase (CAT) antioxidant enzyme activity in three Mytilus mussel congeners sampled in control conditions and after one or three episodic exposures to heat stress in air or in seawater. (A) Mytilus californianus, (B) Mytilus galloprovincialis and (C) Mytilus trossulus. SWC, control temperature seawater; SWX and AirX indicate the number (X=1 or 3) of heat stress episodes in seawater or in air (33±1°C seawater; 32±1°C air). Full details of treatments are given in Materials and Methods. The number of points (corresponding to individual mussels) equals the sample size. In this and all subsequent boxplots, the central line indicates the median, the lower and upper bounds of the box indicate the 25th and 75th percentiles, respectively, and the whiskers indicate 1.5× the interquartile range.

SOD activity

We also predicted a lower and/or less-inducible activity of the antioxidant SOD in M. trossulus. The data only partially supported this prediction. The linear model indicated significant effects of species (P<0.001), treatment (P=0.033) and a species-by-treatment interaction (P=0.017) (Table S2). The three species had statistically indistinguishable SOD activity in control conditions, but they responded differently to heat stress episodes (Fig. 2). SOD activity in M. galloprovincialis exceeded that of M. trossulus only in SW1. Similar to CAT activity, SOD activity in M. galloprovincialis exceeded that of M. californianus in SW1, Air1 and Air3. Contrary to the prediction, SOD activity in heat-tolerant M. californianus (Fig. 2A) and heat-sensitive M. trossulus (Fig. 2C) remained similar across the five treatments.

Fig. 2.

Superoxide dismutase (SOD) antioxidant enzyme activity in three Mytilus mussel congeners sampled in control conditions and after one or three episodic exposures to heat stress in air or in seawater. (A) Mytilus californianus, (B) M. galloprovincialis and (C) M. trossulus. See Fig. 1 and Materials and Methods for details. The number of points (corresponding to individual mussels) equals the sample size.

Fig. 2.

Superoxide dismutase (SOD) antioxidant enzyme activity in three Mytilus mussel congeners sampled in control conditions and after one or three episodic exposures to heat stress in air or in seawater. (A) Mytilus californianus, (B) M. galloprovincialis and (C) M. trossulus. See Fig. 1 and Materials and Methods for details. The number of points (corresponding to individual mussels) equals the sample size.

Lipid oxidative damage (TBARS)

If accumulation of oxidative damage to lipids contributes to reduced heat tolerance, we would predict greater levels of this damage marker in heat-sensitive M. trossulus. There was a significant effect of species on levels of TBARS (P=0.002; Table S3), but not in the expected direction. On average, levels of lipid oxidation in heat-tolerant M. galloprovincialis were roughly 1.5–2 times higher than those in the other two congeners (Fig. 3). We confirmed the constitutive excess of TBARS in M. galloprovincialis relative to M. trossulus in an independent experiment (W.W.D., unpublished data). Furthermore, there was no accumulation of TBARS lipid peroxidation products in any of the species after 3 days of elevated body temperature in either air or seawater (Ptreatment=0.442; Ptreatment×species=0.411).

Fig. 3.

Lipid oxidation in three Mytilus mussel congeners sampled in control conditions and after three episodic exposures to heat stress in air or in seawater. (A) Mytilus californianus, (B) M. galloprovincialis and (C) M. trossulus. Lipid oxidation was measured using a thiobarbituric acid reactive substances (TBARS) assay (MDA, malondialdehyde). See Fig. 1 and Materials and Methods for details. The number of points (corresponding to individual mussels) equals the sample size.

Fig. 3.

Lipid oxidation in three Mytilus mussel congeners sampled in control conditions and after three episodic exposures to heat stress in air or in seawater. (A) Mytilus californianus, (B) M. galloprovincialis and (C) M. trossulus. Lipid oxidation was measured using a thiobarbituric acid reactive substances (TBARS) assay (MDA, malondialdehyde). See Fig. 1 and Materials and Methods for details. The number of points (corresponding to individual mussels) equals the sample size.

DNA oxidative damage (8-OHdG)

Limited nucleic acid material precluded replication of this assay, but using pooled samples, we found that gill 8-OHdG levels in heat-tolerant M. californianus were more than double those of both heat-sensitive M. trossulus and heat-tolerant M. galloprovincialis in control conditions. Levels of 8-OHdG in M. californianus remained higher in the seawater and air heat stress treatments (Fig. 4). There was little indication of accumulation of oxidative DNA damage following multiple episodes of heat stress in any of the three congeners. These results contradict the prediction that heat-sensitive M. trossulus should suffer more DNA oxidative damage.

Fig. 4.

Levels of DNA oxidation in three Mytilus mussel congeners sampled in control conditions and after three episodic exposures to heat stress in air or in seawater. (A) Mytilus californianus, (B) M. galloprovincialis and (C) M. trossulus. DNA oxidation was measured using 8-hydroxy-2′-deoxyguanosine (8-OHdG). See Fig. 1 and Materials and Methods for details. Samples were pooled for one replicate per species and treatment.

Fig. 4.

Levels of DNA oxidation in three Mytilus mussel congeners sampled in control conditions and after three episodic exposures to heat stress in air or in seawater. (A) Mytilus californianus, (B) M. galloprovincialis and (C) M. trossulus. DNA oxidation was measured using 8-hydroxy-2′-deoxyguanosine (8-OHdG). See Fig. 1 and Materials and Methods for details. Samples were pooled for one replicate per species and treatment.

Protein abundance and protein oxidative damage

We quantified protein abundance and signals of protein oxidative damage (amino acid side-chain carbonylation) in control conditions and after three episodes of heat stress in air and in seawater. We predicted that a greater number of proteins, and/or specific proteins important in thermal defenses, would suffer increased levels of carbonylation following repeated episodes of heat stress in the heat-sensitive congener M. trossulus. Both protein abundance gels and carbonylation blots were generally highly reproducible across biological replicates (Fig. S1). One M. trossulus abundance gel (in the seawater heat treatment) and one M. trossulus carbonylation blot (in the control treatment) were excluded from analyses because of poor protein separation and image quality.

Protein abundance patterns

The air and seawater heat treatments induced largely distinct changes in protein abundance in both M. galloprovincialis and M. trossulus (abundance data were not collected for M. californianus). A total of 51 of the 1340 protein spots on M. galloprovincialis abundance gels were significantly differentially expressed among at least two of the three treatments at P<0.01 (20 of these with a q-value<0.10) (Table 1; Table S4). Of the 51, only 17 proteins increased in abundance relative to controls in both heat treatments, and 6 proteins decreased in abundance in both heat treatments (including spot 31357 identified as 60 s ribosomal protein l13a; Table 2). A larger proportion of the differentially expressed proteins (33 of 51) varied between the two heat treatments. Notably, four protein spots for M. galloprovincialis were strongly upregulated in the seawater heat treatment and were identified with mass spectrometry as variants of Hsp70 or glucose-regulated protein 78 (Table 2); each of these chaperone proteins was approximately half as abundant after heat stress in air. Spot 124369, identified as M. edulis heavy metal binding protein (c1q domain containing), was similarly much more abundant after heat stress in seawater than in air (Table 2). There was clear treatment separation of protein abundance profiles within principal component (PC) space, particularly along the first dimension for which the repeated seawater heat stress treatment had elevated scores (Fig. 5B).

Fig. 5.

Principal component analysis coordinates for proteome-wide protein abundance (left) and protein carbonylation signal (right) for individual mussels in each of three treatments. (A) Mytilus californianus, (B,C) M. galloprovincialis and (D,E) M. trossulus. In each analysis, the first three principal components (PCs) explained more than 64% of the variance in the data (up to 84.3%). Stems below points connect to coordinates in the PC1–PC2 plane. Number of points of each color equals the sample size.

Fig. 5.

Principal component analysis coordinates for proteome-wide protein abundance (left) and protein carbonylation signal (right) for individual mussels in each of three treatments. (A) Mytilus californianus, (B,C) M. galloprovincialis and (D,E) M. trossulus. In each analysis, the first three principal components (PCs) explained more than 64% of the variance in the data (up to 84.3%). Stems below points connect to coordinates in the PC1–PC2 plane. Number of points of each color equals the sample size.

Table 1.

Proteomics summary of differentially expressed proteins (left) or differentially carbonylated proteins (right) in three Mytilus mussel congeners in each of three pairwise treatment comparisons

Proteomics summary of differentially expressed proteins (left) or differentially carbonylated proteins (right) in three Mytilus mussel congeners in each of three pairwise treatment comparisons
Proteomics summary of differentially expressed proteins (left) or differentially carbonylated proteins (right) in three Mytilus mussel congeners in each of three pairwise treatment comparisons
Table 2.

Mytilus mussel proteins identified via mass spectrometry and their pairwise treatment comparison ratios

Mytilus mussel proteins identified via mass spectrometry and their pairwise treatment comparison ratios
Mytilus mussel proteins identified via mass spectrometry and their pairwise treatment comparison ratios

Of the 1374 spots on M. trossulus protein abundance gels, 57 were significantly differentially expressed in any comparison at P<0.01 (27 of these with a q-value<0.10) (Table 1; Table S4). More proteins increased in abundance in heat stress treatments relative to controls than decreased; 18 proteins were elevated in both seawater and air heat stress treatments (including 47924 Hps70; 47944 Hsp70-1; 47906 heat shock protein 70 B2-like; Table 2). Unlike for M. galloprovincialis, the M. trossulus proteins that changed significantly in both heat treatments tended to be more highly upregulated in air than in seawater (6 of 18 instances) or upregulated to equal degrees in the two heat treatments (10 of 18 proteins) (Table S4). Another 6 proteins decreased in abundance in both heat treatments relative to controls. As for M. galloprovincialis, a larger proportion of differentially expressed proteins (41 of 57) varied between the two heat treatments than differed between either heat treatment and the control. Accordingly, there was clear separation of proteome-wide abundance patterns among the three treatments in PC space, indicating distinct expression patterns in each treatment (Fig. 5D).

Protein carbonylation patterns

Overall, there was little evidence for induction of protein carbonylation by either repeated heat treatment. There was also little variation in overall sensitivity to this form of damage among the three congeners (Table 1). Instead, the data indicate a broadly similar pattern of protein carbonylation in control conditions and after repeated heat exposures for each Mytilus congener.

Only 2 of the 112 consistently detected carbonylated protein spots on M. californianus blots were significantly different in any comparison (P<0.01). One of these increased in signal relative to the control following seawater heat stress, and the other decreased in the same treatment (Table 1; Table S4). Because we did not have independent protein abundance gels from which to pick spots, the identity of these two proteins is unknown. The first three PCs explained 84.35% of the variance but revealed no separation of the treatments in PC space (Fig. 5A).

Similar to M. californianus, protein carbonyl signals for M. galloprovincialis were largely unchanged following air or seawater heat treatment. Only 10 of 507 carbonylated protein spots were significantly different in any comparison (P<0.01) (Table 1; Table S4). In 12 of the 15 total pairwise instances of significant changes, carbonylation signals of individual proteins decreased in the repeated heat stress treatments relative to the controls (Table S4). No protein carbonyl signals increased in both seawater and air heat treatments; 5 significantly decreased by similar amounts in both heat treatments (including spot IDs 67231, 70506 and 67341, the last of which was identified as actin A3; Table 2). Thus, there is no evidence of a significant rise in carbonylation. These carbonylation signal changes were independent of changes in protein abundance. In 9 of 10 cases, we were able to map the differentially carbonylated spots to a corresponding spot on protein abundance gels, and none of those 9 proteins changed significantly in abundance between any two treatments (Table S4). As for M. californianus, there was little evidence of treatment separation of proteome-wide carbonylation patterns in PC space for M. galloprovincialis (Fig. 5C).

The most heat-sensitive congener, M. trossulus, also exhibited very little change in protein carbonylation signals following either heat treatment. Only 2 of 422 carbonylated protein spots were significantly different in any comparison for M. trossulus (Table 1; Table S4). One spot (ID 46774) had 4.7- and 7.3-fold increases in carbonylation signal after seawater and air heat treatment, respectively. However, the matching spot on protein abundance gels (ID 47920) was 8.72 and 12.01 times more abundant in the heat treatments than in controls (Table S4). As a result, the ratio of carbonyl signal to total protein decreased for this spot in both heat treatments. Based on the location of this spot adjacent to M. trossulus ID 47924, which was identified as a Hsp70, we believe but have not confirmed that this protein is also a Hsp70. Mytilus galloprovincialis protein spots 30950 and 30954 were also found in the same gel positions and identified as Hsp70 variants. The other significant protein (ID 46705) had decreased carbonylation signal in the seawater heat treatment. In agreement with the few statistically significant changes and in line with the other congeners, there was no separation of the three treatments in PC space for M. trossulus (Fig. 5E).

Differential resistance to oxidative stress has received attention in the literature as a possible mediator of variation in thermal tolerance (Kassahn et al., 2009; Pörtner, 2002), including in the Mytilus mussel congeners studied here (Popovic and Riginos, 2020; Tomanek and Zuzow, 2010). Behavioral observations had also led us to speculate that M. trossulus' relative avoidance of periods of valve closure during recovery from episodes of heat stress was a means to avoid additional hypoxia–reoxygenation cycles and concomitant oxidative stress (Dowd and Somero, 2013). However, using repeated heat treatments that are known to induce differential mortality among the three species (Dowd and Somero, 2013), we found little evidence to support the hypothesis that oxidative damage contributes to differences in acute heat tolerance among three Mytilus congeners.

Species differences in susceptibility to heat-induced oxidative damage

Species differences in these data are generally inconsistent with the hypothesis that oxidative stress limits thermal tolerance among Mytilus congeners, at least under the realistic, episodic exposure regimes we employed. The primary pieces of evidence in support of the hypothesis are higher (CAT) or more inducible (SOD) antioxidant capacities in heat-tolerant M. galloprovincialis than in heat-sensitive M. trossulus. Other work has revealed similar heat inducibility of antioxidant defenses in M. galloprovincialis (Georgoulis et al., 2021; Wang et al., 2018). However, these elevated antioxidant capacities, which are combined with the genomic features and expression differences described above, do not protect M. galloprovincialis from having equal or sometimes higher (in the case of TBARS) levels of oxidative macromolecular damage. This pattern is somewhat surprising, given the fact that M. trossulus has a faster aerobic ‘pace of life’ than M. galloprovincialis based on citrate synthase enzyme capacity (Lockwood and Somero, 2011; Fig. S2), which might be expected to lead to greater cumulative levels of ROS formation in M. trossulus. However, ROS formation and aerobic metabolism are often only weakly coupled (Koch et al., 2021; Monaghan et al., 2009). Our conclusion that oxidative stress is not a major contributor to differences in heat tolerance among Mytilus congeners is further supported by the results from the second relatively heat-tolerant species M. californianus, which exhibits antioxidant capacities on a par with those of heat-sensitive M. trossulus yet suffers similar (TBARS) or greater (8-OHdG) quantities of oxidative damage.

Furthermore, there was little evidence for a pronounced increase in any of the measures of oxidative damage to macromolecules, or for induction of antioxidant defenses, in the more thermally sensitive M. trossulus after exposure to repeated (three) episodes of heat stress in air or seawater. Again, our temperature treatments are known to induce substantial, but delayed, mortality in M. trossulus (Dowd and Somero, 2013). Thus, our data strongly suggest that accumulation of oxidative damage is not a key contributor to this mortality. Metrics of antioxidant capacity and oxidative damage also were little changed in the other two Mytilus congeners. Similar results indicating no effect of different thermal treatments on antioxidant parameters or total protein carbonyls were observed within a species of intertidal limpet (Drake et al., 2017). The 2D protein carbonyl blots similarly provide no quantitative evidence for a broad increase in oxidative protein damage in the heat-sensitive M. trossulus relative to its more heat-tolerant congeners (Table 1).

In lieu of a widespread increase in protein oxidation, we predicted that the individual proteins that are susceptible to heat-induced carbonylation in M. trossulus might be important in thermal defenses. Protein carbonylation is known to be limited to a subset of the proteome (Cabiscol et al., 2000), with the most susceptible proteins including those with catalytic motifs that contain transition metals as well as enzymes associated with energy production and conversion (Chaudhuri et al., 2006; Maisonneuve et al., 2009). For example, SOD and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were selectively carbonylated in yeast treated with hydrogen peroxide (Costa et al., 2002); GAPDH, aconitase and the alpha subunit of ATP synthase were carbonylated in yeast mutants lacking functional MnSOD (O'Brien et al., 2004). Other studies have similarly observed oxidative damage to protein components of the antioxidant system (Bruno et al., 2009), which could have the undesirable effect of dampening a cell's ability to defend itself against further oxidative damage (O'Brien et al., 2004).

The one protein in M. trossulus for which the carbonylation signal increased significantly after repeated episodes of heat stress in air and in seawater was putatively identified as the molecular chaperone hsp70. This raises the intriguing possibility that one of the most important proteins for maintenance of proteome integrity during heat stress is itself susceptible to heat-induced oxidative damage. However, the abundance of the corresponding protein spot increased more following heat stress than did the carbonylation signal, indicating that the relative amount of protein oxidation actually decreased following heat stress in this inducible chaperone. In vitro studies of purified M. trossulus hsp70 chaperone function before and after oxidation would be needed to evaluate the functional impacts of carbonylation. Other oxidative stress studies using similar methods have identified molecular chaperones as specific targets of carbonylation, including age-related damage to the protein-folding chaperone BiP in mammalian liver (Rabek et al., 2003), hsp70 in a copper toxicity study in yeast (Shanmuganathan et al., 2004) and hsp70 in a neurocognitive impairment model (Sultana et al., 2010).

Rather than implying a role for oxidative damage in differential heat tolerance, our protein abundance data confirm that air or seawater exposure to heat stress induces changes in protein expression patterns whose primary role may be to promote the conformational integrity of the proteome. Indeed, we observed strong upregulation of several proteins identified as heat shock proteins in both M. galloprovincialis and M. trossulus, as observed in previous Mytilus transcript and protein expression studies (Feidantsis et al., 2020; Fields et al., 2012; Lockwood et al., 2010; Tomanek and Zuzow, 2010). These results offer additional support for the notions that susceptibility to protein destabilization/denaturation is a primary consequence of heat stress and that variation in this susceptibility is a key driver of Mytilus species differences in heat tolerance (Gracey et al., 2008; Hofmann and Somero, 1995; Lockwood and Somero, 2012). Temperature effects on macromolecular stability loom large as a selective force in the evolution of thermal tolerance (Liao et al., 2021; Somero et al., 2017). While oxidative stress clearly plays a role in some aspects of life-history evolution and tradeoffs (Janssens and Stoks, 2018; Koch et al., 2021; Monaghan et al., 2009), its role in acute thermal tolerance remains less clear.

Caveats for oxidative stress measures

Despite considerable evidence that heat stress increases ROS production and/or oxidative damage in marine organisms (Abele et al., 1998, 2001; Heise et al., 2003; Jimenez et al., 2016; Wang et al., 2018), our measurements indicate a lack of pronounced changes in any of the oxidative damage metrics following repeated episodes of acute heat stress. We turn our attention to potential caveats to this observation.

First, it is possible that three ∼4 h episodes of acute heat stress constitute insufficient time for the manifestation of changes in the antioxidant system or in oxidative damage metrics. However, previous studies have observed abrupt changes (within 1 h) in oxidative damage metrics over single episodes (e.g. Bagnyukova et al., 2007; Wang et al., 2018), including acute rises in gill lipid hydroperoxides during heat exposures in M. californianus (Jimenez et al., 2016) and M. galloprovincialis (Wang et al., 2018). It may well be that chronic exposure to high temperature induces fundamentally different changes in these measures to those following the episodic regime we employed. For example, the activity of antioxidants including CAT and SOD increased and TBARS accumulated in response to chronic thermal stress (>5 days) in seawater in M. galloprovincialis (Feidantsis et al., 2020). However, in the region from which we sampled mussels, episodic exposure to heat stress during low-tide periods is the most likely scenario to challenge the thermal tolerance of these species.

Second, some oxidative damage products might be dynamically cleared during the ∼20 h period between the end of the final heat episode and our sampling time, precluding their detection in our design. For example, lipid oxidation products increased transiently and recovered within 24 h in a fish (Bagnyukova et al., 2007), and lipid peroxidation and DNA 8-OHdG damage accumulated in the mussel Perna perna when exposed to air but returned to baseline levels within 3 h of re-immersion (Almeida et al., 2005). In contrast, 8-OHdG peaked 24 h after reperfusion in a mammalian ischemia model (Liu et al., 2004). However, we observed no significant changes in TBARS over a shorter-term time series in heat-sensitive M. trossulus in seawater; we sampled at the peak of the heat ramp, immediately at the end of the ramp, and after 7 h of recovery (data not shown; n=5, Psampling-time=0.541).

The situation for protein carbonylation is also complicated, in part because it is challenging to distinguish between carbonylation as a sign of damage (McDonagh and Sheehan, 2006) and carbonylation as the result of ROS-dependent signaling. Because they form via multiple mechanisms, protein carbonyls were long thought to be a comprehensive and stable index of oxidative protein damage (e.g. in disease states; Hecker and Wagner, 2018; Korolainen et al., 2007). However, reversal of protein carbonylation (decarbonylation) over short time scales of less than 30 min has been reported (Wong et al., 2008), putatively involving a thiol-dependent mechanism (Wong et al., 2013), although this process is not widely documented. Carbonylated proteins might also be degraded by the proteasome (Grune et al., 1997), preventing their accumulation. For example, we observed a decrease in carbonylation of an actin protein in M. galloprovincialis in both heat treatments, whereas in earlier work a mussel actin was deemed highly susceptible to carbonylation caused by pollution and hydrogen peroxide (McDonagh et al., 2005). Furthermore, carbonylation of specific amino acid side-chains is only one of many possible oxidative modifications of proteins, and some of these such as surface methionine oxidation are known to be reversible (Stadtman and Levine, 2003). There is some evidence that methionine oxidation and protein carbonylation may be linked (Moskovitz and Oien, 2010). Interestingly, in M. californianus, recent work suggested a peak in expression of a methionine-R-sulfoxide reductase protein during high tide (Elowe and Tomanek, 2021), perhaps indicative of a repair phase for methionine oxidation when oxygen and food are abundant.

With these caveats in mind, we cannot rule out transient rises and falls in the other measures of oxidative damage. Nonetheless, the near-complete absence of accumulation of any of the measured oxidative damage products after repeated episodes of thermal stress that are known to induce differential mortality strongly suggests that oxidative stress is not differentiating Mytilus congeners in their tolerance of acute heat stress. Instead, the intermittently warm and hypoxic lifestyle of all Mytilus mussels appears to have selected over evolutionary time either for resistance to ROS production and excess oxidative damage (Sokolov et al., 2021) or, as in the case of lipid oxidation, for favoring infrequent repair over the costly maintenance of a constitutive defense (Jimenez et al., 2016). To examine the generality of this conclusion that oxidative stress is not a major player in differential acute heat tolerance of Mytilus congeners, it will be necessary to gather similar comparative data from other genera that do not routinely experience episodic oxidative challenges in the same manner.

Is all heat stress the same?

Our results support the assertion that air and seawater heat stress impose unique challenges for Mytilus mussels and other intertidal animals (Dowd and Somero, 2013). The distinction is particularly evident in the protein abundance results (Fig. 5B,D). Interestingly, these contextual differences are also most prominent in the more heat tolerant of the blue mussels, M. galloprovincialis, for which repeated exposure to heat stress in seawater is most distinct among the three treatments for global protein expression (Fig. 5B) and is associated with reduced antioxidant enzyme capacities (Figs 1 and 2). The same SW3 treatment induced the greatest increase in mortality for both M. galloprovincialis and M. trossulus (Dowd and Somero, 2013).

It is possible that these context-dependent differences in the response to heat stress have played a role in the invasive spread of M. galloprovincialis on the Pacific coast of North America. For example, our protein abundance patterns (including stronger induction of heat shock proteins after seawater heat stress than after air heat stress) and antioxidant data provide some functional evidence supporting the field observation that M. galloprovincialis fares worse in subtidal locations (Schneider and Helmuth, 2007). Interestingly, in other parts of the world, M. galloprovincialis has invaded intertidal habitats first (Robinson and Griffiths, 2002). Successful invasions arise from numerous contributions, including the genetic and physiological features of the populations' founding individuals and different suites of abiotic and ecological constraints in different regions. For example, prior work on the Pacific coast of North America has shown that the two blue mussels differ in their responses to the combined effects of temperature and salinity. The native M. trossulus is generally more successful in areas where salinity is lower on average, even if those habitats have warmer seawater temperatures (Braby and Somero, 2006a). It would be informative to repeat these experiments at a range of salinities to further examine the physiology underlying this interaction. Teasing apart the physiological consequences of such interactions remains a substantial challenge in the study of species invasions.

Oxidative damage and mussel pace of life

Lastly, the relationships between aerobic metabolic capacity, as indexed by citrate synthase activity, and either heat tolerance or susceptibility to acute oxidative damage are not straightforward among these Mytilus congeners. Mytilus californianus appears to have the highest aerobic pace of life of the three congeners studied here (Fig. S2), despite being the most heat tolerant of the three congeners we tested, with a median lethal temperature above 38°C in our source population (Denny et al., 2011). By some measures, M. californianus also suffers higher baseline levels of oxidative damage (8-OHdG) than the other congeners. However, lipid TBARS values are highest in M. galloprovincialis despite it having the lowest gill aerobic capacity. Interestingly, M. californianus evolved at a lower average temperature in the eastern Pacific (driven by the cool California current) than did M. galloprovincialis in the Mediterranean Sea. Thus, adaptation of aerobic and antioxidant enzymatic capacities may simply reflect long-term average oxidative rate, rather than predict the consequences of acute thermal stress.

Conclusion

In closing, we found little support for the hypothesis that accumulation of oxidative damage products contributes to higher mortality following acute episodes of heat stress in the mussel M. trossulus than in its more heat-tolerant congeners M. galloprovincialis and M. californianus. In several other respects, these mussels' antioxidant and oxidative damage profiles differ following exposure to heat stress in air versus in seawater, perhaps contributing to context-dependent patterns of invasiveness. Further comparative work in other genera is needed to examine the generality of these conclusions, given the expected increases in frequency and intensity of acute thermal stress events with climate change.

Laurie Kost collected citrate synthase data for M. californianus and coordinated the 2011 PISCO ‘Biomechanics and Ecological Physiology of Intertidal Communities' course for blue mussel genotyping. Jerod Sapp and Stephen Gosnell assisted with collection and shipment of bay mussels. Brent Lockwood provided unpublished data on M. galloprovincialis and M. trossulus gill citrate synthase activities. Mass spectrometry samples were processed in the laboratory of Dietmar Kültz at the University of California, Davis. Mark Denny provided constructive comments on the manuscript.

Author contributions

Conceptualization: W.W.D., G.N.S.; Methodology: W.W.D.; Formal analysis: W.W.D.; Investigation: W.W.D.; Writing - original draft: W.W.D.; Writing - review & editing: G.N.S.; Visualization: W.W.D.; Supervision: G.N.S.; Project administration: G.N.S.; Funding acquisition: G.N.S.

Funding

This is contribution number 532 from PISCO, the Partnership for Interdisciplinary Studies of Coastal Oceans, funded by the David and Lucile Packard Foundation, and the Gordon and Betty Moore Foundation.

Data availability

Data are publicly available from the WSU Research Exchange repository: https://doi.org/10.7273/000004807.

Abele
,
D.
,
Burlando
,
B.
,
Viarengo
,
A.
and
Pörtner
,
H.-O.
(
1998
).
Exposure to elevated temperatures and hydrogen peroxide elicits oxidative stress and antioxidant response in the Antarctic intertidal limpet Nacella concinna
.
Comp. Biochem. Physiol.
120
,
425
-
435
.
Abele
,
D.
,
Tesch
,
C.
,
Wencke
,
P.
and
Pörtner
,
H. O.
(
2001
).
How does oxidative stress relate to thermal tolerance in the Antarctic bivalve Yoldia eightsi?
Antarct. Sci.
13
,
111
-
118
.
Almeida
,
E. A.
,
Bainy
,
A. C. D.
,
Dafre
,
A. L.
,
Gomes
,
O. F.
,
Medeiros
,
M. H. G.
and
Di Mascio
,
P.
(
2005
).
Oxidative stress in digestive gland and gill of the brown mussel (Perna perna) exposed to air and re-submersed
.
J. Exp. Mar. Biol. Ecol.
318
,
21
-
30
.
Baek
,
S.-H.
,
Kim
,
J.-Y.
,
Choi
,
J.-H.
,
Park
,
E.-M.
,
Han
,
M.-Y.
,
Kim
,
C.-H.
,
Ahn
,
Y.-S.
and
Park
,
Y.-M.
(
2000
).
Reduced glutathione oxidation ratio and 8 OHdG accumulation by mild ischemic pretreatment
.
Brain Res.
856
,
28
-
36
.
Bagnyukova
,
T. V.
,
Danyliv
,
S. I.
,
Zin'ko
,
O. S.
and
Lushchak
,
V. I.
(
2007
).
Heat shock induces oxidative stress in rotan Perccottus glenii tissues
.
J. Thermal Biol.
32
,
255
-
260
.
Bayne
,
B. L.
,
Bayne
,
C. J.
,
Carefoot
,
T. C.
and
Thompson
,
R. J.
(
1976
).
The physiological ecology of Mytilus californianus Conrad 1. Metabolism and energy balance
.
Oecologia
22
,
211
-
228
.
Beers
,
R. F.
and
Sizer
,
I. W.
(
1952
).
A spectrophotometric method for measuring the breakdown of hydrogen peroxide by catalase
.
J. Biol. Chem.
195
,
133
-
140
.
Bickler
,
P. E.
and
Buck
,
L. T.
(
2007
).
Hypoxia tolerance in reptiles, amphibians and fishes: life with variable oxygen availability
.
Annu. Rev. Physiol.
69
,
145
-
170
.
Braby
,
C. E.
and
Somero
,
G. N.
(
2006a
).
Ecological gradients and relative abundance of native (Mytilus trossulus) and invasive (Mytilus galloprovincialis) blue mussels in the California hybrid zone
.
Mar. Biol.
148
,
1249
-
1262
.
Braby
,
C. E.
and
Somero
,
G. N.
(
2006b
).
Following the heart: temperature and salinity effects on heart rate in native and invasive species of blue mussels (genus Mytilus)
.
J. Exp. Biol.
209
,
2554
-
2566
.
Branch
,
G. M.
and
Steffani
,
C. N.
(
2004
).
Can we predict the effects of alien species? A case-history of the invasion of South Africa by Mytilus galloprovincialis (Lamarck)
.
J. Exp. Mar. Biol. Ecol.
300
,
189
-
215
.
Bruno
,
M.
,
Moore
,
T.
,
Nesnow
,
S.
and
Ge
,
Y.
(
2009
).
Protein carbonyl formation in response to propiconazole-induced oxidative stress
.
J. Proteome Res.
8
,
2070
-
2078
.
Buckley
,
B. A.
,
Owen
,
M.-E.
and
Hofmann
,
G. E.
(
2001
).
Adjusting the thermostat: the threshold induction temperature for the heat-shock response in intertidal mussels (genus Mytilus) changes as a function of thermal history
.
J. Exp. Biol.
204
,
3571
-
3579
.
Cabiscol
,
E.
,
Piulats
,
E.
,
Echave
,
P.
,
Herrero
,
E.
and
Ros
,
J.
(
2000
).
Oxidative stress promotes specific protein damage in Saccharomyces cerevisiae
.
J. Biol. Chem.
275
,
27393
-
27398
.
Campbell
,
S. A.
(
1970
).
The carotenoid pigments of Mytilus edulis and Mytilus californianus
.
Comp. Biochem. Physiol.
32
,
97
-
115
.
Carlton
,
J. T.
and
Geller
,
J. B.
(
1993
).
Ecological roulette: the global transport of nonindigenous marine organisms
.
Science
261
,
78
-
82
.
Castro
,
J. S.
,
Braz-Mota
,
S.
,
Campos
,
D. F.
,
Souza
,
S. S.
and
Val
,
A. L.
(
2020
).
High temperature, pH, and hypoxia cause oxidative stress and impair the spermatic performance of the Amazon fish Colossoma macropomum
.
Front. Physiol.
11
,
772
.
Chaudhuri
,
A. R.
,
De Waal
,
E. M.
,
Pierce
,
A.
,
Van Remmen
,
H.
,
Ward
,
W. F.
and
Richardson
,
A.
(
2006
).
Detection of protein carbonyls in aging liver tissue: a fluorescence-based proteomic approach
.
Mech. Ageing Dev.
127
,
849
-
861
.
Colombo
,
G.
,
Clerici
,
M.
,
Garavaglia
,
M. E.
,
Giustarini
,
D.
,
Rossi
,
R.
,
Milzani
,
A.
and
Dalle-Donne
,
I.
(
2016
).
A step-by-step protocol for assaying protein carbonylation in biological samples
.
J. Chromatogr. B
1019
,
178
-
190
.
Costa
,
V. M. V.
,
Amorim
,
M. A.
,
Quintanilha
,
A.
and
Moradas-Ferreira
,
P.
(
2002
).
Hydrogen peroxide-induced carbonylation of key metabolic enzymes in Saccharomyces cerevisiae: the involvement of the oxidative stress response regulators Yap1 and Skn7
.
Free Radic. Biol. Med.
33
,
1507
-
1515
.
Crapo
,
J. D.
,
McCord
,
J. M.
and
Fridovich
,
I.
(
1978
).
Preparation and assay of superoxide dismutases
.
Meth. Enzymol.
53
,
382
-
393
.
Daguin
,
C.
and
Borsa
,
P.
(
1999
).
Genetic characterisation of Mytilus galloprovincialis Lmk. in North West Africa using nuclear DNA markers
.
J. Exp. Mar. Biol. Ecol.
235
,
55
-
65
.
Dalle-Donne
,
I.
,
Rossi
,
R.
,
Giustarini
,
D.
,
Milzani
,
A.
and
Colombo
,
R.
(
2003
).
Protein carbonyl groups as biomarkers of oxidative stress
.
Clin. Chim. Acta
329
,
23
-
38
.
Denny
,
M. W.
and
Harley
,
C. D. G.
(
2006
).
Hot limpets: predicting body temperature in a conductance-mediated thermal system
.
J. Exp. Biol.
209
,
2409
-
2419
.
Denny
,
M. W.
,
Dowd
,
W. W.
,
Bilir
,
L.
and
Mach
,
K. J.
(
2011
).
Spreading the risk: Small-scale body temperature variation among intertidal organisms and its implications for species persistence
.
J. Exp. Mar. Biol. Ecol.
400
,
175
-
190
.
Dong
,
Y.
and
Somero
,
G. N.
(
2009
).
Temperature adaptation of cytosolic malate dehydrogenases of limpets (genus Lottia): differences in stability and function due to minor changes in sequence correlate with biogeographic and vertical distributions
.
J. Exp. Biol.
212
,
169
-
177
.
Dong
,
Y.
,
Miller
,
L. P.
,
Sanders
,
J. G.
and
Somero
,
G. N.
(
2008
).
Heat-shock protein 70 (Hsp70) expression in four limpets of the genus Lottia: interspecific variation in constitutive and inducible synthesis correlates with in situ exposure to heat stress
.
Biol. Bull.
215
,
173
-
181
.
Dong
,
Y.-W.
,
Liao
,
M.-L.
,
Han
,
G.-D.
and
Somero
,
G. N.
(
2022
).
An integrated, multi-level analysis of thermal effects on intertidal molluscs for understanding species distribution patterns
.
Biol. Rev.
97
,
554
-
581
.
Dowd
,
W. W.
(
2012
).
Challenges for biological interpretation of environmental proteomics data in non-model organisms
.
Integr. Comp. Biol.
52
,
705
-
720
.
Dowd
,
W. W.
and
Somero
,
G. N.
(
2013
).
Behavior and survival of Mytilus congeners following episodes of elevated body temperature in air and seawater
.
J. Exp. Biol.
216
,
502
-
514
.
Dowd
,
W. W.
,
Harris
B. N.
,
Cech
J. J.
Jr
and
Kültz
D.
(
2010a
).
Proteomic and physiological responses of leopard sharks (Triakis semifasciata) to salinity change
.
J. Exp. Biol.
213
,
210
-
224
.
Dowd
,
W. W.
,
Renshaw
,
G. M. C.
,
Cech
,
J. J.
Jr
and
Kültz
,
D.
(
2010b
).
Compensatory proteome adjustments imply tissue-specific structural and metabolic reorganization following episodic hypoxia or anoxia in the epaulette shark
.
Physiol. Genomics
42
,
93
-
114
.
Dowd
,
W. W.
,
Felton
,
C. A.
,
Heymann
,
H. M.
,
Kost
,
L. E.
and
Somero
,
G. N.
(
2013
).
Food availability, more than body temperature, drives correlated shifts in ATP-generating and antioxidant enzyme capacities in a population of intertidal mussels (Mytilus californianus)
.
J. Exp. Mar. Biol. Ecol.
449
,
171
-
185
.
Drake
,
M. J.
,
Miller
,
N. A.
and
Todgham
,
A. E.
(
2017
).
The role of stochastic thermal environments in modulating the thermal physiology of an intertidal limpet, Lottia digitalis
.
J. Exp. Biol.
220
,
3072
-
3083
.
Dutton
,
J. M.
and
Hofmann
,
G. E.
(
2008
).
Spatial and temporal variation in distribution and protein ubiquitination for Mytilus congeners in the California hybrid zone
.
Mar. Biol.
154
,
1067
-
1075
.
Elowe
,
C.
and
Tomanek
,
L.
(
2021
).
Circadian and circatidal rhythms of protein abundance in the California mussel (Mytilus californianus)
.
Mol. Ecol.
30
,
5151
-
5163
.
Feidantsis
,
K.
,
Giantsis
,
I. A.
,
Vratsistas
,
A.
,
Makri
,
S.
,
Pappa
,
A.-Z.
,
Drosopoulou
,
E.
,
Anestis
,
A.
,
Mavridou
,
E.
,
Exadactylos
,
A.
,
Vafidis
,
D.
et al. 
(
2020
).
Correlation between intermediary metabolism, Hsp gene expression, and oxidative stress-related proteins in long-term thermal-stressed Mytilus galloprovincialis
.
Am. J. Physiol. Regul. Integr. Comp. Physiol.
319
,
R264
-
R281
.
Fields
,
P. A.
,
Rudomin
,
E. L.
and
Somero
,
G. N.
(
2006
).
Temperature sensitivities of cytosolic malate dehydrogenases from native and invasive species of marine mussels (genus Mytilus): sequence-function linkages and correlations with biogeographic distribution
.
J. Exp. Biol.
209
,
656
-
667
.
Fields
,
P. A.
,
Zuzow
,
M. J.
and
Tomanek
,
L.
(
2012
).
Proteomic responses of blue mussel (Mytilus) congeners to temperature acclimation
.
J. Exp. Biol.
215
,
1106
-
1116
.
Garland
,
T.
Jr
and
Adolph
,
S. C.
(
1994
).
Why not to do two-species comparative studies: Limitations on inferring adaptation
.
Physiol. Zool.
67
,
797
-
828
.
Georgoulis
,
I.
,
Feidantsis
,
K.
,
Giantsis
,
I. A.
,
Kakale
,
A.
,
Bock
,
C.
,
Pörtner
,
H. O.
,
Sokolova
,
I. M.
and
Michaelidis
,
B.
(
2021
).
Heat hardening enhances mitochondrial potential for respiration and oxidative defence capacity in the mantle of thermally stressed Mytilus galloprovincialis
.
Sci. Rep.
11
,
17098
.
Gleason
,
L. U.
,
Miller
,
L. P.
,
Winnikoff
,
J. R.
,
Somero
,
G. N.
,
Yancey
,
P. H.
,
Bratz
,
D.
and
Dowd
,
W. W.
(
2017
).
Thermal history and gape of individual Mytilus californianus correlate with oxidative damage and thermoprotective osmolytes
.
J. Exp. Biol.
220
,
4292
-
4304
.
Gracey
,
A. Y.
,
Chaney
,
M. L.
,
Boomhower
,
J. P.
,
Tyburczy
,
W. R.
,
Connor
,
K.
and
Somero
,
G. N.
(
2008
).
Rhythms of gene expression in a fluctuating intertidal environment
.
Curr. Biol.
18
,
1501
-
1507
.
Grimsrud
,
P. A.
,
Picklo
,
M. J.
, Sr
,
Griffin
,
T. J.
and
Bernlohr
,
D. A.
(
2007
).
Carbonylation of adipose proteins in obesity and insulin resistance: identification of adipocyte fatty acid-binding protein as a cellular target of 4-hydroxynonenal
.
Mol. Cell. Proteomics
6
,
624
-
637
.
Grune
,
T.
,
Reinheckel
,
T.
and
Davies
,
K. J. A.
(
1997
).
Degradation of oxidized proteins in mammalian cells
.
FASEB J.
11
,
526
-
534
.
Halliwell
,
B.
and
Gutteridge
,
J. M. C.
(
2007
).
Free Radicals in Biology and Medicine
.
Oxford
:
Oxford University Press
.
Harley
,
C. D. G.
(
2008
).
Tidal dynamics, topographic orientation, and temperature-mediated mass mortalities on rocky shores
.
Mar. Ecol. Progr. Ser.
371
,
37
-
46
.
Harris
,
R. M. B.
,
Beaumont
,
L. J.
,
Vance
,
T. R.
,
Tozer
,
C. R.
,
Remenyi
,
T. A.
,
Perkins-Kirkpatrick
,
S. E.
,
Mitchell
,
P. J.
,
Nicotra
,
A. B.
,
Mcgregor
,
S.
,
Andrew
,
N. R.
et al. 
(
2018
).
Biological responses to the press and pulse of climate trends and extreme events
.
Nat. Clim. Change
8
,
579
-
587
.
Heath
,
D. D.
,
Rawson
,
P. D.
and
Hilbish
,
T. J.
(
1995
).
PCR-based nuclear markers identify alien blue mussel (Mytilus spp.) genotypes on the west coast of Canada
.
Can. J. Fish. Aquat. Sci.
52
,
2621
-
2627
.
Hecker
,
M.
and
Wagner
,
A. H.
(
2018
).
Role of protein carbonylation in diabetes
.
J. Inherit. Metab. Dis.
41
,
29
-
38
.
Heise
,
K.
,
Puntarulo
,
S.
,
Pörtner
,
H. O.
and
Abele
,
D.
(
2003
).
Production of reactive oxygen species by isolated mitochondria of the Antarctic bivalve Laternula elliptica (King and Broderip) under heat stress
.
Comp. Biochem. Physiol. C Toxicol. Pharmacol.
134
,
79
-
90
.
Heise
,
K.
,
Puntarulo
,
S.
,
Nikinmaa
,
M.
,
Abele
,
D.
and
Pörtner
,
H.-O.
(
2006
).
Oxidative stress during stressful heat exposure and recovery in the North Sea eelpout Zoarces viviparus L
.
J. Exp. Biol.
209
,
353
-
363
.
Helmuth
,
B. S. T.
(
1998
).
Intertidal mussel microclimates: predicting the body temperature of a sessile invertebrate
.
Ecol. Monogr.
68
,
51
-
74
.
Hermes-Lima
,
M.
and
Storey
,
K. B.
(
1995
).
Antioxidant defenses and metabolic depression in a pulmonate land snail
.
Am. J. Physiol.
268R
,
1386
-
1393
.
Hertzberg
,
S.
,
Partali
,
V.
and
Liaaen-Jensen
,
S.
(
1988
).
Animal Carotenoids. 32. Carotenoids of Mytilus edulis (Edible Mussel)
.
Acta Chem. Scand.
B42
,
495
-
503
.
Hilbish
,
T. J.
,
Brannock
,
P. M.
,
Jones
,
K. R.
,
Smith
,
A. B.
,
Bullock
,
B. N.
and
Wethey
,
D. S.
(
2010
).
Historical changes in the distributions of invasive and endemic marine invertebrates are contrary to global warming predictions: the effects of decadal climate oscillations
.
J. Biogeogr.
37
,
423
-
431
.
Hofmann
,
G. E.
and
Somero
,
G. N.
(
1995
).
Evidence for protein damage at environmental temperatures: seasonal changes in levels of ubiquitin conjugates and hsp70 in the intertidal mussel Mytilus trossulus
.
J. Exp. Biol.
198
,
1509
-
1518
.
Hofmann
,
G. E.
and
Somero
,
G. N.
(
1996
).
Interspecific variation in thermal denaturation of proteins in the congeneric mussels Mytilus trossulus and M. galloprovincialis: evidence from the heat-shock response and protein ubiquitination
.
Mar. Biol.
126
,
65
-
75
.
Janssens
,
L.
and
Stoks
,
R.
(
2018
).
Rapid larval development under time stress reduces adult life span through increasing oxidative damage
.
Func. Ecol.
32
,
1036
-
1045
.
Jimenez
,
A. G.
,
Alves
,
S.
,
Dallmer
,
J.
,
Njoo
,
E.
,
Roa
,
S.
and
Dowd
,
W. W.
(
2016
).
Acclimation to elevated emersion temperature has no effect on susceptibility to acute, heat-induced lipid peroxidation in an intertidal mussel (Mytilus californianus)
.
Mar. Biol.
163
,
55
.
Jørgensen
,
C. B.
(
1974
).
On gill function in the mussel Mytilus edulis L
.
Ophelia
13
,
187
-
232
.
Kassahn
,
K. S.
,
Crozier
,
R. H.
,
Pörtner
,
H. O.
and
Caley
,
M. J.
(
2009
).
Animal performance and stress: responses and tolerance limits at different levels of biological organisation
.
Biol. Rev.
84
,
277
-
292
.
Koch
,
R. E.
,
Buchanan
,
K. L.
,
Casagrande
,
S.
,
Crino
,
O.
,
Dowling
,
D. K.
,
Hill
,
G. E.
,
Hood
,
W. R.
,
Mckenzie
,
M.
,
Mariette
,
M. M.
,
Noble
,
D. W. A.
et al. 
(
2021
).
Integrating mitochondrial aerobic metabolism into ecology and evolution
.
Trends Ecol. Evol.
36
,
321
-
332
.
Koehn
,
R. K.
(
1991
).
The genetics and taxonomy of species in the genus Mytilus
.
Aquaculture
94
,
125
-
145
.
Korolainen
,
M. A.
,
Nyman
,
T. A.
,
Nyyssönen
,
P.
,
Hartikainen
,
E. S.
and
Pirttilä
,
T.
(
2007
).
Multiplexed proteomic analysis of oxidation and concentrations of cerebrospinal fluid proteins in Alzheimer disease
.
Clin. Chem.
53
,
657
-
665
.
Li
,
C.
and
Jackson
,
R. M.
(
2002
).
Reactive species mechanisms of cellular hypoxia-reoxygenation injury
.
Am. J. Physiol.
282
,
C227
-
C241
.
Liao
,
M.-L.
,
Dong
,
Y.-W.
and
Somero
,
G. N.
(
2021
).
Thermal adaptation of mRNA secondary structure: stability versus lability
.
Proc. Natl. Acad. Sci. USA
118
,
e2113324118
.
Liu
,
H.
,
Uno
,
M.
,
Kitazato
,
K. T.
,
Suzue
,
A.
,
Manabe
,
S.
,
Yamasaki
,
H.
,
Shono
,
M.
and
Nagahiro
,
S.
(
2004
).
Peripheral oxidative biomarkers constitute a valuable indicator of the severity of oxidative brain damage in acute cerebral infarction
.
Brain Res.
1025
,
43
-
50
.
Lockwood
,
B. L.
and
Somero
,
G. N.
(
2011
).
Invasive and native blue mussels (genus Mytilus) on the California coast: the role of physiology in a biological invasion
.
J. Exp. Mar. Biol. Ecol.
400
,
167
-
174
.
Lockwood
,
B. L.
and
Somero
,
G. N.
(
2012
).
Functional determinants of temperature adaptation in enzymes of cold- versus warm-adapted mussels (Genus Mytilus)
.
Mol. Biol. Evol.
29
,
3061
-
3070
.
Lockwood
,
B. L.
,
Sanders
,
J. G.
and
Somero
,
G. N.
(
2010
).
Transcriptomic responses to heat stress in invasive and native blue mussels (genus Mytilus): molecular correlates of invasive success
.
J. Exp. Biol.
213
,
3548
-
3558
.
Logan
,
C. A.
,
Kost
,
L. E.
and
Somero
,
G. N.
(
2012
).
Latitudinal differences in Mytilus californianus thermal physiology
.
Mar. Ecol. Prog. Ser.
450
,
93
-
105
.
Luo
,
S.
and
Wehr
,
N. B.
(
2009
).
Protein carbonylation: avoiding pitfalls in the 2,4-dinitrophenylhydrazine assay
.
Redox Rep.
14
,
159
-
166
.
Lushchak
,
V. I.
,
Bagnyukova
,
T. V.
,
Lushchak
,
O. V.
,
Storey
,
J. M.
and
Storey
,
K. B.
(
2005
).
Hypoxia and recovery perturb free radical processes and antioxidant potential in common carp (Cyprinus carpio) tissues
.
Int. J. Biochem. Cell Biol.
37
,
1319
-
1330
.
Maisonneuve
,
E.
,
Ducret
,
A.
,
Khoueiry
,
P.
,
Lignon
,
S.
,
Longhi
,
S.
,
Talla
,
E.
and
Dukan
,
S.
(
2009
).
Rules governing selective protein carbonylation
.
PLoS ONE
4
,
e7269
.
McCord
,
J. M.
and
Fridovich
,
I.
(
1969
).
Superoxide dismutase: an enzymic function for erythrocuprein (hemocuprein)
.
J. Biol. Chem
244
,
6049
-
6055
.
McDonagh
,
B.
and
Sheehan
,
D.
(
2006
).
Redox proteomics in the blue mussel Mytilus edulis: carbonylation is not a pre-requisite for ubiquitination in acute free radical-mediated oxidative stress
.
Aquat. Toxic.
79
,
325
-
333
.
McDonagh
,
B.
,
Tyther
,
R.
and
Sheehan
,
D.
(
2005
).
Carbonylation and glutathionylation of proteins in the blue mussel Mytilus edulis detected by proteomic analysis and Western blotting: Actin as a target for oxidative stress
.
Aquat. Toxic.
73
,
315
-
326
.
McDonald
,
J. H.
and
Koehn
,
R. K.
(
1988
).
The mussels Mytilus galloprovincialis and M. trossulus on the Pacific coast of North America
.
Mar. Biol.
99
,
111
-
118
.
McDonald
,
J. H.
,
Seed
,
R.
and
Koehn
,
R. K.
(
1991
).
Allozymes and morphometric characters of three species of Mytilus in the Northern and Southern Hemispheres
.
Mar. Biol.
111
,
323
-
333
.
Miller
,
L. P.
and
Dowd
,
W. W.
(
2017
).
Multimodal in situ datalogging quantifies inter-individual variation in thermal experience and persistent origin effects on gaping behavior among intertidal mussels (Mytilus californianus)
.
J. Exp. Biol.
220
,
4305
-
4319
.
Mislan
,
K. A. S.
and
Wethey
,
D. S.
(
2015
).
A biophysical basis for patchy mortality during heat waves
.
Ecology
96
,
902
-
907
.
Monaghan
,
P.
,
Metcalfe
,
N. B.
and
Torres
,
R.
(
2009
).
Oxidative stress as a mediator of life history trade–offs: mechanisms, measurements and interpretation
.
Ecol. Lett.
12
,
75
-
92
.
Moskovitz
,
J.
and
Oien
,
D. B.
(
2010
).
Protein carbonyl and the methionine sulfoxide reductase system
.
Antioxid. Redox Signal.
12
,
405
-
415
.
Mueller
,
I.
,
Hoffman
,
M.
,
Dullen
,
K.
and
O'brien
,
K.
(
2014
).
Moderate elevations in temperature do not increase oxidative stress in oxidative muscles of Antarctic notothenioid fishes
.
Polar Biol.
37
,
311
-
320
.
O'Brien
,
K. M.
,
Dirmeier
,
R.
,
Engle
,
M.
and
Poyton
,
R. O.
(
2004
).
Mitochondrial protein oxidation in yeast mutants lacking manganese-(MnSOD) or copper- and zinc-containing superoxide dismutase (CuZnSOD)
.
J. Biol. Chem.
279
,
51817
-
51827
.
Perkins
,
D. N.
,
Pappin
,
D. J. C.
,
Creasy
,
D. M.
and
Cottrell
,
J. S.
(
1999
).
Probability-based protein identification by searching sequence databases using mass spectrometry data
.
Electrophoresis
20
,
3551
-
3567
.
Pierce
,
A.
,
Mirzaei
,
H.
,
Muller
,
F.
,
De Waal
,
E.
,
Taylor
,
A. B.
,
Leonard
,
S.
,
Van Remmen
,
H.
,
Regnier
,
F.
,
Richardson
,
A.
and
Chaudhuri
,
A.
(
2008
).
GAPDH is conformationally and functionally altered in association with oxidative stress in mouse models of amyotrophic lateral sclerosis
.
J. Mol. Biol.
382
,
1195
-
1210
.
Pikul
,
J.
,
Leszczynski
,
D. E.
and
Kummerow
,
F. A.
(
1983
).
Elimination of sample autoxidation by butylated hydroxytoluene additions before thiobarbituric acid assay for malonaldehyde in fat from chicken meat
.
J. Agric. Food Chem.
31
,
1338
-
1342
.
Popovic
,
I.
and
Riginos
,
C.
(
2020
).
Comparative genomics reveals divergent thermal selection in warm- and cold-tolerant marine mussels
.
Mol. Ecol.
29
,
519
-
535
.
Pörtner
,
H. O.
(
2002
).
Climate variations and the physiological basis of temperature dependent biogeography: systemic to molecular hierarchy of thermal tolerance in animals
.
Comp. Biochem. Physiol. A Mol. Integr. Physiol.
132
,
739
-
761
.
Pörtner
,
H.-O.
(
2010
).
Oxygen- and capacity-limitation of thermal tolerance: a matrix for integrating climate-related stressor effects in marine ecosystems
.
J. Exp. Biol.
213
,
881
-
893
.
Rabek
,
J. P.
,
Boylston Iii
,
W. H.
and
Papaconstantinou
,
J.
(
2003
).
Carbonylation of ER chaperone proteins in aged mouse liver
.
Biochem. Biophys. Res. Commun.
305
,
566
-
572
.
Rawson
,
P. D.
,
Joyner
,
K. L.
,
Meetze
,
K.
and
Hilbish
,
T. J.
(
1996
).
Evidence for intragenic recombination within a novel genetic marker that distinguishes mussels in the Mytilus edulis species complex
.
Heredity
77
,
599
-
607
.
Reinheckel
,
T.
,
Körn
,
S.
,
Möhring
,
S.
,
Augustin
,
W.
,
Halangk
,
W.
and
Schild
,
L.
(
2000
).
Adaptation of protein carbonyl detection to the requirements of proteome analysis demonstrated for hypoxia/reoxygenation in isolated rat liver mitochondria
.
Arch. Biochem. Biophys.
376
,
59
-
65
.
Rivera-Ingraham
,
G. A.
,
Rocchetta
,
I.
,
Meyer
,
S.
and
Abele
,
D.
(
2013
).
Oxygen radical formation in anoxic transgression and anoxia-reoxygenation: Foe or phantom? Experiments with a hypoxia tolerant bivalve
.
Mar. Envir. Res.
92
,
110
-
119
.
Robinson
,
T. B.
and
Griffiths
,
C. L.
(
2002
).
Invasion of Langebaan Lagoon, South Africa, by Mytilus galloprovincialis – effects on natural communities
.
African Zool.
37
,
151
-
158
.
Schneider
,
K. R.
(
2008
).
Heat stress in the intertidal: comparing survival and growth of an invasive and native mussel under a variety of thermal conditions
.
Biol. Bull.
215
,
253
-
264
.
Schneider
,
K. R.
and
Helmuth
,
B.
(
2007
).
Spatial variability in habitat temperature may drive patterns of selection between an invasive and native mussel species
.
Mar. Ecol. Progr. Ser.
339
,
157
-
167
.
Seed
,
R.
(
1992
).
Systematics, evolution and distribution of mussels belonging to the genus Mytilus: an overview
.
Am. Malacol. Bull.
9
,
123
-
137
.
Seuront
,
L.
,
Nicastro
,
K. R.
,
Zardi
,
G. I.
and
Goberville
,
E.
(
2019
).
Decreased thermal tolerance under recurrent heat stress conditions explains summer mass mortality of the blue mussel Mytilus edulis
.
Sci. Rep.
9
,
17498
.
Shanmuganathan
,
A.
,
Avery
,
S. V.
,
Willetts
,
S. A.
and
Houghton
,
J. E.
(
2004
).
Copper-induced oxidative stress in Saccharomyces cerevisiae targets enzymes of the glycolytic pathway
.
FEBS Lett.
556
,
253
-
259
.
Shick
,
J. M.
,
Widdows
,
J.
and
Gnaiger
,
E.
(
1988
).
Calorimetric studies of behavior, metabolism and energetics of sessile intertidal animals
.
Amer. Zool.
28
,
161
-
181
.
Sokolov
,
E. P.
,
Adzigbli
,
L.
,
Markert
,
S.
,
Bundgaard
,
A.
,
Fago
,
A.
,
Becher
,
D.
,
Hirschfeld
,
C.
and
Sokolova
,
I. M.
(
2021
).
Intrinsic mechanisms underlying hypoxia-tolerant mitochondrial phenotype during hypoxia-reoxygenation stress in a marine facultative anaerobe, the blue mussel Mytilus edulis
.
Front. Mar. Sci.
8
.
Somero
,
G. N.
,
Lockwood
,
B. L.
and
Tomanek
,
L.
(
2017
).
Biochemical Adaptation: Response to Environmental Challenges from Life's Origins to the Anthropocene
.
Sunderland, MA
:
Sinauer Associates, Inc
.
Stadtman
,
E. R.
(
1993
).
Oxidation of free amino acids and amino acid residues in proteins by radiolysis and by metal-catalyzed reactions
.
Ann. Rev. Biochem.
62
,
797
-
821
.
Stadtman
,
E. R.
and
Levine
,
R. L.
(
2003
).
Free radical-mediated oxidation of free amino acids and amino acid residues in proteins
.
Amino Acids
25
,
207
-
218
.
Stillman
,
J. H.
(
2019
).
Heat waves, the new normal: Summertime temperature extremes will impact animals, ecosystems, and human communities
.
Physiology
34
,
86
-
100
.
Stillman
,
J. H.
and
Somero
,
G. N.
(
2000
).
A comparative analysis of the upper thermal tolerance limits of Eastern Pacific porcelain crabs, genus Petrolisthes: Influences of latitude, vertical zonation, acclimation, and phylogeny
.
Physiol. Biochem. Zool.
73
,
200
-
208
.
Storey
,
J. D.
and
Tibshirani
,
R.
(
2003
).
Statistical significance for genomewide studies
.
Proc. Natl. Acad. Sci. USA
100
,
9440
-
9445
.
Sultana
,
R.
,
Perluigi
,
M.
,
Newman
,
S. F.
,
Pierce
,
W. M.
,
Cini
,
C.
,
Coccia
,
R.
and
Butterfield
,
D. A.
(
2010
).
Redox proteomic analysis of carbonylated brain proteins in mild cognitive impairment and early Alzheimer's disease
.
Antioxid. Redox Signal.
12
,
327
-
336
.
Sun
,
J.-L.
,
Zhao
,
L.-L.
,
Liao
,
L.
,
Tang
,
X.-H.
,
Cui
,
C.
,
Liu
,
Q.
,
He
,
K.
,
Ma
,
J.-D.
,
Jin
,
L.
,
Yan
,
T.
et al. 
(
2020
).
Interactive effect of thermal and hypoxia on largemouth bass (Micropterus salmoides) gill and liver: Aggravation of oxidative stress, inhibition of immunity and promotion of cell apoptosis
.
Fish Shellfish Immunol.
98
,
923
-
936
.
Thomsen
,
J.
,
Himmerkus
,
N.
,
Holland
,
N.
,
Sartoris
,
F. J.
,
Bleich
,
M.
and
Tresguerres
,
M.
(
2016
).
Ammonia excretion in mytilid mussels is facilitated by ciliary beating
.
J. Exp. Biol.
219
,
2300
-
2310
.
Tomanek
,
L.
(
2011
).
Environmental proteomics: Changes in the proteome of marine organisms in response to environmental stress, pollutants, infection, symbiosis, and development
.
Annu. Rev. Mar. Sci.
3
,
14.1
-
14.27
.
Tomanek
,
L.
and
Helmuth
,
B.
(
2002
).
Physiological ecology of rocky intertidal organisms: a synergy of concepts
.
Integr. Comp. Biol.
42
,
771
-
775
.
Tomanek
,
L.
and
Somero
,
G. N.
(
2000
).
Time course and magnitude of synthesis of heat-shock proteins in congeneric marine snails (Genus Tegula) from different tidal heights
.
Physiol. Biochem. Zool.
73
,
249
-
256
.
Tomanek
,
L.
and
Zuzow
,
M. J.
(
2010
).
The proteomic response of the mussel congeners Mytilus galloprovincialis and M. trossulus to acute heat stress: implications for thermal tolerance limits and metabolic costs of thermal stress
.
J. Exp. Biol.
213
,
3559
-
3574
.
Wang
,
J.
,
Ren
,
R.-M.
and
Yao
,
C.-L.
(
2018
).
Oxidative stress responses of Mytilus galloprovincialis to acute cold and heat during air exposure
.
J. Mollus. Stud.
84
,
285
-
292
.
Welker
,
A. F.
,
Moreira
,
D. C.
,
Campos
,
É. G.
and
Hermes-Lima
,
M.
(
2013
).
Role of redox metabolism for adaptation of aquatic animals to drastic changes in oxygen availability
.
Comp. Biochem. Physiol. A
165
,
384
-
404
.
Widdows
,
J.
and
Shick
,
J. M.
(
1985
).
Physiological responses of Mytilus edulis and Cardium edule to aerial exposure
.
Mar. Biol.
85
,
217
-
232
.
Wilkins
,
N. P.
,
Fujino
,
K.
and
Gosling
,
E. M.
(
2008
).
The mediterranean mussel Mytilus galloprovincialis Lmk. in Japan
.
Biol. J. Linn. Soc.
20
,
365
-
374
.
Williams
,
C. M.
,
Buckley
,
L. B.
,
Sheldon
,
K. S.
,
Vickers
,
M.
,
Pörtner
,
H.-O.
,
Dowd
,
W. W.
,
Gunderson
,
A. R.
,
Marshall
,
K. E.
and
Stillman
,
J. H.
(
2016
).
Biological impacts of thermal extremes: Mechanisms and costs of functional responses matter
.
Integr. Comp. Biol.
56
,
73
-
84
.
Wong
,
C. M.
,
Cheema
,
A. K.
,
Zhang
,
L.
and
Suzuki
,
Y. J.
(
2008
).
Protein carbonylation as a novel mechanism in redox signaling
.
Circ. Res.
102
,
310
-
318
.
Wong
,
C.-M.
,
Marcocci
,
L.
,
Das
,
D.
,
Wang
,
X.
,
Luo
,
H.
,
Zungu-Edmondson
,
M.
and
Suzuki
,
Y. J.
(
2013
).
Mechanism of protein decarbonylation
.
Free Rad. Biol. Med.
65
,
1126
-
1133
.
Wonham
,
M. J.
(
2004
).
Mini-review: Distribution of the Mediterranean mussel Mytilus galloprovincialis (Bivalvia: Mytilidae) and hybrids in the northeast Pacific
.
J. Shellfish Res.
23
,
535
-
543
.

Competing interests

The authors declare no competing or financial interests.

Supplementary information