Intraspecific differences in thermal tolerance of the diamondback watersnake (Nerodia rhombifer): effects of ontogeny, latitude, and sex

Sex differences in antioxidant ability and energy metabolism level resulting in the difference of hypoxia tolerance in red swamp crayfish (Procambarus clarkii)

Sexual dimorphism widely exists in crustaceans. However, sex differences in the hypoxia tolerance of crayfish have rarely been reported. In this study, the differences in hypoxia tolerance between the two sexes of crayfish were assessed according to mortality, pathological features of hepatopancreas, antioxidant enzyme activity and differentially expressed genes (DEGs) analysis using transcriptome. The results showed that male crayfish displayed significantly higher mortality than the female under hypoxia stress (p < 0.05). Furthermore, female crayfish demonstrated higher levels of antioxidant enzyme activity. Hematoxylin-eosin staining analysis revealed that the damage of hepatopancreas was more severe in the male crayfish compared to the female crayfish. Additionally, there was higher expression level of the DEGs in hypoxia-inducible factor (HIF) pathway and higher energy metabolism level in the female compared to the male. Together, these findings suggest that the female crayfish with higher antioxidant ability and energy metabolism level exhibits stronger hypoxia tolerance than the male crayfish, providing the theoretical support for investigating sex differences in hypoxia tolerance among crustaceans.

High-temperature stress will put the thermo-sensitive teleost yellow catfish (Tachysurus fulvidraco) in danger through reducing reproductivity

Recently, concerns for species that sex differentiation is influenced by temperature in the context of global warming have increased because disrupted operational sex ratios could threaten population maintenance. In contrast, little attention has been given to the reproductive ability of populations that experienced elevated temperatures. In this study, we demonstrated that high temperature (HT) would decrease population size via three different aspects of reproductive ability for the first time. We show that, in a thermo-sensitive teleost yellow catfish, a short period of HT (+3 °C) exposure during the critical period of sex differentiation leads to a different percentage of masculinization of XX genotypic females (1-23%) in wet-lab and natural water bodies. Combining the results of gonadal appearance, histology, sperm parameters, and fertilization rate, we found that XX pseudo-males induced by HT display significantly discounted fertility and reproductive performance compared to XY normal males. We demonstrate that the survival of the XY genotype is lower than XX genotype under environmental stress, including HT, hypoxia, and parasite infection, and the differential survival seems unrelated to male-biased sexual size dimorphism. The mathematical model predicts that the phenotypic female percent will be stabilized at 50% and the population will be sustainably maintained when masculinizing force is less than 0.5, while HT will put the population in danger when the masculinizing force exceeds 0.5. However, when we combine the real-world data of reproductive ability and mathematic model, our results suggest the population size decreases and the long-term survival of the studied species are threatened under the projected pace of increasing temperature. These findings will be useful for understanding the long-term effects of increasing temperature on sex ratio, reproduction and population maintenance in teleost.

Variation in upper thermal tolerance among 19 species from temperate wetlands

Communities usually possess a multitude of interconnected trophic interactions within food webs. Their regulation generally depends on a balance between bottom-up and top-down effects. However, if sensitivity to temperature varies among species, rising temperatures may change trophic interactions via direct and indirect effects. We examined the critical thermal maximum (CT_max) of 19 species from temperate wetlands (insect predators, amphibian larvae, zooplankton and amphipods) and determined if they vary in their sensitivity to warming temperatures. CT_max differed between the groups, with predatory insects having higher CT_max than amphibians (both herbivorous larval anurans and predatory larval salamanders), amphipods and zooplankton. In a scenario of global warming, these differences in thermal tolerance may affect top-down and bottom-up processes, particularly considering that insect predators are more likely to maintain or improve their performance at higher temperatures, which could lead to increased predation rates on the herbivores in the food web. Further studies are needed to understand how the energy flows through communities, how species' energy budgets may change and whether other physiological and behavioral responses (such as phenotypic plasticity and thermoregulation) can buffer or increase these changes in the top-down regulation of wetland food webs.

Thermal performance curves reveal shifts in optima, limits and breadth in early life

