# Climate-Related Variation in Body Dimensions within Four Lacertid Species

Authors: Stanislav Volynchik
Journal: International Journal of Zoology
 (2014)

---

## Abstract

A close relationship between habitat and external morphology is widespread among many animals, including reptiles. Here, I studied the relationship between abiotic environmental conditions and body size of four lacertid species (Phoenicolacerta laevis, Ophisops elegans, Acanthodactylus boskianus, and Mesalina guttulata) occurring in Israel. I examined the effect of average annual temperature and average annual precipitation on body and limb dimensions, using linear statistical models. Temperature- and precipitation-related geographic clines in body size showed the same trend among all species. Females displayed stronger phenotypic response to temperature gradient than conspecific males, suggesting a sex-specific effect of natural selection. Snout-vent length (SVL) was negatively correlated with temperature, supporting Bergmann’s rule in O. elegans and in female P. laevis and A. boskianus, but not in M. guttulata. Precipitation was positively related to SVL in O. elegans and M. guttulata, and in female P. laevis and A. boskianus. The relative extremity lengths, especially hind limb segments, generally increase towards hot and dry locations, following Allen’s rule. Among the Mediterranean region species (P. laevis, O. elegans) the morphological-environmental link with temperature was stronger than in desert dwellers (A. boskianus, M. guttulata), for which precipitation was the major determinant of spatial variation.


---

## Body

## 1. Introduction

Latitudinal variation in corporeal dimensions has often been explained as a response to temperature gradients, which are strong predictors of animals’ external morphology. Zoogeographic principles such as Bergmann’s and Allen’s rules, in classical (thermal) interpretation, express the relationship between body size or shape and environmental temperature [1]. Therefore, the most renowned and generally accepted explanation for these rules is based on the heat-conservation hypothesis [2–4]: reduced surface area-to-volume ratio provides better heat retention in order to reduce body heat loss at low temperatures (Bergmann’s rule). Bergmann [5] stated that in cooler climates, large-bodied endotherms, with a low body surface area to body volume ratio (i.e., a greater thermal inertia), typically retain body heat better than smaller ones. Similar principles contributing to reduce heat loss underlie Allen’s rule [6], which predicted that animals living in colder environments tend to have shorter protruding body parts, such as ears, tails, and limbs. Furthermore, several authors [7, 8] have shown that the geographical intraspecific variation in corporeal dimensions of homeotherms has been related to either ambient temperature or moisture, or both. On the whole, there is evidence showing the link between phenotypic plasticity and food availability, that is, the “resource rule” [9–11], depending, in turn, on a complex of climatic (mainly water- and temperature-related) parameters. Yom-Tov and Nix [12] concluded that in some Australian mammals precipitation, which mostly contributes to primary production in hot and dry habitats, is often better correlated with body size than is temperature. Therefore, food resources constitute the major factor determining body size variation. Recently, Pincheira-Donoso and Meiri [13] based on an intercontinental dataset (65 reptile species from Israel, Argentina, and Chile) found that for both Israeli and South American squamates the resource rule is better predictor of body size than temperature.Multiple studies have noted that many tetrapods show size-latitude gradients. Clines in body dimensions have been well documented among many mammals, reptiles, and amphibians [14–16]. In reptiles, even earlier investigations [9] based on species assemblages have revealed that, in general, they do not exhibit noticeable latitudinal trends, although some snake taxa show a weak gradient. However, the appearance of spatial trends in body size for ectotherms [2], including vertebrates [17–19], is still disputed [13, 20–23]. In recent years, many studies have tested Bergmann’s rule for various ectotherm taxa and have found a wide range of patterns depicting the connection between body size and environmental gradient [16, 18, 24–27]. Atkinson and Sibly [28] and Cruz et al. [29] have claimed that some ectotherms follow Bergmann’s rule, whereas for other species the opposite pattern is true [30, 31]. Ashton and Feldman [17] found that turtles follow Bergmann’s rule whereas squamates contrast it. Cruz et al. [29] demonstrated a strong negative temperature-body size gradient amongLiolaemus lizard species, reflecting a profound effect of thermal inertia at higher latitudes or altitudes, in accordance with Bergmann’s rule. However, Pincheira-Donoso et al. [20, 32] within the same lineage, but based on an enlarged sample size, failed to find statistical support for Bergmann’s clines and, claiming that in these animals the evolution of larger bodies in colder areas is disadvantageous, they suggested that it applied exclusively to endotherms. It is possible that the spatial trends in body dimensions are influenced by a complex set of both abiotic (e.g., seasonality and temperature) and biotic (e.g., intra- and interspecific competition) factors and their interactions [17, 18, 33] and that no single mechanism can explain them.However, to my knowledge, some important issues remain obscure. For instance, tests of intraspecific body size or shape patterns in squamates accounting for gender-specific aspects are scarce because the papers often provide data for the pooled male and female samples. Taking into account that reptiles are often sexually dimorphic ([34, 35] see also the results) it is possible that conclusions of some ecogeographic studies may suffer from biased sampling. In this context, my primary goal is to compare the sexes in respect to their relationship with climate variables. A nonphylogenetic approach was used here to test the main hypothesis that males and females facing different thermal and metabolic demands follow different patterns in their phenotypic responses to environmental gradient. Further, I evaluate the combined effect of temperature and precipitation on body dimensions within all the species studied. Finally, I also focus on species-specific peculiarities in the appearance of environmentally related phenotypic variation. Diurnal lizards including arid-zone dwellers bask in large amounts of time, especially in the morning hours (personal observations). Under this scenario, I hypothesize that in the desertM. guttulata andA. boskianus, which are affected by exposure to extreme daily temperature fluctuations, a smaller size is likely to be preferable. Small-bodied individuals may have better heating and cooling abilities, heating faster after the cold nights, and can thus remain active for longer periods. It is generally accepted [36, 37] that temperature is the most important factor influencing lizards’ activities and that thermal demands play an important role in habitat selection by reptiles. In the wild, lizards apply a variety of strategies to achieve their preferred body temperature. Thermoregulation is realized behaviourally through shuttling between warm and cold sites [38], as well as physiologically by altering peripheral blood flow and modifying heart rates [39, 40]. In addition, some species apply ventilatory heat loss [41, 42] via evaporation through their mouth or cloaca.In arid areas, where the present study mainly took place, rainfall, as well as temperature, is an important factor in animal ecology, including determining food availability and thus also growth rate [11]. The dispute over which geographic component determines body size most makes an examination of the respective influence of temperature- and water-related variables on morphological features timely. In the present study, I used a range of linear measurements and ratios in four oviparous lizard species in order to determine the correlation between morphology, temperature, and precipitation. I examined environmentally induced phenotypic variability and the validity of ecomorphological patterns, such as Bergmann’s and Allen’s rules, in four Middle-Eastern wide-ranging lacertid lizards:Phoenicolacerta laevis,Ophisops elegans,Acanthodactylus boskianus, andMesalina guttulata.Acanthodactylus andMesalinaspecies are mainly desert inhabitants, while thePhoenicolacerta andOphisops species mostly inhabit the Mediterranean region (a more mesic area). In Israel these species occur along pronounced temperature and rainfall gradients ([43] and see below), which makes this region a convenient model for examining intraspecific ecomorphological variation. Thus, in Israel mean annual rainfall ranges from ~1000 mm in the north to 25 mm in the south. There is a negative north-south gradient in ambient temperature with annual averages of about 15°C and 25°C in the north and south, respectively. Another strong climatic feature is a U-shaped west/east temperature gradient from the temperate Mediterranean coast to the colder central mountain range, and the very hot Jordan Valley in the east (Figure 1). The daily and annual temperature fluctuations also differ greatly across Israel.Climatic maps of the studied areas (excluding Turkish localities). (a) Average annual temperature and (b) average annual precipitation. The ranges of temperature and precipitation for each species are given in brackets.
(a)
(b)Israel, which has a wide range of climates and landscapes, including mountain regions, coastal plains, humid valleys, rocky deserts, and sand dunes, is a hotspot of reptile species richness [44]. Environmental temperatures (e.g., extreme annual/diurnal values) and rainfall intensity (Figures 1(a) and 1(b)) are markedly different across the distribution ranges of both desert and the Mediterranean region species. Hence, the appearance of morphological variation in ectotherms is expected to be adjusted by species-specific ecological interactions leading to interspecific differences in body size (or shape) gradient. In this context, I expect a different link between both abiotic predictors and body/limb proportions within ecologically different lizard species. For example, temperature may be a stronger ecological determinant of body or extremities size variation forP. laevis andO. elegans living in cooler environments than for the desertA. boskianus andM. guttulata. The present study is the first including analyses of both Bergmann’s and Allen’s rules to investigate a functional connection between abiotic habitat conditions and the external morphology of the Middle East lacertid lizards; it also discusses the functional link between corporeal measurements and geographic variables. The findings are expected to contribute to our understanding of the potential importance of abiotic components as sources of morphological intraspecific variability in ectothermic vertebrates, whose metabolic rates and life history are greatly affected by environments.


