Abstract

Herbivorous turtle, Chelonia mydas, inhabiting the south China Sea and breeding in Peninsular Malaysia, and Natator depressus, a carnivorous turtle inhabiting the Great Barrier Reef and breeding at Curtis Island in Queensland, Australia, differ both in diet and life history. Analysis of plasma metabolites levels and six sex steroid hormones during the peak of their nesting season in both species showed hormonal and metabolite variations. When compared with results from other studies progesterone levels were the highest whereas dihydrotestosterone was the plasma steroid hormone present at the lowest concentration in both C. mydas and N. depressus plasma. Interestingly, oestrone was observed at relatively high concentrations in comparison to oestradiol levels recorded in previous studies suggesting that it plays a significant role in nesting turtles. Also, hormonal correlations between the studied species indicate unique physiological interactions during nesting. Pearson correlation analysis showed that in N. depressus the time of oviposition was associated with elevations in both plasma corticosterone and oestrone levels. Therefore, we conclude that corticosterone and oestrone may influence nesting behaviour and physiology in N. depressus. To summarise, these two nesting turtle species can be distinguished based on the hormonal profile of oestrone, progesterone, and testosterone using discriminant analysis.

1. Introduction

Hormones have various functions across animal taxa. In birds, for example, testosterone is linked to mating success, while corticosterone has been shown to support energetically demanding processes [1]. In vertebrates glucocorticoid hormones can assist in the modulation of life-history events, survival probabilities, and fecundity [2]. Similarly steroid hormones play an essential role in the reproductive behaviour and sexual development in marine turtles [3]. Sex determination in sea turtles is regulated by incubation temperatures and by exogenous steroid hormones [4, 5]. In addition to determining sex, steroids have also been causally related to various physiological and behavioural actions in sea turtles. For example, corticosterone and testosterone have been associated with hepatic regulation and ovarian function during nesting activities in Caretta caretta [68]. Corticosterone concentration decreases when a turtle is subjected to stress or drops to low levels when a female turtle is breeding [911]. Nesting Lepidochelys olivacea during arribadas (mass nesting behaviour) maintain low corticosterone level as a mechanism of controlling stress [12]. However, a low adrenocortical response may assist reproduction in sea turtles [7, 8]. Fatty acids and glucose, derived from fat and glycogen, respectively, are characterised as indicators of the required energy for nesting [7, 8, 13]. In Chelonia mydas corticosterone increases glycerol release from adipose tissue suggesting that this hormone induces the mobilization of triglycerides accumulated in the adipocytes [14]. Consequently, it causes the release of free fatty acid into the blood during nesting [7, 8]. Because sea turtles are largely aphagic during reproductive and nesting activities [15, 16], corticosterone induces free fatty acids release from adipose tissue to compensate for the lack of food and the energy demands during courtship and nesting activities in sea turtles, indicating that corticosterone is an essential steroid hormone for energetic homeostasis of sea turtles during reproduction [13].

Thyroxine is another hormone responsible for regulating energy availability in animals. High plasma thyroxine concentrations increase the metabolic rate in animal tissues [17], reflecting thyroxine levels as indicators of metabolic activity. In the rainbow trout, Salmo gairdneri, thyroid activity increases thyroxine secretion, which affects body metabolism and it is associated with yolk deposition, mating, or ovulation [18].

