Heat stress does not induce wasting symptoms in the giant California sea cucumber (Apostichopus californicus) Declan Dawson Taylor1 ,2 ,* , Jonathan J. Farr2 ,3 ,* , Em G. Lim4 , Jenna L. Fleet2 ,5 , Sara J. Smith Wuitchik2 ,6 ,7 ,8 and Daniel M. Wuitchik2 ,6 1 Earth, Ocean and Atmospheric Sciences, University of British Columbia, Vancouver, British Columbia, Canada Bamfield Marine Sciences Center, Bamfield, British Columbia, Canada 3 Biological Sciences, University of Alberta, Edmonton, Alberta, Canada 4 Biological Sciences, Simon Fraser University, Vancouver, British Columbia, Canada 5 Biological Sciences, University of Winnipeg, Winnipeg, MB, Canada 6 Biological Sciences, Boston University, Boston, MA, United States of America 7 Informatics Group, Harvard University, Cambridge, MA, United States of America 8 Biology, Mount Royal University, Calgary, Alberta, Canada * These authors contributed equally to this work. 2 ABSTRACT Submitted 10 June 2022 Accepted 20 November 2022 Published 7 February 2023 Corresponding authors Sara J. Smith Wuitchik, sara.wuitchik@gmail.com, ssmith6@mtroyal.ca Daniel M. Wuitchik, wuitchik@bu.edu Academic editor Eric Ward Additional Information and Declarations can be found on page 11 Oceanic heatwaves have significant impacts on disease dynamics in marine ecosystems. Following an extreme heatwave in Nanoose Bay, British Columbia, Canada, a severe sea cucumber wasting event occurred that resulted in the mass mortality of Apostichopus californicus. Here, we sought to determine if heat stress in isolation could trigger wasting symptoms in A. californicus. We exposed sea cucumbers to (i) a simulated marine heatwave (22 ◦ C), (ii) an elevated temperature treatment (17 ◦ C), or (iii) control conditions (12 ◦ C). We measured the presence of skin lesions, mortality, posture maintenance, antipredator defences, spawning, and organ evisceration during the 79-hour thermal exposure, as well as 7-days post-exposure. Both the 22 ◦ C and 17 ◦ C treatments elicited stress responses where individuals exhibited a reduced ability to maintain posture and an increase in stress spawning. The 22 ◦ C heatwave was particularly stressful, as it was the only treatment where mortality was observed. However, none of the treatments induced wasting symptoms as observed in the Nanoose Bay event. This study provides evidence that sea cucumber wasting may not be triggered by heat stress in isolation, leaving the cause of the mass mortality event observed in Nanoose unknown. Subjects Ecology, Marine Biology, Zoology Keywords Wasting disease, Skin ulceration syndrome, Sea cucumber, Thermal stress, Echino- derm DOI 10.7717/peerj.14548 INTRODUCTION Copyright 2022 Dawson Taylor et al. Climate change is increasing the intensity, duration, range, and frequency of marine heatwaves across the globe with potentially catastrophic effects on organism fitness, ecosystems, and human economies (Frölicher, Fischer & Gruber, 2018; Allan et al., 2021). In marine ecosystems these extreme climatic events often cause immediate and mass mortality at all trophic levels from thermal stress, starvation, toxicity, and hypoxia (Cavole Distributed under Creative Commons CC-BY 4.0 OPEN ACCESS How to cite this article Dawson Taylor D, Farr JJ, Lim EG, Fleet JL, Smith Wuitchik SJ, Wuitchik DM. 2022. Heat stress does not induce wasting symptoms in the giant California sea cucumber (Apostichopus californicus). PeerJ 11:e14548 http://doi.org/10.7717/peerj.14548 et al., 2016; Di Lorenzo & Mantua, 2016; von Biela et al., 2019; Suryan et al., 2021). Marine heatwaves may exacerbate disease because thermal stress can compromise the immune response, and warmer temperatures can increase the virulence of pathogens (Harvell et al., 1999; Marcogliese, 2008; Branco et al., 2012; Matozzo et al., 2012; Oliver et al., 2017). This temperature-associated increase in virulence and infectivity has been associated with stimulating and intensifying sea star wasting disease to devastating effects (Harvell et al., 1999; Bates, Hilton & Harley, 2009; Eisenlord et al., 2016; Hewson et al., 2018; Aquino et al., 2021). A notorious example of the impacts of a wasting disease outbreak was the 2013-current sea star wasting epidemic, which has affected more than 20 sea star species in the Northeast Pacific Ocean over the last decade (Hewson et al., 2018). Sea star wasting disease encompasses a broad set of symptoms including twisted arms, lesions, deflation/loss of turgor, loss of arms, lack of grip strength in tube feet, and liquefaction (Bates, Hilton & Harley, 2009; Menge et al., 2016; Hewson et al., 2018). In most sea star species, the driving cause of wasting disease remains largely uncertain. It has often been assumed that wasting in echinoderms is driven by infectious agents (Hewson et al., 2014; Bucci et al., 2017; Miner et al., 2018; Hewson, Johnson & Tibbetts, 2020). Warm temperatures have also been linked to several mass mortality events (Bates, Hilton & Harley, 2009; Eisenlord et al., 2016; Harvell et al., 2019). Because of the population-level impacts on ecologically important species, understanding the pathogenic and environmental drivers of sea star wasting remains an area of active research (Aalto et al., 2020; Aquino et al., 2021; Hewson, 2021). Wasting is not limited to sea stars and is an emerging concern across closely related taxa. For example, sea urchins have faced a variety of disease-linked mass mortality events and epizootics (Feehan & Scheibling, 2014). Several of these bacterial and amoebic diseases have been linked to warm water anomalies, from climate events or storms (Sweet, 2020). Red sea urchins (Mesocentrotus franciscanus) and purple sea urchins (Strongylocentrotus purpuratus) have suffered from ‘‘bald sea urchin disease’’ and ‘‘sea urchin wasting disease’’ along the Pacific coast of British Columbia (B.C.). Pathology of these diseases includes lesions to the body wall and a shortening or loss of spines, and both have been associated with mass mortality events (Feehan & Scheibling, 2014; Sweet, 2020). While these diseases are believed to have bacterial origins, a single bacterial strain has not been identified as the cause (Sweet, 2020). Disease-mediated mass mortality events in urchins and sea stars have caused trophic cascades and dramatic ecosystem shifts, highlighting the importance of understanding these ecological phenomena (Schultz, Cloutier & Côté, 2016; Suryan et al., 2021; Traiger et al., 2022). Recent evidence has emerged that wasting may occur in sea cucumbers as well, promoting further concerns about the impacts of marine diseases in shallow, near-shore ecosystems. Since 2014, there have been reports of giant California sea cucumbers (Apostichopus californicus) with wasting symptoms similar to those of sea stars throughout their range in the northeast Pacific Ocean (Hewson, Johnson & Tibbetts, 2020). Sea cucumber wasting disease is poorly defined, with symptoms typically including non-focal lesions and fissures across the body wall, epidermal tissue sloughing, and rapid liquefaction (Hewson, Johnson & Tibbetts, 2020; Fig. 1C). Much like in sea stars, there are numerous potential causal agents that could elicit this suite of symptoms. Heat stress may be particularly relevant Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 2/17 A D N Heatwave (22oC) Average Temperature (ºC) 20 Warm (17oC) 15 B C Control (12oC) 10 5 January April July October January Figure 1 Temperature anomalies and sea cucumber wasting. (A) Map of western North America with (B) an inset of Vancouver Island showing the data collection locations of Bamfield Inlet (green), Nanoose Bay (purple). (C) A. californicus in Nanoose Bay exhibiting wasting symptoms. (D) Eight-day rolling average seasonal sea surface temperatures in Nanoose Bay from 2010 to 2020 (black line with grey SE ribbon) and 2021 (colour gradient line), with experimental temperature treatments highlighted. Full-size DOI: 10.7717/peerj.14548/fig-1 given that after a series of heatwaves from June 25 to July 07, 2021 (Environment and Climate Change Canada, 2021; Ocean Networks Canada) a wasting event with very high mortality occurred in the Strait of Georgia near Nanoose, B.C., from late August to October 2021 (Fig. 1D; E Lim, pers. comm., 2021). At its peak, up to 94% of observed A. californicus at a single site showed evidence of wasting. On average, across all affected sites (n = 7), 50% of observed A. californicus exhibited wasting symptoms. The six to ten week delay between peak near-surface sea temperatures and widespread mortality may be indicative of a potential link between thermal stress and wasting in A. californicus, perhaps mediated by the progression of disease symptoms and immune failure. Temperature is an especially compelling mechanism considering that healthy individuals were abundant below depths of 19 metres where the water was cooler (E Lim, pers. comm., 2021). The effect of temperature-related stressors on wasting in sea cucumbers has yet to be tested. Here, we assess whether acute heat stress rapidly induces wasting in A. californicus. We also evaluate behavioural stress response signatures through stiffening, spawning, and evisceration to explore how they are affected by extreme elevated temperature. Sea cucumbers are important benthic detritivores that play an important role in organic matter decomposition and sediment aeration (Purcell, Conand & Byrne, 2016). Evidence of a connection between temperature and mass wasting events in sea stars (Bates, Hilton & Harley, 2009; Eisenlord et al., 2016; Harvell et al., 2019) necessitates a proactive assessment of whether heat stress can induce similar symptoms in sea cucumbers. Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 3/17 MATERIALS & METHODS Collection and acclimation We collected Apostichopus californicus via SCUBA from a depth of 7–12 m in July 2021 from Scott’s Bay and Bamfield Inlet in Barkley Sound, British Columbia (48◦ 500 0200 N, 125◦ 080 4500 W). The sea cucumbers were initially collected and used in an unrelated short term tagging experiment comparing T-bar and coded wire tags at the end of July 2021, after which they were all tagged with coded wire tags. Tagging had no measurable impact on sea cucumber behaviour or physiology in the short or long term (Leedham et al., unpublished data; Cieciel, Pyper & Eckert, 2009). After this experiment sea cucumbers were then maintained in an open flow-through system at the Bamfield Marine Sciences Centre, with inflow from Bamfield Inlet for four months prior to the start of the thermal stress experiment. Experimental design To determine an average wet weight of each individual, sea cucumbers were weighed twice at the beginning of the experiment. Individuals (N = 56) were randomly divided into one of three temperature treatments: control (12 ◦ C, n = 19), warm (17 ◦ C, n = 19), or heatwave (22 ◦ C, n = 18). Our experimental temperatures reflect average summer temperatures and those recorded in the Strait of Georgia when wasting of A. californicus was observed (Fig. 1D). Specifically, the average temperature across July and August 2020 and 2021 was 17.5 ◦ C at 5 m depth and 11.9 ◦ C at 20 m depth (Fig. S1). In the summer 2021 heatwave, a maximum temperature at 5 m depth of 21.6 ◦ C was reached on August 04, 2021 (Pawlowicz, 2017; Xuereb et al., 2018; Chen et al., 2021; Ocean Networks Canada Data Archive, 2022). Our control temperature (12 ◦ C) was that of the Bamfield Inlet at the time of the temperature experiment and was the temperature the cucumbers were acclimated to. We separated the sea cucumbers into 27 isolated experimental containers (61×41×22.2 cm) with two individuals per container, aside from two smaller containers (33 ×45.7 ×11.4 cm) that housed one individual each. A plastic mesh divider was used to separate sea cucumbers within each container to allow for tracking individuals throughout the experiment. The containers were randomly placed into sea tables which acted as water baths. The target water bath temperatures were achieved by having a constant flow of cooler inlet sea water to act as the control (12 ◦ C), allowed to reach ambient room temperature for the warm treatment (17 ◦ C), or heated using two 800 W aquarium heaters for the heatwave treatment (22 ◦ C). The 22 ◦ C treatment temperature was gradually increased over a 24-hour period in the sea table water baths (day 1; Fig. S2). During the temperature treatment period, there was no seawater flow to the experimental containers in which sea cucumbers were housed for any of the temperature treatments. To prevent stagnation and maintain high oxygen levels, all containers were constantly aerated with airstones. The water remained at target temperatures for 79 h, after which they were lowered back to the control temperature of 12 ◦ C over 9-hours. Days five through 12 were a recovery period where all sea cucumber holding containers were maintained on the flow through system. Throughout the experiment, the sea cucumbers were closely monitored. To avoid water Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 4/17 fouling, water changes were completed in individual containers as needed to maintain nitrate and ammonium levels below 0.5 ppm, and fresh sea water was heated to the appropriate treatment temperature prior to water changes. Water quality was measured using API R saltwater ammonium and nitrate kits, and food was withheld during the thermal stress window. We monitored if sea cucumbers were defecating three times each day throughout the experiment, based on the presence of fecal pellets, as cessation of defecation is considered to be indicative of a loss or atrophy of digestive organs (Swan, 1961; Fankboner & Cameron, 1985). Phenotyping sea cucumber stress We assessed sea cucumbers for wasting symptoms throughout the experiment and counted skin lesions on days 5, 6, 7, and 12. Epidermal damage was classified as either minor lesions or major ulcers based on their size and appearance. We defined minor lesions as regions of epidermal damage at the spine tips without signs of considerable discoloration, where the dermis was not damaged across its full thickness (Fig. 2B). We classified major ulcers as epidermal ulceration exposing underlying white mutable connective tissues beneath (Fig. 2D). We assessed body stiffness using a 3-point ordinal scale to measure heat-stress induced deviations from normal behaviour. To do so, first we removed each individual from their bin and gently palpated them for 10 s to encourage them to stiffen. Then we placed them on an elevated platform to measure their ability to maintain their posture over 5 s. The stiffness testing platform was constructed by cutting a piece of plastic pipe into an eight cm long, eight cm in diameter, and four cm tall section, which we glued onto an eight cm tall block. We assigned a score of 0 if the organism failed to stiffen at all (full droop), a score of 1 if it did not remain stiff for the full 5 s when placed on the platform (partial droop), and a score of 2 if it maintained its posture for the entire 5 s (no droop). Stiffness was measured on days 1–5 of the experiment (as a baseline and throughout the treatment), on day 7 (48 h after experimental endpoint) and day 12 (7 days after experimental endpoint). Frequency of handling individuals occurred equally across treatments. As some sea cucumbers spawned, we recorded whether there was any evidence of spawning every 12 h. Spawning was assigned to each container rather than individual, since some individuals were co-housed (Nbin = 30). We also evaluated whether specimens had eviscerated based on the presence or absence of ejected internal organs during these routine checks. Statistical analyses All statistical analyses were conducted in R (v4.0.3, R Core Team, 2020). We tested the distribution of the number of lesions using the fitDist function (Rigby & Stasinopoulos, 2005) and determined that minor lesion count best fit a geometric distribution. We modelled the maximum number of minor lesions as a function of treatment, weight, evisceration (binary), and defecation status (binary). We included the sea table and bin as random effects. We then used backwards selection with the stepGAIC function from GAMLSS (Rigby & Stasinopoulos, 2005) to determine the most parsimonious combination of variables that best explained the maximum number of lesions. Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 5/17 Total Lesions / Indiv. 12 A B 8 4 0 Control Major Lesions / Indiv. 12 Warm Heat Wave C D 8 4 0 Control Warm Heat Wave Figure 2 Skin ulcers observed in Apostichopus californicus across treatments. (A) Minor ulcer count, defined as small lesions on the ends of spines (B). (C) Major ulcer count, defined as open wounds that expose white subdermal tissue (D). Full-size DOI: 10.7717/peerj.14548/fig-2 We examined the covariates that affected the likelihood of mortality with a logistic regression model. Sea cucumber mortality was a binary measure, with treatment, evisceration, defecation status, initial droop score, and initial weight as explanatory variables. Terms for the final model were selected via backwards model selection using stepGAIC (Rigby & Stasinopoulos, 2005) to produce the most parsimonious model of variation in sea cucumber mortality. To assess changes in stiffness, we constructed full ordinal regression models using the clmm function (Christensen, 2019) with temperature treatment, date, and the interaction between treatment and date as explanatory variables. We restricted stiffness measurements to those taken before the heat treatment began (day 1), during the treatment (days 2–4) Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 6/17 and immediately after the heat treatment (day 5). We included individual identity as a random effect to account for repeated measures on the same individuals over time. We also included bin and sea table as random effects to account for the paired (co-housing) and blocked (five bins per sea table) experimental design. We used backwards model selection to determine the most parsimonious models using dredge (Bartoń, 2020). To assess the drivers of evisceration, we constructed a logistic regression model with treatment, weight, and defecation status as explanatory variables, along with sea table as a random effect. We determined the most parsimonious model through backwards selection. To compare the incidence of stress spawning across temperature treatments, we conducted a Dunn’s Kruskal-Wallace (K-W) test using the dunnTest function from the FSA package (Ogle et al., 2022). All data and annotated code for the analyses described above are publicly available at https://github.com/declan-taylor/sea_cucumber_wasting (DOI: 10.5281/zenodo.7259325). RESULTS The temperature treatments varied slightly from the target temperatures during the experimental heatwaves. During the 79 h treatment, the mean temperature of the control treatment was 12.4 ◦ C and varied from 10.8 - 14.0 ◦ C; the mean of the warm treatment was 16.6 ◦ C and ranged from 14.8 ◦ C to 17.9 ◦ C; the mean of the heatwave treatment was 21.7 ◦ C and varied from 19.6 ◦ C to 23.3 ◦ C. Temperature treatments were significantly different from each other (K-W χ2 = 463.32, df = 2, p < 2.2e−16). Minor skin lesions occurred during the experiment and were not statistically different between treatments (ncontrol = 17, nwarm = 15, nheatwave = 10; Fig. 2). Major ulcers were only observed in the warm (n = 1) and heatwave treatments (n = 2). These major ulcers appeared to heal throughout the recovery period and were re-classified as minor lesions on day 12. The maximum number of minor lesions per individual was not significantly explained by treatment, weight, evisceration, or defecation status. Mortality was only observed in the 22 ◦ C heatwave treatment (Fig. 3; n = 5), while there were no mortalities observed in the control or warm treatments. Based on backwards model selection, mortality was driven by treatment and weight (Table S1). Body stiffness was lower in the warm and heatwave treatments compared to the control treatment (Fig. 3). Backwards-selected models indicated that temperature treatment and day affected stiffness (Table S2). Sea cucumbers were significantly less likely to have high stiffness scores relative to the control treatment in the warm (p = 1.99e−7) and heatwave (p = 2.44e−11) treatments. Structural stiffness values were significantly likely to be lower than day 1 on day 3 (p = 1.37e−5), day 4 (p = 2.50e−5) and day 5 (p = 8.66e−5), but not on day 2 (p = 0.0627; Table S2). Stress spawning was observed during the temperature treatment (n = 11 bins), with spawning occurring in the 17 ◦ C warm (n = 5 bins) and the 22 ◦ C heatwave (n = 4 bins) treatments. However, there was no significant difference in spawning between temperature treatments (K-W χ 2 = 1.94, df = 2, p = 0.379). Evisceration was observed across all treatments (ncontrol = 2, nwarm = 5, nheatwave = 5). Treatment temperature did not explain a significant amount of the variance in the Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 7/17 Ramp Temperature Stress Recovery Full Droop 15 Heat Wave 10 5 Number of Sea Cucumbers 0 Partial Droop 15 Warm 10 5 0 No Droop 15 Control 10 5 0 1 2 3 4 Experiment Day 5 6 7 Mortality Figure 3 Stiffness and mortality of Apostichopus californicus. Heat ramp (day 1), during temperature stress (days 2–4), and recovery (days 5–7) from the temperature treatment. Full-size DOI: 10.7717/peerj.14548/fig-3 occurrence of evisceration (Table S3), however, weight (p = 0.0383) and defecation status (p = 0.0163) were significant drivers increasing likelihood of evisceration (Table S3). Along with evisceration of internal organs, we also observed the expulsion of the respiratory tree in two individuals. Both individuals were in the 22 ◦ C heatwave treatment and mortality shortly followed the expulsion of the respiratory tree. DISCUSSION The objective of our study was to determine if heat stress can induce wasting symptoms in Apostichopus californicus. While we saw minor skin lesions at all treatment levels, and major ulcers in the warm and heatwave treatments, these were not characteristic of wasting symptoms observed previously (Hewson, Johnson & Tibbetts, 2020; Fig. 1). Neither the minor lesions nor major ulcers that we observed matched the wasting symptoms reported in A. californicus in Nanoose Bay, B.C. (E Lim, pers. comm., 2021), or the isolated wasting events reported along the Pacific coast (Hewson, Johnson & Tibbetts, 2020). We also did not see any sloughing of body tissues or liquefaction, as has been anecdotally reported in previous wasting events (Schroeder, 2017; Hewson, Johnson & Tibbetts, 2020). Despite Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 8/17 the major ulcers sharing some resemblance in colour, texture, and location to wasting disease symptoms, the sea cucumbers did not exhibit the full suite of symptoms that are typical from a wasting cucumber, and, within the treatment period, these lesions did not develop into fissures covering the dorsal body wall as seen in the wild (Fig. 1C). In Nanoose, wasting symptoms developed into a mass outbreak and mortality event over several months. However, unlike in reports of widespread wasting-driven mortality in wild A. californicus, the major ulcers in our specimens healed within the 7-day recovery period, and there was no evidence of these symptoms persisting or spreading to co-housed individuals. As such, we could not conclude that the sea cucumbers in our experiment were afflicted by the fatal wasting disease that has been previously reported (Schroeder, 2017; Hewson, Johnson & Tibbetts, 2020). Future experiments could address the multi-week delay between the heat dome and maximum die-off through prolonged exposure to elevated temperatures and utilise histopathology to further investigate the etiology of epidermal ulceration. While we are uncertain of the ultimate cause of the lesions that we observed, they were likely associated with handling while assessing their posture maintenance, which increased the frequency of potential abrasions on the epidermal surface tissue. Major ulcers may have begun as minor lesions that were then exacerbated by the high physiological stress caused by the 22 ◦ C treatment. White skin ulcerations, like the major ulcers we observed, are a recognized condition in other Holothuroidea, described as a Skin Ulceration Disease or Skin Ulceration Syndrome (SUS; Delroisse et al., 2020). SUS has been documented in commercially farmed Apostichopus japonicus and Holothuria scabra, and has been characterized by white ulcers on both sides of the body wall (Wang et al., 2007; Deng et al., 2009; Li et al., 2012; Zhang et al., 2018). Minor SUS symptoms in commercially farmed A. japonicus and the major ulcers in our A. californicus specimens are visually similar (Deng et al., 2009; Zhang et al., 2018). Unlike the SUS symptoms reported in A. japonicus, we did not see any indication of swelling or discolouration of the peristomes, and we did not see an initial abundance of ulcers around the mouth or cloaca (Becker et al., 2004; Wang et al., 2007; Delroisse et al., 2020). Extreme cases of SUS in farmed A. japonicus also bear resemblance to the wasting symptoms in wild A. californicus, raising further questions about the causes of skin ulceration in Holothuroidea (Delroisse et al., 2020). Despite not seeing evidence of wasting, we observed compelling evidence of stress symptoms in response to thermal stress. We observed five mortality events in the 22 ◦ C heatwave treatment. None of these mortalities were co-housed together, suggesting that it is unlikely a fatal disease spread between cohoused individuals. Therefore, we suggest that the heatwave treatment is close to the upper critical thermal tolerance of A. californicus, but the warm treatment temperature does not cause sufficient stress for mortality to occur. Our findings align with previous work on larval life stages, where Ren, Liu & Pearce (2018) found that at 22 ◦ C, larval A. californicus experienced reduced survival which was not observed at 16 ◦ C or 18 ◦ C. Beyond lethal effects, we observed reduced stiffening behaviour associated with both the 22 ◦ C heatwave and 17 ◦ C warm water treatment. A temperature-induced loss in stiffness may have implications for sea cucumber fitness under warming sea temperatures, Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 9/17 as a limited ability to stiffen may inhibit their ability to avoid predation or maintain posture while feeding and spawning. Thermal stress may have reduced stiffness by causing muscular fatigue and relaxation (Dowd & Somero, 2013) in the circular and longitudinalambulacral muscles (Gao & Yang, 2015). Stiffening is also caused by protein-mediated changes in mutable collagenous tissue within the dermis of sea cucumbers (Takehana et al., 2014). Therefore, heat stress may have reduced stiffening by denaturing or decreasing the production of tensilin, a stiffening protein, or increasing the production of the de-stiffening protein softenin (Yamada et al., 2010; Takehana et al., 2014; Tamori et al., 2016). Unlike stiffening behaviour, we did not observe any significant trends in spawning behaviour or evisceration. Of the 11 spawning events, nine occurred in elevatedtemperature treatments (warm or heatwave) and this trend, although statistically insignificant, was expected as stress spawning has been previously reported in other sea cucumber species (Battaglene et al., 2002; Rakaj et al., 2018; Schagerström et al., 2021). We were unable to associate spawning with individual phenotypes due to the co-housing experimental design. Evisceration appeared random across treatments, potentially because A. californicus in all treatments may have been reacting to handling stimulation and stress during the experiment (Ding et al., 2019). However, we found evidence that biological mechanisms, weight and defecation status, partially explain the non-treatment related variation in evisceration (Table S1). Defecation status in particular showed that A. californicus that were not defecating were more likely to eviscerate. This may have occurred because the energetic cost of eviscerating digestive organs would be lower for A. californicus that had already ceased using their organs, either because they were preparing to eviscerate (Swan, 1961) or undergoing viscera atrophy (Fankboner & Cameron, 1985). As such, when stressed by handling, A. californicus that had already begun seasonal reductions in digestive function may have been more likely to eviscerate from stress and overstimulation. Unlike digestive tract evisceration, we do not believe that the expulsion of the respiratory tree in the 22 ◦ C heatwave treatment was linked to seasonal senescence. For both individuals, respiratory evisceration was followed by mortality, suggesting that this is an indication of extreme physiological distress from thermal exposure. Since we observed skin lesions under thermal stress that were likely unassociated with wasting, other factors are likely causing recent wasting outbreaks in A. californicus. Like wasting, severe cases of SUS are highly transmissible and result in mortality (Delroisse et al., 2020); bacterial and viral pathogenic causal agents have been previously linked to SUS and wasting-like symptoms, both in other sea cucumbers (Deng et al., 2008; Deng et al., 2009; Liu et al., 2010; Zhang et al., 2018; Delroisse et al., 2020) and sea stars (Hewson et al., 2014; Hewson et al., 2018; Work et al., 2021). A study examining a single wasting A. californicus specimen found a high viral load but was unable to identify a specific pathogen causing wasting symptoms (Hewson, Johnson & Tibbetts, 2020). Avenues for future research on wasting diseases in A. californicus could address the potential for shared etiology with SUS, given the symptomatic similarities, as well as the possibility of pathogenic origins. Such investigations would benefit from histopathological examinations of potentially diseased tissue to reveal more about the symptom’s pathogenesis. Infectious factors (viral, bacterial) and non-infectious factors (chemical pollution, hypoxia, eutrophication) should both be Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 10/17 investigated because widespread environmental degradation and anthropogenic climate change are shifting pathogenic dynamics globally (Marcogliese, 2008; Allan et al., 2021). These investigations would add valuable insight to the field of wasting diseases, given the scarcity of published information on these events in echinoderms. In this study, we exposed A. californicus to extreme thermal stress as measured by mortality, degraded stiffening behaviour, and the development of skin lesions. Despite this, we found no evidence that wasting is induced by temperature stress alone. Therefore, the August 2021 mass wasting event in Nanoose, British Columbia, was likely not triggered solely by the anomalous heatwave. Determining the factors that cause and exacerbate wasting in A. californicus is essential for predicting and managing mass mortality events. Sea cucumbers are ecologically important benthic detritivores, which break down organic matter, recycle nutrients, and maintain sediment health (Wheeling, Verde & Nestler, 2007; Purcell, Conand & Byrne, 2016). Efforts to protect, manage, and sustainably harvest A. californicus in the face of global environmental change will require a comprehensive understanding of their stress responses, disease dynamics, and the novel threat of sea cucumber wasting. ACKNOWLEDGEMENTS The species collections and experiments took place on the traditional territories of the Huu-ay-aht First Nations, a Nuu-chah-nulth Nation and signatory to the Maa-nulth First Nations Final Agreement, and we are grateful for the opportunity to conduct research in protected and sacred areas. We would like to thank Chloe Curry, Arya Horon, and Juliane Jones for providing lab assistance; Payton Arthur, Mike Chung, Gabrielle Languedoc, Sammie Foley, Juliane Jones, and Carter Burtlake for feedback on early versions of this manuscript. ADDITIONAL INFORMATION AND DECLARATIONS Funding Funding was provided by Bamfield Marine Sciences Centre to the Fall Program and by Mount Royal University to Sara J Smith Wuitchik. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Grant Disclosures The following grant information was disclosed by the authors: Bamfield Marine Sciences Centre to the Fall Program. Mount Royal University. Competing Interests The authors declare there are no competing interests. Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 11/17 Author Contributions • Declan Dawson Taylor conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft. • Jonathan J. Farr conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft. • Em G. Lim conceived and designed the experiments, prepared figures and/or tables, and approved the final draft. • Jenna L. Fleet conceived and designed the experiments, prepared figures and/or tables, and approved the final draft. • Sara J. Smith Wuitchik conceived and designed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft. • Daniel M. Wuitchik conceived and designed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft. Data Availability The following information was supplied regarding data availability: The data and code is available at GitHub: https://github.com/declan-taylor/sea_ cucumber_wasting. Supplemental Information Supplemental information for this article can be found online at http://dx.doi.org/10.7717/ peerj.14548#supplemental-information. REFERENCES Aalto EA, Lafferty KD, Sokolow SH, Grewelle RE, Ben-Horin T, Boch CA, Raimondi PT, Bograd SJ, Hazen EL, Jacox MG, Micheli F, De Leo GA. 2020. Models with environmental drivers offer a plausible mechanism for the rapid spread of infectious disease outbreaks in marine organisms. Scientific Reports 10:5975 DOI 10.1038/s41598-020-62118-4. Allan RP, Arias PA, Berger S, Canadell JG, Cassou C, Chen D, Cherchi A, Connors SL, Coppola E, Cruz FA, Diongue-Niang A, Doblas-Reyes FJ, Douville H, Driouech F, Edwards TL, Engelbrecht F, Eyring V, Fischer E, Flato GM, Forster P, Fox-Kemper B, Fuglestvedt JS, Fyfe JC, Gillett NP, Gomis MI, Gulev SK, Gutiérrez JM, Hamdi R, Harold J, Hauser M, Hawkins E, Hewitt HT, Johansen TG, Jones C, Jones RG, Kaufman DS, Klimont Z, Kopp RE, Koven C, Krinner G, Lee J-Y, Lorenzoni I, Marotzke J, Masson-Delmotte V, Maycock TK, Meinshausen M, Monteiro PMS, Morelli A, Naik V, Notz D, Otto F, Palmer MD, Pinto I, Pirani A, Plattner G-K, Raghavan K, Ranasinghe R, Rogelj J, Rojas M, Ruane AC, Sallée J-B, Samset BH, Seneviratne SI, Sillmann J, Sörensson AA, Stephenson TS, Storelvmo T, Szopa S, Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 12/17 Thorne PW, Trewin B, Vautard R, Vera C, Yassaa N, Zaehle S, Zhai P, Zhang X, Zickfeld K. 2021. IPCC: climate Change 2021: the physical science basis. Cambridge: Cambridge University PressAvailable at https://www.ipcc.ch/report/ar6/wg1/#SPM . Aquino CA, Besemer RM, De Rito CM, Kocian J, Porter IR, Raimondi PT, Rede JE, Schiebelhut LM, Sparks JP, Wares JP, Hewson I. 2021. Evidence that microorganisms at the animal-water interface drive sea star wasting disease. Frontiers in Microbiology 11:610009 DOI 10.3389/fmicb.2020.610009. Bartoń KA. 2020. MuMIn: multi-model inference. R package version 1.47.1. Bates AE, Hilton BJ, Harley CDG. 2009. Effects of temperature, season and locality on wasting disease in the keystone predatory sea star Pisaster Ochraceus. Diseases of Aquatic Organisms 86:245–251 DOI 10.3354/dao02125. Battaglene SC, Seymour JE, Ramofafia C, Lane I. 2002. Spawning induction of three tropical sea cucumbers, Holothuria Scabra, H. Fuscogilva and Actinopyga Mauritiana. Aquaculture 207:29–47 DOI 10.1016/S0044-8486(01)00725-6. Becker P, Gillan D, Lanterbecq D, Jangoux M, Rasolofonirina R, Rakotovao J, Eeckhaut I. 2004. The skin ulceration disease in cultivated juveniles of Holothuria scabra (Holothuroidea, Echinodermata). Aquaculture 242:13–30 DOI 10.1016/j.aquaculture.2003.11.018. Branco PC, Pressinotti LN, Borges JCS, Iunes RS, Kfoury JR, da Silva MO, Gonzalez M, dos Santos MF, Peck LS, Cooper EL, da Silva JRMC. 2012. Cellular biomarkers to elucidate global warming effects on Antarctic sea urchin Sterechinus neumayeri. Polar Biology 35:221–229 DOI 10.1007/s00300-011-1063-5. Bucci C, Francoeur M, McGreal J, Smolowitz R, Zazueta-Novoa V, Wessel GM, GomezChiarri M. 2017. Sea star wasting disease in asterias forbesi along the Atlantic Coast of North America. PLOS ONE 12:e0188523 DOI 10.1371/journal.pone.0188523. Cavole L, Demko A, Diner R, Giddings A, Koester I, Pagniello C, Paulsen M-L, Ramirez-Valdez A, Schwenck S, Yen N, Zill M, Franks P. 2016. Biological impacts of the 2013–2015 warm-water anomaly in the northeast pacific: winners, losers, and the future. Oceanography 29:273–285 DOI 10.5670/oceanog.2016.32. Chen Z, Shi J, Liu Q, Chen H, Li C. 2021. A Persistent and intense marine heatwave in the northeast pacific during 2019–2020. Geophysical Research Letters 48:e2021GL093239 DOI 10.1029/2021GL093239. Christensen RHB. 2019. Ordinal–regression models for ordinal data. R package version 2019.12-10. Available at https://CRAN.R-project.org/package=ordinal. Cieciel K, Pyper B, Eckert G. 2009. Tag retention and effects of tagging on movement of the giant red sea cucumber Parastichopus californicus. North American Journal of Fisheries Management 29:288–294 DOI 10.1577/M07-194.1. Delroisse J, Van Wayneberghe K, Flammang P, Gillan D, Gerbaux P, Opina N, Todinanahary GGB, Eeckhaut I. 2020. Epidemiology of a SKin Ulceration Disease (SKUD) in the sea cucumber Holothuria scabra with a review on the SKUDs in Holothuroidea (Echinodermata). Scientific Reports 10:22150 DOI 10.1038/s41598-020-78876-0. Deng H, He C, Zhou Z, Liu C, Tan K, Wang N, Jiang B, Gao X, Liu W. 2009. Isolation and pathogenicity of pathogens from skin ulceration disease and viscera ejection Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 13/17 syndrome of the sea cucumber Apostichopus japonicus. Aquaculture 287:18–27 DOI 10.1016/j.aquaculture.2008.10.015. Deng H, Zhou Z, Wang N, Liu C. 2008. The syndrome of sea cucumber (Apostichopus japonicus) infected by virus and bacteria. Virologica Sinica 23:63–67 DOI 10.1007/s12250-008-2863-9. Ding K, Zhang L, Sun L, Lin C, Feng Q, Zhang S, Yang H, Brinkman R, Lin G, Huang Z. 2019. Transcriptome analysis provides insights into the molecular mechanisms responsible for evisceration behavior in the sea cucumber Apostichopus japonicus. Comparative Biochemistry and Physiology Part D: Genomics and Proteomics 30:143–157 DOI 10.1016/j.cbd.2019.02.008. Dowd WW, Somero GN. 2013. Behavior and survival of Mytilus congeners following episodes of elevated body temperature in air and seawater. Journal of Experimental Biology 216:502–514 DOI 10.1242/jeb.076620. Environment and Climate Change Canada. 2021. Warns western Canadians about dangerous record-high temperatures. Available at https://www.canada.ca/en/ environment-climate-change/news/2021/06/environment-and-climate-change-canadawarns-western-canadians-about-dangerous-record-high-temperatures.html. Eisenlord ME, Groner ML, Yoshioka RM, Elliott J, Maynard J, Fradkin S, Turner M, Pyne K, Rivlin N, Van Hooidonk R, Harvell CD. 2016. Ochre star mortality during the 2014 wasting disease epizootic: role of population size structure and temperature. Philosophical Transactions of the Royal Society B: Biological Sciences 371:20150212 DOI 10.1098/rstb.2015.0212. Fankboner PV, Cameron JL. 1985. Seasonal atrophy of the visceral organs in a sea cucumber. Canadian Journal of Zoology 63:2888–2892 DOI 10.1139/z85-432. Feehan CJ, Scheibling RE. 2014. Effects of sea urchin disease on coastal marine ecosystems. Marine Biology 161:1467–1485 DOI 10.1007/s00227-014-2452-4. Frölicher TL, Fischer EM, Gruber N. 2018. Marine heatwaves under global warming. Nature 560:360–364 DOI 10.1038/s41586-018-0383-9. Gao F, Yang H. 2015. Chapter 4—anatomy. Developments in Aquaculture and Fisheries Science 39:53–76 DOI 10.1016/B978-0-12-799953-1.00004-0. Harvell CD, Kim K, Burkholder JM, Colwell RR, Epstein PR, Grimes DJ, Hofmann EE, Lipp EK, Osterhaus ADME, Overstreet RM, Porter JW, Smith GW, Vasta GR. 1999. Emerging marine diseases—climate links and anthropogenic factors. Science 285:1505–1510 DOI 10.1126/science.285.5433.1505. Harvell CD, Montecino-Latorre D, Caldwell JM, Burt JM, Bosley K, Keller A, Heron SF, Salomon AK, Lee L, Pontier O, Pattengill-Semmens C, Gaydos JK. 2019. Disease epidemic and a marine heatwave are associated with the continental-scale collapse of a pivotal predator (Pycnopodia helianthoides). Science Advances 5:eaau7042 DOI 10.1126/sciadv.aau7042. Hewson I. 2021. Microbial respiration in the asteroid diffusive boundary layer influenced sea star wasting disease during the 2013-2014 northeast Pacific Ocean mass mortality event. Marine Ecology Progress Series 668:231–237 DOI 10.3354/meps13735. Hewson I, Bistolas KSI, Quijano Cardé EM, Button JB, Foster PJ, Flanzenbaum JM, Kocian J, Lewis CK. 2018. Investigating the complex association between viral Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 14/17 ecology, environment, and northeast Pacific sea star wasting. Frontiers in Marine Science 5:77 DOI 10.3389/fmars.2018.00077. Hewson I, Button JB, Gudenkauf BM, Miner B, Newton AL, Gaydos JK, Wynne J, Groves CL, Hendler G, Murray M, Fradkin S, Breitbart M, Fahsbender E, Lafferty KD, Kilpatrick AM, Miner CM, Raimondi P, Lahner L, Friedman CS, Daniels S, Haulena M, Marliave J, Burge CA, Eisenlord ME, Harvell CD. 2014. Densovirus associated with sea-star wasting disease and mass mortality. Proceedings of the National Academy of Sciences of the United States of America 111:17278–17283 DOI 10.1073/pnas.1416625111. Hewson I, Johnson MR, Tibbetts IR. 2020. An unconventional flavivirus and other RNA viruses in the sea cucumber (Holothuroidea; Echinodermata) virome. Viruses 12:1057 DOI 10.3390/v12091057. Li C, Feng W, Qiu L, Xia C, Su X, Jin C, Zhou T, Zeng Y, Li T. 2012. Characterization of skin ulceration syndrome associated microRNAs in sea cucumber Apostichopus japonicus by deep sequencing. Fish & Shellfish Immunology 33:436–441 DOI 10.1016/j.fsi.2012.04.013. Liu H, Zheng F, Sun X, Hong X, Dong S, Wang B, Tang X, Wang Y. 2010. Identification of the pathogens associated with skin ulceration and peristome tumescence in cultured sea cucumbers Apostichopus japonicus (Selenka). Journal of Invertebrate Pathology 105:236–242 DOI 10.1016/j.jip.2010.05.016. Lorenzo EDi, Mantua N. 2016. Multi-year persistence of the 2014/15 North Pacific marine heatwave. Nature Climate Change 6:1042–1047 DOI 10.1038/nclimate3082. Marcogliese D. 2008. The impact of climate change on the parasites and infectious diseases of aquatic animals. Revue scientifique et technique 27:267–484 DOI 10.1016/j.jip.2010.05.016. Matozzo V, Chinellato A, Munari M, Finos L, Bressan M, Marin MG. 2012. First evidence of immunomodulation in bivalves under seawater acidification and increased temperature. PLOS ONE 7:e33820 DOI 10.1371/journal.pone.0033820. Menge BA, Cerny-Chipman EB, Johnson A, Sullivan J, Gravem S, Chan F. 2016. Sea star wasting disease in the keystone predator pisaster ochraceus in oregon: insights into differential population impacts, recovery, predation rate, and temperature effects from long-term research. PLOS ONE 11:e0153994 DOI 10.1371/journal.pone.0153994. Miner CM, Burnaford JL, Ambrose RF, Antrim L, Bohlmann H, Blanchette CA, Engle JM, Fradkin SC, Gaddam R, Harley CDG, Miner BG, Murray SN, Smith JR, Whitaker SG, Raimondi PT. 2018. Large-scale impacts of sea star wasting disease (SSWD) on intertidal sea stars and implications for recovery. PLOS ONE 13:e0192870 DOI 10.1371/journal.pone.0192870. Ocean Networks Canada Data Archive. 2022. Temperature data from 14 Jun 20191 Jan 2022. Victoria: University of Victoria, Canada. Available at https://data. oceannetworks.ca/DataSearch?locationCode=BSM.J1&deviceCategoryCode=CTD (accessed on 14 February 2022). Ogle DH, Doll JC, Wheeler P, Dinno A. 2022. FSA: fisheries stock analysis. R package version 0.9.3. GitHub. Available at https://github.com/fishR-Core-Team/FSA. Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 15/17 Oliver ECJ, Benthuysen JA, Bindoff NL, Hobday AJ, Holbrook NJ, Mundy CN, PerkinsKirkpatrick SE. 2017. The unprecedented 2015/16 Tasman Sea marine heatwave. Nature Communications 8:16101 DOI 10.1038/ncomms16101. Pawlowicz R. 2017. Seasonal cycles, hypoxia, and renewal in a coastal Fjord (Barkley Sound, British Columbia). Atmosphere-Ocean 55:264–283 DOI 10.1080/07055900.2017.1374240. Purcell SW, Conand C, Byrne SUM. 2016. Ecological roles of exploited sea cucumbers. In: Oceanography and marine biology. Boca Raton, Florida, USA: CRC Press. R Core Team. 2020. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing. Available at https://www.r-project.org . Rakaj A, Fianchini A, Boncagni P, Lovatelli A, Scardi M, Cataudella S. 2018. Spawning and rearing of Holothuria tubulosa: a new candidate for aquaculture in the Mediterranean region. Aquaculture Research 49:557–568 DOI 10.1111/are.13487. Ren Y, Liu W, Pearce CM. 2018. Effects of stocking density, ration and temperature on growth, survival and metamorphosis of auricularia larvae of the california sea cucumber, Parastichopus Californicus. Aquaculture Research 49:517–525 DOI 10.1111/are.13482. Rigby RA, Stasinopoulos DM. 2005. Generalized additive models for location, scale and shape. Journal of the Royal Statistical Society: Series C (Applied Statistics) 54:507–554 DOI 10.1111/j.1467-9876.2005.00510.x. Schagerström E, Christophersen G, Sunde J, Bakke S, Matusse NR, Dupont S, Sundell KS. 2021. Controlled spawning and rearing of the sea cucumber, Parastichopus Tremulus. Journal of the World Aquaculture Society 53:224–240 DOI 10.1111/jwas.12816. Schroeder L. 2017. Wasting-like lesions occurring on California sea cucumbers. The Dredgings 57:3. Schultz JA, Cloutier RN, Côté IM. 2016. Evidence for a trophic cascade on rocky reefs following sea star mass mortality in British Columbia. PeerJ 4:e1980 DOI 10.7717/peerj.1980. Suryan RM, Arimitsu ML, Coletti HA, Hopcroft RR, Lindeberg MR, Barbeaux SJ, Batten SD, Burt WJ, Bishop MA, Bodkin JL, Brenner R, Campbell RW, Cushing DA, Danielson SL, Dorn MW, Drummond B, Esler D, Gelatt T, Hanselman DH, Hatch SA, Haught S, Holderied K, Iken K, Irons DB, Kettle AB, Kimmel DG, Konar B, Kuletz KJ, Laurel BJ, Maniscalco JM, Matkin C, McKinstry CAE, Monson DH, Moran JR, Olsen D, Palsson WA, Pegau WS, Piatt JF, Rogers LA, Rojek NA, Schaefer A, Spies IB, Straley JM, Strom SL, Sweeney KL, Szymkowiak M, Weitzman BP, Yasumiishi EM, Zador SG. 2021. Ecosystem response persists after a prolonged marine heatwave. Scientific Reports 11:6235 DOI 10.1038/s41598-021-83818-5. Swan EF. 1961. Seasonal evisceration in the sea cucumber, Parastichopus Californicus (Stimpson). Science 133:1078–1079 DOI 10.1126/science.133.3458.1078. Sweet M. 2020. Sea urchin diseases: effects from individuals to ecosystems. Developments in Aquaculture and Fisheries Science 43:219–226 DOI 10.1016/B978-0-12-819570-3.00012-3. Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 16/17 Takehana Y, Yamada A, Tamori M, Motokawa T. 2014. Softenin, a novel protein that softens the connective tissue of sea cucumbers through inhibiting interaction between collagen fibrils. PLOS ONE 9:e85644 DOI 10.1371/journal.pone.0085644. Tamori M, Ishida K, Matsuura E, Ogasawara K, Hanasaka T, Takehana Y, Motokawa T, Osawa T. 2016. Ultrastructural changes associated with reversible stiffening in catch connective tissue of sea cucumbers. PLOS ONE 11:e0155673 DOI 10.1371/journal.pone.0155673. Traiger SB, Bodkin JL, Coletti HA, Ballachey B, Dean T, Esler D, Iken K, Konar B, Lindeberg MR, Monson D, Robinson B, Suryan RM, Weitzman BP. 2022. Evidence of increased mussel abundance related to the Pacific marine heatwave and sea star wasting. Marine Ecology 43:e12715 DOI 10.1111/maec.12715. von Biela V, Arimitsu M, Piatt J, Heflin B, Trowbridge SK, Schoen SK, Trowbridge J, Clawson C. 2019. Extreme reduction in nutritional value of a key forage fish during the Pacific marine heatwave of 2014–2016. Marine Ecology Progress Series 613:171–182 DOI 10.3354/meps12891. Wang P, Chang Y, Yu J, Li C, Xu G. 2007. Acute peristome edema disease in juvenile and adult sea cucumbers Apostichopus japonicus (Selenka) reared in North China. Journal of Invertebrate Pathology 96:11–17 DOI 10.1016/j.jip.2007.03.001. Wheeling RJ, Verde EA, Nestler JR. 2007. Diel cycles of activity, metabolism, and ammonium concentration in tropical holothurians. Marine Biology 152:297–305 DOI 10.1007/s00227-007-0683-3. Work TM, Weatherby TM, De Rito CM, Besemer RM, Hewson I. 2021. Sea star wasting disease pathology in Pisaster ochraceus shows a basal-to-surface process affecting color phenotypes differently. Diseases of Aquatic Organisms 145:21–33 DOI 10.3354/dao03598. Xuereb A, Benestan L, Normandeau É, Daigle RM, Curtis JMR, Bernatchez L, Fortin M-J. 2018. Asymmetric oceanographic processes mediate connectivity and population genetic structure, as revealed by RADseq, in a highly dispersive marine invertebrate (Parastichopus Californicus). Molecular Ecology 27:2347–2364 DOI 10.1111/mec.14589. Yamada A, Tamori M, Iketani T, Oiwa K, Motokawa T. 2010. A novel stiffening factor inducing the stiffest state of holothurian catch connective tissue. Journal of Experimental Biology 213:3416–3422 DOI 10.1242/jeb.044149. Zhang Z, Xing R, Lv Z, Shao Y, Zhang W, Zhao X, Li C. 2018. Analysis of gut microbiota revealed Lactococcus garviaeae could be an indicative of skin ulceration syndrome in farmed sea cucumber Apostichopus japonicus. Fish & Shellfish Immunology 80:148–154 DOI 10.1016/j.fsi.2018.06.001. Dawson Taylor et al. (2022), PeerJ, DOI 10.7717/peerj.14548 17/17