We quantified four gait characteristics (stride length, stride frequency,step length and floating distance) over a range of running speeds in 11 lacertid lizard species known to vary in maximal sprint speed and microhabitat use. For each species, we measured snout—vent length (SVL), body mass and hindlimb length. We tested which variables determine sprint speed, how each species modulates sprint speed and whether species occurring in different microhabitats differ in speed modulation strategy. In doing so, we aimed to test the assumption that sprint speed is correlated with hindlimb length through its effect on stride length. Variation in maximal sprint speed is determined by variation in both stride length and frequency, but those species that modulate their sprint speed mainly by altering stride length appear to attain the highest maximal speeds. At maximal sprint speed, long-limbed species take larger strides because of the positive effect of hindlimb length on step length and floating distance. However, when hindlimb length is statistically controlled for, mass has a negative effect on step length. None of the three morphological variables explained the interspecific variation in stride frequency at maximal sprint speed. Possibly, differences in physiological properties (e.g. muscle contraction speed) underlie the variation in stride frequency. The 11 species modulate their speed in different ways. Lacertids often seen in vertical microhabitats do not seem to be either pronounced stride length or frequency modulators. These species alter their speed by combining the two strategies. However, species occurring mostly in open and vegetated microhabitats are, respectively, stride length and stride frequency modulators. This difference in running style is substantiated by interspecific morphological differences.

Locomotor performance varies considerably within and among species. In some cases, variation in locomotor performance has been linked to ecological differences (e.g. Garland et al.,1988; Garland and Losos,1994; McPeek et al.,1996; Shipley et al.,1996; Garland,1999; Perry, 1999; Melville and Swain, 2000). For instance, Anolis lizards that differ in sprint performance on dowels of different diameters use differently sized perches in nature(Irschick and Losos, 1999). The relationship between locomotor performance and ecological factors is generally interpreted in evolutionary terms, i.e. differences in selective pressures among species are considered to shape the tight correlation between the two variables. However, the variation in locomotor performance may also be explained by proximate factors. A large body of literature exists on the influence of morphological differences on the variation in performance (e.g. Garland, 1983; Heglund and Taylor, 1988; Emerson, 1991; Srygley and Dudley, 1993; Garland and Losos, 1994; Miles, 1994; Rayner, 1995; Kolok, 1999; Melville and Swain, 2000) and,for terrestrial locomotion, the positive correlation between (hind)limb length and speed seems self-evident (e.g. Garland and Janis, 1993; Garland and Losos, 1994; Bauwens et al.,1995; Bonine and Garland,1999; Melville and Swain,2000). In all these cases, the underlying biomechanical assumption is that long limbs allow animals to take large strides, resulting in high sprint speeds.

We have three problems with the cause/effect assumption that hindlimb length affects stride length, which affects sprint speed. First, this often-cited assumption has never been explicitly tested. Second, stride length is not the sole determinant of sprint speed. Stride length is made up of two components: step length (or the distance the body moves forward during stance)and floating distance [or the distance the body moves forward when the(hind)limbs are not in contact with the substratum]. Moreover, sprint speed is the result of combining a particular stride length with a particular stride frequency. This implies that, in theory, species that are unable to lengthen their stride (for whatever reason) can ultimately attain identical sprint speeds solely by increasing their stride frequency. Moreover, interspecific variation in stride length and stride frequency at maximal sprint speed could originate in different ways. One possibility is that different species employ different speed modulation strategies, i.e. they alter stride length and frequency with speed to different extents. Ultimately, maximal sprint speed will differ between species that alter stride length and frequency with speed in diverse ways because different speed modulation strategies place different limits on the maximal attainable speed. However, all species might employ the same speed modulation strategy, but the increase in stride length or frequency with speed might be limited at different points among species. How far each species can increase its stride length or frequency might be determined by its design (e.g. morphology, physiology, biochemistry). In this case, species will differ in maximal sprint speed. To elucidate the mechanistic basis of the variation in sprint speed, we therefore need information on how each species realises its maximal sprint speed (i.e. speed modulation strategy) and on the values of the gait characteristics at maximal sprint speed.

The third problem with the hindlimb length→stride length→sprint speed assumption is that the way in which animals realise their maximal sprint speed might differ among species occupying different microhabitats. To our knowledge, no quantitative data exist on the relationship between ecological variables and speed modulation. White and Anderson(1994) hypothesised that differences in gait characteristics with speed might reflect differences in foraging strategy (active versus ambush predators) among teiid lizards, and Van Damme et al.(1998) suggested that the way in which lizards realise their ultimate sprint speed may be influenced by predation pressure.

One of the most obvious potential ecological factors likely to influence locomotor style is microhabitat use because different substrata (e.g. open,sandy, patches of vegetation, inclined substrata) will impose different or even conflicting demands on steady-running style. Species from open areas can increase sprint speed by increasing both stride length and frequency; the former by including large floating distances in the stride cycle. Species living among patches of vegetation, in contrast, would be expected to modulate speed by altering stride frequency. Because these species will often have to move around obstacles, they must be able to change direction fast and frequently. In this context, it seems beneficial to increase the number of ground contact phases (and thus stride frequency) as these present the only possibilities for altering the direction of motion. Species occurring on inclined surfaces can increase speed by increasing stride frequency or stride length: in the latter case, only if the change in stride length is due to an increase in step length alone. No floating phases can be incorporated in the stride cycle because losing contact with the surface may have disastrous consequences, for climbing animals

In this study, we aimed to test the hindlimb length→stride length→sprint speed assumption for 11 lacertid lizard species known to vary considerably in maximal sprint speed and microhabitat use. Specifically,we test the following hypotheses: (i) that differences in both stride length and frequency contribute to variation in maximal sprint speed; (ii) that differences in both step length and floating distance contribute to differences in stride length at maximal sprint speed; (iii) that differences in morphology (snout—vent length, SVL, mass, hindlimb length) explain differences in step length, floating distance and stride frequency at maximal sprint speed; (iv) that the variation in maximal sprint speed is due to differences in modulation strategies (stride length versusfrequency); (v) that the way in which stride length changes with speed depends on the way in which both step length and floating distance change with speed;(vi) that morphological differences explain the way in which step length,floating distance and stride frequency change with speed and (vii) that species occurring in different microhabitats use different speed modulation strategies.

