Animals can respond to extreme climates by behaviourally avoiding it or by physiologically coping with it. We understand behavioural and physiological thermoregulation, but water balance has largely been neglected. Climate change includes both global warming and changes in precipitation regimes, so improving our understanding of organismal water balance is increasingly urgent. We assessed the hydric physiology of US federally endangered blunt-nosed leopard lizards (Gambelia sila) by measuring cutaneous evaporative water loss (CEWL), plasma osmolality and body condition. Measurements were taken throughout their active season, the short period of year when these lizards can be found aboveground. Compared to a more mesic species, G. sila had low CEWL which is potentially desert-adaptive, and high plasma osmolality that could be indicative of dehydration. We hypothesized that throughout the G. sila active season, as their habitat got hotter and drier, G. sila would become more dehydrated and watertight. Instead, CEWL and plasma osmolality showed minimal change for females and non-linear change for males, which we hypothesize is connected to sex-specific reproductive behaviours and changes in food availability. We also measured thermoregulation and microhabitat use, expecting that more dehydrated lizards would have lower body temperature, poorer thermoregulatory accuracy and spend less time aboveground. However, we found no effect of CEWL, plasma osmolality or body condition on these thermal and behavioural metrics. Finally, G. sila spends considerable time belowground in burrows, and burrows may serve not only as essential thermal refugia but also hydric refugia.

Testing acclimation plasticity informs our understanding of organismal physiology and applies to conservation management amidst our rapidly changing climate. Although there is a wealth of research on the plasticity of thermal and hydric physiology in response to temperature acclimation, there is a comparative gap for research on acclimation to different hydric regimes, as well as the interaction between water and temperature. We sought to fill this gap by acclimating western fence lizards (Sceloporus occidentalis) to experimental climate conditions (crossed design of hot or cool, dry or humid) for 8 days, and measuring cutaneous evaporative water loss (CEWL), plasma osmolality, hematocrit and body mass before and after acclimation. CEWL changed plastically in response to the different climates, with lizards acclimated to hot humid conditions experiencing the greatest increase in CEWL. Change in CEWL among individuals was negatively related to treatment vapor pressure deficit and positively related to treatment water vapor pressure. Plasma osmolality, hematocrit and body mass all showed greater changes in response to temperature than to humidity or vapor pressure deficit. CEWL and plasma osmolality were positively related across treatment groups before acclimation and within treatment groups after acclimation, but the two variables showed different responses to acclimation, suggesting that they are interrelated but governed by different mechanisms. This study is among few that assess more than one metric of hydric physiology and that test the interactive effects of temperature and humidity. Such measurements will be essential for predictive models of activity and survival for animals under climate change.

The climate crisis necessitates predicting how organisms respond to changing environments, but this requires understanding the mechanisms underlying thermal tolerance. The Hierarchical Mechanisms of Thermal Limitation (HMTL) hypothesis proposes that respiratory capacity and marginal stability of proteins and membranes interact hierarchically to determine thermal performance and limits. An untested prediction of the HMTL hypothesis is that behavioral anapyrexia (i.e., reduced body temperature in hypoxia) is exacerbated when metabolic demand is high. We tested this prediction by manipulating western fence lizards' (Sceloporus occidentalis) metabolic demand and oxygen environment, then measuring selected body temperatures. Lizards with elevated metabolic demand selected lower body temperatures at higher oxygen concentrations than resting lizards, but this occurred only at oxygen concentrations <12% O-2, suggesting thermal limits are unaffected by naturally-occurring oxygen variation. Given our results and the ubiquity of behavioral anapyrexia, the HMTL hypothesis may generally explain how oxygen and temperature interactively affect reptile performance.

Accurately predicting the responses of organisms to novel or changing environments requires the development of ecologically-appropriate experimental methodology and process-based models. For ectotherms, thermal performance curves (TPCs) have provided a useful framework to describe how organismal performance is dependent on temperature. However, this approach often lacks a mechanistic underpinning, which limits our ability to use TPCs predictively. Furthermore, thermal dependence varies across traits, and performance is also limited by additional abiotic factors, such as oxygen availability. We test a central prediction of our recent Hierarchical Mechanisms of Thermal Limitation (HMTL) Hypothesis which proposes that natural hypoxia exposure will reduce maximal performance and cause the TPC for whole-organism performance to become more symmetrical. We quantified TPCs for two traits often used as fitness proxies, sprint speed and aerobic scope, in lizards under conditions of normoxia and high-elevation hypoxia. In line with the predictions of HMTL, anaerobically fuelled sprint speed was unaffected by acute hypoxia while the TPC for aerobic scope became shorter and more symmetrical. This change in TPC shape resulted from both the maximum aerobic scope and the optimal temperature for aerobic scope being reduced in hypoxia as predicted. Following these results, we present a mathematical framework, which we call Temperature-Oxygen Performance Surfaces, to quantify the interactive effects of temperature and oxygen on whole-organism performance in line with the HMTL hypothesis. This framework is transferrable across traits and levels of organization to allow predictions for how ectotherms will respond to novel combinations of temperature and other abiotic factors, providing a useful tool in a time of rapidly changing environmental conditions. Read the free Plain Language Summary for this article on the Journal blog.

