Exploring the mechanisms behind swimming performance limits to ocean warming and acidi ﬁ cation in the Atlantic king scallop, Pecten maximus

Recently, we could show that scallops show limitations of muscular performance like a reduced force under ocean warming and acidi ﬁ cation. However, the underlying mechanisms at the cellular level are not completely understood. Metabolomics has become a valuable tool to evaluate the responses of marine organisms to various stressors. In the present study we therefore used a semi-targeted, multi tissue NMR based metabolomic approach to analyze metabolite patterns in the Atlantic king scallop, Pecten maximus , that were long-term acclimated to different end of century conditions of ocean warming (OW), ocean acidi ﬁ cation (OA) and their combination (OWA). We investigated tissue speci ﬁ c metabolic pro ﬁ les and metabolite concentrations in frozen tissues from gills, mantle and phasic and tonic adductor muscle of P. maximus under present conditions using 1 H-HR-MAS NMR spectroscopy. A set of 33 metabolites revealed a clear tissue-speci ﬁ c pattern which can be attributed to the individual functions of the respective tissue type. We then evaluated the impact of OW, OA and OWA on the metabolic pro ﬁ les of the different tissues. OW was the main driver of the changes in metabolites. In particular, energy-related metabolites seem to play an important role in the physiological response of scallops to OW and OWA. In combination with pathway analysis and network exploration we propose a possible correlation between metabolic changes in the adductor muscle and limited swimming performance of P. maximus under future climate. While the metabolic response of the phasic muscle seems to mainly depend on net consumption of energy related metabolites such as ATP and phospho-L-arginine, the tonic muscle seems to rely on metabolizing speci ﬁ c amino acids and beta-oxidation to account for the elevated energetic requirements under ocean warming and acidi ﬁ cation.


Introduction
Increasing greenhouse gas emissions are not only responsible for the ongoing global warming trend, but as anthropogenic CO 2 has been taken up by the oceans (Bindoff et al., 2019;Licker et al., 2019), they also lead to a decrease in ocean pH through the marine carbonate system (Zeebe, 2012).This process known as "ocean acidification" (OA) is a global phenomenon and called the evil twin of warming, since both ocean warming and acidification are occurring in parallel.It is therefore important to investigate the influence of both factors on the physiology of marine organisms in order to be able to make realistic statements about climate-induced changes in marine ecosystems (Pörtner et al., 2014;Gattuso et al., 2018;Landrigan et al., 2020).Attempts to gain an understanding of the combined effects of elevated temperatures and PCO 2 levels on marine species in their ecosystem have been underway for some time (Doney et al., 2012;Byrne and Przeslawski, 2013), and the effects of these climate variables turned out to be context and species-dependent (Harvey et al., 2013).The combination of both factors has recently been reviewed (Baag and Mandal, 2022) and literature shows that ocean warming is the main driver affecting an organism's metabolism while the effects of acidification play more of an additive role (e.g.Pinsky et al., 2019;Matoo et al., 2021).
For water-breathing poikilotherms whose body temperature varies with environmental temperature, the concept of oxygen and capacity limited thermal tolerance highlights that, with increasing temperature, marine ectotherms experience progressive hypoxemia and finally, when reaching the critical temperature, rely on anaerobic ATP production to fuel the rising gap in energy demand, which is no longer supported by solely aerobic mitochondrial ATP production (Pörtner, 2002;Pörtner et al., 2017).Accordingly, since the aerobic power budget (excess of aerobic energy after maintenance costs are covered; Guderley and Pörtner, 2010;Pörtner et al., 2017) is limited and an organism cannot rely long-term on anaerobiosis, there are trade-offs between an organism's performance parameters (growth, exercise, reproduction) and there are bottlenecks in an organism's performance capacity under stressful (energy-demanding) environmental conditions such as ocean warming and acidification (OWA).
The impact of ocean acidification on marine life is manyfold (see reviews by e.g.Melzner et al., 2020;Doney et al., 2020;Findlay et al., 2022) and various calcifiers suffer from reduced growth and calcification rates even though the previously high level of concern about the fate of calcifiers in a future ocean no longer holds (see reviews by Zhao et al., 2020;Leung et al., 2022).Bivalves have a limited capacity for acid base regulation and are experiencing a drop in extracellular pH when exposed to OA (Lannig et al., 2010;Heinemann et al., 2012;Gazeau et al., 2014;Zhao et al., 2017) including the scallop Pecten maximus (Schalkhausser et al., 2013;Schalkhausser et al., 2014).Such studies indicate that OA impact and sensitivity differ between wild and aquaculture animals (Richards et al., 2015;Stapp et al., 2018).Among bivalves, scallops are unique due to their swimming (escape) behavior via jet propulsion through adductor muscle contractions (Guderley and Tremblay, 2016).The adductor muscle consists of two types that interact during swimming exercise (Peŕez et al., 2008;Tremblay et al., 2012;Guderley and Tremblay, 2016).The striated phasic muscle, which is essentially responsible for swimming through claps, and the smaller smooth tonic muscle, which is mainly responsible for closing the shells and holding them closed (Chantler, 2006;Sun et al., 2018) During swimming, muscle contraction cannot solely be fueled by mitochondrial oxidative phosphorylation and is highly dependent on ATP generation from phosphate reserves like phospho-L-arginine (PLA) and anaerobic glycolysis (Guderley and Tremblay, 2016).In comparison to other bivalves like oysters or mussels that are able to colonize the intertidal zone, their active lifestyle gives scallops less flexibility to cope with changing environmental factors (Guderley and Pörtner, 2010;Ivanina and Sokolova, 2016;Götze et al., 2020).

