Intrinsic Mechanisms Underlying Hypoxia-Tolerant Mitochondrial Phenotype During Hypoxia-Reoxygenation Stress in a Marine Facultative Anaerobe, the Blue Mussel Mytilus edulis

Hypoxia is common in marine environments and a major stressor for marine organisms inhabiting benthic and intertidal zones. Several studies have explored the responses of these organisms to hypoxic stress at the whole organism level with a focus on energy metabolism and mitochondrial response, but the instrinsic mitochondrial responses that support the organelle’s function under hypoxia and reoxygenation (H/R) stress are not well understood. We studied the effects of acute H/R stress (10 min anoxia followed by 15 min reoxygenation) on mitochondrial respiration, production of reactive oxygen species (ROS) and posttranslational modifications (PTM) of the proteome in a marine facultative anaerobe, the blue mussel Mytilus edulis. The mussels’ mitochondria showed increased OXPHOS respiration and suppressed proton leak resulting in a higher coupling efficiency after H/R stress. ROS production decreased in both the resting (LEAK) and phosphorylating (OXPHOS) state indicating that M. edulis was able to prevent oxidative stress and mitochondrial damage during reoxygenation. Hypoxia did not lead to rearrangement of the mitochondrial supercomplexes but impacted the mitochondrial phosphoproteome including the proteins involved in OXPHOS, amino acid- and fatty acid catabolism, and protein quality control. This study indicates that mussels’ mitochondria possess intrinsic mechanisms (including regulation via reversible protein phosphorylation) that ensure high respiratory flux and mitigate oxidative damage during H/R stress and contribute to the hypoxia-tolerant mitochondrial phenotype of this metabolically plastic species.

Hypoxia is common in marine environments and a major stressor for marine organisms inhabiting benthic and intertidal zones. Several studies have explored the responses of these organisms to hypoxic stress at the whole organism level with a focus on energy metabolism and mitochondrial response, but the instrinsic mitochondrial responses that support the organelle's function under hypoxia and reoxygenation (H/R) stress are not well understood. We studied the effects of acute H/R stress (10 min anoxia followed by 15 min reoxygenation) on mitochondrial respiration, production of reactive oxygen species (ROS) and posttranslational modifications (PTM) of the proteome in a marine facultative anaerobe, the blue mussel Mytilus edulis. The mussels' mitochondria showed increased OXPHOS respiration and suppressed proton leak resulting in a higher coupling efficiency after H/R stress. ROS production decreased in both the resting (LEAK) and phosphorylating (OXPHOS) state indicating that M. edulis was able to prevent oxidative stress and mitochondrial damage during reoxygenation. Hypoxia did not lead to rearrangement of the mitochondrial supercomplexes but impacted the mitochondrial phosphoproteome including the proteins involved in OXPHOS, amino acid-and fatty acid catabolism, and protein quality control. This study indicates that mussels' mitochondria possess intrinsic mechanisms (including regulation via reversible protein phosphorylation) that ensure high respiratory flux and mitigate oxidative damage during H/R stress and contribute to the hypoxia-tolerant mitochondrial phenotype of this metabolically plastic species.

