Thermal preference, tolerance and temperature-dependent respiration in the California sea hare Aplysia californica ()
1. INTRODUCTION
In the case of ectothermic organisms the growth rate is temperature-dependent. Several laboratory studies of A. californica have documented these rates at different rearing temperatures. It has been a standard hatchery practice to maintain a post-metamorphosis rearing temperature of 15˚C, with the objective to prevent the onset of sexual maturity before 9 months of age [1-4].
Recently, environmental factors involving influence on growth of A. californica have been studied in laboratoryreared under controlled conditions [5]. The decision to analyze Aplysia’s behavior in this particular ambiance has its logic in the similar conduct Aplysia expresses in the field or in a regulated context. Nonetheless, there are differences caused by the greater environmental factors in the field compared to laboratory environments [6].
Aplysia californica is found between Oregon and the Gulf of California, and Aplysia vaccaria is extended from California to the Gulf of California [7]. [8] found only one specimen of Aplysia californica in Las Conchas Ecuador. [3,4,9] varied on the following factors such as temperature, food or stocking density in A. californica, maintained in laboratory conditions, while keeping the others constant; these studies demonstrated that each factor had a strong effect on growth rate and size at sexual maturity. One factor that does not appear to influence growth and maturation, however, is the time of year the life cycle begins [2,10,11].
[5] correlated that rearing temperatures would influence the growth rate and longevity at which A. californica becomes sexually mature in the laboratory. They observed that A. californica lives more days in low temperatures than higher, but the size is larger in high temperatures. Probably there are multiple effects in its metabolism, unfortunately there is not enough information in order to know them or to determine which would be its best performance at a given temperature.
Behavioral thermoregulation can be adaptive in two complementary ways: 1) it can help an organism avoid extreme heat or cold that could be damaging or lethal [12]; 2) it can increase the time an animal spending at optimal physiological temperatures [13]. The thermal preference of a species corresponds closely with a temperature that maximizes growth, and other physiological processes [14,15]. The temperature intervals in which the organism congregates or spends time are defined operationally as preferred temperature or thermal preference [16]. In order to determine these two processes two methods are utilized, acute and gravitational by bipartite definition [17]. As for the gravitational method it takes 24 hours and the preferendum to determine organisms gravitates freely towards their preferred temperature [18, 19].
[20] mentions that critical thermal limits are typically determined experimentally within a controlled laboratory environment and are often based on dynamic non-lethal or static lethal methods. The Critical Thermal Method involves a dynamic changing temperature at a constant rate until a Critical Thermal End Point (predefined sublethal endpoint), is reached. In order to measure thermal tolerance in marine environment, where organisms are often exposed to relatively high temperatures, the CTMax previously mentioned is likely to be one of the principal methods in terms of ecological significance [21].
Thermal stress in aquatic organisms is likely to increase with the effect of global climate change. CTE is a behavioral stress response and is defined as the arithmetic mean of collecting thermal points at which mobile activity becomes disorganized to the point that organisms lose their ability to escape conditions that will lead to death [22]. The analogous behavioral stress response of aquatic macroinvertebrates varies amongst taxa but generally includes inability to remain attached to a substrate or increased movement, followed by immobility with a lack of response when stimulated [23].
In the intertidal zone, organisms can spend part of their lives immersed in stable and benign seawater and the remaining time in extreme conditions during emersion. Therefore, the relationship between temperature and metabolism of animals living in the intertidal zone can provide a suitable model to clarify their physiological adaptation to temperature fluctuation [24]. Among the physiological parameters that correlate with environmental changes such as temperature, the measurement of the metabolic rate allows us to determine the energetic cost that combinations have on an organism [25-27]. According to [28,29] the oxygen consumption rate (OCR) is intimately associated with metabolic work and the energy flow that an organism can use for homeostatic control mechanisms.
The aim of this study is to know in Aplysia californica, the preferred temperature in which it can be correlated to the optimal growth, the effect of critical temperature in order to know its thermal tolerance, and in which temperature it has its best metabolic response.
2. MATERIALS AND METHODS
Sea Hare and Experimental Design
Adult Aplysia californica were collected in nets in the intertidal zone in the area of Punta Morro Baja California, México (32.25˚ Lat N 116.983 Long W). The organisms were kept for thirty days in reservoirs of 2000 L, 35‰ salinity and 16˚C ± 1˚C as the same conditions measure in the moment of capture in a flow-through water system. The photoperiod was maintained at 12 - 12 h light-dark with a 30 minute transition periods in between.
We used 90 organisms acclimated to 16˚C, 19˚C and 22˚C ± 1˚C experimental condition. The Aplysia californica was kept for 21 days in acclimated time. Afterwards, a group of 6 organisms was introduced to the horizontal thermal gradient described by [30]. A Neslab thermoregulator (HX model 150) was connected at one end to cool the water gradient. A 1000 watts heater was placed at the opposite end connected to a temperature controller (Moeller CI-K3-125-M), generating a steady gradient of 10˚C to 30˚C, (y = 6.20 + 1.60x, r2 = 0.98 where x = gradient segments, and each segment = temperature gradient). Twenty aeration stones were placed along the gradient to maintain a dissolved oxygen concentration of 5.8 - 7.2 mg·L−1 and thus eliminate the effect of stratification of the water column. The depth of the water column in the gradient was 9 cm and to maintain good water quality turnover rate was 10.0 - 12.0 L·h−1. Every organism was tagged with a nylon tread, a circular plastic tag allowed us to identify the experimental organism behavior; now we could determine the preferred temperature gravitational method as applied and described by [18].
Six tagged Aplysia were introduced to the gradient in the segment where the water temperature was the same as that of the experimental acclimation temperatures. They were exposed at day to a light intensity of 0.55 × 1015 quanta·sec−1·cm−2 and at night to a light intensity of 0.09 × 1015 quanta·sec−1·cm−2. For each temperature there were four trials performed (N = 60). Nudibranches were not fed 24 h before the trial period to avoid interference from the diet [31,32].
The velocity displacement of the sea hare in the gradient was determined by observing the tagged organism in each segment of the trough every hour during the 24 h period, expressed in cm·h−1. A Mann-Whitney test was used to determine the differences in preferred temperature and the displacement velocity of sea hare in the light and dark cycle [33]. Temperature preference data was processed with Exploratory Data Analysis [34] and plotted as parallel boxes.
To determine the Critical Thermal Maxima (CTMax), 45 sea hare were acclimated in three 2000 L tanks at temperatures of 16˚C, 19˚C and 22˚C ± 1˚C, and a salinity of 35‰, for 21 days to allow acclimation to experimental temperatures. Twenty four hours before trials initiated the sea hare were not fed. The organisms were placed on acrylic plates (30 × 25 × 0.4 cm) and left undisturbed overnight; the acrylic plates were suspended 12 h in 40 L glass aquarium with aeration. Each aquarium was thermostatically controlled with a 1000 W heater and an air stone (15 × 5 × 5 cm). The temperature was raised from the acclimation temperatures of 16˚C, 19˚C and 22˚C ± 1˚C at a rate of 1˚C every 15 minutes [30,35,36].
Determinations of CTMax of the sea hare responses were made between 10:00 and 13:00 h. visual monitoring was performed and recorded the time and temperature at which they showed a curled mantle and started to lose attachment, the organisms immediately were transferred into acclimation tanks for recovery. The organisms were used only once. The data from sea hare that did not recover after returning to acclimation temperatures was discarded. The Kruskall-Wallis test was used to determine differences between groups [33]. The temperature, at which 50% remain attached, was calculated using linear regression, previous transformation to arcsine SigmaStat, version 3.1 (Systat Software Inc., San Jose, CA) on all points below 100%.
Oxygen consumption was measured in a semi open system, 15 organisms from each acclimation temperature were placed individually in a 1.0 L Erlenmeyer flask. Sea hare were placed in the respirometric system 24 h before the beginning of the measurements in order to allow their acclimation to the system. The dissolved oxygen in the Erlenmeyer flasks was measured with an YSI 52 oxymeter equipped with a polarographic sensor; organisms were maintained without food for 24 h. A water sample was taken to measure the initial concentration of dissolved oxygen, the water flow in the Erlenmeyer flasks were kept closed for an hour, since according to [37] this time is sufficient to avoid the dissolved oxygen decrease to below 30% not causing stress in the organisms. At the end, new water samples were taken from each Erlenmeyer flask to measure the final concentration of dissolved oxygen. The difference between the initial and final gas concentration was the oxygen consumed by the organisms expressed in mg O2 h−1·g−1 d.w. The oxygen consumption rate (OCR) was calculated according to Equation (1) [38,39].
(1)
where Ct is the change in the oxygen content in the respirometric chambers before and after the test; C0, is the change of the oxygen content in the blank (control); V is the volume of the chamber; and W and T are the wet sea hare biomass in kg and the time duration in hours, respectively. During the experimental phase, a control Erlenmeyer flask was maintained without organisms in order to measure the oxygen consumption and ammonium production of the microorganisms present in the respirometric system, where necessary corrections were made. The Kruskall-Wallis test was used to determine differences between groups [33].
3. RESULT AND DISCUSSION
The thermoregulatory behavior of sea hare Aplysia californica was determined in a horizontal thermal gradient, 18.3˚C for the day cycle and the 20.8˚C for the night cycle (Figure 1); the organisms were never found below 11˚C or above 31.0˚C in trough temperatures. There was a significant difference in preferred temperatures of sea hare between day and night cycles (P < 0.05). [40] also found a diel rhythm in the freshwater gastropod Helisoma trivolis, which had a preferred temperature of 21˚C - 22˚C during night cycle and 17˚C - 18˚C during day cycle. [36] obtained in (Megastrea Lithopoma) undosa an interval of preferred temperature of 16.2˚C in the day cycle and 14.4˚C in night cycle.
Diel differences in preferred temperatures between day and night cycles have been observed in other aquatic ectotherms, including Amia calva [16], Procambarus

