Genetically Encoded Lizard Color Divergence for Camouflage and Thermoregulation.

Local adaptation is critical in speciation and evolution, yet comprehensive studies on proximate and ultimate causes of local adaptation are generally scarce. Here, we integrated field ecological experiments, genome sequencing, and genetic verification to demonstrate both driving forces and molecular mechanisms governing local adaptation of body coloration in a lizard from the Qinghai-Tibet Plateau. We found dark lizards from the cold meadow population had lower spectrum reflectance but higher melanin contents than light counterparts from the warm dune population. Additionally, the colorations of both dark and light lizards facilitated the camouflage and thermoregulation in their respective microhabitat simultaneously. More importantly, by genome resequencing analysis, we detected a novel mutation in Tyrp1 that underpinned this color adaptation. The allele frequencies at the site of SNP 459# in the gene of Tyrp1 are 22.22% G/C and 77.78% C/C in dark lizards and 100% G/G in light lizards. Model-predicted structure and catalytic activity showed that this mutation increased structure flexibility and catalytic activity in enzyme TYRP1, and thereby facilitated the generation of eumelanin in dark lizards. The function of the mutation in Tyrp1 was further verified by more melanin contents and darker coloration detected in the zebrafish injected with the genotype of Tyrp1 from dark lizards. Therefore, our study demonstrates that a novel mutation of a major melanin-generating gene underpins skin color variation co-selected by camouflage and thermoregulation in a lizard. The resulting strong selection may reinforce adaptive genetic divergence and enable the persistence of adjacent populations with distinct body coloration.

Sex-specific splicing of Z- and W-borne nr5a1 alleles suggests sex determination is controlled by chromosome conformation.

Pogona vitticeps has female heterogamety (ZZ/ZW), but the master sex-determining gene is unknown, as it is for all reptiles. We show that nr5a1 (Nuclear Receptor Subfamily 5 Group A Member 1), a gene that is essential in mammalian sex determination, has alleles on the Z and W chromosomes (Z-nr5a1 and W-nr5a1), which are both expressed and can recombine. Three transcript isoforms of Z-nr5a1 were detected in gonads of adult ZZ males, two of which encode a functional protein. However, ZW females produced 16 isoforms, most of which contained premature stop codons. The array of transcripts produced by the W-borne allele (W-nr5a1) is likely to produce truncated polypeptides that contain a structurally normal DNA-binding domain and could act as a competitive inhibitor to the full-length intact protein. We hypothesize that an altered configuration of the W chromosome affects the conformation of the primary transcript generating inhibitory W-borne isoforms that suppress testis determination. Under this hypothesis, the genetic sex determination (GSD) system of P. vitticeps is a W-borne dominant female-determining gene that may be controlled epigenetically.

A reverse genetic approach in geckos with the CRISPR/Cas9 system by oocyte microinjection.

Reptiles are important model organisms in developmental and evolutionary biology, but are used less widely than other amniotes such as mouse and chicken. One of the main reasons for this is that has proven difficult to conduct CRISPR/Cas9-mediated genome editing in many reptile species despite the widespread use of this technology in other taxa. Certain features of reptile reproductive systems make it difficult to access one-cell or early-stage zygotes, which represents a key impediment to gene editing techniques. Recently, Rasys and colleagues reported a genome editing method using oocyte microinjection that allowed them to produce genome-edited Anolis lizards. This method opened a new avenue to reverse genetics studies in reptiles. In the present article, we report the development of a related method for genome editing in the Madagascar ground gecko (Paroedura picta), a well-established experimental model, and describe the generation of Tyr and Fgf10 gene-knockout geckos in the F0 generation.

Paranannizziopsis spp. Infection in Wild Vipers, Europe.

We describe the detection of Paranannizziopsis sp. fungus in a wild population of vipers in Europe. Fungal infections were severe, and 1 animal likely died from infection. Surveillance efforts are needed to better understand the threat of this pathogen to snake conservation.

Forelimb muscle activation patterns in American alligators: Insights into the evolution of limb posture and powered flight in archosaurs.

