An interaction between climate change and infectious disease drove widespread amphibian declines

Climate change might drive species declines by altering species interactions, such as host–parasite interactions. However, few studies have combined experiments, field data, and historical climate records to provide evidence that an interaction between climate change and disease caused any host declines. A recently proposed hypothesis, the thermal mismatch hypothesis, could identify host species that are vulnerable to disease under climate change because it predicts that cool‐ and warm‐adapted hosts should be vulnerable to disease at unusually warm and cool temperatures, respectively. Here, we conduct experiments on Atelopus zeteki, a critically endangered, captively bred frog that prefers relatively cool temperatures, and show that frogs have high pathogen loads and high mortality rates only when exposed to a combination of the pathogenic chytrid fungus (Batrachochytrium dendrobatidis) and high temperatures, as predicted by the thermal mismatch hypothesis. Further, we tested various hypotheses to explain recent declines experienced by species in the amphibian genus Atelopus that are thought to be associated with B. dendrobatidis and reveal that these declines are best explained by the thermal mismatch hypothesis. As in our experiments, only the combination of rapid increases in temperature and infectious disease could account for the patterns of declines, especially in species adapted to relatively cool environments. After combining experiments on declining hosts with spatiotemporal patterns in the field, our findings are consistent with the hypothesis that widespread species declines, including possible extinctions, have been driven by an interaction between increasing temperatures and infectious disease. Moreover, our findings suggest that hosts adapted to relatively cool conditions will be most vulnerable to the combination of increases in mean temperature and emerging infectious diseases.     

Impacts of thermal mismatches on chytrid fungus Batrachochytrium dendrobatidis prevalence are moderated by life stage,body size,elevation and latitude

Global climate change is increasing the frequency of unpredictable weather conditions; however, it remains unclear how species‐level and geographic factors, including body size and latitude, moderate impacts of unusually warm or cool temperatures on disease. Because larger and lower‐latitude hosts generally have slower acclimation times than smaller and higher‐latitude hosts, we hypothesised that their disease susceptibility increases under ‘thermal mismatches’ or differences between baseline climate and the temperature during surveying for disease. Here, we examined how thermal mismatches interact with body size, life stage, habitat, latitude, elevation, phylogeny and International Union for Conservation of Nature (IUCN) conservation status to predict infection prevalence of the chytrid fungus Batrachochytrium dendrobatidis (Bd) in a global analysis of 32 291 amphibian hosts. As hypothesised, we found that the susceptibility of larger hosts and hosts from lower latitudes to Bd was influenced by thermal mismatches. Furthermore, hosts of conservation concern were more susceptible than others following thermal mismatches, suggesting that thermal mismatches might have contributed to recent amphibian declines.     

Infection risk decreases with increasing mismatch in host and pathogen environmental tolerances

The fungal pathogen Batrachochytrium dendrobatidis (Bd) has caused the greatest known wildlife pandemic, infecting over 500 amphibian species. It remains unclear why some host species decline from disease‐related mortality whereas others persist. We introduce a conceptual model that predicts that infection risk in ectotherms will decrease as the difference between host and pathogen environmental tolerances (i.e. tolerance mismatch) increases. We test this prediction using both local‐scale data from Costa Rica and global analyses of over 11 000 Bd infection assays. We find that infection prevalence decreases with increasing thermal tolerance mismatch and with increasing host tolerance of habitat modification. The relationship between environmental tolerance mismatches and Bd infection prevalence is generalisable across multiple amphibian families and spatial scales, and the magnitude of the tolerance mismatch effect depends on environmental context. These findings may help explain patterns of amphibian declines driven by a global wildlife pandemic.     

Host‐specific thermal profiles affect fitness of a widespread pathogen