Previous observations on scallops swimming performance
In previous studies, we investigated the influence of ocean warming and acidification on swimming performance of the king scallop, Pecten maximus (Schalkhausser et al., 2013(Schalkhausser et al., , 2014)).According to Schalkhausser et al. (2014) P. maximus' swimming ability was affected by warming more than by acidification (PCO 2 0.112 kPa).While the number of claps to exhaustion was similar for P. maximus between 10°C and 20°C, the force of the phasic adductor muscle was significantly reduced under long-term acclimation to 20°C.The time to exhaustion was reduced and the recovery period prolonged in warm-exposed P. maximus.Warming resulted in a significant decrease in hemolymph oxygen levels (extracellular partial pressure of O 2 : P e O 2 ).In combination with a significant increase in respiration rate at 20°C, the nearly 50% reduction in P e O 2 indicated that P. maximus had surpassed its optimal temperature range and was exhibiting a progressive warming-induced mismatch between aerobic energy supply and demand (Pörtner et al., 2017).In a subsequent study on P. maximus exposed long-term to OW (20°C, 0.04 kPa PCO 2 ) and OWA (20°C, 0.112 kPa PCO 2 ; Bock et al., 2019), we investigated the availability of energy-rich phosphates, in particular PLA and inorganic phosphate (Pi) using in vivo 31 P NMR spectroscopy and calculated the mitochondrial maximum surplus oxidative flux of the phasic muscle.In both OW-and OWA-exposed scallops, PLA and maximum surplus oxidative flux of phasic muscle were reduced, which was attributed to the lower P e O 2 in the haemolymph found in OW-and OWA-exposed scallops described in Schalkhausser et al. (2014).
The aim of the present study was to use state-of-the-art metabolomics and thereby elucidate the mechanisms at the cellular level which underpin the performance limitations on the organism-level in P. maximus.Metabolomics enables the holistic tracking of changing metabolic processes in response to an external stimulus such as environmental changes (Viant, 2008, for a recent review see Wishart, 2019).Metabolomics is believed to unravel the phenotype of an organism from the response to an external stimulus (Fiehn, 2002;Wishart et al., 2022).Using tissue samples from our previous experiments (see Schalkhausser et al., 2014), we applied a semi-untargeted metabolomics approach based on NMR spectroscopy to a multi-tissue dataset for obtaining tissue specific insights into the metabolic response of Pecten maximus.In addition, metabolic responses of the adductor muscle were then related to whole-animal performance and capacity parameters (clapping response and oxygen consumption of P. maximus as described in Schalkhausser et al., 2014) to gain a more integrative understanding.To further support our interpretation of the data, a pathway analysis was complemented by a network exploration.This combination enables a better interpretation of the metabolite changes and the generation of testable experimental hypotheses.
In particular, we focused on potential pathway modifications of the intermediary metabolism of the two adductor muscle types (phasic and tonic muscle).We propose a cellular mechanism that might explain the observed limitations in whole-animal performance of P. maximus under future climate change (Schalkhausser et al., 2014;Bock et al., 2019).

Experimental design and tissue collection
All investigations were performed on tissue samples that were taken from our experiment described in Schalkhausser et al. (2014).Briefly, wild living Atlantic king scallop Pecten maximus with similar size were collected by SCUBA divers at Morlaix Bay (Baie de Morlaix, Les Grandes Fourches, France; 48°42′33.6″N,3°55′59.30″W)and transported to the Alfred Wegener Institute (Bremerhaven, Germany).Animals were considered as typical from the sampling spot.After recovery from transportation for 2 weeks, scallops were cleaned from epibionts and randomly placed in recirculation aquarium systems and exposed to two different temperatures (10°C; C vs. 20°C; OW) and two different PCO 2 values (~0.04 kPa (400 matm); OA vs. ~0.112kPa (1,120 matm); OWA) for at least 50 days.Three times per week scallops were fed live phytoplankton (DT´s Premium Reef Blend).Mortality was 9% for the C group, 0% for OA, 20% for OW and 9% in OWA group at the end of the experiment (Nov.2011), but not statistically different between groups (see Schalkhausser et al., 2014).Scallops were dissected on ice, tissues (gills, mantle, phasic and tonic muscle of the adductor muscle) were freeze-clamped, shock-frozen in liquid nitrogen and stored at −80°C until present metabolite analyses were performed.

Metabolite analysis via NMR spectroscopy
Metabolite profiling of the different tissues was conducted using high resolution magic angle spinning nuclear magnetic resonance spectroscopy ( 1 H-HR-MAS NMR spectroscopy) using a wide-bore 400 MHz NMR spectrometer (9.4 T WB with Avance III HD electronics, Bruker Biospin, Germany) equipped with a triple tuned 1 H-13 C-31 P-HR-MAS NMR probe.All measurements were performed on frozen tissue samples similar to Podbielskie et al. (2016).Briefly, a sample was taken of the frozen organ of interest using a biopsy punch needle (3 mm diameter) and immediately placed in a standard zirconium HR-MAS rotor for untargeted metabolite profiling based on 1 H-NMR spectroscopy.A drop of D 2 O (containing 0.05% TSP as standard) was added to the sample for lock and calibration purposes.All samples were measured at 4°C at a spinning rate of 3000 Hz using the Bruker acquisition software TopSpin 3.5pl.NMR parameters were as follows (see also Schmidt et al., 2017;Rebelein et al., 2018): pulse program: 1D-Carr-Purcell-Meiboom-Gill (CPMG) pulse train including f1 presaturation (Bruker protocol cpmgpr1d), 90°bp pulse length 8.4 ms, time domain 70,656, sweep width of 8802 Hz (22 ppm), acquisition time 4.01 s, relaxation delay 4 s, four dummy scans and 64-256 number of scans depending on the signal-to-noise ratio.

