Colonization history affects heating rates of invasive cane toads

Amphibians in hot climates may be able to avoid high temperatures by controlling their rates of heating. In northern Australia, invasive cane toads (Rhinella marina) experience hot dry conditions in newly-colonized (western) sites but milder conditions in longer-occupied (eastern) sites. Under standardized conditions, toads from western sites heated less rapidly than did conspecifics from an eastern site. The availability of free water slowed heating rates of eastern but not western toads. Thus, the colonization of climatically extreme sites has been accompanied by a rapid shift in the toads’ ability to remain cool under hot conditions, even when free water is not available.


Methods
We hand-collected adult toads (n = 8 individuals per site) from four sites across the toads' tropical range within www.nature.com/scientificreports/ toad mass (after gently squeezing the animal in a standardized manner to induce it to empty its bladder) and snout-vent length (SVL) immediately before conducting the trials. Toads were not fed for three days prior to experiments, to ensure they would not defecate during the experiment and minimize variability in mass due to stomach contents. Toads from all four populations were housed in a room kept at 18 °C, then moved concurrently to a temperature-controlled room set at 37 °C. All toads were in separate containers (ventilated plastic boxes of 1-L capacity), half of which had dry paper towel as substrate whereas the other half had 40 mL of water, enough to keep the ventral portion of the body moist but not the rest of the body.
We measured toad body temperatures at the beginning of the trial, and after 20 min and 40 min, using an infrared thermometer (Digitech QM7215) held ~ 10 cm from the toad's dorsal surface. At the beginning and end of the experiment we measured internal temperatures with a cloacal probe (Digitech QM7215 with probe attachment), to check that our measurements of external body temperature offer robust estimates of internal temperature also. Cloacal temperatures were taken within 10 s of each toad's removal from the container. After a trial, toads were kept at a temperature of 25 °C, allowed to fully hydrate and monitored for wellbeing during recovery. No adverse effects of the trials were evident. We used mixed model repeated measures analysis to identify factors affecting body temperatures of cane toads during the 40-min heating trials. Sex and body mass were used as covariates in the analysis with climate at each collection site (# consecutive months per year with average maximum temperature > 30 °C and with < 25 mm rainfall (see above)) and wet/dry treatment as factors and time (0, 20, 40 min) as the repeated measure. Individual ID nested within population was used as a random effect. Our measure of climate characteristics at each site (0, 3, 5, or 6 hot dry months per year) was analysed as a continuous variable. However, for graphical purposes, our figures present data for each site depicted as a categorical variable (by location name). The data conformed to assumptions of normality and variance homogeneity and analyses were run using JMP Pro 14.0 (SAS Institute, Cary, NC).
The statistical analysis revealed two significant two-way interactions involving climatic conditions of the sites from which toads were collected. The first interaction term indicated that although all toads commenced trials at the same temperatures, individuals from the population with the coolest moistest climate (Townsville) heated faster than did toads from other populations (climate × time interaction P < 0.015, Fig. 3). Second, the effect of water treatment on body temperature depended on the climate from which the toads originated. When access to water was removed, toads from the coolest moistest climate (Townsville) exhibited higher body temperatures than did conspecifics from hotter drier locations (climate × water interaction P < 0.02, Tables 1, 2, Fig. 4).

Discussion
The heating rates of toads largely conformed to our predictions. Toads from a population that experiences a relatively brief period of hot dry conditions each year (Townsville) heated faster than did conspecifics from locations that experienced hotter drier conditions for extended periods, and this disparity was greatest if the toads did not have access to free water (Figs. 2, 4). Those patterns are consistent with the hypothesis that anurans can exert some level of control over heating rates, and that this ability is most developed in animals from climates where hot and dry conditions imperil organismal viability. The disparity in heating rates between Townsville toads with versus without access to free water suggests that if water is available, these animals use it to increase evaporation rates from the skin surface and hence, slow down their rate of heating. Interestingly, however, no such disparity was evident in conspecifics from hotter drier climates (Figs. 2, 4). In the seasonally dry landscapes of northwestern Australia, water reserves are critical to toad survival 12,13 and thus, toads in hot dry regions may benefit from an ability to curtail heating rates without expending water in the process.
The challenges of finding cool moist diurnal retreat sites may be greater for individuals at the expanding range edge of a biological invasion, than it is for philopatric individuals in long-established populations. Living in the same area for a long period allows an organism to learn the location of suitable shelter-sites, whereas a recentlyarrived organism at an invasion front has no such opportunity. As a result, individuals that are at the forefront of an invasion (or are translocated) may be forced to utilize retreats that provide less-than-optimal thermal regimes. Radiotelemetric studies of cane toads have shown that individuals encountering novel habitats tend to shelter in more superficial (exposed) sites 14 , and hence exhibit higher body temperatures than do philopatric  www.nature.com/scientificreports/ conspecifics 15 . That situation may exacerbate the thermal and hydric challenges that a colonizing toad encounters as the invasion front moves into hot dry climates (as in north-western Australia). Although geographic variation in the seasonal duration of hot dry conditions across our collection sites was significantly correlated with patterns of variation in body temperatures during the heating experiment, our data are not sufficient to demonstrate causation. Other factors also vary among toads from the four sites that we studied, including physiology 19 , behavioural syndromes 20 and morphology 21 . Thus, the geographic pattern that we found in rates of heating might be influenced by such characteristics, rather than (or as well as) by climatic  www.nature.com/scientificreports/ variation. Future work could usefully attempt to tease apart the roles of such factors, by employing larger sample sizes and looking for other correlates of rates of heating. It would also be instructive to examine commongarden-raised offspring of field-collected toads, to tease apart the extent to which rates of heating are affected by phenotypically plastic responses to the conditions experienced during an individual's life 20 .
Our data do not reveal the proximate mechanisms by which toads control rates of heating. Because the small size of the experimental enclosures gave minimal opportunities to manipulate heating rates by behavioural adjustments, morphological and/or physiological mechanisms are likely to be involved. Our measurements of skin thickness have revealed strong geographic heterogeneity in this trait, but with no close correspondence to the geographic variation in heating rates documented in the current paper 16 . Physiological mechanisms may be more important. For example, decreases in heartrate or shifts in blood flow (away from peripheral vessels) might reduce rates of heat uptake 17 . Future work could usefully explore the mechanisms by which toads can control their own rates of heating; and also, whether the geographic divergence in heating rate, and in the effect of water availability Table 2. Parameter estimates and standard errors from a mixed model repeated measures analysis on factors affecting surface temperatures of cane toads. Sex and body mass were used as covariates in the analysis which focused on interactive effects of climate at each capture site (# consecutive months per year with average maximum temperature > 30 °C and < 25 mm rainfall), treatment (presence or absence of water) and time (0, 20, 40 min). Individual ID nested within population was used as a random effect.  Figure 4. The interactive effect of climate at sites of collection (# consecutive months > 30 °C and < 25 mm rain) and presence of water during the trial, on body temperatures of cane toads during 40-min heating trials. When trials were conducted without access to water, toads from the population that experiences the coolest and wettest conditions (Townsville, 0 consecutive hot dry months per year) had higher average body temperatures than did toads from populations exposed to hotter drier climates. Data on toad temperatures are averaged across readings taken at 0, 20 and 40 min.
Scientific RepoRtS | (2020) 10:12553 | https://doi.org/10.1038/s41598-020-69529-3 www.nature.com/scientificreports/ on heating rate, are driven by adaptation (via genetic and/or epigenetic changes) or phenotypic plasticity (induced by exposure to harsh conditions). Despite the brief timespan of the toads' colonization of northwestern Australia (decades only), a complex mix of heritability and plasticity underlie other geographically-divergent traits in this species [18][19][20][21] . Although mechanisms remain unknown, the capacity of cane toads in hot dry climates to reduce heating rates, even if they have no access to water, may well have contributed to the remarkable success of this invasive anuran in the harsh climatic conditions of the Australian outback.