Skip to main content

Environmental causes of between-population difference in growth rate of a high-altitude lizard



Ectothermic animals living in cold (high latitude or high elevation) regions are predicted to grow slower due to limited thermal opportunities for activity and food resources than those living in warm regions. However, the Qinghai toad-headed lizards (Phrynocephalus vlangalii) grow faster and reach a larger adult size at a high-elevation site than at a low-elevation site. In this study, we aimed to identify the genetic and environmental causes of this between-population difference in growth rate by conducting mark-recapture and common garden experiments on juvenile growth rate, and investigating the thermal environment, lizard body temperature, potential prey availability at the two elevation sites.


Compared with low-elevation individuals, high-elevation juvenile lizards had higher growth rates in the field, but grew at similar rates in the laboratory. High-elevation lizards had higher active body temperatures than low-elevation lizards despite similar air temperatures in the period of field investigation. The high-elevation site had relatively more and larger preys than the low-elevation site.


Inter-population difference in growth rate of P. vlangalii may primarily result from developmental plasticity in response to the difference in environmental resources, rather than genetic differentiation. The higher growth rate of high-elevation lizards is likely associated with higher potential food availability and higher active body temperatures.


Growth rate is a critical life history trait influencing individual fitness [1, 2]. For example, growth rate may determine the age and body size at maturity, survival or reproductive success [3,4,5,6]. Growth rate of organisms often vary along a geographic (e.g. latitudinal and altitudinal) gradient [7,8,9]. Such geographic variation in growth rate may not only stem from genotypic difference among populations, but also be attributable to environmental factors such as temperature and food availability, particularly in ectotherms [5, 10,11,12,13,14]. To better understand the role of growth rate in generating geographic variations in other life history traits, such as body size, we need to identify the relative roles of extrinsic versus intrinsic factors generating geographic variation in these traits [15, 16].

Altitudinal variation in growth rate has been reported in a number of species [9, 17,18,19,20]. However, the underlying mechanisms have received much less attention [but see 2123]. Compared with the ectotherms living at low-elevation sites, the high-elevation ectotherms may grow slower due to the limited thermal opportunity for activity [6, 13, 24] or resource limitations [17, 25]. However, sometimes high-elevation ectotherms can grow faster than low-elevation ones due to greater food availability [21, 26, 27]. Those species in high-elevation (cold) regions may grow faster through the following ways that are not mutually exclusive. First, they may evolve fast-growing genotypes to compensate the disadvantage of a low growth rate due to cold environment [3, 15]. Second, environmental factors, including precipitation, food availability and sunshine duration, facilitate animals at colder sites to have higher growth rates than those at warmer sites [21]. Of those environmental factors that potentially induced growth rate variation, thermal regime and food availability are especially important because they may have strong proximate effects on the growth rate of animals [6, 14, 15, 24, 28,29,30]. High food availability and warm temperatures may enable ectotherms to grow fast [26, 31, 32]. In addition, most studies of altitudinal variation in reptile life history are within the altitudinal range of 150–2000 m [6, 13, 14, 21, 24, 27, 33,34,35], with few studies carried out in the regions above 2000 m in elevation [26, 36,37,38]. With an average elevation exceeding 4500 m, the Qinghai-Tibetan plateau region offers an opportunity for studying altitudinal variation in life history traits of organisms.

The Qinghai toad-headed lizard (Phrynocephalus vlangalii), which is widely distributed in the Qinghai-Tibetan Plateau with an elevation range from 2000 to 4500 m, provides an excellent model to study altitudinal variation of growth rate in reptiles. Our previous study indicated that P. vlangalii individuals at a high-elevation site (Maduo, 4250 m elevation) can grow faster and reach a larger adult size than those at a low-elevation site (Maqu, 2930 m elevation) [39]. In this study, we conducted a mark-recapture experiment in the field and a common-garden experiment in the laboratory to determine genetic vs environmental components of this between-population difference in juvenile growth rate. To identify the environmental causes of the between-population difference in growth rate, we further examined between-population differences in thermal environment, body temperature, diet composition, and potential food availability at these two sites separated by c.a. 1300 m elevation. We hypothesized that lizards from the high-elevation population would grow faster than those from the low-elevation population under identical laboratory conditions, if high-elevation lizards evolve fast-growing genotypes to compensate the disadvantage of harsh environment; or lizards from the high-elevation population would grow faster than those from the low-elevation population in the field rather than in the laboratory, if the between-population difference in lizard growth rate is induced primarily by environmental differences between the two sites (specifically, more food or thermal opportunities at the high-elevation site than at the low-elevation site).


Study species and study areas

Phrynocephalus vlangalii is a small ground-dwelling viviparous agamid lizard (up to 80 mm snout-vent length, SVL), and typically found in open spaces in arid or semi-arid regions covered by sparse vegetation in the Qinghai-Tibetan plateau. Many life-history traits of this species show large amounts of geographic variation [40,41,42]. Courtship occurs in May, and parturition occurs between mid-July and late August [43].