Host behavior can interact with environmental context to influence outcomes of pathogen exposure and the impact of disease on species and populations. Determining whether the thermal behaviors of individual species influence susceptibility to disease can help enhance our ability to explain and predict how and when disease outbreaks are likely to occur. The widespread disease chytridiomycosis (caused by the fungal pathogen Batrachochytrium dendrobatidis, Bd) often has species‐specific impacts on amphibian communities; some host species are asymptomatic, whereas others experience mass mortalities and population extirpation. We determined whether the average natural thermal regimes experienced by sympatric frog species in nature, in and of themselves, can account for differences in vulnerability to disease. We did this by growing Bd under temperatures mimicking those experienced by frogs in the wild. At low and high elevations, the rainforest frogs Litoria nannotis, L. rheocola, and L. serrata maintained mean thermal regimes within the optimal range for pathogen growth (15–25°C). Thermal regimes for L. serrata, which has recovered from Bd‐related declines, resulted in slower pathogen growth than the cooler and less variable thermal regimes for the other two species, which have experienced more long‐lasting declines. For L. rheocola and L. serrata, pathogen growth was faster in thermal regimes corresponding to high elevations than in those corresponding to low elevations, where temperatures were warmer. For L. nannotis, which prefers moist and thermally stable microenvironments, pathogen growth was fastest for low‐elevation thermal regimes. All of the thermal regimes we tested resulted in pathogen growth rates equivalent to, or significantly faster than, rates expected from constant‐temperature experiments. The effects of host body temperature on Bd can explain many of the broad ecological patterns of population declines in our focal species, via direct effects on pathogen fitness. Understanding the functional response of pathogens to conditions experienced by the host is important for determining the ecological drivers of disease outbreaks.     

Mutualism fails when climate response differs between interacting species

Successful species interactions require that both partners share a similar cue. For many species, spring warming acts as a shared signal to synchronize mutualist behaviors. Spring flowering plants and the ants that disperse their seeds respond to warming temperatures so that ants forage when plants drop seeds. However, where warm‐adapted ants replace cold‐adapted ants, changes in this timing might leave early seeds stranded without a disperser. We investigate plant seed dispersal south and north of a distinct boundary between warm‐ and cold‐adapted ants to determine if changes in the ant species influence local plant dispersal. The warm‐adapted ants forage much later than the cold‐adapted ants, and so we first assess natural populations of early and late blooming plants. We then transplant these plants south and north of the ant boundary to test whether distinct ant climate requirements disrupt the ant–plant mutualism. Whereas the early blooming plant's inability to synchronize with the warm‐adapted ant leaves its populations clumped and patchy and its seedlings clustered around the parents in natural populations, when transplanted into the range of the cold‐adapted ant, effective seed dispersal recovers. In contrast, the mutualism persists for the later blooming plant regardless of location because it sets seed later in spring when both warm‐ and cold‐adapted ant species forage, resulting in effective seed dispersal. These results indicate that the climate response of species interactions, not just the species themselves, is integral in understanding ecological responses to a changing climate. Data linking phenological synchrony and dispersal are rare, and these results suggest a viable mechanism by which a species' range is limited more by biotic than abiotic interactions – despite the general assumption that biotic influences are buried within larger climate drivers. These results show that biotic partner can be as fundamental a niche requirement as abiotic resources.     

Too hot to die? The effects of vegetation shading on past,present, and future activity budgets of two diurnal skinks from arid Australia

Behavioral thermoregulation is an important mechanism allowing ectotherms to respond to thermal variations. Its efficiency might become imperative for securing activity budgets under future climate change. For diurnal lizards, thermal microhabitat variability appears to be of high importance, especially in hot deserts where vegetation is highly scattered and sensitive to climatic fluctuations. We investigated the effects of a shading gradient from vegetation on body temperatures and activity timing for two diurnal, terrestrial desert lizards, Ctenotus regius, and Morethia boulengeri, and analyzed their changes under past, present, and future climatic conditions. Both species’ body temperatures and activity timing strongly depended on the shading gradient provided by vegetation heterogeneity. At high temperatures, shaded locations provided cooling temperatures and increased diurnal activity. Conversely, bushes also buffered cold temperature by saving heat. According to future climate change scenarios, cooler microhabitats might become beneficial to warm‐adapted species, such as C. regius, by increasing the duration of daily activity. Contrarily, warmer microhabitats might become unsuitable for less warm‐adapted species such as M. boulengeri for which midsummers might result in a complete restriction of activity irrespective of vegetation. However, total annual activity would still increase provided that individuals would be able to shift their seasonal timing towards spring and autumn. Overall, we highlight the critical importance of thermoregulatory behavior to buffer temperatures and its dependence on vegetation heterogeneity. Whereas studies often neglect ecological processes when anticipating species’ responses to future climate change the strongest impact of a changing climate on terrestrial ectotherms in hot deserts is likely to be the loss of shaded microhabitats rather than the rise in temperature itself. We argue that conservation strategies aiming at addressing future climate changes should focus more on the cascading effects of vegetation rather than on shifts of species distributions predicted solely by climatic envelopes.     