INTRODUCTION
Oxygen deficiency (hypoxia) is a major stressor for aerobic eukaryotes that depend on oxygen to drive mitochondrial ATP synthesis. For these organisms, oxygen fluctuations are stressful due to the limitations of energy supply during hypoxia, as well as the oxidative stress and high costs of recovery during reoxygenation (Garbarino et al., 2015;Lesnefsky et al., 2017;Sokolova et al., 2019). Oxygen is rarely limiting in terrestrial habitats (except underground burrows and high altitudes) but in aquatic environments oxygen fluctuations are common. Many benthic marine organisms -particularly sessile species from the habitats where oxygen availability varies regularly and predictably such as the intertidal zone -have adapted to survive under the fluctuating oxygen conditions. The ecological success of these species is often linked with their ability to transition to anaerobiosis and enter the metabolically depressed state during hypoxia (Hochachka and Mustafa, 1972;Guppy et al., 1994;Hochachka and Somero, 2002). The metabolic rate depression in facultative anaerobes involves complex regulatory mechanisms that ensure a concomitant decrease in ATP consumption and production, onset of anaerobic glycolysis, and changes in membrane properties to maintain ion homeostasis with minimum energy input (Hochachka and Mustafa, 1972;Hochachka and Mommsen, 1983;Hochachka and Guppy, 1987;Buck and Hochachka, 1993;Guppy et al., 1994;Storey and Storey, 2004;Storey, 2015). However, the metabolically arrested state allows only time-limited survival, and to continue the life cycle, the organisms depend on the return of oxygen. Therefore, tolerance to frequent oxygen fluctuations requires maintenance of the mitochondrial integrity and rapid restoration of the energetic and redox homeostasis during reoxygenation (Honda et al., 2005;Pamenter, 2014;Pell et al., 2016;Sokolova et al., 2019).
Recent studies indicate that mitochondria of hypoxia-tolerant organisms are capable of preserving the mitochondrial integrity and function during hypoxia-reoxygenation (H/R) stress via rapid suppression and restoration of oxidative phosphorylation (OXPHOS) subject to oxygen availability, mitochondrial proteome reorganization to limit production of reactive oxygen species (ROS), and upregulation of mitochondrial quality control pathways (Buck and Pamenter, 2006;Pamenter, 2014;Pamenter et al., 2016;Bundgaard et al., 2019;Steffen et al., 2020). These changes are observed during long-term (hours to days) exposures of the organisms to hypoxia and subsequent reoxygenation in vivo (Kurochkin et al., 2008;Ivanina et al., 2010Ivanina et al., , 2011Ivanina et al., , 2012Ivanina et al., , 2016Pamenter et al., 2016;Steffen et al., 2020;Ouillon et al., 2021) and likely reflect a complex integrated response encompassing the direct impact of the oxygen fluctuations on the mitochondria, indirect effects of the changes in intracellular milieu, responses driven by retrograde signaling to the nucleus, as well as endocrine and/or neural regulation (Lee et al., 2020;Du et al., 2021). This complexity, while physiologically relevant and important, makes it difficult to disentangle the contributions of different mechanisms to the functional changes of the mitochondria during H/R stress. Therefore, a reductionist approach using a simplified and thus more tractable system (such as the isolated mitochondria) could provide useful insights. Mitochondria are semi-autonomous organelles that possess their own DNA and protein translation machinery, as well as a dedicated system of quality control and post-translational regulatory mechanisms (McBride et al., 2006) making them capable of direct responses to H/R transitions independently of the rest of the cell. The functional changes during acute H/R exposures of isolated mitochondria provide evidence for such intrinsic mitochondrial regulation (Onukwufor et al., 2014(Onukwufor et al., , 2016(Onukwufor et al., , 2017Sappal et al., 2015), but the molecular mechanisms of this regulation have not yet been investigated.
The intrinsic mitochondrial responses to H/R transitions must rely on post-translational regulatory mechanisms such as the changing levels of the metabolic intermediates, shifts in the state of the electron carriers of the ETS, and posttranslational modifications (PTM) of mitochondrial proteins. Among the PTM mechanisms, assembly of the respiratory supercomplexes and reversible protein phosphorylation (catalyzed by mitochondrial phosphatases and kinases) are proposed to play a key role in the modulation of mitochondrial metabolism (Pagliarini and Dixon, 2006;Lapuente-Brun et al., 2013;Chaban et al., 2014;Mathers and Staples, 2019). Mitochondrial supercomplexes (SC) involve a close association between complexes CI, CIII, and CIV as a potential mechanism to channel electron flux, increase the efficiency of OXPHOS and mitigate ROS generation (Guerra-Castellano et al., 2018). The functional consequences of SC rearrangement are debated, but recent studies showed links between SC stability and plasticity and metabolic responses to stress in animals (Bundgaard et al., 2018(Bundgaard et al., , 2020aFalfushynska et al., 2020). Furthermore, reversible phosphorylation of the proteins of the mitochondrial complexes CI and CIV can contribute to the regulation of the mitochondrial function during H/R stress by modulating ETS activity, protein-protein and protein-lipid interactions, and ROS production (Helling et al., 2012;Gowthami et al., 2019;Kalpage et al., 2019;Mathers and Staples, 2019). Stress-induced shifts in phosphoproteome and modulation of the activity of protein phosphatases and kinases have been shown in marine mollusks exposed to H/R stress Falfushynska et al., 2020); however, these studies have been conducted during the whole-organism exposures to H/R stress, which does not allow discerning the intrinsic mitochondrial mechanisms of the observed changes. Overall, the potential contribution of the supercomplex rearrangement and reversible protein phosphorylation to hypoxia-tolerant mitochondrial phenotype has not been well studied in facultative anaerobes and warrants further investigation.
The mussels of the genus Mytilus are widespread facultative anaerobes with ecological and economical importance (Astorga, 2014). Their habitats (including the intertidal zone, eutrophicated estuaries, and coastal habitats) are characterized by recurrent oxygen fluctuations, hence they are highly adapted to these conditions (Gosling, 1992). The energy metabolism pathways of Mytilus spp. are well characterized (Gosling, 1992;Babarro et al., 2007;Connor and Gracey, 2011), and the genomes and transcriptomes of several Mytilus species have been recently published (Martino et al., 2019;Paterno et al., 2019;Li et al., 2020) making these organisms an excellent model to study the molecular and physiological mechanisms of mitochondrial responses to H/R stress. Our present study aimed to elucidate the potential autonomous regulatory mechanisms (i.e., those independent of the cellular metabolic context and regulatory mechanisms) of mitochondrial responses to H/R stress in the blue mussels Mytilus edulis. We hypothesized that isolated mitochondria of a hypoxia-tolerant facultative anaerobe, M. edulis, possess intrinsic mechanisms that permit maintaining the respiratory flux and ATP synthesis and minimize ROS production during acute H/R stress. We also hypothesized that posttranslational mechanisms such as the rearrangement of supercomplexes and reversible protein phosphorylation would be involved in the intrinsic mitochondrial response to H/R transitions. To test these hypotheses, we investigated the direct effects of hypoxia and reoxygenation on the respiration rate under phosphorylating conditions (measured as the OXPHOS activity), proton leak and ROS generation in the isolated mitochondria from the digestive gland of M. edulis. Furthermore, we determined the abundance of the supercomplexes (formed by associations between CI, CIII, and CIV) and shifts in the mitochondrial phosphoproteome in the isolated mitochondria of M. edulis exposed to H/R stress in vitro. The results of this study shed light on the metabolic plasticity of mitochondrial function and demonstrate autonomous mitochondrial mechanisms involved in adaptive metabolic responses to oxygen fluctuations in the hypoxiatolerant facultative anaerobe M. edulis.

Animal Collection and Maintenance
Experimental blue mussels Mytilus edulis were collected in Warnemünde harbor in the Baltic Sea near Rostock, Germany. The animals were acclimated at salinity 15 and 15 • C for a week before the experiment. The acclimation temperature and salinity were close to the habitat conditions at the time of collection. During the acclimation, the mussels were continuously fed ad libitum with Premium Reef Blend (CoralSands, Wiesbaden, Germany).

MitochondrialṀO 2 and ROS Production
Mitochondria were isolated from the mussels' digestive gland as described elsewhere (Haider et al., 2018;. The digestive gland was chosen as one of the most metabolically active organs and a main site of energy storage in mussels (Gosling, 1992). The osmolarity of isolation and assay buffers were optimized in pilot experiments to achieve high quality of mitochondrial isolations (assessed by the respiratory control ratio and cytochrome c test of mitochondrial integrity) for the studied M. edulis population (data not shown). Briefly, 0.5-1.2 g tissue was homogenized on ice in a glass-Teflon homogenizer containing 10 ml of isolation buffer (∼630 mOsm) containing 100 mM sucrose, 200 mM KCl, 50 mM NaCl, 30 mM HEPES pH 7.5, and 8 mM EGTA supplemented with protease inhibitors (1 mM phenylmethylsulfonyl fluoride (PMSF) and 2 µg/ml aprotinin). Homogenate was centrifuged at 4 • C and 2000 × g for 5 min, the supernatant was collected and centrifuged at 8000 × g for 8 min. The resulting mitochondrial pellet was resuspended in isolation buffer and kept on ice. All assays commenced within 5 min of mitochondrial isolation.
Mitochondrial respiration rate (ṀO 2 ) and ROS production were measured simultaneously in the high-resolution respirometer Oxygraph 2 k (Oroboros, Innsbruck, Austria) at 15 and 25 • C using the DatLab6 software. The assay temperatures were chosen to reflect the acclimation temperature (15 • C) typical for early summer at the study site (June average monthly surface water temperature ∼14.8 • C), and the elevated temperature close to the summer maximum (25 • C). Mitochondrial suspension was added to the instrument's chamber filled with 2 ml of assay buffer (∼525 mOsm) containing 30 mM HEPES pH 7.2, 185 mM sucrose, 130 mM KCl, 10 mM NaCl, buffer, 10 mM glucose, 1 mM MgCl 2 , 10 mM KH 2 PO 4 , and 1% fatty acid-free bovine serum albumin (BSA). The final mitochondrial protein concentrations in the assay chamber were ∼1-1.5 mg ml −1 . Because earlier studies in bivalves Haider et al., 2020) and mammals (Hansen et al., 2010;Wang et al., 2016) showed that a common intracellular osmolyte, taurine, might improve mitochondrial function, we tested the potential mitoprotective effect of taurine in the mussel mitochondria. For this, assay buffer was supplemented either with 50 mM taurine or 50 mM glycine. Glycine is the most abundant free amino acid in oyster tissues (Haider et al., 2020) and was used to distinguish taurine-specific effects on mitochondria (if any) from non-specific effects of amino acid addition. After addition of mitochondrial suspension and stabilization of the polarographic oxygen sensor, the following substrate-uncoupler-inhibitor titration (SUIT) protocol was used: 1) 5 mM pyruvate, 0.5 mM malate to spark Complex I-driven respiration (LEAK I); 2) 10 mM succinate to stimulate Complex II (LEAK I + II); 3) 2.5 mM ADP to stimulate OXPHOS; 4) 2.5 µM oligomycin to inhibit mitochondrial F O ,F 1 -ATPase followed by titration (in 0.5 µM steps) with the uncoupler carbonyl cyanide-chlorophenyl hydrazine (CCCP) to measure maximum ETS activity. ROS production was measured in parallel withṀO 2 using the Oxygraph-2K LED2 module with Fluorescent Sensor Green (525 nm) and Amplex Red (Thermo Fisher Scientific, Waltham, MA, United States) as described elsewhere (Ouillon et al., 2021). Addition of cytochrome C (0.8 µM) in the OXPHOS state (Gnaiger, 2012) increased mitochondrial respiration by <10% indicating integrity of isolated mitochondria (data not shown). The LEAK stateṀO 2 indicates the proton leak rate of resting (non-phosphorylating) mitochondria,ṀO 2 of ADP-stimulated mitochondria corresponds to the OXPHOS capacity, and CCCP-uncoupledṀO 2 reflects the maximum ETS capacity. Respiratory control ratio was calculated by dividing OXPHOS by LEAKI + II respiration (Estabrook, 1967).
To test for the effects of H/R stress on mitochondrial respiration and ROS production, steps 1-3 (up to and including ADP addition) of the above-described SUIT protocol were carried out to measure the pre-hypoxic (baseline) LEAK I + II and OXPHOS oxygen consumption and ROS production rates. The OXPHOS respiration continued until all oxygen in the chamber was consumed, and mitochondria were held in anoxia for 10 min. The chamber was then opened for reoxygenation (2.5-5 min to achieve >80% of air saturation), closed and the post-reoxygenation OXPHOS oxygen consumption rate was measured. Next, 2.5 µM of oligomycin was added to record oxygen consumption in the post-reoxygenation LEAK state (LEAK I + II). At the end of the assay, 0.5 µM of rotenone and 2.5 µM of antimycin A were added to assess non-mitochondrial respiration (<10% of the total oxygen consumption rate, data not shown). The effects of H/R stress on the oxygen consumption rates were tested at 15 and 25 • C, and ROS production at 25 • C due to low ROS generation rates at 15 • C.
The impact of H/R stress on mitochondrial ETS activity was measured in CCCP-uncoupled mitochondria at 15 • C using pyruvate (5 mM), succinate (10 mM), or the mixture of pyruvate and succinate as substrates. Malate (0.5 mM) was added to spark the respiration. Mitochondria were allowed to respire until all oxygen was consumed in the chamber, kept in anoxia for 10 min, and re-oxygenated to measure the post-hypoxic ETS activity.
After completion of the assays, an aliquot of each mitochondrial suspension was lysed in hypoosmotic solution using three freezing-thawing cycles followed by sonication, and mitochondrial protein content was measured according to Bradford using a commercial kit (Merck KGaA, Darmstadt, Germany) with BSA as a standard. Respiration rates were expressed in µmol O 2 min −1 g −1 protein and ROS production in µmol H 2 O 2 min −1 g −1 protein.

Isolation of Phosphoproteins for LC/MS Analysis
Isolation of mitochondrial phosphoproteins was done using a chromatography approach with the Pro-Q R Diamond Phosphoprotein Enrichment Kit (Invitrogen, Carlsbad, CA, United States) according to the manufacturer's recommendations. Briefly, 40-60 mg of mitochondrial suspension was lysed with a proprietary Invitrogen's buffer, diluted to 0.1-0.2 mg ml −1 protein and applied to a prewashed chromatographic column. After the lysate passed through the column, the column was intensively washed to remove unbound proteins, followed by elution of an enriched fraction of phosphoproteins. The phosphoproteins were concentrated using several rounds of centrifugation with Vivaspin concentrators (Sartorius, Göttingen, Germany). Samples were stored at −80 • C in 25 mM Tris, pH 7.5, 0.25% CHAPS until further analyses.

LC/MS Analysis
Sample preparation for mass spectrometric analysis of enriched M. edulis mitochondrial phosphoproteins was performed as described in detail previously . Five biological replicates per condition including normoxic controls (i.e., isolated mitochondria prior to H/R stress) and mitochondria exposed to reoxygenation following 10 min of experimental anoxia in vitro were prepared for MS analysis. In short, protein extracts were concentrated to 15 µl, mixed with sample buffer [200 mM Tris-HCl pH 6.8, 50 mM EDTA pH 8.0, 40% glycerol, 8% SDS, 0.08% bromophenol blue, 1 mM dithiothreitol (DTT)], denatured at 90 • C for 10 min, and briefly separated in 10% SDS mini gels. After Coomassie staining, the proteincontaining gel areas (upper 10 mm of each gel band) were excised, dried in a vacuum centrifuge, reduced (20 mM DTT, 30 min, 55 • C, 800 rpm), alkylated (40 mM iodoacetamide, 30 min, room temperature in the dark), and dehydrated (100% acetonitrile HPLC grade, 2 min, room temperature). Proteins were then digested (4 µg/µl trypsin solution, 37 • C, 16 h), and peptides were purified by Stage Tipping, adjusted to 2 µg/µl and frozen at −80 • C until analysis. Purified peptides were analyzed by reversed phase liquid chromatography (LC) electrospray ionization (ESI) MS/MS. Thus, an EASY nLC 1000 (Thermo Fischer Scientific, Bremen, Germany) was coupled to a QExactive mass spectrometer (Thermo Fischer Scientific, Bremen, Germany). Peptides were loaded onto in-house selfpacked fused silica columns (100 µm × 20 cm) containing reverse-phase C18 material (ReproSil-Pur 120-AQ 3.0 µm, Dr. Maisch GmbH, Ammerbuch, Germany). Peptides were eluted using a non-linear binary gradient of 166 min from 2 to 99% solvent B [0.1% (v/v) acetic acid in acetonitrile] in solvent A [0.1% (v/v) acetic acid] at a constant flow rate of 300 nl/min. The full scan was recorded with a mass range from 300 -1650 m/z and a resolution of 70,000 at 200 m/z. The 15 most abundant ions were isolated with an isolation width of 2 Th and fragmented by higher-energy collisional dissociation (HCD) at a normalized collision energy (NCE) of 27.5. Fragment ion spectra were recorded with a resolution of 17,500 at 200 m/z. Ions with unassigned charge states as well as charge 1 and higher than 6 were excluded from fragmentation. Fragmented ions were excluded from fragmentation for 30 s. Lock mass correction was enabled.

MS Data Analysis
MS/MS spectra were extracted and analyzed using the Sorcerer Sequest software (v5.1.1, Sage-N Research) and a combined target-decoy database. As a dedicated complete M. edulis genome sequence was not yet available at the time, we constructed the database from all available related Mytilus species. To this end, all NCBI protein entries listed for the taxon "Mytilus" (77,824 proteins on 09/08/2020; Supplementary Table 1) were retrieved and clustered at 97% identity with Cd-Hit (Li and Godzik, 2006) to remove redundant proteins. A set of common laboratory contaminants was added and all sequence entries were reversed and appended to the database as decoy sequences (to allow for false discovery rate (FDR) calculation), resulting in a total of 127,890 sequences in the database. MS/MS-based peptide and protein identifications of all 10 samples were merged and evaluated with Scaffold version 4.8.4 1 , with the clustering feature enabled. Protein level FDR and peptide level FDR were set to 1%. A minimum of 2 peptides was required for protein identification. For semi-quantitative analysis and comparison of relative protein abundances between samples, normalized spectral abundance factors (NSAF) were calculated using total spectrum counts (Zybailov et al., 2006). The mass spectrometry proteomics data were deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Vizcaíno et al., 2016) with the dataset identifier PXD027955 and 10.6019/PXD027955.

Blue Native PAGE
Mitochondrial membrane proteins from isolated M. edulis mitochondria were separated by Blue Native PAGE (BN-PAGE) as previously described (Bundgaard et al., 2018(Bundgaard et al., , 2020a. Membrane proteins were extracted with 1% dodecyl maltoside (DDM) or 3 g g −1 protein digitonin on ice for 15 or 5 min, respectively. Mouse mitochondrial membrane proteins extracted with 1% DDM were used as a size marker to identify the M. edulis bands. In-gel complex I activity was assessed as a purple stain after incubation for 30 min at room temperature with 1 mg/mL nitroblue tetrazolium with 0.1 µM NADH (Greggio et al., 2017;Bundgaard et al., 2020a). The gels were scanned using a Bio-Rad ChemiDoc MP imager and analyzed using the Bio-Rad Image Lab band detection software (Bio-Rad, Hercules, CA, United States). The intensity of the supercomplex band was expressed relative to the intensity of the strong complex V (ATP synthase) band to control for loading differences.

Statistics
The significance of the effects of amino acid addition or H/R stress on the studied mitochondrial traits was tested using paired Student's t-test. The effects of H/R stress on the abundance of mitochondrial supercomplexes were assessed using unpaired t-test. Statistical analyses were conducted using GraphPad Prism v. 7.02 (GraphPad Software Inc., La Jolla, CA, United States) software. All differences were considered significant if the probability of Type II error was <0.05. Testing for statistically significant differences in protein abundance between hypoxia and control samples in the phosphoproteome NSAF data set was performed with Perseus (Tyanova and Cox, 2018) using a p-value-based Welch's t-test (p-value threshold 0.05) after imputation of missing values.

Effects of Amino Acids and Temperature on Mitochondrial Functions
Addition of 50 mM taurine or glycine had no significant effect on mitochondrial oxygen consumption or ROS production regardless of the temperature (Table 1). Therefore, amino acid addition has not been considered in the further analyses of the effects of temperature and H/R stress on the studied mitochondrial traits.
Elevated temperature (25 • C) led to an increase in mitochondrial respiration and ROS emission in the LEAK and OXPHOS states (Figure 1). Likewise, ETS activity was higher at 25 • C than at 15 • C (Figure 1). Q 10 for the mitochondrial oxygen consumption was 2.2-2.5 regardless of the mitochondrial activity state. Q 10 for ROS emission was 2.4-2.5 in the LEAK state and 1.6 in the OXPHOS state.
The rate of the ROS emission was higher in the LEAK state compared with the OXPHOS state regardless of the temperature (

Effects of H/R Stress on Mitochondria Functional Traits and Supercomplex Formation
Effects of H/R stress on mitochondrial oxygen consumption were similar at 15 and 25 • C showing suppression of LEAK respiration, elevated OXPHOS activity and increased RCR after H/R stress (Figure 2). ETS activity was elevated after H/R stress in mitochondria respiring on succinate but not in those using pyruvate or pyruvate-succinate mixture (Figure 3).
To improve the assay sensitivity and decrease the time for aerobic-anaerobic transition, all subsequent analyses (including ROS emission, supercomplex formation and phosphoproteomics) were conducted in mitochondria respiring at 25 • C where the mitochondrial activity was higher. At 25 • C, ROS emission in the LEAK state was strongly (by 56%) suppressed after the H/R exposure, whereas ROS emission in the OXPHOS state was not significantly affected by the H/R stress (Figure 4). The rate of the electron leak (measured as a ratio of H 2 O 2 production to O 2 consumption) decreased by 49 and 28% in the LEAK and OXPHOS state mitochondria, respectively, after the H/R stress ( Figure 4C). The LEAK state mitochondria converted ∼3.6% of consumed O 2 to H 2 O 2 under the control conditions, but only ∼1.9% after the H/R stress ( Figure 2C). In the OXPHOS state, the mussels' mitochondria converted ∼1.3% and <1% of consumed O 2 to H 2 O 2 under before and after the H/R stress, respectively ( Figure 4C).
Supercomplex (CI-CIII) formation was detected in mussels but it was not affected by exposure of isolated mussel mitochondria to H/R stress (Figure 5). The SC abundance relative to that of complex V was instead dependent on the use of detergent, as found previously in several species (Bundgaard et al., 2020a,b).

Mitochondrial Phosphoproteome Shifts During the H/R Stress
Proteomic analysis identified 169 mitochondrial phosphorylated proteins in the isolated mitochondria of M. edulis (Supplementary Table 2) out of which 25 (15%) significantly changed in abundance in response to the H/R stress (Welch's t-test, P < 0.05; Table 2; Supplementary Table 2). Of these, 16 phosphorylated proteins showed an increase in abundance, whereas abundance of 9 phosphorylated proteins decreased.
Among the phosphoproteins upregulated during the H/R stress, two (the NDUFAB1 subunit of the mitochondrial Complex I and A-kinase anchoring protein 1) showed an "off-on" response and were found exclusively in the samples exposed to the H/R stress but not in normoxia (Table 2). Furthermore, eight more phosphorylated proteins showed a large increase in abundance Paired t-test with Bonferroni correction showed no significant differences between the control and amino acid-treated mitochondria (P > 0.05). Respiration rates were expressed in µmol O 2 min −1 g −1 protein, and H 2 O 2 production in µmol H 2 O 2 min −1 g −1 protein.  LEAK respiration, (B) OXPHOS respiration, and (C) respiratory control ratios. The percentage change following H/R stress was calculated by standardizing the control values at the respective temperatures to 100%. Significant differences in each mitochondrial trait between normoxia (C) and reoxygenation (R) are denoted by asterisks (*P < 0.05, ***P < 0.001). n.s. indicates differences that were not statistically significant (P > 0.05).
after the H/R stress including kynurenine/α-aminoadipate aminotransferase (7.1-fold increase), propionyl-CoA carboxylase β chain (4.9-fold), cytochrome c oxidase assembly protein (3.8-fold), Golgi phosphoprotein 3 (2.9-fold), Bax inhibitor ETS, TCA cycle, amino acid metabolism, and protein processing) showed a 10-70% increase in abundance in response to H/R exposure ( Table 2). Among the phosphorylated mitochondrial proteins that were significantly less abundant after the H/R stress, two proteins (mitochondrial metalloendopeptidase OMA1 and NADH:ubiquinone oxidoreductase subunit AB1) showed an "on-off " response (i.e. were not detected in the phosphorylated fraction after the H/R stress), and two (glutaryl-CoAdehydrogenase and hydroxyacid-oxoacid transhydrogenase) showed a 2.5-fold decrease in the abundance after the H/R stress ( Table 2). Six phosphorylated proteins including the NDUFB11 subunit of the mitochondrial Complex I, F-type H + -transporting ATPase subunits γ and f, mitochondrial ribosomal protein S22, and two enzymes involved in the intermediary metabolism [ (3,5)-(2,4)-dienoyl-CoA isomerase, and 2-methoxy-6polyprenyl-1,4-benzoquinol methylase] showed a modest but statistically significant decrease (by 30-40%) in abundance after the H/R stress ( Table 2). Two more proteins showed large but statistically non-significant changes in the abundance of the phosphorylated form after H/R stress, including translocase of outer membrane 20 (TOM20) that was not found in the mitochondrial phosphoproteome after H/R stress, and reticulon-4-interacting protein that showed a >5-fold decrease in phosphorylated form abundance (Supplementary Table 2).

DISCUSSION
Mussels' mitochondria show robust coupling and 2.2-2.5-fold increase in both LEAK and OXPHOS respiration during a temperature increase from 15 to 25 • C. Elevated temperature led to higher ROS efflux in the mussels mitochondria proportional to the increase inṀO 2 so that electron leak rate (i.e., the fraction of consumed oxygen converted to ROS) did not change during warming. Interestingly, the ROS emission in the OXPHOS state was less temperature sensitive (Q 10 = 1.6) than ROS emission of the resting (LEAK state) mitochondria (Q 10 = 2.4-2.5). Such thermally robust mitochondrial function might reflect adaptations to fluctuating temperature regimes in this typical intertidal species that colonized the subtidal zone of the Baltic Sea about 10,000 years ago (Śmietanka et al., 2014).

H/R Stress Improves Mitochondrial Efficiency and Suppresses ROS Formation
Mitochondrial oxygen consumption is a driving force of the electron flux in the ETS and ATP production by OXPHOS, as well as a protective mechanism against oxygen toxicity (Sung et al., 2010). As the major oxygen consumers and ROS producers, mitochondria are key targets of H/R stress. Mitochondrial tolerance to intermittent hypoxia requires the ability to maintain high respiration flux while avoiding excessive ROS production during oxygen fluctuations . Our present study indicates that the mitochondria of the blue mussels M. edulis show features consistent with the hypoxia-tolerant phenotype. Thus, isolated mitochondria of M. edulis exposed to H/R stress in vitro showed a 5-18% decrease in proton leak rate and a 5-14% increase in OXPHOS activity resulting in higher coupling efficiency of mitochondria. Mitochondrial proton leak causes uncoupling of oxygen consumption from ATP synthesis and must be compensated by consumption of additional oxygen and substrates to maintain the activity of ETS proton pumps and prevent depolarization of the mitochondrial membrane. Therefore, lower proton leak during reoxygenation in the mussel mitochondria can contribute to improved ATP synthesis efficiency and decrease energy costs of mitochondrial maintenance. Furthermore, an increase in the OXPHOS rate after reoxygenation indicates improved ATP synthesis capacity in the mussel mitochondria. This increase appears largely driven by improved capacity of the mitochondrial ETS to utilize succinate whereas the ability to oxidize NADH-linked substrates such as pyruvate is not affected. Succinate is a common anaerobic end product in marine bivalves that accumulates in high concentrations during hypoxia (Hochachka and Mustafa, 1972;Bayne, 2017). High succinate oxidation capacity might be metabolically beneficial during post-hypoxic recovery in mussels helping to generate ATP with a highly abundant substrate (succinate) and restore normal levels of metabolic intermediates. Elevated OXPHOS and ETS capacity after H/R stress might thus assist in restoring of energy homeostasis after hypoxia and contribute to the increase in oxygen consumption during posthypoxic recovery (so called "oxygen debt") commonly found in hypoxia-tolerant invertebrates and fish (Ellington, 1983;Lewis et al., 2007;Vismann and Hagerman, 2008).
Exposure of isolated mussels' mitochondria to H/R stress strongly modulated mitochondrial ROS emission, which is FIGURE 4 | Effects of hypoxia and reoxygenation (H/R) on oxygen consumption rates (ṀO 2 ) and ROS efflux of isolated digestive gland mitochondria of M. edulis. (A) oxygen consumption rates during LEAK I+II (PLI+II) and OXPHOS respiration, (B) ROS efflux during LEAK I+II (PLI+II) and OXPHOS respiration, (C) Electron leak rate during LEAK I+II (PLI+II) and OXPHOS respiration. All activities were assessed at 25 • C. The percentage change following H/R stress was calculated by standardizing the control values at the respective temperatures to 100%. Significant differences in each mitochondrial trait between normoxia (C) and reoxygenation (R) are denoted by asterisks (**P < 0.01). considered a main driver of H/R-induced mitochondrial stress and injury (Honda et al., 2005;Cadenas, 2018). In hypoxiasensitive organisms such as terrestrial mammals, post-hypoxic reoxygenation commonly leads to a burst of ROS production causing damage to mitochondrial membrane and proteins and suppression the ETS activity and ATP production (Honda et al., 2005;Venditti et al., 2013;Cadenas, 2018). The damage to ETS further enhances electron leak and ROS production resulting in a vicious cycle of ROS-induced ROS release (Zorov et al., 2006). Unlike this detrimental pattern commonly observed in mammals, mitochondria of M. edulis suppressed electron leak in response to H/R stress. The overall rate of ROS emission and the electron leak rate were strongly (by ∼56 and ∼49%, respectively) suppressed in the LEAK state mitochondria of mussels after H/R stress. The ROS emission rate in the OXPHOS state did not significantly decrease after H/R stress, but the opposing effects of H/R stress on the oxygen consumption (i.e., a slight increase) and ROS emission (a slight decrease) combined to produce a significant (by ∼28%) suppression of the electron leak rate in the OXPHOS state mitochondria. These changes in the ROS emission and electron leak rates induced by the H/R stress can prevent oxidative stress and mitochondrial damage during reoxygenation and might reflect adaptation of the mussels to frequent fluctuation in the oxygen conditions experienced in intertidal and estuarine habitats. From the perspective of oxidative stress mitigation, suppression of the ROS emission in the LEAK state, where electron leak rates and ROS generation rates are high due to the high mitochondrial membrane potential (Miwa and Brand, 2003;Lambert and Brand, 2004), might play an especially important protective role during H/R stress in mussel mitochondria. Fold changes were calculated from NSAF values in Scaffold. INF -indicates an "on-off" response where a phosphorylated protein was exclusively found in the normoxic samples but not in those exposed to H/R stress, or vice versa. Accession numbers refer to the sequences from the NCBI reference genomes (M. coruscus or M. edulis) used for protein identification.
Earlier studies using exposures to H/R in vivo (e.g., in whole organisms) or in situ (e.g., during ischemia-reperfusion of organs, tissues, or intact cells) also show robust mitochondrial ETS and OXPHOS function in hypoxia-tolerant species contrasting the loss of OXPHOS capacity and suppressed ETS activity commonly observed in hypoxia-tolerant species (Pamenter, 2014;Sokolova, 2018;Cheng et al., 2021). Thus, OXPHOS respiration and ETS activity were greatly impaired after H/R exposure in hypoxia-sensitive organisms like rainbow trout (Oncorhynchus mykiss) (Onukwufor et al., 2014(Onukwufor et al., , 2016(Onukwufor et al., , 2017, bay scallops (Argopecten irradians) (Ivanina et al., 2016) and mammals (Paradis et al., 2016;Lesnefsky et al., 2017). A decrease in the mitochondrial OXPHOS rates was also found in isolated rainbow trout mitochondria (Onukwufor et al., 2014) and in permeabilized heart fibers of a hypoxia-sensitive shovelnose ray (Aptychotrema rostrata) (Hickey et al., 2012) after exposure to H/R stress in vitro. Mitochondrial ETS and OXPHOS rates of hypoxia-tolerant organisms such as freshwater turtles, fish, and benthic and intertidal marine bivalves (clams, oysters, mussels) remained constant or increased after H/R exposure in vivo (Hickey et al., 2012;Cook et al., 2013;Galli and Richards, 2014;Du et al., 2016;Ivanina et al., 2016;Gerber et al., 2019;Napolitano et al., 2019;Ouillon et al., 2021) but not after anoxic exposure in freshwater turtles, where mitochondrial respiration was suppressed (Galli et al., 2013;Bundgaard et al., 2018Bundgaard et al., , 2020a. Similarly, an increase in OXPHOS activity was found in isolated mitochondria of hypoxia-tolerant zebrafish (Danio rerio) after H/R exposure in vitro (Napolitano et al., 2019) but a decrease in OXPHOS was found after in vivo prolonged hypoxia of zebrafish in permeabilized skeletal muscle fibers (Cadiz et al., 2019). Taken together, our present findings and earlier published research indicate that, although different responses to hypoxia may take place in some vertebrate species (Pamenter, 2014;Pamenter et al., 2016Pamenter et al., , 2018Bundgaard et al., 2020b), in bivalves the mitochondrial phenotype associated with hypoxia tolerance involves the ability to maintain or enhance OXPHOS capacity, mitigate ROS emission and prevent collapse of the ETS during reoxygenation, and that these mitochondrial traits are partially independent of cellular or organismal regulatory mechanisms.

Mitochondrial Proteome Rearrangements During H/R Stress: Dynamics of Supercomplexes
Supercomplexes are dynamic structures consisting of different combinations of complexes involved in the ETS system that can undergo rearrangement in response to varying physiological or environmental conditions including hypoxia (Ramírez-Aguilar et al., 2011;Timón-Gómez and Barrientos, 2020). In hypoxiaexposed mammalian mitochondria, CI was reversibly deactivated (Galkin et al., 2009) likely due to dissociation of the large supercomplexes containing CI, CIII, and CIV, although the capacity of the individual respiratory complexes was intact (Rosca et al., 2008). In plants, SCs were observed to dissociate and reassemble during H/R stress as a result of cellular acidification (Ramírez-Aguilar et al., 2011). However, our present study in mussels as well as recent research in hypoxia-tolerant freshwater turtles (Bundgaard et al., 2018) shows that the dynamic rearrangement of SCs is not involved in response to hypoxia in these species. In M. edulis, H/R stress did not influence the SCs' abundance (Figure 5), and stability of CICIII 2 SC during H/R stress was reported in freshwater turtles (Bundgaard et al., 2018). The functional significance of SC rearrangement remains debatable (Signes and Fernandez-Vizarra, 2018;Timón-Gómez and Barrientos, 2020), and the role of SC stability in maintaining robust ETS activity or in limiting ROS during oxygen fluctuations in hypoxia-tolerant organisms requires further investigation. It is worth noting that SC of ectotherms appear generally more stable than those of mammals, regardless of the species-specific hypoxia tolerance (Bundgaard et al., 2020a,b). Thus, it is possible that the stability of SCs in the mitochondria of mussels observed in our present study is a shared ectotherm trait rather than adaptation to hypoxic stress.

Phosphoproteome Shifts
Reversible phosphorylation is emerging as a pivotal posttranslational modification (PTM) mechanism involved in the regulation of mitochondrial functions (Pagliarini and Dixon, 2006;Helling et al., 2012). The steadily increasing numbers of reported mitochondrial kinases, phosphatases, and phosphoproteins imply an important role of protein phosphorylation in different mitochondrial processes, including responses to H/R stress (Goldenthal and Marín-García, 2004;Horbinski and Chu, 2005). Mitochondrial proteome phosphorylation patterns and their role in the responses to H/R stress are less well studied in hypoxia-tolerant organisms including mollusks. Differential phosphorylation of proteins associated with OXPHOS, Krebs cycle, and intermediary metabolism has been reported during H/R exposure in vivo in the hypoxia-tolerant Pacific oyster, C. gigas . Furthermore, H/R exposure in vivo modulated activities of protein kinases A and C in the mitochondria of M. edulis, C. gigas, and Arctica islandica . Our present study demonstrates that reversible protein phosphorylation is also involved in the intrinsic mitochondrial responses to H/R stress that act independently of the cellular context during in vitro exposures to fluctuating oxygen conditions. The pathways regulated by reversible phosphorylation during H/R stress in isolated mussel mitochondria include ETS and OXPHOS (10 out of 25 differentially phosphorylated proteins), amino acid and fatty acid metabolism (6 proteins), mitochondrial protein and mtDNA homeostasis (5 proteins), and substratelevel phosphorylation (2 proteins). It is worth noting that in our present study, the change in the protein abundance in the phosphoprotein fraction after H/R stress was interpreted as a change in the phosphorylation state of the respective protein due to PTM. In principle, changes in the abundance of a certain phosphoprotein might also reflect changes in the total protein abundance in the mitochondria due to de novo synthesis and/or protein degradation. However, these alternative explanations appear less likely given the cell-free system and the short exposure times (10 min in anoxia followed by immediate sample collection after reoxygenation) used in the present study.
Of note, H/R exposure of isolated mussels' mitochondria led to a major increase in abundance of phosphorylated A-kinase anchoring protein 1 (AKAP1). AKAP1 is a scaffold protein that recruits protein kinase A (PKA) and other signaling proteins, as well as RNA, to the outer mitochondrial membrane (Liu et al., 2020). AKAP1 acts as a signaling hub regulating metabolic homeostasis and mitochondrial quality control and is involved in response to H/R stress in mitochondria (Merrill and Strack, 2014;Liu et al., 2020). In mammalian models, AKAP1 becomes phosphorylated in response to energy stress via AMPactivated protein kinase (AMPK)-dependent mechanism, and this phosphorylation is essential for activation of mitochondrial respiration (Liu et al., 2020). Although the mechanisms of AKAP1-dependent mitochondrial regulation are not yet fully unraveled, multiple studies indicate that functional AKAP1 is essential for cell survival during H/R stress and is involved in regulation of PKA activity in mitochondria [review in Merrill and Strack (2014); Liu et al. (2020)]. If similar mechanisms exist in mussel mitochondria, phosphorylation of AKAP1 might indicate an adaptive response to a decrease in ATP levels during H/R stress and positively modulate mitochondrial respiration and PKA recruitment. The latter hypothesis is also consistent with an important role of PKA in regulating the mitochondrial ETS activity in mussels .
Mitochondrial Complex I is an important target for PTMdependent regulation during H/R stress in animals including mammals (Gowthami et al., 2019) and mollusks . Complex I activity is often down-regulated in response to H/R stress as a mechanism to suppress mitochondrial ROS generation by reversible electron transport (RET) (Pell et al., 2016); however, the role of reversible protein phosphorylation in this regulation is not known. In marine bivalves including mussels, Complex I activity is responsive to experimental manipulation of the phosphorylation status . In M. edulis, PKA-dependent phosphorylation led to an increase in Complex I activity whereas PKC-dependent phosphorylation had no effect . Unspecific dephosphorylation (by cerium oxide) suppressed Complex I activity . Our present study strongly indicates that H/R stress affects phosphorylation status of several Complex I subunits in M. edulis. Thus, the subunit NDUFAB1 of Complex I was exclusively phosphorylated in the mitochondrial samples exposed to H/R, and the abundance of the phosphorylated forms of NDUFS2 and NDUFB6 increased by 30-60% after the H/R stress. In contrast, the abundance of the phosphorylated NDUFB11 decreased by 40% after the H/R stress. Functional implications of PTM of these subunits are not known in M. edulis. Interestingly, the subunit NDUFAB1 that is exclusively phosphorylated in the mussels after H/R stress is known to play a protective role during H/R stress in mammalian mitochondria supporting high respiratory flux, stabilizing supercomplexes and mitigating ROS generation (Hou et al., 2019). It is tempting to speculate that a similar protective mechanism mediated by the PTM of NDUFAB1 subunit might be in play in the mussels' mitochondria during H/R stress, but this hypothesis requires further investigation.
It is worth noting that mitochondrial Complex IV was not affected by PTM during H/R stress in the mussels' mitochondria. In mammals, Complex IV is considered a key regulatory target during hypoxia whose activity and oxygen sensitivity is modulated by PTM (including reversible phosphorylation) (Prabu et al., 2006;Kadenbach, 2021) and differential expression of hypoxia-specific subunits (Kocha et al., 2015). In M. edulis, Complex IV can also be modulated by experimental phosphorylation/dephosphorylation, albeit to a lesser extent than Complex I . In marine mollusks including M. edulis, cytochrome c oxidase is not considered rate-limiting due to the high apparent excess capacity of this enzyme relative to the maximum ETS activity Ouillon et al., 2021). Furthermore, oxygen affinity of cytochrome c oxidase does not appear modulated by H/R stress in marine bivalves (Ouillon et al., 2021). This indicates that cytochrome c oxidase might be less important as a regulatory target of ETS during H/R stress in mollusks. Consistent with this notion, none of the Complex IV subunits were identified as differentially phosphorylated during H/R stress in our present analysis. Interestingly, H/R stress led to a major increase (by ∼3.8-fold) of the phosphorylated form of cytochrome c oxidase assembly protein COX11. This protein acts as copper chaperone delivering Cu 2+ to the active center of cytochrome c oxidase and might be involved in mitigation of ROS generation in mitochondria (Radin et al., 2015). It is possible that elevated phosphorylation levels of this protein might affect redox properties or Cu 2+ binding affinity of this protein; however, at present the functional implications of COX11 PTM remain speculative.
H/R stress resulted in a modest (by ∼30%) but statistically significant decrease in the abundance of phosphorylated subunits γ and f of the mitochondrial F o , F 1 -ATPase in the mussels and a ∼70% increase in abundance of phosphorylated mitochondrial protease ATP23 responsible for F 0 , F 1 -ATPase assembly. These findings indicate that F o , F 1 -ATPase is another target of regulatory PTM during H/R stress in bivalve mitochondria. Phosphorylation of ATPase subunits (particularly subunit γ) plays an important role in F o , F 1 -ATPase assembly (Reinders et al., 2007). In yeast, dephosphorylation of subunit γ facilitates oligomerization of F o , F 1 -ATPase enhancing its activity (Reinders et al., 2007). If a similar mechanism exists in mussels, dephosphorylation of subunit γ of F o , F 1 -ATPase would activate F o , F 1 -ATPase through enhanced oligomerization in the inner mitochondrial membrane and might contribute to the observed increase in the OXPHOS activity after the H/R stress in mussels.
Amino acid and fatty acid catabolism in mitochondria represents another key target for regulation through PTM in bivalve mitochondria as shown by differential phosphorylation levels induced by the H/R stress in multiple enzymes involved in the respective pathways. Overall, six enzymes involved in amino acid and fatty acid catabolism appeared differentially phosphorylated in the control and H/R exposed mitochondria of mussels ( Table 2). Most of these phosphorylated enzymes showed a modest change in abundance (typically, ∼20-40% change compared with the controls), except glutaryl-CoA dehydrogenase (showed a 2.5 decrease in the phosphorylation level) and phosphorylated kynurenine/alpha-aminoadipate aminotransferase and propionyl-CoA carboxylase beta chain that showed a ∼−5-7-fold higher abundance after H/R stress. In oysters (C. gigas), phosphorylated forms of enoyl-CoA isomerase and acyl-CoA dehydrogenase family enzymes (involved in amino acid and fatty acid catabolism) decreased in abundance after H/R stress in vivo . Presently, the functional consequences of reversible phosphorylation of these enzymes are not known but their high representation in differentially expressed mitochondrial phosphoproteome of bivalves after H/R stress indicates the importance of modulation of amino acid and fatty acid metabolism, possibly to facilitate the use of diverse mitochondrial substrates for ATP production during post-hypoxic recovery.
Mitochondrial protein quality control mechanisms play a key role in maintaining the integrity of the mitochondrial proteome and are commonly regulated during H/R stress at the transcriptional, translational and post-translational levels Falfushynska et al., 2020;Steffen et al., 2020). Among these proteins, mitochondrial metalloendopeptidase OMA1 appears an important target for PTM. In M. edulis, the abundance of phosphorylated OMA1 form was strongly down-regulated in H/R stress, so that no phosphorylated form of OMA1 was found in the mitochondrial proteome after H/R stress. Similarly, in C. gigas exposed to hypoxia in vivo the abundance of phosphorylated OMA1 decreased ∼6-fold in hypoxia, and no phosphorylated OMA1 was found after 1 h of reoxygenation . These findings indicate that dephosphorylation of OMA1 plays an important role in response to H/R stress, particularly the post-hypoxic recovery, in facultative anaerobes such as marine bivalves, and can occur via intrinsic mitochondrial mechanisms. Dephosphorylation of OMA1 in response to low mitochondrial ATP levels might activate this enzyme, thereby fostering mitochondrial fragmentation and suppressing mitochondrial fusion in stress (Baker et al., 2014;Bohovych et al., 2015;. Other important proteins involved in the maintenance of mitochondrial proteome showed an apparent change in phosphorylation status after the H/R stress (including the mitochondrial ribosomal protein S22, TOM20 and mitochondrial processing peptidase subunit beta), consistent with the key role of protein homeostasis in mitochondrial stress resistance (Sokolova, 2018). However, the functional implications of PTM of these proteins are not yet known even in model organisms and therefore, no conclusions can presently be drawn with regard to their possible effects on mussel mitochondria.
In conclusion, our study shows that isolated digestive gland mitochondria of the blue mussel M. edulis possess intrinsic regulatory mechanisms that improve mitochondrial efficiency and limit ROS generation during hypoxia-reoxygenation transition. These adaptive functional changes are triggered by oxygen fluctuations in isolated mitochondria within minutes and thus reflect autonomous mitochondrial mechanisms that do not require an intact cellular environment or retrograde nuclear signaling. The intrinsic mitochondrial responses to H/R stress are at least in part dependent on the shifts in the mitochondria phosphoproteome including such key pathways as OXPHOS, amino acid and fatty acid catabolism, and protein quality control, but do not appear to involve a change in SC formation. The findings of our present study complement earlier research that showed an important role of PTM in regulating glycolytic metabolic fluxes and substrate-level phosphorylation Storey, 1995, 1997;Russell and Storey, 1995;Dawson et al., 2013;Lama et al., 2013;Storey, 2015) and highlight the importance of reversible protein phosphorylation as a global metabolic regulatory mechanism during oxygen fluctuations in facultative anaerobes.

DATA AVAILABILITY STATEMENT
The mass spectrometry proteomics data were deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Vizcaíno et al., 2016) with the dataset identifier PXD027955 and doi: 10.6019/PXD027955.

AUTHOR CONTRIBUTIONS
ES and AB: methodology, investigation, data curation, validation, formal analysis, and writing -review and editing. LA: data curation, visualization, writing -original draft, and writing -review and editing. SM: data curation, validation, formal analysis, resources, and writing -review and editing. AF: methodology, resources, supervision, and writingreview and editing. DB: methodology, resources, supervision, validation, and writing -review and editing. CH: methodology, investigation, data curation, validation, and writing -review and editing. IS: conceptualization, methodology, validation, formal analysis, resources, supervision, project administration, funding acquisition, writing -original draft, and writing -review and editing. All authors contributed to the article and approved the submitted version.