Taking into consideration the important roles of steroids and metabolites (i.e., glucose and triglycerides) during the reproductive cycles, we hypothesised that despite similarities during nesting, the steroid hormone and metabolite profiles in sea turtles is specific for each population and species. Hence, two turtle species (C. mydas versus Natator depressus) with different diet (herbivorous versus carnivorous) and distribution (Peninsular, Malaysia, versus Queensland, Australia) were likely to have unique hormone and metabolite profiles. We studied also the role of oestrone in nesting C. mydas and N. depressus and measured the plasma concentration profile of thyroxine, a hormone that has not been investigated before in nesting turtles. Blood metabolites (i.e., glucose and triglycerides) in conjunction with steroid hormones assist turtles to conserve energy during breeding and nesting season (e.g., C. mydas [7, 8]). Hence, understanding metabolites role and their interactions with hormones during nesting activities in sea turtles will provide a more comprehensive knowledge on sea turtle biology. Consequently, we also examined the interactions among hormone and metabolite levels in both species as well as the influence of nesting history, time of oviposition, and size (curved carapace length measurements) to hormone and metabolite levels in N. depressus. Similarly, we studied a potential link between egg size (egg dimensions) and hormone and/or metabolite levels in N. depressus. Last but not least, we used discriminant analysis to separate the two species based on their hormones and metabolite levels. We argued that this method would provide an additional insight on the nesting biology of the studied species.

2. Materials and Methods

2.1. Field Sites

This study investigated hormonal and metabolite patterns collected at the peak of the respective nesting seasons for C. mydas (Ibrahim, personal observations) and N. depressus [19, 20]. Chelonia mydas sampling was conducted between June and September 2004 at Ma’Daerah beach (4°34′15.50′N, 103°27′51.93′E), Kuala Terengganu State, in Peninsular Malaysia. Natator depressus sampling was conducted between November and December 2006 at Curtis Island (23°45′S, 151°17′E), Gladstone, in central Queensland, Australia.

2.2. Blood Collection

Blood samples (2 mL) were collected from the dorsal cervical sinus from 14 nesting C. mydas and 20 nesting N. depressus using a 21 G needle and a 5 mL syringe immediately after oviposition as previously been described [21]. Blood samples were initially stored in heparinized tubes on ice (2 hr) prior to being centrifuged at 6,000 rpm for 10 minutes. For C. mydas in Peninsular Malaysia, plasma was decanted and frozen in liquid nitrogen and then transferred on dry ice to Queensland and stored at −25°C until analysed. The N. depressus plasma samples were kept for <2 hr on ice and then stored at −25°C until analysed.

2.3. Biometric Measurements in Natator depressus

Midline curved carapace length (CCL) was measured with a flexible tape by two independent turtle researchers (the author MI always conducted one of the measurements) as previously been described [22]. Egg size was mean product of the maximum and minimum egg diameter of three eggs from each nesting turtle. Time of oviposition (TO) was recorded as a continuous variable from 16:00 to 24:00. Nesting history was acquired from the Queensland Turtle Research Project database. The year each individual turtle was reported nesting for the first time was grouped in three categories: first time nesting turtles in 2006; younger remigrant turtles, which began breeding during 1998–2002; and older remigrant turtles, which began breeding during 1993–1997.

2.4. Steroid Hormone Analysis

Plasma samples from C. mydas and N. depressus were analysed directly using specific enzyme linked immunoassays (ELISAs) for testosterone (T), dihydrotestosterone (DHT), progesterone (P), thyroxine (T4), oestrone (E1) (dbc Diagnostics, Canada), and corticosterone (B) (Abacus Diagnostics, Australia). Oestrone and B ELISA assays have been validated previously for use in marine turtle and crocodilian samples, respectively [23, 24]. Still, in this study, all ELISA assays were validated by using pooled samples from nesting C. mydas () and N. depressus (–4) over a total period of four months and with estimates comparable to those provided in the kits by the manufacturers. Pooled samples in duplicates were used for intra-assay, inter-assay variation, and minimum assay sensitivity (Table 1). For E1, T, DHT, T4, P, and B ELISA kits the published intra-assay precision values, based upon three human serum samples (except B, which was based upon two human serum samples) were assayed ten times (Table 1). Inter-assay variations of T4, P, and DHT ELISA kits were not possible because only one kit of these steroids was used. However, because no major intra-assay variations were observed in these ELISA kits there was a confidence of the kit validation in marine turtle plasma. The percent coefficient of variation (% CV) was determined as equal to SD/mean 100(%).