Populations of insects can differ in how sensitive their development, growth, and performance are to environmental conditions such as temperature and daylength. The environmental sensitivity of development can alter phenology (seasonal timing) and ecology. Warming accelerates development of most populations. However, high-elevation and season-limited populations can exhibit developmental plasticity to either advance or prolong development depending on conditions. We examine how diurnal temperature variation and daylength interact to shape growth, development, and performance of several populations of the montane grasshopper, Melanoplus boulderensis, along an elevation gradient. We then compare these experimental results to observed patterns of development in the field. Although populations exhibited similar thermal sensitivities of development under long-day conditions, development of high-elevation populations was more sensitive to temperature under short-day conditions. This developmental plasticity resulted in rapid development of high elevation populations in short-day conditions with high temperature variability, consistent with their observed capacity for rapid development in the field when conditions are permissive early in the season. Notably, accelerated development generally did not decrease body size or alter body shape. Developmental conditions did not strongly influence thermal tolerance but altered the temperature dependence of performance in difficult-to-predict ways. In sum, the high-elevation and season-limited populations exhibited developmental plasticity that enables advancing or prolonging development consistent with field phenology. Our results suggest these patterns are driven by the thermal sensitivity of development increasing when days are short early in the season compared to when days are long later in the season. Developmental plasticity will shape phenological responses to climate change with potential implications for community and ecosystem structure.

Two themes emerging from the special issue "Beyond CTMAX and CTMIN : Advances in Studying the Thermal Limits of Reptiles and Amphibians" are: (1) the need to identify mechanisms that determine the shape of thermal performance curves and (2) how these curves can be best used predictively.

BACKGROUND:High-quality genomic resources facilitate investigations into behavioral ecology, morphological and physiological adaptations, and the evolution of genomic architecture. Lizards in the genus Sceloporus have a long history as important ecological, evolutionary, and physiological models, making them a valuable target for the development of genomic resources. FINDINGS:We present a high-quality chromosome-level reference genome assembly, SceUnd1.0 (using 10X Genomics Chromium, HiC, and Pacific Biosciences data), and tissue/developmental stage transcriptomes for the eastern fence lizard, Sceloporus undulatus. We performed synteny analysis with other snake and lizard assemblies to identify broad patterns of chromosome evolution including the fusion of micro- and macrochromosomes. We also used this new assembly to provide improved reference-based genome assemblies for 34 additional Sceloporus species. Finally, we used RNAseq and whole-genome resequencing data to compare 3 assemblies, each representing an increased level of cost and effort: Supernova Assembly with data from 10X Genomics Chromium, HiRise Assembly that added data from HiC, and PBJelly Assembly that added data from Pacific Biosciences sequencing. We found that the Supernova Assembly contained the full genome and was a suitable reference for RNAseq and single-nucleotide polymorphism calling, but the chromosome-level scaffolds provided by the addition of HiC data allowed synteny and whole-genome association mapping analyses. The subsequent addition of PacBio data doubled the contig N50 but provided negligible gains in scaffold length. CONCLUSIONS:These new genomic resources provide valuable tools for advanced molecular analysis of an organism that has become a model in physiology and evolutionary ecology.

The stress phenotype is multivariate. Recent advances have broadened our understanding beyond characterizing the stress response in a single dimension. Simultaneously, the toolbox available to ecophysiologists has expanded greatly in recent years, allowing the measurement of multiple biomarkers from an individual at a single point in time. Yet these advances-in our conceptual understanding and available methodologies-have not yet been combined in a unifying multivariate statistical framework. Here, we offer a brief review of the multivariate stress phenotype and describe a general statistical approach for analysis using nonparametric multivariate analysis of variance with residual randomization in permutation procedures (RRPP) implemented using the "RRPP" package in R. We also provide an example illustrating the novel insights that can be gained from a holistic multivariate approach to stress and provide a tutorial for how we analyzed these data, including annotated R code and a guide to interpretation of outputs (Online Appendix 1). We hope that this statistical methodology will provide a quantitative framework facilitating the unification of our theoretical understanding and empirical observations into a quantitative, multivariate theory of stress.