## 2. Material and Methods

### 2.1. Samples and Site Locations

The observed specimens were adult males and females of four lizard species (more than 100 individuals per species), with large geographic ranges. I measured 679 preserved specimens of four lacertid species in the Natural History Museum of the Department of Zoology, Faculty of Life Sciences, Tel Aviv University (TAUM). The specimens had been collected between 1941 and 2012, mostly 86.6% from Israel but some 71 from Egypt (the Sinai Peninsula) and 20 specimens from Turkey.The following specimens were collected:Phoenicolacerta laevis, 104 specimens were collected from Mt. Hermon, the Golan Heights, the Galilee, the foothills of Judea, and the coastal plain (Israel);Ophisops elegans, 127 specimens were collected from Israel (Mt. Hermon, the Golan Heights, the Galilee, the Judean foothills, and the northern and central Negev) and 20 specimens from Turkey;Acanthodactylus boskianus, 164 specimens were collected from Israel (the A’rava Valley and the northern, central, and southern Negev) and 50 specimens from Egypt (northern Sinai, Sinai mountains, and the western Sinai coastal plain);Mesalina guttulata, 193 specimens were collected from Israel (the Judean Hills, the Judean Desert, Dead Sea area, and the northern, central, and southern Negev) and 21 specimens from Egypt (the northern Sinai and Sinai mountains).


### 2.2. Study Species

Phoenicolacerta laevis (Gray, 1838) is a medium-sized lizard, SVL up to 77 mm (according to our records), generally distributed in the Mediterranean landscape from Turkey to Israel and Jordan. This fairly adaptable species is found in relatively humid biotopes, such as temperate forests, valleys, Mediterranean-type shrub vegetation, cultivated fields, and gardens ([45, 46] pers. obs.). It is found from 200 m below sea level up to 2200 m above sea level (TAUM records).Ophisops elegans (Ménétriés, 1832) is a small, slender lizard. The largest specimen in our sample has a SVL of 59 mm. It occurs in open arid plains of the Mediterranean region and Central Asia and was found in sparse forests, vegetated rocky hillsides, and wadis [46–48]. Its vertical distribution ranges from 240 m below sea level to 2000 m above sea level (TAUM records).Acanthodactylus boskianus (Daudin, 1802) is a medium-sized lizard, with adults up to 88 mm SVL. It is widespread in desert and semidesert areas of the Saharo-Arabian region from northern Africa through the eastern Mediterranean and the Arabian Peninsula to Iraq. It inhabits rather hard substrate and occupies a wide range of sandy and gravel habitats. It was collected from 400 m below sea level to 1300 m above sea level [46, 48, 49].Mesalina guttulata (Lichtenstein, 1823) is a rather small, slender lizard, with a maximum SVL of 57 mm. It is widely distributed in North Africa and south-west Asia. It usually inhabits desert and semidesert regions with various terrains, prefers hard substrate and scattered rocks, and avoids sand dunes [46–48]. It is found from 350 m below sea level to 2000 m above sea level (TAUM records).


### 2.3. Morphological Measurements

For each specimen, I measured the following: snout-vent length (SVL); humerus length (HmL), distance from axilla to apex of elbow; forearm length (FaL), distance from elbow apex to apex of wrist; femur length (FeL), distance from groin to apex of knee; and tibia length (TbL), distance from apex of knee to heel. From these, total fore limb length (FLL) and total hind limb length (HLL), which were obtained as the sums of the corresponding segment lengths (i.e., HmL + FaL and FeL + TbL), as well as the limbs-body ratios (HmL/SVL, FaL/SVL, FeL/SVL, TbL/SVL) were calculated. Tail length was not included in the analysis because many individuals had broken or regenerated tails. For each specimen I also recorded collection data (region and locality with geographic coordinates and elevation), as well as body weight (M), which was taken shortly after capture to the nearest 0.1 g. Sex was determined by tail base shape and hemipenis eversion or, in doubtful cases, by dissection. Measurements were made with a digital caliper (Mitutoyo CD-6′′CX) to the nearest 0.01 mm, except for SVL, which was measured to the nearest 0.5 mm using a millimeter ruler. I could not record all measurements from each specimen due to the bad preservation conditions of some reptiles, so the sample sizes vary among features (Table 1). All measurements were made by the author.Table 1
Morphological characteristics of lizard species studied and their intersexual comparisons (one-way ANOVA with multiple testing corrections; Duncan’s multiple range test, MRT).
Species
Features
Range
Mean ± SD
n
♂
♂
/
♀
♀
P
F
♂
♂
♀
♀
♂
♂
♀
♀
Phoenicolacerta laevis
M (g)SVL (mm)HmL (mm)FaL (mm)FeL (mm)TbL (mm)
2.5–10.550.5–775.3–95.9–9.327.62–11.848.75–13
3–11.450–765.22–7.345.9–7.87.7–10.58.5–11.4
5.85
±
1.91 63.53±6.65 6.99±0.797.61±0.779.66±1.0311.04±1.02
5.13
±
1.94 62.97±6.016.34±0.476.96±0.468.78±0.639.87±0.62
48/5648/5647/5647/5647/5647/56
0.0590.6540.00010.00010.00010.0001
3.620.2026.3628.4927.4051.01
Ophisops elegans
M (g)SVL (mm)HmL (mm)FaL (mm)FeL (mm)TbL (mm)
0.9–4.240–594.45–6.14.78–6.725.3–8.727.15–10.53
1.2–3.540–53.53.9–5.44.23–6.225.2-7.476.46–9.22
2.5
±
0.72 48.07±4.275.3±0.45.82±0.57.22±0.719±0.81
1.98
±
0.51 46.47±3.444.66±0.355.15±0.416.29±0.577.78±0.63
85/6285/6284/6184/6184/6184/61
<0.00010.016<0.0001<0.0001<0.0001<0.0001
23.225.8598.3873.2069.1195.38
Acanthodactylus boskianus
M (g)SVL (mm)HmL (mm)FaL (mm)FeL (mm)TbL (mm)
1.7–13.852–885.2–9.185.92–9.928.45–13.3810.8–16.4
3–12.150–80.54.62–7.35.8–8.488.3–12.2810–14.54
6.99
±
2.15 65.9±6.466.9±0.757.96±0.7911.2±1.0413.77±1.19
5.35
±
1.63 61.11±5.836.26±0.567.12±0.629.81±0.7512.12±0.87
121/93121/93113/85113/85114/85114/85
<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001
37.2331.3143.3365.57108.62116.06
Mesalina guttulata
M (g)SVL (mm)HmL (mm)FaL (mm)FeL (mm)TbL (mm)
0.8–640–574.13–5.94.5–6.245.95–8.27.17–9.67
1.16–3.740–563.9–5.44.2–5.65.57–7.766.64–8.8
2.04
±
0.63 45.61±3.015±0.315.34±0.37.07±0.468.39±0.45
2.01
±
0.56 46.9±3.734.62±0.294.95±0.36.52±0.417.63±0.44
98/79119/95118/95118/95118/95118/95
0.7480.005<0.0001<0.0001<0.0001<0.0001
0.107.8078.2284.0979.61152.51


### 2.4. Meteorological Data

The climate elements considered were average annual temperature and average annual precipitation within the natural locations of the animals. Spatially interpolated climate data [50] for the 1950–2000 period, on grids with spatial resolution of 30 arc-second (often referred to as 1-km spatial resolution), were taken from http://www.worldclim.org and imported into a GIS application.


### 2.5. Data Analyses

