pmc logo imageJournal ListSearchpmc logo image
Logo of bmcebBioMed Central web siteReference to the article.Search.Manuscript submission.Registration.Journal front page.
BMC Evol Biol. 2008; 8: 68.
Published online 2008 February 27. doi: 10.1186/1471-2148-8-68.
PMCID: PMC2268677
The evolution of body size under environmental gradients in ectotherms: why should Bergmann's rule apply to lizards?
Daniel Pincheira-Donoso,1 David J Hodgson,1 and Tom Tregenzacorresponding author1
1Centre for Ecology and Conservation, School of Biosciences, University of Exeter, Cornwall Campus, Penryn TR10 9EZ, Cornwall, UK
corresponding authorCorresponding author.
Daniel Pincheira-Donoso: D.PincheiraDonoso/at/exeter.ac.uk; David J Hodgson: D.J.Hodgson/at/exeter.ac.uk; Tom Tregenza: T.Tregenza/at/exeter.ac.uk
Received August 13, 2007; Accepted February 27, 2008.
Abstract
Background
The impact of environmental gradients on the evolution of life history traits is a central issue in macroecology and evolutionary biology. A number of hypotheses have been formulated to explain factors shaping patterns of variation in animal mass. One such example is Bergmann's rule, which predicts that body size will be positively correlated with latitude and elevation, and hence, with decreasing environmental temperatures. A generally accepted explanation for this phenotypic response is that as body mass increases, body surface area gets proportionally smaller, which contributes to reduced rates of heat-loss. Phylogenetic and non-phylogenetic evidence reveals that endotherms follow Bergmann's rule. In contrast, while previous non-phylogenetic studies supported this prediction in up to 75% of ectotherms, recent phylogenetic comparative analyses suggest that its validity for these organisms is controversial and less understood. Moreover, little attention has been paid to why some ectotherms conform to this rule, while others do not. Here, we investigate Bergmann's rule in the six main clades forming the Liolaemus genus, one of the largest and most environmentally diverse genera of terrestrial vertebrates. A recent study conducted on some species belonging to four of these six clades concluded that Liolaemus species follow Bergmann's rule, representing the only known phylogenetic support for this model in lizards. However, a later reassessment of this evidence, performed on one of the four analysed clades, produced contrasting conclusions.
Results
Our results fail to support Bergmann's rule in Liolaemus lizards. Non-phylogenetic and phylogenetic analyses showed that none of the studied clades experience increasing body size with increasing latitude and elevation.
Conclusion
Most physiological and behavioural processes in ectotherms depend directly upon their body temperature. In cold environments, adaptations to gain heat rapidly are under strong positive selection to allow optimal feeding, mating and predator avoidance. Therefore, evolution of larger body size in colder environments appears to be a disadvantageous thermoregulatory strategy. The repeated lack of support for Bergmann's rule in ectotherms suggests that this model should be recognized as a valid rule exclusively for endotherms.
Background
Geographical variation in environmental conditions is a major ecological factor involved in evolutionary diversification [1,2]. Since thermal regimes are particularly sensitive to latitude and altitude, geographical location imposes profound selection on organisms' metabolism, morphology and behaviour [3-5], leading to covariation between phenotypic traits and geographical gradients [5-7]. Body size is known to exhibit substantial variation in relation to thermal differences among habitats [4,8,9]. However, factors other than environmental temperature (e.g. sexual selection, predation) may also impose selection on body size [2,10-16]. As a result, models predicting patterns of evolutionary change in body mass in response to thermal variation are controversial [17-20]. One such example is Bergmann's rule [18], which suggests that species body size increases with increasing latitude and elevation, and hence, with decreasing environmental temperatures [6,21-23].
Different hypotheses have attempted to elucidate the causal factors promoting the pattern predicted by Bergmann's rule [18,24-26]. Potential explanations have focused on heat-conservation strategies [4,6,20,27], later maturation to larger body size [28,29], phylogenetic constraints [4,30], interspecific variation in migration [4], differential resistance to starvation [31-33], and effects on somatic cell sizes [18,26,34,35]. So far, the heat-conservation hypothesis, based on the fact that increases in body volume lead to decreases in relative surface area appears to be the most likely explanation [27].
Since the formulation of Bergmann's rule more than 150 years ago [21], several studies have explored the universality of its predictions across different lineages. In general, the validity of this rule appears to depend on the thermoregulatory physiology of the studied model systems. For endotherms, most evidence has supported Bergmann's rule [36-38]. This is presumably because the benefits of reducing heat-loss rates through larger body size are advantageous for cold climate organisms that maintain optimal body temperatures by metabolic generation of heat. However, the relevance of Bergmann's rule for ectotherms is less obvious, and supportive evidence is elusive [6,20,31,39].
Bergmann's rule in ectotherms
More than 99% of species are ectotherms [40]. Consequently, no biological prediction can be considered universal if it is not supported by these organisms. The first tests to evaluate Bergmann's rule in ectotherms [41,42] claimed that up to 75% of studied species support its predictions. However, these studies were based on simplistic statistical approaches, when the importance of conducting phylogenetic comparative studies was not yet appreciated [6,31,43,44]. Unsurprisingly, the development of phylogenetic comparative analyses to test Bergmann's rule over the last decade has provided divergent lines of evidence. While a series of studies confirmed previous findings [26,27,31,39,45], others were more equivocal [18,24,26,29-31]. For example, Ray [42] gathered data from previously published studies to evaluate Bergmann's rule. Regarding fishes, this author concluded that it is "obeyed by a great number of fishes as shown by numerous reports in the literature". Curiously, no citations were provided to support this claim [29,42]. More recently however, in a rigorous study conducted on more than 600 fish populations of different freshwater species from North America, Belk and Houston [29] observed that these ectotherms tend to reverse Bergmann's rule.
With regard to tetrapod ectotherms (i.e. amphibians and reptiles), evidence is similarly controversial [31]. A number of phylogenetic studies have shown that groups of primarily aquatic and tropical-subtropical species follow Bergmann's rule. For example, some anurans, and most urodeles (salamanders and newts) and turtles exhibit a negative relationship between body mass and environmental temperatures [30,31,41]. On the other hand, debate has progressively intensified in relation to squamate reptiles (lizards, snakes). Non-phylogenetic [41] and phylogenetic [30,31] studies have revealed that squamates exhibit only a weak tendency to conform to this model. Nevertheless, research conducted on some widespread groups has concluded that these reptiles can exhibit body size trends predicted by Bergmann's rule [20,27,46]. For example, a recent phylogenetic study conducted on Liolaemus lizards supported this temperature-size model [27]. Cruz et al. [27] observed that species of the clade Liolaemus boulengeri [11] show larger body size at higher latitudes and elevations. An additional integrative analysis including species belonging to congeneric clades led Cruz et al. [27] to conclude that "the strong, positive, size-latitude relationship in the L. boulengeri clade apparently accounted for the pattern observed for the entire dataset". However, after enlarging the sample of taxa belonging to one of the studied clades (clade boulengeri), Pincheira-Donoso et al. [39] observed that these species do not support Bergmann's rule.
Among factors that might explain the disparity observed across studies testing Bergmann's rule in ectotherms, the limited availability of phylogenetic studies is perhaps the most obvious [29,31,47]. Only scarce or generalized tests have been conducted on ectotherms, and most of this research focuses on the analyses of few species per clade (e.g. a single species representing clades consisting of more than 100 taxa) [20,31]. Moreover, studies focused on analysing explicit patterns of body size variation in the species of the same clade in response to continuous geographical gradients are still rare [27,29,46]. Recent criticism suggests that a rigorous test of Bergmann's rule should focus on species belonging to a monophyletic clade exhibiting substantial variation in patterns of body size and occurring in a broad geographical area encompassing a wide range of environmental conditions [6,39,46].
Lizards of the South American genus Liolaemus offer unique opportunities to test predictions concerning the impact of selection on traits recognized as labile to variation in environmental temperatures, such as body size [11,27,48]. This clade represents one of the most diverse known amniote lineages. Consisting of more than 190 named species, Liolaemus occurs in the widest range of environmental conditions recorded for any lizard genus [49-53]. These iguanians range from tropical-subtropical areas in Brazil and Peru, and the Atacama Desert in Chile, to austral Patagonia in Tierra del Fuego, the southernmost place where reptiles have been recorded [53-59]. The altitudinal distribution of Liolaemus is also one of the broadest known among squamate reptiles, occurring from sea level to over 5000 m in the Andes range [48,50,57,60,61]. These biological features satisfy all of the requirements recognized as essential for a model group employed to test predictions concerning evolutionary radiations [6,39,46].
Here, we investigate the effect that continuous variation in environmental temperatures imposes on the evolution of body size among species of the Liolaemus genus, using a comparative approach. We studied a set of more than 120 species (see additional file 1: Supplementary table), 63 of which are included in an explicit phylogenetic hypothesis (Fig. 1). This species sample represents almost the entire biogeographical, ecological and morphological diversity known for these lizards. We aim to test the hypothesis that increasing latitude and elevation (and therefore, decreasing environmental temperatures) are associated with larger body size [6,21]. We suggest as an alternative hypothesis that large body size is disadvantageous for ectotherms in cold-climates, because it demands longer time basking to achieve optimal metabolic temperatures. Therefore, we expect to find weak or no evidence in support of Bergmann's rule.
Figure 1Figure 1
Phylogenetic relationships for 63 of the 126 Liolaemus taxa included in this study. The clades lineomaculatus (lnm) and montanus (mnt) are abbreviated.
Results
Body size patterns under environmental gradients
We found non significant relationships between environmental gradients and body size in most of the tests performed. In none of the studied clades was Bergmann's rule supported (Table 1). For both methods of data analysis, significantly high or low standardized residual (z-scores) values revealed the presence of outliers which might affect the accuracy of the models. For independent clade analyses, extreme z-scores (see methods) were found in the clades montanus (one case with z-scores > 2.05) and Donosolaemus-magellanicus (one case with z-scores < -2.26). The relationship between independent contrasts for body size and adjusted latitudinal midpoint (ALM, latitude adjusted for altitude, see below) revealed the presence of five outliers affecting the model (over 8.06% of the cases with z-scores > 2.12, and < -2.21). The observed differences in the correlation and regression analyses including and excluding outliers are detailed below.
Table 1Table 1
Results of least squares regression analyses conducted on clades separately (non-phylogenetic) and on independent contrasts (phylogenetic) in the genus Liolaemus. In the clades Donosolaemus-magellanicus and montanus, these analyses were conducted including (more ...)
Linear bivariate regression analyses showed that body size does not vary predictably with adjusted latitudinal midpoint (ALM) in the genus Liolaemus. For analyses conducted on clades separately, low proportions of variance were explained in the groups chiliensis, fitzingerii, wiegmannii and lineomaculatus (Table 1). Regression analyses including outliers showed non-significant relationships between SVL and ALM in the clades montanus and Donosolaemus-magellanicus. When excluding outliers, we found a significant negative relationship between SVL and ALM for the clade montanus and a non-significant relationship in the clade Donosolaemus-magellanicus (Table 1; Fig. 2).
Figure 2Figure 2
Linear regressions of raw data for ln(snout-vent length) against adjusted latitudinal midpoint in the six main clades forming the Liolaemus genus. Top: regressions for the clades chiliensis (black), Donosolaemus-magellanicus including outliers (blue; (more ...)
Linear regression analyses conducted on phylogenetic independent contrasts (through the origin) revealed that ALM cannot predict body size when including outliers (Table 1; Fig. 3a), and when excluding outliers (Table 1; Fig. 3b).
Figure 3Figure 3
Linear regression analyses (through the origin) of phylogenetically independent contrasts (IC) for ln(snout-vent length) against adjusted latitudinal midpoint in the entire dataset of Liolaemus species for which phylogenetic information was available. (more ...)
Effects of sample size reduction
Non-phylogenetic linear regressions conducted on the 63 species for which phylogenetic information was available (see Additional file 1: Supplementary table), revealed that reduction of sample size for phylogenetic analyses (from 126 to 63 species; see methods) does not produce qualitative differences in the results. As observed in the whole dataset (126 species), body size does not vary predictably with adjusted latitudinal midpoint (ALM) in any of the studied groups. For these 63 species, non-significant proportions of variance were explained in the clades chiliensis (R2 = 0.063, F1,23 = 1.55, P = 0.23), Donosolaemus-magellanicus (R2 = 0.14, F1,5 = 0.78, P = 0.42), fitzingerii (R2 = 0.17, F1,16 = 3.22, P = 0.092), lineomaculatus (R2 = 0.79, F1,1 = 3.66, P = 0.31), montanus (R2 = 0.17, F1,1 = 0.21, P = 0.73) and wiegmannii (R2 = 0.16, F1,5 = 0.92, P = 0.38). The analysis on standardized residuals confirmed that none of the regression models suffered from oultiers.
Discussion
This analysis of Bergmann's rule is unique in including a representative proportion of the geographical, ecological, morphological and phylogenetic diversity of one of the most species rich terrestrial ectotherm vertebrate genera, Liolaemus lizards [see [6,39,46]]. Our results fail to support Bergmann's rule. Non-phylogenetic and phylogenetic analyses showed that increasing latitude and elevation do not predict increasing body size in Liolaemus species. These results further contradict the only known phylogenetic evidence in favour of Bergmann's rule reported in squamate reptiles [[27], see also [39]], and a recent large-scale analysis conducted on assemblages of European lizards supporting this macroecological model [20].
Thermal gradients, thermoregulatory physiology and body size
Research on Bergmann's rule has inspired debate about why increasing body size should be advantageous for cold climate species, and why its predictions are observed in some groups but not in others [6,31,39]. For example, it has been suggested that the anatomical characteristics of the skin (e.g. density, or epidermal covering such as feathers and fur) in groups as different as birds, mammals and turtles make large body size advantageous in colder environments and that the absence of these skin characteristics in squamate reptiles explains why they generally fail to follow Bergmann's rule [6,31,62]. In contrast, much less reference has been made to the role that thermoregulatory physiology of organisms (endothermy and ectothermy) can play during the process of body size evolution [39]. Indeed, an important part of the available discussions focuses on the benefits of increasing body size for heat retention in cold climate species. However, since conservation of body heat in ectotherms is only possible once they have optimally achieved it (which is metabolically resolved in endotherms), it appears to be simplistic to exclude the fact that increasing body mass also reduces the rates of heat gain [31], which may be as critical as the need for retaining heat [6,32]. This is illustrated, for example, by the observation that gravid females in several lizard species tend to bask more often or for longer than non-gravid females [63-66].
It might be argued that lizard body size and heating rates do not covary in a linear fashion, as these organisms may adjust (plastically or genetically) some physiological and behavioural traits to increase heating rates in cold environments [5,7,20,67,68]. For example, it has been observed that physiological adaptations may increase heating rates up to 17% in cold climate lizards [20,68]. Likewise, behavioural adjustments such as selection of basking sites less exposed to wind [3] or modifications in the postures and body orientation to the sun may contribute to gain heat more rapidly [3,69,70]. However, even though these adaptations contribute to increase heating rates in lizards inhabiting cold environments, their efficiency has limits. Indeed, a number of studies have shown that annual and daily activity, and thermoregulatory processes in cold climate lizards may vary substantially in comparison to warm climate lizards. For example, it has been observed that high elevation (e.g. over 4000 m) lizards tend to exhibit considerably shorter daily activities, often no more than four or five hours a day [71-73], than species distributed in lower latitudes and elevations, which may remain active for over twelve hours per day [57,61,72,74-76]. One of the most plausible factors to explain this pattern is that mean daily temperatures are lower, hot hours per day are fewer, and warm seasons are shorter in cold climates [22,71,77-79]. Therefore, even in presence of the above mentioned physiological and behavioural adjustments, cold environments restrict severely the patterns of activity in lizards, as a consequence of their failure to overcome the selective pressures imposed by low temperatures. This suggests that the thermoregulatory limitations determined by the ectothermal condition of lizards might be one of the main (if not the main) factors constraining the maximum attainable limits of body size in cold climate species.
Large body mass in cold climate lizards would have dramatic consequences for ecological performance (imposed by natural selection) and reproductive success (imposed by sexual selection), because in conditions of suboptimal metabolic temperatures most physiological functions (and hence, behavioural responses) occur at suboptimal rates [3,80]. For example, performance at prey capture, predator evasion, endurance, digestion, mate courtship, sperm production and conversion of lipids is substantially reduced at suboptimal body temperatures [3,66,80-83]. Also, basking for longer periods may increase the risk of predation by diurnal hunters [66]. Consequently, large body mass in cold climate lizards would be disadvantageous. Hence, the hypothesis of heat conservation cannot be accepted (at least in these organisms), as it predicts increasing body size in low-temperature environments. These claims are supported by a number of previous studies conducted in lizards [30,31,39].
Finally, it might also be argued that large ectotherms can maintain higher constant minimum body temperatures than do smaller species. However, it is unlikely that even relatively large cold climate lizards are able to retain heat overnight. Indeed, it has been observed in several species [84] that during the initial period of basking, the body temperature of lizards is very similar to the environmental temperature. Additionally, field observations conducted at high latitudes and elevations in South America (DPD, unpublished data) reveal that different sized Liolaemus lizards exhibit similar temperatures before initiating basking.
Trophic niche evolution and body size
The evolution of trophic niches in lizards may provide additional evidence to support the idea that larger body size is disadvantageous for species living in cold environments. Although most lizards feed on other animals, omnivory and herbivory have evolved independently many times among these squamates [80,85]. Plant matter is a low energy food for lizards [3]. Therefore, omnivory and herbivory are advantageous for these reptiles only when certain morphological and thermoregulatory conditions are met. In general, it has been observed that plant consumption has evolved primarily in species occurring in warm habitats, where vegetation is extremely abundant (e.g. tropical forests), or where abundance of animal prey is extremely low (e.g. deserts, small islands) [11,80,86,87]. In these environments, lizards can attain large body size and high body temperatures, both traits considered essential requirements for efficient digestion of plant matter [3,88-90]. A large body allows for a more voluminous gut, and high body temperature provides optimal conditions for endosymbionts (bacteria and protozoa) specialized in breaking down otherwise indigestible plant matter [3,9,91-93]. Consequently, according to these rules of herbivory [3,11], most lizards in which plant consumption has evolved are expected to reverse Bergmann's rule, because larger species tend to occur in warm environments.
Contrary to predictions, in a recent study Espinoza et al. [11] found that omnivory and herbivory have not only evolved in warm-climate and large lizards, but also in small species occurring in cold environments [94]. These authors observed that a large number of small Liolaemidae species (which includes Liolaemus) restricted to high latitudes or elevations, have substantial proportions of plant matter in their diets (i.e. 11%–100% of digestive volume content). Although the body size pattern observed in omnivorous and herbivorous Liolaemidae species differs significantly from the pattern known for plant-consumer lizards occurring in warm climates, the explanation may be exactly the same. Since digestion of plant matter requires the highest possible body temperatures [11,93], evolution of large body size in cold-climate Liolaemidae species would decrease the thermoregulatory efficiency needed to achieve optimal metabolic temperatures. This in turn would affect negatively the digestion of plant matter. Therefore, in the case of Liolaemidae species occurring at high latitudes and elevations, evolution of large size would be disadvantageous, challenging patterns predicted by Bergmann's rule.
Why do some ectotherms appear to follow Bergmann's rule?
Thermal environment differs from some other sources of selection in that it is a dominant selective force at low temperatures [31,39], but at higher temperatures, ectotherms are essentially freed from thermal constraints and hence other factors may become dominant in determining body size. These include fecundity selection [10,13,14], reduction of mass-specific energy requirements [3], greater physiological homeostasis [9], and reduction of mortality rates by predation [2,95-97]. Therefore, ectotherm clades whose geographical distributions are primarily restricted to tropical or subtropical ecosystems, such as continental turtles [30,31], may exhibit increasing body size with increasing latitude and moderate elevations. This pattern probably reflects a response to selective factors other than the thermal environment which is essentially benign throughout the geographic range of the group. Alternatively, observations of Bergmann's clines in some ectotherm vertebrates (i.e. turtles, salamanders) may be explained by the large taxonomic scale at which studies have been conducted [30,31] which could mean that different patterns of body size variation at an interspecific scale might be obscured by patterns observed at inter-clade scales [27]. Consequently, for example, although Ashton [45] observed that salamanders tend to follow Bergmann's rule, we could expect that a test incorporating a substantially larger number of species from a wider diversity of areas find that these amphibians do not follow Bergmann's clines (the situation might be different for turtles, as these animals are more restricted to tropical zones).
Conclusion
Physiological processes in ectotherms are strongly dependent upon their body temperature [3,80] with numerous consequences for survival and reproduction. Since heating rates in ectotherms are critically determined by body mass, large body size in species occurring at high latitudes or elevations is likely to be disadvantageous. We found no evidence to support Bergmann's rule in Liolaemus lizards: increasing latitude and elevation are not associated with larger body size in these reptiles (Fig. 2). We suggest that Bergmann's rule should be recognized as a macroecological prediction employed to investigate patterns of body size evolution only in endotherms, as it was originally proposed by Bergmann [21].
Methods
Phylogenetic scale and study species
Even though many studies testing Bergmann's rule have focused on intraspecific comparisons [19,34,37,73,98,99], the rule [6,21,23] was originally developed on the basis of multiple species comparisons [21] which is the approach that we took. Meiri & Thomas [100] provided a detailed discussion on the taxonomical scale of Bergmann's rule in an historical context.
We gathered data on body size with respect to latitude and elevation from 4942 adult Liolaemus specimens of both sexes representing a total of 126 species from museum collections (see appendix and Additional file 1: Supplementary table). Since body size of lizards may be affected by distribution in islands [101-103], we excluded insular Liolaemus taxa from our dataset (e.g. L. brattstroemi, L. melaniceps). Collection numbers and localities for most of the studied species and specimens can be found in Pincheira-Donoso & Núñez [48]. The species forming our sample cover a representative proportion of the total biogeographical, ecological, morphological and taxonomical diversity of the genus Liolaemus. The studied taxa belong to the clades chiliensis, Donosolaemus-magellanicus, fitzingerii, lineomaculatus, montanus and wiegmannii [11,48,57]. These groups represent the six main lineages known for this genus. The studied species encompass almost the entire geographical range known for Liolaemus [57,104]. The total dataset comprises individuals coming from austral Patagonia in Argentina and Chile, the high Andes plateau (> 4500–5000 m), the Atacama Desert, tropical areas in eastern Brazil, austral forests in southern Chile and several intermediate and temperate areas [48,50,51,54,57,60,104].
Environmental estimations
Since environmental temperatures decrease with increasing latitude and altitude [27,45,78], it is necessary to account for the combined effect of latitude and elevation when testing thermal dependence of traits [27]. Latitudinal and altitudinal data were used as estimators of the thermal conditions under which species live by converting them to a single Adjusted Latitudinal Midpoint (ALM) for each species. This scale, recently calibrated by Cruz et al. [27], and similar versions have been used for estimations of species' environmental conditions in comparative studies [11,39,105,106]. The ALM is obtained assuming that environmental temperature in altitudinal transects declines 0.65°C for each 100 m of increased elevation [27,78]. Cruz et al. [27] obtained a corrected latitudinal value for latitude and altitudinal thermal covariation using the formula y = 0.009x - 6.2627, where x represents the altitudinal midpoint for each species, and y the corrected temperature for latitude which is added to the latitudinal midpoint for each species (the details of this formula were provided personally by F.B. Cruz, as they are not published in the same format in Cruz et al. [27]). The final value is referred to as the adjusted latitudinal midpoint for South American areas where Liolaemus occurs, calculated for each species [27].
Statistical analysis and phylogenetic control
To investigate the effects of environmental variation on Liolaemus species size, we used snout-vent length (SVL) as a proxy for body size. Snout-vent length is a widely used measure of body size in squamate reptiles. This parameter is positively correlated with other body variables, such as mass, and with ecological and life-history traits [3,27,86]. Several recent studies have selected the largest recorded values for SVL as an estimation of species size [27,57,107]. However, utilization of single extreme values may lead to misleading results when testing these kinds of hypotheses [108,109]. Recent evidence [110] shows that estimation of asymptotic size can be reasonably obtained using the largest individual per sample, but only in lineages that follow asymptotic growth curves, such as Anolis lizards [110]. Since it is unknown whether Liolaemus follow asymptotic growth curves [39], the use of the largest recorded individual per sample may bias analyses. An alternative method that offers greater statistical power is the estimation of body size by obtaining intermediate percentiles between the maximum record and the sample mean. For example, it has been found in agamid species that the use of mean sample values may lead to substantial underestimation of asymptotic size, whereas the use of the largest individual per sample may lead to overestimation [108]. In contrast, these intermediate percentiles (between mean and maximum record) provide the most accurate estimate of asymptotic size. Moreover, such percentiles have low variance, low dependence on sample size, and are amenable to bootstrap estimations of confidence intervals, when compared to the largest individual per sample [108]. Consequently, we used the mean SVL of the largest two-thirds of the adults for each studied species [111]. Since species of the genus Liolaemus exhibit various patterns of sexual size dimorphism (i.e. larger males, larger females or no sexual dimorphism) [48,50,54,57,60,61,112,113], SVL mean values were calculated for males and females, to estimate a mean value for the species [27,39]. Whenever possible (117 of 126 species) samples comprised a similar number of males and females.
Prior to analyses SVL values were ln-transformed. This ln-transformation reduces skewness of the original measurements, and helps to homogenize variance [114-116]. After ln-transformation, SVL met the statistical assumptions required for parametric analyses (according to Shapiro-Wilk tests).
It is now well appreciated that samples consisting of species sharing phylogenetic histories cannot be considered independent evolutionary entities in comparative analyses [43,117]. Consequently, simple statistical analyses might provide phylogenetically biased evidence, necessitating the development of explicit approaches to control the impact of shared ancestry [43,44,118]. On the other hand, recent studies have shown that conventional regression and correlation analyses on raw data may perform better than independent contrasts data under certain conditions of phenotypic evolution [119-122]. For example, Carvalho et al. [121] suggested that it is not always necessary to perform phylogenetic simulations for statistical analyses when there is little phylogenetic effect, and that comparative analyses in some cases should be applied as a conservative approach. Hence, providing results from both conventional statistical procedures and explicit phylogenetic analyses may be a more robust approach [2,119]. In the case of body size patterns observed in Liolaemus species, the distribution of this trait across the phylogeny may reveal only a partial effect of shared ancestry. Therefore, we analyzed our data in two different ways. First, using conventional ordinary least squares bivariate regression analyses on raw ln-transformed data, separately for each of the main clades [43,44]. We identified these clades as detailed above following the latest available nomenclature [11,27,48,57] (Fig. 1). Each of the recognized clades has previously been supported by different phylogenetic hypotheses based on molecular, morphological and combined molecular and morphological evidence [11,27,51,57,59,123].
Second, we analyzed the dataset using an explicit phylogenetic approach, by calculating independent contrasts [43] as implemented in COMPARE version 4.6b [124]. We examined species relationships using a phylogenetic hypothesis of the genus Liolaemus derived from different recent studies [11,27,51,59,123]. Because this phylogeny (Fig. 1) is based on both molecular and morphological data, we performed phylogenetic analyses under a speciational Brownian motion model of evolutionary change, assuming branch lengths equal to 1.0 [11,125-127]. Since phylogenetic information was only available for 63 of the 126 species included in our dataset (see additional file 1: Supplementary table), we conducted all analyses in three different ways to test for the potential effect that reduction of sample size can have on the results: (i) ordinary least squares regressions on raw data separately for each of the six clades, including the entire species sample (i.e. 126 species), (ii) the same regression analyses (on raw data) for separate clades conducted in the entire dataset, but using only the 63 species included in the phylogeny, and (iii) ordinary least squares regression analyses (through the origin) on phylogenetically independent contrast for those species for which phylogenetic information was available. Since the existence of outliers may have large effects on statistical analyses [116], we performed all the regressions including and excluding outliers. Data points with standardized residuals (i.e. z-scores) greater than 2 or less than -2 were considered outliers [114]. Except for one analysis (non-phylogenetic regression on the montanus clade), regressions including and excluding outliers did not differ qualitatively.
Authors' contributions
DPD participated in the design of the study, conducted all the field work, performed the statistical and phylogenetic analyses and wrote the manuscript.
DJH participated in the design of the study, in the statistical analyses and in the writing of the manuscript.
TT participated in the design of the study, in the statistical analyses and in the writing of the manuscript.
Appendix
The studied material is housed in the herpetological collections of the following institutions. Collections identified with an asterisk (*) indicate the existence of specimens with collection data, but without an official collection number at the moment of our study. Division Reptiles and Amphibians, Museo Nacional de Historia Natural de Chile (MNHN*), Zoological Museum, Facultad de Ciencias Naturales y Oceanograficas, Universidad de Concepcion, Chile (MZUC*), Division Zoology, Museo de Historia Natural de Concepcion, Chile (MHNC*), Department of Cell Biology and Genetics, Facultad de Medicina, Universidad de Chile (DBCGUCH*), Instituto de Biologia Animal, Facultad de Ciencias Agrarias, Universidad Nacional de Cuyo, Argentina (IBAUNC*), Instituto Argentino de Investigaciones de las Zonas Aridas, CRICYT, Argentina (IADIZA), Natural History Museum of London (NHML), Jose Miguel Cei Diagnostic Collection (JMC-DC), Jose Alejandro Scolaro Diagnostic Collection (JAS-DC), and in the Herpetological Collection of the senior author, D. Pincheira-Donoso (CHDPD*).
Supplementary Material
Additional file 1
Supplementary Table. Snout-vent length (SVL), latitudinal and altitudinal range of the species included in this study. Taxa for which phylogenetic information is available are indicated in bold.
Acknowledgements
We thank H. Núñez (MNHN), J. Artigas, M. Contreras and E. Solar (MZUC), J. Navarro (DBCGUCH), E. Pereyra (IBAUNC), F. Videla (IADIZA), C. McCarthy (NHML), J.M. Cei (JMC-DC) and J.A. Scolaro (JASDC) for permission to study collections under their direction or for providing data from these samples. F.M. Jaksic, H. Núñez, J.A. Scolaro provided essential literature. Shai Meiri and an anonymous referee provided valuable observations to improve significantly the manuscript. DP-D thanks the financial support provided by Universities UK through an Overseas Research Student Award, the University of Exeter for an Exeter Research Student Award and a School of Biological and Chemical Sciences PhD Scholarship, Phrynosaura Chile and Oxford University Press. TT is funded by a Royal Society Research Fellowship. This work was supported by the European Social Fund.
References
  • Coyne, JA; Orr, HA. Speciation. Massachusetts , Sinauer Associates; 2004.
  • Schluter, D. The ecology of adaptive radiation. Oxford , Oxford University Press; 2000. p. 288.
  • Pough, FH; Andrews, RM; Cadle, JE; Crump, ML; Savitzky, AH; Wells, KD. Herpetology. Third . New Jersey , Pearson, Prentice Hall; 2004.
  • Gaston, KJ; Blackburn, TM. Pattern and process in macroecology. Massachusetts , Blackwell Science; 2000.
  • Huey, RB; Hertz, PE; Sinervo, B. Behavioral drive versus behavioural inertia in evolution: a null model approach. American Naturalist. 2003;161:357–366. doi: 10.1086/346135. [PubMed]
  • Blackburn, TM; Gaston, KJ; Loder, N. Geographic gradients in body size: a clarification of Bergmann’s rule. Diversity and Distributions. 1999;5:165–174. doi: 10.1046/j.1472-4642.1999.00046.x.
  • Bogert, CM. Thermoregulation in reptiles, a factor in evolution. Evolution. 1949;3:195–211. doi: 10.2307/2405558. [PubMed]
  • Brown, JH; Gillooly, JF; Allen, AP; Savage, VM; West, GB. Toward a metabolic theory of ecology. Ecology. 2004;85:1771–1789. doi: 10.1890/03-9000.
  • Brown, JH; Sibly, RM. Life-history evolution under a production constraint. Proceedings of the National Academy of Sciences, USA. 2006;103:17595–17599. doi: 10.1073/pnas.0608522103.
  • Andersson, M. Sexual selection. Princeton , Princeton University Press; 1994.
  • Espinoza, RE; Wiens, JJ; Tracy, CR. Recurrent evolution of herbivory in small, cold-climate lizards: Breaking the ecophysiological rules or reptilian herbivory. Proceedings of the National Academy of Sciences, USA. 2004;101:16819–16824. doi: 10.1073/pnas.0401226101.
  • Reznick, DN; Shaw, FH; Rodd, FH; Shaw, RG. Evaluation of the rate of evolution in natural populations of guppies (Poecilia reticulata). Science. 1997;275:1934–1937. doi: 10.1126/science.275.5308.1934. [PubMed]
  • Shine, R. Vertebral numbers in male and female snakes: the roles of natural, sexual and fecundity selection. Journal of Evolutionary Biology. 2000;13:455–465. doi: 10.1046/j.1420-9101.2000.00181.x.
  • Shine, R. Life-history evolution in reptiles. Annual Reviews of Ecology, Evolution and Systematics. 2005;36:23–46. doi: 10.1146/annurev.ecolsys.36.102003.152631.
  • Darwin, C. The descent of man and selection in relation to sex. London , John Murray; 1871.
  • Peters, RH. The ecological implications of body size. Cambridge , Cambridge University Press; 1983.
  • Angilletta, MJ; Dunham, AE. The temperature-size rule in ectotherms: simple evolutionary explanations may not be general. American Naturalist. 2003;162:332–342. doi: 10.1086/377187. [PubMed]
  • Partridge, L; Coyne, JA. Bergmann’s rule in ectotherms: Is it adaptive? Evolution. 1997;51:632–635. doi: 10.2307/2411139.
  • Walters, RJ; Hassall, M. The temperature-size rule in ectotherms: May a general explanation exist after all? American Naturalist. 2006;167:510–523. doi: 10.1086/501029. [PubMed]
  • Olalla-Tarraga, MA; Rodriguez, MA; Hawkins, BA. Broad-scale patterns of body size in squamate reptiles of Europe and North America. Journal of Biogeography. 2006;33:781–793. doi: 10.1111/j.1365-2699.2006.01435.x.
  • Bergmann, C. Ueber die Verhaltnisse der warmeokonomie der thiere zu ihrer grosse. Gottinger Studien. 1847;3:595–708.
  • Blackburn, TM; Ruggiero, A. Latitude, elevation and body mass variation in Andean passerine birds. Global Ecology and Biogeography. 2001;10:245–259. doi: 10.1046/j.1466-822X.2001.00237.x.
  • James, FC. Geographic size variations in birds and its relationship with climate. Ecology. 1970;51:365–390. doi: 10.2307/1935374.
  • Arnett, AE; Gotelli, NJ. Geographic variation in life-history traits of the ant lion, Myrmeleon immaculatus: evolutionary implications of Bergmann’s rule. Evolution. 1999;53:1180–1188. doi: 10.2307/2640821.
  • Kaspari, M; Vargo, EL. Colony size as a buffer against seasonality: Bergmann’s rule in social insects. American Naturalist. 1995;145:610–632. doi: 10.1086/285758.
  • Mousseau, TA. Ectotherms follow the converse to Bergmann’s rule. Evolution. 1997;51:630–632. doi: 10.2307/2411138.
  • Cruz, FB; Fitzgerald, LA; Espinoza, RE; Schulte, JA. The importance of phylogenetic scale in tests of Bergmann’s and Rapoport’s rules: lessons from a clade of South American lizards. Journal of Evolutionary Biology. 2005;18:1559–1574. [PubMed]
  • Atkinson, D. Temperature and organism size: a biological law for ectotherms. Advances in Ecological Research. 1994;25:1–58.
  • Belk, MC; Houston, DD. Bergmann’s rule in ectotherms: A test using freshwater fishes. American Naturalist. 2002;160:803–808. doi: 10.1086/343880.
  • de Queiroz, A; Ashton, KG. The phylogeny of a species-level tendency: Species heritability and possible deep origins of Bergmann's rule in tetrapods. Evolution. 2004;58:1674–1684. [PubMed]
  • Ashton, KG; Feldman, CR. Bergmann’s rule in nonavian reptiles: turtles follow it, lizards and snakes reverse it. Evolution. 2003;57:1151–1163. [PubMed]
  • Cushman, JH; Lawton, JH; Manly, BFJ. Latitudinal patterns in Europe ant assemblages: variation in species richness and body size. Oecologia. 1993;95:30–37.
  • Lindstet, SL; Boyce, MS. Seasonality, fasting endurance, and body size in mammals. American Naturalist. 1985;125:873–878. doi: 10.1086/284385.
  • Van Voorhies, WA. Bergmann size clines: A simple explanation for their occurrence in ectotherms. Evolution. 1996;50:1259–1264. doi: 10.2307/2410666.
  • Van Voorhies, WA. On the adaptive nature of Bergmann size cline: a reply to Mousseau, Partridge and Coyne. Evolution. 1997;51:635–640. doi: 10.2307/2411140.
  • Ashton, KG; Tracy, MC; de Queiroz, A. Is Bergmann’s rule valid for mammals? American Naturalist. 2000;156:390–415. doi: 10.1086/303400.
  • Graves, GR. Bergmann’s rule near the equator: latitudinal clines in body size of an Andean passerine bird. Proceedings of the National Academy of Sciences, USA. 1991;88:2322–2325. doi: 10.1073/pnas.88.6.2322.
  • Meiri, S; Dayan, T. On the validity of Bergmann’s rule. Journal of Biogeography. 2003;30:331–351.
  • Pincheira-Donoso, D; Tregenza, T; Hodgson, DJ. Body size evolution in South American Liolaemus lizards of the boulengeri clade: a contrasting reassessment. Journal of Evolutionary Biology. 2007;20:2067–2071. doi: 10.1111/j.1420-9101.2007.01394.x. [PubMed]
  • Atkinson, D; Sibly, RM. Why are organisms usually bigger in colder environments? Making sense of a life history puzzle. Trends in Ecology and Evolution. 1997;12:235–239. doi: 10.1016/S0169-5347(97)01058-6.
  • Lindsey, CC. Body sizes of poikilotherm vertebrates at different latitudes. Evolution. 1966;20:456–465. doi: 10.2307/2406584.
  • Ray, C. The application of Bergmann’s and Allen’s rules to the poikilotherms. Journal of Morphology. 1960;106:85–108. doi: 10.1002/jmor.1051060104. [PubMed]
  • Felsenstein, J. Phylogenies and the comparative method. American Naturalist. 1985;125:1–15. doi: 10.1086/284325.
  • Harvey, PH; Pagel, MD. The comparative method in evolutionary biology. Oxford , Oxford University Press; 1991.
  • Ashton, KG. Do amphibians follow Bergmann’s rule? Canadian Journal of Zoology. 2002;80:708–716. doi: 10.1139/z02-049.
  • Ashton, KG. Body size variation among mainland populations of the western rattlesnake (Crotalus viridis). Evolution. 2001;55:2523–2533. [PubMed]
  • Ashton, KG. Are ecological and evolutionary rules being dismissed prematurely? Diversity and Distributions. 2001;7:289–295. doi: 10.1046/j.1366-9516.2001.00115.x.
  • Pincheira-Donoso, D; Núñez, H. Las especies chilenas del género Liolaemus. Taxonomía, sistemática y evolución. Santiago , National Museum of Natural History Press; 2005. pp. 1–487.
  • Etheridge, R; Espinoza, RE. Taxonomy of the Liolaeminae (Squamata: Iguania: Tropiduridae) and a semi-annotated bibliography. Smithsonian Herpetological Information Service. 2000;126:1–64.
  • Pincheira-Donoso, D. Anfibios y reptiles de la Provincia de El Loa. In: Ramírez GM, Pincheira-Donoso D. , editor. Fauna del Altiplano y Desierto de Atacama Vertebrados de la Provincia de El Loa. Calama , Phrynosaura Ediciones; 2005. pp. 93–150.
  • Schulte, JA; Losos, JB; Cruz, FB; Núñez, H. The relationship between morphology, escape behaviour and microhabitat occupation in the lizard clade Liolaemus (Iguanidae: Tropidurinae: Liolaemini). Journal of Evolutionary Biology. 2004;17:408–420. doi: 10.1046/j.1420-9101.2003.00659.x. [PubMed]
  • Sura, P. Encyclopedia of extant amphibians and reptiles. Nowy Sacz , Wydawnictwo Fundacja; 2005. p. 544.
  • Pincheira-Donoso, D; Scolaro, JA; Sura, P. A monographic catalogue on the systematics and phylogeny of the South American iguanian lizard family Liolaemidae (Squamata, Iguania). Zootaxa. 2008.
  • Cei, JM. Reptiles del centro, centro-oeste y sur de la Argentina. Herpetofauna de las zonas áridas y semiáridas. Torino , Museo Regionale di Scienze Naturali di Torino; 1986. p. 527.
  • Donoso-Barros, R. Liolaemus. In: Peters JA, Donoso-Barros R. , editor. Catalogue of the Neotropical Squamata Part II Lizards and Amphisbaenians. Bulletin United States National Museum 297; 1970. pp. 170–195.
  • Jaksic, FM; Schwenk, K. Natural history observations on Liolaemus magellanicus, the southernmost lizard in the world. Herpetologica. 1983;39:457–461.
  • Etheridge, R. A review of lizards of the Liolaemus wiegmannii group (Squamata, Iguania, Tropiduridae), and a history of morphological change in the sand-dwelling species. Herpetological Monographs. 2000;14:293–352. doi: 10.2307/1467049.
  • Laurent, RF. New forms of lizards of the subgenus Eulaemus of the genus Liolaemus (Reptilia: Squamata: Tropiduridae) from Perú and northern Chile. Acta Zoológica Lilloana. 1998;44:1–26.
  • Pincheira-Donoso, D; Scolaro, JA; Schulte II, JA. The limits of polymorphism in Liolaemus rothi: molecular and phenotypic evidence for a new species of the Liolaemus boulengeri clade (Iguanidae, Liolaemini) from boreal Patagonia of Chile. Zootaxa. 2007;1452:25–42.
  • Cei, JM. Reptiles del noroeste, nordeste y este de la Argentina. Herpetofauna de las selvas subtropicales, puna y pampas. Torino , Museo Regionale di Scienze Naturali di Torino; 1993. p. 947.
  • Donoso-Barros, R. Reptiles de Chile. Santiago , Ediciones Universidad de Chile; 1966.
  • Calder, WA. Size, function and life history. Massachusetts , Harvard University Press; 1984.
  • Fitch, HS. Autecology of the copperhead. University of Kansas Museum of Natural History, Miscellaneous Publications. 1960;13:85–288.
  • Hirth, HF; King, AC. Body temperatures of snakes in different seasons. Journal of Herpetology. 1969;3:101–102. doi: 10.2307/1563230.
  • Kurfess, JF. Mating, gestation and growth rate in Lichanura r. roseofusca. Copeia. 1967;1967:477–479. doi: 10.2307/1442149.
  • Shine, R. "Costs" of reproduction in reptiles. Oecologia. 1980;46:92–100. doi: 10.1007/BF00346972.
  • Bartholomew, GA. Physiological control of body temperature. In: Gans C, Pough FH. , editor. Biology of the Reptilia Vol 12. New York , Academic Press; 1982. pp. 167–213.
  • Díaz, JA; Bauwens, D; Asensio, B. A comparative study of the relation between heating rates and ambient temperatures in lacertid lizards. Physiological Zoology. 1996;69:1359–1383.
  • Heath, JE. Temperature regulation and diurnal activity in horned lizards. University of California Publications in Zoology. 1965;64:97–136.
  • Muth, A. Thermoregulatory postures and orientation to the sun: a mechanistic evaluation for the zebra-tailed lizard, Callisaurus draconoides. Copeia. 1977;1977:710–720. doi: 10.2307/1443171.
  • Hellmich, W. On ecotypic and autotypic characters, a contribution to the knowledge of the evolution of the genus Liolaemus (Iguanidae). Evolution. 1951;5:359–369. doi: 10.2307/2405681.
  • Jaksic, FM. Ecología de los vertebrados de Chile. Santiago , Ediciones de la Universidad Católica de Chile; 1998.
  • Sears, MW; Angilletta, MJ. Body size clines in Sceloporus lizards: proximate mechanisms and demographic constraints. Integrative and Comparative Biology. 2004;44:433–442. doi: 10.1093/icb/44.6.433.
  • Fox, SF; Shipman, PA. Social behavior at high and low elevations. Environmental release and phylogenetic effects in Liolaemus. In: Fox SF, McCoy JK, Baird TA. , editor. Lizard social behavior. Baltimore & London , Johns Hopkins University Press; 2003. pp. 310–355.
  • Fuentes, ER; Jaksic, FM. Activity temperatures of eight Liolaemus (Iguanidae) species in central Chile. Copeia. 1979;1979:546–548. doi: 10.2307/1443242.
  • Avila-Pires, TCS. Lizards of Brazilian Amazonia (Reptilia: Squamata). Zoologische Verhandelingen. 1995;299:1–706.
  • Conrad, V; Pollak, LW. Methods in climatology. Massachusetts , Harvard University Press; 1950.
  • Lutgens, FK; Tarbuck, EJ. The atmosphere, and introduction to meteorology. New Jersey , Prentice Hall; 1998.
  • Spellerberg, IF. Adaptations of reptiles to cold. In: Bellairs A, Barry-Cox C. , editor. Morphology and biology of reptiles. Vol. 3. Linnean Society Symposium Series; 1976. pp. 261–285.
  • Pianka, ER; Vitt, LJ. Lizards. A windows to the evolution of diversity. Berkeley, Los Angeles & London , University of California Press; 2003. p. 333.
  • Hailey, A; Davies, PMC. Activity and thermoregulation of the snake Natrix maura. 2. A synoptic model of thermal biology and the physiological ecology of performance. Journal of Zoology. 1988;214:325–342.
  • Ibargüengoytía, NR. Field, selected body temperature and thermal tolerance of the syntopic lizards Phymaturus patagonicus and Liolaemus elongatus (Iguania: Liolaemidae). Journal of Arid Environments. 2005;62:435–448. doi: 10.1016/j.jaridenv.2005.01.008.
  • Martin-Vallejo, J; Garcia-Fernandez, J; Perez-Mellado, V; Vicente-Vellardon, JL. Habitat selection and thermal ecology of the sympatric lizards Podarcis muralis and Podarcis hispanica in a mountain region of central Spain. Herpetological Journal. 1995;5:181–188.
  • Brattstrom, BH. Body temperatures of reptiles. American Midland Naturalist. 1965;73:376–422. doi: 10.2307/2423461.
  • Cooper, WE; Vitt, LJ. Distribution, extent, and evolution of plant consumption by lizards. Journal of Zoology. 2002;257:487–517. doi: 10.1017/S0952836902001085.
  • Pough, FH. Lizard energetics and diet. Ecology. 1973;54:837–844. doi: 10.2307/1935678.
  • Olesen, JM; Valido, A. Lizards as pollinators and seed dispersers: an island phenomenon. Trends in Ecology and Evolution. 2003;18:177–181. doi: 10.1016/S0169-5347(03)00004-1.
  • Alexander, RM. Energy for animal life. Oxford , Oxford University Press; 1999.
  • King, G. Reptiles and herbivory. New York , Chapman and Hall; 1996.
  • Wilson, KJ; Lee, AK. Energy expenditure of a large herbivorous lizard. Copeia. 1974;1974:338–348. doi: 10.2307/1442528.
  • Stevens, CE; Hume, ID. Comparative physiology of the vertebrate digestive system. Cambridge , Cambridge University Press; 1995.
  • Troyer, K. Transfer of fermentative microbes between generations in a herbivorous lizard. Science. 1982;216:540–542. doi: 10.1126/science.216.4545.540. [PubMed]
  • Schall, JJ; Dearing, MD. Body temperature of the herbivorous bonaire island whiptail lizard (Cnemidophorus murinus). Journal of Herpetology. 1994;28:526–528. doi: 10.2307/1564976.
  • Vitt, LJ. Shifting paradigms: herbivory and body size in lizards. Proceedings of the National Academy of Sciences, USA. 2004;101:16713–16714. doi: 10.1073/pnas.0407439101.
  • Losos, JB; Schoener, TW; Langerhans, RB; Spiller, DA. Rapid temporal reversal in predator-driven natural selection. Science. 2006;314:1111. doi: 10.1126/science.1133584. [PubMed]
  • Losos, JB; Schoener, TW; Spiller, DA. Predator-induced behaviour shifts and natural selection in field-experimental lizard populations. Nature. 2004;432:505–508. doi: 10.1038/nature03039. [PubMed]
  • Brown, JH. Macroecology. Chicago , University of Chicago Press; 1995.
  • Mayr, E. Geographical character gradients and climatic adaptation. Evolution. 1956;10:105–108. doi: 10.2307/2406103.
  • Rensch, B. Some problems of geographical variation and species-formation. Proceedings of the Linnean Society of London. 1938;150:275–285.
  • Meiri, S; Thomas, GH. The geography of body size - challenges of the interspecific approach. Global Ecology and Biogeography. 2007;16:689–693. doi: 10.1111/j.1466-8238.2007.00343.x.
  • Case, TJ. A general explanation for insular body size trends in terrestrial vertebrates. Ecology. 1978;59:1–18. doi: 10.2307/1936628.
  • Dunham, AE; Tinkle, DW; Gibbons, JW. Body size in island lizards: a cautionary tale. Ecology. 1978;59:1230–1238. doi: 10.2307/1938236.
  • Stone, PA; Snell, HL; Snell, HM. Island biogeography of morphology and social behaviour in the lava lizards of the Galapagos islands. In: Fox SF, McCoy JK, Baird TA. , editor. Lizard social behavior. Baltimore and London , Johns Hopkins University Press; 2003. pp. 190–239.
  • Donoso-Barros, R. Catálogo herpetológico chileno. Boletín del Museo Nacional de Historia Natural de Chile. 1970;31:49–124.
  • Wiens, JJ; Parra-Olea, G; García-París, M; Wake, DB. Phylogenetic history underlies elevational biodiversity patterns in tropical salamanders. Proceedings of the Royal Society of London, Biological Sciences. 2007;274:919–928. doi: 10.1098/rspb.2006.0301.
  • Pincheira-Donoso, D; Hodgson, DJ; Tregenza, T. Comparative evidence for strong phylogenetic inertia in precloacal signalling glands in a species-rich lizard clade. Evolutionary Ecology Research. 2008;10:11–28.
  • Etheridge, R; Christie, MI. Two new species of the lizard genus Liolaemus (Squamata: Liolaemidae) from northern Patagonia, with comments on Liolaemus rothi. Journal of Herpetology. 2003;37:325–341. doi: 10.1670/0022-1511(2003)037[0325:TNSOTL]2.0.CO;2.
  • Brown, RP; Znari, M; El Mouden, ELH; Harris, P. Estimating asymptotic body size and testing geographic variation in Agama impalearis. Ecography. 1999;22:277–283. doi: 10.1111/j.1600-0587.1999.tb00503.x.
  • Meiri, S. Size evolution in island lizards. Global Ecology and Biogeography. 2007;16:702–708. doi: 10.1111/j.1466-8238.2007.00327.x.
  • Stamps, JA; Andrews, RM. Estimating asymptotic size using the largest individuals per sample. Oecologia. 1992;92:203–512. doi: 10.1007/BF00317842.
  • Losos, JB; Butler, M; Schoener, TW. Sexual dimorphism in body size and shape in relation to habitat use among species of Caribbean Anolis lizards. In: Fox SF, McCoy JK, Baird TA. , editor. Lizard social behaviour. Baltimore and London , John Hopkins University Press; 2003. pp. 356–380.
  • Etheridge, R. Lizards of the Liolaemus darwinii Complex (Squamata: Iguania: Tropiduridae) in Northern Argentina. Bollettino del Museo Regionale di Scienze Naturali di Torino. 1993;11:137–199.
  • Scolaro, JA; Cei, JM. Systematic status and relationships of Liolaemus species of the archeforus and kingii groups: a morphological and taxonumerical approach (Reptilia: Tropiduridae). Bollettino del Museo Regionale di Scienze Naturali di Torino. 1997;15:369–406.
  • Field, A. Discovering statistics using SPSS. Second Edition. London , Sage; 2006. p. 779.
  • Miles, DB; Ricklefs, RE. The correlation between ecology and morphology in deciduous forest passerine birds. Ecology. 1984;65:1629–1640. doi: 10.2307/1939141.
  • Zar, JH. Biostatistical analysis. New Jersey , Prentice-Hall; 1999.
  • Cheverud, JM; Dow, MM; Leutenegger, W. The quantitative assessment of phylogenetic constraints in comparative analyses: sexual dimorphism in body weight among primates. Evolution. 1985;39:1335–1351. doi: 10.2307/2408790.
  • Brooks, DR; McLennan, DA. Phylogeny, ecology and behaviour. A research programme in comparative biology. Chicago , University of Chicago Press; 1991.
  • Harvey, PH; Rambaut, A. Comparative analyses for adaptive radiations. Philosophical Transactions of the Royal Society of London, B. 2000;355:1599–1605. doi: 10.1098/rstb.2000.0715.
  • Björklund, M. Are "comparative methods" always necessary? Oikos. 1997;80:607–612. doi: 10.2307/3546637.
  • Carvalho, P; Diniz-Filho, JAF; Bini, LM. Factors influencing changes in trait correlations across species after using phylogenetic independent contrasts. Evolutionary Ecology. 2006;20:591–602. doi: 10.1007/s10682-006-9119-7.
  • Ricklefs, RE; Starck, JM. Applications of phylogenetically independent contrasts: a mixed progress report. Oikos. 1996;77:167–172. doi: 10.2307/3545598.
  • Schulte, JA; Macey, JR; Espinoza, RE; Larson, A. Phylogenetic relationships in the iguanid lizard genus Liolaemus: multiple origins of viviparous reproduction and evidence for recurring Andean vicariance and dispersal. Biological Journal of the Linnean Society. 2000;69:75–102. doi: 10.1006/bijl.1999.0346.
  • Martins, EP. COMPARE, version 4.6b. Computer programs for the statistical analysis of comparative data. Distributed by the author at http://compare.bio.indiana.edu/. Indiana , Department of Biology, Indiana University; 2004.
  • Garland, T; Dickerman, AW; Janis, CM; Jones, JA. Phylogenetic analysis of covariance by computer simulation. Systematic Biology. 1993;42:265–292. doi: 10.2307/2992464.
  • Martins, EP; Garland, T. Phylogenetic analyses of the correlated evolution of continuous characters: a simulation study. Evolution. 1991;45:534–557. doi: 10.2307/2409910.
  • Rohlf, FJ; Chang, WS; Sokal, RR; Kim, J. Accuracy of estimated phylogenies: effects of tree topology and evolutionary model. Evolution. 1990;44:1671–1684. doi: 10.2307/2409346.