Sugar-coated survival: N-glycosylation as a unique bearded dragon venom resistance trait within Australian agamid lizards.

In the ongoing evolutionary arms race between predators and prey, adaptive innovations often trigger a reciprocal response. For instance, the emergence of α-neurotoxins in snake venom has driven prey species targeted by these snakes to evolve sophisticated defense mechanisms. This study zeroes in on the particular motifs within the orthosteric sites of post-synaptic nicotinic acetylcholine receptors (nAChR) that confer resistance to α-neurotoxins, often through structural alterations of nAChR. This research examined Australian agamid lizards, a primary prey group for Australian elapid snakes, which are subject to predatory selection pressures. We previously showed that Pogona vitticeps (Central bearded dragon) was resistant to α-neurotoxic snake venoms through a steric hindrance form resistance evolving within the nAChR orthosteric, specifically through the 187-189NVT motif resulting in the presence of N-glycosylation, with the branching carbohydrate chains impeding the binding by the neurotoxins. This adaptive trait is thought to be a compensatory mechanism for the lizard's limited escape capabilities. Despite the significance of this novel adaptation, the prevalence and evolutionary roots of such venom resistance in Australian agamids have not been thoroughly investigated. To fill this knowledge gap, we undertook a comprehensive sequencing analysis of the nAChR ligand-binding domain across the full taxonomical diversity of Australian agamid species. Our findings reveal that the N-glycosylation resistance mechanism is a trait unique to the Pogona genus and absent in other Australian agamids. This aligns with Pogona's distinctive morphology, which likely increases vulnerability to neurotoxic elapid snakes, thereby increasing selective pressures for resistance. In contrast, biolayer interferometry experiments with death adder (Acanthophis species) venoms did not indicate any resistance-related binding patterns in other agamids, suggesting a lack of similar resistance adaptations, consistent with these lineages either being fast-moving, covered with large defensive spines, or being arboreal. This research not only uncovers a novel α-neurotoxin resistance mechanism in Australian agamids but also highlights the complex dynamics of the predator-prey chemical arms race. It provides a deeper understanding of how evolutionary pressures shape the interactions between venomous snakes and their prey.

Systematic assessment of the brown tree frog (Anura: Pelodryadidae: Litoria ewingii) reveals two endemic species in South Australia.

The brown tree frog (Litoria ewingii) is a relatively widespread, commonly encountered pelodryadid frog from south-eastern Australia, known for its characteristic whistling call. The distribution of Litoria ewingii spans over more than 350,000 km2, encompassing a range of moist temperate habitats, and is fragmented by well-known biogeographic barriers. A preliminary analysis of mitochondrial DNA sequences revealed evidence for deep phylogenetic structure between some of these fragmented populations. In this study, we sought to re-evaluate the systematics and taxonomy of Litoria ewingii sensu lato by analysing variation in nuclear and mitochondrial DNA, adult morphology and male advertisement calls throughout the species range. Our analyses reveal two additional, deeply divergent and allopatric lineages in South Australia. We herein re-describe Litoria ewingii from Tasmania, southern New South Wales, Victoria and south-eastern South Australia, resurrect the name Litoria calliscelis for a species occurring in the Mount Lofty Ranges and Fleurieu Peninsula in South Australia, and describe a new species, Litoria sibilus sp. nov., endemic to Kangaroo Island.

A Russian Doll of Resistance: Nested Gains and Losses of Venom Immunity in Varanid Lizards.