Animals

We collected nine of the 11 species during the summers of 1994-1999 at different sites in Western Europe and the Canary Islands (see Table 1 for a list of species). Latastia longicaudata and Takydromus sexlineatus were obtained commercially. Adult males and non-gravid females were used in this study. We transported all lizards to the laboratory at the University of Antwerp, Belgium. Small species (<80 mm SVL) were housed in groups of four or five, with species kept apart, in glass terraria (100 cm×50 cm) with a sandy substratum, some stones and patches of vegetation. Multiple refuges and basking spots were available in each terrarium, and two light bulbs (100 W each) provided light and heat for 10 h per day, enabling the lizards to thermoregulate. The larger species, Gallotia galloti and Lacerta bilineata, were housed individually in identical terraria. We fed the animals live crickets dusted with calcium daily. Water was always present.

Table 1.

Species used in this study, number of individuals for which we obtained data on gait characteristics, number of digitised sequences and maximal sprint speed

SpeciesNumber of individualsNumber of sequencesVmax (m s-1)
Gallotia galloti (Oudart) 41 1.93 
Lacerta oxycephala (Dumeril & Bibron) 47 2.02 
Lacerta bedriagae (Camerano) 10 46 1.79 
Podarcis sicula (Rafinesque) 43 1.67 
Podarcis muralis (Laurenti) 10 48 2.14 
Podarcis hispanica Steindachner 40 2.03 
Podarcis atrata Steindachner 40 1.09 
Podarcis tiliguerta Gmelin 10 47 1.55 
Lacerta bilineata Daudin 26 2.68 
Takydromus sexlineatus Daudin 24 1.33 
Latastia longicaudata (Reuss) 31 3.34 
SpeciesNumber of individualsNumber of sequencesVmax (m s-1)
Gallotia galloti (Oudart) 41 1.93 
Lacerta oxycephala (Dumeril & Bibron) 47 2.02 
Lacerta bedriagae (Camerano) 10 46 1.79 
Podarcis sicula (Rafinesque) 43 1.67 
Podarcis muralis (Laurenti) 10 48 2.14 
Podarcis hispanica Steindachner 40 2.03 
Podarcis atrata Steindachner 40 1.09 
Podarcis tiliguerta Gmelin 10 47 1.55 
Lacerta bilineata Daudin 26 2.68 
Takydromus sexlineatus Daudin 24 1.33 
Latastia longicaudata (Reuss) 31 3.34 

Data are taken from Vanhooydonck and Van Damme(2001) and Vanhooydonck et al.(2001).

Morphometrics

We took the following external measurements with digital callipers(Mitutoyo, CD-15DC) to the nearest 0.01 mm: snout—vent length (SVL),measured from the tip of the snout to the posterior edge of the anal scale;femur length, measured from the hip to the knee joint; tibia length, measured from the knee joint to the ankle; and metatarsal length, measured from the ankle to the base of the third toe. All animals were weighed on an electronic balance (A&D Instruments, FX-3200) to the nearest 0.01 g. The sum of femur length, tibia length and metatarsal length (hereafter referred to as hindlimb length) was used in subsequent analyses. Thus, our measurement of hindlimb length does not correspond to the measure of hindlimb span (e.g. Bonine and Garland, 1999),which takes variation in the width of the pelvis into account.

We calculated the mean values for each species of the three morphological variables (SVL, mass, hindlimb length) and logarithmically(log10)-transformed them before statistical analyses.

Performance

Data on maximal sprint speed of the 11 species are taken from our previous studies (Vanhooydonck and Van Damme,2001; Vanhooydonck et al.,2001) and given in Table 1. We calculated the logarithm (log10) of the mean maximal sprint speed for each species and used these values in subsequent analyses.

Microhabitat use

We observed lizards in a large (4 m×4 m) terrarium containing the following eight types of microhabitat: (i) a dense stand of purple moor grass(Molinia caerulea), with a total area of 3 m2; (ii) a strip of grass (Poa annua), with an area of 3 m2; (iii) a patch of moss with a surface area of 3 m2; (iv) an area of approximately 3 m2 closely planted with shrubs (Loniceraspp.); (v) a strip of sand, also 3 m2 in surface area; (vi) a wall of piled stones, 0.5 m high and with a top surface area of 0.5 m2;(vii) a series of logs, loosely stacked to a height of 0.3 m; and (viii) a tree trunk of 20 cm diameter, cut off at a height of 40 cm. Eight 500 W lamps,one above each type of microhabitat, provided heat and light throughout the experiment.

The smaller species (SVL<80 mm) were introduced into the test terrarium in groups of no more than eight individuals, the larger species in groups of no more than five individuals. Species were always kept apart. To avoid disturbance, all observations were made from behind a blind. After an initial habituation period of 1 h, we recorded the microhabitat occupied by each individual lizard in the terrarium at 30 min intervals. To facilitate interpretation, we decided to pool the types of microhabitat described above into three classes that clearly differ in structural features: an `open'microhabitat [a combination of low grass (ii), moss (iii) and sandy (iv)patches], a `vegetated' microhabitat [a combination of moor grass (i) and shrub (v) patches) and a `vertical' microhabitat [stone wall (vi), logs (vii)and tree trunk (viii)].

For each species, we measured the proportion of time spent in each of the three microhabitats (open, vertical or vegetated). We arcsine-transformed all proportions before subsequent statistical analyses. We used the arcsine transformation because it is the most appropriate when dealing with proportions (Sokal and Rohlf,1995).

We were unable to make observations on microhabitat use in Podarcis hispanica and Podarcis atrata because the experimental arrangement described above was not available when we performed the morphological, performance and kinematic measurements on these two species. All analyses concerning microhabitat use are thus for nine species only.

Gait analysis

