Haematological and biochemical reference intervals for wild green turtles (Chelonia mydas): a Bayesian approach for small sample sizes

Abstract Animal health is directly linked to population viability, which may be impacted by anthropogenic disturbances and diseases. Reference intervals (RIs) for haematology and blood biochemistry are essential tools for the assessment of animal health. However, establishing and interpreting robust RIs for threatened species is often challenged by small sample sizes. Bayesian predictive modelling is well suited to sample size limitations, accounting for individual variation and interactions between influencing variables. We aimed to derive baseline RIs for green turtles (Chelonia mydas) across two foraging aggregations in North Queensland, Australia, using Bayesian generalized linear mixed-effects models (n = 97). The predicted RIs were contained within previously published values and had narrower credible intervals. Most analytes did not vary significantly with foraging ground (76%, 22/29), body mass (86%, 25/29) or curved carapace length (83%, 24/29). Length and body mass effects were found for eosinophils, heterophil:lymphocyte ratio, alkaline phosphatase, aspartate transaminase and urea. Significant differences between foraging grounds were found for albumin, cholesterol, potassium, total protein, triglycerides, uric acid and calcium:phosphorus ratio. We provide derived RIs for foraging green turtles, which will be helpful in future population health assessments and conservation efforts. Future RI studies on threatened species would benefit from adapting established veterinary and biomedical standards.