The interplay between predator and prey has catalyzed the evolution of venom systems, with predators honing their venoms in response to the evolving resistance of prey. A previous study showed that the African varanid species Varanus exanthematicus has heightened resistance to snake venoms compared to the Australian species V. giganteus, V. komodoensis, and V. mertensi, likely due to increased predation by sympatric venomous snakes on V. exanthematicus. To understand venom resistance among varanid lizards, we analyzed the receptor site targeted by venoms in 27 varanid lizards, including 25 Australian varanids. The results indicate an active evolutionary arms race between Australian varanid lizards and sympatric neurotoxic elapid snakes. Large species preying on venomous snakes exhibit inherited neurotoxin resistance, a trait potentially linked to their predatory habits. Consistent with the 'use it or lose it' aspect of venom resistance, this trait was secondarily reduced in two lineages that had convergently evolved gigantism (V. giganteus and the V. komodoensis/V. varius clade), suggestive of increased predatory success accompanying extreme size and also increased mechanical protection against envenomation due to larger scale osteoderms. Resistance was completely lost in the mangrove monitor V. indicus, consistent with venomous snakes not being common in their arboreal and aquatic niche. Conversely, dwarf varanids demonstrate a secondary loss at the base of the clade, with resistance subsequently re-evolving in the burrowing V. acanthurus/V. storri clade, suggesting an ongoing battle with neurotoxic predators. Intriguingly, within the V. acanthurus/V. storri clade, resistance was lost again in V. kingorum, which is morphologically and ecologically distinct from other members of this clade. Resistance was also re-evolved in V. glebopalma which is terrestrial in contrast to the arboreal/cliff dwelling niches occupied by the other members of its clade (V. glebopalma, V. mitchelli, V. scalaris, V. tristis). This 'Russian doll' pattern of venom resistance underscores the dynamic interaction between dwarf varanids and Australian neurotoxic elapid snakes. Our research, which included testing Acanthophis (death adder) venoms against varanid receptors as models for alpha-neurotoxic interactions, uncovered a fascinating instance of the Red Queen Hypothesis: some death adders have developed more potent toxins specifically targeting resistant varanids, a clear sign of the relentless predator-prey arms race. These results offer new insight into the complex dynamics of venom resistance and highlight the intricate ecological interactions that shape the natural world.

The Eastern Bandy Bandy Vermicella annulata, expresses high abundance of SVMP, CRiSP and Kunitz protein families in its venom proteome.

The Australian elapid snake radiation (Hydrophiinae) has evolved in the absence of competition from other advanced snakes. This has resulted in ecological specialisation in Australian elapids and the potential for venom proteomes divergent to other elapids. We characterised the venom of the Australian elapid Vermicella annulata (eastern bandy bandy). The venom was analysed using a two-dimensional fractionation process consisting of reverse-phase high-performance liquid chromatography then sodium dodecyl sulphate polyacrylamide gel electrophoresis, followed by bottom-up proteomics. Resulting peptides were matched to a species-specific transcriptome and 87% of the venom was characterised. We identified 11 toxins in the venom from six families: snake venom metalloproteinases (SVMP; 24.2%; two toxins) that are class P-III SVMPs containing a disintegrin-like domain, three-finger toxins (3FTx; 21.6%; five toxins), kunitz peptides (KUN; 19.5%; one toxin), cysteine-rich secretory proteins (CRiSP; 18%; one toxin), and phospholipase A2 (PLA2; 4%; two toxins). The venom had low toxin diversity with five protein families having one or two toxins, except for 3FTx with five different toxins. V. annulata expresses an unusual venom proteome, with high abundances of CRiSP, KUN and SVMP, which are not normally highly expressed in elapid venoms. This unusual venom composition could be an adaptation to its specialised diet. BIOLOGICAL SIGNIFICANCE: Although the Australian elapid radiation represents the most extensive speciation event of elapids on any continent, with 100 terrestrial species, the venom composition of these snakes has rarely been investigated, with only five species currently characterised. Here we provide the venom proteome of a sixth species, Vermicella annulata. The venom of this species could be particularly informative from an evolutionary perspective, as it is an extreme dietary specialist, only preying on blind snakes (Typhlopidae). We show that V. annulata expresses a highly unusual venom for an elapid, due to the high abundance of the protein families SVMP, CRiSP, and KUN, which together make up 61% of the venom. When averaged across all species, a typical elapid venom is 82% PLA2 and 3FTx. This is the second recorded instance of an Australian elapid having evolved highly divergent venom expression.

The Unusual Metalloprotease-Rich Venom Proteome of the Australian Elapid Snake Hoplocephalus stephensii.