To obtain estimates of gait characteristics, we filmed all lizards with a high-speed video camera (NAC 1000) at 500 frames s-1 in dorsal view while they were running on a moving treadmill belt. To facilitate subsequent digitisation, we marked the lizards with small white dots (of non-toxic paint)mid-dorsally on the pelvis and on the dorsal and ventral sides of the second toe of the hindlimb. After marking, the lizards were placed in an incubator set at the different species-selected body temperatures (38 °C for L. longicaudata, 35 °C for the other species) (see Castilla et al., 1999) for at least 1 h. The lizards were then placed, one at a time, on the belt moving at a low speed. Once the animals started running, the speed of the belt was increased slowly but steadily. We therefore obtained video sequences of each individual lizard at a wide range of running velocities.

In an initial analysis, the speed attained during each recorded running bout was assessed by digitising (using a NAC-1000 XY coordinator) a marker on the moving belt, as a reference, and the pelvic marker on the first and last frame of potentially useful sequences. A sequence was considered `potentially useful' only if the lizard took at least five sequential steps at a constant speed and moved in a straight line. For each individual, we then aimed to select (at least) five sequences evenly spread throughout the range of speeds for further detailed analyses. For practical reasons (e.g. narrow range of speeds), this was not always possible. On these sequences, we digitised a marker on the moving belt and the pelvic and toe markers every third frame(i.e. with a time interval of 0.006 s). We chose to digitise a hindfoot because the posterior limbs generate most of the propulsive power in these animals (Snyder, 1952). The numbers of individuals and sequences used for complete digitisation and determination of maximal sprint speed (see Vanhooydonck and Van Damme,2001; Vanhooydonck et al.,2001) are given in Table 1. The speeds at which lizards ran for each digitised sequence and species are presented in Fig. 1.

Fig. 1.

Speeds at which the lizards ran on the treadmill for each species. Each dot represents a digitised sequence. Full species names are given in Table 1.

Fig. 1.

Speeds at which the lizards ran on the treadmill for each species. Each dot represents a digitised sequence. Full species names are given in Table 1.

We obtained the following variables from the digitised images using the method described by Van Damme et al.(1998): (i) stride length, the distance travelled by the centre of mass during one complete cycle of limb movements (cf. Alexander, 1977a,b, 1982), here measured as the distance travelled from the time the hindfoot touches the ground until the time when that same hindfoot touches the ground in the following cycle; (ii)stride frequency, the number of complete cycles per second; (iii) step length,the distance travelled by the centre of mass during the stance of a particular leg (Alexander, 1977a; McMahon, 1984), here measured as the distance travelled from the time the hindfoot touches the ground until the time that same hindfoot loses contact with the substratum in the same cycle; (iv) floating distance, the distance travelled by the centre of mass when neither of the hindfeet is in contact with the ground [(calculated as stride length minus 2(step length)], here measured as the distance travelled from the time both hindfeet lose contact with the substratum until the time one of the hindfeet touches the ground again.

We logarithmically (log10)-transformed running speed and the kinematic variables. For each species, we then regressed each of the four gait characteristics against running speed (see Table 2 for regression equations). We used the regression equation to estimate the values of the four gait characteristics (i.e. stride length, stride frequency, step length and floating distance) at maximal sprint speed (as measured on the moving belt)for each species. These estimated values at maximal sprint speed and the slopes (a in Table 2)of each regression equation were used in subsequent analyses, the former to test which variables determine maximal sprint speed and the latter to test how species realise their maximal sprint speed and whether speed modulation strategies differ among species occupying different microhabitats.

Table 2.

Regression equations describing the relationships between the gait characteristics and running speed of lizards running on the moving belt

Stride lengthStride frequencyStep lengthFloating distance
Speciesabababab
Gallotia galloti 0.40 -1.01 0.60 1.01 0.26 -1.34 0.009 0.001 
Lacerta oxycephala 0.35 -1.09 0.65 1.09 0.27 -1.38 0.006 -0.001 
Lacerta bedriagae 0.42 -0.94 0.58 0.94 0.33 -1.19 0.003 -0.006 
Podaras sicula 0.48 -0.95 0.52 0.95 0.36 -1.24 0.009 0.0001 
Podaras muralis 0.46 -1.11 0.54 1.11 0.43 -1.34 -0.003 -0.006 
Podaras hispanica 0.33 -1.15 0.67 1.15 0.29 -1.41 0.00001 -0.004 
Podaras atrata 0.25 -1.17 0.75 1.17 0.25 -1.40 -0.001 -0.005 
Podaras tiliguerta 0.53 -0.98 0.47 0.98 0.41 -1.23 0.005 -0.006 
Lacerta bilineata 0.38 -0.96 0.62 0.96 0.24 -1.26 0.010 -0.001 
Takydromus sexlineatus 0.31 -1.11 0.69 1.11 0.18 -1.33 0.004 -0.007 
Latastia longicaudata 0.53 -0.97 0.47 0.97 0.37 -1.26 0.015 -0.001 
Stride lengthStride frequencyStep lengthFloating distance
Speciesabababab
Gallotia galloti 0.40 -1.01 0.60 1.01 0.26 -1.34 0.009 0.001 
Lacerta oxycephala 0.35 -1.09 0.65 1.09 0.27 -1.38 0.006 -0.001 
Lacerta bedriagae 0.42 -0.94 0.58 0.94 0.33 -1.19 0.003 -0.006 
Podaras sicula 0.48 -0.95 0.52 0.95 0.36 -1.24 0.009 0.0001 
Podaras muralis 0.46 -1.11 0.54 1.11 0.43 -1.34 -0.003 -0.006 
Podaras hispanica 0.33 -1.15 0.67 1.15 0.29 -1.41 0.00001 -0.004 
Podaras atrata 0.25 -1.17 0.75 1.17 0.25 -1.40 -0.001 -0.005 
Podaras tiliguerta 0.53 -0.98 0.47 0.98 0.41 -1.23 0.005 -0.006 
Lacerta bilineata 0.38 -0.96 0.62 0.96 0.24 -1.26 0.010 -0.001 
Takydromus sexlineatus 0.31 -1.11 0.69 1.11 0.18 -1.33 0.004 -0.007 
Latastia longicaudata 0.53 -0.97 0.47 0.97 0.37 -1.26 0.015 -0.001 

We pooled all sequences per species.

All regressions were significant, except the regressions of floating distance against running speed in some species (given in italic type).

Equations are of the form logy=alogx+b,where x is running speed and y is the gait characteristic. On the basis of these regression equations, we estimated stride length, stride frequency, step length and floating distance at maximal sprint speed.

Statistical analyses

Because species share parts of their evolutionary history, they cannot be regarded as independent data points in statistical analyses (Felsenstein, 1985, 1988; Harvey and Pagel, 1991; Garland et al., 1993). A fundamental assumption in traditional statistical analyses, however, is that data points are independent of each other. In recent years, various computer programs have been developed to take phylogenetic relationships into account in statistical analyses (Harvey and Pagel,1991; Garland et al.,1993; Losos and Miles,1994). In the present study, we used the phylogenetic independent contrasts approach (Felsenstein, 1985, 1988).

The independent contrasts method requires information on the topology of the phylogenetic tree. The phylogeny of the Lacertidae is incompletely resolved, and we had to use a combination of results from mitochondrial DNA(Harris et al., 1998; Harris and Arnold, 1999) and morphological studies (Arnold, 1989, 1998) to compile a tree that is the best approximation currently available(Fig. 2). Some unresolved polytomies remain. However, all attempts at phylogeny reconstruction in lacertids, even those using molecular techniques, have been unable to unravel these relationships, and the polytomies may, in fact, represent explosive speciation events (Arnold,1989; Harris et al.,1998; Harris and Arnold,1999; Fu, 2000). We therefore considered the unresolved nodes to be `hard' polytomies (see Purvis and Garland, 1993). All branch lengths were set to unity, because few data are available on divergence times within lacertids. It has been shown that the actual length of the branches does not usually greatly affect the outcome of phylogenetic statistical analyses (Martins and Garland,1991; Walton,1993; Irschick et al.,1996; Díaz-Uriarte and Garland, 1998). Moreover, checks of branch lengths with the PDTREE program (Garland et al., 1999)did not show any significant correlation between the absolute values of the standardised contrasts and their standard deviations, implying that branch lengths were adequate for all traits.

