Ontogenetic reduction in thermal tolerance is not alleviated by earlier developmental acclimation in Rana temporaria

Drivers of Intraspecific Variation in Thermal Traits and Their Importance for Resilience to Global Change in Amphibians

Intraspecific variation can be as great as variation across species, but the role of intraspecific variation in driving local and large-scale patterns is often overlooked, particularly in the field of thermal biology. In amphibians, which depend on environmental conditions and behavior to regulate body temperature, recognizing intraspecific thermal trait variation is essential to comprehensively understanding how global change impacts populations. Here, we examine the drivers of micro- and macrogeographical intraspecific thermal trait variation in amphibians. At the local scale, intraspecific variation can arise via changes in ontogeny, body size, and between the sexes, and developmental plasticity, acclimation, and maternal effects may modulate predictions of amphibian performance under future climate scenarios. At the macrogeographic scale, local adaptation in thermal traits may occur along latitudinal and elevational gradients, with seasonality and range-edge dynamics likely playing important roles in patterns that may impact future persistence. We also discuss the importance of considering disease as a factor affecting intraspecific variation in thermal traits and population resilience to climate change, given the impact of pathogens on thermal preferences and critical thermal limits of hosts. Finally, we make recommendations for future work in this area. Ultimately, our goal is to demonstrate why it is important for researchers to consider intraspecific variation to determine the resilience of amphibians to global change.

Thermal sensitivity of Rhinella arenarum tadpole at low concentrations of dimethoate pesticides

One of the main causes of contamination of aquatic environments, which affects biotic communities, is the use of pesticides in agricultural regions. Amphibians are considered good bio-indicators of aquatic pollution, because they are one of the most susceptible groups to pollution. Several studies suggest that both pollution and climate change produce synergistic effects in amphibians which amplify the toxicity afecting survival, and malformations with an increase in temperature. We studied the sensitivity of sublethal concentrations of dimethoate in Rhinella arenarum tadpoles on two fitness related thermal traits including locomotor swimming performance and thermal tolerance limits (CTmax = critical thermal maximum and CTmin = critical thermal minimum). The locomotor performance of R. arenarum tadpoles decreased with increasing sublethal dimethoate concentrations up to ∼60 % at intermediates dimethoate concentration. The tadpoles showed a tendency to decrease their tolerance to high temperatures (CTmax) with increasing dimethoate concentration around ∼0.5 °C, however no significant differences were found among treatments. Similarly, tadpoles showed decreases in their cold resistance (CTmin) with dimethoate concentrations, around 1 °C the high concentrations of dimethoate. The increase of atypical climatic events, such as heat waves may put R. arenarum tadpoles at greater risk when exposed to dimethoate. Our results show that the sublethal concentrations of the dimethoate pesticide may affect the fitness and survival of the larvae of R. arenarum in natural, and seminatural enviroments.

Photoperiod effects in a freshwater community: Amphibian larvae develop faster and zooplankton abundance increases under an early-season photoperiod

Organisms that shift their phenologies in response to global warming will experience novel photic environments, as photoperiod (daylength) continues to follow the same annual cycle. How different organisms respond to novel photoperiods could result in phenological mismatches and altered interspecific interactions. We conducted an outdoor mesocosm experiment exposing green frog (Rana clamitans) larvae, gray treefrog (Hyla versicolor) larvae, phytoplankton, periphyton, and zooplankton to a three-month shift in photoperiod: an early-season photoperiod (simulating April) and a late-season photoperiod (simulating July). We manipulated photoperiod by covering and uncovering tanks with clear or light-blocking lids to mimic realistic changes in daylength. We assessed amphibian life history traits and measured phytoplankton, periphyton, and zooplankton abundances. Green frog larvae and gray treefrog metamorphs were more developed under the early-season photoperiod. Gray treefrog total length was also reduced, but photoperiod did not affect green frog total length. Although phytoplankton and periphyton abundances were not affected by photoperiod, copepod nauplii were in greater abundance under the early-season photoperiod. Overall, this simplified aquatic community did not exhibit significant changes to structure when exposed to a three-month shift in photoperiod. Temperate amphibians that breed earlier in the year may develop faster, which may have long-term costs to post-metamorphic growth and performance. Asynchronous shifts in zooplankton abundances in response to altered photoperiods could subsequently affect freshwater community structure. While photoperiod has been shown to individually affect freshwater organisms, our study using replicated outdoor wetland communities shows that the comprehensive effects of photoperiod may be less important than other cues such as temperature and precipitation.