This study was conducted at a low-elevation site, Maqu (Gansu, western China, 34°00′N, 102°04′E, 2930 m elevation), and a high-elevation site, Maduo (Qinghai, western China, 34°55′N, 98°12′E, 4250 m elevation). Over a linear distance of 370 km, these two sites have distinct mean annual air temperatures (1.4 °C in Maqu vs − 1.7 °C in Maduo), rainfall (550 mm in Maqu vs 379 mm in Maduo), and sunshine duration (2594 h in Maqu vs 2842 h in Maduo). Toad-headed lizards are abundant in both study sites, and previous phylogenetic analysis has indicated that these two populations belong to a single lineage [44].

Growth of juvenile lizards in the field and laboratory

In mid-June (10th‒18th) and late August (22nd‒30th) of 2011 and 2012, a study plot of 4000 m2 in each site (80 × 50 m2) was visited 4 times. At each visit, active lizards were captured by hand whenever possible, weighed and measured for SVL, noting sex and reproductive condition. In the years of 2011 and 2012, a total of 370 juveniles (low-elevation: 58 in 2011 and 75 individuals in 2012; high-elevation: 107 in 2011 and 130 individuals in 2012) were captured and marked individually by toe-clipping upon first capture. After collecting the data, lizards were released immediately at their site of capture. Data from the mark-recapture study were used to estimate the growth rate in the field.

In mid-July of 2011, some gravid females captured in the field (29 from the low-elevation and 30 from the high-elevation) were transferred to our laboratory in Hangzhou Normal University. Females gave birth from late July to mid-August in the laboratory. A total of 79 newborns (51 from the low-elevation and 28 from the high-elevation) were collected, and housed in fifteen 50 × 40 × 30 cm terraria (3–4 individuals per terrarium) that placed in a temperature-controlled room set at 20 ± 2 °C, with individuals from different populations being assigned evenly to each aquarium. A 60 W light bulb was suspended above each terrarium allowing thermoregulation for 10 h daily. Small mealworms (larvae of Tenebrio molitor), and water enriched with vitamins and minerals were provided ad libitum throughout the experiment. Newborns were weighed and measured for SVL at birth, and at an age of 30 days to determine postnatal growth in a month.

Ambient temperature and lizard body temperature

From early June to late August of 2011, we placed 10 Thermocron iButton temperature loggers (DS1921, 711C, MAXIM Integrated Products/Dallas Semiconductor Ltd, Sunnyvale, USA) in each study site to record the ground surface temperature (Tg) in the natural habitat at an interval of 3 h. Each iButton was pressed into in a silvery-white balloon and placed at a 0.5 cm depth from the surface. Additionally, during a 1-week period in early August 2011 (1st‒7th), the active lizards at both study sites were captured by hand, and their body (cloacal) temperatures (Tb), Tg and air temperatures at 10 cm above the substrate (Ta) were recorded immediately after capture using a UT325 electronic thermometer (Uni-trend Group, Shanghai, China) to investigate between-population difference in active body temperature of lizards.

Diet composition and prey availability

In early August 2011 (1st‒7th), sixty-three adult lizards captured from the two study sites (15♀♀, 23 ♂♂ in the low-elevation population; 9♀♀, 16 ♂♂ in the high-elevation population) were immediately sacrificed for diet analysis. We dissected specimens, removed the stomachs, and placed the animal and its stomach in 10% formalin. Specimens were transferred to our laboratory in Hangzhou Normal University. Stomach contents were spread in a Petri dish and all prey items were identified to the taxonomic level of Order with the aid of a dissecting microscope in the laboratory.

Potential prey availability in the two sites was estimated by means of pitfall traps. Despite having a risk of biased estimation of prey availability, pitfall traps are considered as an effective technique to estimate ground-dwelling arthropods, which predominate in the diet composition of lizards. In the years of 2011 and 2012, we visited the study sites four times (early September 2011, mid-June, late July and late August 2012), and randomly placed thirty pitfalls (diameter = 20 cm, depth = 15 cm) throughout the study sites. All the trapped arthropods were collected after 48 h, and placed into wide-necked bottles with ethyl acetate and transferred to the laboratory for analysis. Arthropods captured in pitfall traps were identified to the taxonomic level of Order, counted and measured individually for body length and width. The volume (V) of each prey item was calculated using the equation of an ellipsoid: V = (4/3)πab2, where a = 1/2 body length and b = 1/2 body width. We used the total volume of prey as a measure of prey availability.

Data analysis