I applied one-way ANOVA with multiple testing corrections (Duncan’s multiple range test, MRT) to examine the presence of sexual size dimorphism (SSD) in body and limb dimensions in all species studied. I used simple ordinary least squares (OLS) regression model, treating temperature and precipitation as continuous predictors, to assess the relationship between abiotic environmental components and morphological traits. In this analysis, the examined samples suppose a great dispersion of the data that cannot be addressed with log-transformations. In fact, both body size ranges within the samples and the sample sizes across the observed locations vary among the species. However, prior to analyses the data were tested for normality (Shapiro-Wilk’s test) and homogeneity of variances (Levene’s test). As argued in the Introduction, testing the possible presence of intraspecific clines requires that sexual size dimorphism within the species studied to be taken into account. In view of these considerations, relationships between abiotic variables and the measurements/ratios were tested for each sex separately. I also controlled for the effect of temperature and precipitation on limb length (for each sex separately) by MANCOVA using SVL as a covariate. In this analysis, both abiotic factors treated separately as categorical predictors and total fore limb length (FLL) and total hind limb length (HLL) as dependent variables. Multiple ordinary least squares (OLS) regression analysis, using SVL as a dependent variable and treating climatic factors as predictors, was applied to examine the combined effect of temperature and precipitation on body size. Here, in order to produce a more meaningful index of the observed effects, the slope coefficients were transformed into standardized “beta” coefficients (β). All statistical analyses were done using Statistica 8 (StatSoft, Inc., USA).


## 2.1. Samples and Site Locations

The observed specimens were adult males and females of four lizard species (more than 100 individuals per species), with large geographic ranges. I measured 679 preserved specimens of four lacertid species in the Natural History Museum of the Department of Zoology, Faculty of Life Sciences, Tel Aviv University (TAUM). The specimens had been collected between 1941 and 2012, mostly 86.6% from Israel but some 71 from Egypt (the Sinai Peninsula) and 20 specimens from Turkey.The following specimens were collected:Phoenicolacerta laevis, 104 specimens were collected from Mt. Hermon, the Golan Heights, the Galilee, the foothills of Judea, and the coastal plain (Israel);Ophisops elegans, 127 specimens were collected from Israel (Mt. Hermon, the Golan Heights, the Galilee, the Judean foothills, and the northern and central Negev) and 20 specimens from Turkey;Acanthodactylus boskianus, 164 specimens were collected from Israel (the A’rava Valley and the northern, central, and southern Negev) and 50 specimens from Egypt (northern Sinai, Sinai mountains, and the western Sinai coastal plain);Mesalina guttulata, 193 specimens were collected from Israel (the Judean Hills, the Judean Desert, Dead Sea area, and the northern, central, and southern Negev) and 21 specimens from Egypt (the northern Sinai and Sinai mountains).


## 2.2. Study Species

Phoenicolacerta laevis (Gray, 1838) is a medium-sized lizard, SVL up to 77 mm (according to our records), generally distributed in the Mediterranean landscape from Turkey to Israel and Jordan. This fairly adaptable species is found in relatively humid biotopes, such as temperate forests, valleys, Mediterranean-type shrub vegetation, cultivated fields, and gardens ([45, 46] pers. obs.). It is found from 200 m below sea level up to 2200 m above sea level (TAUM records).Ophisops elegans (Ménétriés, 1832) is a small, slender lizard. The largest specimen in our sample has a SVL of 59 mm. It occurs in open arid plains of the Mediterranean region and Central Asia and was found in sparse forests, vegetated rocky hillsides, and wadis [46–48]. Its vertical distribution ranges from 240 m below sea level to 2000 m above sea level (TAUM records).Acanthodactylus boskianus (Daudin, 1802) is a medium-sized lizard, with adults up to 88 mm SVL. It is widespread in desert and semidesert areas of the Saharo-Arabian region from northern Africa through the eastern Mediterranean and the Arabian Peninsula to Iraq. It inhabits rather hard substrate and occupies a wide range of sandy and gravel habitats. It was collected from 400 m below sea level to 1300 m above sea level [46, 48, 49].Mesalina guttulata (Lichtenstein, 1823) is a rather small, slender lizard, with a maximum SVL of 57 mm. It is widely distributed in North Africa and south-west Asia. It usually inhabits desert and semidesert regions with various terrains, prefers hard substrate and scattered rocks, and avoids sand dunes [46–48]. It is found from 350 m below sea level to 2000 m above sea level (TAUM records).


## 2.3. Morphological Measurements

For each specimen, I measured the following: snout-vent length (SVL); humerus length (HmL), distance from axilla to apex of elbow; forearm length (FaL), distance from elbow apex to apex of wrist; femur length (FeL), distance from groin to apex of knee; and tibia length (TbL), distance from apex of knee to heel. From these, total fore limb length (FLL) and total hind limb length (HLL), which were obtained as the sums of the corresponding segment lengths (i.e., HmL + FaL and FeL + TbL), as well as the limbs-body ratios (HmL/SVL, FaL/SVL, FeL/SVL, TbL/SVL) were calculated. Tail length was not included in the analysis because many individuals had broken or regenerated tails. For each specimen I also recorded collection data (region and locality with geographic coordinates and elevation), as well as body weight (M), which was taken shortly after capture to the nearest 0.1 g. Sex was determined by tail base shape and hemipenis eversion or, in doubtful cases, by dissection. Measurements were made with a digital caliper (Mitutoyo CD-6′′CX) to the nearest 0.01 mm, except for SVL, which was measured to the nearest 0.5 mm using a millimeter ruler. I could not record all measurements from each specimen due to the bad preservation conditions of some reptiles, so the sample sizes vary among features (Table 1). All measurements were made by the author.Table 1
Morphological characteristics of lizard species studied and their intersexual comparisons (one-way ANOVA with multiple testing corrections; Duncan’s multiple range test, MRT).
Species
Features
Range
Mean ± SD
n
♂
♂
/
♀
♀
P
F
♂
♂
♀
♀
♂
♂
♀
♀
Phoenicolacerta laevis
M (g)SVL (mm)HmL (mm)FaL (mm)FeL (mm)TbL (mm)
2.5–10.550.5–775.3–95.9–9.327.62–11.848.75–13
3–11.450–765.22–7.345.9–7.87.7–10.58.5–11.4
5.85
±
1.91 63.53±6.65 6.99±0.797.61±0.779.66±1.0311.04±1.02
5.13
±
1.94 62.97±6.016.34±0.476.96±0.468.78±0.639.87±0.62
48/5648/5647/5647/5647/5647/56
0.0590.6540.00010.00010.00010.0001
3.620.2026.3628.4927.4051.01
Ophisops elegans
M (g)SVL (mm)HmL (mm)FaL (mm)FeL (mm)TbL (mm)
0.9–4.240–594.45–6.14.78–6.725.3–8.727.15–10.53
1.2–3.540–53.53.9–5.44.23–6.225.2-7.476.46–9.22
2.5
±
0.72 48.07±4.275.3±0.45.82±0.57.22±0.719±0.81
1.98
±
0.51 46.47±3.444.66±0.355.15±0.416.29±0.577.78±0.63
85/6285/6284/6184/6184/6184/61
<0.00010.016<0.0001<0.0001<0.0001<0.0001
23.225.8598.3873.2069.1195.38
Acanthodactylus boskianus
M (g)SVL (mm)HmL (mm)FaL (mm)FeL (mm)TbL (mm)
1.7–13.852–885.2–9.185.92–9.928.45–13.3810.8–16.4
3–12.150–80.54.62–7.35.8–8.488.3–12.2810–14.54
6.99
±
2.15 65.9±6.466.9±0.757.96±0.7911.2±1.0413.77±1.19
5.35
±
1.63 61.11±5.836.26±0.567.12±0.629.81±0.7512.12±0.87
121/93121/93113/85113/85114/85114/85
<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001
37.2331.3143.3365.57108.62116.06
Mesalina guttulata
M (g)SVL (mm)HmL (mm)FaL (mm)FeL (mm)TbL (mm)
0.8–640–574.13–5.94.5–6.245.95–8.27.17–9.67
1.16–3.740–563.9–5.44.2–5.65.57–7.766.64–8.8
2.04
±
0.63 45.61±3.015±0.315.34±0.37.07±0.468.39±0.45
2.01
±
0.56 46.9±3.734.62±0.294.95±0.36.52±0.417.63±0.44
98/79119/95118/95118/95118/95118/95
0.7480.005<0.0001<0.0001<0.0001<0.0001
0.107.8078.2284.0979.61152.51


## 2.4. Meteorological Data

