Dataset: Nurse shark respirometry: tank estimates of respiration and Overall Dynamic Body Acceleration (ODBA) from Mote Marine Lab, 2014 (Shark Metabolic Rate project)

These methods describe all those associated with the project.

Capture and maintenance:

Respirometry experiments were conducted on juvenile nurse, lemon, and blacktip sharks. Nurse sharks (n=8, 53-132 cm TL) were captured via rod and reel from the Florida Keys, USA. Lemon sharks (n= 30, 69-100 cm TL) were captured with cast nets from Cape Canaveral, FL, and the Florida Keys, USA. Blacktip sharks (n = 8, 53-64 cm TL) were captured with rod and reel from Terra Ceia Bay, FL, USA. All animals were transported to Mote Marine Laboratory in Sarasota, FL, and held in a 150,000 liter indoor, recirculating tank for the duration of experiments. Sharks were fed a diet consisting of herring, squid, and shrimp every other day to satiation, but were fasted prior to the beginning of trials to achieve a post-absorptive state. Nurse sharks were fasted for at least 72 h prior to trials, and lemon and blacktip sharks were fasted for at least 48 h prior to trials. All sharks were kept on a constant 12 h light:dark cycle.

Respirometry trials were run at two temperature groups representing the low (~20°C) and high (~30°C) end of the temperature range these species are likely to experience in the wild. Sharks were acclimated to trial temperatures in the holding tank for at least a week prior to experimentation. Trials with lemon and blacktip sharks were all run within 1 –2 months of initial capture. Nurse shark trials were run with individuals that had been maintained in captivity for at least one year.

Accelerometry:

During trials, sharks were equipped with Cefas G6A+ acceleration data loggers (Cefas, Inc., Lowestoft, UK), which recorded triaxial acceleration at 25 Hz. Tags were attached to the first dorsal fin of sharks at two points using monofilament (Fig. 4) at least 24 h in advance of the start of a trial. Since sharks tagged in the field would also need to be tracked acoustically in order to retrieve the data loggers, the loggers used in captive trials were paired with a mock acoustic tag (Vemco Inc., Nova Scotia, model V9) in order to maintain the same weight and drag as tags used in field studies (see Fig. 4). The paired tag package measured 37 x 36 x 15 mm and weighed 23 g in air, representing 0.02-0.002% of the body mass of the study animals.

Respirometry:

Trials were conducted in a circular, closed respirometer constructed from a modified fiberglass holding tank, as described in Whitney et al. (2016). Briefly, the respirometer was sealed using a lid constructed from a PVC frame with translucent plastic sheeting stretched across it, and dissolved oxygen (DO) and temperature were measured using a handheld YSI (model Pro 2030, Yellow Springs, OH). In order to ensure even water mixing in the respirometer and create water flow past the YSI probe for accurate DO measurements, a pump was set up in the centre of the tank facing into a T-shaped pipe made of PVC which housed the YSI probe, providing sufficient water movement to mix water throughout the system without creating a current in the tank. The pump and YSI were enclosed in a circular cage made of PVC and rigid plastic mesh during lemon and nurse shark trials to encourage the sharks to swim in full circles around the outer edge of the tank. This mesh cage was not used during blacktip trials because it appeared to induce stress in blacktip sharks. See Fig. 4 for a picture of the respirometer set-up.

Lemon and nurse sharks were placed into the respirometry system at least 12 h prior to the start of trials to allow them to acclimate to the system overnight. Blacktip sharks appeared to fatigue after extended periods in the smaller tank system, and were acclimated to the respirometer for 1 h prior to the start of trials rather than overnight. After the acclimation period, the respirometer tank was isolated from its flow-through system and sealed off with the lid. The tank was surrounded by a curtain to limit extraneous disturbances, and the trials monitored remotely using a live digital video feed. DO and water temperature were recorded every five minutes, and behavior monitored constantly throughout the trials. Trials began with the DO near 100% saturation and were run until the DO reached 80% saturation. To assess background respiration, a blank respirometer (without an animal) was measured for 4 h during each group of trials.

Data analysis:

Periods of the trials where sharks displayed consistent behavior (either constant swimming or resting) for at least 20 min, were used to analyze metabolic rate. Oxygen consumption rate (MO2, mgO2 kg-1 h-1) was calculated for each one of these periods using equation 1:

 (3)   MO2 = (S-b)(60 min)(V)/(W)

Where S is the slope of the oxygen degradation curve (in min), b is the slope of the background respiration curve, V is the volume of the respirometer in liters, and W is the weight of the shark in kg. The volume of the shark (<10 l) was considered to be negligible relative to the respirometer volume (2494 l), representing an error of <0.5%, and was thus not incorporated into our model.

Accelerometer data were analyzed using Igor Pro (Wavemetrics, Lake Oswego, OR). Static acceleration (indicating animal body position) was separated from dynamic acceleration (indicating animal movement) using a 3 s box smoother (Shepard, 2008), and ODBA calculated as the sum of the three dynamic acceleration axes. A mean ODBA value was calculated for each time interval where MO2 was analyzed during respirometry trials, to provide paired ODBA-MO2 points for analyses.

SMR's (Standard Metabolic Rate) were calculated for lemon and nurse sharks at each experimental temperature by averaging metabolic rates during all resting intervals recorded. Since blacktip sharks are a ram-ventilating species, SMR was not directly calculated, but was estimated using the intercept of the ODBA-MO2 relationship. RMRs were calculated for each species and temperature grouping as the mean metabolic rate of all periods where the study animal showed consistent swimming activity. Rest periods for lemon and nurse sharks were not included in RMR calculations.

Due to difficulties in keeping blacktip sharks in captivity for extended periods, a full ODBA -MO2 calibration was not conducted at temperatures near 20°C. However, a pair of trials was run at 21.6°C on one individual, allowing for the calculation of a preliminary Q10 value. Since there were insufficient data at the low temperatures to extrapolate an SMR for blacktip sharks, this Q10 value was calculated based on differences in RMR. Only RMR data that overlapped in ODBA levels between the two temperatures were used in the Q10 calculation to ensure the comparison was made between metabolic rates from similar activity levels. Q10 values for nurse and lemon sharks were calculated using SMR data. Q10 values were calculated according to the Van’t Hoff equation

 (4)   Q10 = (R2/R1)^^10/T2-T1

Where R1 is the metabolic rate at temperature T1, and R2 the metabolic rate at temperature T2.

Error estimation and modeling:

Error estimation and modeling were performed in R (R Core Team, 2010), using the lme4 (Bates et al. 2015) and MuMIn (Barton, 2016) packages. Linear mixed models were constructed to describe the relationship between ODBA and oxygen consumption for each species, with ODBA, activity state (resting or swimming), and temperature group as predictor variables, and individual included as a random effect. All combinations of variables and interaction effects were assessed using a comprehensive model selection table produced by the MuMIn package. Models were compared against each other using a corrected Akaikes Information Criterion (AICc), residuals, log likelihood, and R2 of the models. Using the regression line produced by the best-fit model, the predicted MO2 for each analysis interval was calculated, and compared against the measured MO2 to calculate the standard error of the estimate for each species [(predicted MO2 – observed MO2) / predicted MO2 x 100], which was examined as a percentage of the measured MO2. Normality of the residuals of the optimal models was tested using an Anderson-Darling test (Wright et al., 2014).]