Fig. 2.

Hypothesised phylogenetic relationships among the 11 lacertid lizard species under study. The phylogeny is incompletely resolved, and this tree should be considered a best approximation. Polytomies are treated as `hard'for the purposes of phylogenetic analyses (see text).

Fig. 2.

Hypothesised phylogenetic relationships among the 11 lacertid lizard species under study. The phylogeny is incompletely resolved, and this tree should be considered a best approximation. Polytomies are treated as `hard'for the purposes of phylogenetic analyses (see text).

For each species, we calculated the independent contrasts of the log-transformed mean values of the morphological (SVL, mass, hindlimb length)and performance (maximal sprint speed) variables, of the arcsine-transformed proportions of the observations on microhabitat use (open, vertical, vegetated microhabitat) and of the species-specific slopes of the regressions of the four gait characteristics (stride length, stride frequency, step length,floating distance) against running speed and of the same gait characteristics at maximal sprint speed (extrapolations based on regression equations, see above).

The independent contrasts were then used as input for a series of backward stepwise multiple regressions through the origin (see Garland et al., 1992). A graphical representation of the different regression analyses, apart from those concerning microhabitat use, is shown in Fig. 3A. In the regressions in which we used the contrasts of the slopes of the gait characteristics against speed, we never introduced the contrasts of the slopes of stride length against speed and the contrasts of the slopes of stride frequency against speed simultaneously in the analyses because these two variables are, by definition, inversely related. To test for the context-dependency of speed modulation among the 11 lizard species, we correlated the contrasts of the three microhabitat variables with the contrasts of the slopes of stride frequency against speed. By definition, introducing the contrasts of the slopes of stride length against speed in the analyses gives the opposite result.

Fig. 3.

Graphical representations of the multiple regression analyses we performed on the morphological, gait characteristics and performance data in the 11 lacertid lizard species. (A) All potential relationships. Each arrow represents one analysis. Independent variables introduced into the same multiple regression are boxed. Arrows are drawn from the independent variables box to the dependent variable box. We performed two series of analyses: one that included the (contrasts of the) values for the four gait characteristics at maximal sprint speed, and the second on the (contrasts of the) slopes of each of the four gait characteristics against speed. (B) Significant relationships in the multiple regression analyses including the contrasts of the estimated values of the four gait characteristics at maximal sprint speed.(C) Significant relationships in the multiple regression analyses including the contrasts of the slopes of each gait characteristic against speed. By definition, the contrasts of the slopes of stride length and stride frequency against speed are inversely related so we never introduced them simultaneously into the analyses. Instead, we included only the contrasts of the slopes of stride frequency against speed in the analysis; including the contrasts of the slopes of stride length against speed gives the opposite result(r=0.56; P=0.075).

Fig. 3.

Graphical representations of the multiple regression analyses we performed on the morphological, gait characteristics and performance data in the 11 lacertid lizard species. (A) All potential relationships. Each arrow represents one analysis. Independent variables introduced into the same multiple regression are boxed. Arrows are drawn from the independent variables box to the dependent variable box. We performed two series of analyses: one that included the (contrasts of the) values for the four gait characteristics at maximal sprint speed, and the second on the (contrasts of the) slopes of each of the four gait characteristics against speed. (B) Significant relationships in the multiple regression analyses including the contrasts of the estimated values of the four gait characteristics at maximal sprint speed.(C) Significant relationships in the multiple regression analyses including the contrasts of the slopes of each gait characteristic against speed. By definition, the contrasts of the slopes of stride length and stride frequency against speed are inversely related so we never introduced them simultaneously into the analyses. Instead, we included only the contrasts of the slopes of stride frequency against speed in the analysis; including the contrasts of the slopes of stride length against speed gives the opposite result(r=0.56; P=0.075).

Gait characteristics at maximal sprint speed

The variation in contrasts of maximal sprint speed is explained, by definition, by the variations in contrasts of stride length and frequency(Fig. 3B). Species that attain a high sprint speed do so by taking larger (standardised regression coefficient r=0.90) and more frequent (r=0.69) strides.

The variation in contrasts of stride length, in turn, is largely (97%)explained by the variations in contrasts of step length and floating distance(Fig. 3B). In other words,large strides are accomplished by taking large steps (r=0.80) and having long floating phases (r=0.37).