The evolution of archosaurs provides an important context for understanding the mechanisms behind major functional transformations in vertebrates, such as shifts from sprawling to erect limb posture and the acquisition of powered flight. While comparative anatomy and ichnology of extinct archosaurs have offered insights into musculoskeletal and gait changes associated with locomotor transitions, reconstructing the evolution of motor control requires data from extant species. However, the scarcity of electromyography (EMG) data from the forelimb, especially of crocodylians, has hindered understanding of neuromuscular evolution in archosaurs. Here, we present EMG data for nine forelimb muscles from American alligators during terrestrial locomotion. Our aim was to investigate the modulation of motor control across different limb postures and examine variations in motor control across phylogeny and locomotor modes. Among the nine muscles examined, m. pectoralis, the largest forelimb muscle and primary shoulder adductor, exhibited significantly smaller mean EMG amplitudes for steps in which the shoulder was more adducted (i.e., upright). This suggests that using a more adducted limb posture helps to reduce forelimb muscle force and work during stance. As larger alligators use a more adducted shoulder and hip posture, the sprawling to erect postural transition that occurred in the Triassic could be either the cause or consequence of the evolution of larger body size in archosaurs. Comparisons of EMG burst phases among tetrapods revealed that a bird and turtle, which have experienced major musculoskeletal transformations, displayed distinctive burst phases in comparison to those from an alligator and lizard. These results support the notion that major shifts in body plan and locomotor modes among sauropsid lineages were associated with significant changes in muscle activation patterns.

Vomeronasal organ volume increases with body size and is dissociated with the loss of a visual signal in Sceloporus lizards.

Many organisms communicate using signals in different sensory modalities (multicomponent or multimodal). When one signal or component is lost over evolutionary time, it may be indicative of changes in other characteristics of the signalling system, including the sensory organs used to perceive and process signals. Sceloporus lizards predominantly use chemical and visual signals to communicate, yet some species have lost the ancestral ventral colour patch used in male-male agonistic interactions and exhibit increased chemosensory behaviour. Here, we asked whether evolutionary loss of this sexual signal is associated with larger vomeronasal organ (VNO) volumes (an organ that detects chemical scents) compared with species that have retained the colour patch. We measured VNO coronal section areas of 7-8 adult males from each of 11 Sceloporus species (4 that lost and 7 that retained the colour patch), estimated sensory and total epithelium volume, and compared volumes using phylogenetic analysis of covariance, controlling for body size. Contrary to expectations, we found that species retaining the ventral patch had similar relative VNO volumes as did species that had lost the ancestral patch, and that body size explains VNO epithelium volume. Visual signal loss may be sufficiently compensated for by increased chemosensory behaviour, and the allometric pattern may indicate sensory system trade-offs for large-bodied species.

Neural thyroid hormone metabolism integrates seasonal changes in environmental temperature with the neuroendocrine reproductive axis.

We asked if environmental temperature alters thyroid hormone metabolism within the hypothalamus, thereby providing a neuroendocrine mechanism by which temperature could be integrated with photoperiod to regulate seasonal rhythms. We used immunohistochemistry to assess the effects of low-temperature winter dormancy at 4 °C or 12 °C on thyroid-stimulating hormone (TSH) within the infundibulum of the pituitary as well as deiodinase 2 (Dio2) and 3 (Dio3) within the hypothalamus of red-sided garter snakes (Thamnophis sirtalis). Both the duration and, in males, magnitude of low-temperature dormancy altered deiodinase immunoreactivity within the hypothalamus, increasing the area of Dio2-immunoreactivity in males and females and decreasing the number of Dio3-immunoreactive cells in males after 8-16 weeks. Reciprocal changes in Dio2/3 favor the accumulation of triiodothyronine within the hypothalamus. Whether TSH mediates these effects requires further study, as significant changes in TSH-immunoreactive cell number were not observed. Temporal changes in deiodinase immunoreactivity coincided with an increase in the proportion of males exhibiting courtship behavior as well as changes in the temporal pattern of courtship behavior after emergence. Our findings mirror those of previous studies, in which males require low-temperature exposure for at least 8 weeks before significant changes in gonadotropin-releasing hormone immunoreactivity and sex steroid hormones are observed. Collectively, these data provide evidence that the neuroendocrine pathway regulating the reproductive axis via thyroid hormone metabolism is capable of transducing temperature information. Because all vertebrates can potentially use temperature as a supplementary cue, these results are broadly applicable to understanding how environment-organism interactions mediate seasonally adaptive responses.

High sugar diet alters immune function and the gut microbiome in juvenile green iguanas (Iguana iguana).