NMR processing
All NMR data were processed and analyzed using Chenomx NMR suite 8.4 professional suite (Chenomx Inc., Canada).Data were automatically zero-filled, processed with an exponential multiplication of 0.3 Hz, phase and base-line corrected using the automatic functions within Chenomx.The spectra were automatically shim-corrected and calibrated to the TSP standard.After completion of the processing procedure, all spectra from each individual tissue of the single experimental groups were overlaid within the Chenomx profiler.The superposition of the individual spectra allows a quick overview of the quality of the individual NMR spectra and the assessment of comparability.In this way, differences in the line widths and shifts of individual NMR signals can be specifically identified and different baselines between the spectra can be recognized.These factors are essential for an automatic evaluation of metabolite profiles.Shifts of NMR signals, just like different baselines between individual spectra, would lead to an erroneous determination of metabolite concentrations.Interestingly, there was no need for alignment shifts of NMR signals in individual spectra.Spectra that differed in baseline in comparison to the other NMR spectra were marked and later identified as potential outliers within MetaboAnalyst (MetaboAnalyst 5.0., see below).A standard metabolite profile from scallop gill tissue from our previous study (Götze et al., 2020) was overlaid to an example spectrum of each tissue and checked for consistency of the metabolite specific NMR peaks.A set of 33 metabolites was identified and was used for a targeted binning procedure on batches of spectra from each tissue within Chenomix.The water region (4.7-5.25 ppm) was excluded and removed from the spectra.Integrals of the metabolite specific bins were then automatically calculated in arbitrary units, normalized to integers and exported into a csv file.For further comparative statistical analysis, the csv files of the individual spectra were merged.For this, only those bins from the individual metabolites were used for the additional analysis that were distinct for the identification and quantification of the metabolite profiles, or showed a significant content, to get an unambiguous assignment of the metabolites.This included both the integral and the position of the bins in the spectrum that could be uniquely identified.

Statistics
All spectral data were analyzed for differences in spectral patterns and individual metabolites using the web-based analysis tool MetaboAnalyst (MetaboAnalyst 5.0).The entire analytic procedure followed the recommendation for metabolomic analysis as described in Chong et al. (2019) and corresponding YouTube tutorials of the same group.In a first step, the csv files of the metabolite specific bins were imported into MetaboAnalyst and normalized using pareto scaling as recommended.The data were checked for outliers using an operator controlled visual inspection of the particular heatmaps in combination with an unsupervised principal component analysis (PCA).Identified outliers were removed from the list before further analysis (4 samples for gill and 1 sample for tonic muscle).The identified outliers matched the spectra that had already been marked, as they had different baselines compared to the other NMR spectra (see above).Dendrograms of the data excluding the outliers were used to observe clusters of specific sample spectra between groups.Differences in patterns and classification discrimination between groups were determined by Partial Least Square-discriminant analysis (PLS-DA).Metabolites of interest were classified using the VIP (variance of importance) score of the PLS-DA.Classification and cross validation by permutation tests based on separation distance were performed in view of possible overfitting and meaning fullness of the PLS-DA.A significance analysis of microarray/metabolites (SAM) analysis was performed to identify significant different metabolites between groups.The delta value to control the false discovery rate (FDR) was set conservatively to the highest possible level (leading to a maximum FDR of 0.1).For the discussion we only considered those metabolites that showed definite significant changes in SAM.In addition, ANOVA combined with a Post-Hoc test (Fisher's) was used to underpin the results from the SAM analysis where possible.To better understand the interdependencies of the altered metabolites and a possible link to the observed limitations at the whole-animal level, a pathway analysis (PW analysis) was performed in combination with a network (NW) exploration in the muscle samples of control and warming exposed scallops.While the PW analysis establishes a connection to metabolic pathways on the basis of the altered metabolites, the NW exploration looks at the interaction of metabolites and their exchange and thus points to novel interactions.However, the two methods only allow two groups to be compared with each other.Due to its specific role in swimming performance, these analyses were therefore performed on the adductor muscle between the control group and warming.For the network exploration, the metabolite interactions between the control group and the respective experimental groups were also included for comparison.The pathway analysis was performed within MetaboAnalyst, similar to (Götze et al., 2020), to get an idea of the long-term induced changes in metabolic pathways limiting the swimming performance of the adductor muscle.Only those pathways with a significant impact >0.1 were taken into account.For network exploration a metabolite-metabolite interaction network based on STITCH ("search tool for interactions of chemicals") was used to elucidate potential linkages.STITCH is a searchable database that includes interactions for over 300,000 molecules and 1.6 million proteins.The similarity of chemical structures in combination with text mining are used to make predictions about relationships between molecules.For the network exploration the filter for nodes based on degree and a filter for nodes based on betweenness were adjusted for each exploration individually to reduce the complexity of the network and to identify those metabolites that can act as specific bottlenecks in the metabolic network (see figure legends for exact values).Both analyses were only performed on phasic and tonic muscle tissues to underline the particular interaction of the two tissues in relation to the observed changes in swimming performance of scallops under warming conditions.

Results
Figure 1 shows a comparative overview of 1 H-HR-MAS NMR spectra of the four tissues examined from Pecten maximus under control incubations.Although signal intensity varied between investigated tissue types, tissue-specific spectral patterns were quite similar.All large and medium sized signals are seen in each spectrum and there were no obvious additional signals that were observed in one or some tissues only.A set of 33 metabolites was automatically aligned to all spectra of tissue type.A closer look at the individual spectra, however, revealed tissue-specific differences between metabolite profiles with respect to signal intensity, and thus metabolite concentration.The spectrum from the mantle tissue biopsy showed the broadest line widths in relation to the standard (TSP).However, the quality was sufficient for a comparative approach between tissues (see Supplementary for example spectra of the different tissues).Figure 1 displays a multivariate comparison between the tissues by an unsupervised principal component analysis (PCA) and a supervised partial least square discriminant analysis (PLS-DA).Both analyses highlight significant differences in the metabolite profiles between the individual tissues, which were also evident in the visual comparison of the individual NMR spectra (see Supplementary), including differences between the two muscle types.The three principal components of the PCA explained 94.3% of the total variance (Figure 1A) and the three components of the PLS-DA produced a clear separation between the four tissue types.Interestingly, more than 70% were explained by the first component (Figure 1B).The clear separation is confirmed by the corresponding scree plot for the PCA.The cross validation as well as the permutation tests for the PLS-DA indicated a good prediction capability and meaningfulness (see Supplementary).The significantly different metabolite concentrations and the corresponding correlations between tissues after Anova analysis and post hoc test are shown in Table 1.Gills and mantle showed the highest concentrations of osmotically relevant metabolites such as imidazole, homarine and glycine, as well as of the two sugars (glucose and UDP glucose).In contrast, the two muscle tissues showed increased concentrations of the membrane-relevant metabolites choline and O-phosphocholine.L-arginine and carnitine were also significantly more concentrated in the two muscle types compared to gill and mantle.Interestingly, the tonic muscle showed metabolite concentrations that were not observed as prominently in the other tissues.These were high concentrations of hypotaurine, tyrosine, methionine, succinate and lactate, whereas  ATP concentrations were lowest.The phasic muscle contained the highest concentrations of the amino acids alanine, glutamate and glutamine and of acetate, whereas the gill tissue exhibited the highest trigonelline levels.The mantle showed high concentrations of threonine and valine, but very low betaine concentrations compared to the other tissues.
The metabolic response to ocean acidification (OA), warming (OW) and the combination of both (OWA) did partially differ between tissue types.Figure 2 shows a graphical representation of the significant treatment-dependent metabolite changes of gill and phasic muscle tissues.According to the SAM technique, tonic and mantle tissue did not show any significant changes.In gill tissue only the energy related metabolite ATP changed in comparison to control conditions, with ATP levels showing a decrease under OW and OWA (Figure 2A).Most of the changes in concentration were observed in the phasic muscle.Here, eleven metabolites showed concentration changes.These were the essential amino acids leucine and threonine, but also amino acids involved in energy metabolism, such as L-arginine, alanine and glutamine and glutamate, which were lower in the warm groups compared to the control.Acetate and O-phosphocholine, which are utilized in membrane lipid metabolism, were also lower in the warming groups.Ophosphocholine was lowest in the OA group.A similar pattern was observed for lactate.In addition, succinylacetone showed decreased levels in the warming groups.However, glycolate showed increased concentrations under OW and OWA (Figure 2B).

Pathway and network analysis
A combination of pathway analysis and network correlation were performed on the phasic and tonic tissues from the adductor muscle in particular to get more insights of the metabolic mechanisms behind the swimming performance limits observed under ocean warming in comparison to control conditions.The metabolic pathways with the greatest impact identified in the PW analysis are listed in Table 2.The comparison between the metabolic profiles of the adductor muscle of control and OW- Changes in gill tissue (A) and phasic muscle (B) metabolite concentrations in P. maximus exposed to different treatments, controls (C: 10°C, 0.04 kPa PCO 2 ), ocean acidification (OA: 10°C, 0.112 kPa PCO 2 ), ocean warming (OW: 20°C, 0.04 kPa PCO 2 ) and combination of both (OWA: 20°C, 0.112 kPa PCO 2 ).Significant differences to control are marked with "a".3).A reduction of the network to the most important metabolites (nodes) and their connections (betweenness), highlights the central role of L-arginine, in addition to ATP, in the network.L-arginine is directly linked to oxygen, ATP and their degradation products, indicating a central role of L-arginine in energy supply during long-term acclimatization to warming in the adductor muscle.In addition, it is noteworthy that the amino acid L-methionine is in close proximity to L-arginine and ATP and that L-arginine is directly connected to carbon dioxide and sodium.In addition, network explorations from the phasic muscle under OA, OW and OWA were performed for comparison.The reduced networks confirmed the central role of L-arginine under OW and OWA, whereas under OA the metabolite phosphorylcholine was most connected (see Supplementary).