Design‐ and model‐based recommendations for detecting and quantifying an amphibian pathogen in environmental samples

Accurate pathogen detection is essential for developing management strategies to address emerging infectious diseases, an increasingly prominent threat to wildlife. Sampling for free‐living pathogens outside of their hosts has benefits for inference and study efficiency, but is still uncommon. We used a laboratory experiment to evaluate the influences of pathogen concentration, water type, and qPCR inhibitors on the detection and quantification of Batrachochytrium dendrobatidis (Bd) using water filtration. We compared results pre‐ and post‐inhibitor removal, and assessed inferential differences when single versus multiple samples were collected across space or time. We found that qPCR inhibition influenced both Bd detection and quantification in natural water samples, resulting in biased inferences about Bd occurrence and abundance. Biases in occurrence could be mitigated by collecting multiple samples in space or time, but biases in Bd quantification were persistent. Differences in Bd concentration resulted in variation in detection probability, indicating that occupancy modeling could be used to explore factors influencing heterogeneity in Bd abundance among samples, sites, or over time. Our work will influence the design of studies involving amphibian disease dynamics and studies utilizing environmental DNA (eDNA) to understand species distributions.     

Upward ant distribution shift corresponds with minimum,not maximum,temperature tolerance

Rapid climate change may prompt species distribution shifts upward and poleward, but species movement in itself is not sufficient to establish climate causation. Other dynamics, such as disturbance history, may prompt species distribution shifts resembling those expected from rapid climate change. Links between species distributions, regional climate trends and physiological mechanism are needed to convincingly establish climate‐induced species shifts. We examine a 38‐year shift (1974–2012) in an elevation ecotone between two closely related ant species, Aphaenogaster picea and A. rudis. Even though A. picea and A. rudis are closely related with North American distributions that sometimes overlap, they also exhibit local‐ and regional‐scale differences in temperature requirements so that A. rudis is more southerly and inhabits lower elevations whereas A. picea is more northerly and inhabits high elevations. We find considerable movement by the warm‐habitat species upward in elevation between 1974 and 2012 with A. rudis, replacing the cold‐habitat species, A. picea, along the southern edge of the Appalachian Mountain chain in north Georgia, USA. Concomitant with the distribution shifts, regional mean and maximum temperatures remain steady (1974–2012), but minimum temperatures increase. We collect individuals from the study sites and subject them to thermal tolerance testing in a controlled setting and find that maximum and minimum temperature acclimatization occurs along the elevation gradient in both species, but A. rudis consistently becomes physiologically incapacitated at minimum and maximum temperatures 2 °C higher than A. picea. These results indicate that rising minimum temperatures allow A. rudis to move upward in elevation and displace A. picea. Given that Aphaenogaster ants are the dominant woodland seed dispersers in eastern deciduous forests, and that their thermal tolerances drive distinct differences in temperature‐cued synchrony with early blooming plants, these climate responses not only impact ant‐ant interactions, but might have wide implications for ant‐plant interactions.     

Elevated temperatures alter an ant‐aphid mutualism

Ant‐hemipteran mutualisms are keystone interactions that can be variously affected by warming: these mutualisms can be strengthened or weakened, or the species can transition to new mutualist partners. We examined the effects of elevated temperatures on an ant‐aphid mutualism in the subalpine zone of the Rocky Mountains in Colorado, USA. In this system, inflorescences of the host plant, Ligusticum porteri Coult. & Rose (Apiaceae), are colonized by the ant‐tended aphid Aphis asclepiadis Fitch or less frequently by the non‐ant tended aphid Cavariella aegopodii (Scopoli) (both Hemiptera: Aphididae). Using an 8‐year observational study, we tested for two key mechanisms by which ant‐hemipteran mutualisms may be altered by climate change: shifts in species identity and phenological mismatch. Whereas the aphid species colonizing the host plant is not changing in response to year‐to‐year variation in temperature, we found evidence that a phenological mismatch between ants and aphids could occur. In warmer years, colonization of host plant inflorescences by ants is decreased, whereas for A. asclepiadis aphids, host plant colonization is mostly responsive to date of snowmelt. We also experimentally established A. asclepiadis colonies on replicate host plants at ambient and elevated temperatures. Ant abundance did not differ between aphid colonies at ambient vs. elevated temperatures, but ants were less likely to engage in tending behaviors on aphid colonies at elevated temperatures. Sugar composition of aphid honeydew was also altered by experimental warming. Despite reduced tending by ants, aphid colonies at elevated temperatures had fewer intraguild predators. Altogether, our results suggest that higher temperatures may disrupt this ant‐aphid mutualism through both phenological mismatch and by altering benefits exchanged in the interaction.     