The climate elements considered were average annual temperature and average annual precipitation within the natural locations of the animals. Spatially interpolated climate data [50] for the 1950–2000 period, on grids with spatial resolution of 30 arc-second (often referred to as 1-km spatial resolution), were taken from http://www.worldclim.org and imported into a GIS application.


## 2.5. Data Analyses

I applied one-way ANOVA with multiple testing corrections (Duncan’s multiple range test, MRT) to examine the presence of sexual size dimorphism (SSD) in body and limb dimensions in all species studied. I used simple ordinary least squares (OLS) regression model, treating temperature and precipitation as continuous predictors, to assess the relationship between abiotic environmental components and morphological traits. In this analysis, the examined samples suppose a great dispersion of the data that cannot be addressed with log-transformations. In fact, both body size ranges within the samples and the sample sizes across the observed locations vary among the species. However, prior to analyses the data were tested for normality (Shapiro-Wilk’s test) and homogeneity of variances (Levene’s test). As argued in the Introduction, testing the possible presence of intraspecific clines requires that sexual size dimorphism within the species studied to be taken into account. In view of these considerations, relationships between abiotic variables and the measurements/ratios were tested for each sex separately. I also controlled for the effect of temperature and precipitation on limb length (for each sex separately) by MANCOVA using SVL as a covariate. In this analysis, both abiotic factors treated separately as categorical predictors and total fore limb length (FLL) and total hind limb length (HLL) as dependent variables. Multiple ordinary least squares (OLS) regression analysis, using SVL as a dependent variable and treating climatic factors as predictors, was applied to examine the combined effect of temperature and precipitation on body size. Here, in order to produce a more meaningful index of the observed effects, the slope coefficients were transformed into standardized “beta” coefficients (β). All statistical analyses were done using Statistica 8 (StatSoft, Inc., USA).


## 3. Results

### 3.1. Sexual Dimorphism

Summary data of basic morphological features for males and females of each species are given in Table1. Simple intersexual comparisons (one-way ANOVA) of the lizards’ external morphology indicate that maleO. elegansandA. boskianus are significantly (P<0.0001) longer and heavier than females (Table 1). InM. guttulata females are longer but not heavier than males. InP. laevis males are marginally, nonsignificantly heavier than females but not longer.Moreover, there is a pronounced sexual dimorphism (P<0.0001) in extremity length. Within all the species studied in absolute terms (mm) HmL, FaL, FeL, and TbL in males were much greater (from 7.3 to 13.6%) than those in conspecific females. Assessing the ratios of fore and hind limbs’ segments relative to SVL (i.e., HmL/SVL, FaL/SVL, FeL/SVL, and TbL/SVL) made very little difference (data not shown).


### 3.2. Environment-Related Variation