2.5. Enzyme Linked Immunoassays

Detailed information on ELISA assays is provided in Ikonomopoulou [25]. In brief, all ELISA kits (E1, T, DHT, T4, P, and B) used an antisteroid antibody, a steroid-biotin conjugate solution, and a steroid-avidin-horse radish peroxide conjugate. Appropriate calibration standards, positive control, and plasma samples were added to the plate and incubated for 1-2 hrs, depending on the steroid and the manufacturer’s protocol. After the incubation period, the plate was washed and incubated with TMB (3,3′,5,5′-tetramethylbenzidine) for 15 minutes on a Benchmark Microwell Plate Reader. The steroid concentrations (E1, T, DHT, T4, P, and B) of the unknown samples were interpolated (using an Excel program for PC) based on the steroid calibrations curves.

All samples were within the range of the calibration curves enabling the steroid protocols given by the manufacturers to be employed without modifications (e.g., sample dilutions). The ELISA kits were tested for precision, sensitivity, accuracy, and steroid specificity.

2.6. Precision (CV%)

Calibration curves for each steroid ELISA kit for C. mydas and N. depressus were completed by measuring the absorbance of six calibrators as described in Ikonomopoulou [25]. Two calibrators (for each measured steroid) were run as unknown samples to test intra-assay variation, determining accuracy in each ELISA kit.

2.7. Sensitivity

The sensitivity of the ELISA kits was determined by measuring the mean absorbance (OD) of the lowest standard concentrations as described in Ikonomopoulou [25] and comparing this with the absorbance of samples containing no hormone. The sensitivity was reported to be E1: 10.0 pg/mL; T: 0.022 ng/mL; (dbc, Inc.); DHT: 6.0 pg/mL; T4: 0.6 μg/100 μL; P: 0.1 ng/mL; and B: 0.17 ng/mL.

2.8. Accuracy

ELISA accuracy for marine turtles (C. mydas and N. depressus) was determined using a positive control sample containing the appropriate steroid (E1, T, DHT, P, B, and T4) in human serum. The expected concentration of the control was divided by the measured concentration and multiplied by 100 to obtain a percentage. The determined accuracy in each ELISA kit for marine turtles was as follows: E1 (89.3%); T (95.2%); DHT (90.70%); T4 (100.61%); P (94.72%); and B (100.85%). All the % accuracy values for the steroids were within the manufacturers reported levels E1 (88–125%); T (80.5–110.1%); DHT (88–97.9%); T4 (93.9–116.4%); P (78–124%); and B (94–116%).

2.9. Specificity (% CV)

The ELISA kits were highly specific to the steroids they measured. All the ELISA kits (T, E1, P, T4, and DHT) were tested for specificity in the commercial laboratory (unpublished data, dpc, Inc) and B ELISA kit were tested in the laboratories of Abacus Diagnostics (unpublished data, Abacus Diagnostics). Information on all steroid ELISA kits that crossed reactivity with other steroids is given in Ikonomopoulou [25].

2.10. Metabolite Analysis

Glucose and triglyceride in the plasma of nesting C. mydas and N. depressus females were measured using an automated Olympus AU400 analyser and Olympus reagents (Olympus Diagnostica, Clare Island) by The University of Queensland, Department of Veterinary Pathology, Unit Diagnostic Service. Glucose and triglyceride assays have been previously validated for marine turtles by the Department of Veterinary Pathology Unit (Table 2). Glucose analysis was performed according to HKG6P-DH normal application with NIST SRM 965 traceability (Veterinary Pathology Unit, internal methods). The triglyceride analysis was done according to GPO-POD normal application with Isotope dilution mass spectrometry reference method traceability (Veterinary Pathology Unit, internal methods).

2.11. Hormone and Metabolites Statistical Analysis