The seesaw effect of winter temperature change on the recruitment of cotton bollworms Helicoverpa armigera through mismatched phenology

Knowing how climate change affects the population dynamics of insect pests is critical for the future of integrated pest management. Rising winter temperatures from global warming can drive increases in outbreaks of some agricultural pests. In contrast, here we propose an alternative hypothesis that both extremely cold and warm winters can mismatch the timing between the eclosion of overwintering pests and the flowering of key host plants. As host plants normally need higher effective cumulative temperatures for flowering than insects need for eclosion, changes in flowering time will be less dramatic than changes in eclosion time, leading to a mismatch of phenology on either side of the optimal winter temperature. We term this the “seesaw effect.” Using a long‐term dataset of the Old World cotton bollworm Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) in northern China, we tested this seesaw hypothesis by running a generalized additive model for the effects of the third generation moth in the preceding year, the winter air temperature, the number of winter days below a critical temperature and cumulative precipitation during winter on the demography of the overwintering moth. Results confirmed the existence of the seesaw effect of winter temperature change on overwintering populations. Pest management should therefore consider the indirect effect of changing crop phenology (whether due to greenhouse cultivation or to climate change) on pest outbreaks. As arthropods from mid‐ and high latitudes are actually living in a cooler thermal environment than their physiological optimum in contrast to species from lower latitudes, the effects of rising winter temperatures on the population dynamics of arthropods in the different latitudinal zones should be considered separately. The seesaw effect makes it more difficult to predict the average long‐term population dynamics of insect pests at high latitudes due to the potential sharp changes in annual growth rates from fluctuating minimum winter temperatures.     

Improving methods to measure critical thermal limits in phloem‐feeding pest insects

The ability to accurately assess thermal tolerance in the laboratory without compromising ecological relevance is essential to predict the impacts of global climate change on phytophagous pest insects such as the phloem‐feeding aphids. One method to study thermal tolerance employs a temperature‐controlled column to measure critical thermal limits. However, assessments are commonly made with little relation to the natural environment of the study species. This study measured critical thermal minima (CT_min) for three cereal aphids – Sitobion avenae (Fabricius), Metopolophium dirhodum (Walker), and Rhopalosiphum padi (L.) (all Hemiptera: Aphididae) – in the absence and presence of host plant material to determine the best experimental design. Results revealed that CT_min measured in the presence of the host plant was significantly lower, suggesting that performing the measurement in the absence of the host plant could result in an underestimation of insect thermal tolerance. In addition, the study highlights the importance of understanding how an insect interacts with its environment, as this can reveal behavioural variation integral to differential survival at unfavourable temperatures.     

The capacity to cope with climate warming declines from temperate to tropical latitudes in two widely distributed Eucalyptus species

As rapid climate warming creates a mismatch between forest trees and their home environment, the ability of trees to cope with warming depends on their capacity to physiologically adjust to higher temperatures. In widespread species, individual trees in cooler home climates are hypothesized to more successfully acclimate to warming than their counterparts in warmer climates that may approach thermal limits. We tested this prediction with a climate‐shift experiment in widely distributed Eucalyptus tereticornis and E. grandis using provenances originating along a ~2500 km latitudinal transect (15.5–38.0°S) in eastern Australia. We grew 21 provenances in conditions approximating summer temperatures at seed origin and warmed temperatures (+3.5 °C) using a series of climate‐controlled glasshouse bays. The effects of +3.5 °C warming strongly depended on home climate. Cool‐origin provenances responded to warming through an increase in photosynthetic capacity and total leaf area, leading to enhanced growth of 20–60%. Warm‐origin provenances, however, responded to warming through a reduction in photosynthetic capacity and total leaf area, leading to reduced growth of approximately 10%. These results suggest that there is predictable intraspecific variation in the capacity of trees to respond to warming; cool‐origin taxa are likely to benefit from warming, while warm‐origin taxa may be negatively affected.     