The present work aimed to study whether a high sugar diet can alter immune responses and the gut microbiome in green iguanas. Thirty-six iguanas were split into four treatment groups using a 2×2 design. Iguanas received either a sugar-supplemented diet or a control diet, and either a lipopolysaccharide (LPS) injection or a phosphate-buffered saline (PBS) injection. Iguanas were given their respective diet treatment through the entire study (∼3 months) and received a primary immune challenge 1 and 2 months into the experiment. Blood samples and cloacal swabs were taken at various points in the experiment and used to measure changes in the immune system (bacterial killing ability, lysis and agglutination scores, LPS-specific IgY concentrations), and alterations in the gut microbiome. We found that a sugar diet reduces bacterial killing ability following an LPS challenge, and sugar and the immune challenge temporarily alters gut microbiome composition while reducing alpha diversity. Although sugar did not directly reduce lysis and agglutination following the immune challenge, the change in these scores over a 24-h period following an immune challenge was more drastic (it decreased) relative to the control diet group. Moreover, sugar increased constitutive agglutination outside of the immune challenges (i.e. pre-challenge levels). In this study, we provide evidence that a high sugar diet affects the immune system of green iguanas (in a disruptive manner) and alters the gut microbiome.

Glucocorticoid and glycemic responses to immune challenge in a viviparous snake afflicted with an emerging mycosis.

Disease may be both a cause and consequence of stress, and physiological responses to infectious disease may involve stress coping mechanisms that have important fitness consequences. For example, glucocorticoid and glycemic responses may affect host fitness by altering resource allocation and use in hosts, and these responses may be affected by competing stressors. To better understand the factors that affect host responses to infection, we challenged the immune system of field acclimatized pygmy rattlesnakes, Sistrurus miliarius, with a sterile antigen, lipopolysaccharide (LPS), and measured the glucocorticoid and glycemic response in healthy non-reproductive snakes, snakes afflicted with an emerging mycosis (ophidiomycosis), and pregnant snakes. We hypothesized that LPS challenge would result in a glucocorticoid and glycemic response typical of the vertebrate acute phase response (APR), and therefore predicted that LPS challenge would result in an acute increase in plasma corticosterone (CORT) and a decline in plasma glucose in all individuals. Additionally, we hypothesized that the APR would be attenuated in individuals simultaneously coping with additional challenges to homeostasis (i.e., disease or reproduction). As predicted, immune challenge elicited an acute increase in plasma CORT and a decrease in plasma glucose. Snakes coping with ophidiomycosis and pregnant snakes were able to mount a robust glucocorticoid and hypoglycemic response to LPS challenge, which was contrary to our hypothesis. Our findings clarify directions of causality linking infection, glucocorticoids, and glucose, and emphasize the importance of future research examining the fitness consequences of interactions between stress and disease in wildlife threatened by emerging pathogens.

The relevance of vascular adjustments to hemodynamic control in the face of temperature change in Crotalus durissus.

The presence of cardiac shunts in ectothermic tetrapods is thought to be consistent with active vascular modulations for proper hemodynamic support. Local control of blood flow modulates tissue perfusion and thus systemic conductance (Gsys) is assumed to increase with body temperature (Tb) to accommodate higher aerobic demand. However, the general increase of Gsys presses for a higher right-to-left (R-L) shunt, which reduces arterial oxygen concentration. In contrast, Tb reduction leads to a Gsys decrease and a left-to-right shunt, which purportedly increases pulmonary perfusion and plasma filtration in the respiratory area. This investigation addressed the role of compensatory vascular adjustments in the face of the metabolic alterations caused by Tb change in the South American rattlesnake (Crotalus durissus). Cardiovascular recordings were performed in decerebrated rattlesnake preparations at 10, 20 and 30°C. The rise in Tb increased metabolic demand, and correlated with an augmentation in heart rate. Although cardiac output increased, systemic stroke volume reduced while pulmonary stroke volume remained stable. Although that resulted in a proportionally higher increase in pulmonary blood flow, the R-L shunt was maintained. While the systemic compliance of large arteries was the most relevant factor in regulating arterial systemic blood pressure, peripheral conductance of pulmonary circulation was the major factor influencing the final cardiac shunt. Such dynamic adjustment of systemic compliance and pulmonary resistance for shunt modulation has not been demonstrated before and contrasts with previous knowledge on shunt control.