The variation in contrasts of step length is largely explained by the variations in contrasts of mass and hindlimb length(Fig. 3B). The variation in contrasts of SVL does not contribute significantly to the variation in contrasts of step length. Species that take large steps have relatively long hindlimbs (r=1.92) and are relatively light (r=-1.45).

Only the variation in contrasts of hindlimb length contributes significantly to the variation in contrasts of floating distance(Fig. 3B). The variations in contrasts of SVL and mass do not contribute significantly to the variation in contrasts of floating distance. Thus, species with long hindlimbs have long floating phases.

Variation in none of the contrasts of the three morphological variables(SVL, mass and hindlimb length) explains the variation in contrasts of stride frequency (r=0.50, F3,7=0.76, P=0.55).

Speed modulation

The variation in contrasts of maximal sprint speed is marginally explained(P=0.075) by the variation in contrasts of the slopes of stride length against speed (Fig. 3C). Because the contrasts of the slopes of stride frequency and stride length against running speed are inversely related (by definition), the opposite relationship is obtained for the contrasts of the slopes of stride frequency and maximal sprint speed. This means that species attain a higher maximal sprint speed if they modulate their speed mainly by increasing stride length rather than stride frequency.

The variation in contrasts of the slopes of stride length, in turn, is explained by the variations in contrasts of the slopes of step length and floating distance against speed (Fig. 3C). Lizards modulate their stride length by changing both step length (r=0.83) and, to a lesser extent, floating distance(r=0.37).

The variation in contrasts of SVL and hindlimb length explains the variation in contrasts of the slopes of step length against speed(Fig. 3C). The variation in contrasts of mass does not contribute to this variation. Species in which step length changes with speed have relatively long hindlimbs (r=2.30) and a relatively small SVL (r=-2.16).

The variation in contrasts of the slopes of floating distance against speed is explained by the variations in contrasts of mass and hindlimb length(Fig. 3C). The variation in contrasts of SVL does not contribute significantly to the variation in contrasts of the slopes of floating distance against speed. Species with relatively long hindlimbs (r=1.77) and a relatively small mass(r=-1.12) alter their floating distance with speed.

The variations in contrasts of mass and hindlimb length explain a significant proportion of the variation in contrasts of the slopes of stride frequency against speed (Fig. 3C): relatively heavy species (r=1.78) with relatively short hindlimbs (r=-2.14) alter stride frequency with speed. The variation in contrasts of SVL does not contribute significantly to the variation in contrasts of the slopes of stride frequency against speed.

The variation in contrasts of time spent in the open microhabitat is explained by the variation in contrasts of the slopes of stride frequency against speed (Fig. 4A). As mentioned above, the contrasts of the slopes of stride length and stride frequency are inversely related, so the opposite is true for the relationship between the contrasts of time spent in the open microhabitat and the contrasts of the slopes of stride length against running speed. This means that species occurring mostly in open microhabitats increase their speed mainly by increasing stride length and less by increasing stride frequency.

Fig. 4.

Bivariate correlations between the contrasts of microhabitat use and the contrasts of the slopes of stride frequency against speed (forced through the origin). (A) Relationship between the contrasts of time spent in the open microhabitat and the contrasts of the slopes of stride frequency against speed. (B) Correlation between the contrasts of time spent in the vegetated microhabitat and the contrasts of the slopes of stride frequency against speed. In both cases, including the contrasts of the slopes of stride length against speed gives the opposite result (r=0.76, P=0.017; r=-0.62, P=0.074, respectively).

Fig. 4.

Bivariate correlations between the contrasts of microhabitat use and the contrasts of the slopes of stride frequency against speed (forced through the origin). (A) Relationship between the contrasts of time spent in the open microhabitat and the contrasts of the slopes of stride frequency against speed. (B) Correlation between the contrasts of time spent in the vegetated microhabitat and the contrasts of the slopes of stride frequency against speed. In both cases, including the contrasts of the slopes of stride length against speed gives the opposite result (r=0.76, P=0.017; r=-0.62, P=0.074, respectively).

The variation in contrasts of time spent in the vertical microhabitat is not explained by the variation in contrasts of the slopes of stride frequency against speed (r=0.20, F1,7=0.30, P=0.60). The variation in contrasts of time spent in the vegetated microhabitat is, however, explained by the variation in contrasts of the slopes of stride frequency against speed(Fig. 4B). In other words,species occurring mostly amongst vegetation increase their speed mainly by increasing stride frequency and less by increasing stride length.

Among the 11 lacertid lizards under study here, it is clear that the variation in maximal sprint speed is determined by variation in both stride length and stride frequency. The highest sprint speeds are, however, achieved by those species that modulate their speed mainly by increasing stride length. This is not surprising because stride frequency cannot be increased indefinitely. The dynamics, such as the force/velocity relationships, of the muscles involved and the inertial properties of the musculo-skeletal system at the moment of maximal performance (i.e. maximal activation of muscles) set the potential maximal frequency. Because of the great overall resemblance between lacertid lizards, it can be assumed that the maximal attainable stride frequencies for the different species will be of the same order of magnitude. When comparing stride frequencies at maximal sprint speed among species, the lowest stride frequency is more than half the highest (the lowest stride length is one-third of the highest). If species were to employ a stride frequency modulation strategy, they would reach their respective stride frequency limit earlier (i.e. at a lower speed) than if they were to employ a stride length modulation strategy. In the latter case, they can therefore continue to increase sprint speed for longer, resulting in a higher maximal sprint performance.

Stride length, in turn, is determined by step length and floating distance. In this case, speed-dependent changes in step length have a greater effect on changes in stride length than does variation in floating distance with speed. Hindlimb length correlates positively (for a given body mass) both with step length and with floating distance. The positive effect of hindlimb length on step length is expected if the joint angles at touch-down and take-off are equal in short-limbed and long-limbed lizards. In this case, the distance the body moves forward during stance will be greater in species with longer limbs. Long-legged birds and mammals appear to be able to take longer steps(Hoyt et al., 2000), but Irschick and Jayne (1999)found no correlation between limb length and step length among five species of lizards.