We calculated size-specific and mass-specific growth rates during the experiments using the formula ln (measurement2/measurement1)/(date2 − date1). Prior to analysis, the normality of distributions and homogeneity of variances in the data were tested using the Kolmogorov–Smirnov test and Bartlett’s test, respectively. Analysis of variance (ANOVA) or covariance (ANCOVA) was used to determine between-population differences in body and ambient temperature, food availability and individual growth rate of lizards. Log-likelihood ratio test (G-test) was used to determine between-site difference in recapture rate and the proportion of the individuals having empty stomachs. Values are presented as mean ± standard error (SE), and the significance level is set at α = 0.05.


Growth rate of juveniles in the field and laboratory

A total of 81 juveniles (30 from the low-elevation population and 51 from the high-elevation population) were recaptured at least once during the mark-recapture experiment. The proportion of recaptured individuals in total marked lizards did not differ between the two study sites (G = 0.50, df = 1, P = 0.817). Initial body size of recaptured lizards did not differ between sites (low-elevation vs high-elevation: SVL, 41.43 ± 0.73 vs 41.66 ± 0.39 mm, F1, 79 = 0.09, P = 0.762; mass, 2.78 ± 0.14 vs 2.80 ± 0.08 g, F1, 79 = 0.01, P = 0.934). Lizards at the high-elevation site increased their SVL (F1, 79 = 4.41, P = 0.039) and body mass (F1, 74 = 8.35, P < 0.01) more rapidly than those at the low-elevation site (Fig. 1).

Fig. 1

Mean values (+SE) of size- and mass-specific growth rates of juveniles in the field and newborns in the laboratory from two different elevation sites. The asterisks indicate significant differences (*P < 0.05, **P < 0.01)

Thirty-nine newborns (22 from the low-elevation population, 17 from the high-elevation population) died during the growth experiment. The proportion of survived individuals did not differ between the two study populations (G = 2.25, df = 1, P = 0.134). Initial body size of survived individuals from low-elevation population was greater than that from high-elevation population (mixed-model ANOVA with site of origin as a fixed factor and clutch as a random factor, low-elevation vs high-elevation: SVL, 28.8 ± 0.4 vs 26.7 ± 0.4 mm, F1, 18.3 = 6.03, P = 0.024; mass, 0.93 ± 0.05 vs 0.73 ± 0.05 g, F1, 19.1 = 3.15, P = 0.092). In the first month, the newborns did not differ in either size-specific growth rate (F1, 8.6 = 0.81, P = 0.394), or mass-specific growth rate (F1, 14.3 = 0.22, P = 0.645) between populations (Fig. 1).

Ambient temperature and lizard body temperature

The records from iButton temperature loggers showed that the mean daily ground surface temperatures (Tg) in the low-elevation and high-elevation site were 16.0, and 13.1°C, with ranges of 9.7–24.8 and 8.4–20.0 °C, respectively (Fig. 2). During the period of field investigation (from August 1st to 7th in 2011), lizard body temperatures (Tb) were significantly higher than air temperatures (Ta) at both low- and high-elevation sites (paired t-test: low-elevation, t = 20.65, df = 138, P < 0.001; high-elevation, t = 27.85, df = 147, P < 0.001). Lizard Tbs were higher than Tgs at high-elevation site (t = 7.58, df = 147, P < 0.001), but not at low-elevation site (t = 1.33, df = 138, P = 0.185) (Fig. 3). In the same time period, mean lizard Tbs (F1, 263 = 21.57, P < 0.001) and Tgs (F1, 263 = 8.05, P < 0.01) at high-elevation site were higher than those at low-elevation site, but mean Tas did not differ between the two sites (F1, 263 = 0.06, P = 0.802) (Fig. 3). Lizard Tbs were positively correlated with Tgs in both populations (low-elevation, r2 = 0.48, F1, 137 = 124.15, P < 0.001, Tb = 14.87 + 0.52Tg; high-elevation, r2 = 0.63, F1, 146 = 244.03, P < 0.001, Tb = 14.28 + 0.58Tg). Between-site difference in Tb was still evident after removing the effect of Tg (ANCOVA with Tg as the covariate, F1, 262 = 14.02, P < 0.001).

Fig. 2

Ground surface temperatures (from early June to late August) in natural habitats used by Phrynocephalus vlangalii at two different elevation sites

Fig. 3

Body temperatures of Phrynocephalus vlangalii (a), ground surface temperatures (b) and air temperatures at 10 cm above the substrate (c) during a 1-week period in August at two different elevation sites. Numbers above (or below) the error bars represent sample sizes

Diet composition and prey availability

