Geographic variation in life-history traits: growth season affects age structure, egg size and clutch size in Andrew’s toad (Bufo andrewsi)
© Liao et al. 2016
Received: 17 October 2015
Accepted: 28 January 2016
Published: 9 February 2016
Environmental variation associated with season length is likely to promote differentiation in life-history traits, but has been little studied in natural populations of ectotherms. We investigated patterns of variation in egg size, clutch size, age at sexual maturity, maximum age, mean age, growth rate and adult body size in relation to growth season length among 17 populations of Andrew’s toad (Bufo andrewsi) at different latitudes and altitudes in the Hengduan Mountains, western China.
We found that egg size, age at sexual maturity, and mean age increased with decreasing length of the growth season, whereas clutch size showed a converse cline. Body size did not increase with decreasing length of the growth season, but was tightly linked to lifetime activity (i.e. the estimated number of active days during lifetime). Males and females differed in their patterns of geographic variation in growth rates, which may be the result of forces shaping the trade-off between growth and reproduction in different environments.
Our findings suggest that growth season plays an important role in shaping variation in life-history traits in B. andrewsi across geographical gradients.
A central aim in life-history research is to explain patterns of growth, development, reproductive investment and survival at the level of species, populations and individuals [1, 2]. Due to physical and physiological constraints and common dependence on limited resources, investment in one life-history trait is often traded off against investment in other traits . Optimal resource allocation to specific life-history traits for animals across broad environmental ranges depends on differential adaptation to local conditions such as temperature, food supply, predation, and competition, and on the ability for plastic responses [4–7]. Most studies which have examined environmental effects on life-history traits such as age, body size, growth, egg size and clutch size in ectotherms consider that temperature exerts a strong effect [5, 8–18]. For example, many ectotherms delay their reproduction with increasing latitude and/or altitude, while simultaneously decreasing offspring number and increasing offspring size likely as a strategy to enhance offspring survival [10, 16, 19, 20]. Body size is another important life-history trait which has a strong impact on individual fitness, and Bergmann’s rule states that body size and temperature are inversely correlated . Geographic variation in body size predominantly follows Bergmann’s rule in endotherms [22–25], whereas both Bergmann’s rule as well as the reverse of it are frequently reported for ectotherms [14, 26, 27].
In amphibians, latitudinal and/or altitudinal patterns of body size variation may be complex . The influence of temperature on amphibian growth during the larval and juvenile phases affects the timing of maturity and, subsequently, the size of adults . Time constraints associated with the length of the activity season may play an important role in shaping allocation patterns between investment into growth and longevity [26, 29]. For instance, shorter growth seasons should increase the time taken to reach adulthood, leading to a higher age at maturation [8, 30]. This is for example also observed in reptiles, for which longevity often increases with decreasing growth season [7, 31, 32].
Bufo andrewsi is a medium-sized anuran with female-biased sexual size dimorphism. The taxonomic validity of B. andrewsi as a species distinct from B. gargarizans has been debated in recent years [33, 34], but we regard B. andrewsi as a distinct species for the present study. Bufo andrewsi is widely distributed in subtropical forests of the Hengduan Mountains in China, at elevations ranging from 750 to 3500 m . The species has a relatively long spawning period, with egg laying taking place from early February to mid-June. During the breeding season, males clasp females when they enter the breeding pond, and females produce a single clutch. Egg volume and total reproductive output increase with altitude, whereas clutch size decreases with it . Moreover, in accordance with Bergmann’s rule B. andrewsi attain larger adult sizes at metamorphosis, higher average age and larger body sizes despite slower growth rates at higher altitudes . However, differences in average body size among anuran populations not only depend on average age but also on the effective number of days available for growth . Hence, our understanding of observed patterns of interpopulational variation in life-history traits across geographical gradients in amphibians would benefit from further insights into the role of growth season length.
Here, we document the effects of the length of the growth season on geographical variation in life-history traits in B. andrewsi. To this end, we analyzed a large dataset derived from 17 distinct populations from the Hengduan Mountains, China. In particular, we investigated how egg size, clutch size, age at maturity, mean age, maximum age and growth rates in males and females vary with variation in the length of the growth season for both sexes. We asked whether differences in effective number of days available for growth can explain the observed variation in body size along geographical gradients.
Sampling of populations
Following the protocol of Hjernquist et al. , we compiled data on the length of the growth season based on average ambient temperatures between 2007 and 2013 obtained from the Chinese Meteorological Administration (http://www.cma.gov.cn). Field observations revealed that the toads begin to be active when the average daily ambient temperature is approximately 6 °C, a threshold which does not vary among sites (unpublished data); we assume that the length of the growth season closely follows the length of the activity period. In our study area, the length of the growth season (i.e. number of days with a mean daily temperature above 6 °C) differs between localities by a maximum of approximately 150 days (Lingguan and Muyangchang, see Additional file 1: Appendix S1). Low latitude populations in the southern Yunnan province become active in early February, continuing activity until early December. In high latitude populations in the northern Sichuan province, on the other hand, the activity period starts in early May and finishes by mid-October. The length of the growth season can also be used to calculate the lifetime activity period of an individual, by multiplying age (in years) with the length of the growth season in days (following Hjernquist et al. ).
We collected a total of 208 amplectant toads from ten populations (see details in Additional file 1: Appendix S1) between 2008 and 2013, and transported them to laboratories close to spawning sites. We kept pairs separately in tanks (40 × 50 × 60 cm) filled with pond water and allowed them to oviposit. Once oviposition was completed, we counted the total number of eggs in each clutch, and measured SVL of males and females to the nearest 0.1 mm using callipers. We placed 100 randomly selected eggs from each clutch on a glass plate to take a digital image, and measured egg diameter using these photographs (to the nearest 0.01 mm, not including the jelly capsules) at 400× magnification using a Motic Images 3.1 digital camera mounted on a Moticam 2006 light microscope. After the experiments, we returned all individuals and egg strings to their place of collection.
Skeletrochronology is considered the most reliable method of age estimation in wild amphibians, and is based on counting lines of arrested growth (LAGs) in long bones [38, 39]. Paraffin sectioning and Harris’s haematoxylin staining was used to produce histological sections (for details see ). The collected digits were washed in water for 2 h and then decalcified in 5 % nitric acid for 48 h, before being washed in running tap water overnight and stained with Ehrlich’s haematoxylin for 75 min. The stained bones were dehydrated and embedded in small paraffin blocks. We selected the cross-section (approximately 13 μm thick) of the phalanx which had the smallest medullar cavity and the thickest cortical bone, and mounted it on glass slides. We recorded the number of lines of arrested growth from mid-diaphyseal sections under a light microscope, having previously confirmed that the number of LAGs reflects actual age based on four marked and recaptured individuals . We assessed endosteal resorption of LAGs by comparing the diameter of the smallest cross-section of 1-year old toads to the diameter of the resorption line in adults . Double lines and false lines could be readily distinguished from the true LAGs and therefore did not create problems for age estimation. In total, we determined the age of 1563 individuals (1166 males and 397 females).
The minimum age of adult toads at a spawning location was used as an estimate of the age at first reproduction in a population, and the maximum age as an estimate of longevity. This is reasonable because the toads captured at the spawning sites represent the age distribution of the reproductive population. We did not find a significant effect of sample size on minimum and maximum age in each population (both P > 0.12).
Von Bertalanffy’s  model, a standard method for describing the growth of animals with asymptotic growth after maturity [42–44], was used to estimate growth parameters. The equation has the form of SVLt = SVLmax (1 – e-kt+b), where SVLt is body size at age t for per individual, SVLmax is the estimated asymptotic maximum size in a population, k is a growth coefficient and b is a constant. Growth rate can then be presented as R = dSVL/dt = k (SVLmax –SVLt).
All analyses were conducted using Type III sums of squares tests in SPSS 21.0 (Statistical Product and Service Solutions Company, Chicago, USA). Differences in minimum age (or mean age or maximum age or lifetime activity) among the 17 populations were tested for using Generalized Linear Mixed Models (GLMMs), with minimum age (or age or maximum age or lifetime activity) as the dependent variable, growth season as fixed effect, and sex, latitude (meters from the equator) and altitude as covariates; population was treated as a random factor. To investigate variation in body size, we treated body size as the dependent variable, growth season as fixed effect, and age, sex, latitude and altitude as covariates, and population as a random factor. To investigate variation in egg size, we treated egg size as the dependent variable, growth season as fixed effect, and female size, latitude and altitude as covariates. Finally, to test variation in egg size and a trade-off between clutch size and egg size, we treated clutch size as the dependent variable, growth season as fixed effect, and female size, latitude, altitude, egg size as covariates.
Egg size and clutch size
Response of egg size (1), clutch size (2), minimum age (3), maximum age (4), average age (5), lifetime activity (6), and adult body length (7) to length of the growth season for 17 Bufo andrewsi populations along a geographical gradient as revealed by by General Linear Mixed Models (GLMM)
Source of variation
1. Egg size
2. Clutch size
3. Minimum age
4. Maximum age
6. Lifetime activity
7. Body size
Age and lifetime activity
Body size and growth rate
Although mean body size was not predicted by the length of the growth season, altitude or latitude, the effects of sex and age on body size were significant across populations (Table 1; Fig. 3e; sex: F 1,1806.435 = 3035.206, P < 0.001; age: F 1,1810.021 = 189.421, P < 0.001). The relationship between growth rate and length of the growth season tended to be sex-specific. Male growth rates increased with decreasing length of the growth season, whereas female growth rate tended to decrease with it (Fig. 3f).
Our analysis of life-history traits across 17 B. andrewsi populations provides strong support for the notion that life-history strategies respond to environmental conditions. Specifically, we find that egg size, minimum age at first reproduction, longevity and mean age increase with decreasing length of the growth season, whereas clutch size displays a converse cline. In other words, populations with longer growth seasons mature earlier, have a smaller mean age and longevity, and produce smaller eggs and larger clutch sizes. However, body size and lifetime activity in B. andrewsi does not increase with decreasing length of the growth season.
Like in most other ectotherms [7, 45], life-history traits and strategies of anurans are influenced by varying environmental conditions [5, 8, 13–16]. The observed increase in egg size with decreasing length of the growth season may be linked to adaptive beneficial effects on larval and metamorphic performance  and, hence, can also influence lifetime fitness [46, 47]. Larger eggs contain a higher amount of yolk , which is assumed to be advantageous under the high energetic requirements posed by cold environments. Deviations from the general pattern of egg size by the YL population might relate to differences in temperature regimes between breeding sites, for which no local temperature measurements are available. Maternal condition can influence investment in clutch size , and lower resource availability possibly leads to resource limitation on clutch production at shorter activity periods. We observed a trade-off between egg size and clutch size, suggesting investment in larger eggs comes at a cost to egg numbers in accordance with the basic tenets of life-history theory .
The larger age at sexual maturity and mean age at shorter growth season in B. andrewsi is consistent with previous studies [7, 16, 51]. In populations with longer growth seasons such as Lingguan, Muping and Qibie, toads matured at the age of 1-2 years, and the average age of reproducing toads was 5.1 years. At populations with shorter growth season such as Caopuo, Baozhigou and Muyangchang, individuals matured at the age of 3-4 years, and the mean age of reproducing individuals was 7.1 years. Earlier maturation and lower age in populations with longer growth seasons may reflect higher predation risks [52, 53]. In addition, variation in environmental temperature is likely correlated with fluctuation of food availability. Invertebrates, the major food resources of anurans, are less abundant in low-temperature environments . Low food availability at shorter growth seasons might increase juvenile mortality and juveniles might need longer to reach adulthood, resulting in a higher age at sexual maturity and mean age [45, 49]. In B. andrewsi, predation risk and food availability of juveniles and adults are expected to be lower in populations with shorter growth seasons . Thus, temperature-dependent metabolic rates may contribute to the observed variation in sexual maturation and mean age [10, 55]. In addition to minimum and mean age, maximum age in ectotherms often increases with decreasing length of the growth season [7, 16, 55].
As has been observed in previous studies [16, 29], mean body size did not increase with decreasing length of the growth season, despite a strong negative correlation between age at maturity and growth season length. For B. andrewsi, lifetime activity therefore cannot explain the observed variation in body size among populations, probably due to a non-significant correlation between the length of the growth season and body size variation. A negative correlation between growth rate and maximum age across geographic gradients has previously been suggested to explain observed body size clines . For instance, later reproduction and higher maximum age at long growth seasons play a more important role in increasing body size in anurans than slower growth does in reducing size, resulting in a negative correlation between body size and growth season [15, 43]. A lifetime growth period, on the other hand, will fail to compensate for the effects of slow growth on body size, resulting in a positive correlation between body size and growth season [14, 27]. In this study, a non-significant variation in body size may be attributed to different contributions of growth rates and age structures in determining body size at varying growth season lengths, in agreement with other studies on life-history traits and seasonality in ectotherms [29, 39].
Growth was more rapid in females compared to males for shorter growth seasons, while males had higher growth rates than females in longer growth seasons (interaction of both sex and lifetime activity: F 1,1513.102 = 2.204, P < 0.01; Fig. 3f). A previous study suggested that male B. andrewsi mature earlier at lower altitudes and latitudes, also showing a more male-biased operational sex ratio . Therefore, stronger male-male competition at longer growth seasons might also select for accelerated growth in males, whereas more balanced sex ratios and weakened male–male competition at shorter growth seasons may result in slower growth [18, 56]. Life-history theory describes a trade-off between reproductive investment and growth . For female B. andrewsi, lower reproductive investment and faster growth are observed in longer growth seasons under warmer conditions [15, 37]. In the present study, we also found a trade-off between reproductive investment and growth through slower growth and higher reproductive investment under environments with shorter activity periods.
Geographical variation in life-history traits in anurans are driven by environmental and genetic differences [5, 26, 57, 58], and plastic responses to environmental conditions are common [26, 59, 60]. Approaches such as common garden experiments are needed to verify the genetic basis for population differences in life-history traits in B. andrewsi in the future.
Taken together, our study shows that individuals from populations with longer seasonal activity periods produce smaller eggs, mature earlier and have larger mean age, whereas individuals with shorter activity periods produce smaller clutches. However, the body sizes do not increase with decreasing length of the growth season across geographical populations. In addition, males and females exhibit differential variation in growth rates across geographical gradients, which can be explained by a trade-off in resource allocation between growth and reproduction in different environments.
Financial support was provided by the National Natural Sciences Foundation of China (31471996), Sichuan Province Outstanding Youth Academic Technology Leaders Program (2013JQ0016) and Sichuan Province Department of Education Innovation Team Project (15TD0019). All fieldwork performed in Fengtongzhai Nature Reserve, Wanglang Nature Reserve, Baima Snow Mountains Nature Reserve and Yele Nature Reserve was permitted by the State Forestry Administration of China. The field studies did not involve endangered or protected species. All animal work was approved by China West Normal University Institutional Animal Care Use Committee (IACUC) and Ethics Committee this study.
Open AccessThis 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.
- Gadgil M, Bossert W. Life historical consequences of natural selection. Am Nat. 1971;104:1–24.View ArticleGoogle Scholar
- Begon M, Harper JL, Townsend CR. Ecology: individuals, populations and communities. London: Blackwell Scientific Publications; 1990.Google Scholar
- Roff DA. The evolution of life-histories: theory and analysis. New York: Chapman and Hall; 1992.Google Scholar
- Sorci G, Clobert J, Belichon S. Phenotypic plasticity of growth and survival in the common lizard Lacerta vivipara. J Anim Ecol. 1996;65:781–90.View ArticleGoogle Scholar
- Berven KA. The genetic basis of altitudinal variation in the wood frog Rana sylvatica. I. An experimental analysis of life-history traits. Evolution. 1982;36:962–83.View ArticleGoogle Scholar
- Berven KA, Gill DE. Interpreting geographic variation in life-history traits. Am Zool. 1983;23:85–97.View ArticleGoogle Scholar
- Adolph SC, Porter WP. Growth, seasonality and lizard life histories: age and size at maturity. Oikos. 1996;77:267–78.View ArticleGoogle Scholar
- Miaud C, Guyétant R, Elmberg J. Variation in life-history traits in the common frog Rana temporaria (Amphibia: Anura): a literature review and new data from the French Alps. J Zool (Lond). 1999;249:61–73.View ArticleGoogle Scholar
- Storz JF, Balasingh J, Bhat HR, Nathan PT, Doss DPS, Prakash AA, et al. Clinal variation in body size and sexual dimorphism in an Indian fruit bat, Cynopterus sphinx (Chiroptera: Pteropodidae). Biol J Linn Soc. 2011;72:17–31.View ArticleGoogle Scholar
- Morrison C, Hero JM. Geographic variation in life-history characteristics of amphibians: a review. J Anim Ecol. 2003;72:270–9.View ArticleGoogle Scholar
- Angilletta MJ, Oufiero CE, Leaché AD. Direct and indirect effects of environmental temperature on the evolution of reproductive strategies: an information-theoretic approach. Am Nat. 2006;168:123–35.View ArticleGoogle Scholar
- Gvoždík V, Moravec J, Kratochvíl L. Geographic morphological variation in parapatric Western Palearctic tree frogs, Hyla arborea and Hyla savignyi: are related species similarly affected by climatic conditions? Biol J Linn Soc. 2008;95:539–56.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Sinsch U, Marangoni F, Oromí N, Leskovar C, Sanuy D, Tejedo M. Proximate mechanisms determining size variability in natterjack toads. J Zool (Lond). 2010;281:272–81.Google Scholar
- Liao WB, 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.View ArticleGoogle Scholar
- Hjernquist MB, Söderman F, Jönsson KI, Herczeg G, Laurila A, Merilä J. Seasonality determines patterns of growth and age structure over a geographic gradient in an ectothermic vertebrate. Oecologia. 2012;170:641–9.View ArticlePubMedGoogle Scholar
- Liao WB. Evolution of sexual size dimorphism in a frog obeys the inverse of Rensch’s rule. Evol Biol. 2013;40:493–9.View ArticleGoogle Scholar
- Liao WB, Liu WC, Merilä J. Andrew meets Rensch: Sexual size dimorphism and the inverse of Rensch’s rule in Andrew’s toad (Bufo andrewsi). Oecologia. 2015;177:389–99.View ArticlePubMedGoogle Scholar
- Tinkle DW, Wilbur HM, Tilley SG. Evolutionary strategies in lizard reproduction. Evolution. 1970;24:55–74.View ArticleGoogle Scholar
- Räsänen K, Söderman F, Laurila A, Merilä J. Geographic variation in maternal investment: acidity affects egg size and fecundity in Rana arvalis. Ecology. 2008;89:2553–62.View ArticlePubMedGoogle Scholar
- Bergmann C. Über die Verhältnisse der Wärmeökonomie der Thiere zu ihrer Grösse. Gött Stud. 1847;1:595–708.Google Scholar
- Olalla-Tárraga MA, 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.View ArticleGoogle Scholar
- Meiri S, Dayan T. On the validity of Bergmann’s rule. J Biogeogr. 2003;30:331–51.View ArticleGoogle Scholar
- Ashton KG, Feldman CR. Bergmann’s rule in nonavian reptiles: turtles follow it, lizards and snakes reverse it. Evolution. 2003;57:1151–63.View ArticlePubMedGoogle Scholar
- Pincheira-Donoso D, Hodgson DJ, Tregenza T. The evolution of body size under environmental gradients in ectotherms: why should Bergmann’s rule apply to lizards? BMC Evol Biol. 2008;8:68.PubMed CentralView ArticlePubMedGoogle Scholar
- Mousseau TA. Ectotherms follow the converse to Bergmann’s rule. Evolution. 1997;51:630–2.View ArticleGoogle Scholar
- Ma X, Lu X, Merilä J. Altitudinal decline of body size in a Tibetan frog. J Zool (Lond). 2009;279:364–71.View ArticleGoogle Scholar
- Wells KD. The ecology and behaviour of amphibians. Chicago and London: University of Chicago; 2007.View ArticleGoogle Scholar
- Laugen AT, Laurila A, Jonsson KI, Soderman F, Merilä J. Do common frogs (Rana temporaria) follow Bergmann’s rule? Evol Ecol Res. 2005;7:717–31.Google Scholar
- Karl I, Fischer K. Altitudinal and environmental variation in lifespan in Copper butterfly Lycaena tityrus. Funct Ecol. 2009;23:1132–8.View ArticleGoogle Scholar
- Horvátová TH, Cooney CR, Fitze PS, Oksanen TA, Jelić D, Ghira I, et al. Length of activity season drives geographic variation in body size of a widely distributed lizard. Ecol Evol. 2013;3:2423–42.Google Scholar
- Roitberg ES, Kuranova VN, Bulakhova NA, Orlova VF, Eplanova GV, Zinenko OI, et al. Variation of reproductive traits and female body size in the most widely-ranging reptile species: testing the effects of reproductive mode, lineage, and climate. Evol Biol. 2013;40:420–38.PubMed CentralView ArticlePubMedGoogle Scholar
- Macey JR, Shulte II JA, Larson A, Fang Z, Wang Y, Tuniyev BS, et al. Phylogenetic relationships of toads in the Bufo bufo species group from the eastern escarpment of the Tibetan Plateau: a case of vicariance and dispersal. Mol Phylogenet Evol. 1998;9:80–7.View ArticlePubMedGoogle Scholar
- Frost DR. Amphibian species of the world: an online reference. New York: American Museum of Natural History; 2013.Google Scholar
- Fei L, Ye CY. The colour handbook of amphibians of Sichuan. Beijing: China Forestry Publishing House; 2001.Google Scholar
- Liao WB, Lu X, Jehle R. Altitudinal variation in maternal investment and trade-off between egg size and clutch size in the Andrew’s Toad. J Zool (Lond). 2014;293:84–91.View ArticleGoogle Scholar
- Liao WB. Elevational variation in the life-history of anurans in a subtropics montane forest of Sichuan, southwestern China. PhD thesis. Wuhan: Wuhan University; 2009.Google Scholar
- Castanet J, Smirina E. Introduction to the skeletochronological method in amphibians and reptiles. Ann Sci Nat Zool. 1990;11:191–6.Google Scholar
- Sinsch U. Review: Skeletochronological assessment of demographic life-history traits in amphibians. Herpetol J. 2015;25:5–13.Google Scholar
- Li ST, Wu X, Li DY, Lou SL, Mi ZP, Liao WB. Body size variation of Odorous Frog (Odorrana grahami) across altitudinal gradients. Herpetol J. 2013;23:187–92.Google Scholar
- Von Bertalanffy L. Quantitative laws in metabolism and growth. Quar Rev Biol. 1957;32:217–31.View ArticleGoogle Scholar
- Hemelaar AMS. Age, growth and other population characteristics of Bufo bufo from different latitudes and altitudes. J Herpetol. 1988;22:369–88.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.Google Scholar
- Angilletta MJ, Steury T, Sears M. Temperature, growth rate, and body size in ectotherms: fitting pieces of a life-history puzzle. Integr Comp Biol. 2004;44:498–509.View ArticlePubMedGoogle Scholar
- Smith DC. Adult recruitment in chorus frogs: effects of size and date of metamorphosis. Ecology. 1987;68:344–50.View ArticleGoogle Scholar
- Berven KA. Factors affecting population fluctuations in larval and adult stages of the wood frog (Rana sylvatica). Ecology. 1990;71:1599–606.View ArticleGoogle Scholar
- Komoroski MJ, Nagle RD, Congdon JD. Relationships of lipids to ovum size in amphibians. Phys Zool. 1998;71:633–41.View ArticleGoogle Scholar
- Gliwicz ZM, Guisande C. Family-planning in Daphnia: resistance to starvation in offspring born to mothers grown at different food levels. Oecologia. 1992;91:463–7.View ArticleGoogle Scholar
- Roff DA. Life-history evolution. Sunderland, Massachusetts, USA: Sinauer Associates; 2002.Google Scholar
- 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.View ArticleGoogle Scholar
- Charlesworth B. Evolution in age-structured populations. Cambridge: Cambridge University Press; 1994.View ArticleGoogle Scholar
- Schemske DW, Mittelbach GC, Cornell HW, Sobel JM, Kaustuv R. Is there a latitudinal gradient in the importance of biotic interactions? Ann Rev Ecol Evol Syst. 2009;40:245–69.View ArticleGoogle Scholar
- Lessard JP, Sackett TE, Reynolds WN, Fowler DA, Sanders NJ. Determinants of the detrital arthropod community structure: the effects of temperature, resources, and environmental gradients. Oikos. 2011;120:333–43.View ArticleGoogle Scholar
- Brown JH, Gillooly JF, Allen AP, van Savage M, West GB. Toward a metabolic theory of ecology. Ecology. 2004;85:1771–89.View ArticleGoogle Scholar
- Liao WB, Lu X. Proximate mechanisms leading to large male-mating advantage in the Andrew’s toad Bufo andrewsi. Behaviour. 2011;148:1087–102.View ArticleGoogle Scholar
- Laugen AT, Laurila A, Merilä J. Latitudinal and temperature dependent variation in embryonic development and growth in Rana temporaria. Oecologia. 2003;135:548–54.View ArticlePubMedGoogle Scholar
- Lindgren B, Laurila A. Proximate causes of adaptive growth rates: growth efficiency variation among latitudinal populations of Rana temporaria. J Evol Biol. 2005;18:820–8.View ArticlePubMedGoogle Scholar
- Alho JS, Herczeg G, Laugen AT, Räsänen K, Laurila A, Merilä J. Allen’s rule revisited: quantitative genetics of extremity length in the common frog along a latitudinal gradient. J Evol Biol. 2011;24:59–70.View ArticlePubMedGoogle Scholar
- Muir AP, Biek R, Thomas R, Mable BK. Local adaptation with high gene flow: temperature parameters drive adaptation to altitude in the common frog (Rana temporaria). Mol Ecol. 2014;23:561–74.PubMed CentralView ArticlePubMedGoogle Scholar