The positive effect of hindlimb length on floating distance is in agreement with biomechanical predictions. A stride cycle that includes a floating phase can be regarded as two single leg jumps in sequence. The distance an animal can jump is determined, amongst other factors, by the distance through which it accelerates during push-off (Alexander,1972). In this case, the latter distance relates to hindlimb length. Thus, long-limbed animals are predicted to perform the longest jumps. Quantitative data on jumping performance in Anolis lizards(Losos, 1990) and tree frogs(Emerson, 1991) substantiate this prediction.

Among the 11 lacertid species under study, mass has a negative effect on step length for lizards with equally long hindlimbs. If we assume that muscle mass scales proportionally with body mass (m) and that force(F) scales proportionally with the cross-sectional area of the muscles, thus with m2/3, the acceleration (a) of the centre of mass will be lower for a given muscle activation pattern(F=ma). This will then result in a smaller forward displacement of the body and, thus, a shorter step length. Hoyt et al.(2000) also found that relative step length (i.e. expressed in units of leg length) is smaller in larger mammals and birds than in smaller ones.

None of the morphological variables that we introduced in the regression analysis explained a significant proportion of the interspecific variation in stride frequency at maximal sprint speed. Among seven teiid lizard species,however, stride frequency appeared to be negatively affected by body mass at body-length-equivalent speeds (White and Anderson, 1994), and stride frequency at burst running speed decreased with body mass in the lizard Dipsosaurus dorsalis(Johnson et al., 1993). In contrast, mass had no effect on stride frequency in painted turtles Chrysemys picta (Zani and Claussen, 1995). The physiological properties of the limb muscles may determine stride frequency. For example, limb cycling frequency increases with increasing body temperature in D. dorsalis(Swoap et al., 1993). It has been suggested that the positive relationship between frequency and temperature is due to the thermal dependence of muscle contraction speed (e.g. Bauwens et al., 1995). However,the seemingly straightforward link between muscle contractile properties and maximal stride frequency remains obscure because muscle cycling frequency for maximal power output largely exceeds maximal limb cycling frequency in D. dorsalis (Johnson et al.,1993; Swoap et al.,1993).

Speed modulation strategy during steady locomotion differs among lacertid species occurring in different microhabitats. Species that occur mainly in open microhabitats increase sprint speed mainly by increasing stride length,while species occurring in patches of vegetation mainly increase stride frequency. This is as expected. In an open microhabitat, species can include a floating phase in their stride cycle to increase stride length. A vegetated microhabitat, in contrast to an open one, poses conflicting demands on running style. Lizards occurring among patches of vegetation need to change direction fast and frequently. They can do so only when their feet are touching the ground. Therefore, it seems logical for these animals to increase speed mainly by increasing stride frequency. This reasoning is based on our data on steady-speed locomotion. In contrast, during unsteady locomotion, such as when accelerating, some lizard species are known to increase speed by increasing stride frequency initially and by subsequently increasing stride length(Irschick and Jayne,1998).

Each steady-running style (stride length versus stride frequency modulation) seems to place specific demands on the locomotor apparatus of the lizards. While stride-length-modulating species appear to be long-limbed (for a given mass or SVL) and light or small (for a given hindlimb length),stride-frequency-modulating species are short-limbed (for a given mass) and heavy (for a given hindlimb length). Long-limbed lizards are able to change stride length over a greater range because limb length sets the upper limit on step length (in theory, maximal step length equals twice leg length) and floating distance (as argued above). Because, at a given speed, stride length and frequency are inversely related, long-limbed lizards have lower stride frequencies at any speed than short-limbed animals.

It is unclear why stride-length-modulating species should be light or small and stride frequency modulators heavy (all relative to hindlimb length). In Mammalia, the opposite effect was found(Strang and Steudel, 1990):larger mammals change speed mainly by altering stride length, while smaller mammals do so mainly by increasing stride frequency. However, Strang and Steudel (1990) state that it is not clear how the interrelationship between size and hindlimb length affects their results. Therefore, the morphological basis for the differences in how species change their sprint speed is far from clear. Data are needed on the physiology and architecture of the musculoskeletal system to determine the factors that affect the way in which animals modulate sprint speed.

Species occurring mainly in vertical microhabitats do not seem to be either pronounced stride length or frequency modulators. They alter their speed by combining the two strategies. Neither stride length nor frequency can be indefinitely increased in this context. When climbing up vertical surfaces,animals need to establish a secure grip on the surface because slipping could be disastrous. This is in sharp contrast to locomotion on the level, where slipping does not constitute such a major hazard. Taking time to get a good grip and to position the legs precisely potentially constrains stride frequency. Taking larger strides in the same context can only be done by increasing step length; floating phases will cause the animals to fall from a vertical surface. We tested this idea by correlating the contrasts of the slopes of floating distance against speed with the contrasts of time spent in the vertical microhabitat. The relationship between the two variables was negative (r=-0.65; P=0.057) and therefore appears to substantiate our hypothesis. The inability to include floating phases into the stride cycles greatly limits stride length in lizards living on steep inclines.

By employing a stride length modulation strategy, species from open microhabitats are able to attain faster maximal sprint speeds than species occurring mainly in the other two microhabitats. Being able to escape fast in an open environment appears to be `ecologically relevant', while high sprint speeds might be less important in enclosed environments (see B. Vanhooydonck and R. Van Damme, in preparation). On the basis of these results, it is tempting to conclude that the speed modulation strategy of lizards has been shaped by natural selection to fulfil the specific requirements set by the environmental context. However, we realise that sprint speed might be confounding the relationship between microhabitat use and speed modulation strategy. The sole objective of a lizard living in open areas might be to attain high sprint speeds. The way to maximise level sprinting speed is to modulate speed mainly by increasing stride length (see above). Because sprint capacity is correlated with occurrence in an open microhabitat (B. Vanhooydonck and R. Van Damme, in preparation), the latter might be correlated indirectly to speed modulation strategy. Theoretically, this could be tested by comparing species with the same maximal sprint speeds but from different microhabitats. In practice, this is impossible because the lacertid lizards studied here seem to differ in locomotor capacity depending, amongst other factors, on their microhabitat use (B. Vanhooydonck and R. Van Damme, in preparation).

Fundamental to this study are the assumptions that natural selection moulds the way in which lizards adapt their speed to their life style (open,vegetated or vertical microhabitat) and that running style is fixed within species. In other words, we assumed that, if species use a running style`adapted' to life in their natural surroundings, the same animal will use the same speed modulation strategy in different contexts. Data on locomotor styles of level-running and climbing geckoes substantiate this assumption(Zaaf et al., 2001). However,to test explicitly whether running style in the lacertid lizards studied here is plastic, we need quantitative data on their speed modulation strategies on different, ecologically relevant, substrata (e.g. level versusinclined surfaces, open areas versus patches of vegetation).