Crested geckos (Correlophus ciliatus, formerly Rhacodactylus ciliatus) were rediscovered in New Caledonia 25 years ago and despite being common in the pet trade, there is no published information on their physiology. We measured thermoregulation (preferred body temperature, thermal set-point range, and voluntary limits) and performance (thermal performance curves [TPC] for 25 cm sprint speed and 1 m running speed) of adult and juvenile crested geckos in the laboratory to describe their thermal tolerances, differences among life stages, correlations between behavior and performance, and correlations with natural temperatures. Despite lacking special lighting or heating requirements in captivity, crested geckos displayed typical thermal biology for a lizard with no difference among life stages. They thermoregulated to a narrow set-point range (TSET , 24-28°C), that broadly overlaps natural air temperatures in New Caledonia, during activity. Somewhat surprisingly, the optimal temperature for performance (TOPT , 32°C) was substantially above preferred body temperatures and approximated the average maximum temperature voluntarily experienced (VTMAX , 33°C). Preferred body temperatures, by contrast, corresponded to the lower threshold temperature (Td ) where the TPC deviated from exponential, which we suggest is the temperature where performance is optimized after accounting for the costs of metabolic demand and overheating risk. Our results demonstrate that despite their lack of specific requirements when housed in human dwellings, crested geckos actively thermoregulate to temperatures that facilitate performance, and have thermal biology typical of other nocturnal or shade-dwelling species. Additionally, crested geckos appear at little risk of direct climate change-induced decline because increased temperatures should allow increased activity.

High-quality genomic resources facilitate population-level and species-level comparisons to answer questions about behavioral ecology, morphological and physiological adaptations, as well as the evolution of genomic architecture. Squamate reptiles (lizards and snakes) are particularly diverse in characteristics that have intrigued evolutionary biologists, but high-quality genomic resources for squamates are relatively sparse. Lizards in the genus Sceloporus have a long history as important ecological, evolutionary, and physiological models, making them a valuable target for the development of genomic resources. We present a high-quality chromosome-level reference genome assembly, SceUnd1.0, (utilizing 10X Genomics Chromium, HiC, and PacBio data) and tissue/developmental stage transcriptomes for the Eastern Fence Lizard, Sceloporus undulatus. We performed synteny analysis with other available squamate chromosome-level assemblies to identify broad patterns of chromosome evolution including the fusion of micro- and macrochromosomes in S. undulatus. Using this new S. undulatus genome assembly we conducted reference-based assemblies for 34 other Sceloporus species to improve draft nuclear genomes assemblies from 1% coverage to 43% coverage on average. Across these species, typically >90% of reads mapped for species within 20 million years divergence from S. undulatus, this dropped to 75% reads mapped for species at 35 million years divergence. Finally we use RNAseq and whole genome resequencing data to compare the three assemblies as references, each representing an increased level of sequencing, cost and assembly efforts: Supernova Assembly with data from10X Genomics Chromium library; HiRise Assembly that added data from HiC library; and PBJelly Assembly that added data from PacBio sequencing. We found that the Supernova Assembly contained the full genome and was a suitable reference for RNAseq, but the chromosome-level scaffolds provided by the addition of the HiC data allowed the reference to be used for other whole genome analysis, including synteny and whole genome association mapping analyses. The addition of PacBio data provided negligible gains. Overall, these new genomic resources provide valuable tools for advanced molecular analysis of an organism that has become a model in physiology and evolutionary ecology.

Research on the thermal ecology and physiology of free-living organisms is accelerating as scientists and managers recognize the urgency of the global biodiversity crisis brought on by climate change. As ectotherms, temperature fundamentally affects most aspects of the lives of amphibians and reptiles, making them excellent models for studying how animals are impacted by changing temperatures. As research on this group of organisms accelerates, it is essential to maintain consistent and optimal methodology so that results can be compared across groups and over time. This review addresses the utility of reptiles and amphibians as model organisms for thermal studies by reviewing the best practices for research on their thermal ecology and physiology, and by highlighting key studies that have advanced the field with new and improved methods. We end by presenting several areas where reptiles and amphibians show great promise for further advancing our understanding of how temperature relations between organisms and their environments are impacted by global climate change.

It is frequently hypothesized that animals employ a generalized "stress response," largely mediated by glucocorticoid (GC) hormones, such as corticosterone, to combat challenging environmental conditions. Under this hypothesis, diverse stressors are predicted to have concordant effects across biological levels of an organism. We tested the generalized stress response hypothesis in two complementary experiments with juvenile and adult male Eastern fence lizards (Sceloporus undulatus). In both experiments, animals were exposed to diverse, ecologically-relevant, acute stressors (high temperature or red imported fire ants, Solenopsis invicta) and we examined their responses at three biological levels: behavioral; physiological (endocrine [plasma corticosterone and blood glucose concentrations] and innate immunity [complement and natural antibodies]); and cellular responses (gene expression of a panel of five heat-shock proteins in blood and liver) at 30 or 90 min post stress initiation. In both experiments, we observed large differences in the cellular response to the two stressors, which contrasts the similar behavioral and endocrine responses. In the adult experiment for which we had innate immune data, the stressors affected immune function independently, and they were correlated with CORT in opposing directions. Taken together, these results challenge the concept of a generalized stress response. Rather, the stress response was context specific, especially at the cellular level. Such context-specificity might explain why attempts to link GC hormones with life history and fitness have proved difficult. Our results emphasize the need for indicators at multiple biological levels and whole-organism examinations of stress.