Developmental environments do not affect thermal physiological traits in reptiles: an experimental test and meta-analysis

On a global scale, organisms face significant challenges due to climate change and anthropogenic disturbance. In many ectotherms, developmental and physiological processes are sensitive to changes in temperature and resources. Developmental plasticity in thermal physiology may provide adaptive advantages to environmental extremes if early environmental conditions are predictive of late-life environments. Here, we conducted a laboratory experiment to test how developmental temperature and maternal resource investment influence thermal physiological traits (critical thermal maximum: CTmax and thermal preference: Tpref) in a common skink (Lampropholis delicata). We then compared our experimental findings more broadly across reptiles (snakes, lizards and turtles) using meta-analysis. In both our experimental study and meta-analysis, we did not find evidence that developmental environments influence CTmax or Tpref. Furthermore, the effects of developmental environments on thermal physiology did not vary by age, taxon or climate zone (temperate/tropical). Overall, the magnitude of developmental plasticity on thermal physiology appears to be limited across reptile taxa suggesting that behavioural or evolutionary processes may be more important. However, there is a paucity of information across most reptile taxa, and a broader focus on thermal performance curves themselves will be critical in understanding the impacts of changing thermal conditions on reptiles in the future.

Metamorphosis Imposes Variable Constraints on Genome Expansion through Effects on Development

Genome size varies ∼100,000-fold across eukaryotes and has long been hypothesized to be influenced by metamorphosis in animals. Transposable element accumulation has been identified as a major driver of increase, but the nature of constraints limiting the size of genomes has remained unclear, even as traits such as cell size and rate of development co-vary strongly with genome size. Salamanders, which possess diverse metamorphic and non-metamorphic life histories, join the lungfish in having the largest vertebrate genomes-3 to 40 times that of humans-as well as the largest range of variation in genome size. We tested 13 biologically-inspired hypotheses exploring how the form of metamorphosis imposes varying constraints on genome expansion in a broadly representative phylogeny containing 118 species of salamanders. We show that metamorphosis during which animals undergo the most extensive and synchronous remodeling imposes the most severe constraint against genome expansion, with the severity of constraint decreasing with reduced extent and synchronicity of remodeling. More generally, our work demonstrates the potential for broader interpretation of phylogenetic comparative analysis in exploring the balance of multiple evolutionary pressures shaping phenotypic evolution.

A comprehensive database of amphibian heat tolerance

Rising temperatures represent a significant threat to the survival of ectothermic animals. As such, upper thermal limits represent an important trait to assess the vulnerability of ectotherms to changing temperatures. For instance, one may use upper thermal limits to estimate current and future thermal safety margins (i.e., the proximity of upper thermal limits to experienced temperatures), use this trait together with other physiological traits in species distribution models, or investigate the plasticity and evolvability of these limits for buffering the impacts of changing temperatures. While datasets on thermal tolerance limits have been previously compiled, they sometimes report single estimates for a given species, do not present measures of data dispersion, and are biased towards certain parts of the globe. To overcome these limitations, we systematically searched the literature in seven languages to produce the most comprehensive dataset to date on amphibian upper thermal limits, spanning 3,095 estimates across 616 species. This resource will represent a useful tool to evaluate the vulnerability of amphibians, and ectotherms more generally, to changing temperatures.

Thermal tolerance and acclimation capacity in the European common frog (Rana temporaria) change throughout ontogeny