We would like to thank Jan Scholliers for taking care of the animals. Ignacio Barandalla and Jan Scholliers helped with the digitisations. Duncan J. Irschick and an anonymous reviewer provided constructive comments on an earlier draft of the manuscript. B.V. was a scholarship student at the IWT(grant no. 951359) when most of this work was carried out. R.V.D. is a senior research assistant of the Science Fund — Flanders (FWO-VI). This research was supported financially by a GOA-BOF project of the University of Antwerp (1999-2002).

Alexander, R. McN. (
1972
).
Animal Mechanics.
Seattle: University of Washington Press.
Alexander, R. McN. (
1977a
). Terrestrial locomotion. In
Mechanics and Energetics of Animal Locomotion
(ed. R. McN. Alexander and G. Goldspink), pp.
168
-203. London: Chapman & Hall.
Alexander, R. McN. (
1977b
). Mechanics and scaling of terrestrial locomotion. In
Scale Effects of Animal Locomotion
(ed. T. J. Pedley), pp.
93
-110. New York, London: Academic Press.
Alexander, R. McN. (
1982
).
Locomotion of Animals.
Glasgow: Blackie.
Arnold, E. N. (
1989
). Towards a phylogeny and biogeography of the Lacertidae: relationships within an old-world family of lizards derived from morphology.
Bull. Brit. Mus. (Nat. Hist.)Zool.
55
,
209
-257.
Arnold, E. N. (
1998
). Structural niche, limb morphology and locomotion in lacertid lizards (Squamata, Lacertidae); a preliminary survey.
Bull. Brit. Mus. (Nat. Hist.)Zool.
64
,
63
-89.
Bauwens, D., Garland, T., Jr, Castilla, A. M. and Van Damme,R. (
1995
). Evolution of sprint speed in lacertid lizards:morphological, physiological and behavioral covariation.
Evolution
49
,
848
-863.
Bonine, K. E. and Garland, T., Jr (
1999
). Sprint performance of phrynosomatid lizards, measured on a high-speed treadmill, correlates with hindlimb length.
J. Zool.,Lond.
248
,
255
-265.
Castilla, A. M., Van Damme, R. and Bauwens, D.(
1999
). Field body temperatures, mechanisms of thermoregulation and evolution of thermal characteristics in lacertid lizards.
Natura Croatica
8
,
253
-274.
Díaz-Uriarte, R. and Garland, T., Jr(
1998
). Effects of branch lengths errors on the performance of phylogenetically independent contrasts.
Syst. Biol.
47
,
654
-672.
Emerson, S. B. (
1991
). The ecomorphology of Bornean tree frogs (family Rhacophoridae).
Zool. J. Linn. Soc.
101
,
337
-357.
Felsenstein, J. (
1985
). Phylogenies and the comparative method.
Am. Nat.
125
,
1
-15.
Felsenstein, J. (
1988
). Phylogenies and quantitative characters.
Annu. Rev. Ecol. Syst.
19
,
445
-471.
Fu, J. (
2000
). Toward the phylogeny of the family Lacertidae — Why 4708 base pairs of mtDNA sequences cannot draw the picture.
Biol. J. Linn. Soc.
71
,
203
-217.
Garland, T., Jr (
1983
). The relation between maximal running speed and body mass in terrestrial mammals.
J. Zool., Lond.
199
,
157
-170.
Garland, T., Jr (
1999
). Laboratory endurance capacity predicts variation in field locomotor behaviour among lizard species.
Anim. Behav.
57
,
77
-83.
Garland, T., Jr, Dickerman, A. W., Janis, C. M. and Jones, J. A. (
1993
). Phylogenetic analysis of covariance by computer simulation.
Syst. Biol.
42
,
265
-292.
Garland, T., Jr, Geiser, F. and Baudinette, R. V.(
1988
). Comparative locomotor performance of marsupial and placental mammals.
J. Zool., Lond.
215
,
505
-522.
Garland, T., Jr, Harvey, P. H. and Ives, A. R.(
1992
). Procedures for the analysis of comparative data using phylogenetically independent contrasts.
Syst. Biol.
41
,
18
-32.
Garland, T., Jr and Janis, C. M. (
1993
). Does metatarsal/femur ration predict maximal running speed in cursorial mammals?
J. Zool., Lond.
229
,
133
-151.
Garland, T., Jr and Losos, J. B. (
1994
). Ecological morphology of locomotor performance in squamate reptiles. In
Ecological Morphology. Integrative Organismal Biology
(ed. P. C. Wainwright and S. M. Reilly), pp.
240
-302. Chicago: The University of Chicago Press.
Garland, T., Jr., Midford, P. E. and Ives, A. R.(
1999
). An introduction to phylogenetically based statistical methods, with a new method for confidence intervals on ancestral values.
Am. Zool.
39
,
374
-388.
Harris, D. J. and Arnold, E. N. (
1999
). Relationships of wall lizards, Podarcis (Reptilia: Lacertidae) based on mitochondrial DNA sequences.
Copeia
1999
,
749
-754.
Harris, D. J., Arnold, E. N. and Thomas, R. H.(
1998
). Relationships of lacertid lizards (Reptilia: Lacertidae)estimated from mitochondrial DNA sequences and morphology.
Proc. R. Soc. Lond. B
265
,
1939
-1948.
Harvey, P. H. and Pagel, M. D. (
1991
).
The Comparative Method in Evolutionary Biology.
Oxford: Oxford University Press.
Heglund, N. C. and Taylor, C. R. (
1988
). Speed,stride frequency and energy cost per stride: how do they change with body size and gait?
J. Exp. Biol.
138
,
301
-318.
Hoyt, D. F., Wickler, S. J. and Cogger, E. A.(
2000
). Time of contact and step length: the effect of limb length, running speed, load carrying and incline.
J. Exp. Biol.
203
,
221
-227.
Irschick, D. J., Austin, C. C., Petren, K., Fisher, R. N.,Losos, J. B. and Ellers, O. (
1996
). A comparative analysis of clinging ability among pad-bearing lizards.
Biol. J. Linn. Soc.
59
,
21
-35.
Irschick, D. J. and Jayne, B. C. (
1998
). Effects of incline on speed, acceleration, body posture and hindlimb kinematics in two species of lizard Callisaurus draconoides and Uma scoparia.
J. Exp. Biol.
201
,
273
-287.
Irschick, D. J. and Jayne, B. C. (
1999
). Comparative three-dimensional kinematics of the hindlimb for high-speed bipedal and quadrupedal locomotion of lizards.
J. Exp. Biol.
202
,
1047
-1065.
Irschick, D. J. and Losos, J. B. (
1999
). Do lizards avoid habitats in which performance is submaximal? The relationship between sprinting capabilities and structural habitat use in Caribbean anoles.
Am. Nat.
154
,
293
-305.
Johnson, T. P., Swoap, S. J., Bennett, A. F. and Josephson, R. K. (
1993
). Body size, muscle power output and limitations on burst locomotor performance in the lizard Dipsosaurus dorsalis.
J. Exp. Biol.
174
,
199
-213.
Kolok, A. S. (
1999
). Interindividual variation in the prolonged locomotor performance of ectothermic vertebrates: a comparison of fish and herpetofaunal methodologies and a brief review of the recent fish literature.
Can. J. Fish. Aquat. Sci.
56
,
700
-710.
Losos, J. B. (
1990
). The evolution of form and function: morphology and locomotor performance in West Indian Anolislizards.
Evolution
44
,
1189
-1203.
Losos, J. B. and Miles, D. B. (
1994
). Adaptation, constraint and the comparative method: phylogenetic issues and methods. In
Ecological Morphology. Integrative Organismal Biology
(ed. P. C. Wainwright and S. M. Reilly), pp.
60
-98. Chicago: The University of Chicago Press.
Martins, E. P. and Garland, T., Jr (
1991
). Phylogenetic analyses of the correlated evolution of continuous characters: a simulation study.
Evolution
45
,
534
-557.
McMahon, T. A. (
1984
).
Muscles,Reflexes and Locomotion.
Princeton: Princeton University Press.
McPeek, M. A., Schrot, A. K. and Brown, J. M.(
1996
). Adaptation to predators in a new community: swimming performance and predator avoidance in damselflies.
Ecology
77
,
617
-629.
Melville, J. and Swain, R. (
2000
). Evolutionary relationships between morphology, performance and habitat openness in the lizard genus Niveoscincus (Scincidae: Lygosominae).
Biol. J. Linn. Soc.
70
,
667
-683.
Miles, D. B. (
1994
). Covariation between morphology and locomotory performance in Sceloporine lizards. In
Lizard Ecology. Historical and Experimental Perspectives
(ed. L. J. Vitt and E. R. Pianka), pp.
207
-236. Princeton: Princeton University Press.
Perry, G. (
1999
). The evolution of search modes: ecological versus phylogenetic perspectives.
Am. Nat.
153
,
98
-109.
Purvis, A. and Garland, T., Jr (
1993
). Polytomies in comparative analyses of continuous characters.
Syst. Biol.
42
,
569
-575.
Rayner, J. M. V. (
1995
). Flight mechanisms and constraints on flight performance.
Isr. J. Zool.
41
,
321
-342.
Shipley, L. A., Spallinger, D. E., Gross, J. E., Thompson, H. N. and Wunder, B. A. (
1996
). The dynamics and scaling of foraging velocity and encounter rate in mammalian herbivores.
Funct. Ecol.
10
,
234
-244.
Snyder, R. C. (
1952
). Quadrupedal and bipedal locomotion of lizards.
Copeia
1952
,
64
-70.
Sokal, R. R. and Rohlf, F. J. (
1995
).
Biometry. The Principles and Practice of Statistics in Biological Research.
Third edition. New York: W. H. Freeman &Company.
Srygley, R. B. and Dudley, R. (
1993
). Correlations of the position of center of body mass with butterfly escape tactics.
J. Exp. Biol.
174
,
155
-166.
Strang, K. T. and Steudel, K. (
1990
). Explaining the scaling of transport costs: the role of stride frequency and stride length.
J. Zool., Lond.
221
,
343
-358.
Swoap, S. J., Johnson, T. P., Josephson, R. K. and Bennett, A. F. (
1993
). Temperature, muscle power output and limitations on burst locomotor performance of the lizard Dipsosaurus dorsalis.
J. Exp. Biol.
174
,
185
-197.
Van Damme, R., Aerts, P. and Vanhooydonck, B.(
1998
). Variation in morphology, gait characteristics and speed of locomotion in two populations of lizards.
Biol. J. Linn. Soc.
63
,
409
-427.
Vanhooydonck, B. and Van Damme, R. (
2001
). Evolutionary trade-offs in locomotor capacities in lacertid lizards: are splendid sprinters clumsy climbers?
J. Evol. Biol.
14
,
46
-54.
Vanhooydonck, B., Van Damme, R. and Aerts, P.(
2001
). Speed and stamina trade-off in lacertid lizards.
Evolution
55
,
1040
-1048.
Walton, B. M. (
1993
). Physiology and phylogeny:the evolution of locomotor energetics in hylid frogs.
Am. Nat.
141
,
26
-50.
White, T. D. and Anderson, R. A. (
1994
). Locomotor patterns and costs as related to body size and form in teiid lizards.
J. Zool., Lond.
233
,
107
-128.
Zaaf, A., Van Damme, R., Herrel, A. and Aerts, P.(
2001
). Spatio-temporal gait characteristics of level and vertical locomotion in a ground-dwelling and a climbing gecko.
J. Exp. Biol.
204
,
1233
-1246.
Zani, P. A. and Claussen, D. L. (
1995
). Effects of extrinsic load on locomotion in painted turtles (Chrysemys picta).
Copeia
1995
,
735
-738.