Stable isotopes reveal sex- and context-dependent amino acid routing in green anole lizards (Anolis carolinensis).

Allocation of acquired resources to phenotypic traits is affected by resource availability and current selective context. While differential investment in traits is well documented, the mechanisms driving investment at lower levels of biological organization, which are not directly related to fitness, remain poorly understood. We supplemented adult male and female Anolis carolinensis lizards with an isotopically labelled essential amino acid (13C-leucine) to track routing in four tissues (muscle, liver, gonads and spleen) under different combinations of resource availability (high- and low-calorie diets) and exercise training (sprint training and endurance capacity). We predicted sprint training should drive routing to muscle, and endurance training to liver and spleen, and that investment in gonads should be of lower priority in each of the cases of energetic stress. We found complex interactions between training regime, diet and tissue type in females, and between tissue type and training, and tissue type and diet in males, suggesting that males and females adjust their 13C-leucine routing strategies differently in response to similar environmental challenges. Importantly, our data show evidence of increased 13C-leucine routing in training contexts not to muscle as we expected, but to the spleen, which turns over blood cells, and to the liver, which supports metabolism under differing energetic scenarios. Our results reveal the context-specific nature of long-term trade-offs associated with increased chronic activity. They also illustrate the importance of considering the costs of locomotion in studies of life-history strategies.

There is no limitation for CO2 excretion across the lung in exercising American alligators (Alligator mississippiensis).

Vertebrates utilize various respiratory organs such as gills, lungs and skin in combination with diverse cardiovascular structures, including single-, three- and four-chambered hearts, to enable oxygen delivery and carbon dioxide removal. They also exhibit differences in aerobic and anaerobic metabolism during exertion, but the cardiorespiratory gas transport of all vertebrates is a four-step process governed by Fick's Principle and Fick's Law of Diffusion over the entire range of metabolic rates. Hillman et al. (2013) suggested that previous exercise studies have focused too narrowly on mammals and proposed that the cardiorespiratory system's excess capacity serves an evolutionary role in enhancing CO2 excretion in non-mammalian vertebrates. In contrast, an analysis by Hicks and Wang (2021) concluded that vertebrates maintain effective gas exchange even at peak activity, finding no evidence of arterial hypercapnia at maximal oxygen consumption and thus challenging the proposal of significant limitations to pulmonary or branchial CO2 efflux. In the present study, we investigated the limits for CO2 exchange in exercising American alligators (Alligator mississippiensis) and provide evidence that the cardiorespiratory system is adequately built to sustain CO2 excretion during strenuous exercise and maintain arterial PCO2, with no evidence of diffusion limitation for pulmonary CO2 excretion.

Sex and early-life conditions shape telomere dynamics in an ectotherm.

Telomeres, the repetitive DNA regions that protect the ends of chromosomes, and their shortening have been linked to key life history trade-offs among growth, reproduction and lifespan. In contrast to most endotherms, many ectotherms can compensate for telomere shortening throughout life by upregulation of telomerase in somatic tissues. However, during development, marked by rapid growth and an increased sensitivity to extrinsic factors, the upregulation of telomerase may be overwhelmed, resulting in long-term impacts on telomere dynamics. In ectotherms, one extrinsic factor that may play a particularly important role in development is temperature. Here, we investigated the influence of developmental temperature and sex on early-life telomere dynamics in an oviparous ectotherm, Lacerta agilis. While there was no effect of developmental temperature on telomere length at hatching, there were subsequent effects on telomere maintenance capacity, with individuals incubated at warm temperatures exhibiting less telomere maintenance compared with cool-incubated individuals. Telomere dynamics were also sexually dimorphic, with females having longer telomeres and greater telomere maintenance compared with males. We suggest that selection drives this sexual dimorphism in telomere maintenance, in which females maximise their lifetime reproductive success by investing in traits promoting longevity such as maintenance, while males invest in short-term reproductive gains through a polygynous mating behaviour. These early-life effects, therefore, have the potential to mediate life-long changes to life histories.

Trained quantity discrimination in invasive red-eared slider and a comparison with the native stripe-necked turtle.