Phenotypic plasticity may allow ectotherms with complex life histories such as amphibians to cope with climate-driven changes in their environment. Plasticity in thermal tolerance (i.e., shifts of thermal limits via acclimation to higher temperatures) has been proposed as a mechanism to cope with warming and extreme thermal events. However, thermal tolerance and, hence, acclimation capacity, is known to vary with life stage. Using the common frog (Rana temporaria) as a model species, we measured the capacity to adjust lower (CT_min ) and upper (CT_max ) critical thermal limits at different acclimation temperatures. We calculated the acclimation response ratio as a metric to assess the stage-specific acclimation capacity at each of seven consecutive ontogenetic stages and tested whether acclimation capacity was influenced by body mass and/or age. We further examined how acclimation temperature, body mass, age, and ontogenetic stage influenced CT_min and CT_max . In the temperate population of R. temporaria that we studied, thermal tolerance and acclimation capacity were affected by the ontogenetic stage. However, acclimation capacity at both thermal limits was well below 100% at all life stages tested. The lowest and highest acclimation capacity in thermal limits was observed in young and late larvae, respectively. The relatively low acclimation capacity of young larvae highlights a clear risk of amphibian populations to ongoing climate change. Ignoring stage-specific differences in thermal physiology may drastically underestimate the climate vulnerability of species, which will hamper successful conservation actions.

A state-space approach to understand responses of organisms, populations and communities to multiple environmental drivers

Understanding the response of biotic systems to multiple environmental drivers is one of the major concerns in ecology. The most common approach in multiple driver research includes the classification of interactive responses into categories (antagonistic, synergistic). However, there are situations where the use of classification schemes limits our understanding or cannot be applied. Here, we introduce and explore an approach that allows us to better appreciate variability in responses to multiple drivers. We then apply it to a case, comparing effects of heatwaves on performance of a cold-adapted species and a warm-adapted competitor. The heatwaves had a negative effect on the native (but not on the exotic) species and the approach highlighted that the exotic species was less responsive to multivariate environmental variation than the native species. Overall, we show how the proposed approach can enhance our understanding of variation in responses due to different driver intensities, species, genotypes, ontogeny, life-phases or among spatial scales at any level of biological organization.

Giménez et al. explore a “state-space” approach (SSEA) to examine variation in effects of multiple environmental drivers on biological systems. They illustrate the SSEA with a case study where larvae of an exotic crab were less responsive to an experimental heatwave than those of a native species.

Predators like it hot: Thermal mismatch in a predator-prey system across an elevational tropical gradient

Climate change may have dramatic consequences for communities through both direct effects of peak temperatures upon individual species and through interspecific mismatches in thermal sensitivities of interacting organisms which mediate changes in interspecific interactions (i.e. predation). Despite this, there is a paucity of information on the patterns of spatial physiological sensitivity of interacting species (at both landscape and local scales) which could ultimately influence geographical variation in the effects of climate change on community processes. In order to assess where these impacts may occur, we first need to evaluate the spatial heterogeneity in the degree of mismatch in thermal tolerances between interacting organisms. We quantify the magnitude of interspecific mismatch in maximum (CT_max ) and minimum (CT_min ) thermal tolerances among a predator-prey system of dragonfly and anuran larvae in tropical montane (242-3,631 m) and habitat (ponds and streams) gradients. To compare thermal mismatches between predator and prey, we coined the parameters maximum and minimum predatory tolerance margins (PTM_max and PTM_min ), or difference in CT_max and CT_min of interacting organisms sampled across elevational and habitat gradients. Our analyses revealed that: (a) predators exhibit higher heat tolerances than prey (~4°C), a trend which remained stable across habitats and elevations. In contrast, we found no differences in minimum thermal tolerances between these groups. (b) Maximum and minimum thermal tolerances of both predators and prey decreased with elevation, but only maximum thermal tolerance varied across habitats, with pond species exhibiting higher heat tolerance than stream species. (c) Pond-dwelling organisms from low elevations (0-1,500 m a.s.l.) may be more susceptible to direct effects of warming than their highland counterparts because their maximum thermal tolerances are only slightly higher than their exposed maximum environmental temperatures. The greater relative thermal tolerance of dragonfly naiad predators may further increase the vulnerability of lowland tadpoles to warming due to potentially enhanced indirect effects of higher predation rates by more heat-tolerant dragonfly predators. However, further experimental work is required to establish the individual and population-level consequences of this thermal tolerance mismatch upon biotic interactions such as predator-prey. ​.

Low desiccation and thermal tolerance constrains a terrestrial amphibian to a rare and disappearing microclimate niche