All values are reported as mean hormone and metabolite levels (ng/mL) ± standard error (SE). R statistical programming language was used for the analysis [26]. All traits were analysed using a one way ANOVA, using “lm” function of the “stats” package [26], with the metabolite/hormone level as the outcome variable and turtle species as the explanatory variables. Residual analysis was performed to assess potential violations of the assumptions of ANOVA. No such violations were detected (results not shown). As an extra safeguard, and for assessing the robustness of such estimates, a bootstrap procedure was used. This procedure gave similar results to the simple ANOVA and thus only the ANOVA results are presented. Least square means were estimated using the Ismeans R package [27] and are presented here. DHT was not included in the statistical analysis since the values were found to be below the assays sensitivity level. In addition, we explored the associations between metabolites and hormones in the two species. We performed Pearson’s correlation analysis, using “Hmisc” package [28], within each species testing two different hypothesis. (i) Are the correlations in each species in our sample different from zero? (ii) Are the two sets of correlations different from each other (i.e., between species comparison)?

We also assessed whether a combination of the above hormones and metabolites could be used to differentiate the two species. This analysis would answer a question such as the following. “Given a blood sample from a turtle of these species, could we use a combination of hormones/metabolites to identify from which species the sample came?” This question was approached using linear discriminant analysis, in R using the MASS package [29].

2.12. The Association of Hormones and Metabolites with Natator depressus Biometric Measurements

Data on nesting history, turtle size (CCL), egg size, and the time of oviposition were also collected. These biometric measurements could be associated with the hormone and metabolite profile of turtles and thus a Pearson correlation analysis was performed. It should be noted that the analysis of time was performed both on the original scale and with decimal time. The results were identical and thus only these for decimal time are presented. It should be noted that in the statistical analysis for the evaluation of nesting history on hormones/metabolites profile only two categories were used because there was only one turtle reported as first time nesting turtle in 2006.

3. Results

3.1. Species-Specific Steroid Hormonal Profile in Sea Turtles

Steroid hormonal profiling from plasma samples obtained from C. mydas and N. depressus while nesting showed that among all of the analysed steroid hormones the mean levels (±SE) of P were the highest ( versus  ng/mL) in both turtle species. Progesterone was followed by T4 ( versus  ng/mL), T ( versus  ng/mL), E1 ( versus  ng/mL), and B (: below assay sensitivity versus  ng/mL) and DHT values for both species were below assay sensitivity (~0.002 ± 0.001). Hence, DHT data obtained from the statistical analysis should be interpreted with caution.

As it can be seen in Table 3, statistically significant differences were observed for all hormones except thyroxine (ANOVA: F1,32 = 0.23). Chelonia mydas had higher levels of progesterone but lower levels for all other hormones and metabolites.

The within species correlations are shown in Tables 4 and 5. The correlation between progesterone and testosterone was significantly different from zero for both species (Pearson correlation: 0.80 and 0.63 for C. mydas and N. depressus, resp.), whereas the correlations between corticosterone-progesterone, corticosterone-oestrone, and testosterone-thyroxine were statistically significant from zero for N. depressus (Pearson correlation: 0.47, 0.64, and 0.51, resp.). Interestingly the thyroxine-testosterone correlation closely approached statistical significance for difference between the two species (Pearson correlations of −0.19 and 0.51, resp.) (; Table 6).

The linear discriminant analysis using a stepwise procedure resulted in 3 variables being selected: oestrone, progesterone, and testosterone. The optimal function for discrimination was D = (−4.35 × oestrone) + (0.0975 × progesterone) − (2.915 × testosterone). This function was able to separate 100% the two species with negative values being for N. depressus and positive values being for C. mydas (Figure 1).

Average inter- and intra-assay coefficients of variation based on the absorbance of standards for E1, T, DHT, B, and T4 and the sensitivity for each assay are shown in Table 1.