Figure 1. Thermoregulatory behaviour of adults of Aplysia californica in a day-night cycle. The shadow zone bordered by triangles represents the 95% confidence interval of the median. The bars include 50% of the distribution and the vertical lines represent the quartiles.
spiculifer [41,42], and Procambarus clarkii. According to [43] this behavior may be explained with the bioenergetics hypothesis, suggesting that diel changes in preferred temperatures may reflect differences in the temperatures favorable for feeding versus temperatures favorable when food is limited. In adults A. californica the diel thermoregulatory behavior obtained indicates that in the habitat of sea hare the food is a limiting factor.
This study concludes that the selection of A. californica of these specific preferred temperatures, which has a correlation with its optimum growth rate and sexual maturity. [5] observed that at 205 days of age the fastest growth rates occurred between 18˚C and 21˚C, and animals reared in the same thermal interval reached sexual maturity at ages 190 to 197 days.
The displacement velocity of sea hare in the gradient showed an initial velocity of 30 cm by hour and gradually diminished to 18 cm per hour, with several fluctuations. We can see a small effect at 18 hours (moment of transition from day to night); the displacement velocity had a very slow activity <5 cm per hour, with an increasing activity at 2:00 a.m. (Figure 2). According to [44], sea hare used an orthothermokinesis orientation mechanism, because it offered the greatest scope of activity [17]. Thus, sea hare remains within a relatively narrow temperature range, decreasing the displacement velocity to 18.0 cm·h−1. [6] observed A. californica in the intertidal zone, where they moved over a distance of 10 m or more per day. This behaviour is correlated with the displacement velocity observed in the horizontal thermal gradient during the diel cycle.
CTMax was evaluated at three temperatures (16˚C, 19˚C, 22˚C), and we observed attachment percentages changed in relation to the acclimated temperature (Figure 3). For A. californica the relationship between percentage of attached organisms on the plate and acclimation temperature is represented as follows.