Drier and hotter conditions caused by climate change threaten species that exist close to their physiological limits, as well as those with limited ability to move. Habitat specialists may also be particularly vulnerable if they have specific abiotic requirements. Here we assess whether thermal and hydric constraints can explain the highly restricted and declining distributions of the critically endangered terrestrial-breeding frog, Geocrinia alba. We also evaluate the species' vulnerability to climate change based on the similarity of current microclimatic conditions to their physiological limits. We found that G. alba had low thresholds of thermal and desiccation tolerance relative to other anuran species. The estimated thermal optimum (T_opt ) and critical thermal maxima (CT_max ) were 23.3°C and 29.6°C, respectively, and adult frogs had an absorption threshold (AT, the lowest water potential at which water can be absorbed from a substrate) of -50 kPa, the lowest recorded for an amphibian. Comparing environmental conditions and water loss in the field using agar models showed that riparian habitats where frogs occur provide a unique microclimate in the landscape, offering significantly lower desiccation risk during extreme summer conditions compared to immediately adjacent riparian and terrestrial habitats. Monitoring of microclimate conditions within occupied frog habitats over 2 years showed that in extreme dry and hot years the AT was exceeded at six of eight sites, and T_opt was exceeded at two of eight sites. Given their specific physiological limits, the apparent rarity of suitable microclimates and a regional drying-warming trend, we suggest that G. alba occupies a potentially disappearing niche and may be indicative of other habitat specialists that rely on ephemeral drainages. More broadly, this study highlights that desiccation thresholds may tightly constrain amphibian distributions and need to be considered along with thermal tolerance thresholds when predicting the impacts of climate change.

Metamorphosis in an Era of Increasing Climate Variability

Most animals have complex life cycles including metamorphosis or other discrete life stage transitions, during which individuals may be particularly vulnerable to environmental stressors. With climate change, individuals will be exposed to increasing thermal and hydrologic variability during metamorphosis, which may affect survival and performance through physiological, behavioral, and ecological mechanisms. Furthermore, because metamorphosis entails changes in traits and vital rates, it is likely to play an important role in how populations respond to increasing climate variability. To identify mechanisms underlying population responses and associated trait and life history evolution, we need new approaches to estimating changes in individual traits and performance throughout metamorphosis, and we need to integrate metamorphosis as an explicit life stage in analytical models.

Critical Thermal Limits Do Not Vary between Wild-caught and Captive-bred Tadpoles of Agalychnis spurrelli (Anura: Hylidae)

Captive-bred organisms are widely used in ecology, evolution and conservation research, especially in scenarios where natural populations are scarce or at risk of extinction. Yet, it is still unclear whether captivity may alter thermal tolerances, crucial traits to predict species resilience to global warming. Here, we study whether captive-bred tadpoles of the gliding treefrog (Agalychnis spurrelli) show different thermal tolerances than wild-caught individuals. Our results show that there are no differences between critical thermal limits (CTmax and CTmin) of captive-bred and wild-caught tadpoles exposed to three-day acclimatization at 20 °C. Therefore, we suggest that the use of captive-bred amphibians is valid and may be appropriate in experimental comparisons to thermal physiological studies of wild populations.

Egg incubation temperature does not influence adult heat tolerance in the lizard Anolis sagrei

Extreme heat events are becoming more common as a result of anthropogenic global change. Developmental plasticity in physiological thermal limits could help mitigate the consequences of thermal extremes, but data on the effects of early temperature exposure on thermal limits later in life are rare, especially for vertebrate ectotherms. We conducted an experiment that to our knowledge is the first to isolate the effect of egg (i.e. embryonic) thermal conditions on adult heat tolerance in a reptile. Eggs of the lizard Anolis sagrei were incubated under one of three fluctuating thermal regimes that mimicked natural nest environments and differed in mean and maximum temperatures. After emergence, all hatchlings were raised under common garden conditions until reproductive maturity, at which point heat tolerance was measured. Egg mortality was highest in the warmest treatment, and hatchlings from the warmest treatment tended to have greater mortality than those from the cooler treatments. Despite evidence that incubation temperatures were stressful, we found no evidence that incubation treatment influenced adult heat tolerance. Our results are consistent with a low capacity for organisms to increase their physiological heat tolerance via plasticity, and emphasize the importance of behavioural and evolutionary processes as mechanisms of resilience to extreme heat.

