Behavioural Stress Propagation in Benthic Invertebrates Caused by Acute pH Drop-Induced Metabolites

Studies on pH stress in marine animals typically focus on direct or species-specific aspects. We here test the hypothesis that a drop to pH = 7.6 indirectly affects the intra- and interspecific interactions of benthic invertebrates by means of chemical communication. We recorded fitness-relevant behaviours of small hermit crabs Diogenes pugilator, green shore crabs Carcinus maenas, and harbour ragworms Hediste diversicolor in response to short-term pH drop, and to putative stress metabolites released by conspecifics or gilt-head sea bream Sparus aurata during 30 min of acute pH drop. Not only did acute pH drop itself impair time to find a food cue in small hermit crabs and burrowing in harbour ragworms, but similar effects were observed under exposure to pH drop-induced stress metabolites. Stress metabolites from S. aurata, but not its regular control metabolites, also induced avoidance responses in all recipient species. Here, we confirm that a short-term abrupt pH drop, an abiotic stressor, has the capacity to trigger the release of metabolites which induce behavioural responses in conspecific and heterospecific individuals, which can be interpreted as a behavioural cost. Our findings that stress responses can be indirectly propagated through means of chemical communication warrant further research to confirm the effect size of the behavioural impairments caused by stress metabolites and to characterise their chemical nature.


INTRODUCTION
Compared to the open ocean, coastal areas, and particularly intertidal zones are highly variable environments characterised by abrupt changes in water parameters. This includes fluctuations in pH beyond 0.3 units, the levels predicted for average change related to the process of ocean acidification toward the end of the century (Caldeira and Wickett, 2003;Sabine et al., 2004;IPCC, 2019). While this may mean that intertidal species are more resilient to climate change due to their acquired tolerances of pH fluctuations, it also forces them to live more frequently near their physiological tolerance limits (Truchot and Duhamel-Jouve, 1980;Hofmann et al., 2011;Sokolova, 2018;Wolfe et al., 2020). Superimposing ocean acidification on natural pH fluctuations may further increase the variability of organisms' responses (Eriander et al., 2015). Therefore, organisms inhabiting intertidal areas are interesting models to study the effects of short-term pH fluctuations within the IPCC predicted range. A low environmental pH can directly alter animal behaviour through several pathways, which include (i) deviation of energy budgets toward the stress response (Pörtner, 2008), (ii) fleeing to avoid the sources of stress (Pörtner and Peck, 2010;Abreu et al., 2016), (iii) disrupted information detection and processing leading to impaired decision-making Porteus et al., 2018;Thomas et al., 2020;Cothran et al., 2021), and (iv) alteration of the chemical signals themselves impacting the sensory environment and the transfer of information (Wyatt et al., 2014;Roggatz et al., 2016Roggatz et al., , 2019. Behavioural effects triggered by lowered pH are known to occur in different taxonomic groups such as crustaceans (de la Haye et al., 2011), marine ragworms (Bond, 2018), and fish (Munday et al., 2009), although recent research debates both their ubiquitousness and effect size (Clark et al., 2020a;Clements et al., 2020a,b, but see Clark et al., 2020b;Munday et al., 2020;Williamson et al., 2021 for the importance of which comparisons are made). In addition, it should be mentioned that although there is little literature disentangling the behavioural effects of pH drop vs. high CO2, an increasing body of recent research evidenced the role of pH-dependent altered chemical communication (Roggatz et al., 2016;Schirrmacher et al., 2020;Porteus et al., 2021;Velez et al., 2021).
In aquatic environments, where visibility can be limited, infochemicals and chemosensory functions, such as pheromones used for mating (Karlson and Lüscher, 1959) are crucial for communication (Hardege, 1999;Jordão and Volpato, 2000;Ashur et al., 2017;Saha et al., 2019). However, infochemicals are also a mechanism to signal stress-but stress communication has mostly been researched within the context of biotic stressors. For example, alarm cues released following mechanical damage from fish epidermal club cells can trigger panic reactions in conspecifics and heterospecifics (Toa et al., 2004;Júnior et al., 2010). Such mechanisms of chemical communication also occur in invertebrates, as evidenced by the reduced out-ofburrow activity of the marine polychaete Alitta virens exposed to damaged conspecifics (Watson et al., 2005;Ende et al., 2017). Disturbance cues refer to chemicals that may be stored in gill epithelium or urine, and are released voluntarily or involuntarily by disturbed or stressed prey, to induce early antipredator responses in recipients to anticipate potential threats (Bairos-Novak et al., 2017;Goldman et al., 2020a). The central role of chemical communication in the aquatic environment and the recent evidence of its impairment by ocean acidification pinpoints the need for a deeper understanding of such potential environmental modulation of chemical signalling (Chivers et al., 2013).
Overall, though, it is not well understood yet whether abiotic stressors such as pH drop can also induce the release of chemical cues and whether these can propagate the stress response to other individuals. For such cues, we here introduce the term "stress metabolites" as defining such compounds released, voluntarily or not, by an organism in response to abiotic stressors. These can be detected and processed by conspecifics and/or heterospecifics, leading to the induction of a behavioural stress response (Hazlett, 1985). Of note, Draper and Weissburg (2019) defined the "direct effects" of elevated CO 2 on predator-prey interactions as the effects modifying the sensory pathways (cue production and emission by the donor and cue reception and processing by the receiver). Here, we term the effects of acute pH drop as "direct" when the animal is exposed to the stressor itself and as "indirect" when the animal is exposed to the stress metabolites released upon pH stress. Both direct and indirect effects may be guided by sensory pathways as well as metabolism in senders and receivers (as defined in Draper and Weissburg, 2019). We term "stress propagation" the signalling of stress responses within or between species by means of stress metabolites ( Figure 1A). Detecting stress metabolites may allow other individuals to modify their behaviour to avoid a change in the abiotic environment, and communities to coordinate or potentiate their response to ensure survival (Mothersill et al., 2007;Giacomini et al., 2015;Abreu et al., 2016). However, to investigate the function of any pH drop-induced chemical cues, the original stressor (pH drop) must first be removed from the experimental design.
In this study, we investigate, within and among species inhabiting the intertidal zone, the indirect effects of acute pH drop through stress metabolites induced by it on fitness-relevant behaviours. Using a full factorial design, we test the hypothesis that an acute pH drop to end-of-century level (7.6) will both directly (i.e., when the animal is exposed to the abiotic stressor) and indirectly (i.e., when the animal is exposed to the stress metabolites) affect behaviour through the induction of stress metabolites. We observe three marine invertebrates (small hermit crab Diogenes pugilator, green shore crab Carcinus maenas, and harbour ragworm Hediste diversicolor) exposed to control pH and pH drop in combination with conditioned water from both conspecifics and their potential vertebrate predator, the gilt-head sea bream S. aurata.

