Biology Letters
You have access Research articles

An extreme cold event leads to community-wide convergence in lower temperature tolerance in a lizard community

Published:https://doi.org/10.1098/rsbl.2020.0625

    Abstract

    Extreme climate events are predicted to increase in frequency and severity due to contemporary climate change. Recent studies have documented the evolutionary impacts of extreme events on single species, but no studies have yet investigated whether such events can drive community-wide patterns of trait shifts. On 22 January 2020, subtropical south Florida experienced an extreme cold episode during which air temperatures dropped below the lower thermal limit of resident lizard populations. In the week immediately after the cold event, we documented decreased lower thermal limits (CTmin) of six co-occurring lizard species that vary widely in ecology, body size and thermal physiology. Although cold tolerance of these species differed significantly before the cold snap, lizards sampled immediately after had converged on the same new, lower limit of thermal tolerance. Here, we demonstrate that extreme climate events can drive substantial and synchronous community-wide trait changes and provide evidence that tropical and subtropical ectotherms—often characterized as unable to withstand rapid changes in climatic conditions—can endure climatic conditions that exceed their physiological limits. Future studies investigating the mechanisms driving these trait shifts will prove valuable in understanding the ability of ectotherm communities to mitigate climate change.

    1. Introduction

    Extreme climate events can be powerful agents of change [1]. Sudden shifts in temperature and rainfall regimes, or extreme weather events such as tropical cyclones or cold spells [27], can lead to dramatic changes in population-level traits [2,8]. However, the extent to which co-occurring species respond similarly to the same exceptional climate conditions is unclear. Do species-specific differences in ecology, morphology or physiology lead to differential phenotypic responses to shared extreme conditions, or can such events drive community-wide trait shifts, with responses comparable among species? While the evolutionary consequences of extreme climate events on populations of single species is a burgeoning area of research [1], no studies have yet explored the community-wide responses of multiple co-occurring species exposed to the same extreme climate event.

    Under contemporary climate change, temperature fluctuations are expected to increase in frequency and magnitude leading to exceptional heat waves and cold snaps [9]. Two non-mutually exclusive classes of biodiversity are considered highly vulnerable to increased temperature variation: (i) ectotherms, those species without internal physiological temperature regulation [10], and (ii) tropical species, which typically have narrow thermal tolerances owing to low climate variability [11,12]. Tropical lizards, therefore, present a valuable model to understand responses of thermally vulnerable species to rapid changes in climate.

    On 22 January 2020, south Florida experienced an extreme short-term cold event, with temperatures dropping overnight to 4.4°C—the lowest recorded temperature in the previous decade. Prior to this severe cold snap, we had measured the lower limits of thermal tolerance (critical thermal minima, CTmin) of six co-occurring tropical and subtropical lizard species in Miami, south Florida, USA [13]. Previous studies have suggested that exposure to exceptionally high or low temperatures can lead to changes in physiological traits associated with thermal biology in lizards, including CTmin [1416]. It is possible that such changes could result from adaptation via natural selection [15,16]—only those individuals in the population able to tolerate the extreme conditions survive—or acclimation via phenotypic plasticity [17], in which individual physiologies may be labile and shift in concert with abiotic conditions.

    In the week following the extreme cold event, we tested whether the surviving lizards were able to tolerate colder temperatures than those measured previously. We also investigated whether post-cold snap thermal tolerance limits were comparable across species despite substantial interspecific differences in thermal physiology, body size and ecology. Finally, we re-sampled thermal limits 10 weeks after the cold event to test whether differences in thermal tolerance observed directly after the event represented a short-term plastic response to cold temperature exposure, or if the shifts remained consistent, more indicative of an extreme selection event.

    2. Methods

    (a) Species sampling and field collection

    We sampled the lower thermal tolerance (CTmin) of six lizard species in Miami, south Florida, USA, between 21 October and 29 November 2016 [13]. Following the extreme cold event on 22 January 2020, we re-sampled all species from 29 January to 6 February 2020. The lizard species in this community are ecologically (figure 1b) and morphologically diverse (electronic supplementary material, table S1) and live 1–2 years in nature ([18,19]; J. T. Stroud 2015–2020, unpublished data). Hemidactylus mabouia were not sampled in the second post-cold snap sampling session owing to nocturnal curfews imposed in the City of Miami in response to the spread of SARS-CoV-2 (COVID-19 coronavirus; in place from 26 March 2020).

    Figure 1.

    Figure 1. Following an extreme cold event in south Florida, air temperatures (a) hit a decade-low of 4.4°C and six lizard species experienced short-term exposure to temperatures below their lower thermal tolerance limits. (b) Representation of the lizard community and their ecological niches drawn to scale with respect to relative size differences between species and not tree size. (c) Comparing lower thermal limits (CTmin) of species before (pre) and after (post) the extreme cold event, we observed shifts towards tolerance of lower temperatures. New, lower thermal limits remained 10 weeks after the extreme cold event. Horizontal blue-dashed line at 7.2°C illustrates the immediate post-cold event community-wide mean CTmin: the most recent cold event experienced temperatures below this limit for a substantially longer period of time than previous events (circled in (a)). *p < 0.05; n.s., not significant.

    (b) Assessing lower thermal tolerance limits (CTmin)

    Lower critical thermal limits (CTmin) were measured as previously in [13]. In short, lizards were acclimated overnight to room temperature before all thermal trials. A digital thermocouple was placed in the cloaca of each lizard and internal body temperature was monitored during gradual cooling until a righting response was lost (see electronic supplementary material).

    (c) Statistical analyses and calculating selection differentials

    We tested for a significant difference in CTmin between pre- and post-event populations (hereafter treatment) using analysis of covariance (ANCOVA), with species, treatment (and their interaction) and sex as model factors, and log-transformed body length (snout–vent length, SVL), body mass and starting body temperature as covariates. As the interaction between treatment and species in the full term was significant, we ran independent models per species to obtain p-values, including treatment as a model factor and log-transformed SVL as a covariate. To test if the variance in CTmin was greater in populations before or after the cold event, we conducted upper one-tailed F-tests for each species independently. A community-wide p-value was calculated via a Fisher's test of combined probability. Selection differentials were calculated by standardizing pre-cold and (plus 1 week) post-cold event samples combined to have a zero mean and unit variance for each species independently. Post-cold event means were then subtracted from pre-cold event means to calculate independent standardized differentials for each species [20]. Finally, we then tested for convergence in CTmin by fitting full models independently on pre- and post-cold event lizards; species, treatment (and their interaction) and sex were model factors, with log-transformed body length, mass and starting body temperature as covariates.

    3. Results

    (a) Species shift to tolerate colder temperatures

    Following the extreme cold event in south Florida, we observed shifts to lower thermal tolerance (CTmin, i.e. increased tolerance of colder temperatures) in all species of lizard tested (table 1 and figure 1). Species-level decreases in CTmin ranged from 1.09°C (Anolis cristatellus) to 4.04°C (Basiliscus vittatus). Variances in CTmin were smaller for all species after the cold event, except B. vittatus (table 2). Standardized selection differentials ranged from 0.92 to 1.80 (table 2).

    Table 1. (a) Following an extreme cold climate event in Miami, south Florida USA, we observed shifts to lower thermal tolerance limits (CTmin) in six lizard species, which remained consistent 10 weeks after the cold episode. Mean CTmin values per species are shown (±1 s.e.); sample sizes are given in parentheses. (b) Full community model indicates that CTmin of all species decreased following the cold event, including body size, mass, sex and starting temperature of CTmin trial as covariates. (c) Final model drops all non-significant covariates. **p < 0.01, ***p < 0.001.

    pre-cold event CTmin (°C) 1 week post-cold event CTmin (°C) F-value (d.f.) p-value 10 weeks post-cold event CTmin (°C) F-value (d.f.) p-value ecological niche: structural microhabitat
    (a) pre-/post-cold event differences in CTmin
    Anolis carolinensis 9.75 ± 0.4 (12) 7.43 ± 0.3 (11) 18.63 (1,21) <0.001*** 6.92 ± 0.2 (11) 32.96 (1,21) <0.001*** canopy branches, leaves, trunks
    Anolis cristatellus 8.04 ± 0.3 (10) 6.95 ± 0.2 (13) 8.85 (1,21) 0.007** 6.43 ± 0.2 (16) 19.05 (1,24) <0.001*** low tree trunks and ground
    Anolis distichus 9.6 ± 0.5 (10) 7.15 ± 0.2 (12) 19.96 (1,20) <0.001*** 6.75 ± 0.4 (11) 42.41 (1,19) <0.001*** broad tree trunks
    Anolis sagrei 9.05 ± 0.3 (12) 6.91 ± 0.5 (10) 14.17 (1,20) 0.001** 7.25 ± 0.3 (15) 23.18 (1,25) <0.001*** low tree trunks and ground
    Basiliscus vittatus 11.29 ± 0.3 (11) 7.25 ± 0.2 (11) 111.20 (1,20) <0.001*** 6.81 ± 0.4 (10) 85.84 (1,19) <0.001*** ground and riparian areas
    Hemidactylus mabouia 8.57 ± 0.5 (6) 7.38 ± 0.5 (6) 3.03 (1,10) 0.113 trees and branches; nocturnal
    (b) full community: CTmin shifts (full model)
    cold event 20.79 (1,88) <0.001***
    species 5.05 (5,88) <0.001***
    log(body size, SVL) 1.06 (1,88) 0.307
    log(mass) 0.78 (1,88) 0.379
    sex 1.33 (3,88) 0.269
    log(starting temp.°C) 1.52 (1,88) 0.221
    cold event × species 3.49 (5,88) 0.006**
    (c) full community: CTmin shifts (reduced model)
    cold event 23.27 (1,112) <0.001***
    species 9.77 (5,112) <0.001***
    cold event × species 4.33 (5,112) 0.001**

    Table 2. Minimum and maximum lower thermal limit (CTmin) values (range in parentheses), upper one-tailed F-tests for equality of variance and standardized selection differentials for observed shifts in CTmin in six lizard species sampled before and one week after an extreme cold event in Miami, south Florida, USA. The F-value represents an F : 1 ratio of variances between pre-cold event versus post-cold event samples. Community-wide shift in CTmin variance, p = 0.047 (Fisher's combined probability test).

    species pre-event CTmin range (°C) post-event CTmin range (°C) F p-value selection differential, s
    Anolis carolinensis 7.0–12.0 (5.0) 5.4–9.1 (3.7) 2.16 0.118 1.34
    Anolis cristatellus 6.0–9.4 (3.4) 5.2–8.0 (2.8) 1.33 0.315 1.07
    Anolis distichus 7.1–12.1 (5.0) 5.7–8.6 (2.9) 4.67 0.010 1.39
    Anolis sagrei 7.8–10.7 (2.9) 4.9–9.3 (4.4) 0.38 0.935 1.26
    Basiliscus vittatus 10.0–13.1 (3.1) 5.9–8.5 (2.6) 1.81 0.177 1.80
    Hemidactylus mabouia 6.8–9.7 (2.9) 6.3–9.2 (2.9) 1.27 0.401 0.92

    (b) Community-wide convergence on lower thermal limits

    Prior to the cold event, there was a significant difference in CTmin between species (F5,55 = 8.14, p < 0.001; figure 1). Interspecific differences also existed when excluding B. vittatus, which had the highest CTmin (F4,45 = 2.98, p = 0.029). By contrast, there were no interspecific differences in CTmin after the January 2020 cold event (F5,57 = 0.45, p = 0.813). Full model output is available in the electronic supplementary material.

    (c) Lower thermal limits maintained after the extreme cold event

    Re-sampling ca 10 weeks after the cold event revealed that CTmin did not differ from those observed immediately after the cold event (plus 1 week versus plus 10 weeks post-cold event: F1,110 = 1.42, p = 0.236), but significantly differed from those prior to the cold event (pre-cold event versus plus 10 weeks post-cold event: F1,108 = 36.67, p < 0.001; individual species differences in table 1). There also remained no interspecific differences in CTmin after 10 weeks (F4,56 = 1.90, p = 0.123; community-wide CTmin average = 6.28 ± 0.12°C (s.e.), mean CTmin range per species from 6.43 to 7.25°C).

    4. Discussion

    Extreme climate events such as unusual cold snaps, heat waves, droughts or exceptional rainfall are expected to increase in frequency and severity due to human-mediated climate change [2124]. Such atypical conditions can drastically affect populations and may drive long-lasting and large-scale effects on species and communities [3,25]. Following an extreme cold episode in subtropical south Florida, we observed shifts to increased tolerance of lower temperatures in six species of ecologically, physiologically and morphologically diverse lizards, resulting in community-wide convergence of lower thermal limits.

    Two non-exclusive mechanisms may have driven this observed community-wide shift in thermal tolerance: directional selection and physiological plasticity. A previous study demonstrated the role selection may play in driving the evolution of lizard cold tolerance [16]. During the winter of 2013–14, the south-central United States—north-west of south Florida—experienced record-setting extreme cold weather during the so-called ‘polar vortex.’ Similar to our observations, populations of the American green anole (Anolis carolinensis) had significantly lower thermal tolerance limits (CTmin) immediately after this extreme cold event [16]. Resampling the same A. carolinensis populations five months later revealed that the lower thermal tolerance limits of the populations following the cold event were maintained [16]. In concert with evidence of evolutionary change in genes relevant to thermal physiology, the maintenance of lower thermal limits suggested a plastic response was unlikely and so the observed genetic and phenotypic trait shifts most likely reflected the operation of natural selection [25].

    Here, we also observed that lower thermal limits recorded immediately after the cold event remain consistent over two months later in all species tested, suggesting there may have been community-wide directional selection towards lower thermal tolerance. Examination of the pre- and post-cold snap distributions and selection differentials suggests that, at least for some species, mortality would have had to be very high for directional selection to explain these patterns (table 2). While local reports of both perished and cold-stunned lizards following the cold event are well documented (see electronic supplementary material; albeit, all larger Iguana spp.), it is difficult to determine from our field observations whether sufficient mortality occurred for such extreme selection to alone explain these patterns. Another signature of selection would be a decrease in trait variance as only a subset of the initial population would have survived [26]. While we observed less variation in CTmin in post- versus pre-cold event sampling in all species except Anolis sagrei (table 2), we only observed a statistically significant decrease in one species (Anolis distichus; table 2). It is unclear, therefore, whether directional selection favouring the survival of more cold-tolerant individuals is responsible for the replicated CTmin shifts we observed across this community.

    Thermal trait shifts can also arise through phenotypic plasticity (i.e. acclimation; [17,2729]). For example, seasonal acclimation can lead to lower thermal limits during colder winter temperatures relative to hotter summers in some ectotherms (insects [30], freshwater and marine invertebrates [31], freshwater fish [32]), including temperate lizards [3336]. Tropical lizards, however, have a much weaker ability to acclimate to colder temperature regimes [3739]. In Miami specifically, species which exhibit seasonal variation in CTmin also acclimate during experimental exposure to weeks-long colder ambient temperatures (e.g. A. carolinensis and A. sagrei; [40,41]), while those with no seasonal variation do not (e.g. A. cristatellus; [14,40,42]). While some seasonal variation may exist, it remains unclear whether any tropical lizards are able to acclimate in response to short-term, rapid and extreme dips to low temperatures, as observed here. Plastic responses to such ‘cold shocks' have been observed in a range of insects [43,44], although the physiological mechanisms which underlie these responses remain unknown [30].

    It is also possible that directional selection and physiological plasticity have both operated in our study. If species within this lizard community varied in their acclimation ability, as indicated by previous studies [14,40,42], then it is possible that plasticity may underlie extreme shifts in trait distributions of CTmin in some species, whereas selection is more plausible for others. For example, the post-cold event CTmin levels of brown basilisks (B. vittatus) were substantially below those of any individual measured prior to the cold event. Such a pattern suggests that it was highly unlikely that natural selection was responsible, and plasticity may have played a greater role than in other species where selective mortality of existing variation could explain the shift (e.g. Anolis spp.). We also observed high interspecific variation in CTmin prior to the cold event, as well as differences in habitat use and body size (figure 1), suggesting that species may have been expected to respond differently to an extreme event. However, we observed convergence among all species on a new, lower CTmin. While differences in ecology and behaviour translate into different thermal ecologies of diurnal species during the day [45], species differ little in nocturnal behaviour: all species in this community sleep perched among above-ground vegetation ([18,46,47]; figure 2). As such, extreme cold temperatures which occurred in the middle of the night were likely experienced similarly across species.

    Figure 2.

    Figure 2. All diurnal lizards in this study sleep exposed on above-ground vegetation during periods of nocturnal inactivity. For example, here are typical sleeping site locations of Anolis sagrei in Miami, FL. Photo: J. T. Stroud, taken at 21.30 on 15 October 2019.

    Here, we observed that six co-occurring species of lizards converged on a new, lower limit of thermal tolerance following an extreme cold climate event. Multi-species field studies of thermal physiology remain rare, yet are likely to provide valuable insights into both our forecasting of the community-wide effects of contemporary climate change, and in developing our understanding of how climate and physiology interact in community assembly and associated dynamics. It is clear that extreme climate events pose particular risks to those species or life stages with limited ability to regulate internal body temperatures (e.g. ectotherms; [4851]) or those adapted to a narrow range of climatic conditions (e.g. thermal niche specialists; [12,52]); future physiological, genetic and demographic studies are needed to clarify the mechanisms that underlie both the immediate and long-term effects on the physiology of extreme climate events on thermally vulnerable species.

    Ethics

    Research was conducted under Florida Fish and Wildlife permit LSSC 16-00013 and University of Miami IACUC protocol 17-061.

    Data accessibility

    All raw data and code are included as electronic supplementary material.

    Authors' contributions

    J.T.S. conceived and designed the study. J.T.S., C.C.M. and W.B. acquired lizards and thermal trait data. J.T.S. analysed the data. All authors contributed to the interpretation of the data. J.T.S. drafted the manuscript and all authors provided critical revisions. All authors approved the final version of the manuscript and agreed to be held accountable for the content therein.

    Competing interests

    We declare we have no competing interests.

    Funding

    We received no funding for this study.

    Acknowledgements

    We would like to thank S. Clements, H. Howell, A. Messerman, C. Searcy and L. Stemle for field assistance. We would like to thank Sean Tomlinson and two anonymous reviewers for comments that improved the manuscript.

    Footnotes

    Lead contact.

    Joint first authors.

    Electronic supplementary material is available online at https://doi.org/10.6084/m9.figshare.c.5172233.

    Published by the Royal Society. All rights reserved.

    References