P. laevis. The studied locations are characterized by a highly significant correlation between latitude and both temperature (r=-0.607, P<0.0001) and precipitation (r=0.879, P<0.0001), and between temperature and precipitations as well (r=-0.732, P<0.0001). Considering body size gradient, the scatter plots shown in Figure 2 indicate that although the individuals exhibit a notable dispersion (it is also true for other species studied) the major trends become evident. Male SVL does not vary significantly in relation to climate (temperature: OLS slope ± SE = -0.25±0.45, r=-0.081, P=0.585; precipitation: slope = -0.001±0.006, r=-0.032, P=0.827), whereas females are greatly affected by it, growing significantly heavier and longer (Figures 2(a) and 2(b)) in cool and wet habitats (temperature: slope = -1.47±0.31, r=-0.543, P<0.0001; precipitation: slope = 0.01±0.004, r=0.369, P=0.005). Continuous sampling (Table 2) revealed a stronger effect of climate on corporeal proportions than on absolute sizes. Analysis of latitudinal variability in body proportions showed a significant relationship between temperature and relative extremities length. In both sexes, relatively larger limbs are associated with hot habitats; among females, the recorded temperature-related correlation was highly significant for these morphological variables. The influence of precipitation was less obvious, although in females it was still significant for all ratios. In males, only HmL/SVL varied in relation to precipitation. With SVL as a covariate, however, a MANCOVA indicates a statistically significant effect of temperature on both FLL and HLL in females (Table 3). The results of multiple regression with SVL as a response (Table 4) indicate that the body size-environment link is more obvious in respect to temperature than to precipitation. Thus, in females SVL is much more affected by the former variable: P=0.0009, β=-0.61±0.17, than by the latter: P=0.588, β=-0.09±0.17.Table 2
Climate-morphology relationship between average annual temperature (AAT)/average annual precipitation (AAP) and body measurements/ratios (regression analysis with temperature and precipitation as the predictors and features as the responses). Positive correlation (+) and negative correlation (−); single sign: significant difference (P < 0.05) and double sign: highly significant difference (P < 0.0001).
Features
Phoenicolacerta laevis
Ophisops elegans
Acanthodactylus boskianus
Mesalina guttulata
AAT
AAP
AAT
AAP
AAT
AAP
AAT
AAP
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
M
−
+
−
+
+
+
+
SVL
− −
+
−
−
+
+
−
+
+
+
HmL
−
−
+
++
FaL
−
−
+
+
+
FeL
−
−
++
+
TbL
−
−
+
+
HmL/SVL
+
++
−
−
+
−
−
+
FaL/SVL
+
++
−
+
− −
−
+
FeL/SVL
+
++
− −
+
+
+
−
−
+
−
−
TbL/SVL
+
++
−
+
−
−
+
−
−Table 3
Results of a MANCOVA examining the effect of temperature and precipitation on total fore limb length (FLL) and total hind limb length (HLL) inPhoenicolacerta laevis with SVL as a covariate. Significant effects are marked in bold letter.
Sourceof variation
Dependent variable
df
MS
F
P
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
Temperature
FLLHLL
77
88
0.521.11
0.791.52
1.451.49
2.433.33
0.2140.199
0.028 0.004
Covariate (SVL)
FLLHLL
11
11
73.21127.5
23.4947.34
201.7170.7
71.64103.8
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
0.332.76
4.047.66
0.933.71
12.3216.79
0.3410.062
<0.001<0.001
Error
FLLHLL
3737
4646
0.360.74
0.330.45
Precipitation
FLLHLL
88
88
0.470.99
0.641.50
1.261.30
1.873.27
0.2940.275
0.0880.005
Covariate (SVL)
FLLHLL
11
11
71.08120.8
21.6246.77
191.3158.2
62.75101.7
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
0.302.96
4.937.98
0.813.88
14.3317.38
0.3740.056
<0.001<0.001
Error
FLLHLL
3636
4646
0.370.76
0.340.45Table 4
The effect of climate variables on snout-vent length (SVL). Results of multiple regression with SVL as a dependent variable and climatic factors as predictors.
Effect of temperature
Effect of precipitation
β  ± SE
t
P
β  ± SE
t
P
Phoenicolacerta laevis
♂♂
-
0.21
±
0.21
−0.98
0.328
-
0.18
±
0.21
−0.85
0.399
♀♀
-
0.615
±
0.17
−3.51
0.0009
-
0.09
±
0.17
−0.54
0.588
Ophisops elegans
♂♂
-
0.27
±
0.11
−2.51
0.013
0.26
±
0.11
2.38
0.019
♀♀
-
0.32
±
0.12
−2.66
0.009
0.25
±
0.12
2.08
0.041
Acanthodactylus boskianus
♂♂
-
0.03
±
0.09
−0.34
0.727
0.136
±
0.09
1.37
0.207
♀♀
-
0.13
±
0.11
−1.15
0.251
0.21
±
0.11
−1.77
0.079
Mesalina guttulata
♂♂
-
0.03
±
0.09
−0.41
0.68
0.21
±
0.09
2.36
0.019
♀♀
-
0.14
±
0.09
−1.51
0.134
0.33
±
0.09
3.39
0.001The relationship between SVL and (a) average annual temperature and (b) average annual precipitation inPhoenicolacerta laevis: solid line with circles: males; dashed line with triangles: females.
(a)
(b)O. elegans. Examination of the environmental relationship across the observed area of this mostly Mediterranean species demonstrates a highly significant correlation between latitude and both temperature (r=-0.801, P<0.0001) and precipitation (r=0.433, P<0.0001). Temperature is negatively correlated with precipitation (r=-0.388, P<0.0001). Snout-vent lengths (Figures 3(a) and 3(b)) of both sexes increase with temperature (males: slope = -0.72±0.18, r=-0.402, P=0.0001; females: slope = -0.49±0.14, r=-0.408, P=0.001) and with precipitation (males: slope = 0.007±0.001, r=0.395, P=0.0002; females: slope = 0.005±0.001, r=0.374, P=0.002).The relationship between SVL and (a) average annual temperature and (b) average annual precipitation inOphisops elegans: solid line with circles: males; dashed line with triangles: females.
(a)
(b)The corporeal ratios and, especially, lengths (Table2) are noticeably influenced by the local climate conditions. Regression analysis indicates that the sexes generally follow similar trends in their corporal dimensions and become significantly larger at higher latitudes. By analogy with SVL, appendages were positively correlated with precipitation and were negatively correlated with temperature. In addition, in both sexes the longer femora may be associated with increased precipitation. However, after accounting for SVL, a MANCOVA reveals a nonsignificant effect of both abiotic predictors on limb size (Table 5). Multiple regression demonstrates the equivalency of two climate influences in their relationships with SVL (Table 4). Furthermore, males and females, in general, did not differ from one another in their responses to either temperature (β=-0.27±0.11 versus β=-0.32±0.12) or precipitation (β=0.26±0.11 versus β=0.25±0.12) gradient.Table 5
Results of a MANCOVA examining the effect of temperature and precipitation on total fore limb length (FLL) and total hind limb length (HLL) inOphisops elegans with SVL as a covariate.
Sourceof variation
Dependent variable
df
MS
F
P
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
Temperature
FLLHLL
1717
1616
0.120.48
0.260.71
0.811.21
1.401.47
0.6720.278
0.1860.155
Covariate (SVL)
FLLHLL
11
11
25.2256.96
5.2915.25
163.3143.2
28.6431.83
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
1.011.38
2.512.65
6.513.47
13.595.52
0.0130.067
0.0010.023
Error
FLLHLL
6565
4242
0.150.39
0.180.47
Precipitation
FLLHLL
2323
1717
0.130.56
0.240.64
0.831.56
1.291.30
0.6730.086
0.2430.241
Covariate (SVL)
FLLHLL
11
11
21.8649.02
5.3414.81
140.6136.6
28.2229.71
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
1.232.01
2.633.05
7.965.59
13.936.12
0.0060.021
0.0010.018
Error
FLLHLL
5959
4141
0.150.35
0.180.49A. boskianus. In the desert habitats of this species temperature was uncorrelated with latitude (r=-0.082, P=0.247) across the observed geographic range. There is, however, a pronounced latitudinal precipitation gradient (r=0.815, P<0.0001) and a negative relationship between annual averages of temperature and precipitation: r=-0.367, P<0.0001. In this lizard, the sexes differ in respect to their environmental-morphological connection (Figures 4(a) and 4(b)). The SVL of males did not correlate consistently with temperature (slope = −0.21 ± 0.29, r=-0.067, P=0.485) and precipitation (slope = 0.01±0.008, r=0.134, P=0.160), whereas the SVL of females varied in relation to these variables: slope = -0.81±0.34, r=-0.243, P=0.021, and slope = 0.02±0.007, r=0.280, P=0.007, for temperature and precipitation, respectively.The relationship between SVL and (a) average annual temperature and (b) average annual precipitation inAcanthodactylus boskianus: solid line with circles: males; dashed line with triangles: females.
(a)
(b)Corporeal measurements (Table2) showed a close relationship with rainfall in both sexes and increased precipitation is correlated to shorter extremities relative to SVL. Despite the fact that the appendages tend to lengthen relative to their SVL with rise in temperature, this abiotic component would seem to affect only hind limb segments (FeL and TbL). The MANCOVA results (Table 6) also suggest that within both sexes, the temperature gradient significantly contributes only to HLL variation. Precipitation had a significant effect on hind extremity in both sexes; furthermore, within males it influenced the variation in FLL as well. Analysis of combined effect of environmental variables when both temperature and precipitation were regressed against SVL (Table 4) reveals the role of moisture as a major predictor of intraspecific variability. Moreover, males and females were similar in this trend: β=-0.03±0.09 (temperature), β=0.136±0.09 (precipitation) for males and β=0.13±0.11 versus β=0.21±0.11 for females, respectively.Table 6
Results of a MANCOVA examining the effect of temperature and precipitation on total fore limb length (FLL) and total hind limb length (HLL) inAcanthodactylus boskianus with SVL as a covariate. Significant effects are marked in bold letter.
Sourceof variation
Dependent variable
df
MS
F
P
♂♂
♀♀
♂♂
♀♀
♂♂
♀♀
♂♂
♀♀
Temperature
FLLHLL
3535
2828
0.891.51
0.531.31
1.411.91
1.251.75
0.1180.012
0.2390.040
Covariate (SVL)
FLLHLL
11
11
60.44173.1
34.3661.25
95.26218.6
80.1881.76
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
2.034.45
4.9529.11
3.205.62
11.5638.85
0.0780.021
0.001<0.001
Error
FLLHLL
6565
5151
0.630.79
0.420.74
Precipitation
FLLHLL
3333
3131
0.861.51
0.461.16
1.781.86
1.011.71
0.023 0.016
0.4810.049
Covariate (SVL)
FLLHLL
11
11
64.15179.2
23.9442.22
131.8220.1
51.4361.51
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
1.614.22
7.0434.82
3.325.18
15.1450.73
0.0730.026
<0.001<0.001
Error
FLLHLL
6767
4848
0.480.81
0.460.68M. guttulata. The studied geographic range of this desert lacertid is characterized by a highly significant correlation between latitude and both temperature (r=0.657, P<0.0001) and precipitation (r=0.682, P<0.0001), whereas temperature was uncorrelated with precipitation (r=0.030, P=0.654). In both sexes SVL (Figures 5(a) and 5(b)) increased along a rainfall gradient: slope = 0.006±0.002, r=0.215, P=0.018, for males; slope = 0.01 ± 0.004, r=0.317, P=0.001, for females. Continuous sampling also revealed an insignificant relationship between temperature and SVL in this species. Regression analysis (Table 2) indicated an important role of precipitation in the creation of clines in body and limb dimensions. Moreover, spatial changes of both climatic components may induce some allometric shifts in latitudinal space. Thus, temperature rise has led to the appearance of longer (relative to SVL) extremities among females. This observation was confirmed by MANCOVA; a significant association was observed between either FLL or HLLand temperature, as well as between HLL and precipitation. In males, this approach, however, indicates a nonsignificant effect of both climatic variables on limb dimensions (Table 7).Table 7
Results of a MANCOVA examining the effect of temperature and precipitation on total fore limb length (FLL) and total hind limb length (HLL) inMesalina guttulata with SVL as a covariate. Significant effects are marked in bold letter.
Sourceof variation
Dependent variable
df
MS
F
P
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
Temperature
FLLHLL
2525
2424
0.210.45
0.240.54
1.561.43
1.811.91
0.0650.110
0.029 0.019
Covariate (SVL)
FLLHLL
11
11
20.0844.05
9.8720.11
147.5140.7
72.7971.14
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
4.5010.82
8.9220.73
33.1134.55
65.8373.31
<0.001<0.001
<0.001<0.001
Error
FLLHLL
9191
6969
0.130.31
0.130.28
Precipitation
FLLHLL
2828
2525
0.190.41
0.220.56
1.421.29
1.622.11
0.1070.181
0.0590.008
Covariate (SVL)
FLLHLL
11
11
19.5041.81
9.6319.08
140.9130.7
68.6670.94
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
3.569.01
8.0219.18
25.7528.17
57.1771.33
<0.001<0.001
<0.001<0.001
Error
FLLHLL
8888
6868
0.130.32
0.140.26The relationship between SVL and (a) average annual temperature and (b) average annual precipitation inMesalina guttulata: solid line with circles: males; dashed line with triangles: females.
(a)
(b)Similarly to other desert dwellers (A. boskianus), in this species too precipitation level is a stronger ecological determinant of body length (Table 4) than temperature gradient: P=0.019, β=0.21±0.09, versus P=0.680, β=-0.03±0.09, for males; P=0.001, β=0.33±0.09, versus P=0.134, β=-0.14±0.09, for females, respectively.


## 3.1. Sexual Dimorphism