Understanding thermal performance at life stages that limit persistence is necessary to predict responses to climate change, especially for ectotherms whose fitness (survival and reproduction) depends on environmental temperature. Ectotherms often undergo stage-specific changes in size, complexity and duration that are predicted to modify thermal performance. Yet performance is mostly explored for adults, while performance at earlier stages that typically limit persistence remains poorly understood. Here, we experimentally isolate thermal performance curves at fertilization, embryo development and larval development stages in an aquatic ectotherm whose early planktonic stages (gametes, embryos and larvae) govern adult abundances and dynamics. Unlike previous studies based on short-term exposures, responses with unclear links to fitness or proxies in lieu of explicit curve descriptors (thermal optima, limits and breadth), we measured performance as successful completion of each stage after exposure throughout, and at temperatures that explicitly capture curve descriptors at all stages. Formal comparisons of descriptors using a combination of generalized linear mixed modelling and parametric bootstrapping reveal important differences among life stages. Thermal performance differs significantly from fertilization to embryo development (with thermal optimum declining by ∼2°C, thermal limits shifting inwards by ∼8-10°C and thermal breadth narrowing by ∼10°C), while performance declines independently of temperature thereafter. Our comparisons show that thermal performance at one life stage can misrepresent performance at others, and point to gains in complexity during embryogenesis, rather than subsequent gains in size or duration of exposure, as a key driver of thermal sensitivity in early life.

How repeatable is CTmax within individual brook trout over short- and long-time intervals?

As stream temperatures increase due to factors such as heated runoff from impervious surfaces, deforestation, and climate change, fish species adapted to cold water streams are forced to move to more suitable habitat, acclimate or adapt to increased thermal regimes, or die. To estimate the potential for adaptation, a (within individual) repeatable metric of thermal tolerance is imperative. Critical thermal maximum (CT_max) is a dynamic test that is widely used to measure thermal tolerance across many taxa and has been used in fishes for decades, but its repeatability in most species is unknown. CT_max tests increase water temperature steadily over time until loss of equilibrium (LOE) is achieved. To determine if CT_max is a consistent metric within individual fish, we measured CT_max on the same lab-held individually-marked adult brook trout Salvelinus fontinalis at three different times (August & September 2016, September 2017). We found that CT_max is a repeatable trait (Repeatability ± S.E.: 0.48 ± 0.14). CT_max of individuals males was consistent over time, but the CT_max of females increased slightly over time. This result indicates that CT_max is a robust, repeatable estimate of thermal tolerance in a cold-water adapted fish.

Thermal performance responses in free-ranging elasmobranchs depend on habitat use and body size

Temperature is one of the most influential drivers of physiological performance and behaviour in ectotherms, determining how these animals relate to their ecosystems and their ability to succeed in particular habitats. Here, we analysed the largest set of acceleration data compiled to date for elasmobranchs to examine the relationship between volitional activity and temperature in 252 individuals from 8 species. We calculated activation energies for the thermal performance response in each species and estimated optimum temperatures using an Arrhenius breakpoint analysis, subsequently fitting thermal performance curves to the activity data. Juveniles living in confined nursery habitats not only spent substantially more time above their optimum temperature and at the upper limits of their performance breadths compared to larger, less site-restricted animals, but also showed lower activation energies and broader performance curves. Species or life stages occupying confined habitats featured more generalist behavioural responses to temperature change, whereas wider ranging elasmobranchs were characterised by more specialist behavioural responses. The relationships between the estimated performance regimes and environmental temperature limits suggest that animals in confined habitats, including many juvenile elasmobranchs within nursery habitats, are likely to experience a reduction of performance under a warming climate, although their flatter thermal response will likely dampen this impact. The effect of warming on less site-restricted species is difficult to forecast since three of four species studied here did not reach their optimum temperature in the wild, although their specialist performance characteristics may indicate a more rapid decline should optimum temperatures be exceeded.

Oxygen Limitation Does Not Drive the Decreasing Heat Tolerance of Grasshoppers during Development

Thermal physiology changes as organisms grow and develop, but we do not understand what causes these ontogenetic shifts. According to the theory of oxygen- and capacity-limited thermal tolerance, an organism's heat tolerance should change throughout ontogeny as its ability to deliver oxygen varies. As insects grow during an instar, their metabolic demand increases without a proportional increase in the size of tracheae that supply oxygen to the tissues. If oxygen delivery limits heat tolerance, the mismatch between supply and demand should make insects more susceptible to heat and hypoxia as they progress through an instar. We tested this hypothesis by measuring the heat tolerance of grasshoppers (Schistocerca americana) on the second and seventh days of the sixth instar, in either a normoxic or a hypoxic atmosphere (21% or 10% O₂, respectively). As expected, heat tolerance decreased as grasshoppers grew larger. Yet contrary to expectation, hypoxia had no effect on heat tolerance across all stages and sizes. Although heat tolerance declines as grasshoppers grow, this pattern must stem from a mechanism other than oxygen limitation.