Introduction
Blood analyses are routinely used for conservation, ecology and rehabilitation purposes and can indicate population declines and long-term survival challenges (Seminoff and Shanker, 2008;Hamann et al., 2010;Perrault et al., 2017;Stacy and Innis, 2017;Perrault et al., 2021). Indicator species, such as green turtles (Chelonia mydas), help in assessing the threats a particular ecosystem or habitat is facing Lutz, 2004, De Cáceres et al., 2010). Ecosystem functions and services can be assessed by species morphology, behaviour, demography, physiology, biogeochemical composition and socioeconomic importance (Castro Tavares et al., 2019). Green turtles contribute to ecosystem functioning in foraging grounds and nesting beaches by transporting significant amounts of nutrients from nutrient-rich foraging grounds to nutrient-poor nesting beaches (Bjorndal and Jackson, 2002). Although considered migratory in their early life stages and during breeding seasons, green turtles show strict fidelity to foraging grounds as small as 2 km 2 (Musick and Limpus, 1997;Shimada, 2015) and can reflect the qualitative status of their local habitat. Seagrass meadows are considered essential carbon storage and sequestration sites, and green turtles maintain nutrient-rich areas and contribute to the biodiversity of seagrass species through grazing and seed dispersal (Duarte et al., 2010, Fourqurean et al., 2012, Scott et al., 2020. A comprehensive understanding of a species' baseline information is required for correctly interpreting haematological and biochemical data and includes assessing physiology and anatomy, reproductive biology, behaviour, food habits and nutritional requirements, home range and expected parasite fauna (Ryser-Degiorgis, 2013). Interpretation of health status is also dependent on the comparison of blood analyte values with suitable reference intervals (RIs). The utility of RIs for each analyte relies on methodological, physiological and environmental factors. General guidelines for the development of RIs in healthy animals are available for species commonly encountered in the veterinary profession (Friedrichs et al., 2012), and a recent systematic analysis is available for nondomestic species (Moore et al., 2020). These guidelines recommend as selection criteria to include a sufficiently large sample size, typically 100-200 individuals for normally distributed data, and to provide statistical reliability when using frequentist statistics (Jennen-Steinmetz and Wellek, 2005;Geffré et al., 2009;Friedrichs et al., 2012;Wellek et al., 2014;Klee et al., 2018;Moore et al., 2020). Specific guidelines for establishing RIs in sea turtles are also available (Page-Karjian and ; however, some aspects (e.g. the recommended sample size) differ from standards approved for human or veterinary medicine, such as those of the Clinical Laboratory Standards Institute (CLSI) guidelines or the American Society of Veterinary Clinical Pathology Guidelines (CLSI, 2010;Friedrichs et al., 2012;Wilkinson et al., 2016).
Frequentist statistical methods, which are traditionally used in wildlife research, are best suited to large sample sizes (n > 120) and low variability between samples (Friedrichs et al., 2012;Moore et al., 2020). By nature, threatened species research is limited by small sample sizes due to population size and distribution, species biology and behaviour, restricted funding and resources, permitting limitations or short project timeframes (Steidl et al., 1997, Bissonette, 1999, Lloyd-Smith et al., 2005, Kéry, 2010, Harden et al., 2018. Only a limited number of studies on vertebrate wildlife have collected samples from >120 individuals (see reviews on the topic by Cray, 2015, Moore et al., 2020, which is also true for green turtle biochemical and haematological studies (Supplementary Table S1). In addition, blood values in reptiles may be highly variable across species, populations, sex and life stages, and this has been related to the lack of robust RIs (Stahl, 2006;Mitchell and Tully, 2008;Stacy and Innis, 2017). Following guidelines for developing RIs in threatened species is therefore problematic, and unreliable RIs may lead to false interpretations on population health. This may directly hinder conservation, management and rehabilitation efforts by, for example, overseeing abnormal findings within a population that may lead to falsenegative diagnoses and inadequate enforcement measures (Deem and Harris, 2017;Sacchi et al., 2020). As a result, an unnoticed decline in population health will reveal a decline in reproductive output and/or population viability, and will lead to increased caseloads in rehabilitation centres (Hamann et al., 2010;Commonwealth of Australia, 2017;Deem and Harris, 2017). Statistical approaches that account for small sample sizes would therefore provide an increased reliability and clinical utility in the determination of RIs (Steidl et al., 1997, Harden et al., 2018, Sacchi et al., 2020.
While comprehensive and representative sampling and prioritizing large sample sizes remain important considerations, alternative statistical approaches such as Bayesian statistics effectively account for small sample sizes (van de Schoot et al., 2021). Bayesian statistics are very popular in the biomedical and ecological sciences, as they update the probability for a hypothesis as more data becomes available (van de Schoot et al., 2021). Bayesian models are based on the available data and account for individual variation in the calculation of predictive credible intervals, which are equivalent to frequentist 95% confidence intervals (CIs; i.e. mean ± two standard deviations) (Hespanhol et al., 2019). As a result, Bayesian models can be used to minimize the number of individuals included in a study (Katki et al., 2005, Sottas et al., 2011. These characteristics make Bayesian statistics highly suitable for threatened species research. The aims of this study were to measure biochemical and haematological analytes of green turtle foraging aggregations (C. mydas, n = 97), to calculate RIs that could be used for the examined regions and to compare the predicted intervals against other studies. The turtles were sourced from two geographically and ecologically distinct foraging  grounds in North Queensland, Australia (Howick Group of Islands and Townsville region; Fig. 1). The examined foraging grounds were in a marine-protected area (MPA) with very little anthropogenic impacts (Howick Group of Islands) and in an MPA located in an industrialized region currently experiencing a port expansion (Townsville region) , Queensland Government, 2021. Green turtle grazing has been reported at both sites , Flint et al., 2019. Howick Group of Islands is found in the northern Great Barrier Reef Marine Park (−14.416695 • S 144.880484 • E), ∼30 km from the Cape York region catchment and consists of mid-shelf, unpopulated reefs. The area is considered to be free from chemical pollutants, fishing pressure and coastal development (Villa et al., 2017, Flint et al., 2019. In contrast, Townsville region has an estimated population of >230 000 (Australian Bureau of Statistics, 2022) and is influenced by anthropogenic impacts such as industrial runoff, urbanization and coastal dredging (Villa et al., 2017). Based on past studies and the threatened status of this species, we anticipated sampling limitations, and were interested in using statistical methods suitable for small sample sizes (n < 120). We aimed to develop Bayesian linear mixed-effects models that would account for the effects of low sample size, geographical location, length and mass on the selected analytes. Further objectives of this study were to compare the resulting intervals in wild turtles between industrialized versus offshore foraging grounds.

Study sites
This study was conducted in two major foraging grounds in North Queensland, Australia: (i) Cleveland Bay (19O13 05 S, 146O55 19 E)  are separated by over 500 km. Sampling was conducted exclusively in winter (between June and August 2019) to avoid travelling during the cyclone season.

Animals and sampling protocol
Haematological and biochemical analyte values were determined from plasma obtained from wild turtles (n = 97) captured from a boat using rodeo technique (n = 85) as described in Limpus and Reed (1985) or hand-captured in shallow water (n = 12). Sampling was opportunistic and tidedependent and predominantly took place in the mornings. A general health assessment was conducted by performing a physical examination to record any injuries, epibiota or presence of tumoral lesions (e.g. fibropapillomatosis) following standard procedures outlined in Deem and Harris (2017) and Harris et al. (2017). Only assessed healthy turtles, without macroscopic anomalies, were selected for further examination. Where possible, the animal's eyes were covered with a cloth to minimize stress. Turtles were tagged with approved titanium identification tags and curved carapace length (CCL) from notch to tip to the nearest 2 mm was measured. Turtles were allocated into life stages based on CCL as per Chaloupka and Limpus (2001) with juveniles (immature) CCL < 65 cm, sub-adults (immature) 65 cm > CCL < 90 cm, and adults (mature) CCL > 90 cm. Body temperature was measured using a thermocouple (8402-20 Thermistor 237 Thermometer, Cole-Palmer Instruments, Vernon Hills, IL, USA), and by inserting the probe 5 cm into the cloaca (Flint, 2013;Stacy and Innis, 2017). Total body mass measurements were recorded using a specially designed harness, which secured each animal around the base of each limb. The harness was then attached to a digital scale where the mass was measured to the nearest 0.1 kg while the animal remained suspended. The harness was removed immediately after weighing. Blood samples were taken from all turtles as described below. Once sampling was completed, the turtles were released in the same area they were captured. Randomly selected juvenile turtles (14.4%, n = 14/97) were also assessed by laparoscopic examination to determine their sex. This standard procedure was conducted last and was part of a longitudinal monitoring study conducted on a yearly basis . All turtles were tagged, measured and weighed following standard operating procedures (DBCA, 2017, DES, 2018. The protocol related to measurements and health assessments was standardized and occurred in the following order: capture, physical examination, measurements, blood sampling and laparoscopies (on selected animals). All procedures and protocols were approved by the Great Barrier Reef Marine Park Authority (permit number G19/42769.1) and the Department of Environment and Science, Queensland Government (permit numbers SPP18-001167 and PTU18-001419-2).

Blood sampling and processing
Blood was sampled from the external jugular vein, which is located on the superficial, lateral regions of the neck. Prior to venepuncture, the skin was disinfected using 70% ethanol swabs (Liv-Wipe, Livingstone, Livingstone Int., Mascot NSW, Australia). Blood samples (2 ml) were collected using a 10ml syringe (Shandong Hapool Medical Technology Co., Heze, China) with a 22-gauge × 1 1 / 2 inch needle (Terumo, Japan). No expected or unexpected adverse events occurred. Sample quality was assured by immediate visual inspection of each blood sample. Any sample suspected of contamination with lymph fluid was discarded and an additional sample was collected.
Packed cell volume was determined as an indicator for hydration state and anaemia and was measured twice to determine the average value (Livingstone Microhaematocrit Capillary Tubes, Livingstone Int., Mascot, NSW, Australia; and Pico 17 Microcentrifuge, Thermo Fisher Scientific, Waltham, MA, USA). Duplicate blood smears were prepared using a clean glass slide for the smear and as the spreader slide (Thermo Scientific Menzel-Gläser, Thermo Fisher Scientific, Waltham, MA, USA). The remaining blood was transferred to a sodium-heparin-coated blood collection tube (BD Vacutainer LH 34 I.U., BD Vacutainer Systems, Plymouth, UK), which was gently rocked to ensure proper mixing of blood components. Smears were initially fixated with methanol and were stained once these were returned to the laboratory. Blood smears were interpreted from turtles captured in Townsville region (blood smear quality from turtles captured at Howick Group of Islands was insufficient). Pre-analytical errors may influence analyte values and therefore the RIs determined. Pre-and post-analytical procedures were standardized and followed recommendations for field sampling techniques of reptilian blood (Fullarton, 2012, Eshar et al., 2018. Samples were kept at 4 • C (39.2 • F), either using refrigeration when available or a cooler box with ice packs for up to 12 hours before centrifugation, and blood tubes were prevented from direct contact with ice packs. Blood smears were stained (Diff Quick and Wright's stain) and examined using a light microscope (Olympus BX43, Olympus Corp., Tokyo, Japan) at 40× magnification following standard procedures. Blood smears were blindly assessed (JCU Veterinary Diagnostic Pathology Laboratory, Townsville, Queensland, Australia), and leukocyte identification was determined upon consensus ( Fig. 2). A white blood cell (i.e. leucocytes) differential count was performed on at least 150 cells, and the cells classified as heterophils, lymphocytes, monocytes, eosinophils or basophils Ebanks, 1984, Samour et al., 1998). The heterophil:lymphocyte (H:L) ratio was also determined.

Statistical analyses
All statistical analyses were produced with R statistical software, using the package ggplot2 for data visualization (Hadley, 2016; R Core Team, 2019). The statistical approaches used in our study were based on the methods used by Logan (2020) Sacchi et al. (2020) and Spinks et al. (2021). Distribution of the response variables (i.e. the biochemical and haematological analytes) were either Gaussian or Gamma, and log-transformed models were considered (normality or non-normality results for each analyte are described in Table 3 and were based on the best model fit). Bayesian generalized linear mixed-effects models were developed for all biochemical and haematological variables, except for sodium and chloride, which were assessed using Bayesian generalized additive models (best model fit). The models were fit using uninformative normal priors or with weak informative priors to allow for regularization whenever a more informative prior was required (Korner-Nievergelt et al., 2015). The posterior prior was derived from the prior distribution, and suitability was confirmed with visual posterior checks. Models were run with the No-U-Turn sampler, using three chains and 5000 iterations. The first 1000 iterations were discarded to converge the model to the correct posterior distribution.
Models were fitted separately for each response variable (i.e. biochemical or haematological analytes). The response variables were first explored graphically and were then statistically analysed by fitting the models previously mentioned. We included location, mass, and CCL as fixed effects, and animal ID as a random effect to account for inter-animal variability. Collinearity between mass and CCL in the studied locations is very common, especially in mature turtles . We included both variables into our calculations to account for exceptions where collinearity might not be the case (e.g. young turtles with an increased growth rate, or turtles that vary in body condition for the same CCL) (Eckert et al., 1999). In addition, the uninformative and weak informative priors used in Bayesian statistics help reparametrizing the model, accounting for collinearity within the data (Ogle and Barber, 2020). The resulting predictions were then back transformed, when applicable, to obtain the final RIs in their original scale. The predicted values for each parameter are reported as estimated marginal mean (EMM), and as lower and upper highest posterior density credible intervals (HPD-CIs) (Table 3), which are analogous to frequentist CIs (Lee, 1989). HPDCI and CI only differ in the way the predicted parameter is treated, i.e. Bayesian HPDCI treats the predicted parameter as a random variable, whereas frequentist CI treats it as a fixed variable.
All models were fit in a Bayesian analytical framework available in the packages rstanarm (Goodrich, 2020), brms (Bürkner, 2017(Bürkner, , 2018 and gamm4 (Wood et al., 2017). Model assumptions (e.g. linearity and homogeneity of variance) were visually confirmed with diagnostic residual plots, all of which were satisfactory, using the packages coda (Plummer et al., 2006), bayesplot (Gabry and Mahr, 2021), ggmcmc (Fernández-i-Marín, 2016) and DHARMa (Hartig, 2020). The final model selection was based on diagnostic residual plots (e.g. DHARMa residual plotting, Hartig, 2019), on the fit of the data to the selected model and on the corrected Akaike Information Criterion for small sample sizes (AICc, Barton and Barton, 2015). Outlier identification and exclusion was performed with residual plotting using the package DHARMa (Hartig, 2019), and negative analyte values were excluded prior to running the models. Sample size estimates using G * Power analysis revealed a total sample size of 159 turtles to achieve a Power of 0.8 (effect size 0.25, α 0.05, three groups). However, this estimate relates specifically to frequentist statistical approaches, since Bayesian methods do not assume fixed/known effect sizes. Posterior prior distributions were derived instead, all of which were satisfactory.
Specific contrasts were conducted for comparisons across locations with the package emmeans (α = 0.05) (Lenth, 2016). Posterior probability distributions using the Markov Chain Monte Carlo (MCMC) estimation assessed the effects of location, mass and CCL on the measured analytes (Fernández-i-Marín, 2016). The differences in the parameter intervals were based on 95% Bayesian Uncertainty Intervals (UIs) for modelled higher posterior density (HPD) median effects. Statistical significance (P < 0.05) was inferred when the 95% UIs did not overlap. Whenever referring to location differences throughout the manuscript, it should be noted that mass and CCL were accounted for in the specific contrasts.
Correlations between variables were assessed using Pearson's and Spearman's correlation coefficient analyses (strong correlation assumed when P < 0.05 and r > 0.5, Supplementary Appendix S1). Effect size indexes (Hedges' g) were calculated where possible for comparison with other studies (Supplementary Table S3). Additional body condition indices (BCI) were determined by converting straight carapace length from the measured CCL values (Bjorndal and Bolten, 1989;Bjorndal et al., 2000;Norton and Wyneken, 2015). This study followed recommended human and veterinary guidelines (Fig. 3) during the data collection and analysis process to ensure the reliability of the established RIs, i.e. CLSI guidelines, FAIR principles, ARRIVE Guidelines (Supplementary Appendix S3) and American Society of Veterinary Clinical Pathology Guidelines (CLSI, 2010, McGrath et al., 2010, Friedrichs et al., 2012, Wilkinson et al., 2016, Percie du Sert et al., 2020.

Animal characteristics
A total of 97 wild turtles were captured. Of these, 26% were adults, 8% subadults and 66% juveniles (Table 2). Laparoscopic examination of a subset of the sample revealed 12 female and 2 male juvenile turtles (n = 14/97). The different distribution of life stages across the two sites was reflected in sample distribution, with adults (n = 25) and subadults (n = 8) only caught at Howick Group of Islands, and the Townsville region group consisting entirely of juveniles (n = 40). All turtles appeared healthy on physical examination, were in good body condition and had no apparent external lesions, except from one animal which was missing a front limb (ID number QA94686). Mean CCL of all turtles was 65.5 cm (range, 36.8-115.2 cm), and mean mass was 43.4 kg (range, 6.15-147.2 kg). Mean BCI was 1.16 (range, 0.88-1.53). No significant differences in animal characteristics (i.e. CCL, mass, BCI and cloacal temperature) between sites were identified (Table 2). Mass and CCL were strongly and positively correlated (Pearson's correlation coefficient 0.97, P < 0.05, t = 43, df = 103). Animal data for each study site and life stage are provided in Table 2.

RIs for blood biochemical and haematological parameters
RIs (EMM, 95% upper and lower HPDCI limits) are reported in Table 3 (measured, original data are reported in Supplementary Table S2). The majority of blood analytes were not statistically different in the turtles across locations, mass or CCL (76%, 86% and 83%, respectively) ( Table 3).
Location was associated with significant (P < 0.05) differences in 7/29 (24%) blood analytes: albumin, cholesterol, potassium, total protein, triglycerides, uric acid and calcium:phosphorus ratio (Table 3). Mass and CCL were associated with eosinophil percentage and H:L ratio (P < 0.05). Plasma levels of alkaline phosphatase and urea were also influenced by CCL and mass (Table 3). An exception was that CCL had a significant effect on aspartate transaminase (P < 0.05), whereas mass did not. These findings

Research article
Conservation Physiology • Volume 10 2022 North Queensland, Australia (n = 97), which were derived using Bayesian predictive modelling. The Bayesian generalized linear mixed-effects models accounted for the effects of low sample size (n < 120), geographical location, length and mass. RIs for wild turtles were predicted by including both locations into the Bayesian model, and by accounting for potential differences across them (Table 3).
Our predicted intervals were narrower and within previously reported values or intervals that had been calculated using frequentist statistics (Aguirre and Balazs, 2000, Hamann et al., 2006, Whiting et al., 2007, Arthur et al., 2008, Flint et al., 2010 (Supplementary Table S1). A wide range of factors are known to affect blood analytes in sea turtles, including geographical location, diet composition, sex, age of maturity (i.e. mass and length), captivity, season or weather conditions and sample handling and processing (Herbst and Jacobson, 2002, Hamann et al., 2006, Drake et al., 2017, Stacy and Innis, 2017, Harden et al., 2018, Sacchi et al., 2020. In our study, most haematological and biochemical analytes had no significant association (P > 0.05) with location, mass or CCL (76%, 86% and 83%, respectively), with some exceptions detailed below (Table 3). Other studies, most of which predominantly sampled immature green turtles as well (31/37 studies, 84%, Supplementary Table S1), reported significant effects of CCL and/or mass on the measured blood analytes (Bolten and Bjorndal, 1992, Hasbun et al., 1998, Labrada-Martagon et al., 2010. The difference in the impact of mass and CCL on blood analytes between the current and previous studies may be related to differences in the statistical treatment of our data, location, seasonality and/or diet (Stacy and Innis, 2017).

Statistical approach
Clinical guidelines recommend establishing RIs with a large enough sample size (n > 120) and using predictive statistical models (e.g. linear mixed effects models or Bayesian statistics) to minimize variability within and between analytes (Katki et al., 2005, CLSI, 2010, Sottas et al., 2011, Friedrichs et al., 2012, Ozarda, 2016, Harden et al., 2018, Sacchi et al., 2020. Most veterinary studies refer to the American Society of Veterinary Clinical Pathology (ASVCP) Guidelines (Friedrichs et al., 2012), which have also been promoted by sea turtle researchers (Page-Karjian et al., 2015;Stacy and Innis, 2017;Stacy et al., 2019;de Mello and Alvarez, 2020;. Other recommendations include sampling a minimum of 20 animals to establish RIs, with larger sample sizes preferred to calculate more reliable results (Page-Karjian and . From a clinical perspective, however, Bayesian models have the advantage of accounting for small sample sizes and overcome important limitations of frequentist likelihood models, such as biassed maximum likelihood estimates (Katki et al., 2005, van de Schoot et al., 2021. Bayesian statistics also have the ability to incorporate independent information about both fixed and random factors and to fit models when complex and multiple interactions exist between variables (van de Schoot et al., 2021). For example, Bayesian models have been used to establish haematological RIs in lizards (Sacchi et al.,   Results for tests of significance (P) are displayed for the effects of location (i.e. Howick Group of Islands versus Townsville region), mass and CCL. Analyte, unit, sample size (n), distribution (P: parametric, NP: non-parametric), EMM, % higher posterior density credible intervals (HPDCI) lower and upper limits (analogous to mean +/− 2 SD) and results of the tests of significance are reported. The differences in the parameter ranges were based on 95% Bayesian UIs for modelled HPDCI. Statistical significance was assessed with posterior probability distributions using the MCMC estimation. Statistical significance (sig) was inferred when the 95% UIs did not overlap (P < 0.05). 2020), identify abnormal biochemical analytes in veterinary medicine (Knox et al., 1998) and predict wildlife population declines over time (King et al., 2009). We present our Bayesian modelling as an example for establishing robust RIs in green turtle studies limited by small sample sizes.
A paradigm shift to develop standardized procedures for sea turtles specifically, and for threatened species research in general, has been called out by several authors and organizations (Lawson et al., 2021;Stacy and Innis, 2017;Mashkour et al., 2020;, Ryser-Degiorgis, 2013, Stokes et al., 2010. Failure to achieve this strategic priority will result in increased false-positive and falsenegative diagnoses and unreliable population health estimates. Ultimately, evidence-informed rehabilitation and conservation efforts will be enhanced by accurate and representative RIs. Threatened species studies would therefore benefit from adapting established veterinary and biomedical standards, such as the ASVCP Guidelines (Friedrichs et al., 2012). However, if the recommended sample size (n > 120) cannot be reached, using alternative predictive approaches such as Bayesian statistics is strongly encouraged. Previous studies using Bayesian modelling frameworks used or recommended sample sizes ranging from 20 for RIs predictions in box turtles (Terrapene ornata) (Harden et al., 2018) and 36 for estimating mortality rates in alligator snapping turtles (Macrochelys temminckii) (Steen and Robinson Jr, 2017) to 100-140 for sex ratio predictions in loggerhead turtles (Caretta caretta) (Shertzer et al., 2018). The minimum sample size to be used with Bayesian models for the determination of RIs should be confirmed with prior predictive checking, which is particularly relevant in complex models with small sample sizes (van de Schoot et al., 2021).

10
RI studies should always be carefully designed and aimed for collecting randomized data from as many individuals as possible. Sampling bias might still occur in Bayesian statistics, and the sampled individuals might not accurately represent the population if the sample size is too small. For further information and recommendations on Bayesian modelling, prior selection and sample size, we refer the interested reader to a recent review by van de Schoot et al. (2021) and to a comprehensive RI study on lizards (Sacchi et al., 2020).

Haematological analyses: White blood cell differential counts
Reptilian leukocytes are considered indicators for systemic stressors, with heterophils fulfilling the surrogate role of neutrophils in lower vertebrates (Campbell, 2006, Campbell, 2015, Stacy and Innis, 2017, Flint et al., 2019. Heterophils seem to have similar functions to those found in avian blood, as they rely on oxygen-independent mechanisms to combat microorganisms (Stacy et al., 2011, Campbell, 2015. We compared our leukocyte percentages (%) with those previously reported for green turtles (Samour et al., 1998, Flint et al., 2010, Lewbart et al., 2014, March et al., 2018 and calculated effect size indexes (Hedges' g) and % difference in mean values where possible (Supplementary Table S3). The H:L ratio is a reliable method to estimate stress responses in vertebrates (Davis et al., 2008, Krams et al., 2012. Elevated H:L ratios may reflect physiological differences between green turtle aggregations, or indicate a sub-clinical, inflammatory response (Davis et al., 2008, Goessling et al., 2015. Globally, green turtles are reported to have a low H:L ratio (Lewbart et al., 2014;Muñoz-Pérez et al., 2017); but the opposite trend has also been found in the United Arab Emirates and Australia (Samour et al., 1998, March et al., 2018. High H:L ratio has also been found in loggerhead turtles (Caretta caretta) in the Atlantic Ocean (Casal et al., 2009, Deem et al., 2009, Kelly et al., 2015. The H:L ratio was approximately 1:1 in our turtles. Only two Australian studies to date reported white blood cell counts, both in locations >1000 km South of our field sites; Flint et al. (2010) reported a 1:3 H:L ratio for Southern Queensland turtles (Australia) and March et al. (2018) reported 4:1 and 2:1 ratios for rehabilitating green turtles in New South Wales (Australia). In our study, mass was found to have a negative effect on the H:L ratio and a positive effect on the eosinophil percentages (%) ( Table 3). The various relationships observed suggests age-related changes and increasing exposure to environmental stressors or infectious agents to influence H:L ratios (Aguirre et al., 1995, Deem et al., 2009, Oh and Hustead, 2011, Muñoz et al., 2013. This finding is also observed in humans who may experience a dominance shift of lymphocytes to neutrophils with ageing (Li et al., 2015). The increasing eosinophilia in the turtles may reflect a decrease in heterophils, or an increasing parasite burden as the turtles age (Aguirre et al., 1995, Deem et al., 2009, Muñoz et al., 2013. Eosinophil percentages (%) of the sampled turtles were lower than those previously reported (Samour et al., 1998, Lewbart et al., 2014. PCV of the turtles included in our study were within previously reported intervals for green turtles (Flint et al., 2010;Lewbart et al., 2014).

Biochemical analytes
The established biochemical intervals for wild turtles fell within previously reported blood values, ranges or intervals for green turtles in Australian waters and elsewhere (Bolten and Bjorndal, 1992, Aguirre and Balazs, 2000, Hamann et al., 2006, Whiting et al., 2007, Arthur et al., 2008, Flint et al., 2010, March et al., 2018, with some exceptions (effect size and mean % difference are outlined in Supplementary Table  S3). Our RIs were narrower than those previously defined for the same regions (Flint et al., 2010), which we attributed to the differences in the statistical methodology used. In this section, we decided to focus on the analytes that differed across studies and refer the interested reader to Stacy and Innis (2017) for a detailed summary on clinical pathology in sea turtles.
Glucose-Sick animals often present hypoglycaemia or hyperglycaemia, which is usually associated with a stress response (Innis et al., 2009;Stacy and Innis, 2017). Hypoglycaemia has been associated with exhaustion and prolonged fasting (Deem, 2009, Stacy andInnis, 2017). We found higher plasma glucose levels than reported by Hamann et al. (2006), which we associated to methodological differences in the assays used, since the turtles in both studies were deemed to be healthy. Our study utilized a glucose hexokinase method, which on average, has fewer known interferences than the more commonly used glucose oxidase methods (Link et al., 2015, Dickson et al., 2019. Interferents with the glucose oxidase method could also reflect the rapidity with which the plasma was separated from the red blood cells (Kunze et al., 2020).
Enzymes-Aspartate transaminase was significantly associated with CCL, which suggests age-related changes, i.e. growth (Oh and Hustead, 2011). Alkaline phosphatase, which was influenced by mass and CCL, is an enzyme related to bone formation and osteoblast activity (van Straalen et al., 1991). This enzyme has been shown to be higher in juvenile and subadult turtles (Bolten and Bjorndal, 1992), which could have been the case in our study. Our turtles had lower plasma creatine kinase levels than those reported by March et al. (2018) in rehabilitating turtles. Elevated creatine kinase could be related to muscle catabolism (e.g. cachectic animals), capture methods and acute stress responses; however, further research is needed to confirm these hypotheses in reptiles (Anderson et al., 2013, Petrosky et al., 2015.

Nitrogenous compounds-
The results from our study demonstrated higher urea values than those reported previously (Hamann et al., 2006, Whiting et al., 2007, which may be related to a higher-protein diet (Singer, 2003 al., 2007). Creatinine concentrations measured in this study were also lower than those previously reported (Bolten and Bjorndal, 1992, Aguirre and Balazs, 2000, Hamann et al., 2006, Whiting et al., 2007, Flint et al., 2010. Analytical differences cannot be excluded either, as creatinine levels in this study were analysed using Jaffe-based chemistry (Supplementary Appendix S4), whereas other laboratories may use an enzymatic method (Delanghe and Speeckaert, 2011). However, creatinine is of minimal clinical relevance in sea turtles, and we decided to disregard this finding (Manire et al., 2002;Innis et al., 2009).
Electrolytes and minerals-Sick turtles often present elevated electrolytes (usually sodium, potassium, chloride and phosphorus), which has been linked to dehydration, renal disease, hyperaldosterism or salt gland dysfunction (Innis et al., 2009, Keller et al., 2012, Stacy and Innis, 2017. Electrolytes and minerals may be influenced by diet (in particular calcium, magnesium, sodium or phosphorus) or by the reproductive physiology of nesting females (e.g. calcium), and do not necessarily reflect pathological disorders (Raphael, 2003, Stacy and Innis, 2017, Bloodgood et al., 2019. Calcium and magnesium, for example, are associated with skeletal formation, contribute to the activation of other enzymes and can be found in high concentrations in vegetation (Bloodgood et al., 2019). The Ca:P ratio is a strong indicator for UVB deficiency, metabolic bone disease and nutritional secondary hyperparathyroidism in captive reptiles (Perrault et al., 2012;Stacy and Innis, 2017). In this study, the examined turtles exhibited a normal Ca:P ratio. Mildly to markedly inverted Ca:P ratios have also been reported in healthy turtles and might be related to life stage, diet or metabolic imbalances (Stringer et al., 2010, Kelly et al., 2015, Stacy and Innis, 2017. The turtles from this study had lower phosphorus levels than reported by Flint et al. (2010), higher magnesium levels than reported by Whiting et al. (2007) and lower sodium levels than reported by Bolten and Bjorndal (1992) (Supplementary Table S3). Our findings could be attributed to dietary differences, as no concurrent abnormalities or methodological differences in the assays used were found (Supplementary Appendix S4).

Location effects
Significant differences between Townsville region and Howick Group of Islands (P < 0.05) were found for albumin, total protein, potassium, cholesterol, triglycerides, uric acid and Ca:P ratio (Tables 3 and Supplementary Table S2). All values were still within previously reported blood values, ranges or intervals (Bolten and Bjorndal, 1992, Aguirre and Balazs, 2000, Hamann et al., 2006, Whiting et al., 2007, Arthur et al., 2008, Flint et al., 2010, March et al., 2018. Interestingly, neither mass nor CCL affected any of these parameters. We hypothesize that diet composition contributed to the analyte differences across the two sites (Whiting et al., 2007, Stacy et al., 2018, Bloodgood et al., 2019, Putillo et al., 2020. For example, Whiting et al. (2007) found that green turtles that consumed mainly seagrass had higher plasma protein levels than turtles that consumed algae. Total protein levels were higher at Howick Group of Islands than at Townsville region. It is likely that the foraging grounds at Howick Group of Islands are richer in protein sources due to higher food availability and/or nutritional content. Location differences across the same two capture sites were also reported by Flint et al. (2019), who assessed the effects of catastrophic weather events on green turtle blood analytes in, 2014-2015. Unfortunately, Flint et al. (2019 did not provide information on the statistical analyses performed, which prevented comparison of statistical methodologies.
Further, there appears to be a lack of research detailing seagrass protein content in these foraging grounds. Other factors to consider that influence total protein in sea turtle species are debilitation or malnutrition (Aguirre et al., 1995, Deem et al., 2009Innis et al., 2009;March et al., 2018) or depletion of energy during nesting (Stacy and Innis, 2017;. None of the examined turtles, however, was deemed to be unhealthy based on physical examination and on the clinical analyses. With regards to the other analytes, triglycerides and cholesterol were lowest in the Townsville region. Uric acid levels were low in comparison to other studies and were also lowest in the Townsville region. Uric acid tends to be low in healthy sea turtles (Hamann et al., 2006, Innis et al., 2009, is likely related to dietary influences (Jones and Seminoff, 2013, Barajas-Valero et al., 2021 and is sometimes found to be increased in unhealthy and/or stranded turtles (Deem et al., 2009, Innis et al., 2009, March et al., 2018. Since neither mass nor CCL influenced triglycerides, cholesterol or uric acid (Table 3), the observed location differences may be related to the nutritional composition of the foraging grounds and prey availability, rather than to dietary shifts across life stages.

Study limitations
A number of study limitations should be acknowledged. Although by comparison with other studies in green turtles, the present study's sample size was large, it was still below the recommended threshold considered to be adequate for the generation of RIs when using frequentist statistics (n < 120). This limitation was moderated by using statistical procedures (Bayesian methods) that mitigate the weakened statistical power associated with conventional frequentist statistical analysis. Despite the variation in mass and CCL in our study, two-thirds of turtles sampled were juveniles (66%, n = 64/97) and one third (34%, 33/97) were subadult and adult turtles (n = 8/97 subadult animals of undetermined sex, n = 21/97 adult females, and n = 4/97 adult males). From our sample, subadult and adult animals were mainly found in Howick Group of Islands and juvenile animals were mainly found in Townsville region (Table 2). To address this imbalance, our models also accounted for the effects of mass and CCL in the predictions.

Conclusions
Our study provides biochemical and haematological RIs for wild green turtles foraging in North Queensland, Australia, determined using Bayesian statistics that accounted for the effects of small sample sizes. Our estimated RIs fell within existing intervals and had narrower credible intervals. Location, mass and CCL effects were found for 24%, 14% and 17% of analytes, respectively. We recommend that population-specific RIs are produced with predictive statistical approaches that account for small sample sizes and for the effects of geographical location, length and mass; if they are to be used with confidence to evaluate sea turtle health. Randomized and representative sampling of the target population is essential for the determination of RIs. This is particularly important in threatened species research, which is often subject to sample size limitations. Unreliable predictions may result in false-negative or false-positive diagnoses, which can result in inadequate enforcement measures that may threaten population viability. Evidence-based sea turtle conservation and rehabilitation efforts will be enhanced by using accurate and precise RIs.

Funding
The Department of Environment and Science (Queensland Government) provided in-kind support for conducting the field trip to Howick Group of Islands. Funding for S.