Summary data of basic morphological features for males and females of each species are given in Table1. Simple intersexual comparisons (one-way ANOVA) of the lizards’ external morphology indicate that maleO. elegansandA. boskianus are significantly (P<0.0001) longer and heavier than females (Table 1). InM. guttulata females are longer but not heavier than males. InP. laevis males are marginally, nonsignificantly heavier than females but not longer.Moreover, there is a pronounced sexual dimorphism (P<0.0001) in extremity length. Within all the species studied in absolute terms (mm) HmL, FaL, FeL, and TbL in males were much greater (from 7.3 to 13.6%) than those in conspecific females. Assessing the ratios of fore and hind limbs’ segments relative to SVL (i.e., HmL/SVL, FaL/SVL, FeL/SVL, and TbL/SVL) made very little difference (data not shown).


## 3.2. Environment-Related Variation

P. laevis. The studied locations are characterized by a highly significant correlation between latitude and both temperature (r=-0.607, P<0.0001) and precipitation (r=0.879, P<0.0001), and between temperature and precipitations as well (r=-0.732, P<0.0001). Considering body size gradient, the scatter plots shown in Figure 2 indicate that although the individuals exhibit a notable dispersion (it is also true for other species studied) the major trends become evident. Male SVL does not vary significantly in relation to climate (temperature: OLS slope ± SE = -0.25±0.45, r=-0.081, P=0.585; precipitation: slope = -0.001±0.006, r=-0.032, P=0.827), whereas females are greatly affected by it, growing significantly heavier and longer (Figures 2(a) and 2(b)) in cool and wet habitats (temperature: slope = -1.47±0.31, r=-0.543, P<0.0001; precipitation: slope = 0.01±0.004, r=0.369, P=0.005). Continuous sampling (Table 2) revealed a stronger effect of climate on corporeal proportions than on absolute sizes. Analysis of latitudinal variability in body proportions showed a significant relationship between temperature and relative extremities length. In both sexes, relatively larger limbs are associated with hot habitats; among females, the recorded temperature-related correlation was highly significant for these morphological variables. The influence of precipitation was less obvious, although in females it was still significant for all ratios. In males, only HmL/SVL varied in relation to precipitation. With SVL as a covariate, however, a MANCOVA indicates a statistically significant effect of temperature on both FLL and HLL in females (Table 3). The results of multiple regression with SVL as a response (Table 4) indicate that the body size-environment link is more obvious in respect to temperature than to precipitation. Thus, in females SVL is much more affected by the former variable: P=0.0009, β=-0.61±0.17, than by the latter: P=0.588, β=-0.09±0.17.Table 2
Climate-morphology relationship between average annual temperature (AAT)/average annual precipitation (AAP) and body measurements/ratios (regression analysis with temperature and precipitation as the predictors and features as the responses). Positive correlation (+) and negative correlation (−); single sign: significant difference (P < 0.05) and double sign: highly significant difference (P < 0.0001).
Features
Phoenicolacerta laevis
Ophisops elegans
Acanthodactylus boskianus
Mesalina guttulata
AAT
AAP
AAT
AAP
AAT
AAP
AAT
AAP
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
M
−
+
−
+
+
+
+
SVL
− −
+
−
−
+
+
−
+
+
+
HmL
−
−
+
++
FaL
−
−
+
+
+
FeL
−
−
++
+
TbL
−
−
+
+
HmL/SVL
+
++
−
−
+
−
−
+
FaL/SVL
+
++
−
+
− −
−
+
FeL/SVL
+
++
− −
+
+
+
−
−
+
−
−
TbL/SVL
+
++
−
+
−
−
+
−
−Table 3
Results of a MANCOVA examining the effect of temperature and precipitation on total fore limb length (FLL) and total hind limb length (HLL) inPhoenicolacerta laevis with SVL as a covariate. Significant effects are marked in bold letter.
Sourceof variation
Dependent variable
df
MS
F
P
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
Temperature
FLLHLL
77
88
0.521.11
0.791.52
1.451.49
2.433.33
0.2140.199
0.028 0.004
Covariate (SVL)
FLLHLL
11
11
73.21127.5
23.4947.34
201.7170.7
71.64103.8
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
0.332.76
4.047.66
0.933.71
12.3216.79
0.3410.062
<0.001<0.001
Error
FLLHLL
3737
4646
0.360.74
0.330.45
Precipitation
FLLHLL
88
88
0.470.99
0.641.50
1.261.30
1.873.27
0.2940.275
0.0880.005
Covariate (SVL)
FLLHLL
11
11
71.08120.8
21.6246.77
191.3158.2
62.75101.7
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
0.302.96
4.937.98
0.813.88
14.3317.38
0.3740.056
<0.001<0.001
Error
FLLHLL
3636
4646
0.370.76
0.340.45Table 4
The effect of climate variables on snout-vent length (SVL). Results of multiple regression with SVL as a dependent variable and climatic factors as predictors.
Effect of temperature
Effect of precipitation
β  ± SE
t
P
β  ± SE
t
P
Phoenicolacerta laevis
♂♂
-
0.21
±
0.21
−0.98
0.328
-
0.18
±
0.21
−0.85
0.399
♀♀
-
0.615
±
0.17
−3.51
0.0009
-
0.09
±
0.17
−0.54
0.588
Ophisops elegans
♂♂
-
0.27
±
0.11
−2.51
0.013
0.26
±
0.11
2.38
0.019
♀♀
-
0.32
±
0.12
−2.66
0.009
0.25
±
0.12
2.08
0.041
Acanthodactylus boskianus
♂♂
-
0.03
±
0.09
−0.34
0.727
0.136
±
0.09
1.37
0.207
♀♀
-
0.13
±
0.11
−1.15
0.251
0.21
±
0.11
−1.77
0.079
Mesalina guttulata
♂♂
-
0.03
±
0.09
−0.41
0.68
0.21
±
0.09
2.36
0.019
♀♀
-
0.14
±
0.09
−1.51
0.134
0.33
±
0.09
3.39
0.001The relationship between SVL and (a) average annual temperature and (b) average annual precipitation inPhoenicolacerta laevis: solid line with circles: males; dashed line with triangles: females.
(a)
(b)O. elegans. Examination of the environmental relationship across the observed area of this mostly Mediterranean species demonstrates a highly significant correlation between latitude and both temperature (r=-0.801, P<0.0001) and precipitation (r=0.433, P<0.0001). Temperature is negatively correlated with precipitation (r=-0.388, P<0.0001). Snout-vent lengths (Figures 3(a) and 3(b)) of both sexes increase with temperature (males: slope = -0.72±0.18, r=-0.402, P=0.0001; females: slope = -0.49±0.14, r=-0.408, P=0.001) and with precipitation (males: slope = 0.007±0.001, r=0.395, P=0.0002; females: slope = 0.005±0.001, r=0.374, P=0.002).The relationship between SVL and (a) average annual temperature and (b) average annual precipitation inOphisops elegans: solid line with circles: males; dashed line with triangles: females.
(a)
(b)The corporeal ratios and, especially, lengths (Table2) are noticeably influenced by the local climate conditions. Regression analysis indicates that the sexes generally follow similar trends in their corporal dimensions and become significantly larger at higher latitudes. By analogy with SVL, appendages were positively correlated with precipitation and were negatively correlated with temperature. In addition, in both sexes the longer femora may be associated with increased precipitation. However, after accounting for SVL, a MANCOVA reveals a nonsignificant effect of both abiotic predictors on limb size (Table 5). Multiple regression demonstrates the equivalency of two climate influences in their relationships with SVL (Table 4). Furthermore, males and females, in general, did not differ from one another in their responses to either temperature (β=-0.27±0.11 versus β=-0.32±0.12) or precipitation (β=0.26±0.11 versus β=0.25±0.12) gradient.Table 5
Results of a MANCOVA examining the effect of temperature and precipitation on total fore limb length (FLL) and total hind limb length (HLL) inOphisops elegans with SVL as a covariate.
Sourceof variation
Dependent variable
df
MS
F
P
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
Temperature
FLLHLL
1717
1616
0.120.48
0.260.71
0.811.21
1.401.47
0.6720.278
0.1860.155
Covariate (SVL)
FLLHLL
11
11
25.2256.96
5.2915.25
163.3143.2
28.6431.83
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
1.011.38
2.512.65
6.513.47
13.595.52
0.0130.067
0.0010.023
Error
FLLHLL
6565
4242
0.150.39
0.180.47
Precipitation
FLLHLL
2323
1717
0.130.56
0.240.64
0.831.56
1.291.30
0.6730.086
0.2430.241
Covariate (SVL)
FLLHLL
11
11
21.8649.02
5.3414.81
140.6136.6
28.2229.71
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
1.232.01
2.633.05
7.965.59
13.936.12
0.0060.021
0.0010.018
Error
FLLHLL
5959
4141
0.150.35
0.180.49A. boskianus. In the desert habitats of this species temperature was uncorrelated with latitude (r=-0.082, P=0.247) across the observed geographic range. There is, however, a pronounced latitudinal precipitation gradient (r=0.815, P<0.0001) and a negative relationship between annual averages of temperature and precipitation: r=-0.367, P<0.0001. In this lizard, the sexes differ in respect to their environmental-morphological connection (Figures 4(a) and 4(b)). The SVL of males did not correlate consistently with temperature (slope = −0.21 ± 0.29, r=-0.067, P=0.485) and precipitation (slope = 0.01±0.008, r=0.134, P=0.160), whereas the SVL of females varied in relation to these variables: slope = -0.81±0.34, r=-0.243, P=0.021, and slope = 0.02±0.007, r=0.280, P=0.007, for temperature and precipitation, respectively.The relationship between SVL and (a) average annual temperature and (b) average annual precipitation inAcanthodactylus boskianus: solid line with circles: males; dashed line with triangles: females.
(a)
(b)Corporeal measurements (Table2) showed a close relationship with rainfall in both sexes and increased precipitation is correlated to shorter extremities relative to SVL. Despite the fact that the appendages tend to lengthen relative to their SVL with rise in temperature, this abiotic component would seem to affect only hind limb segments (FeL and TbL). The MANCOVA results (Table 6) also suggest that within both sexes, the temperature gradient significantly contributes only to HLL variation. Precipitation had a significant effect on hind extremity in both sexes; furthermore, within males it influenced the variation in FLL as well. Analysis of combined effect of environmental variables when both temperature and precipitation were regressed against SVL (Table 4) reveals the role of moisture as a major predictor of intraspecific variability. Moreover, males and females were similar in this trend: β=-0.03±0.09 (temperature), β=0.136±0.09 (precipitation) for males and β=0.13±0.11 versus β=0.21±0.11 for females, respectively.Table 6
Results of a MANCOVA examining the effect of temperature and precipitation on total fore limb length (FLL) and total hind limb length (HLL) inAcanthodactylus boskianus with SVL as a covariate. Significant effects are marked in bold letter.
Sourceof variation
Dependent variable
df
MS
F
P
♂♂
♀♀
♂♂
♀♀
♂♂
♀♀
♂♂
♀♀
Temperature
FLLHLL
3535
2828
0.891.51
0.531.31
1.411.91
1.251.75
0.1180.012
0.2390.040
Covariate (SVL)
FLLHLL
11
11
60.44173.1
34.3661.25
95.26218.6
80.1881.76
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
2.034.45
4.9529.11
3.205.62
11.5638.85
0.0780.021
0.001<0.001
Error
FLLHLL
6565
5151
0.630.79
0.420.74
Precipitation
FLLHLL
3333
3131
0.861.51
0.461.16
1.781.86
1.011.71
0.023 0.016
0.4810.049
Covariate (SVL)
FLLHLL
11
11
64.15179.2
23.9442.22
131.8220.1
51.4361.51
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
1.614.22
7.0434.82
3.325.18
15.1450.73
0.0730.026
<0.001<0.001
Error
FLLHLL
6767
4848
0.480.81
0.460.68M. guttulata. The studied geographic range of this desert lacertid is characterized by a highly significant correlation between latitude and both temperature (r=0.657, P<0.0001) and precipitation (r=0.682, P<0.0001), whereas temperature was uncorrelated with precipitation (r=0.030, P=0.654). In both sexes SVL (Figures 5(a) and 5(b)) increased along a rainfall gradient: slope = 0.006±0.002, r=0.215, P=0.018, for males; slope = 0.01 ± 0.004, r=0.317, P=0.001, for females. Continuous sampling also revealed an insignificant relationship between temperature and SVL in this species. Regression analysis (Table 2) indicated an important role of precipitation in the creation of clines in body and limb dimensions. Moreover, spatial changes of both climatic components may induce some allometric shifts in latitudinal space. Thus, temperature rise has led to the appearance of longer (relative to SVL) extremities among females. This observation was confirmed by MANCOVA; a significant association was observed between either FLL or HLLand temperature, as well as between HLL and precipitation. In males, this approach, however, indicates a nonsignificant effect of both climatic variables on limb dimensions (Table 7).Table 7
Results of a MANCOVA examining the effect of temperature and precipitation on total fore limb length (FLL) and total hind limb length (HLL) inMesalina guttulata with SVL as a covariate. Significant effects are marked in bold letter.
Sourceof variation
Dependent variable
df
MS
F
P
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
♂
♂
♀
♀
Temperature
FLLHLL
2525
2424
0.210.45
0.240.54
1.561.43
1.811.91
0.0650.110
0.029 0.019
Covariate (SVL)
FLLHLL
11
11
20.0844.05
9.8720.11
147.5140.7
72.7971.14
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
4.5010.82
8.9220.73
33.1134.55
65.8373.31
<0.001<0.001
<0.001<0.001
Error
FLLHLL
9191
6969
0.130.31
0.130.28
Precipitation
FLLHLL
2828
2525
0.190.41
0.220.56
1.421.29
1.622.11
0.1070.181
0.0590.008
Covariate (SVL)
FLLHLL
11
11
19.5041.81
9.6319.08
140.9130.7
68.6670.94
<0.001<0.001
<0.001<0.001
Factors
FLLHLL
11
11
3.569.01
8.0219.18
25.7528.17
57.1771.33
<0.001<0.001
<0.001<0.001
Error
FLLHLL
8888
6868
0.130.32
0.140.26The relationship between SVL and (a) average annual temperature and (b) average annual precipitation inMesalina guttulata: solid line with circles: males; dashed line with triangles: females.
(a)
(b)Similarly to other desert dwellers (A. boskianus), in this species too precipitation level is a stronger ecological determinant of body length (Table 4) than temperature gradient: P=0.019, β=0.21±0.09, versus P=0.680, β=-0.03±0.09, for males; P=0.001, β=0.33±0.09, versus P=0.134, β=-0.14±0.09, for females, respectively.