Interactive effects of experimental heating rates, ontogeny and body mass on the upper thermal limits of anuran larvae

Biological and methodological factors influence the upper thermal limits (UTL) of ectothermic animals, but most factors have been studied independently. Few studies have integrated variables, so our understanding about sources of UTL variation remains fragmentary. Thereby, we investigated synergic effects of experimental protocols (heating rates, ΔTs) and biological factors (ontogeny and body mass) on the UTL on the larvae of two anuran species (Physalaemus nattereri and Boana pardalis), specifically their Critical Thermal Maximum (CT_max). The species displayed slightly different responses to ΔTs: In B. pardalis tadpoles both average and variance of CT_max increased at a fastest ΔT, the same response happened in P. nattereri tadpoles at slow and moderate ΔTs. Also, the CT_max of P. nattereri declined at the end of metamorphosis independently of ΔT, but tadpoles at all developmental stages still displayed higher heat tolerance at the slow ΔT. Finally, we detected small, synergic effects of body mass and ΔTs on the CT_max of both species. In small B. pardalis tadpoles and premetamorphic P. nattereri tadpoles, body mass had a positive effect on CT_max, but only at slow and moderate ΔTs, probably indicating physiological responses. A similar trend was observed in large B. pardalis tadpoles at the fast ΔT, but this result is likely to be influenced by thermal inertia. Our findings contribute to integrate the understanding of factors influencing UTL in small ectothermic animals. This understanding is critical to discuss the physiological component of vulnerability to climate change that is related to acute temperatures.

Upper thermal tolerances of different life stages, sexes, and species of widow spiders (Araneae, Theridiidae)

Temperature strongly influences the physiology and behavior of ectotherms. Persistence within different environments can be limited by thermal tolerances. These thermal tolerances can also shift through life stages and differ between sexes. The critical thermal maximum (CT_Max) defines the temperature at which animals experience unorganized locomotion or spasms. In this study, we tested if CT_Max varied between a native and an invasive widow species. We separately tested if CT_Max varied by widow life stage and sex. We predicted that the invasive species would have higher CT_Max due to originally inhabiting warmer climates. We also predicted that juveniles and male widows would possess higher CT_Max because they are more mobile and could experience a greater scope of thermal extremes throughout landscapes. We did not find a difference in CT_Max between the species, but we did find differences across development stages. Temperature of spasms and death decreased with developmental stages, which corresponds with previous studies in spiders. Future studies of ontogenic and interspecific comparisons will be crucial for more broadly understanding how upper tolerances shapes species persistence in changing climates or ability to invade new habitats.

Complex interactions between local adaptation, phenotypic plasticity and sex affect vulnerability to warming in a widespread marine copepod

Predicting the response of populations to climate change requires an understanding of how various factors affect thermal performance. Genetic differentiation is well known to affect thermal performance, but the effects of sex and developmental phenotypic plasticity often go uncharacterized. We used common garden experiments to test for effects of local adaptation, developmental phenotypic plasticity and individual sex on thermal performance of the ubiquitous copepod, Acartia tonsa (Calanoida, Crustacea) from two populations strongly differing in thermal regimes (Florida and Connecticut, USA). Females had higher thermal tolerance than males in both populations, while the Florida population had higher thermal tolerance compared with the Connecticut population. An effect of developmental phenotypic plasticity on thermal tolerance was observed only in the Connecticut population. Our results show clearly that thermal performance is affected by complex interactions of the three tested variables. Ignoring sex-specific differences in thermal performance may result in a severe underestimation of population-level impacts of warming because of population decline due to sperm limitation. Furthermore, despite having a higher thermal tolerance, low-latitude populations may be more vulnerable to warming as they lack the ability to respond to increases in temperature through phenotypic plasticity.

Thermal quality influences habitat use of two anole species