The Australasian region is home to the most diverse elapid snake radiation on the planet (Hydrophiinae). Many of these snakes have evolved into unique ecomorphs compared to elapids on other continents; however, their venom compositions are poorly known. The Australian elapid Hoplocephalus stephensii (Stephen's banded snake) is an arboreal snake with a unique morphology. Human envenoming results in venom-induced consumption coagulopathy, without neurotoxicity. Using transcriptomics and a multi-step fractionation method involving reverse-phase high-performance liquid chromatography, sodium dodecyl sulfate polyacrylamide gel electrophoresis and bottom-up proteomics, we characterized the venom proteome of H. stephensii. 92% of the total protein component of the venom by weight was characterized, and included all dominant protein families and 4 secondary protein families. Eighteen toxins made up 76% of the venom, four previously characterized and 14 new toxins. The four dominant protein families made up 77% of the venom, including snake venom metalloprotease (SVMP; 36.7%; three identified toxins), phospholipase A2 (PLA2; 24.0%; five identified toxins), three-finger toxin (3FTx; 10.2%; two toxins) and snake venom serine protease (SVSP; 5.9%; one toxin; Hopsarin). Secondary protein families included L-amino acid oxidase (LAAO; 10.8%; one toxin), natriuretic peptide (NP; 0.8%; two toxins), cysteine-rich secretory protein (CRiSP; 1.7%; two toxins), c-type lectin (CTL; 1.1%; one toxin), and one minor protein family, nerve growth factor (NGF; 0.8%; one toxin). The venom composition of H. stephensii differs to other elapids, with a large proportion of SVMP and LAAO, and a relatively small amount of 3FTx. H. stephensii venom appeared to have less toxin diversity than other elapids, with only 18 toxins making up three-quarters of the venom.

Lineage diversity within a widespread endemic Australian skink to better inform conservation in response to regional-scale disturbance.

Much attention is paid in conservation planning to the concept of a species, to ensure comparability across studies and regions when classifying taxa against criteria of endangerment and setting priorities for action. However, various jurisdictions now allow taxonomic ranks below the level of species and nontaxonomic intraspecific divisions to be factored into conservation planning-subspecies, key populations, evolutionarily significant units, or designatable units. Understanding patterns of genetic diversity and its distribution across the landscape is a key component in the identification of species boundaries and determination of substantial geographic structure within species. A total of 12,532 reliable polymorphic SNP loci were generated from 63 populations (286 individuals) covering the distribution of the Australian eastern three-lined skink, Bassiana duperreyi, to assess genetic population structure in the form of diagnosable lineages and their distribution across the landscape, with particular reference to the recent catastrophic bushfires of eastern Australia. Five well-supported diagnosable operational taxonomic units (OTUs) existed within B. duperreyi. Low levels of divergence of B. duperreyi between mainland Australia and Tasmania (no fixed allelic differences) support the notion of episodic exchange of alleles across Bass Strait (ca 60 m, 25 Kya) during periods of low sea level during the Upper Pleistocene rather than the much longer period of isolation (1.7 My) indicated by earlier studies using mitochondrial sequence variation. Our study provides foundational work for the detailed taxonomic re-evaluation of this species complex and the need for biodiversity assessment to include an examination of cryptic species and/or cryptic diversity below the level of species. Such information on lineage diversity within species and its distribution in the context of disturbance at a regional scale can be factored into conservation planning regardless of whether a decision is made to formally diagnose new species taxonomically and nomenclaturally.

A return-on-investment approach for prioritization of rigorous taxonomic research needed to inform responses to the biodiversity crisis.

Global biodiversity loss is a profound consequence of human activity. Disturbingly, biodiversity loss is greater than realized because of the unknown number of undocumented species. Conservation fundamentally relies on taxonomic recognition of species, but only a fraction of biodiversity is described. Here, we provide a new quantitative approach for prioritizing rigorous taxonomic research for conservation. We implement this approach in a highly diverse vertebrate group-Australian lizards and snakes. Of 870 species assessed, we identified 282 (32.4%) with taxonomic uncertainty, of which 17.6% likely comprise undescribed species of conservation concern. We identify 24 species in need of immediate taxonomic attention to facilitate conservation. Using a broadly applicable return-on-investment framework, we demonstrate the importance of prioritizing the fundamental work of identifying species before they are lost.

Activity of two key toxin groups in Australian elapid venoms show a strong correlation to phylogeny but not to diet.

BACKGROUND: The relative influence of diet and phylogeny on snake venom activity is a poorly understood aspect of snake venom evolution. We measured the activity of two enzyme toxin groups - phospholipase A2 (PLA2), and L-amino acid oxidase (LAAO) - in the venom of 39 species of Australian elapids (40% of terrestrial species diversity) and used linear parsimony and BayesTraits to investigate any correlation between enzyme activity and phylogeny or diet. RESULTS: PLA2 activity ranged from 0 to 481 nmol/min/mg of venom, and LAAO activity ranged from 0 to 351 nmol/min/mg. Phylogenetic comparative methods, implemented in BayesTraits showed that enzyme activity was strongly correlated with phylogeny, more so for LAAO activity. For example, LAAO activity was absent in both the Vermicella and Pseudonaja/Oxyuranus clade, supporting previously proposed relationships among these disparate taxa. There was no association between broad dietary categories and either enzyme activity. There was strong evidence for faster initial rates of change over evolutionary time for LAAO (delta parameter mean 0.2), but no such pattern in PLA2 (delta parameter mean 0.64). There were some exceptions to the phylogenetic patterns of enzyme activity: different PLA2 activity in the ecologically similar sister-species Denisonia devisi and D. maculata; large inter-specific differences in PLA2 activity in Hoplocephalus and Austrelaps. CONCLUSIONS: We have shown that phylogeny is a stronger influence on venom enzyme activity than diet for two of the four major enzyme families present in snake venoms. PLA2 and LAAO activities had contrasting evolutionary dynamics with the higher delta value for PLA2 Some species/individuals lacked activity in one protein family suggesting that the loss of single protein family may not incur a significant fitness cost.

An Integrative Approach Using Phylogenomics and High-Resolution X-Ray Computed Tomography for Species Delimitation in Cryptic Taxa.

Morphologically cryptic taxa have proved to be a long-standing challenge for taxonomists. Lineages that show strong genomic structuring across the landscape but are phenotypically similar pose a conundrum, with traditional morphological analyses of these cryptic lineages struggling to keep up with species delimitation advances. Micro X-ray computed tomography (CT) combined with geometric morphometric analyses provides a promising avenue for identification of morphologically cryptic taxa, given its ability to detect subtle differences in anatomical structures. However, this approach has yet to be used in combination with genomic data in a comparative analytical framework to distinguish cryptic taxa. We present an integrative approach incorporating genomic and geometric morphometric evidence to assess the species delimitation of grassland earless dragons (Tympanocryptis spp.) in north-eastern Australia. Using mitochondrial and nuclear genes (ND2 and RAG1, respectively), along with $>$8500 SNPs (nuclear single nucleotide polymorphisms), we assess the evolutionary independence of target lineages and several closely related species. We then integrate phylogenomic data with osteological cranial variation between lineages using landmark-based analyses of three-dimensional CT models. High levels of genomic differentiation between the three target lineages were uncovered, also supported by significant osteological differences. By incorporating multiple lines of evidence, we provide strong support for three undescribed cryptic lineages of Tympanocryptis in north-eastern Australia that warrant taxonomic review. Our approach demonstrates the successful application of CT with integrative taxonomic approaches for cryptic species delimitation, which is broadly applicable across vertebrates containing morphologically similar yet genetically distinct lineages. Additionally, we provide a review of recent integrative taxonomic approaches for cryptic species delimitation and an assessment of how our approach can value-add to taxonomic research.

Taxonomy and conservation of grassland earless dragons: new species and an assessment of the first possible extinction of a reptile on mainland Australia.

Taxonomic research is of fundamental importance in conservation management of threatened species, providing an understanding of species diversity on which management plans are based. The grassland earless dragon lizards (Agamidae: Tympanocryptis) of southeastern Australia have long been of conservation concern but there have been ongoing taxonomic uncertainties. We provide a comprehensive taxonomic review of this group, integrating multiple lines of evidence, including phylogeography (mtDNA), phylogenomics (SNPs), external morphology and micro X-ray CT scans. Based on these data we assign the lectotype of T. lineata to the Canberra region, restrict the distribution of T. pinguicolla to Victoria and name two new species: T. osbornei sp. nov. (Cooma) and T. mccartneyi sp. nov. (Bathurst). Our results have significant conservation implications. Of particular concern is T. pinguicolla, with the last confident sighting in 1969, raising the possibility of the first extinction of a reptile on mainland Australia. However, our results are equivocal as to whether T. pinguicolla is extant or extinct, emphasizing the immediate imperative for continued surveys to locate any remaining populations of T. pinguicolla. We also highlight the need for a full revision of conservation management plans for all the grassland earless dragons.

A citizen science model for implementing statewide educational DNA barcoding.

Our aim was to develop a widely available educational program in which students conducted authentic research that met the expectations of both the scientific and educational communities. This paper describes the development and implementation of a citizen science project based on DNA barcoding of reptile specimens obtained from the Museums Victoria frozen tissue collection. The student program was run by the Gene Technology Access Centre (GTAC) and was delivered as a "one day plus one lesson" format incorporating a one-day wet laboratory workshop followed by a single lesson at school utilising online bioinformatics tools. The project leveraged the complementary resources and expertise of the research and educational partners to generate robust scientific data that could be analysed with confidence, meet the requirements of the Victorian state education curriculum, and provide participating students with an enhanced learning experience. During two 1-week stints in 2013 and 2014, 406 students mentored by 44 postgraduate university students participated in the project. Students worked mainly in pairs to process ~200 tissue samples cut from 53 curated reptile specimens representing 17 species. A total of 27 novel Cytochrome Oxidase subunit 1 (CO1) sequences were ultimately generated for 8 south-east Australian reptile species of the families Scincidae and Agamidae.

Integrating phylogeography and high-resolution X-ray CT reveals five new cryptic species and multiple hybrid zones among Australian earless dragons.

Cryptic lineages, comprising species complexes with deep genetic structuring across the landscape but without distinct morphological differences, impose substantial difficulties for systematists and taxonomists in determining true species diversity. Here, we present an integrative approach that combines data from phylogeography and geometric morphometric analyses of three-dimensional cranial models to revisit the uncertain taxonomy of earless dragons from southern and central Australia that at one time or another have been included under the name Tympanocryptis lineata. Our approach finds strong support for seven previously described species, and more importantly, five undescribed Tympanocryptis taxa for which we provide a taxonomic treatment. We also find evidence of introgression and hybridization in three discrete contact zones between lineages, supported by mitochondrial and nuclear genes, as well as morphological analyses. With a sampling design that includes at least five individuals for each genetic lineage with corresponding X-ray microcomputed tomography scans, we perform comparative evolutionary analyses to show that there is a significant phylogenetic signal in Tympanocryptis cranial shape. Our results demonstrate the importance of using multiple specimens in each genetic lineage, particularly in cases of potential hybridization, and that geometric morphometrics, when used in an integrative framework, is a powerful tool in species delimitation across cryptic lineages. Our results lay the groundwork for future evolutionary studies in this widespread group across multiple environmental types and identify several species of immediate conservation concern with a focus on T. petersi sp. nov. We suggest that this species has undergone significant population declines and warrants a full conservation assessment.

Genetic structure and diversity of the endangered growling grass frog in a rapidly urbanizing region.

Two pervasive and fundamental impacts of urbanization are the loss and fragmentation of natural habitats. From a genetic perspective, these impacts manifest as reduced genetic diversity and ultimately reduced genetic viability. The growling grass frog (Litoria raniformis) is listed as vulnerable to extinction in Australia, and endangered in the state of Victoria. Remaining populations of this species in and around the city of Melbourne are threatened by habitat loss, degradation and fragmentation due to urban expansion. We used mitochondrial DNA (mtDNA) and microsatellites to study the genetic structure and diversity of L. raniformis across Melbourne's urban fringe, and also screened four nuclear gene regions (POMC, RAG-1, Rhod and CRYBA1). The mtDNA and nuclear DNA sequences revealed low levels of genetic diversity throughout remnant populations of L. raniformis. However, one of the four regions studied, Cardinia, exhibited relatively high genetic diversity and several unique haplotypes, suggesting this region should be recognized as a separate Management Unit. We discuss the implications of these results for the conservation of L. raniformis in urbanizing landscapes, particularly the potential risks and benefits of translocation, which remains a contentious management approach for this species.

Genetic Connectivity among Populations of an Endangered Snake Species from Southeastern Australia (Hoplocephalus bungaroides, Elapidae).

For endangered species that persist as apparently isolated populations within a previously more extensive range, the degree of genetic exchange between those populations is critical to conservation and management. A lack of gene flow can exacerbate impacts of threatening processes and delay or prevent colonization of sites after local extirpation. The broad-headed snake, Hoplocephalus bungaroides, is a small venomous species restricted to a handful of disjunct reserves near Sydney, Australia. Mark-recapture studies have indicated low vagility for this ambush predator, suggesting that gene flow also may be low. However, our analyses of 11 microsatellite loci from 163 snakes collected in Morton National Park, from six sites within a 10-km diameter, suggest relatively high rates of gene flow among sites. Most populations exchange genes with each other, with one large population serving as a source area and smaller populations apparently acting as sinks. About half of the juvenile snakes, for which we could reliably infer parentage, were collected from populations other than those in which we collected their putative parents. As expected from the snakes' reliance on rocky outcrops during cooler months of the year, most gene flow appears to be along sandstone plateaux rather than across the densely forested valleys that separate plateaux. The unexpectedly high rates of gene flow on a landscape scale are encouraging for future conservation of this endangered taxon. For example, wildlife managers could conserve broad-headed snakes by restoring habitats near extant source populations in areas predicted to be least affected by future climate change.

Reproduction in shark-attacked sea turtles is supported by stress-reduction mechanisms.

Vertebrates exhibit varied behavioural and physiological tactics to promote reproductive success. We examined mechanisms that could enable female loggerhead turtles to undertake nesting activities and maintain seasonal reproduction despite recent shark injuries of varying severity. We proposed that endocrinal mechanisms that regulate both a turtle's stress response and reproductive ability are modified to promote successful and continued reproduction. Irrespective of the degree of injury, females did not exhibit increased levels of the stress hormone corticosterone, nor decreased levels of the reproductive steroid testosterone; hormone responses consistent with stress. When exposed to a capture stressor, females with shark injury did not exhibit any greater corticosterone response than controls. In addition, breeding females showed a reduced corticosterone stress response compared to non-breeding females. Reduced endocrinal responses following shark injury, and during breeding in general may, in part, enable females to maintain behavioural and physiological commitment to reproduction.

Limited effect of anthropogenic habitat fragmentation on molecular diversity in a rain forest skink, Gnypetoscincus queenslandiae.

To examine the effects of recent habitat fragmentation, we assayed genetic diversity in a rain forest endemic lizard, the prickly forest skink (Gnypetoscincus queenslandiae), from seven forest fragments and five sites in continuous forest on the Atherton tableland of northeastern Queensland, Australia. The rain forest in this region was fragmented by logging and clearing for dairy farms in the early 1900s and most forest fragments studied have been isolated for 50-80 years or nine to 12 skink generations. We genotyped 411 individuals at nine microsatellite DNA loci and found fewer alleles per locus in prickly forest skinks from small rain forest fragments and a lower ratio of allele number to allele size range in forest fragments than in continuous forest, indicative of a decrease in effective population size. In contrast, and as expected for populations with small neighbourhood sizes, neither heterozygosity nor variance in allele size differed between fragments and sites in continuous forests. Considering measures of among population differentiation, there was no increase in FST among fragments and a significant isolation by distance pattern was identified across all 12 sites. However, the relationship between genetic (FST) and geographical distance was significantly stronger for continuous forest sites than for fragments, consistent with disruption of gene flow among the latter. The observed changes in genetic diversity within and among populations are small, but in the direction predicted by the theory of genetic erosion in recently fragmented populations. The results also illustrate the inherent difficulty in detecting genetic consequences of recent habitat fragmentation, even in genetically variable species, and especially when effective population size and dispersal rates are low.

Interplay between plasma hormone profiles, sex and body condition in immature hawksbill turtles (Eretmochelys imbricata) subjected to a capture stress protocol.

We investigated plasma hormone profiles of corticosterone and testosterone in immature hawksbill turtles (Eretmochelys imbricata) in response to a capture stress protocol. Further, we examined whether sex and body condition were covariates associated with variation in the adrenocortical response of immature turtles. Hawksbill turtles responded to the capture stress protocol by significantly increasing plasma levels of corticosterone over a 5 h period. There was no significant sex difference in the corticosterone stress response of immature turtles. Plasma testosterone profiles, while significantly different between the sexes, did not exhibit a significant change during the 5 h capture stress protocol. An index of body condition was not significantly associated with a turtle's capacity to produce plasma corticosterone both prior to and during exposure to the capture stress protocol. In summary, while immature hawksbill turtles exhibited an adrenocortical response to a capture stress protocol, neither their sex nor body condition was responsible for variation in endocrine responses. This lack of interaction between the adrenocortical response and these internal factors suggests that the inactive reproductive- and the current energetic- status of these immature turtles are important factors that could influence plasma hormone profiles during stress.