Body Size Is Related to Temperature Preference in Hyla chrysoscelis Tadpoles
Temperature plays a critical role in growth and development of anurans. Here, we investigated thermal preference of tadpoles from a wide-ranging anuran, Hyla chrysoscelis, by establishing an ex situ thermal gradient in a laboratory setting. We collected 64 tadpoles with a mean (±SD) body length of 12.2 ± 2.1 mm (range, 7–15 mm). Mean temperature preference of tadpoles ranged from 12.9 to 30.0°C, and we found a positive relationship between temperature preference and body length. Although body length is only one component of amphibian development, its apparent link with temperature preference is an exciting area for further research in this species as well as other species. Furthermore, given the broad distribution of H. chrysoscelis, there is an opportunity to investigate how thermal preferences could limit population edges and how populations may respond to shifts in environmental temperatures.ABSTRACT
Many organisms have adapted successfully to long-term temperature fluctuations, but when temperatures shift quickly, these same organisms often rely instead on behavioral responses as a means of acclimation. For aquatic ectotherms, conduction and convection with the surrounding water are the primary means by which they manipulate their body temperature (Reynolds and Casterlin, 1979). They are capable of moving through their environment to select areas in the water with the most optimal temperature (Hutchison and Dupré, 1992). However, temperature fluctuations tend to be rather minimal in aquatic environments because of water's high specific heat, thereby limiting the range of temperatures these organisms can realistically access (Hillman et al., 2009). This limitation in aquatic thermal range is important given that body temperatures for aquatic ectotherms rarely differ from the surrounding water by more than 2°C (Spotila et al., 1992). Understanding the temperatures that individuals select at different times in life can help illuminate their thermal needs at a particular age or developmental stage as well as how those needs may be impacted by rapidly shifting global temperatures. Improved knowledge of temperature selection shifts throughout ontogeny can also increase the accuracy and utility of species distribution models, as survivability can vary with body size (Riddell et al., 2018; Gouveia et al., 2019).
Temperature preference is of particular interest to species with complex life histories, such as amphibians and aquatic insects. Temperature preference for some anuran species can shift as tadpoles progress through development, with individuals in later developmental stages showing preference for warmer temperatures than those at earlier stages (Floyd, 1984; Dupré and Petranka, 1985). For example, both Rana (Lithobates) catesbeiana and Rhinella marina tadpoles exhibit a preference for higher temperatures at later developmental stages (Hutchison and Hill, 1978; Floyd, 1984; Dupré et al., 1986). Pseudacris triseriata tadpoles do not show distinct temperature preferences earlier in development, but from stage 36 onward, they show a preference for warmer temperatures (Dupré and Petranka, 1985). However, there are some species of anurans—for example, Rana boylii—that exhibit no difference in temperature selection at the various life-history stages (Catenazzi and Kupferberg, 2013).
Larval amphibians, in general, balance growth and metamorphosis, depending on a variety of environmental conditions, by using the release of particular hormones to stimulate and inhibit responses (Wilbur and Collins, 1973; Werner, 1986; Hillman et al., 2009). In addition, Rana arvalis tadpoles exposed to increased temperatures demonstrated shorter larval periods and larger masses at metamorphosis (Richter-Boix et al., 2015), providing evidence that temperature can interact with other drivers, such as hormones, to alter the rate of growth, metamorphosis, or both. In tadpoles, selection of lower temperatures typically leads to a slower rate of development (Smith-Gill and Berven, 1979), likely the result of decreased metabolism, whereas selection of higher temperatures has been shown to increase growth and developmental rate (Hutchison and Hill, 1978; Wu et al., 2007). The ability for tadpoles to influence their developmental rate is critical for survival, especially to escape stressors, such as predation and pond drying (Wilbur, 1997). In addition, swimming performance increases with temperature in Limnodynastes peronii tadpoles (Niehaus et al., 2011), so the selection for higher temperatures could also maximize the ability of tadpoles to escape predation with improved locomotor performance.
Here, we address the temperature preference of tadpoles of a broadly distributed anuran species, Hyla chrysoscelis (Cope, 1880) (Cope's Gray Treefrog). We hypothesized that temperature preference would increase for H. chrysoscelis tadpoles at larger body sizes to maximize the rate of metamorphosis and emergence as metamorphs, a trend similar to that seen in other North American anurans (Hutchison and Hill, 1978; Floyd, 1984; Dupré and Petranka, 1985; Dupré et al., 1986).
Materials and Methods
Study Site and Focal Species.—
The distribution of H. chrysoscelis extends from the coastal plain of the southeastern United States northward through the most of the midwestern and northeastern states on into Canada (Cline, 2005). H. chrysoscelis typically breeds between April and July, depending on the location (Cline, 2005). Individual clutches contain 30–45 eggs, and the larval period (egg to metamorphosis) is usually 45–65 d, which depends upon the thermal regime of the location (Cline, 2005; Mitchell and Pague, 2014). We collected H. chrysoscelis tadpoles at various development stages on 4–5 June 2018 (n = 48) and 10 July 2018 (n = 16). We collected all individuals from an artificial wetland and concrete retention pool at the Tennessee Aquarium Conservation Institute (TNACI), Chattanooga, Tennessee (35°04′29′′N, −85°20′13′′W). We did not collect water quality or temperature data of the collection site. TNACI is located in the temperate zone of the ridge and valley physiographic province in Tennessee. The site is adjacent to the impounded Tennessee River (Nickajack Reservoir) at an elevation of 200 m a.s.l. Once collected, all individuals acclimated to room temperature (22.8°C) in a laboratory at TNACI for a minimum of 1 h, and we recorded body length (in millimeters) and wet mass (in grams). Tadpoles were measured for length by using a ruler attached to the inside of a plastic container (118 ml [16 oz]). This container held enough water to cover the ruler and individual tadpoles, which minimized tadpole time spent out of the water. Next, the tadpoles were blotted dry before recording their mass. The entire measurement process used standardized approaches so that handling was minimal and consistent between all individuals, lasting ∼1–2 min per tadpole from start to finish. To our knowledge, no individuals were adversely affected by the measurement process.
Temperature Preference Experiment.—
We set up a temperature preference trial following Strickland et al. (2016). Sixteen aluminum rain gutters (150 cm × 10 cm) were arranged around a blind that was used to minimize behavioral disturbance of tadpoles caused by the presence of the observer. Each trough was situated with the lengthwise dimension perpendicular to the blind so that the observer could see the entire water space and tadpole. We filled each gutter with 2 cm of refrigerated, dechlorinated water. Water depths of maximum 2.5 cm have been used to reduce convection in similar experimental setups (Herreid and Kinney, 1967). We created a horizontal thermal gradient by placing one end of each gutter on a hot plate (low temperature setting) and the other end on ice. The hot plate was situated near the blind, whereas the ice bath was situated at the far end. In addition, we placed a handful of ice in the cold end of each gutter to speed up the cooling process and removed it before tadpole placement in the gutters. Because we only had 16 gutters, we conducted four separate trials and gave the individual trials a unique identifier (i.e., trial number). We conducted several of these trials on different dates. At the start of trials, the thermal gradient (mean ± SD) was from 13.2 ± 4.0°C to 25.9 ± 3.1°C at each extreme. At the end of trials, the thermal gradient was from 17.7 ± 1.8°C to 28.6 ± 3.3°C. Trials were conducted between 1000 and 1600 h on each given day.
We placed one tadpole in the center of each trough and allowed 15 min of acclimation before starting the trial. Temperature readings were taken every 15 min for 1 h by using an infrared thermometer (LaserGrip 774, Etekcity) for a total of four temperature readings per individual. Although distance from the target can affect infrared thermometer accuracy, we tested our thermometer and noted that temperature readings from 20 to 164 cm away differed on average by 0.1°C, which we deemed an adequate level of accuracy because the resolution of the thermometer was 0.1°C. Each temperature reading was collected at the body of the tadpole to provide an accurate reading of the temperature selected by the tadpole at that point in time. After the completion of each trial, we returned tadpoles to the artificial wetland surrounding the retention pool to ensure we did not recapture and test the same individual multiple times. All descriptive statistics represent means and SDs.
Analyses.—
We used a linear mixed-effects model (package lme4; Bates et al., 2015) to analyze the relationship between tadpole body length and temperature preference. Given that body length and tadpole mass were correlated (r2 = 0.834; F1,62 = 312.32, P < 0.001), we only analyzed body length as an explanatory variable, rather than both body length and mass. We averaged the four recorded temperatures for each individual to calculate their average temperature preference during the trial. Because all gutters were not oriented in the same direction with respect to the room, and not all trials were conducted on the same day, we wanted to control for the possibility that gutter orientation or trial date affected tadpole behavior. For these reasons, the gutter and trial date were tested in as well as the date a trial was conducted were included as random effects. We tested the residuals of the linear model with a Shapiro–Wilks test to ensure we met assumptions of normality. We used Satterthwaite's method, applied using the package lmerTest (Kuznetsova et al., 2017), to calculate df and p values for our model. We calculated conditional r2 for our model using the function rsquared from the R package piecewiseSEM (Nakagawa and Schielzeth, 2013; Lefcheck, 2016). All analyses were conducted in R software (R Core Team, 2018). Raw data can be found in Appendix 1.
Results
The body length of tadpoles was 12.2 ± 2.1 mm (range, 7–15 mm). Tadpoles had average temperature selections ranging from 12.9 to 30.0°C. Our model explained 13.05% of the variance, and we found a significant, positive relationship (F1,62= 9.454, P = 0.0031, β = 1.0606 ± 0.3449 [mean ± SE]) between body length and mean temperature preference (Fig. 1).



Citation: Journal of Herpetology 55, 1; 10.1670/20-020
Discussion
Temperature plays an integral role in anuran life cycles, wherein it can influence both metabolic (Wu et al., 2007) and developmental (Smith-Gill and Berven, 1979; Goldstein et al., 2017) rates, which contribute to the maturation of tadpoles into adults. Our study demonstrates that preferred temperature increases as body size increases in H. chrysoscelis tadpoles. This pattern is similar to that found in several other anuran tadpoles (Hutchison and Hill, 1978; Floyd, 1984; Dupré and Petranka, 1985; Dupré et al., 1986).
There are two main hypotheses to explain the positive relationship between tadpole body size and temperature preference: (1) escape from environmental pressures and (2) escape from ecological pressures. Seeking out warmer temperatures later in development could be an attempt to increase the rate of development, as tadpoles at later stages of development are in a more vulnerable state (Dupré and Petranka, 1985). From an environmental perspective, this vulnerability is often the result of environmental factors, such as pond drying (Wellborn et al., 1996; Wilbur, 1997), that will lead to mortality if tadpoles cannot escape the deteriorating environment in time. Selecting warmer temperatures generally increases the rate at which development proceeds, thus enabling tadpoles to emerge from their aquatic habitats to avoid desiccation (Newman, 1992).
In an ecological sense, increased temperatures can also enable tadpoles to escape ecological pressures, such as predation and competition. Selection of warmer temperatures may decrease the threat of predation by either increasing locomotion of the tadpoles (Niehaus et al., 2011) or by enabling them to metamorphose and leave the pond sooner (Wilbur, 1997). Dupré and Petranka (1985) discuss that the correlation between thermal preference and body size of tadpoles might also spatially segregate size classes and minimize competition and cannibalism in aquatic habitats, although cannibalism is not known in H. chrysoscelis. Differences in thermal preference may even contribute to interspecific niche partitioning, allowing multiple species to coexist in a given habitat by exploiting different temperature regimes.
Although our study demonstrates a relationship between preferred temperature and body size in H. chrysoscelis tadpoles, it is narrow in scale as it focuses on a single population in the southern third of its North American range. The tadpoles used in our study were from an artificial wetland and retention pond, which might not mimic a natural thermal regime. Temperature preferences may vary among anuran populations due to local thermal adaptations (Drakulić et al., 2017). Therefore, our results might represent a local thermal adaptation to artificial habitats that might differ if compared with thermal preferences of tadpoles from natural habitats. Because of its cosmopolitan distribution in North America, temperature preference of H. chrysoscelis could vary considerably among populations depending on local climates, and future studies should investigate temperature preference of H. chrysoscelis at a larger scale.

Positive relationship between average temperature selected (Celsius) by individual H. chrysoscelis tadpoles and their body length (conditional r2 = 0.1305, F1,62= 9.454, P = 0.0031, β = 1.0606 ± 0.3449 [mean ± SE]). Gray represents individual temperature measurements, and black represents average temperature selection calculations for each individual. The regression line is shown with a 95% confidence interval. The histogram on the right margin shows the frequency with which tadpoles selected different temperatures.
Contributor Notes