Regeneration of secondary forests on previously deforested or degraded land is one of the most dominant forms of land-use change in the tropics. However, the response of animal communities to forest regeneration is poorly understood. To evaluate support for thermal quality as a mechanism driving reptile species distributions during secondary forest succession, we measured operative temperatures and occupancy in three successional forest stages (pasture, secondary forest, and old growth forest) for two anole species common in the landscape (Norops humilis and Norops limifrons). We then measured thermal preference in laboratory experiments and used operative temperature and temperature preference measurements to determine how thermal quality of habitat changes over the course of secondary forest succession, and if occupancy varies as a function of thermal quality. We found that thermal quality was lowest in pasture habitat because of a large frequency of temperatures above the thermal preference range. However, in low thermal quality pasture sites, riparian habitats and remnant trees provided a thermal refuge for both lizard species. Our results support thermal quality as a mechanism for reptile species distributions in altered landscapes and highlight the importance of the maintenance of riparian corridors.

Integrating multiple stressors across life stages and latitudes: Combined and delayed effects of an egg heat wave and larval pesticide exposure in a damselfly

To understand the effects of pollutants in a changing world we need multistressor studies that combine pollutants with other stressors associated with global change such as heat waves. We tested for the delayed and combined impact of a heat wave during the egg stage and subsequent sublethal exposure to the pesticide esfenvalerate during the larval stage on life history and physiology in the larval and adult stage of the damselfly Lestes sponsa. We studied this in a common garden experiment with replicated central- and high latitude populations to explore potential effects of local thermal adaptation and differences in life history shaping the multistressor responses. Exposure of eggs to the heat wave had no effect on larval traits, yet had delayed costs (lower fat and flight muscle mass) in the adult stage thereby crossing two life history transitions. These delayed costs were only present in central-latitude populations potentially indicating their lower heat tolerance. Exposure of larvae to the pesticide reduced larval growth rate and prolonged development time, and across metamorphosis reduced the adult fat content and the flight muscle mass, yet did not affect the adult heat tolerance. The pesticide-induced delayed emergence was only present in the slower growing central-latitude larvae, possibly reflecting stronger selection to keep development fast in the more time-constrained high-latitude populations. We observed no synergistic interactions between the egg heat wave and the larval pesticide exposure. Instead the pesticide-induced reduction in fat content was only present in animals that were not exposed to the egg heat wave. Our results based on laboratory conditions highlight that multistressor studies should integrate across life stages to fully capture cumulative effects of pollutants with other stressors related to global change.

Effects of Maternally-Transferred Methylmercury on Stress Physiology in Northern Water Snake (Nerodia sipedon) Neonates

Biomagnification of methylmercury in aquatic systems can cause elevated tissue mercury (Hg) and physiological stress in top predators. Mercury is known to affect stress hormone levels in mammals, birds and fish. In this study, the effects of maternally-transferred methylmercury on the stress physiology of Northern Water Snake (Nerodia sipedon) neonates were tested. Gravid females were dosed via force-fed capsules during late gestation with 0, 0.01, or 10 µg methylmercury per gram of body mass. Plasma corticosterone levels and leukocyte differentials were analyzed in baseline and confinement-stressed neonates from all dose levels. Neither Hg nor confinement stress had a significant effect on leukocyte differentials nor was Hg related to corticosterone levels. However, stress group neonates showed lower heterophil/lymphocyte ratios and this study was the first to show that neonate N. sipedon can upregulate CORT in response to stress. These results indicate that N. sipedon may be somewhat tolerant to Hg contamination.

Heat tolerance during embryonic development has not diverged among populations of a widespread species (Sceloporus undulatus)

The frequency and magnitude of heat waves have increased in recent decades, imposing additional stresses on organisms in extreme environments. Most reptilian embryos are regularly exposed to thermal stress because they develop in shallow, warm soils for weeks to months. We studied cardiac performance during warming to infer lethal temperatures for embryonic lizards in the Sceloporus undulatus complex. Embryos from four populations throughout the geographical range (New Jersey, South Carolina, Colorado, and Arizona) were warmed at a rate observed in natural nests. Embryos from all populations exhibited a similar pattern of thermal sensitivity, as follows: heart rate rose between 34 and 41°C, remained stable between 41 and 44°C, and dropped sharply between 44 and 47°C. No embryos recovered from cardiac arrest, indicating that the upper lethal temperature was ≤47°C. Despite the putative selective pressures, the thermal limit to cardiac performance seems to have been conserved during the evolution of this species.