## 4. Discussion

Despite the inability of ectotherms to produce significant metabolic heat, the interaction between body size and heat balance based on a set of behavioural, anatomical, and physiological mechanisms [51] that contribute to fitness might explain the peculiarity of ecological drivers in size clines among these animals. The relationship between habitat conditions and morphological variables (body, tail, or limb dimensions) and consequent shifts in phenology due to climate change have been observed across a wide range of reptile taxa [52–54]. It is possible that in diurnal squamates environmental conditions may favour selection for smaller-sized individuals [55] which, due to the ability to achieve more accurate adjustment of their body temperature, have better thermoregulatory capacities [56] and thus display high activity levels. In terrestrial ectotherms, thermoregulation is realized behaviourally by habitat selection and daily activity pattern [57], while at the same time, morphological specialization may also essentially facilitate it [55]. For example, Shine and Madsen [58] pointed out that for relatively small reptile species behavioural thermoregulation seems to be more important, especially in conditions of high diurnal thermal heterogeneity. In this context, desert reptiles inhabiting warmer regions with lower thermoregulatory costs may adjust their body temperature mainly behaviourally and, thus, invest less in adjustment of either body size or shape in relation to ambient temperature. Despite the great data dispersion (as shown in the scatter plots) it was found that among the two “desert” species studied (M. guttulataandA. boskianus) only females of the latter were smaller in hot temperatures, whereas in both sexes of the former, as well as inA. boskianus males, SVL did not correlate with temperature. In the Mediterranean regionP. laevis andO. elegans inhabit cooler climates than the desertA. boskianus andM. guttulata. In these environments, larger gravid reptiles, in addition to the buffering effects of the behavioural thermoregulation supplied by the mother [59, 60], can provide greater temperature stability during the preoviposition period. Volynchik [19], however, found that a large-bodied local viper (Vipera palaestinae), with a Mediterranean distribution pattern, does not follow Bergmann’s rule or its converse. At the same time, the most “cold-loving” lizard (O. elegans) in my own sample obeyed Bergmann’s rule, growing larger in cooler environments. On the other hand, in desert habitats, this pattern of greater thermal inertia seems to be less important, and the behavioural component of thermoregulation becomes paramount.The impact of ambient temperature on external morphology is also significant in respect to relative extremity lengths: different degrees of intraspecific variability in limb proportions supporting Allen’s rule were observed in all the species studied. The present findings indicate that temperature decrease across the observed habitats could lead to reduction in size of extremity segments relative to body length, thus constituting a compensatory adaptation to either cooler or hotter areas. I noted an allometric shift in humerus, forearm, femur, and tibia lengths, relative to SVL, along a temperature gradient amongP. laevis. A similar relationship between average annual temperature and corporeal ratios was noted in females ofM. guttulata. AmongA. boskianus total hind limb length significantly correlated with both thermal regime and rainfall gradient. Hence, both Mediterranean and desert reptiles follow a similar pattern in their body shape variation and may benefit from various morphological modifications matched to habitat range.Body temperature changes affect reptile ecology [61], becoming more important to ectotherms living in cooler environments, which complicate attaining proper body temperature [58, 62]. In support of this, my results show the association between biogeographic origins of the species examined and their relationship with environmental variables. On interspecific comparison, the most pronounced distinctions were noted in phenotypic responses to the temperature gradient between Mediterranean and desert inhabitants, while the between-species differences for precipitation-related variation were weaker. As expected, both Mediterranean region species (P. laevis,O. elegans) display significantly stronger responses to temperature than desert dwellers (A. boskianus,M. guttulata). Apparently, the cooler locations of the former force them to exhibit more consistent morphological variability in respect to temperature changes, whereas for desert species such habitat-morphology connection seems to be less important. It is highly likely that the species’ sensitivity to the environment is associated not only with their ecological differentiation, but also with body size ranges within the sample. Thus, this might be one possible explanation for the different, albeit relatively minor, responses to temperature between ecologically similar species (i.e.,P. laevis versusO. elegans andA. boskianus versusM. guttulata) that were noted here. In the examined samples,L. laevis, owing to its greater measurement ranges thanO. elegans, has a higher potential for more consistent matching of its body dimensions to the observed temperature gradient. The same is true for the two desert lizards;A. boskianus due to its extended body and limb ranges, in comparison withM. guttulata, displays a stronger morphological-environmental link.The examination of gender-specific habitat-morphology interactions revealed that within three of the four species (P. laevis, A. boskianus, and M. guttulata) sensitivity to the environment is coupled with sex. Considering the effect of climate variables on SVL, the findings demonstrate that females, in general, show a more consistent link between temperature and body/limb measurements than males, while the effect of precipitation on external characters was similar for both sexes. Moreover, inP. laevis SVL was significantly more strongly correlated with the temperature gradient in females than in males. These observations indicate the sex-associated effect of natural selection, pointing to different thermal and metabolic demands in males and females. Thus, it is likely that in females the increased body size in colder environments is associated with reproduction-induced requirements. This pattern suggests that females, due to increased thermal requirements, particularly during vitellogenesis [63], are more able to generate a range of suitable morphological modifications for maintenance of optimal body temperature in local environments. Overall, my findings indicate that although females show a more consistent correlation between SVL and a temperature gradient than males, the body shape changes, involving both body and limb dimensions, are rather preliminary in explaining the observed intersexual divergence in the habitat-morphology relationship.Environmental conditions within the nest experienced during early ontogenesis, especially temperature and moisture of substrate, can significantly affect development and growth trajectories of individuals in oviparous reptiles [64–68]. As observed in many species, either posthatching phenotype or offspring fitness are closely associated with thermal regime [63, 69–71]. In both the laboratory and in natural nests [64] temperature is traditionally considered to be an important environmental determinant of reptile size, while the effect of water-related factors is more complex. Several recent studies have tested the effect of substrate moisture during incubation of eggs on a diverse range of reptile species. Significant effects are found in certain species [66, 68, 72, 73] but not in others [74–76] that may reflect either different water permeability of the eggshell or interspecific variation in the water content of eggs [77]. On the other hand, the appearance of morphological variation along a precipitation gradient, which, in desert lizards, is even more evident than the temperature-related variability, can be considered as a response to the local habitats. In reptiles, food availability strongly affects growth trajectories and is one of the possible determinants of body size and shape [78–80]. Overall, in arid regions, where vegetation and arthropod fauna are determined primarily by precipitation [11], wetter localities provide better food resources to lizards at the population level. The present findings suggest that among desert species the possible spatial differences in prey availability and, thus, in nutrition levels, may lead to morphometric differentiation between specimens inhabiting dry or wet sites. However, detailed studies on possible dietary variation linked to local environments are necessary to answer the question of whether the habitat conditions, which are dependent on precipitation, have an effect on body dimensions in these lizards.A nonphylogenetic approach suggests that lizard species tend to exhibit strong morphological specialization in relation to their habitat [53, 81]. Variation in body characters related to locomotion which, in turn, is associated with substrate pattern, has been observed across many lizard taxa [54, 82–84]. Unfortunately, based on available datasets, I could not evaluate whether the shift in substrate usage is associated with limb proportion changes, and so a rational explanation for the effect of vegetation on limb lengths, through its dependence upon hydrothermal regime, remains to be determined. My findings show that, generally, the examined species indicate a uniform pattern in extremities size variation; fore and, especially, hind limb lengths were positively correlated with temperature and negatively with precipitation. The results indicate that ecological adaptation to abiotic factors is a powerful source of extremity size variation, which is under the effect of natural selection. In lizards, the morphology-performance relationships suggest that ground-dwelling species from open landscapes should possess relatively long limb segments to their body length [53]. This body shape pattern increases stride length providing better sprinting abilities [85, 86]. Furthermore, shorter extremities may facilitate the locomotion in highly vegetated areas, where long limbs may be disadvantageous [53]. In this context, my findings suggest that various body shape patterns, such as relatively short/long limbs or their segments, which are suited to the environmental requirements, may reflect an adaptation to either dense or sparse vegetation, resulting from different rainfall intensities across the distribution range.Finally, given the potential importance of habitats for limb size variation, lizards may also use their protruding body parts as additional (together with body surface area) heat exchangers. In these reptiles, for optimal balance between heat gain and heat loss, isometric size changes, to a certain extent, are a “two-edged sword,” in which high thermal inertia of large bodies is accompanied by slow heating rates, and vice versa. In this respect, among desert dwellers relatively long extremities and small size enhance heat exchange and seem to be a morphological adaptation for effective thermoregulation. In desert areas, large daily fluctuations of temperature force the ectotherms to quickly respond to changing conditions and rapidly reach the necessary metabolic rates for locomotor activity. Long-limbed individuals which have an increased surface area-to-volume ratio may enhance their heating and cooling abilities through more effective interactions between the body and either solar radiation or surface temperature gradient. Hence, a consistent link between environmental variables and limb lengths can also be explained by their importance for maintenance of optimal body temperature in local environments. As shown above, the fact that the extremity size-environment link is more obvious in respect to hind limbs than to fore ones reinforces this hypothesis; relatively large hind limbs having an increased heat exchange capacity are much more powerful tool for effective thermoregulation than small fore extremities.In summary, the lacertid lizards I studied generally followed the ecogeographic rules of Bergmann and Allen. In addition, the findings indicate both intersexual and interspecific differences in the habitat-morphology relationship, which can be considered within the framework of their functional connection. My data provide evidence that females tend to show stronger relationships between climate variables and body dimensions than males. In desert species temperature has a smaller effect than rainfall on lizard phenotypes. Under a constant dry climate, water availability seems to be the most important abiotic factor in reptile ecology. In this respect, due to the strong association between precipitation and productivity, my findings could also be interpreted as evidence for the hypothesis that different body size patterns are linked to the spatial variation in primary production, that is, the “resource rule.” In the Mediterranean region lizards the thermal regime is apparently the main driving force of morphological variability. The present study suggests a complex link between environmental variables and morphological characters within lacertid species. The physical conditions, such as temperature and precipitation gradient, were shown to be able to significantly affect both body and limb dimensions, and, thus, as a powerful source of variation, they play an important role in lizard morphology. However, general statements about the habitat-morphology relationships and their potential importance for maintenance of optimal body temperature in reptiles are hard to make due to the enormous diversity of variations in body composition and physiology across species, sexes, and even individuals and the difficulty in evaluating enough mechanisms in enough species. Consequently, additional field and experimental studies testing both body size and shape variations and their relationships with environmental variables, especially at the intraspecific level and accounting for sex-specific aspects are needed to help determine whether the observed patterns are common across species and to what extent they are important for effective thermoregulation.



---

*Source: 795387-2014-05-19.xml*