Cold range edges of marine fishes track climate change better than warm edges

Species around the world are shifting their ranges in response to climate change. To make robust predictions about climate‐related colonizations and extinctions, it is vital to understand the dynamics of range edges. This study is among the first to examine annual dynamics of cold and warm range edges, as most global change studies average observational data over space or over time. We analyzed annual range edge dynamics of marine fishes—both at the individual species level and pooled into cold‐ and warm‐edge assemblages—in a multi‐decade time‐series of trawl surveys conducted on the Northeast US Shelf during a period of rapid warming. We tested whether cold edges show stronger evidence of climate tracking than warm edges (due to non‐climate processes or time lags at the warm edge; the biogeography hypothesis or extinction debt hypothesis), or whether they tracked temperature change equally (due to the influence of habitat suitability; the ecophysiology hypothesis). In addition to exploring correlations with regional temperature change, we calculated species‐ and assemblage‐specific sea bottom and sea surface temperature isotherms and used them to predict range edge position. Cold edges shifted further and tracked sea surface and bottom temperature isotherms to a greater degree than warm edges. Mixed‐effects models revealed that for a one‐degree latitude shift in isotherm position, cold edges shifted 0.47 degrees of latitude, and warm edges shifted only 0.28 degrees. Our results suggest that cold range edges are tracking climate change better than warm range edges, invalidating the ecophysiology hypothesis. We also found that even among highly mobile marine ectotherms in a global warming hotspot, few species are fully keeping pace with climate.     

Precipitation Constrains Amphibian Chytrid Fungus Infection Rates in a Terrestrial Frog Assemblage in Jamaica,West Indies

We model Batrachochytrium dendrobatidis (Bd) infection rates in Jamaican frogs—one of the most threatened amphibian fauna in the world. The majority of species we surveyed were terrestrial direct‐developing frogs or frogs that breed in tank bromeliads, rather than those that use permanent water bodies to breed. Thus, we were able to investigate the climatic correlates of Bd infection in a frog assemblage that does not rely on permanent water bodies. We sampled frogs for Bd across all of the major habitat types on the island, used machine learning algorithms to identify climatic variables that are correlated with infection rates, and extrapolated infection rates across the island. We compared the effectiveness of the machine learning algorithms for species distribution modeling in the context of our study, and found that infection rate rose quickly with precipitation in the driest month. Infection rates also increased with mean temperature in the warmest quarter until 22 °C, and remained relatively level thereafter. Both of these results are in accordance with previous studies of the physiology of Bd. Based on our environmental results, we suggest that frogs occupying high‐precipitation habitats with cool rainy‐season temperatures, though zcurrently experiencing low frequencies of infection, may experience an increase in infection rates as global warming increases temperatures in their habitat.     

Climate variation alters the synchrony of host–parasitoid interactions

Observed changes in mean temperature and increased frequency of extreme climate events have already impacted the distributions and phenologies of various organisms, including insects. Although some research has examined how parasitoids will respond to colder temperatures or experimental warming, we know relatively little about how increased variation in temperature and humidity could affect interactions between parasitoids and their hosts. Using a study system consisting of emerald ash borer (EAB), Agrilus planipennis, and its egg parasitoid Oobius agrili, we conducted environmentally controlled laboratory experiments to investigate how increased seasonal climate variation affected the synchrony of host–parasitoid interactions. We hypothesized that increased climate variation would lead to decreases in host and parasitoid survival, host fecundity, and percent parasitism (independent of host density), while also influencing percent diapause in parasitoids. EAB was reared in environmental chambers under four climate variation treatments (standard deviations in temperature of 1.24, 3.00, 3.60, and 4.79°C), while O. agrili experiments were conducted in the same environmental chambers using a 4 × 3 design (four climate variation treatments × 3 EAB egg densities). We found that EAB fecundity was negatively associated with temperature variation and that temperature variation altered the temporal egg laying distribution of EAB. Additionally, even moderate increases in temperature variation affected parasitoid emergence times, while decreasing percent parasitism and survival. Furthermore, percent diapause in parasitoids was positively associated with humidity variation. Our findings indicate that relatively small changes in the frequency and severity of extreme climate events have the potential to phenologically isolate emerging parasitoids from host eggs, which in the absence of alternative hosts could lead to localized extinctions. More broadly, these results indicate how climate change could affect various life history parameters in insects, and have implications for consumer–resource stability and biological control.     