Discussion
In the present study, a retrospective metabolomic analysis using 1 H-HRMAS NMR spectroscopy was conducted using tissue samples of the king scallop, Pecten maximus, from long-term incubation experiments.Our aim was to gain a better understanding of the cellular mechanisms that respond to ocean warming (OW) and acidification (OWA) in scallops, and to relate the observed changes in metabolic pathways to the reduction in swimming performance (Schalkhausser et al., 2013, Schalkhausser et al., 2014;Bock et al., 2019).

Tissue-specific metabolite concentrations
Following our previously published targeted metabolomics approach we could identify a set of 33 metabolites for all tissues, similar to that defined in gill tissue of P. maximus (Götze et al., 2020).Similarly, Cappello et al. (2018) found only small differences between the metabolic profiles of mussel Mytilus galloprovincialis tissues (gills, phasic muscle and digestive gland).While a main set of metabolites is conserved across all tissue types, metabolite concentrations varied strongly between tissues indicating a differential prioritization of metabolic pathways in each organ.
In gills we found high concentrations of methionine and glycolate.Both metabolites have in common that they have antioxidant properties (methionine: Martıńez et al., 2017; glycolate for mitochondria: Diez et al., 2021).The main function of gills is gas exchange making their large surface areas a target for (external and internal) fluctuations of oxygen radicals.In coastal environments, in particular in lagoons or tidal pools, oxygen levels can vary from hyperoxic during day time to hypoxic conditions during night time.Thus, it is not surprising that gills possess a higher concentration of metabolites associated with oxidative stress protection than muscle.Furthermore, gill tissue, which has an important role in osmo-and ion-regulation relies on organic osmolytes, such as compatible solutes like homarine, the sugars glucose and UDP-glucose, but also the amino acid glycine.The high imidazole concentration in gills can be explained by the strong emphasis on effective pH control (Götze et al., 2020).Other amino acids such as isoleucine, serine and valine are highly concentrated in gills due to their importance in protein synthesis and degradation under osmotic stress (Haider et al., 2019).Trigonelline was also found in high concentrations and is known as a plant alkaloid, but has also been found in a variety of marine invertebrates, including bivalves (Cappello et al., 2018;Poulin et al., 2018;Götze et al., 2020;Frizzo et al., 2021;Zhou et al., 2023).The exact function of alkaloids in animals is not yet clear, but trigonelline is supposed to have various biological functions.In crustaceans urinary trigonelline was shown to serve as a chemical sensor for predator-prey interaction (Poulin et al., 2018) but it also to play a role in glucose metabolism and inflammation (see review of Anthoni et al., 1991;Mathur and Kamal, 2012;Moreno et al., 2022).In the hard clam, Mercenaria mercenaria, increased trigonelline levels were found after animals were exposed to hypersaline stress and were discussed as an oxidative stress response (Zhou et al., 2023).Similarly, after acute warming (to 26°C) Götze et al. (2020) found a change in trigonelline levels in gills of P. maximus when additionally exposed to acute hypoxia and to the combination of hypoxia and hypercapnia.Liao et al. (2019) also reported increased levels of trigonelline under oxidative stress in OA-exposed scallop, Patinopecten yessoensis.However, whether the high trigonelline content is related to its putative role in oxidative stress remains speculative and requires further investigation.Interestingly, the tonic muscle shows similar concentrations in "gill-specific" metabolites such as glycolate and methionine, but also a very prominent serine concentration.Since the metabolite matrix of the tonic muscle is only poorly understood we can only assume that such metabolites in the tonic muscle have a function similar to that proposed for gills.
The metabolite spectrum of the mantle tissue is quite comparable to that of the gill tissue, exemplified by high concentrations of osmotically relevant metabolites, which indicate a role in osmo-and ion regulation as expected for a tissue in direct contact to seawater.However, the mantle spectra exhibited a lower quality due to a broader line width of the NMR signals.The generally poorer quality of 1 H-NMR spectra from mantle tissue is a phenomenon also observed in other NMR studies on tissue extracts from marine bivalves (e.g.Lannig et al., 2010).This was explained by the functional heterogeneity of the mantle tissue, e.g. between inner and outer mantle (Clark, 2020).However, despite this heterogeneity, small n numbers and lower NMR signal quality, distinct differences in metabolite concentrations were observed to those of the other tissues.The concentrations of the amino acids threonine, glycine and valine were highest in the mantle tissue compared to the other tissues.Other amino acids such as alanine, leucine and isoleucine also showed high concentrations in the mantle tissue, indicating a particular importance of amino acids in the mantle tissue.The main functions of the mantle lie in biomineralization and calcification.A high proportion of free amino acids, provided for biomineralization is therefore not surprising.It has been shown only recently that amino acids like serine and threonine are of particular importance for biomineralization and calcification processes in bivalves such as scallops (Yarra et al., 2016, Yarra et al., 2021).
The spectra of the two muscle types also show a characteristic picture in their metabolite profiles.In particular, both are characterized by prominent concentrations of the amino acid Larginine.L-arginine is the precursor for PLA, which corresponds to phosphagen in vertebrate muscle.PLA is used as the main rapid energy reserve to provide ATP and for the formation of octopine in the muscle during exercise and recovery (Grieshaber, 1978).It is therefore the essential metabolite for spontaneous locomotion in the muscle of scallops and indeed, it was shown to vary depending on the lifestyle of different scallop species (Tremblay et al., 2012;Tremblay and Guderley, 2014).In addition, fumarate and carnitine could be detected in both muscle types.These metabolites are also important for energy provision in muscle.Fumarate (fumaric acid) is an intermediate in the Krebs cycle for the production of ATP, but is also involved in the biosynthesis of L-arginine.Carnitine is the shuttle for fatty acids into mitochondria where they are subsequently used for energy production through beta-oxidation.The relatively high content of both metabolites may indicate i) a pronounced biosynthesis of L-arginine and ii) a high proportion of betaoxidation in the adductor muscle.Interestingly, Guderley et al. demonstrated that neither palmitoyl carnitine nor aspartate were oxidized in isolated muscle mitochondria of the tropic scallop, Euvola ziczac (Guderley et al., 1995).Therefore, the question arises why such high carnitine levels were found in the adductor muscle of P. maximus if fatty acids are apparently not used for energy production in muscle mitochondria.Yet, the studies of Guderley et al. were done on a tropical species and the rate of substrate use in mitochondria may differ in tropical and boreal species.Temperature is indeed an important factor for the utilization of a specific substrate in mitochondria in all species (Barbe et al., 2023).In addition to Larginine, fumarate and carnitine, choline, glucose, homarine and some other amino acids also showed comparable concentrations in the two muscle types.However, there were also some fundamental differences.For example, acetate and the amino acids alanine, glutamine and glutamate were highest in phasic muscle, whereas in tonic muscle they were at similar levels as in gills.In tonic muscle ATP levels were lowest, but succinate, lactate and O-phosphocholine Plot of the metabolite-metabolite interaction network of phasic and tonic muscle tissue between control and long-term OW exposed P. maximus.The insert in the upper left corner shows the entire network, while the center presents a subnetwork reduced to metabolites with strong connections.Both ATP and L-arginine have a central role in the network.The interaction of L-arginine with ATP, ADP and AMP and with oxygen, carbon dioxide and sodium is marked by bigger blue circles.The filter for nodes based on degree were set to 15 and a filter for nodes based on betweenness were set to 5. levels were higher than those in other tissues.These differences can be explained by the specific function of each muscle type.The scallops' tonic muscle, the so-called slow-catch "smooth muscle", is mainly responsible for shell closure and prolonged contractures with little energy expenditure (Chantler, 2006;Sun et al., 2018).The phasic or glycolytic muscle is a fast twitching, striated muscle and mainly used during swimming.In line with a higher mitochondrial content in invertebrates' tonic compared to phasic muscle (Nguyen et al., 1997) the observed lower ATP content and elevated succinate levels in tonic muscle appear conceivable.O-phosphocholine is an intermediate for phosphatidylcholines, the key building blocks for cellular membranes, and its observed high levels might indicate that the tonic muscle has a high membrane conversion rate due to its higher mitochondrial content.It is also conceivable that the higher mitochondrial content of tonic muscle needs a higher cell membrane flexibility, as suggested in the theory of membranes as metabolic pacemaker (Hulbert and Else, 1999).The membrane pacemaker theory (MPT) developed for mammals has also been tested in bivalves, although the general applicability to bivalves may differ by species or phenotype and is not always obvious (Pernet et al., 2006(Pernet et al., , 2007;;Sukhotin et al., 2017).However, the differences between individual MPT studies on molluscs could also be due to tissue type, lifestyle or different environmental histories of the species studied (Guderley et al., 2009).Further experiments would be necessary to elucidate the role of MPT for both muscle tissues in bivalves.
Finally, the tonic muscle showed the lowest taurine, TMAO and trigonelline concentrations between tissues.In addition, low concentrations were observed also in other osmotically relevant metabolites such as homarine, glycine and glucose, which were comparable to the phasic muscle.This simply shows that involvement in systemic osmoregulation is only of secondary importance for the tonic and partly also for the phasic muscle.For the amino acids the tonic muscle exhibited highest concentrations of serine and tyrosine, but lowest in alanine, glycine, leucine and threonine concentrations in comparison to the other tissues.The question, however, why these two amino acids were elevated in the tonic muscle and the others were deficient cannot be answered here.

Tissue-specific OA, OW and OWA related metabolite changes
Metabolite profiles after long-term acclimation to ocean acidification (OA), ocean warming (OW) and the combination of ocean warming and acidification (OWA) are tissue-specific.In gills, ATP levels decreased in the two warming groups indicating increased ATP demand in the warmth.Indeed, Schalkhausser observed an increased oxygen consumption after long-term acclimation (Schalkhausser et al., 2014) which may be mirrored in the lowering of ATP levels.Compared to other tissues, gills are an aerobic-active and mitochondria-rich tissue with a nearly 2-fold higher metabolism compared to, for example, mantle tissue (Stapp et al., 2017).In line with increased gill respiration rates determined in the blue mussel after long-term exposure to intermediate and high OA scenarios, the increased gill energy demand can be explained by increased cilia movements involved in water and food transport (Stapp et al., 2017).
The phasic muscle showed only a few changes in metabolite concentrations under OA but more prominent changes in both warming groups, OW and OWA.Glycolate, lactate, Ophosphocholine and threonine were all decreased in the phasic muscle of the OA group.The amino acids alanine, glutamine, glutamate, arginine and leucine, as well as acetate and succinylacetone also displayed lower concentrations in the OW and OWA groups.The observed decrease in O-phosphocholine under OA can be explained by its specific function as an essential constituent for membranes and their turnover.An increased turnover of cell membranes as visible in concentration changes of e.g.phosphocholine and its derivatives (also proposed by Rebelein et al., 2018 andGötze et al., 2020) could be induced by warming.In addition, bivalves are osmoconformers and OA as well as warming, similar to salinity fluctuations, were shown to induce osmotic changes in bivalve tissues (e.g.Morabito et al., 2013 for OA;Jiang et al., 2020;Georgoulis et al., 2022 for warming).Since muscle tissues have a lower content of organic osmolytes than gills (see above), the observed changes of phosphocholine in phasic muscle might indicate a remodeling of the cell membrane to compensate for cell volume changes and ensure cellular homeostasis.This is also indicated by the network exploration of phasic muscle under OA, in which the cell membrane-specific metabolite phosphorylcholine is at the center of the network together with lactate (see Supplementary).At the same temperature, a decrease in lactate levels might indicate metabolic depression under OA.A connection between lactate and metabolic depression was already suggested by Pörtner et al. (2014) in 1994.A depressive effect by the OA induced reduction in extracellular pH has been intensively discussed for different bivalve species (e.g.Thomsen and Melzner, 2010;Lesser, 2016 for blue mussels; Lannig et al., 2010 for oyster; Liu and He, 2012 for mussels, oysters and scallops).On the one hand, Liu and He (2012) reported decreased MO 2 levels under OA in the scallop Chlamys nobilis suggesting metabolic depression for scallops.On the other, one of our earlier organismal studies did not observe decreased MO 2 levels in P. maximus under OA, which argues against metabolic depression under the selected OA conditions for this species (Schalkhausser et al., 2014).The degree of metabolic depression will clearly be influenced by warming.In addition, lactate may be used as an extra energy source via the pyruvate shuttle under routine metabolism under OA conditions, when P e O 2 levels are reduced in the hemolymph (Schalkhausser et al., 2014).However, it must be noted that the prominent NMR signals of lactate and threonine are close to each other in the NMR spectrum.In our semi-targeted binning analysis approach, it cannot be ruled out that nearby NMR signals may influence each other due to overlap, resulting in a misleading profile.The similar pattern of the two metabolites between the four groups could indicate this, so the results for lactate and threonine should be interpreted with caution.Indeed, threonine changes are not as expected for a proteinogenic amino acid, where one would expect that it changes in parallel to other proteinogenic amino acids between treatments, which is not the case.
The changes in amino acid levels under OW and OWA can all be explained by the increasing energy demand due to the higher temperatures, as all changing amino acids are involved to a greater or lesser extent in energy metabolism.In the particular case of Larginine, the decrease in phasic muscle can be explained by PLA reduction and octopine formation.As shown above PLA is the main energy source for the fast-twitching phasic muscle to refuel ATP during swimming and escape responses.In our previous study, using in vivo 31 P-NMR spectroscopy we revealed a PLA decrease in resting phasic muscle of P. maximus following long term exposure to warming and OWA (Bock et al., 2019).In the present study, we obtained metabolic profiles via 1 H-HRMAS NMR spectroscopy.As in the case of lactate and threonine, the 1 H-NMR signals of Larginine and PLA are very close to each other in the 1 H-NMR spectrum, therefore it was not possible to distinguish between PLA and L-arginine.The decreased L-arginine content under OW and OWA conditions may reflect the warming-induced increase in energy demand and dephosphorylation of PLA to fuel ATP, if octopine formation absorbed L-arginine.Unfortunately, we have only rudimentary information on the possible 1 H-NMR signals of octopine and have not found any data on octopine in the literature and databases.Therefore, an increase of octopine from L-arginine decomposition could not be verified in the NMR spectra from the phasic muscle of the OW and OWA groups.The decrease in acetate in the two warming groups could indicate that acetate is metabolized to fatty acids in order to obtain the required energy through beta-oxidation or to synthesize phospholipids (Zhukova, 1986) for the warming-induced rebuilding of cell membranes.
Finally, succinylacetone (SA) decreased in the phasic muscle under OW and OWA.This is a finding that to our knowledge has not been reported in the literature.SA is a metabolite in tyrosine metabolism and has signaling effects (Bechara et al., 2007).In vertebrates SA can inhibit the biosynthesis pathway of heme groups in mitochondria and it also can release iron in the liver and brain (Bechara et al., 2007).More importantly, SA can activate the biosynthesis of glutathione peroxidase and SOD (Rocha et al., 2000;Bechara et al., 2007).The production of reactive oxygen species is a well-known phenomenon in isolated mitochondria under heat stress, also in bivalves (e.g.Heise et al., 2003;Sokolova, 2023) and its importance was recently shown for the yesso scallop under OA (Liao et al., 2019).The observed increase of glycolate in phasic muscle of both warming groups would be in line with a strengthened antioxidative function of glycolate, compensating for the decrease in SA under these conditions.As described for gills the function of glycolate in marine animals is unclear and needs further investigations.

Link to organismal performance
The review of Guderley and Pörtner which focuses on scallops and fish, highlights the crucial role of energy availability under climate change for the aerobic power budget in marine ectotherms (Guderley and Pörtner, 2010).Energy availability depends on energy resources via food and/or stored energy reserves such as carbohydrates and lipids which are processed via aerobic and anaerobic metabolic pathways.Thereby, thermally induced hypoxemia which results from the mismatch between rising oxygen demand and increasingly limited supply via uptake (ventilatory capacity) and transport (cardiac circulatory capacity) is a main driver for temperature induced limitation in ectotherm performance (Pörtner, 2002, Pörtner andLannig, 2009;Sokolova et al., 2012;Pörtner et al., 2017).Present observations in tissuespecific responses may indicate coordinated tissue changes balancing functional constraints for the sake of organism survival (Lannig et al., 2010;Tripp-Valdez et al., 2017;Götze et al., 2020).
In the following, we aim to interpret our previous results of P. maximus' swimming performance and oxidative capacity under ocean warming and acidification (Schalkhausser et al., 2014;Bock et al., 2019) together with the present metabolomics findings in tissue samples.The increased energy demand in ocean warming treatments indicated by elevated resting metabolic rates of P. maximus at 20°C compared to 10°C (Schalkhausser et al., 2014) is accompanied by and correlates with the observed decrease of ATP levels in gill tissue (see above).The decreased levels of L-arginine in phasic muscle of OW-and OWA exposed scallops may have followed the lowered levels of PLA (Bock et al., 2019) including absorption of L-arginine in octopine formation (Grieshaber, 1978).In addition, we showed a mainly OW-induced decrease in the phasic muscle's maximum oxidative flux of OW and OWA exposed P. maximus as well as the delayed replenishing of the phasic muscle's PLA pool in the warmth (Bock et al., 2019).All of this can be linked to the warming induced strong reduction in haemolymph P e O 2 levels of OW-and OWA-exposed P. maximus (Schalkhausser et al., 2014) indicating that the warming induced hypoxemia does not allow for maximal and fast recovery of PLA levels, since PLA replenishment occurs upon aerobic recovery (Grieshaber, 1978;Kamp, 1993;Sokolova et al., 2000).Most metabolite changes in the warming groups were observed in the phasic muscle, which were mainly energy related, as confirmed by the pathway analysis and network exploration.A change in energy demand via the TCA cycle was also reported for the adductor muscle of Argopecten irradians under thermal stress (Song et al., 2022).In line with our conclusion, the authors suggested that stored energy reserves were mobilized under thermal stress.In addition, Song et al. reported the degradation of valine, leucine and isoleucine for A. irradians under thermal stress (Song et al., 2022), similar to our observations.According to the pathway analysis and network exploration, L-arginine plays a central role in the metabolite network for energy provision under long term acclimation to warming in the adductor muscle, which is decoupled from the impact of OA.In terms of energy metabolism, this is not surprising, as L-arginine-phosphate is the basis for providing phosphate groups to ADP for building ATP.However, its direct connection to oxygen as determined from the network analysis also shows its dependence on oxygen and supports our hypothesis that PLA cannot be sufficiently rebuilt due to the reduced oxygen content in the hemolymph.This in turn leads to a reduced support for the increased ATP demand under OW.Interestingly, the network exploration for the warming groups also displayed a direct connection of L-arginine and ATP to methionine, sodium and CO 2 .Furthermore, both muscle types exhibited significantly higher concentrations of carnitine in comparison to gills and mantle (see Table 1).
We therefore propose the following model of metabolic connectivity between the two muscle types, in accordance with the work by Peŕez et al. (2008), who demonstrated the recuperation of the phasic muscle by the tonic muscle during swimming in scallops.Compared to the phasic muscle, for which energy must be provided immediately during swimming, the tonic muscle uses energy rather slowly and for a prolonged period.To accomplish their different tasks the phasic muscle requires fast ATP supply under performance and essentially uses PLA for ATP generation.Only during prolonged exercise at low speeds does it switch to aerobic utilization of carbohydrates and fatty acids to supply ATP, as is generally the case for contracting muscles (Hargreaves and Spriet, 2020).The tonic muscle provides ATP for a prolonged period of performance and therefore preferably utilizes beta-oxidation to generate ATP.This is reasonable, as carbohydrate oxidation has a higher rate of ATP synthesis than fatty acid oxidation; the latter, however, can obtain more energy per molecule in the long term (overall ATP generation of 130 ATP vs 36 ATP per molecule; Hargreaves and Spriet, 2020).Fatty acids are used as the predominant substrate, for example during fasting and hypoglycemia, to provide energy via beta-oxidation in the muscle (Longo et al., 2016;Talley and Mohiuddin, 2023).It is therefore quite conceivable that P. maximus make increased use of beta-oxidation for energy supply in times of increased energy requirements such as under OW and exercise in the tonic muscle.This view is supported by the increase of carnitine in the phasic and tonic muscle in comparison to control (see Figure 4), albeit not significant, and the high concentration of methionine in the tonic muscle (Table 1).The higher carnitine concentration may allow a higher shuttle of fatty acids into tonic muscle mitochondria and thus an increased beta-oxidation (see above).The connection to sodium and methionine in the network of L-arginine might confirm this, since sodium is essential for the transport of carnitine across membranes (Talley and Mohiuddin, 2023) and methionine, from which carnitine can be synthesized endogenously (Longo et al., 2016).The link to CO 2 might indicate the importance of beta-oxidation for energy provision in muscles under unfavorable conditions (Longo et al., 2016), such as OA, OW and the combination of both.
The decreased levels of proteinogenic amino acids in both warming groups in the adductor muscle (Figure 2) might indicate a switch to the use of proteins for energy production as also suggested by (Götze et al., 2020), since carbohydrates are limited.An ongoing depletion of carbohydrates during prolonged exercise together with lowered oxygen levels in the haemolymph would lead to a decrease in power output (Hargreaves and Spriet, 2020).Indeed, a decreased clapping force was observed in the warming groups, while the number of claps remained the same compared to control (Schalkhausser et al., 2014).Since the carbohydrate reserves in the phasic muscle of the warming groups (OW & OWA), just like the PLA level (Bock et al., 2019), are already reduced due to the higher energetic costs in the warmth, there is an earlier reliance on the tonic muscle during swimming (Peŕez et al., 2008).Betaoxidation, with its higher capacity but lower power output, can thus maintain the number of claps, but at the expense of a reduced clapping force.In addition, beta-oxidation can provide the energy during recovery from exercise in mammals (Hargreaves and Spriet, 2020), but the processes and pathways that metabolize fatty acids to provide ATP are much slower than carbohydrate pathways.If one assumes that due to the lower P e O 2 levels in the haemolymph, the energy is rather provided by beta-oxidation and octopine is formed from L-arginine (Grieshaber, 1978), it becomes understandable why recovery after fatigue was prolonged in the scallops of the warming groups (OW & OWA) compared to control.To summarize, if the tonic muscle compensates for the reduced performance of the phasic muscle and if beta-oxidation is favored for energy production under systemic hypoxemia during OW and OWA, the observed organismic limitations of the scallops from the previous Changes in carnitine concentrations in phasic (left) and tonic muscle (right) tissue between treatments.Note, the differences are not significant.Bock et al. 10.3389/fevo.2024.1347160Frontiers in Ecology and Evolution frontiersin.orgstudies can be explained.However, additional studies on the transport and flow of substrates fueling ATP production in the two muscle types will be needed to test this.Future experiments must show to what extent the energy is distributed in the body of the scallops under changing environmental conditions and whether additional energy requirements such as gonadal development will lead to a reduction in the energy budget and/or changes in energy distribution.This could have further consequences for the performance level of scallops from coastal and intertidal areas (Feidantsis et al., 2020).

Conclusions
The presented study of metabolite profiles of individual tissues showed distinct metabolite profiles, which are adjusted to the specific tasks of the organs.After long-term incubation under OA, OW and their combination OWA, metabolic profiles were largely different within the adductor muscle, between phasic and tonic muscles.These changes in metabolite concentrations mainly involved metabolic adjustments essential for the increased energy demand induced by warming and the associated hypoxemia.The proposed model for the interaction of the two muscle types of the adductor muscle can explain the observed performance limitations of scallops under OW and OWA.The tonic muscle apparently pursues a different strategy for energy provision than the phasic muscle, which for strong clapping essentially relies on PLA hydrolysis.This interplay is shifted to lower performance rates by organismic changes such as oxygen reduction in the hemolymph, especially under warming, and thus leads to the observed limitations in the swimming behavior of P. maximus.
3D views of a PCA (A) and a PLS-DA (B) across the metabolite profiles of all four tissue types (gills = red; mantle = green; phasic muscle = blue and tonic muscle = cyan) of the king scallop, P. maximus exposed to control conditions (normocapnia at 10°C).The unsupervised PCA showed a clear separation of all four tissues and the supervised PLS-DA resulted in a distinct classification between tissues.The ellipsoids indicate the respective confidence areas (not available for mantle n=3).

TABLE 1
Summary of the results of the ANOVA test presenting the significant differences in metabolite concentrations between tissue types (gills (C_G), mantle (C_M), phasic (C_P) and tonic muscle (C_T)) of the king scallop, P. maximus exposed to control conditions (normocapnia at 10°C).

TABLE 1 Continued
Bock et al.  10.3389/fevo.2024.1347160Frontiers in Ecology and Evolution frontiersin.orgexposedscallopsidentified changes in four important pathways.All of them are directly related to arginine and energy metabolism.A Network (NW) exploration of metabolite-metabolite interactions is presented in Figure3.The entire network consisted of 647 nodes, 2000 edges and 26 seeds (see insert on the left corner of Figure

TABLE 2
Major metabolic pathways identified by PW analysis contributing to the metabolic response of the adductor muscle (phasic and tonic muscle) of P. maximus to long-term exposures to ocean warming (OW: 20°C, 0.112 kPa PCO 2 ).