Little is known about the behavioral and cognitive traits that best predict invasion success. Evidence is mounting that cognitive performance correlates with survival and fecundity, two pivotal factors for the successful establishment of invasive populations. We assessed the quantity discrimination ability of the globally invasive red-eared slider (Trachemys scripta elegans). We further compared it to that of the native stripe-necked turtle (Mauremys sinensis), which has been previously evaluated for its superior quantity discrimination ability. Specifically, our experimental designs aimed to quantify the learning ability as numerosity pairs increased in difficulty (termed fixed numerosity tests), and the immediate response when turtles were presented with varied challenges concurrently in the same tests (termed mixed numerosity tests). Our findings reaffirm the remarkable ability of freshwater turtles to discern numerical differences as close as 9 vs 10 (ratio = 0.9), which was comparable to the stripe-necked turtle's performance. However, the red-eared slider exhibited a moderate decrease in performance in high ratio tests, indicating a potentially enhanced cognitive capacity to adapt to novel challenges. Our experimental design is repeatable and is adaptable to a range of freshwater turtles. These findings emphasize the potential importance of cognitive research to the underlying mechanisms of successful species invasions.

Severe fire regimes decrease resilience of ectothermic populations.

Understanding populations' responses to environmental change is crucial for mitigating human-induced disturbances. Here, we test hypotheses regarding how three essential components of demographic resilience (resistance, compensation and recovery) co-vary along the distinct life histories of three lizard species exposed to variable, prescribed fire regimes. Using a Bayesian hierarchical framework, we estimate vital rates (survival, growth and reproduction) with 14 years of monthly individual-level data and mark-recapture models to parameterize stochastic integral projection models from five sites in Brazilian savannas, each historically subjected to different fire regimes. With these models, we investigate how weather, microclimate and ecophysiological traits of each species influence their vital rates, emergent life history traits and demographic resilience components in varying fire regimes. Overall, weather and microclimate are better predictors of the species' vital rates, rather than their ecophysiological traits. Our findings reveal that severe fire regimes increase populations' resistance but decrease compensation or recovery abilities. Instead, populations have higher compensatory and recovery abilities at intermediate degrees of fire severity. Additionally, we identify generation time and reproductive output as predictors of resilience trends across fire regimes and climate. Our analyses demonstrate that the probability and quantity of monthly reproduction are the proximal drivers of demographic resilience across the three species. Our findings suggest that populations surpass a tipping point in severe fire regimes and achieve an alternative stable state to persist. Thus, higher heterogeneity in fire regimes can increase the reproductive aspects and resilience of different populations and avoid high-severity regimes that homogenize the environment. Despite being more resistant, species with long generation times and low reproductive output take longer to recover and cannot compensate as much as species with faster paces of life. We emphasize how reproductive constraints, such as viviparity and fixed clutch sizes, impact the ability of ectothermic populations to benefit and recover from disturbances, underscoring their relevance in conservation assessments.

Loggerhead Sea Turtles as Hosts of Diverse Bacterial and Fungal Communities.

Research on microbial communities associated with wild animals provides a valuable reservoir of knowledge that could be used for enhancing their rehabilitation and conservation. The loggerhead sea turtle (Caretta caretta) is a globally distributed species with its Mediterranean population categorized as least concern according to the IUCN Red List of Threatened Species as a result of robust conservation efforts. In our study, we aimed to further understand their biology in relation to their associated microorganisms. We investigated epi- and endozoic bacterial and endozoic fungal communities of cloaca, oral mucosa, carapace biofilm. Samples obtained from 18 juvenile, subadult, and adult turtles as well as 8 respective enclosures, over a 3-year period, were analysed by amplicon sequencing of 16S rRNA gene and ITS2 region of nuclear ribosomal gene. Our results reveal a trend of decreasing diversity of distal gut bacterial communities with the age of turtles. Notably, Tenacibaculum species show higher relative abundance in juveniles than in adults. Differential abundances of taxa identified as Tenacibaculum, Moraxellaceae, Cardiobacteriaceae, and Campylobacter were observed in both cloacal and oral samples in addition to having distinct microbial compositions with Halioglobus taxa present only in oral samples. Fungal communities in loggerheads' cloaca were diverse and varied significantly among individuals, differing from those of tank water. Our findings expand the known microbial diversity repertoire of loggerhead turtles, highlighting interesting taxa specific to individual body sites. This study provides a comprehensive view of the loggerhead sea turtle bacterial microbiota and marks the first report of distal gut fungal communities that contributes to establishing a baseline understanding of loggerhead sea turtle holobiont.