Experimental Design
Small hermit crabs Diogenes pugilator and green shore crabs Carcinus maenas were collected in autumn at low tide in the Ria Formosa lagoon (Portugal). Harbour ragworms (Hediste diversicolor) were sourced via a local supplier (Valbaits) from the Setubal lagoon (Portugal). Animals were acclimated for 1 week at pH 8.1 (pH is reported using National Bureau of Standards scale) in large flow-through open-circuit communal tanks (continuous flow of 0.05 L/min) mimicking natural habitats at the Ramalhete Marine Station (CCMAR, Faro, Portugal) in a direct CO 2controlled system (described in Sordo et al., 2016Sordo et al., , 2018Gregório et al., 2019) with pCO 2 constantly measured and adjusted. Two system seawater conditions were used, either control pH of 8.1 ± 0.015 or pH drop of 7.6 ± 0.008. For a measured total alkalinity of 2,500 µmol/kg SW, the pCO 2 was calculated as 537 and 1,922 µatm in the control and pH drop system water, respectively, according to the CO2calc software (v1.2.0, Robbins et al., 2010). pCO 2 was not estimated in any treatment containing reef buffer, HCl, or dilutions, as total alkalinity was not measured . Behavioural effects of short-term acute pH drop and stress metabolites it induced were tested in three species: Diogenes pugilator, Carcinus maenas, and Hediste diversicolor. Metabolites were obtained by conditioning donors in either control pH = 8.1 (releasing putative control metabolites), or pH drop = 7.6 (releasing putative stress metabolites), followed by pH adjustments for factorial design yielding four experimental conditions CM, pH drop, SM, and pH drop + SM. Predictors are binary coded as 0 (control metabolites, control pH) and 1 (stress metabolites, pH drop). Metabolites originated from conspecifics or the heterospecific Sparus aurata. Behaviour assays for crabs consisted in locating a feeding cue in 300 s. Specimens of H. diversicolor were placed on top of sediment and burrowing behaviour was recorded for 300 s. A range of avoidance behaviours were also recorded. Experimental conditions were CM, control metabolites in control pH; SM, control metabolites in control pH; pH drop, control metabolites in pH drop; pH drop + SM, stress metabolites in pH drop. Animal drawings by A. Murcia and KCWV. in those, and pCO 2 could therefore have differed from that in the system water (see Supplementary Information for details). Seawater parameters were measured daily at 2.30 pm (mean temperature: 20.16 • C ± 1.05 • C, mean salinity: 35.67 PSU ± 0.26 PSU, mean dissolved oxygen: 7.55 mg/L ± 0.16 mg/L). A total pool of approx. 100-120 D. pugilator, 150-200 C. maenas, and over 150 H. diversicolor was used during the experimental period, yielding a total of 921 observations. Because no individual marking could be achieved, we aimed to space out the reuse of animals in successive behavioural assays by circulating animals through recovery tanks afterward, before returning them to their stock population tanks once a day. Animal reuse was randomised across treatments and treatments randomised per day to prevent any confounding effects on the measured behaviours. All experiments were conducted in October 2019, except one independent set of observations (n = 40) in C. maenas conducted in September 2018, which was analysed together with the 2019 data with the year treated as a covariate.
A three-way factorial design of pH drop × stress metabolites × donor species was generated to study the effects of acute pH drop, stress metabolites induced by it, and their combination vs. a control containing regular metabolites (control metabolites) ( Figure 1B). Metabolite donor species were D. pugilator, C. maenas, H. diversicolor (called conspecific donors), and the potential predator Sparus aurata (called heterospecific donor). Metabolite donors were conditioned for 30 min in seawater at regular pH (pH 8.1, 538 µatm CO 2 , putatively inducing control metabolite release), or pH drop (pH 7.6, 1,922 µatm CO 2 , putatively inducing stress metabolite release). Recipient species were D. pugilator, C. maenas, and H. diversicolor which received water containing either conspecific or heterospecific metabolites. To achieve the full factorial design, recipient species received either control or stress metabolites and were tested in either control pH (8.1) or pH drop (7.6), by re-adjusting the pH of the conditioned water before each bioassay (with the pH verified with an accuracy of 0.02 units using a portable pH metre). Overall, four different experimental treatments were obtained: (i) putative control metabolites at pH 8.1 (CM), (ii) putative control metabolites at pH 7.6 (pH drop), (iii) putative stress metabolites at pH 8.1 (SM), (iv) and putative stress metabolites at pH 7.6 (pH drop + SM). Testing conspecific and heterospecific metabolites with three recipient species amounted to a total of 12 treatments. Prior to use in recipient species, conditioned seawater from H. diversicolor and S. aurata (0.0067 fish/L) was 10-fold diluted in fresh system seawater at the desired pH. To explore the possibility that stress metabolites are equivalent to highly concentrated control metabolites, we additionally tested the effects of undiluted S. aurata control metabolites at pH 8.