Five lizards from the low-elevation population and one from the high-elevation population had no prey items in their stomachs. The proportion of individuals having empty stomachs did not differ significantly between the high-elevation (4.0%) and the low-elevation (13.2%) populations (G = 1.64, df = 1, P > 0.10). The wet masses of stomach contents were greater in high-elevation lizards (0.21 ± 0.02 g, n = 24) than those in low-elevation ones (0.15 ± 0.02 g, n = 33) (ANCOVA with body mass as the covariate, F1, 54 = 4.56, P = 0.037). The diet of P. vlangalii was primarily made up of arthropods. Coleoptera, Diptera, Lepidoptera and Hymenoptera (ants) were recorded from specimens in both populations. Lithobiomorpha and Opiliones were only recorded from specimens collected at the high-elevation site, while plant material was only found at the low-elevation site. Three of six major taxa (Coleoptera, Diptera and Lepidoptera) accounted for approximately 80% of the total number, and 90% of the volumetric contribution of food ingested by the lizards from both study sites (Table 1).

Table 1 Prey items found in stomach contents of adult Phrynocephalus vlangalii collected from two different elevation sites

The two study sites possess similar types of prey (Table 2). However, the proportion of some prey differed between the two sites. The high-elevation site had higher proportions of Lepidoptera and Opiliones, whereas the low-elevation site had higher proportions of Hymenoptera and Diptera (Table 2). The availability and size of the three major taxa (Coleoptera, Diptera and Lepidoptera) of prey for the lizard differed significantly between the two sites, with more and larger prey at the high-elevation site than at the low-elevation site (Table 3, Fig. 4).

Table 2 Potential prey items from pitfall traps placed in natural habitats used by Phrynocephalus vlangalii at two different elevation sites
Table 3 Results of repeated measures ANOVA (with site of origin as the between-subject factor and season as the within-subject factor) for the availability of potential prey in natural habitats used by Phrynocephalus vlangalii at two different elevation sites
Fig. 4

Mean values (+SE) of potential prey size and availability in natural habitats used by Phrynocephalus vlangalii at two different elevation sites. The asterisks indicate significant differences (***P < 0.001)


The climate of the Qinghai-Tibetan Plateau is characterized by strong solar radiation, low temperatures and limited precipitation [45], which might affect the growth of those animal species inhabiting this area. Our study demonstrated significant variations in thermal environment, food availability and juvenile growth rate between low- and high-elevation populations of P. vlangalii. Under identical laboratory conditions, however, juvenile lizards from the two populations showed a similar growth rate. These results indicate that divergent growth rate of P. vlangalii is primarily derived from plasticity induced by environmental factors rather than genetic differentiation.

Different genotypes might be favored by natural selection in different environments [46, 47]. Fast growing genotypes are often found in those resource-poor environments causing slow growth [7]. For example, temperate juveniles grew faster than tropical ones under identical laboratory conditions in the eastern water skink Eulamprus quoyii [3]. Similarly, genetic differences in growth rate have been documented using the common garden or reciprocal transplant experiment in other reptile populations [2, 14, 48, 49]. However, our laboratory growth experiment revealed that newborns of P. vlangalii from the two study populations had similar growth rates, which was inconsistent with our first prediction. Similarly, no significant genetic differentiation in growth rate was also reported in the common lizard, Zootoca vivipara [6].