Multivariate analysis of adaptive capacity for upper thermal limits in Drosophila simulans

Thermal tolerance is an important factor influencing the distribution of ectotherms, but our understanding of the ability of species to evolve different thermal limits is limited. Based on univariate measures of adaptive capacity, it has recently been suggested that species may have limited evolutionary potential to extend their upper thermal limits under ramping temperature conditions that better reflect heat stress in nature. To test these findings more broadly, we used a paternal half-sibling breeding design to estimate the multivariate evolutionary potential for upper thermal limits in Drosophila simulans. We assessed heat tolerance using static (basal and hardened) and ramping assays. Our analyses revealed significant evolutionary potential for all three measures of heat tolerance. Additive genetic variances were significantly different from zero for all three traits. Our G matrix analysis revealed that all three traits would contribute to a response to selection for increased heat tolerance. Significant additive genetic covariances and additive genetic correlations between static basal and hardened heat-knockdown time, marginally nonsignificant between static basal and ramping heat-knockdown time, indicate that direct and correlated responses to selection for increased upper thermal limits are possible. Thus, combinations of all three traits will contribute to the evolution of upper thermal limits in response to selection imposed by a warming climate. Reliance on univariate estimates of evolutionary potential may not provide accurate insight into the ability of organisms to evolve upper thermal limits in nature.

Thermal tolerance of the crab Pachygrapsus marmoratus: intraspecific differences at a physiological (CTMax) and molecular level (Hsp70)

Temperature is one of the most important variables influencing organisms, especially in the intertidal zone. This work aimed to test physiological and molecular intraspecific differences in thermal tolerance of the crab Pachygrapsus marmoratus (Fabricius, 1787). The comparisons made focused on sex, size, and habitat (estuary and coast) differences. The physiological parameter was upper thermal limit, tested via the critical thermal maximum (CTMax) and the molecular parameter was total heat shock protein 70 (Hsp70 and Hsp70 plus Hsc70) production, quantified via an enzyme-linked imunosorbent assay. Results showed that CTMax values and Hsp70 production are higher in females probably due to different microhabitat use and potentially due to different hormonal regulation in males and females. Among females, non-reproducing ones showed a higher CTMax value, but no differences were found in Hsp70, even though reproducing females showed higher variability in Hsp70 amounts. As reproduction takes up a lot of energy, its allocation for other activities, including stress responses, is lower. Juveniles also showed higher CTMax and Hsp70 expression because they occur in greater shore heights and ageing leads to alterations in protein synthesis. Comparing estuarine and coastal crabs, no differences were found in CTMax but coastal crabs produce more Hsp70 than estuarine crabs because they occur in drier and hotter areas than estuarine ones, which occur in moister environments. This work shows the importance of addressing intraspecific differences in the stress response at different organizational levels. This study shows that these differences are key factors in stress research, climate research, and environmental monitoring.

Temperature tolerance and stress proteins as mechanisms of invasive species success

Invasive species are predicted to be more successful than natives as temperatures increase with climate change. However, few studies have examined the physiological mechanisms that theoretically underlie this differential success. Because correlative evidence suggests that invasiveness is related to the width of a species' latitudinal range, it has been assumed – but largely untested – that range width predicts breadth of habitat temperatures and physiological thermotolerances. In this study, we use empirical data from a marine community as a case study to address the hypotheses that (1) geographic temperature range attributes are related to temperature tolerance, leading to greater eurythermality in invasive species, and (2) stress protein expression is a subcellular mechanism that could contribute to differences in thermotolerance. We examined three native and six invasive species common in the subtidal epibenthic communities of California, USA. We assessed thermotolerance by exposing individuals to temperatures between 14°C and 31°C and determining the temperature lethal to 50% of individuals (LT₅₀) after a 24 hour exposure. We found a strong positive relationship between the LT₅₀ and both maximum habitat temperatures and the breadth of temperatures experience across the species' ranges. In addition, of the species in our study, invasives tended to inhabit broader habitat temperature ranges and higher maximum temperatures. Stress protein expression may contribute to these differences: the more thermotolerant, invasive species Diplosoma listerianum expressed higher levels of a 70-kDa heat-shock protein than the less thermotolerant, native Distaplia occidentalis for which levels declined sharply above the LT₅₀. Our data highlight differences between native and invasive species with respect to organismal and cellular temperature tolerances. Future studies should address, across a broader phylogenetic and ecosystem scope, whether this physiological mechanism has facilitated the current success of invasive species and could lead to greater success of invasives than native species as global warming continues.