Behavioural Bioassays
For each of the four conditions, we measured the time to find a feeding cue (1/10 diluted mussel juice) vs. a mock (seawater) cue in D. pugilator and C. maenas, or to bury the FIGURE 2 | Effects of predictors (pH drop, stress metabolites, donor) on time-to-success in Diogenes pugilator, Carcinus maenas, and Hediste diversicolor. (A) Effects of predictors on success ratio (aka hazard ratio with success as event, arbitrary units). Likelihood ratio tests for overall model fit were: (D. pugilator: P = 0.06; C. maenas: P = 0.009; H. diversicolor: P < 0.0001). Significant predictors from Cox proportional hazard models are shown with asterisks, and significant interaction terms are shown in grey above the plots. (B) Interaction term plot of marginal effects showing predicted success score (aka risk score with success as event, arbitrary units ± confidence interval) split by donor and recipient species. Crossed solid grey lines represent an interacting effect of pH and metabolites. (C) Kaplan-Meier curves to visualise success probability (cumulative event) for each experimental condition over time. *P ≤ 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001. Experimental treatments were: CM, control metabolites in control pH; pH drop, control metabolites in pH drop; SM, stress metabolites in control pH; pH drop + SM, stress metabolites in pH drop.
head entirely in the sand in H. diversicolor, and dubbed this variable "time-to-success." Behavioural assays were terminated once the animal successfully grabbed the ballasted sponge containing the feeding cue with their pincers (D. pugilator and C. maenas) or buried its entire head (H. diversicolor), or at a maximum time of 300 s. Both feeding and burrowing behaviours are tested response variables in crabs Wang et al., 2018) and ragworm (Bhuiyan et al., 2021) exposed to pH drop. Additionally, avoidance behaviour was binary coded. In D. pugilator and C. maenas, these avoidance behaviours included freezing (suddenly retreating into the shell or attempting to burrow, sudden and lasting arrest in locomotion for at least 120 consecutive seconds or more than 150 s in total), and escaping (walking along the walls or retreating into a corner of the tank for at least 120 consecutive seconds or more than 150 s in total). Such freeze and escape behaviours indicate danger in crustaceans (Katz and Rittschof, 1993;Perrot-Minnot et al., 2017;Tomsic et al., 2017). In H. diversicolor, avoidance behaviours consisted of freezing, which might be accompanied by spread of jaws, sideways-undulating behaviour, and formation of a noticeable slime cap whilst outside the sediment. These behaviours indicate extreme stress in marine polychaetes (Mouneyrac et al., 2003;Burlinson and Lawrence, 2007;McBriarty et al., 2018). The feeding cues were randomly assigned left or right to the crabs. See Supplementary Methods for further details on behavioural assays. All data is given as Supplementary Dataset.