Figure 2. Displacement velocity of sea hare Aplysia californica in a horizontal thermal gradient, in a day night cycle.
For organisms acclimated to 16˚C this equation was:
(2)
For acclimated sea hare to 19˚C this relationship was:
(3)
As for sea hare acclimated to 22˚C this relationship was:
(4)

Figure 3. Percentage of sea hare Aplysia californica remaining attached plotted against different acclimation temperatures. The continuous line denotes the critical thermal maxima at 50%.
where y = temperature and x percentage of attached organisms. For sea hare acclimated at 16˚C, the last organism to lose holding was registered at 31˚C. In the case of the animals who were maintained at 22˚C, the sea hare that ultimately had a loss of attachment reaching a temperature of 39˚C (Figure 3). As for the CTMax (defined by the moment where 50% of the sea hare remained attached), the temperatures were 32.7˚C, 33.6˚C and 36.2˚C to 16˚C, 19˚C and 22˚C acclimated temperatures respectively (Figure 3). After the exposure to determine CTMax, organisms were returned to their acclimation temperature for 96 h, and they had a 98% survival rate.
In this study, the CTMax of Aplysia californica was determined by a methodology that involves dynamic changes in temperature. [30] presented a summary of the critical thermal maxima at 50% for various species of gastropods of genus Haliotis, acclimated to different temperatures and salinities. In the determination of CTMax to 50%, sea hare was subjected to a continuous and constant linear increase in temperature until an endpoint was reached; when we changed acclimation temperatures from 16˚C to 19˚C, we observed an attachment increase of 3%, and from 16˚C to 22˚C, thermal tolerance increased in 10%.
When Aplysia went through the process of reaching the end point of CTMax it had several particular responses. It maintained the mantle spread and stacked it to the acrylic plate during the first minutes, later the mantle looked deployed and well spread. When sea hare was about to arrive to the ending point, its mantle started to curl causing the body to rise and looked like a balloon. Finally at the point of losing attachment its mantle displayed many ruffles, this process continued to intensify until the mantle had a complete loss of adhesion; the endpoint of CTMax is marked when Aplysia starts to float.
According to [45,46] the attachment capacity in marine gastropods is the result from the sum of two main components: the proteinaceous or mucopolysaccharide secretions and the muscular activity that maintains the suction force from the foot and its shape. It also has been shown that temperature affects the capacity of attachment of different species such as sea urchins, abalones, sea cucumbers, and wavy turban snail [30,36,47-49].
Oxygen consumption rate (OCR) in A. californica was closely related to acclimation temperature. At 16˚C and 19˚C sea hare had a relatively stable metabolic rate with a slight variation between 5 to 6 mg O2·h−1·kg−1 w.w., and when acclimation temperature increased from 19˚C to 22˚C, OCR increased to 9 mg O2·h−1·kg−1 w.w. (Figure 4). Once the sea hare was in a temperature of 22˚C there was an increase of 150% of oxygen consumption rate, indicating the high need of oxygen supply requirements. In the case of sea cucumber, Apostichopus japonicus, when exposed to higher temperatures [24] it

Figure 4. Oxygen consumption rate of sea hare Aplysia californica acclimated at different temperatures. The shadow zone bordered by circles represents the 95% confidence interval of the median. The bars include 50% of the distribution and the vertical lines represent the quartiles.
obtained a similar response in OCR; according to [24,50] this pattern suggests that high thermal tolerance may come with a significant fitness and energetic cost. In low temperatures (16˚C to 19˚C) Aplysia californica reduced their oxygen consumption and had an increment on physiologically useful energy, enhanced food conversion, and favored the saving of metabolic energy that can be used for growth, eating or reproduction activities [5].
4. ACKNOWLEDGEMENTS
We thank Jose M. Domínguez and Francisco Javier Ponce from the Drawing Department of CICESE for preparing the figures; Asael G. Arroyo Re for the English reviews of the manuscript.