- Research article
- Open Access
Geographic body size variation of a tropical anuran: effects of water deficit and precipitation seasonality on Asian common toad from southern Asia
BMC Evolutionary Biology volume 19, Article number: 208 (2019)
Two previous studies on interspecific body size variation of anurans found that the key drivers of variation are the species’ lifestyles and the environments that they live in. To examine whether those findings apply at the intraspecific level, we conducted a study of the Asian common toad (Duttaphrynus melanostictus), a terrestrial anuran distributed in tropical regions. The body size of toads from 15 locations, covering the majority of their geographic range, and local environmental data were summarized from published literature. We used a model selection process based on an information-theoretic approach to examine the relationship between toad body size and those environmental parameters.
We found a positive correlation between the body size of the Asian common toad and the water deficit gradient, but no linkage between body size and temperature-related parameters. Furthermore, there was a positive correlation between the seasonality of precipitation and body size of females from different sampled populations.
As a terrestrial anuran, the Asian common toad should experience greater pressure from environmental fluctuations than aquatic species. It is mainly distributed in tropical regions where temperatures are generally warm and stable, but water availability fluctuates. Therefore, while thermal gradients are not strong enough to generate selection pressure on body size, the moisture gradient is strong enough to select for larger size in both males and females in dryer regions. Larger body size supports more efficient water conservation, a pattern in accordance with the prediction that lifestyles of different species and their local habitats determine the relationship between body size and environment. In addition, larger females occur in regions with greater seasonality in precipitation, which may happen because larger females can afford greater reproductive output in a limited reproductive season.
Body size is an important life history trait that influences many aspects of an individual’s biology , and identifying the forces (i.e. environmental gradients, interaction with sympatric species) that influence geographic body size variation across different species has implications for understanding how animals adapt to different abiotic and biotic environments [2,3,4,5,6]. A well-known clinal pattern of body size variation is Bergmann’s rule, which describes the tendency for body size to increase with lower temperatures or higher latitudes, because body size evolves under selective pressure from ambient temperature; being larger helps maintain body heat in cold environments (the heat conservation hypothesis [7, 8]). The majority of endothermic vertebrates show a Bergmann’s cline at the intraspecific  as well as interspecific [10, 11] level, and environmental temperature gradient is a decisive abiotic factor in body size diversification [7,8,9,10,11]. For insight into the issue, it was found that some anurans also show a Bergmann’s cline, with environmental temperature acting as the primary determinant of geographic size variation like in endotherms at both interspecific and intraspecific levels [12, 13]. Specially, many intraspecific investigations demonstrated that some life history traits diversify, interact and trade off under the temperature gradients and represent in the form of body size variation along the gradients [3, 4, 14,15,16,17,18,19,20]. Generally, a bigger size due to longer longevity in colder environments also facilitates heat conservation for many anuran species [3, 4, 16]. Nevertheless, there are also many species show a reversed cline [17,18,19], while some others show no geographic trend [20,21,22] even some life history traits show some general clines [17, 19, 21, 22]. Furthermore, aside from temperature as proposed by the heat conservation hypothesis, environmental factors such as water deficits and seasonality of temperature are emphasized in other hypotheses concerning body size variation in ectotherms, e.g. the water availability hypothesis ; the hibernation hypothesis  and the heat balance hypothesis . Moreover, because males and females may suffer different pressures from the same gradient, some species show sex-specific relationships between body size and environmental/geographic gradients [17, 20, 21, 25, 26]. The question of why there are such disparate trends among anurans remains unanswered. Studies of other clades suggest that these different relationships and new relevant factors emerge due to species-specific traits (i.e. habitat preference or thermal niche) [2, 27, 28] and the particular characteristics of each sex . As for anurans, the body size of anurans in Europe and North America are negatively correlated with and environmental energy gradients in accordance with most previous findings [3, 4, 7,8,9,10,11,12,13]. Nevertheless, anurans from the Brazilian Cerrado in South America showed a positive correlation between body size and water deficit, a measurement of environmental dryness . Those findings indicate that body size in anurans is strongly affected by environmental background: in regions where temperature fluctuates more, but precipitation fluctuates less (e.g., Holarctic), the major forces are thermal conditions, but in regions where temperature fluctuates less than precipitation (e.g., Neotropics), precipitation becomes the limiting factor . Based on these findings, Olalla-Tárraga et al. (2009) proposed the water-energy conservation hypothesis, which posits that larger bodies allow greater heat conservation in cold macroclimates and greater water conservation in dry tropical areas . A subsequent study on Chinese anurans showed that terrestrial anurans were generally larger in colder regions because they face temperature stress directly from the ambient environment, and becoming larger helps to conserve body heat. However, size of aquatic species did not change along thermal gradients. This could be due to the buffering effect of water, which reduces the pressures from macroenvironmental fluctuation. In other words, the relatively stable environmental gradients aquatic anurans experience are not strong enough to drive the evolution of their body size. Ecological traits (referring, for example, to the habitat preference of different species, i.e. terrestrial/aquatic/arboreal) of different clades also play an important role in the body size-environment relationship . Furthermore, a recent interspecific study on new world anurans further indicated that clades with different habitat preferences suffer different selective pressures . These studies of anurans on the American continents and in China focused on interspecific variation in body size [30,31,32]. However, the mechanisms that drive interspecific variation in body size should also apply to intraspecific variation [33, 34]. At present, studies considering both the ecological traits and environmental background within single anuran species are rare. A follow-up study found that body size in Darwin’s frogs (Rhinoderma darwinii) is positively correlated with the magnitude of seasonal temperature differences. Such a linkage may arise due to metabolic rate is reduced further and longer during colder, longer winters, leading to decreased energy depletion during hibernation, improved survival and increased longevity (and hence growth) (a mechanism termed the hibernation hypothesis) . Another study on male lesser treefrogs (Dendropsophus minutus) found that body size is positively linked with increased seasonality of precipitation, in accordance with the predictions of the water-energy conservation hypothesis . Nevertheless, neither study addressed the role of ecological traits and the environmental background in the patterns of variation observed within the focal species. The patterns discovered may owe to direct selective pressures from ambient environmental gradients (terrestrial and arboreal) in both species, but selective pressures probably differ because they live in different regions (temperate forests vs. tropical regions) [24, 35]. More work on other species in other regions is needed to clarify the issue. The Asian common toad (Duttaphrynus melanostictus) is a terrestrial anuran distributed widely throughout Southeast Asia and as far west as India and Pakistan (Fig. 1d) . As a terrestrial species, it spends the majority of its life on land and may suffer more extremes of ambient environmental conditions than aquatic species [31, 41,42,43]. Its body size may evolve differentially due to localized environmental pressures, with populations from different locations exhibiting a size cline along some environmental gradients. The Asian common toad faces macroenvironmental conditions similar to those of the Brazilian Cerrado , with warm and stable temperatures that do not fluctuate much across its geographic range (inferred by overlapping its distribution map with environmental data downloaded from WorldClim  in ARCGIS 10.0). However, precipitation conditions across its distribution range are characteristically variable, fluctuating across both space and time [36, 39]. As a result, populations of Asian common toads in different locations may experience more environmental pressure from precipitation than from temperature . The patchy environmental conditions that Asian common toads face allow us to test the validity of the water-energy conservation hypothesis within a single species. A previous study conducted on body size variation in the Asian common toad within a small geographic range (Fig. 1d) found that this species exhibited a reversed Bergmann’s cline . However, sampling a relatively small proportion of territory in a species with a broad geographic range may give an incomplete picture of intraspecific variation . In the present study, we summarized data covering the majority of the Asian common toad’s distribution to explore the relationship between body size and environmental predictors to gain insight into intraspecific responses to variations in environmental background. We hypothesized that toads would be larger in regions where the environmental water deficit is high, but that there would be no relationship between body size and thermal gradients 
Species and environmental data
The dataset was summarized from the literature, which included 15 populations covering the majority of the distribution of the Asian common toad. The mean body size of males ranged from 50.8 mm to 90.6 mm, and the mean body size of females ranged from 58.0 mm to 107.1 mm (Fig. 1d; Additional file 1: Table S1). The annual mean temperature experienced by each population ranged from 15.8 °C to 28.1 °C; the water deficit level ranged from 0.46 mm to 66.28 mm; and the precipitation seasonality ranged from 49 to 104 (Fig. 1a-c; Additional file 1: Table S2). The mean snout-vent length value of both males and females were spatially autocorrelated (Moran’s I = 0.44, P = 0.001; Moran’s I = 0.28, P = 0.026 separately), indicating that nearby populations have similar body sizes.
Correlation between body size and environmental gradients
When sex differences were considered, of the 32 models constructed, results showed that the mean body size of male Asian common toad was unaffected by environmental temperature gradients but was influenced by environmental dryness level (Table 1; Fig. 2a and b). The best model contained water deficit only, with mean male body size positively and significantly correlated with water deficit (β = 0.369 ± 0.088, P = 0.001; Fig. 2b). This factor alone explained 54.2% of the body size variation in males; no linkage with other predictors was detected (Table 1). Similarly, among females, mean female body size was unaffected by environmental temperature gradients (Table 1; Fig. 2a), and water deficit alone was the second-best-fit model, explaining 33.1% of the variation (β = 0.439 ± 0.156, P = 0.015; Fig. 2b). In contrast to the males, however, the body size of female Asian common toads was also influenced by the seasonality of precipitation, as female body size was positively and significantly correlated with environmental precipitation seasonality. For females, the best model contained precipitation seasonality only, which alone explained 33.8% of variation in body size (β = 0.576 ± 0.202, P = 0.014; Fig. 2c). The residuals of the best models were not spatially autocorrelated (Monte Carlo permutation test, all P > 0.1 for mean and both sexes), indicating that spatial autocorrelation did not bias our results. Sex differences were ignored in the analysis in some previous studies. When we combined the data for both sexes, of the 32 models constructed, the model selection process indicated that body size was unaffected by either environmental temperature gradients or by actual evapotranspiration, as in the sex-specific analysis, but the variation also could not be explained by precipitation seasonality (Table 1). The mean body size of Asian common toads was affected by water availability only, with a positive correlation between mean body size and water deficit (β = 0.362 ± 0.121, P = 0.010), and it alone explained 36.5% of the variation. The residuals of the best models were not spatially autocorrelated (Monte Carlo permutation test, P > 0.1). Due to the high collinearity between annual temperature, temperature seasonality and PET (r > 0.8, Additional file 1: Table S3), we conducted an analysis that excluded annual temperature. Of the 64 models constructed, we had the same results as with the analysis that included annual temperature (Additional file 1: Table S4). In addition, when the population from Bangalore was excluded because it had less accurate data (Additional file 1: Table S1), we obtained the same general results as in the analysis that included Bangalore (Additional file 1: Table S5).
Previous studies of the evolution of anuran body size gave insufficient attention to the ecological traits of different animal clades and focused exclusively on the role of environmental temperature gradients [9,10,11,12,13]. In contrast, we found that in one species of terrestrial anuran living in tropical regions, the key drivers of body size variation were water-related gradients. Different relationships between leading factors in the environment that determine body size in anuran clades may arise due to species-specific ecological traits and the particular environmental background they inhabit [30, 31]. Previous studies on Bergmann’s cline among anurans found a variety of correlations with environmental conditions [3, 4, 12,13,14,15,16,17,18,19,20,21,22,23,24,25,26, 30,31,32], and some species showed a coincident cline indicating that temperature seasonality  and absolute temperature  were the critical factors in the evolution of their body sizes. However, other studies found that precipitation gradients had an important mediating role in the relationship between body size and environment. Interspecific body size variation between anurans in the Brazilian Cerrado was influenced by water deficit level , while body size in the male lesser treefrog was influenced by precipitation seasonality . The different environmental backgrounds those species face may result in the diversified clines [30, 45]. On the other hand, many species showed no or a reversed cline along the thermal and precipitation gradients that produce Bergmann’s clines in other species [17,18,19,20,21,22]. We argue that such patterns may be explained by the fact that previous studies ignored the diverse ecological traits of these species. For instance, aquatic species experience relatively smaller selective pressures from fluctuations in the macroenvironment due to the buffering effect of the aquatic microenvironment . Arboreal and terrestrial species would be expected to respond more to fluctuations in the macroenvironment as they face these conditions directly [24, 42, 43]. Darwin’s frog is a terrestrial frog that lives in temperate forests, and the lesser treefrog is an arboreal species that lives in the tropics. In each species, body size was found to covary with environmental gradients, but in those studies little attention was paid to the important role of their ecologies and environmental backgrounds [24, 35]. Similar patterns and mechanisms also exist among other animal clades. Snakes in Australia with different ecological traits (nocturnal/diurnal activity; fossorial/surface lifestyles) showed different body size-environment relationships. In addition, due to the hot Australian macroenvironment, few areas are cold enough to limit daily activity. Low temperature in other regions [46, 47] does not act on their body size, while high summer temperatures might actually limit their daily activities . Furthermore, some endotherms show the same correlations. The body size of European herbivores showed a weaker relationship with temperature relative to that of carnivores and omnivores , and body size in tuco-tucos (Ctenomys), a subterranean rodent clade from South America, showed a reversed cline to Bergmann’s rule both between and within species . This is likely due to their subterranean lifestyle, which protects them from external temperature fluctuations and subjects them to other selective pressures, such as precipitation levels driving body size adaptation to different digging conditions . The roles of ecological traits and environmental background seem to be important in the body size-environment relationship in all animal clades. For terrestrial anurans, larger body size in colder climates can help to conserve body heat . Such a mechanism is advantageous in regions where the temperature fluctuates more (i.e., the Holarctic). However, in regions where temperatures fluctuate less, selective pressures from temperature gradients are not strong enough to drive body size diversification [30, 49]. In tropical regions, temperatures are usually high and stable, but precipitation gradients fluctuate. Anurans need to maintain their water balance while still conserving energy . In other words, the key factor influencing the body size-environment relationship among anurans in tropical regions is the dryness level, whereas the temperature gradient is the key factor in temperate regions . The Asian common toad is a terrestrial species that spends the majority of its life on land, except to reproduce. We predict that it will experience greater stress from macroenvironmental gradients than an aquatic species would [31, 43]. Further, this toad’s habitat is characterized by warm and stable ambient temperatures with little fluctuation, but precipitation varies greatly (Additional file 1: Table S2). The conservation of body water because of evaporative water loss becomes a priority, rather than how to conserve heat like terrestrial anurans in other temperate regions [30, 31, 50]. Therefore, body size in the Asian common toad may evolve with changing water deficit gradients . In accordance with our predictions, our results indicated that both male and female body size did not covary with temperature-related parameters, but there was a significant positive correlation between body size and water deficit for both sexes. In addition to the water deficit gradient, the body size of the female Asian common toad showed a positive correlation with precipitation seasonality. It is possible that no linkage between body size and environmental temperature gradient was observed because such mild gradients in temperature and PET do not generate enough selection pressure to drive body size [30, 49]. Our findings were inconsistent with the predictions of the hibernation hypothesis, which emphasizes the role of seasonality in temperature (Additional file 1: Table S4) . Unlike the ambient environment in the habitat of Darwin’s frog, temperature is not strongly seasonal across the Asian common toad’s range [36, 39], which means that the Asian common toad does not experience such drastic seasonal temperature fluctuations and selection pressure would be weak [30, 36, 49]. The Asian common toad is also far larger than Darwin’s frog (species mean body size 73.52 mm vs 23.7 mm), which allows it to reserve more energy to survive the subsequent year after hibernation; in some regions, the Asian common toad does not hibernate at all [51,52,53]. The hibernation hypothesis may be more applicable to small species living in temperate regions. The water deficit conditions across the distribution of the Asian common toad are similar to fluctuations in water availability in the Brazilian Cerrado [30, 36, 39]. This means that the water deficit gradient is more likely to become a selective pressure on body size [30, 49]. Consistent with our predictions and with the results of the interspecific comparison conducted in the Cerrado , the body size of both male and female Asian common toads showed a positive correlation with water deficit gradient; populations in dryer areas are larger-bodied than populations in wetter areas in accordance with the predictions of the water-energy conservation hypothesis . Becoming larger helps to conserve body water in dry environments. Further, in contrast to previous studies [13, 24, 30, 31, 35], we also found a positive correlation between female body size and precipitation seasonality, meaning that female Asian common toads grow larger in regions where there is high seasonality in precipitation. Seasonal changes in precipitation can affect the reproductive strategies of female toads [51,52,53,54,55]. In regions with strong seasonality, total yearly precipitation is concentrated into just a few months . Although these toads spend the majority of life on land, mating and tadpole development must occur in water, typically a temporal pool [56, 57]. Thus, reproductive activity is constrained by the availability of these pools [53,54,55,56,57]. In regions characterized by strong seasonality in precipitation, the window for the toad to mate and for the tadpoles to develop is relatively short. Under such conditions, larger bodies enable female Asian common toads to give birth to as many offspring as possible in a limited time . Where precipitation is more evenly distributed throughout the year, in contrast, every month is suitable for toads to reproduce. For instance, toads have been observed in amplexus every month in Singapore, where the precipitation seasonality index is low [36, 39, 58]. Under such environmental conditions, female Asian common toads are no longer pressured to reproduce in bursts. Similar relationships between body size and seasonality of precipitation were also observed in female Polypedates leucomystax, an arboreal species that grows bigger and produces larger eggs in central Thailand than in Singapore. Nevertheless, two co-occurring aquatic species, Microhyla heymonsi and Hylarana erythraea, collected from the same locations exhibit neither body size nor egg size variation as seen in P. leucomystax . The male Asian common toad is not affected by precipitation seasonality in the same way as the female. Differences in the strongest predictors of body size between the sexes may arise due to life history traits, since the body size of each sex is under slightly different selection pressures despite experiencing the same environment [17,18,19,20,21,22]. Female Asian common toads are larger than males (t = − 6.42, P < 0.001; 80.0 mm v.s. 67.0 mm). These larger females suffer relatively less pressure from reduced water availability than males, and as a result female body size is not as tightly coupled to the water deficit gradient . In addition to becoming larger in dry regions to conserve water, female Asian common toads need to lay as many eggs as possible in the short reproductive window such environments allow. As a result, the body size of female toads in regions with highly seasonal precipitation should also increase for that reason. As expected, the body size of male Asian common toads did not vary in response to seasonality of precipitation . The reproductive rhythm of male toads is the same as that of females in the same population [52, 57, 61]. In seasonally breeding populations, female toads must lay as many eggs as possible while they can, which is a costly investment, but sperm production for male toads is far less energy-intensive  and therefore would not be expected to affect male body size. These sex-specific reproductive costs are responsible for the different relationships between body size and the precipitation seasonality index observed in males vs. females. In contrast to our findings in the Asian common toad, male lesser treefrogs are bigger in regions with highly seasonal precipitation . Though both terrestrial anurans and arboreal anurans should be more sensitive to the macroenvironment than aquatic anurans [31, 42], there are other ways to conserve heat and energy besides growing larger. Tree frogs typically have thicker, waterproof skin to guarantee their water balance , while terrestrial species may adapt behavioral strategies to achieve the same goal [41, 62, 63]. The Asian common toad is larger than lesser treefrogs (species mean body size for males: 67.0 mm v.s. 23.7 mm) , which may give them more protection from dry conditions. Therefore, the differences in observed relationships between body size and environmental gradients in the lesser treefrog and the Asian common toad may result from their different lifestyles. More generally, both studies indicated a link between body size and environmental gradients, although the specific factor that explained the variation was dependent on environmental background. The unique responses in body size of male and female toads to environmental gradients suggests that we should take sex into account when analysing intraspecific and interspecific variation, using mean values for each sex as the measurement of species-specific body size [24, 64]. Collapsing the values of male and female body sizes into a single mean for the species may result in information loss; for instance, in the present study, we can detect the correlation between females and seasonality of precipitation gradients only by analyzing the mean SVT for females separately. Our study demonstrated that water deficit plays a key role in geographic body size variation of a tropical terrestrial anuran. In addition, it indicated that precipitation seasonality is another important environmental factor driving body size variation of anurans. Together, the results from studies of the Asian common toad, lesser treefrog and Darwin’s frog suggest that this variation is primarily dependent on the environmental background and ecological traits of different species. For terrestrial anurans in general, water-related gradients act as the main selective pressure in tropical regions, while thermal gradients act as the main pressure in temperate regions . Beyond the role of environmental background and lifestyle, we suggest that species-specific body size itself is an important factor determining geographic variation patterns among anurans. A species’ size is an ecological trait [1, 65, 66], and a relatively larger body size facilitates resistance to harsh environmental conditions, while smaller species are under greater pressure from ambient environmental conditions [65, 66]. For example, a comparative study indicated that body size of lesser treefrog becomes bigger while Blacksmith treefrog (Boana faber) enhance its ability of resistance to evaporative water loss with the dehydrating gradient because of their species-specific body size (21.8 mm v.s. 90.8 mm) while enlarge body size for lesser treefrog is more efficient to keep water balance . In the present study, water deficit explains at most 52% of body size variation in male Asian common toads (67.58 mm), while it explains only 31.7% of the variation for females (80.45 mm). However, the seasonality of temperature alone explains 87.6% of the size variation in the much smaller Darwin’s frog (23.7 mm) . As an ecological trait, body size itself will influence the intraspecific body size-environment relationship, with larger species being less sensitive to environmental gradients. It should be noted that all data in this study were summarized from the literature. Since the measurements were conducted by different people (Additional file 1: Table S1), there may be some errors in the dataset. Thus, further study in the field or laboratory will be required to replicate these results. Another limitation is that we explained at most only half of the variation with environmental factors, which means that more than half of the observed patterns are not associated with local climate. Size-independent traits such as behavioral  and psychological  strategies or different waterproof ability  or interaction of those strategies  are employed to survive under shifting environmental conditions. Besides, many previous studies on intraspecific size variation of anurans indicated that life history traits such as growth rate, longevity, age at sexual maturity, clutch size and egg size may all contribute to geographic body size variation and many authors concerned more on the effects of those traits on the anuran body size variation [3, 4, 14,15,16,17,18,19]. However, those intrinsic traits also covary with different selective pressures like body size [3, 4, 16, 17, 67]. The life history traits of different geographical populations represent interaction and trade-offs under particular environmental conditions and are not fully separable from body size, and body size is the representation of those life history traits evolve and interact with each other [3, 4, 16, 67]. Thus, it is improper to emphasize and conclude that those life history traits alone can produce the observed geographic patterns in body size. On the contrary, it is reasonable to explore how body size (the outcome of all other life history traits evolve and interact with each other along different environmental gradients) evolves under the ambient environments. Similar to the Bergmann’s rule [9, 26, 33], though many intraspecific studies have found some general trends in the diversification of life history traits , many species do not show the same trends [19, 20, 44, 68], it is more interesting to explore why some species do not show the same cline with the general trend. Under the framework of intraspecific Bergmann’s rule, it will gain better insight into how body size evolves if we consider how different life history traits evolve and interact with each other under the same selective pressures and how many of them contribute to the size variation alone . Furthermore, interaction with sympatric species in the community is also a driving force act on the size variation of animals [5, 6], however, due to the limitations of our method, we could not account for those life history traits and relationship of the toad with other species in the community in this study. Further studies conducted in the field or laboratory could examine the evolution of those life history traits under environmental gradients and their relationship with body size in greater detail [44, 67,68,69]. Besides, accounting for issues such as interaction with sympatric species [5, 6], time spent in hibernation [50, 51, 54, 70, 71] or behavioral  and physiological  adaptation is also required to gain insight into all of the potential mechanisms driving body size evolution not only in the Asian common toad, but in anurans more generally, across their various habitats.
Using the Asian common toad as a focal species, our study demonstrated that at the intraspecific level, body size variation of terrestrial anurans in the Asian tropics is primarily determined by water-related, rather than temperature-related, environmental parameters, a pattern consistent with the predictions of the water-energy conservation hypothesis. Additionally, we documented that even in the same environmental conditions, females and males of the same species may evolve different adaptive strategies. Combined with findings in Darwin’s frog and the lesser treefrog, our results suggest that selective pressures at the intraspecific level differ with changes in the macroenvironment.
Species data and environmental predictors
Snout-vent length (SVL) was used to measure body size in the Asian common toad. We summarized data on SVL and the central coordinates of each geographic population from the published literature (published research articles and local faunas: Additional file 1: Table S1). Because our goal was to explore whether the body size of different geographic populations covaries with particular environmental gradients, the mean (mean ± SD) SVL of mature individuals of each population was used as the measure of central tendency [24, 72]. Though previous studies indicated that life history traits such as age and growth rate also contribute to geographic variability in body size , those traits also evolve under different selective pressures with interaction and trade off exist among them . In other words, the life history traits of different populations trade off because of environmental gradients and can themselves be represented in the form of body size diversification [3, 4, 16, 19]. Since male and female amphibians may respond differently to climatic variation [17, 18, 20, 21], analyses performed on unsexed populations can be biased ; we employed the mean body size of both sexes together as the estimate of the population’s body size and separately to conduct our analysis. Furthermore, to better compare with the study of Darwin’s frogs  and interspecific studies [13, 30, 64], we conducted an extra analysis of unsexed mean SVL by population. Because body size evolves under direct selection resulting from different environmental pressures, analyses linking geographic parameters with body size probably cannot capture the causal forces behind the observed patterns [28, 32]. In the present study, we used several environmental gradients in lieu of examining toad body size variation by latitude or altitude. Following previous studies [13, 23, 24, 30, 35], seven environmental factors that may serve as evolutionary determinants of anuran body size were selected and extracted according to the geographic coordinates of the sampled populations that were collected (Additional file 1: Table S1 and Table S2). Then, the relationship between these environmental factors and the body size of Asian common toads was analysed to determine whether they work as selective pressures on the evolution of toad body size. The seven factors are: annual mean temperature (°C), a measure of heat in the environment; annual total precipitation (mm), a measure of water availability; temperature seasonality (standard deviation of temperature across months), an indicator of energy predictability; precipitation seasonality (coefficient of variation of precipitation across months), an indicator of water predictability; actual evapotranspiration (AET, mm), a measure of water and energy balance; potential evapotranspiration (PET, mm), a measure of heat and light inputs; and the water deficit (WD, mm), a measure of dryness level. Data on annual mean temperature, annual total precipitation, temperature seasonality and precipitation seasonality were extracted from WorldClim at a resolution of 0.167° × 0.167° grid cells . Data on actual evapotranspiration and potential evapotranspiration were extracted according to Willmott & Matsuura (2001) at a resolution of 0.5° × 0.5° grid cell . Water deficits were calculated by subtracting AET from PET [30, 37, 38].
Ordinary least squares regression based on Akaike’s information criterion (AIC) was used to generate multiple regression models, with body size as the response value and environmental variables as predictors. We compared the AICc scores of models with all possible predictor combinations. The model with the lowest score was selected as the final, best-fitting model, thus omitting uninformative predictors . Due to the high collinearity between annual temperature and temperature seasonality/PET (r > 0.8, Additional file 1: Table S3), we excluded temperature seasonality and PET in the final analysis. As a check, we also conducted an analysis including temperature seasonality and PET but excluding annual temperature, which produced similar results (Additional file 1: Table S4). Two-tailed significance levels were set at P = 0.05. Multiple regression was run in R (version 3.1.1-R Core Team 2015)  using the package ‘MuMIn’ . ARCGIS 10.0 was used to calculate Moran’s I as a measure of spatial autocorrelation of SVL , and the Monte Carlo permutation test (199 permutations) was used in SAM 3.0 to assess spatial autocorrelation of residuals .
Availability of data and materials
The datasets used and analysed during the current study are available in additional file 1 Table S1 and Table S2.
Akaike’s information criterion
- P. Seas:
Annual total precipitation
- T. Seas:
Annual mean temperature
Roff DA. The evolution of life histories: theory and analysis. New York: Chapman and Hall; 1992.
Rodriguez MA, Lopez-Sanudo IL, Hawkins BA. The geographic distribution of mammal body size in Europe. Glob. Ecol. Biogeogr. 2006;15:173–81.
Liao W, Lu X. Adult body size = f (initial size + growth rate × age): explaining the proximate cause of Bergman’s cline in a toad along altitudinal gradients. Evol. Ecol. 2012;26:579–90.
Liao WB, Luo Y, Lou SL, Lu D, Jehle R. Geographic variation in life-history traits: growth season affects age structure, egg size and clutch size in Andrew's toad (Bufo andrewsi). Front. Zool. 2016;13:6.
Adams DC, West ME, Collyer ML. Location-specific sympatric morphological divergence as a possible response to species interactions in West Virginia Plethodon salamander communities. J. Anim. Ecol. 2007;76:289–95.
Kozak KH, Mendyk RW, Wiens JJ. Can parallel diversification occur in sympatry? Repeated patterns of body-size evolution in coexisting clades of north American salamanders. Evolution. 2010;63:1769–84.
Bergmann C. Ueber die verhältnisse der wärmeökonomie der thiere zuihrer grösse. Gottinger Studien. 1874;3:595–8.
Mayr E. Geographical character gradients and climatic adaptation. Evolution. 1956;10:105–8.
Meiri S, Dayan T. On the validity of Bergmann’s rule. J. Biogeogr. 2003;30:331–51.
Olson V, Davies RG, Orme CDL, Thomas GH, Meiri S, Blackburn TM, Gaston KJ, Owens IPF, Bennett PM. Global biogeography and ecology of body size in birds. Ecol. Lett. 2009;12:249–59.
Clauss M, Dittmann MT, Müller DWH, Meloro C, Codron D. Bergmann’s rule in mammals: a cross-species interspecific pattern. Oikos. 2013;122:1465–72.
Ashton KG. Do amphibians follow Bergmann’s rule? Can. J. Zool. 2002;80:708–16.
Olalla-Tárraga MÁ, Rodríguez MÁ. Energy and interspecific body size patterns of amphibian faunas in Europe and North America: anurans follow Bergmann’s rule, urodeles its converse. Glob. Ecol. Biogeogr. 2007;16:606–17.
Liao WB, Lu X, Shen YW, Hu JC. Age structure and body size of two populations of the rice frog Rana limnocharis from different altitudes. Ital. J. Zool. 2011;78:215–21.
Liu YH, Liao WB, Zhou CQ, Mi ZP. Altitudinal variation in body size in the Rice frog (Rana limnocharis) in southwestern China. Acta Herpetol. 2012;7:57–68.
Yu TL, Wang DL, Busan M, Deng YH. Altitudinal variation in body size in Bufo minshanicus supports Bergmann's rule. Evol. Ecol. 2019;33:449–60.
Hsu FH, Hsieh YS, Wu SH, Kam YC. Altitudinal variation in body size and age structure of the Sauter's frog Rana sauteri in Taiwan. Zool. Stud. 2014;53:62.
Liu Q, Feng H, Jin L, Mi ZP, Zhou ZM, Liao WB. Latitudinal variation in body size in Fejervarya limnocharis supports the inverse of Bergmann's rule. Anim. Biol. 2018;68:113–28.
Zhang L, Ma X, Jiang J, Lu X. Stronger condition dependence in female size explains altitudinal variation in sexual size dimorphism of a Tibetan frog. Biol. J. Linn. Soc. 2012;107:558–65.
Li S, Wu X, Li DY, Lou SL, Mi ZP, Liao WB. Body size variation of odorous frogs (Odorrana grahami) across altitudinal gradients. Herpetol. J. 2013;23:187–92.
Lu X, Li B, Liang JJ. Comparative demography of a temperate anuran, Rana chensinensis, along a relatively fine altitudinal gradient. Can. J. Zool. 2006;84:1789–95.
Feng X, Chen W, Hu J, Jiang J. Variation and sexual dimorphism of body size in the plateau brown frog along an altitudinal gradient. Asian Herpetol Res. 2015;6:291–7.
Nevo E. Adaptive variation in size of cricket frogs. Ecology. 1973;54:1271–81.
Valenzuela-Sánchez A, Cunningham AA, Soto-Azat C. Geographic body size variation in ectotherms: effects of seasonality on an anuran from the southern temperate forest. Front. Zool. 2015;12:37.
Liao WB, Lu X. Age structure and body size of the Chuanxi tree frog Hyla annectans chuanxiensis from two different elevations in Sichuan (China). Zool. Anz. 2010;248:255–63.
Liao WB, Zhou CQ, Yang ZS, Hu JC, Lu X. Age, size and growth in two populations of the dark-spotted frog Rana nigromaculata at different altitudes in southwestern China. Herpetol. J. 2010;20:77–82.
Medina AI, Martí DA, Bidau CJ. Subterranean rodents of the genus Ctenomys (Caviomorpha, Ctenomyidae) follow the converse to Bergmann’s rule. J. Biogeogr. 2007;34:1439–54.
Rypel AL. 2014. The cold-water connection: Bergmann’s rule in north American freshwater fishes. Am. Nat. 2014;183:147–56.
Shine R. Ecological cause for the evolution of sexual dimorphism: a review of the evidence. Q. Rev. Biol. 1989;64:419–61.
Olalla-Tárraga MÁ, Diniz-Filho JAF, Bastos P, Rodríguez MÁ. Geographic body size gradients in tropical regions: water deficit and anuran body size in the Brazilian Cerrado. Ecography. 2009;32:581–90.
Guo C, Lu X. Lifestyle-based approaches provide insights into body size variation across environmental gradients in anurans. Evol. Biol. 2016;43:208–14.
Amado TF, Bidau CJ, Olalla-Tárraga MÁ. Geographic variation of body size in new world anurans: energy and water in a balance. Ecography. 2019;42:456–66.
Blackburn TM, Gaston KJ. Scale in macroecology. Glob. Ecol. Biogeogr. 2002;11:185–9.
Adams DC, Church JO. The evolution of large-scale body size clines in Plethodon salamanders: evidence of heat-balance or species-specific artifact? Ecography. 2011;34:1067–75.
Oyamaguchi HM, Oliveira E, Smith TB. Environmental drivers of body size variation in the lesser treefrog (Dendropsophus minutus) across the Amazon-Cerrado gradient. Biol. J. Linn. Soc. 2017;120:363–70.
Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 2005;25:1965–78.
Willmott CJ, Matsuura K. Terrestrial water budget data archive: monthly time series (1950–1999); 2001. Available at: http://climate.geog.udel.edu/~climate/html_pages/README.wb_ts2.html.
Francis AP, Curri DJA. Globally-consistent richness-climate relationship for angiosperms. Am. Nat. 2003;161:523–36.
IUCN. IUCN red list of threatened species. 2018. http://maps.iucnredlist.org/map.html?id=54707. Accessed 6 Jun 2018.
Fan XL. Thermal adaption of Duttaphrynus melanostictus: pattern and mechanism. Ph. D. Thesis. Nanjing: Nanjing normal university; 2014. (in Chinese)
Navas CA. Implications of microhabitat selection and patterns of activity on the thermal ecology of high elevation neotropical anurans. Oecologia. 1996;108:617–26.
Young JE, Christian KA, Donnellan S, Tracy CR, Parry B. Comparative analysis of cutaneous evaporative water loss in frogs demonstrates correlation with ecological habits. Physiol. Biochem. Zool. 2005;78:847–56.
Naya DE, Naya H, Cook J. Climate change and body size trends in aquatic and terrestrial endotherms: does habitat matter? PLoS One. 2017;12:e0183051.
Laugen AT, Laurila A, Jönsson I, Söderman F, Merilä J. Do common frogs (Rana temporaria) follow Bergmann’s rule? Evol. Ecol. Res. 2005;7:717–31.
Hawkins BA, et al. Energy, water, and broad-scale geographic patterns of species richness. Ecology. 2003;84:3105–17.
Olalla-Tárraga MÁ, Rodríguez MÁ, Hawkins BA. Broad-scale patterns of body size in squamate reptiles of Europe and North America. J. Biogeogr. 2006;33:781–93.
Aragón P, Fitze PS. Geographical and temporal body size variation in a reptile: roles of ecology, phylogeny and ecology structured in phylogeny. PLoS One. 2014;9(8):e104026.
Feldman A, Meiri S. Australian snakes do not follow Bergmann’s rule. Evol. Biol. 2014;41:327–35.
Foster F, Collard M. A reassessment of Bergmann’s rule in modern humans. PLoS One. 2013;8:e72269.
Gouveia SF, Correia I. Geographical clines of body size in terrestrial amphibians: water conservation hypothesis revisited. J. Biogeogr. 2016;43:2075–84.
Shieh JN. The breeding ecology of Bufo melanostitus. M.S. thesis. Taiwan: Tunghai University; 1993. (in Chinese)
Jorgensen CB, Shakuntala K, Vijayakumar S. Body size, reproduction and growth in a tropical toad, Bufo melanostictus, with a comparison of ovarian cycles in tropical and temperate zone anurans. Oikos. 1986;46:379–82.
Boyce MS. Seasonality and patterns of natural selection for life histories. Am. Nat. 1979;114:569–83.
Aichinger M. Annual activity patterns of anurans in a seasonal neotropical environment. Oecologia. 1987;71:583–92.
Bertoluci J, Rodrigues MT. Seasonal patterns of breeding activity of Atlantic rainforest anurans at Boracéia. Southeastern Brazil Amphibia-Reptilia. 2002;23:161–7.
Fei L, Hu S, Ye C, Huang Y. Fauna Sinica, Amphibia, Vol. 3. Anura Ranidae. Beijing: Science Press; 2009. (in Chinese)
Invasive Species Council. Biosecurity failures in Australia, The Asian black-spined toad, https://invasives.org.au/wp-content/uploads/2014/09/Biosecurity-failures-asian-black-spined-toads.pdf; retrieved 2017-05-28; 2014.
Berry PY. The breeding patterns of seven species of Singapore anuran. J. Anim. Ecol. 1964;33:227–43.
Sheridan JA. Reproductive variation corresponding to breeding season length in three tropical frog species. J. Trop. Ecol. 2009;25:583–92.
Simmons LW. Sperm competition and its evolutionary consequences in the insects. Princeton University Press: Princeton; 2001.
Gramapurohit NP, Radder RS. Mating pattern, spawning behavior, and sexual size dimorphism in the tropical toad Bufo melanostictus (Schn.). J. Herpetol. 2012;46:412–6.
Tracy CR, Christian KA, O’Connor MP, Tracy CR. Behavioral thermoregulation by Bufo americanus: the importance of the hydric environment. Herpetology. 1993;49:375–82.
Navas CA, Antoniazzi MM, Jared C. A preliminary assessment of anuran physiological and morphological adaptation to the Caatinga, a Brazilian semi-arid environment. Int. Congr. Ser. 2004;1275:298–305.
Slavenko A, Meiri S. Mean body sizes of amphibian species are poorly predicted by climate. J. Biogeogr. 2015;42:1246–54.
Tracy CR, Christian KA, Tracy CR. Not just small, wet, and cold: effects of body size and skin resistance on thermoregulation and arboreality of frogs. Ecology. 2010;91:1477–84.
Gouveia SF, Bovo RP, Rubalcaba JG, Da Silva FR, Maciel NM, Andrade DV, Martinez PA. Biophysical modeling of water economy can explain geographic gradient of body size in anurans. Am. Nat. 2019;193:51–8.
Morrison C, Hero JM. Geographic variation in life-history characteristics of amphibians: a review. J. Anim. Ecol. 2003;72:270–9.
Wang X, Huang Y, Zhong M, Yang S, Yang X, Jiang J, Hu J. Environmental stress shapes life-history variation in the swelled-vented frog (Feirana quadranus). Evol. Ecol. 2019;33:435–48.
Cvetković D, Tomašević N, Ficetola GF, Crnobrnja-Isailović J, Miaud C. Bergmann’s rule in amphibians: combining demographic and ecological parameters to explain body size variation among populations in the common toad Bufo bufo. J. Zool. Syst. Evol. Res. 2009;47:171–80.
Sinsch U, Dehling JM. Tropical anurans mature early and die young: evidence from eight Afromontane Hyperolius species and a meta-analysis. PLoS One. 2017;12:e0171666.
Church G. The invasion of Bali by Bufo melanostictus. Herpetologica. 1960;16:15–21.
Romano A, Ficetola GF. Ecogeographic variation of body size in the spectacled salamanders (Salamandrina): influence of genetic structure and local factors. J. Biogeogr. 2010;37:2358–70.
Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. New York: Springer Science & Business Media; 2002.
R Core Team. R: a language and environment for statistical computing. Vienna, a, Austria. http://www.R-project.org/; 2015.
Barton, K. MuMIn: Multi-model inference. R package version 1.9.0. 2013. Available at: http://CRAN.R-project.org/package=MuMIn. Last accessed 25 Jan 2014.
Lichstein JW, Simons TR, Shriner SA, Franzreb KE. Spatial autocorrelation and autoregressive models in ecology. Ecol. Monogr. 2002;72:445–63.
Rangel TF, Diniz-Filho JAF, Bini LM. SAM: a comprehensive application for spatial analysis in macroecology. Ecography. 2010;33:46–50.
We thank four anonymous reviewers for comments on earlier draft, their comments of the manuscript deepen my understanding of the geographic body size variation of anurans. We thank editors who patiently handled the manuscript during the submission process. We are grateful to Guo-Yue Zhang for advice on statistical analyses and American Journal Experts for language editing of the earlier manuscript.
The editing service and publication fees of the paper work were supported by research start-up funds of Central South University of Forestry and Technology (2016YJ029) to Cheng Guo. The funding body had no role in the design, collection and conclusion of this study.
Ethics approval and consent to participate
No experiments on live Asian common toads were included in this study. All data in the analyses were summarized from open-access information (published research articles and local faunas). All research methods adhered to the Chinese legal requirements and were evaluated and approved by Central South University of Forestry and Technology’s Institutional Animal Care and Use Committee (protocol no. 2012–018).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Mean body size data of Asian common toad from sampled populations across its distribution range (ranked from north to south). Table S2. Data of environmental predictors of sampled populations. Table S3. Correlation coefficients between each environmental variable (statistically significant [P < 0.05] are shown in bold). Table S4. Multiple regression models for mean body size of each sex of Asian common toad and environmental predictors (excluding annual mean temperature due to it is highly correlated with temperature seasonality and potential evapotranspiration). Models are ranked by AICc from the best- to worst-fitting models. Table S5. Multiple regression models for mean body size of each sex of Asian common toad and environmental predictors (excluding temperature seasonality and potential evapotranspiration due to they are highly correlated with annual mean temperature), excluding the population from Bangalore due to its less accurate data. Models are ranked by AICc from the best- to worst-fitting models
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Guo, C., Gao, S., Krzton, A. et al. Geographic body size variation of a tropical anuran: effects of water deficit and precipitation seasonality on Asian common toad from southern Asia. BMC Evol Biol 19, 208 (2019). https://doi.org/10.1186/s12862-019-1531-z
- Bergmann’s rule
- Duttaphrynus melanostictus
- Precipitation seasonality
- Water-energy conservation hypothesis
- Water deficit