Life‐history traits and physiological limits of the alpine fly Drosophila nigrosparsa (Diptera: Drosophilidae): A comparative study

Interspecific variation in life‐history traits and physiological limits can be linked to the environmental conditions species experience, including climatic conditions. As alpine environments are particularly vulnerable under climate change, we focus on the montane‐alpine fly Drosophila nigrosparsa. Here, we characterized some of its life‐history traits and physiological limits and compared these with those of other drosophilids, namely Drosophila hydei, Drosophila melanogaster, and Drosophila obscura. We assayed oviposition rate, longevity, productivity, development time, larval competitiveness, starvation resistance, and heat and cold tolerance. Compared with the other species assayed, D. nigrosparsa is less fecund, relatively long‐living, starvation susceptible, cold adapted, and surprisingly well heat adapted. These life‐history characteristics provide insights into invertebrate adaptations to alpine conditions which may evolve under ongoing climate change.     

Host‐associated genomic differentiation in congeneric butterflies: now you see it,now you do not

Ecotypic variation among populations may become associated with widespread genomic differentiation, but theory predicts that this should happen only under particular conditions of gene flow, selection and population size. In closely related species, we might expect the strength of host‐associated genomic differentiation (HAD) to be correlated with the degree of phenotypic differentiation in host‐adaptive traits. Using microsatellite and Amplified Fragment Length Polymorphism (AFLP) markers, and controlling for isolation by distance between populations, we sought HAD in two congeneric species of butterflies with different degrees of host plant specialization. Prior work on Euphydryas editha had shown strong interpopulation differentiation in host‐adapted traits, resulting in incipient reproductive isolation among host‐associated ecotypes. We show here that Euphydryas aurinia had much weaker host‐associated phenotypic differentiation. Contrary to our expectations, we detected HAD in Euphydryas aurinia, but not in E. editha. Even within an E. aurinia population that fed on both hosts, we found weak but significant sympatric HAD that persisted in samples taken 9 years apart. The finding of significantly stronger HAD in the system with less phenotypic differentiation may seem paradoxical. Our findings can be explained by multiple factors, ranging from differences in dispersal or effective population size, to spatial variation in genomic or phenotypic traits and to structure induced by past histories of host‐adapted populations. Other infrequently measured factors, such as differences in recombination rates, may also play a role. Our result adds to recent work as a further caution against assumptions of simple relationships between genomic and adaptive phenotypic differentiation.     

Temperature alters reproductive life history patterns in Batrachochytrium dendrobatidis,a lethal pathogen associated with the global loss of amphibians

Understanding how pathogens respond to changing environmental conditions is a central challenge in disease ecology. The environmentally sensitive fungal pathogen Batrachochytrium dendrobatidis (Bd), which causes the amphibian disease chytridiomycosis, has spread globally causing amphibian extirpations in a wide variety of climatic regions. To gain an in‐depth understanding of Bd's responses to temperature, we used an integrative approach, combining empirical laboratory experiments with mathematical modeling. First, we selected a single Bd isolate and serially propagated two lineages of the isolate for multiple generations in two stable thermal conditions: 4°C (cold‐adapted lineage) and 23°C (warm‐adapted lineage). We quantified the production of infectious zoospores (fecundity), the timing of zoospore release, and zoospore activity in reciprocal temperature transplant experiments in which both Bd lineages were grown in either high or low temperature conditions. We then developed population growth models for the Bd lineages under each set of temperature conditions. We found that Bd had lower population growth rates, but longer periods of zoospore activity in the low temperature treatment (4°C) compared to the high temperature treatment (23°C). This effect was more pronounced in Bd lineages that were propagated in the low temperature treatment (4°C), suggesting a shift in Bd's response to low temperature conditions. Our results provide novel insights into the mechanisms by which Bd can thrive in a wide variety of temperature conditions, potentially altering the dynamics of chytridiomycosis and thus, the propensity for Bd to cause amphibian population collapse. We also suggest that the adaptive responses of Bd to thermal conditions warrant further investigation, especially in the face of global climate change.