Statistical Analysis
Avoidance behaviours were analysed using generalised linear models for logit regression with binomial distribution using the stats R package (R Core Team, 2020). Time-to-success was TABLE 1 | Results of Cox proportional hazard model for the main effects of predictors (pH drop, stress metabolites, donor) on the time-to-success analysis in small hermit crab (Diogenes pugilator, n = 320 observations with 120 events of finding a food cue, degrees of freedom = 312), green shore crab (Carcinus maenas, n = 189 observations with 122 events of finding a food cue, degrees of freedom = 181), and harbour ragworm (Hediste diversicolor, n = 325 observations with 234 events of burrowing head in sediment, degrees of freedom = 317).

Predictors
Success ( modelled for the time to reach a feeding cue in D. pugilator and C. maenas or to burrow the head in H. diversicolor depending on the three predictors: pH drop, stress metabolites, and donor. Time-to-success data was analysed with a time-to-event analysis (also called survival analysis) using Cox proportional hazard models from the survival R package (Therneau and Grambsch, 2000;Therneau, 2021). The event being the success of food cue location (in crabs) or burrowing the head in sediment (in worms), the analysis was hereafter referred to as "time-to-success analysis." The exponentiated estimates (hazard ratios) from the Cox proportional hazard models were expressed as "success ratio." Success ratios were visualised using the plot_model function from the sjPlot R package (Lüdecke, 2021). The "success probability" (the probability of the feeding or burrowing event occurring at any time over 300 s) over time was represented by Kaplan-Meier curves drawn using the survminer R package (Kassambara et al., 2020). Animals not reaching the feeding cue or not burying their head were censored and assigned the maximum time of "300 +". All statistical models always included the main effects and interaction terms of the binary predictors: pH drop (pH = 8.1 vs. pH = 7.6), stress metabolites (stress metabolites vs. control metabolites), and donor (conspecific vs. heterospecific).
Covariates (year, number of water uses, crab size, where relevant; see section 2.2 in Supplementary Information) were excluded from models when deemed insignificant compared to the null model by at least two model selection methods among the following: an ANOVA Chi-square test from the stats R package (R Core Team, 2020), the Akaike Information Criterion ( AIC), and Bayesian Information Criterion ( BIC) from the AICcmodavg R package (Mazerolle, 2020). Overall model fit p-values were retrieved from the ANOVA Chi-squared test or the Likelihood ratio test for generalised linear models and Cox-proportional models, respectively. Pairwise comparisons between treatments involved in the interaction term between the three predictors were obtained from statistical models using the emmeans R package (Lenth, 2021), wherein false discovery rate p-value adjustments were applied for post hoc term-wide multiple testing.
To investigate the possibility that stress metabolites are equivalent to highly concentrated control metabolites (which were present in CM and pH drop treatments but in 10-fold dilution), additional tests compared the behavioural response of D. pugilator and H. diversicolor to 10-fold diluted stress metabolites and control metabolites to that of undiluted control metabolites from S. aurata (see Supplementary Information, section 2.4). The full dataset may include animals tested at most twice over the course of the experimental period. We therefore performed random bootstrapping of 50% of the dataset (Supplementary Information, section 2.5). All statistical analyses were conducted in RStudio (R Studio Team, 2020) with a significance threshold of P ≤ 0.05. See R script and Supplementary Information for further details.

Response of Small Hermit Crab Diogenes pugilator to pH, Stress Metabolite, and Donor Predictors
The time-to-success analysis was not conclusive for the stress metabolite predictor in D. pugilator (Figure 2A and Table 1, Z = −0.66, P = 0.5073). On the other hand, the Cox proportional hazard model (overall Likelihood ratio test model fit: P = 0.06) found that both pH drop (feeding success ratio = 0.47, Z = −2.14, P = 0.0325) and metabolites from S. aurata (feeding success ratio = 0.35, Z = −2.76, P = 0.0058) had significant negative effects on the feeding success ratio of D. pugilator. In addition, there was a significant interaction between pH drop and metabolite donor terms on the feeding success of D. pugilator (Z = 2.19, P = 0.0289). Overall, pH drop induced a significantly lower success (risk) score in the conspecific treatment ( Figure 2B). Post hoc analyses evidenced that all treatments (CM, pH drop, SM, and pH drop + SM) induced similar feeding time responses in both the conspecific and heterospecific donor groups (Supplementary Table 2).
Avoidance responses of D. pugilator did not depend on pH drop (Z = 0.25, P = 0.8049) nor on stress metabolites (Z = −0.25, P = 0.7995, overall Chi-square test model fit: P < 0.0001, Table 2), but were significantly more pronounced when metabolites originated from S. aurata (69%) instead of conspecifics (26%, Z = 2.46, P = 0.0139, Figure 3 and Table 2). After splitting the display of avoidance behaviour by donor, pairwise comparisons failed to find differences between the four treatments (CM, pH drop, SM, and pH drop + SM), but confirmed that the donor effect existed across all four treatments (Supplementary Table 3 and Figure 3). However, partitioning the avoidance behaviour of D. pugilator into freezing and escaping responses evidenced that escaping significantly increased with S. aurata control metabolites only when tested in pH drop (Supplementary Figures 1, 2

Response of Green Shore Crab Carcinus maenas to pH, Stress Metabolite, and Donor Predictors
The Cox proportional hazard model (overall Likelihood ratio test model fit: P = 0.009) showed that the feeding success ratio of C. maenas did not significantly vary with pH (Z = −0.80 P = 0.4247), metabolites (Z = 0.68, P = 0.4989), nor donor (Z = 1.09, P = 0.2773) terms (Figure 2 and Table 1). Nevertheless, the interaction terms were significant. Post hoc tests revealed that time-to-success responses were similar across treatments CM, pH drop, SM, and pH drop + SM when C. maenas received conspecific metabolites (Supplementary Table 8). This pattern changed when S. aurata was the metabolite donor as evidenced by significantly lower feeding success scores in C. maenas exposed to pH drop (Z = 2.97, P = 0.0179), SM (Z = 2.48, P = 0.0394), and the trend for pH drop + SM (Z = 2.20, P = 0.0563), compared to the control CM. In C. maenas, the predictors pH drop (Z = −1.23, P = 0.2195), stress metabolites (Z = −1.23, P = 0.2195), and donor (Z = −0.56, P = 0.5785) did not alter avoidance responses but their interaction terms were significant (overall Chi-square test model fit: P < 0.0001, Table 2 and Figure 3). Pairwise comparisons showed no differences between the four treatments (CM, pH drop, SM, pH drop + SM) in avoidance patterns of C. maenas receiving conspecific metabolites. Conversely, C. maenas facing S. aurata metabolites while experiencing treatments of pH drop (94%, Z = −3.60, P = 0.0020), SM (75%, Z = −3.28, P = 0.0031), and pH drop + SM (56%, Z = −2.33, P = 0.04) significantly increased their avoidance display compared to control CM (18%). Moreover, pH drop combined with stress metabolites (pH drop + SM treatment) instead of control metabolites (pH drop treatment) significantly lowered avoidance behaviours (Z = 2.21, P = 0.0405, Supplementary Table 9 and Figure 3). Partitioning the avoidance behaviour of C. maenas into freezing and escaping responses evidenced that escaping explained the avoidance response whereas no clear patterns were observed in the freezing behaviour (Supplementary Figures 1, 2 and Supplementary Tables 10-13, see Supplementary Results).

Response of Harbour Ragworm Hediste diversicolor to pH, Stress Metabolite, and Donor Predictors
Burrowing success responses of H. diversicolor did not depend on the donor species (P = 0.1093, Figure 2 and Table 1). However, the Cox proportional hazard model (overall Likelihood ratio test model fit: P < 0.0001) evidenced that both pH drop (burrowing success ratio = 0.30, Z = −3.91, P < 0.0001) and stress metabolites (burrowing success ratio = 0.35, Z = −3.42, P = 0.0006) terms significantly altered burrowing time responses of H. diversicolor. Significant interactions between the main predictors required to decipher the involvement of within-donor treatment effects (Supplementary Table 14). Post hoc tests revealed that pH drop (Z = 3.91, P = 0.0006), SM (Z = 3.43, P = 0.0012), and pH drop + SM (Z = 3.47, P = 0.0012) treatments significantly lowered burrowing success scores compared to the control CM in H. diversicolor exposed to conspecific metabolites. Moreover, H. diversicolor facing S. aurata metabolites had lower burrowing success scores in the pH drop + SM treatment compared to the control CM (Z = 3.35, P = 0.0049).
H. diversicolor displayed more avoidance in response to S. aurata metabolites (47%) compared to conspecific metabolites (20%, Z = 1.09, P = 0.0006, overall Chi-square test model fit: P < 0.0001, Table 2 and Figure 3). H. diversicolor exposed to S. aurata metabolites were not affected by pH drop (Z = −0.01, P = 0.9906) whereas stress metabolites induced significantly more avoidance than control metabolites (Z = 4.31, P < 0.0001, overall Chi-square test model fit: P < 0.0001, Table 2). There were no differences across pH drop, SM, and pH drop + SM treatments in H. diversicolor receiving conspecific metabolites (Supplementary Table 15). Pairwise comparisons revealed that the SM treatment (93%) triggered significantly more avoidance than CM (33%, Z = −4.31, P = 0.0001) and pH drop + SM treatments (51%, Z = 3.92, P = 0.0003). Specimens of H. diversicolor exposed to S. aurata metabolites displayed different avoidance responses in the control treatment CM (Supplementary Video 1) compared to pH drop (slowed-down burrowing response, Supplementary Video 2), SM (sideways-undulating behaviour with body flips, Supplementary Video 3), and pH drop + SM (raised head and spread mouth parts with freezing behaviour in Supplementary Video 4 and slowed-down movement with sideways undulating behaviour in Supplementary Video 5).
Importantly, the behavioural responses of H. diversicolor and D. pugilator remained unchanged even when we used undiluted control metabolites (see Supplementary Information, section  2.4, Supplementary Figure 3 and Supplementary Tables 16, 17). Of note, the avoidance response of H. diversicolor and D. pugilator faded as water was reused ( Supplementary Information, section  2.2, Supplementary Figure 4).
Furthermore, random bootstrapping of 50% of the dataset ( Supplementary Information, section 2

Direct Effects of pH Drop and Donor Type on Assay Success and Avoidance Behaviours
In this study, we performed over 900 behavioural observations, and our results show species-specific responses, donordependent effects, as well as sometimes large within-treatment variability. Sparus aurata-conditioned water, whether it contained control metabolites or pH stress-induced metabolites, elicited avoidance in all recipient species which can be attributed to antipredator behaviours. Despite the donor effect, however, much of the behavioural variance was explained by the effects of pH drop and stress metabolites, their interaction, or their interaction with the donor term (Figure 2).
We found that pH drop lowered success ratios in feeding and burrowing in Diogenes pugilator and Hediste diversicolor, respectively. Such low pH-induced impairments of fitnessrelevant behaviours were previously found in hermit crabs Pagurus bernhardus and P. tanneri as well as in the Japanese stone crab Charybdis japonica and may indicate impaired decisionmaking (de la Haye et al., 2011Kim et al., 2015;Wu et al., 2017). Past literature also reported delayed burrowing responses in the harbour ragworm H. diversicolor exposed for 28 days to pH drop (Bhuiyan et al., 2021), and an inverse relationship between burrowing and pH in the king ragworm Alitta virens (Batten and Bamber, 1996). These behavioural alterations may reflect the cost of acid-base regulation to counteract hypercapnia (Pörtner, 2008), as marine ragworms overexpress the acid-base regulator Carbonic Anhydrase (CA) at the expense of energy reserves in response to acidification (Freitas et al., 2016;Wage et al., 2016). In the longer term, behavioural costs of an altered burrowing activity may extend to crucial ecological functions as shown by the diminished bioturbation activity of H. diversicolor in similar conditions (Bond, 2018). In contrast, Carcinus maenas FIGURE 3 | Effects of predictors (pH drop, stress metabolites, donor) on the percentage of avoidance behaviour in Diogenes pugilator, Carcinus maenas, and Hediste diversicolor. Avoidance behaviour included freezing and escaping (D. pugilator and C. maenas), or freezing, curling, flipping, and slime secretion (H. diversicolor). (A) Split bars represent the presence (dark area) or absence (light area) of avoidance behaviours. Significant main predictors (and their interaction terms) and covariates are shown above plots. Significant pairwise comparisons between treatments in each donor/recipient are shown as horizontal grey lines. *P ≤ 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001. Experimental treatments were: CM, control metabolites in control pH; pH drop, control metabolites in pH drop; SM, stress metabolites in control pH; pH drop + SM, stress metabolites in pH drop. (B) Interaction plot showing the marginal effects on the predicted avoidance percentages (±confidence interval) of stress metabolites and pH drop within the conspecific and heterospecific groups for each species. Crossing solid grey lines represent an interacting effect of pH and metabolites for each donor. The interactive effects of predictors for H. diversicolor/H. diversicolor are not represented due to the missing CM treatment. Overall effects of donor and water use across the response of H. diversicolor (for either donor) are shown with solid lines, whereas the effects of stress metabolites are for S. aurata/H. diversicolor only (grey dashed line).
did not delay its feeding response to sudden pH drop in presence of conspecific metabolites, which is in line with meta-analyses showing that some aquatic arthropods are more robust to ocean acidification conditions (Kroeker et al., 2013;Wittmann and Pörtner, 2013;Clements and Darrow, 2018). The resilience to pH stress in C. maenas may be due to their pH regulatory capacities allowing them to maintain feeding during short-term assays, but they may show altered responses in longer-term assaysfor instance after 10 weeks-as their energetic requirements increase (Appelhans et al., 2012;Clements and Darrow, 2018). Alternatively, the fact that C. maenas did not respond to pH drop may arise from differences in the type of cues that are presented or that not all behaviours are similarly affected by pH drop (Richardson et al., 2021). Short-term hypercapnia does not seem to be a strong stressor in C. maenas, due to its pH-variable intertidal habitat (Fehsenfeld et al., 2011). Other than reduced success ratios, pH drop itself did not induce any avoidance behaviours indicative of stress.
Besides the direct effect of pH drop on behaviour, previous studies on ocean acidification also found contrasting effects of pH drop on predator-prey interactions (Draper and Weissburg, 2019). Here, we observed an increased escaping tendency of C. maenas in response to S. aurata control metabolites in pH drop but not in control pH. Although this was not the main focus of our study, it could mean that pH drop renders predator odour more potent for its prey, for example through pH-dependent changes in odour molecular structure, receptor binding, or information processing (Munday et al., 2009;Roggatz et al., 2016Roggatz et al., , 2019Schirrmacher et al., 2020). Our results contrast with earlier reports of decreased predator avoidance (e.g., in snails; see Jellison et al., 2016;Froehlich and Lord, 2020). Overall, our data suggests that a short-term acute pH drop has the ability to impair decision-making and predator-prey interactions, and in turn alter fitness-relevant behaviours. Nevertheless, future studies should aim to disentangle the effects of pH drop from that of pCO 2 , which ideally also would test short, medium, and long-term responses.

Indirect Effects of pH Drop on Recipients' Behaviour Through Stress Chemical Communication
Despite the above-mentioned species-specific direct effects of pH drop on behaviour success ratios, we could show here that pH drop also has indirect effects on behaviour by altering chemical communication. We found that stress metabolites released by animals exposed to pH drop altered feeding and burrowing, and increased avoidance behaviours, but that the response was donor and recipient species-specific. Hermit crabs D. pugilator did not react to any stress metabolites. Significant responses in the other two recipient species, however, were directionally equivalent to those induced by pH drop itself since it took longer to find a food cue or to burrow the head in the presence of stress metabolites. Although behavioural stress response propagation through chemical communication is well documented for alarm substances and disturbance cues upon predation stress (Giacomini et al., 2015;Abreu et al., 2016;Mathuru, 2016), our results show here that chemical communication induced by abiotic stressors such as pH drop can trigger the release of waterborne chemicals that in turn trigger a behavioural stress response in its recipients. Stress metabolites conditioned by the potential predator S. aurata induced more avoidance behaviours and delayed feeding in C. maenas and burrowing in H. diversicolor, compared to S. aurata control metabolites. Similarly, the presence of conspecific stress metabolites significantly lowered H. diversicolor burrowing success ratio compared to conspecific control metabolites.
On the other hand, C. maenas did not react to conspecific stress metabolites. The different behavioural effect of conspecific stress metabolites between crabs and H. diversicolor has several possible explanations ranging from different defence strategies related to stress risk perception and evaluation (Hazlett, 1985;Bairos-Novak et al., 2017Goldman et al., 2020b) to different types, concentrations, or ratios of the released chemicals that serve as risk cues (Júnior et al., 2010;Morishita and Barreto, 2011). As we have shown here a pH drop to 7.6 is not a strong stressor in crabs (Fehsenfeld et al., 2011) which may explain that we failed to observe responses to water from pH dropconditioned conspecific green shore crabs. However, the presence of stress metabolites conditioned by S. aurata caused a marked drop in the success ratio of C. maenas at low (but not control) pH, evidenced by significant interaction terms.
Whilst we did not investigate their chemical structure, our experiment allows us to characterise novel properties of stress metabolites. Stress metabolites, like disturbance cues, potentially consist of regularly excreted metabolites such as urea and ammonia (Bairos-Novak et al., 2017;Shrivastava et al., 2019). However, the fact that behavioural responses of H. diversicolor and D. pugilator remained unchanged even when we used undiluted control metabolites might mean that stress metabolites are not just up-concentrated control metabolites (at least in the tested range of 10-fold control metabolites), or that the cocktail of chemicals contains much higher ratios of stress cues relative to control cues. We also observed that avoidance behaviours by both C. maenas and H. diversicolor became somewhat less pronounced with subsequent water uses, suggesting that the recipients may be able to discriminate degraded vs. fresh metabolites and react accordingly to the degree of threat they may indicate (Fuselier et al., 2009;Bairos-Novak et al., 2018). This shows that stress metabolites induced by pH drop are either not identical with metabolites released in normal pH, and/or are at least more than 10-fold higher concentrated than control metabolites.
The avoidance behaviours we recorded in this study are known indicators of stress. S. aurata stress metabolites elicited atypical freezing, eversion of the proboscis, mucus secretion, flipping and sideways-undulating behaviours in H. diversicolor which were previously described as indicators for a physiological stress response following exposure to copper sulphate (Burlinson and Lawrence, 2007), antiparasitic drugs (McBriarty et al., 2018), and trace metals (Mouneyrac et al., 2003). We also observed increased freeze and escape responses in crabs, which indicate a stress response in crustaceans (Katz and Rittschof, 1993;Perrot-Minnot et al., 2017;Tomsic et al., 2017). The finding that S. aurata stress metabolites induced such responses in potential prey species suggests that responses to disturbance cues are do not depend on the "audience, " similar as what has been observed in tadpoles (Bairos-Novak et al., 2020).
A directionally similar response to pH drop and to the stress metabolites it induces shows that pH drop can have indirect effects through stress propagation (Figure 1A), by which stressed animals directly experiencing environmental stressors negatively influence fitness-relevant behaviours of community members. The stress propagation toward naive conspecific individuals might constitute a "positive stress feedback loop" (Feugere et al., 2021). This was first observed by Hazlett (1985) after freshwater crayfish exposed to water conditioned with heat-stressed conspecifics displayed increased alertness, after returning the conditioned water to normal temperature. Similar indirect effects through chemical signalling from stressed donors toward unstressed, naive recipients exposed to the water conditioned by the directly stressed donors but not to the stressor itself, were also shown in response to different types of biotic stressors such as mock predator chase (Toa et al., 2004;Giacomini et al., 2015), handling (Barcellos et al., 2011), and acute fasting (Abreu et al., 2016), and to physical injury including irradiation (Mothersill et al., 2007) and predation (Frisch, 1938;Mathuru et al., 2012;Oliveira et al., 2013).

CONCLUSION AND PERSPECTIVES
In this study, we could show that short-term pH drop of a similar magnitude of that experienced within the intertidal zone, but also aligned to end-of-century predicted average values (Chavez et al., 2017;Landschützer et al., 2018), had negative consequences on fitness-relevant behaviours in harbour ragworm H. diversicolor and small hermit crab D. pugilator. Additionally, we confirm that pH drop events also impede the same behaviours in the same way indirectly via chemical communication, albeit these effects depended on donor and recipient species. Sea bream S. aurata and harbour ragworm H. diversicolor stressed by pH drop released stress metabolites which possibly differ from control metabolites and negatively affect fitness-relevant behaviours in metabolite recipients. Stress metabolites induced similar avoidance behaviours as those exhibited under physiological stress, meaning that a stress response was propagated from donor to recipient. These negative indirect behavioural effects of pH drop warrant further study, especially as our results were inconclusive with regards to the combined treatments of low pH and stress metabolites. Short-term pH drops thus involve behavioural additionally to metabolic trade-offs, which is also of interest for better predicting the response of natural and aquaculture systems under ocean acidification combined with tidal pH fluctuations. One potential limitation of our experimental design was that we could not achieve individual marking of animals. Despite circulating animals to recovery tanks to limit animal reuse, this may have led to a small degree of animal reuse wherein the full dataset may include animals tested at most twice over the course of the experimental period. Nevertheless, we show through random bootstrapping of 50% of the dataset, that the outcomes of the time-to-event analysis remain unchanged even if only half of the observations were analysed, which confirms the validity of our experimental data. Whether or not pseudoreplication is present in ecological studies has been subject to debate (Davies and Gray, 2015), but owing to our large number of observations we can show here that animal reuse had no effect on the model outcomes and effect size was not overestimated. Given the potential high ecological significance of indirect negative effects of pH fluctuations on population fitness, studies aiming to replicate our results should aim to further optimise our experimental design with better technical equipment as that available to us (Baker, 2016).

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

ETHICS STATEMENT
The animal study was reviewed and approved by the University of Hull Ethics Committee under the approvals U020 and FEC_2019_81.

AUTHOR CONTRIBUTIONS
KW designed the study. LF, SS, LA, JF, RN, KR, and KW performed the experiments. LF analysed the data and wrote the first manuscript draft with KW, JH, and HB-H. All authors contributed to the final manuscript.

FUNDING
JH and KW acknowledge funding by Natural Environment Research Council (NERC grant #NE/T001577/1), and KW, JH, and LF acknowledge funding by the University of Hull toward the MolStressH2O research cluster.