Integrating mechanistic and correlative niche models to unravel range-limiting processes in a temperate amphibian

Insights into the causal mechanisms that limit species distributions are likely to improve our ability to anticipate species range shifts in response to climate change. For species with complex life histories, a mechanistic understanding of how climate affects different lifecycle stages may be crucial for making accurate forecasts. Here, we use mechanistic niche modeling (NicheMapR) to derive "proximate" (mechanistic) variables for tadpole, juvenile, and adult Rana temporaria. We modeled the hydroperiod, and maximum and minimum temperatures of shallow (30 cm) ponds, as well as activity windows for juveniles and adults. We then used those ("proximate") variables in correlative ecological niche models (Maxent) to assess their role in limiting the species' current distribution, and to investigate the potential effects of climate change on R. temporaria across Europe. We further compared the results with a model based on commonly used macroclimatic ("distal") layers (i.e., bioclimatic layers from WorldClim). The maximum temperature of the warmest month (a macroclimatic variable) and maximum pond temperatures (a mechanistic variable) were the most important range-limiting factors, and maximum temperature thresholds were consistent with the observed upper thermal limit of R. temporaria tadpoles. We found that range shift forecasts in central Europe are far more pessimistic when using distal macroclimatic variables, compared to projections based on proximate mechanistic variables. However, both approaches predicted extensive decreases in climatic suitability in southern Europe, which harbors a significant fraction of the species' genetic diversity. We show how mechanistic modeling provides ways to depict gridded layers that directly reflect the microenvironments experienced by organisms at continental scales, and to reconstruct those predictors without extrapolation under novel future conditions. Furthermore, incorporating those predictors in correlative ecological niche models can help shed light on range-limiting processes, and can have substantial impacts on predictions of climate-induced range shifts.

Immediate and Persistent Effects of Temperature on Oxygen Consumption and Thermal Tolerance in Embryos and Larvae of the Baja California Chorus Frog, Pseudacris hypochondriaca

The developmental environment has significant immediate effects on phenotypes, but it may also persistently or permanently shape phenotypes across life history. This study examined how developmental temperature influenced embryonic and larval phenotypes of Baja California chorus frog (Pseudacris hypochondriaca), an abundant amphibian in southern California and northern Baja California. We collected egg clutches from native ponds in northern San Diego County within 24 h of fertilization, and clutches were separated and distributed between constant temperatures of 10, 15, 20, and 25°C for incubation. Oxygen consumption rate (V˙O₂), developmental stage, and embryo and yolk masses were measured throughout development. Time to 50% hatch, survival at 50% hatch, and hatch duration were determined. Development rate was strongly affected by temperature, with warmer temperatures reducing time to hatch and hatch duration. Survival to hatch was high across all temperatures, >90%. Mass-specific V˙O₂ of embryos either remained constant or increased throughout development, and by hatching energy demand was significantly increased at higher temperatures. There were limited temperature effects on growth, with embryo and yolk dry mass similar between temperatures throughout embryonic development. To examine long-term effects of embryonic temperature, we reared hatchlings from each temperature until onset of larval feeding. Once feeding, larvae were acclimated to 20 or 25°C (>2 weeks). Following acclimation to 20 or 25°C, we measured larval mass-specific V˙O₂ and critical thermal maximum (CTMax) at a common developmental stage (Gosner stages 32–36, “hindlimb toe differentiation”). Embryonic temperature had persistent effects on larval V˙O₂ and CTMax, with warmer temperatures generally resulting in similar or higher V˙O₂ and CTMax. This partially supported a “warmer is better” effect of embryonic incubation temperature. These results suggest that in a thermally robust amphibian species, temperature may program the phenotype during early development to construct traits in thermal tolerance and energy use that may persist. Overall, P. hypochondriaca displays a thermally robust phenotype, and it is possible that amphibians that possess a wider range of phenotypic plasticity will be relatively more successful mitigating effects of climate change.

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.