Temperature and food quality are likely to be the most important environmental factors influencing lizard growth rates [13, 21]. Due to thermal constraints on body temperatures, lizard growth should be positively correlated with the amount of time over which body temperatures are suitable for activity. An increase in thermal opportunity for activity can increase growth rate of lizards [14, 26]. Due to relatively higher ambient temperatures, the low-elevation regions presumably offer longer potential daily and seasonal activity periods, and thereby tend to produce greater growth rates and prolonged growth periods [26]. However, P. vlangalii individuals at the low-elevation site did not necessarily have higher body temperatures than those at the high-elevation site (Fig. 3). In fact, our survey results in August showed that, in the same time period, Tas at the low-elevation site were not higher, but Tbs and Tgs were lower than at the high-elevation site. Unfortunately, we lacked the data on operative temperature here, so the Tgs recorded by iButton temperature loggers [> minimal recorded value of Tg (12 °C) when lizards were active] were used to estimate lizard Tbs over the active season. There was also no significant between-population difference in estimated Tb (t = 0.46, df = 414, P = 0.647). Solar radiation increases with increasing elevation [50]. Hence the higher intensity of solar radiation at the high-elevation site contributes to higher ground surface temperatures, and thus higher lizard body temperatures [51]. Moreover, high-elevation ectothermic animals can increase their active body temperatures in a more efficient way by using solar radiation heat sources, and/or by being active at low temperatures [52, 53]. For example, despite lower air temperature, high-elevation Liolaemus lizards had similar activity temperatures compared with low-elevation lizards [54]. Body size and coloration can affect thermoregulation efficiency of ectothermic animals [55, 56]. For example, some lizards living at high elevations have darker dorsal colorations to facilitate thermoregulation and warm faster than those at low elevations [55]. Phrynocephalus lizards in the Qinghai-Tibetan Plateau have an abdominal black-speckled area, which becomes larger with increased elevation. The black-speckled area allows lizards to gain heat more efficiently probably by absorbing reflected solar radiation from ground surfaces [57]. It would be of great interest to evaluate between-population differences in dorsal coloration and abdominal black-speckled area size and the contribution of the difference to thermoregulation efficiency of P. vlangalii in future studies. Overall, these findings together with a shorter annual sunshine duration at low-elevation site than at high-elevation site (Fig. 5, data from allow us to infer that the daily and seasonal activity period of low-elevation lizards are not necessarily longer than that of high-elevation lizards. Therefore, P. vlangalii from the high-elevation population might be able to use thermal opportunities more efficiently to facilitate individual growth. Nonetheless, we only quantified lizard thermoregulation in a shorter period (1 week in summer), further extensive studies on the thermoregulation of lizards from the two populations are needed to test our hypothesis on altitudinal variation in thermal effects on lizard growth.

Fig. 5

Monthly mean sunshine duration for the two study sites where Phrynocephalus vlangalii were collected (data from

Consistent with our prediction, potential prey availability was more abundant at the high-elevation site than at the low-elevation site, and high-elevation juveniles grew faster than low-elevation ones. Actually, higher growth rates for high-elevation P. vlangalii were also confirmed in another study using field mark-recapture data and von Bertalanffy growth parameters [39]. As suggested by some previous studies, the effect of food availability on lizard growth rate is more immediate [27]. The observed inter-population difference in growth rate of P. vlangalii also results from the variability in habitat food availability. Interestingly, high-elevation newborns were smaller than low-elevation ones, but high-elevation adults were reported to have larger sizes than low-elevation ones [39,40,41]. These results imply that the environmental factors determining growth rate are large enough to override the variation in newborn size. The higher food resources at high-elevation site could allow P. vlangalii juveniles to grow faster and thus reach larger adult size. On the other hand, food scarcity may restrict the growth potential of juveniles at low-elevation site. Larger newborn sizes in this environment should be favored by natural selection to compensate for lower growth rates [58, 59].

Lizard growth rates can be regulated by complex interaction between those environmental factors such as thermal opportunity for activity and food availability [21]. Consequently, it may be difficult to predict the latitudinal and altitudinal patterns of growth rate in different species of lizards. For example, the canyon lizards (Sceloporus merriami) at an intermediate-elevation site grow fastest despite food availability increasing with elevation [24], while the common lizards (Z. vivipara) raised at a low-elevation site grow faster than those at the high-elevation site, resulting from between-site differences in thermal environments [6]. Additionally, other internal and external factors can have a significant impact on growth rates of lizards, such as the length of daily activity and costs of activity [13], and metabolic expenditure or activity constraints imposed by predation risk [26, 38]. Understanding the contribution of these factors to lizard growth would be an interesting topic for future research.


Our study found that inter-population differences in growth rate of P. vlangalii were induced by the differences in environmental resources, rather than resulting from genetic differentiation. Although the slightly colder environments at the high-elevation site are thought to select for larger offspring that might have higher survival probability [35, 60, 61], smaller newborns were produced at the high-elevation site with more food resources. Fast growth, driven by high food availability, appeared to surpass the adverse effects of small body size and cold thermal environments. Contrarily, larger newborns produced at the low-elevation site ensured their survival, and compensated for slower growth rate in a low-food environment [21]. High-elevation lizards may not always have lower body temperatures and shorter activity period than low-elevation ones possibly due to more efficient ways to improve body temperature, higher solar radiation intensity and longer sunshine duration at the high-elevation site. Therefore, temperature might not be the major limiting factor that shapes altitudinal differences in growth rates. Our data indicated that potential food availability at high-elevation site allow juvenile P. vlangalii to have higher field growth rate and thereby larger body size at adulthood than at low-elevation site. Unfortunately, we only conducted a two-population comparison due to logistic difficulties in the current study. Further studies involving multiple populations along the altitudinal cline would be essential to elucidate the proximate and ultimate causes of altitudinal variations in lizard life history.


T b :

body (cloacal) temperature

T g :

ground surface temperature

T a :

air temperature


  1. 1.

    Shine R, Charnov EL. Patterns of survival, growth, and maturation in snakes and lizards. Am Nat. 1992;139:1257–69.

    Article  Google Scholar 

  2. 2.

    Niewiarowski PH. Energy budgets, growth rates, and thermal constraints: toward an integrative approach to the study of life-history variation. Am Nat. 2001;157:421–33.

    CAS  PubMed  Google Scholar 

  3. 3.

    Caley MJ, Schwarzkopf L. Complex growth rate evolution in a latitudinally widespread species. Evolution. 2004;58:862–9.

    Article  Google Scholar 

  4. 4.

    Clobert J, Oppliger A, Sorci G, Ernande B, Swallow JG, Garland T. Trade-offs in phenotypic traits: endurance at birth, growth, survival, predation and susceptibility to parasitism in a lizard, Lacerta vivipara. Funct Ecol. 2000;14:675–84.

    Article  Google Scholar 

  5. 5.

    Du WG, Warner DA, Langkilde T, Robbins T, Shine R. The roles of pre- and post-hatching growth rates in generating a latitudinal cline of body size in the eastern fence lizard (Sceloporus undulatus). Biol J Linn Soc. 2012;106:202–9.

    Article  Google Scholar 

  6. 6.

    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.

    Article  Google Scholar 

  7. 7.

    Conover DO, Duffy TA, Hice LA. The covariance between genetic and environmental influences across ecological gradients: reassessing the evolutionary significance of countergradient and cogradient variation. Ann NY Acad Sci. 2009;1168:100–29.

    Article  Google Scholar 

  8. 8.

    Roff DA. Life history evolution. Sunderland: Sinauer Associates, Inc.; 2002.

    Google Scholar 

  9. 9.

    Stearns SC. The evolution of life histories. Oxford: Oxford University Press; 1992.

    Google Scholar 

  10. 10.

    Andrews RM, Mathies T, Warner DA. Effect of incubation temperature on morphology, growth, and survival juvenile Sceloporus undulatus. Herpetol Monogr. 2000;14:420–31.

    Article  Google Scholar 

  11. 11.

    Goodman RM. Evidence of divergent growth rates among populations of the lizard Anolis carolinensis based on experimental manipulations of egg size. Popul Ecol. 2010;52:113–22.

    Article  Google Scholar 

  12. 12.

    Niewiarowski PH. Effects of supplemental feeding and thermal environment on growth rates of eastern fence lizards, Sceloporus undulatus. Herpetologica. 1995;51:487–96.

    Google Scholar 

  13. 13.

    Sinervo B, Adolph SC. Thermal sensitivity of growth rate in hatchling Sceloporus lizards: environmental, behavioral and genetic aspects. Oecologia. 1989;78:411–9.

    Article  Google Scholar 

  14. 14.

    Sinervo B, Adolph SC. Growth plasticity and thermal opportunity in Sceloporus lizards. Ecology. 1994;75:776–90.

    Article  Google Scholar 

  15. 15.

    Angilletta MJ, Steury TD, Sears MW. Temperature, growth rate, and body size in ectotherms: fitting pieces of a life-history puzzle. Integr Comp Biol. 2004;44:498–509.

    Article  Google Scholar 

  16. 16.

    Ferguson GW, Talent LG. Life-history traits of the lizard Sceloporus undulatus from two populations raised in a common laboratory environment. Oecologia. 1993;93:88–94.

    Article  Google Scholar 

  17. 17.

    Hodkinson ID. Terrestrial insects along elevation gradients: species and community responses to altitude. Biol Rev. 2005;80:489–513.

    Article  Google Scholar 

  18. 18.

    Kivelä SM, Välimäki P, Carrasco D, Mäenpää MI, Oksanen J. Latitudinal insect body size clines revisited: a critical evaluation of the saw-tooth model. J Anim Ecol. 2011;80:1184–95.

    Article  Google Scholar 

  19. 19.

    Sand H, Cederlund G, Danell K. Geographical and latitudinal variation in growth patterns and adult body size of Swedish moose (Alces alces). Oecologia. 1995;102:433–42.

    Article  Google Scholar 

  20. 20.

    Sears MW, Angilletta MJ. Body size clines in Sceloporus lizards: proximate mechanisms and demographic constraints. Integr Comp Biol. 2004;44:433–42.

    Article  Google Scholar 

  21. 21.

    Iraeta P, Monasterio C, Salvador A, Díaz JA. Mediterranean hatchling lizards grow faster at higher altitude: a reciprocal transplant experiment. Funct Ecol. 2006;20:865–72.

    Article  Google Scholar 

  22. 22.

    Karl I, Janowitz SA, Fischer K. Altitudinal life-history variation and thermal adaptation in the copper butterfly Lycaena tityrus. Oikos. 2008;117:778–88.

    Article  Google Scholar 

  23. 23.

    Sears MW, Angilletta MJ. Life history variation in the sagebrush lizard (Sceloporus graciosus): phenotypic plasticity or local adaptation? Ecology. 2003;84:1624–34.

    Article  Google Scholar 

  24. 24.

    Grant BW, Dunham AE. Elevational covariation in environmental constraints and life histories of the desert lizard Sceloporus merriami. Ecology. 1990;71:1765–76.

    Article  Google Scholar 

  25. 25.

    Chown SL, Klok CJ. Altitudinal body size clines: latitudinal effects associated with changing seasonality. Ecography. 2003;26:445–55.

    Article  Google Scholar 

  26. 26.

    Sears MW. Geographic variation in the life history of the sagebrush lizard: the role of thermal constraints on activity. Oecologia. 2005;143:25–36.

    Article  Google Scholar 

  27. 27.

    Iraeta P, Salvador A, Díaz JA. Life-history traits of two Mediterranean lizard populations: a possible example of countergradient covariation. Oecologia. 2013;172:167–76.

    Article  Google Scholar 

  28. 28.

    Angilletta MJ, Niewiarowski PH, Dunham AE, Leaché AD, Porter WP. Bergmann’s clines in ectotherms: illustrating a life-history perspective with Sceloporine lizards. Am Nat. 2004;164:E168–83.

    Article  Google Scholar 

  29. 29.

    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:E123–35.

    Article  Google Scholar 

  30. 30.

    Lourdais O, Bonnet X, Shine R, DeNardo D, Naulleau G, Guillon M. Capital-breeding and reproductive effort in a variable environment: a longitudinal study of a viviparous snake. J Anim Ecol. 2002;71:470–9.

    Article  Google Scholar 

  31. 31.

    Horváthová T, Cooney CR, Fitze PS, Oksanen TA, Jelić D, Ghira I, Uller T, Jandzik D. Length of activity season drives geographic variation in body size of a widely distributed lizard. Ecol Evol. 2013;3:2424–42.

    Article  Google Scholar 

  32. 32.

    Yamahira K, Conover DO. Intra- vs. interspecific latitudinal variation in growth: adaptation to temperature or seasonality? Ecology. 2002;83:1252–62.

    Article  Google Scholar 

  33. 33.

    Buckley CR, Irschick DJ, Adolph SC. The contributions of evolutionary divergence and phenotypic plasticity to geographic variation in the western fence lizard, Sceloporus occidentalis. Biol J Linn Soc. 2010;99:84–98.

    Article  Google Scholar 

  34. 34.

    Rodríguez-Díaz T, Braña F. Altitudinal variation in egg retention and rates of embryonic development in oviparous Zootoca vivipara fits predictions from the cold-climate model on the evolution of viviparity. J Evol Biol. 2012;25:1877–87.

    Article  Google Scholar 

  35. 35.

    Rohr DH. Demographic and life-history variation in two proximate populations of a viviparous skink separated by a steep altitudinal gradient. J Anim Ecol. 1997;66:567–78.

    Article  Google Scholar 

  36. 36.

    Leache AL, Helmer DA, Moritz C. Phenotypic evolution in high-elevation populations of western fence lizards (Sceloporus occidentalis) in the Sierra Nevada Mountains. Biol J Linn Soc. 2010;100:630–41.

    Article  Google Scholar 

  37. 37.

    Mathies T, Andrews RM. Thermal and reproductive biology of high and low elevation populations of the lizard Sceloporus scalaris: implications for the evolution of viviparity. Oecologia. 1995;104:101–11.

    CAS  Article  Google Scholar 

  38. 38.

    Sears MW. Resting metabolic expenditure as a potential source of variation in growth rates of the sagebrush lizard. Comp Biochem Physiol A. 2005;140:171–7.

    Article  Google Scholar 

  39. 39.

    Lu HL, Xu CX, Jin YT, Hero JM, Du WG. Proximate causes of altitudinal differences in body size in an agamid lizard. Ecol Evol. 2018;8:645–54.

    Article  Google Scholar 

  40. 40.

    Jin YT, Liu NF, Li JL. Elevational variation in body size of Phrynocephalus vlangalii in the North Qinghai-Xizang (Tibetan) Plateau. Belg J Zool. 2007;137:197–202.

    Google Scholar 

  41. 41.

    Li JQ, Zhou R, Liu NF. Life-history variation among three populations of the toad-headed lizard Phrynocephalus vlangalii along an elevation gradient on the northeastern Tibetan Plateau. Herpetol J. 2014;24:17–23.

    Google Scholar 

  42. 42.

    Jin YT, Liu NF. Altitudinal variation in reproductive strategy of the toad-headed lizard, Phrynocephalus vlangalii in North Tibet Plateau (Qinghai). Amphibia-Reptilia. 2007;28:509–15.

    Article  Google Scholar 

  43. 43.

    Wu YY, Fu JZ, Yue BS, Qi Y. An atypical reproductive cycle in a common viviparous Asia agamid Phrynocephalus vlangalii. Ecol Evol. 2015;5:5138–47.

    Article  Google Scholar 

  44. 44.

    Jin YT, Brown RP, Liu NF. Cladogenesis and phylogeography of the lizard Phrynocephalus vlangalii (Agamidae) on the Tibetan plateau. Mol Ecol. 2008;17:1971–82.

    CAS  Article  Google Scholar 

  45. 45.

    Lin ZY, Zhao XY. Spatial characteristics of changes in temperature and precipitation of the Qinghai-Xizang (Tibet) Plateau. Sci China Ser D. 1996;39:442–8.

    Google Scholar 

  46. 46.

    Endler J. Natural selection in the wild. Princeton: Princeton University Press; 1986.

    Google Scholar 

  47. 47.

    Irschick DJ, Meyers JJ. An analysis of the relative roles of plasticity and natural selection in the morphology and performance of a lizard (Urosaurus ornatus). Oecologia. 2007;153:489–99.

    Article  Google Scholar 

  48. 48.

    Bronikowski AM. Experimental evidence for the adaptive evolution of growth rate in the garter snake Thamnophis elegans. Evolution. 2000;54:1760–7.

    CAS  Article  Google Scholar 

  49. 49.

    Niewiarowski PH, Roosenburg WM. Reciprocal transplant reveals sources of geographic variation in growth rates of the eastern fence lizard, Sceloporus undulatus. Ecology. 1993;74:1992–2002.

    Article  Google Scholar 

  50. 50.

    Heinl M, Leitinger G, Tappeiner U. Diurnal surface temperature regimes in mountain environments. Phys Geog. 2013;33:344–59.

    Article  Google Scholar 

  51. 51.

    Carrascal LM, López P, Martín J, Salvador A. Basking and antipredator behaviour in a high altitude lizard: implications of heat-exchange rate. Ethology. 1992;92:143–54.

    Article  Google Scholar 

  52. 52.

    Martín J, López P, Carrascal LM, Salvador A. Adjustment of basking postures in the high-altitude Iberian rock lizard (Lacerta monticola). Can J Zool. 1995;73:1065–8.

    Article  Google Scholar 

  53. 53.

    Valdecantos S, Martínez V, Lobo F, Cruz FB. Thermal biology of Liolaemus lizards from the high Andes: being efficient despite adversity. J Therm Biol. 2013;38:126–34.

    Article  Google Scholar 

  54. 54.

    Marquet PA, Ortíz JC, Bozinovié F, Jaksié FM. Ecological aspects of thermoregulation at high altitudes: the case of andean Liolaemus lizards in northern Chile. Oecologia. 1989;81:16–20.

    CAS  Article  Google Scholar 

  55. 55.

    Reguera S, Zamora-Camacho FJ, Moreno-Rueda G. The lizard psammodromus algirus (squamata lacertidae) is darker at high altitudes. Biol J Linn Soc. 2014;112:132–41.

    Article  Google Scholar 

  56. 56.

    Geen MRS, Johnston GR. Coloration affects heating and cooling in three color morphs of the Australian blue tongue lizard, Tiliqua scincoides. J Therm Biol. 2014;43:54–60.

    Article  Google Scholar 

  57. 57.

    Jin YT, Tong HJ, Zhang KL. The impact of phenotypic characteristics on thermoregulation in a cold-climate agamid lizard, Phrynocephalus guinanensis. Asian Herpetol Res. 2016;73:210–9.

    Google Scholar 

  58. 58.

    Ferguson GW, Fox SF. Annual variation of survival advantage of large juvenile side-blotched lizards, Uta stansburiana—its causes and evolutionary significance. Evolution. 1984;38:342–9.

    PubMed  Google Scholar 

  59. 59.

    Lorenzon P, Clobert J, Massot M. The contribution of phenotypic plasticity to adaptation in Lacerta vivipara. Evolution. 2001;55:392–404.

    CAS  Article  Google Scholar 

  60. 60.

    Sinervo B. The evolution of maternal investment in lizards: an experimental and comparative analysis of egg size and its effects on offspring performance. Evolution. 1990;44:279–94.

    Article  Google Scholar 

  61. 61.

    Sinervo B, Doughty P, Huey RB, Zamudio K. Allometric engineering: a causal analysis of natural selection on offspring size. Science. 1992;258:1927–30.

    CAS  Article  Google Scholar 

Download references

Authors’ contributions

WGD conceived and designed the study. HLL and CXX conducted field and laboratory work. HLL, ZGZ and WGD analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.


We are grateful to C.-Q. Jiang and J.-M. Hero for their assistance in the field and comments on the manuscript.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article.

Consent for publication

Not applicable.

Ethics approval and consent to participate

This study was conducted in accordance to the current Chinese laws on animal welfare and research. The protocol and study were approved by the Animal Ethics Committee at the Institute of Zoology, Chinese Academy of Sciences (Permit Number: IOZ14001).


This work was supported by grants from the Strategic Priority Research Program of Chinese Academy of Sciences (XDA20050201, XDB31000000), and China’s Biodiversity Observation Network (Sino-BON).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information



Corresponding author

Correspondence to Wei-Guo Du.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Lu, HL., Xu, CX., Zeng, ZG. et al. Environmental causes of between-population difference in growth rate of a high-altitude lizard. BMC Ecol 18, 37 (2018).

Download citation


  • Phrynocephalus vlangalii
  • Elevational variation
  • Growth rate
  • Food availability
  • Activity period