Insect thermal tolerance: what is the role of ontogeny, ageing and senescence?

Temperature has dramatic evolutionary fitness consequences and is therefore a major factor determining the geographic distribution and abundance of ectotherms. However, the role that age might have on insect thermal tolerance is often overlooked in studies of behaviour, ecology, physiology and evolutionary biology. Here, we review the evidence for ontogenetic and ageing effects on traits of high- and low-temperature tolerance in insects and show that these effects are typically pronounced for most taxa in which data are available. We therefore argue that basal thermal tolerance and acclimation responses (i.e. phenotypic plasticity) are strongly influenced by age and/or ontogeny and may confound studies of temperature responses if unaccounted for. We outline three alternative hypotheses which can be distinguished to propose why development affects thermal tolerance in insects. At present no studies have been undertaken to directly address these options. The implications of these age-related changes in thermal biology are discussed and, most significantly, suggest that the temperature tolerance of insects should be defined within the age-demographics of a particular population or species. Although we conclude that age is a source of variation that should be carefully controlled for in thermal biology, we also suggest that it can be used as a valuable tool for testing evolutionary theories of ageing and the cellular and genetic basis of thermal tolerance.

Northern grass lizards (Takydromus septentrionalis) from different populations do not differ in thermal preference and thermal tolerance when acclimated under identical thermal conditions

We acclimated adults of Takydromus septentrionalis (northern grass lizard) from four localities (populations) under identical thermal conditions to examine whether local thermal conditions have a fixed influence on thermal preference and thermal tolerance in the species. Selected body temperature (Tsel), critical thermal minimum (CTMin), and critical thermal maximum (CTMax) did not differ between sexes and among localities in lizards kept under identical laboratory conditions for approximately 5 months, and the interaction effects between sex and locality on these measures were not significant. Lizards acclimated to the three constant temperatures (20, 25, and 35 degrees C) differed in Tsel, CTMin, and CTMax. Tsel, CTMin, and CTMax all shifted upward as acclimation temperature increased, with Tsel shifting from 32.0 to 34.1 degrees C, CTMin from 4.9 to 8.0 degrees C, and CTMax from 42.0 to 44.5 degrees C at the change-over of acclimation temperature from 20 to 35 degrees C. Lizards acclimated to the three constant temperatures also differed in the range of viable body temperatures; the range was widest in the 25 degrees C treatment (38.1 degrees C) and narrowest in the 35 degrees C treatment (36.5 degrees C), with the 20 degrees C treatment in between (37.2 degrees C). The results of this study show that local thermal conditions do not have a fixed influence on thermal preference and thermal tolerance in T. septentrionalis.

Ontogenetic shifts in thermal tolerance, selected body temperature and thermal dependence of food assimilation and locomotor performance in a lacertid lizard, Eremias brenchleyi

We used Eremias brenchleyi as a model animal to examine differences in thermal tolerance, selected body temperature, and the thermal dependence of food assimilation and locomotor performance between juvenile and adult lizards. Adults selected higher body temperatures (33.5 vs. 31.7 degrees C) and were able to tolerate a wider range of body temperatures (3.4-43.6 vs. 5.1-40.8 degrees C) than juveniles. Within the body temperature range of 26-38 degrees C, adults overall ate more than juveniles, and food passage rate was faster in adults than juveniles. Apparent digestive coefficient (ADC) and assimilation efficiency (AE) varied among temperature treatments but no clear temperature associated patterns could be discerned for these two variables. At each test temperature ADC and AE were both higher in adults than in juveniles. Sprint speed increased with increase in body temperature at lower body temperatures, but decreased at higher body temperatures. At each test temperature adults ran faster than did juveniles, and the range of body temperatures where lizards maintained 90% of maximum speed differed between adults (27-34 degrees C) and juveniles (29-37 degrees C). Optimal temperatures and thermal sensitivities differed between food assimilation and sprint speed. Our results not only show strong patterns of ontogenetic variation in thermal tolerance, selected body temperature and thermal dependence of food assimilation and locomotor performance in E. brenchleyi, but also add support for the multiple optima hypothesis for the thermal dependence of behavioral and physiological variables in reptiles.