Research and application of a novel selective stacking ensemble model based on error compensation and parameter optimization for AQI prediction.

With the ongoing process of industrialization, the issue of declining air quality is increasingly becoming a critical concern. Accurate prediction of the Air Quality Index (AQI), considered as an all-inclusive measure representing the extent of pollutants present in the atmosphere, is of paramount importance. This study introduces a novel methodology that combines stacking ensemble and error correction to improve AQI prediction. Additionally, the reptile search algorithm (RSA) is employed for optimizing model parameters. In this study, four distinct regional AQI data containing a collection of 34864 data samples are collected. Initially, we perform cross-validation on ten commonly used single models to obtain prediction results. Then, based on evaluation indices, five models are selected for ensemble. The results of the study show that the model proposed in this paper achieves an improvement of around 10% in terms of accuracy when compared to the conventional model. Thus, the model introduced in this study offers a more scientifically grounded approach in tackling air pollution.

The genome of the Australian water dragon (Intellagama lesueurii), an agamid model for urban adaptation.

Squamate reptiles are a highly diverse and intriguing group of tetrapods, offering valuable insights into the evolution of amniotes. The Australian water dragon (Intellagama lesueurii) is a member of the Agamidae, and sister to the core mesic Australian endemic radiation (Amphibolurinae). The species is renowned for its urban adaptability and complex social systems. We report a 1.8 Gb chromosome-length genome assembly together with the annotation of 23,675 protein-coding genes. Comparative analysis with other squamate genomes highlights gene family expansions associated with immune function, energetic homeostasis, and wound healing. This reference genome will serve as a valuable resource for studies of evolution and environmental resilience in lizards.

Hemodynamics During Development and Postnatal Life.

A well-developed heart is essential for embryonic survival. There are constant interactions between cardiac tissue motion and blood flow, which determine the heart shape itself. Hemodynamic forces are a powerful stimulus for cardiac growth and differentiation. Therefore, it is particularly interesting to investigate how the blood flows through the heart and how hemodynamics is linked to a particular species and its development, including human. The appropriate patterns and magnitude of hemodynamic stresses are necessary for the proper formation of cardiac structures, and hemodynamic perturbations have been found to cause malformations via identifiable mechanobiological molecular pathways. There are significant differences in cardiac hemodynamics among vertebrate species, which go hand in hand with the presence of specific anatomical structures. However, strong similarities during development suggest a common pattern for cardiac hemodynamics in human adults. In the human fetal heart, hemodynamic abnormalities during gestation are known to progress to congenital heart malformations by birth. In this chapter, we discuss the current state of the knowledge of the prenatal cardiac hemodynamics, as discovered through small and large animal models, as well as from clinical investigations, with parallels gathered from the poikilotherm vertebrates that emulate some hemodynamically significant human congenital heart diseases.

Developmental energetics in the oviparous corn snake, Pantherophis guttatus, confirms a conservative evolutionary pattern in snakes.

Eggs of oviparous reptiles are ideal models for studying evolutionary patterns of embryonic metabolism since they allow tracking of energy allocation during development. Analyzing oxygen consumption of whole eggs throughout development indicates three patterns among reptiles. Embryos initially grow and consume oxygen exponentially, but oxygen consumption slows, or drops before hatching in some species. Turtles, crocodilians, and most lizards follow curves with initial exponential increases followed by declines, whereas embryonic snakes that have been studied exhibit a consistently exponential pattern. This study measured oxygen consumption of corn snake, Pantherophis guttatus, embryos to determine if this species also exhibits an exponential increase in oxygen consumption. Individual eggs, sampled weekly from oviposition to hatching, were placed in respirometry chambers for 24-h during which oxygen consumption was recorded. Embryos were staged and carcasses and yolk were weighed separately. Results indicate steady inclines in oxygen consumption during early stages of development, with a rapid increase prior to hatching. The findings support the hypothesis that embryonic oxygen consumption of snakes differs from most other non-avian reptiles. Total energy required for development was determined based on calorimetry of initial yolk compared to hatchlings and residual yolk and by integration of the area under the curve plotting oxygen consumption versus age of embryos. The cost of development estimates based on these two methods were 6.4 and 10.0 kJ, respectively. Our results emphasize the unique physiological aspects of snake embryogenesis and illustrate how the study of physiological characteristics can contribute to the broader understanding of reptilian evolution.