Metabolites, blood glucose, and triglycerides are important sources of energy during responses to physiological and stress stimuli including nesting. No significant differences were observed in blood glucose and triglyceride concentrations between C. mydas and N. depressus (ANOVA: F1,32 = 2.23; versus  mmol/L; F1,32 = 3.00; versus  mmol/L, ).

3.2. Correlations among Hormones and Metabolites

In C. mydas, a positive correlation was observed between T and P (Pearson correlation: , ).

On the other hand in N. depressus, T was positively correlated with P (, ) but also with T4 (, ). In addition, B was positively correlated with both P (, ) and E1 (, ).

3.3. Correlation of Nesting History, Carapace Length, Time of Oviposition, and Egg Size with Metabolites and Hormones in Natator depressus

We also analysed whether nesting history, carapace length, time of oviposition, and egg size might be correlated with the hormonal or metabolites plasma profile. Due to technical limitations we recorded measurements only for N. depressus. The results for this Pearson correlation analysis are given in Table 7. Significantly, negative correlations were observed between corticosterone and oestrone with time of oviposition ( and , , resp.), suggesting that the time of oviposition might influence corticosterone and oestrone levels in N. depressus. No statistically significant correlations were estimated for any other traits.

4. Discussion

Marine turtle species are globally distributed and hence exposed to varying environmental conditions. In this study, we observed the hormone and metabolite profiles of C. mydas and N. depressus. The reasons for the noted hormone and metabolite variations remain unknown. Critical information was missing on the nesting status (i.e., when a turtle was nested for the first time and the number of clutches laid in a nesting season) of the studied individuals at the time of sampling. In general, C. mydas nesting in Heron Island lays 5+ clutches [30] while N. depressus in Curtis Island lays approximately 3 clutches per season [31]. Knowing the deposited number of clutches is important because it could influence hormone levels in nesting turtles [6]. For example, hormones such as P, T, E2, and B are at higher levels at the beginning than at the end of a nesting season [7, 8, 32]. Nevertheless, C. mydas sampling occurred during a 3-month period; hence, it was highly unlikely that all of the individuals were at the same nesting status. In addition, hormone level variations as observed from the given SE within individuals of both species were minimum.

This study revealed a different hormone and metabolite profile for C. mydas and N. depressus. These discrepancies might have been influenced by diet and/or geographical location parameters. Chelonia mydas, as a herbivorous species, is predominantly feeding on sea grass, seaweeds, and algae [33] and, depending on the foraging site, its diet can vary, characterising this species as an opportunistic forager [34]. On the contrary, N. depressus is carnivorous and is feeding primarily on soft-bodied invertebrates, including jellyfish, sea cucumbers, sea pens, and soft corals [35]. Indeed, progesterone levels have been reported to be high (C. mydas and C. caretta [36]) or low (C. mydas and C. caretta [37] and Dermochelys coriacea [38]) or to decline after a clutch deposition (C. caretta [6] and Lepidochelys kempii [39]) in different nesting turtle populations worldwide. We observed that progesterone, as the major female sex steroid hormone during pregnancy, had the highest plasma concentration whereas DHT, a major testosterone derivative that is mainly observed in males, showed the lowest plasma concentration of the six analysed steroid hormones in both turtle species. In a previous study it has been shown that progesterone binds with low affinity and capacity in the plasma of nesting C. mydas [40], explaining the high progesterone levels found for C. mydas in this study.

In contrast to progesterone, both testosterone and oestradiol bind with high affinity and low capacity to sex steroid binding protein (SSBP), indicating that most of this hormone remains bound to a high affinity carrier in the body of nesting C. mydas [40]. Prior to migration to nesting sites, female turtles increase their T concentrations [32]. Whittier et al. [6] previously reported a correlation between testosterone and corticosterone in nesting C. caretta but such a correlation was not observed either in N. depressus or in C. mydas (C. mydas [7, 8]; this study). All these observations together suggest a pivotal role of testosterone that is specific to each female turtle species enabling completion of nesting. In addition, the parallel analysis of steroid profiles in C. mydas and N. depressus in this study for the first time supports that steroid plasma hormone concentrations are species and possibly genetic populations specific and they might be influenced by various factors.

Nevertheless, the link between oestrone levels and time of oviposition in N. depressus also supports the hypothesis that hormone profiles, as well as their formed correlations with various factors (e.g., time of oviposition), are interactive for the successful nesting activities. We postulate that oestrone, together with oestradiol, might play an essential function during nesting in turtles. Despite the important role played by androgenic hormones, oestrogenic hormones are the most vital hormones during nesting in sea turtles. Previous studies have been focused predominantly on oestradiol measurements in nesting turtles. Oestradiol 17β remains low or is at undetectable levels in nesting C. mydas and C. caretta turtles [68, 36]. However, other researchers have previously suggested that there must be another major oestrogenic hormone besides oestradiol in nesting turtles (e.g., [6, 36]). Indeed, oestrone plasma levels were higher than those of oestradiol in nesting C. mydas on Raine Island [23], supporting our conclusion about the major role of oestrone in nesting turtles. Nevertheless, we should emphasise that due to insufficient plasma quantity we were unable to measure oestradiol and verify whether its concentration was indeed lower than oestrone as hypothesised. Finally, the fact that N. depressus produces the largest eggs among turtle species [41] and that oestrone levels are influenced by the oviposition time in this species supports the need for further studies in a larger range of sea turtle species and populations that will provide conclusive insights into the role of this hormone in the aforementioned nesting parameters.

In respect to corticosterone levels, as found for the levels of oestrone, our data suggest that they were influenced by the oviposition time in N. depressus. Corticosterone levels were low in both nesting species. However, because corticosterone concentrations in C. mydas were below assay sensitivity levels, caution should be exercised when interpreting these results. Having said this, the low corticosterone concentrations in C. mydas have been characterised as an adaptive mechanism for optimising the reproductive success during stressful situations [9, 10]. Even at such low concentrations as these measured in this study, corticosterone may regulate and supply the necessary energy for the nesting cycle for C. mydas as it has previously been suggested [7, 8].

It should also be highlighted that in this study we identified the two turtle species 100% correctly, using discriminant analysis, based on three hormone levels (i.e., oestrone, progesterone, and testosterone). We suggest that during nesting these hormones are very distinguishable in C. mydas and N. depressus.

5. Conclusions

In conclusion, this study provides some fundamental baseline insights into the hormonal and metabolite patterns of two nesting turtle species with different diet preferences, geographical locations, and evolution pathways. We postulate that the hormone and metabolite differences as well as their interactions were meaningful and necessary for the successful completion of nesting activities, reflecting the species evolved physiological uniqueness and their complexity during nesting. We propose that oestrone is a major oestrogenic hormone in nesting turtles modulating, along with corticosterone at least in N. depressus, time of oviposition. Understanding the hormonal patterns of these species could lead us to comprehend better the effect of turtle species adaptation to different environments.

Finally, we were able to successfully distinguish nesting C. mydas and N. depressus based on the hormone profile of oestrone, progesterone, and testosterone, leading to potential new investigation pathways for identifying turtle species at particular life stages.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by an Australian Research Council Linkage grant (LP0455513) to JMW and HO (Queensland Health Scientific Services; QHSS). All experiments were performed in accordance with approval by the Animal Ethics Committee of the University of Queensland (nos. ANAT/238/06/UQ-ARC/Earthwatch/URG and ANAT/352/07/UQ-ARC/Earthwatch/URG) and the Queensland Environmental Protection Agency (QPWS permit numbers WISP03952406). The authors thank the Malaysian Fisheries Department Team, the rangers in K.T. They thank also the Environmental Protection Agency Team, the project staff, and the volunteers for the 2006 Curtis Island nesting census project.